Reproductive character displacement and environmental filtering shape floral variation between sympatric sister taxa

Divergence in reproductive traits between closely related species that co‐occur contributes to speciation by reducing interspecific gene flow. In flowering plants, greater floral divergence in sympatry than allopatry may reflect reproductive character displacement (RCD) by means of divergent pollinator‐mediated selection or mating system evolution. However, environmental filtering (EF) would prevail for floral traits under stronger selection by abiotic factors than pollination, and lead to sympatric taxa being more phenotypically similar. We determine whether floral UV pigmentation and size show signatures of RCD or EF using a biogeographically informed sister taxa comparison. We determine whether 35 sister pairs in the Potentilleae tribe (Rosaceae) are allopatric or sympatric and confirm that sympatric sisters experience more similar bioclimatic conditions, an assumption of the EF hypothesis. We test whether interspecific differences are greater in allopatry or sympatry while accounting for divergence time. For UV pigmentation, sympatric sisters are more phenotypically similar than allopatric ones. For flower size, sympatric sisters show increased divergence with time since speciation but allopatric ones do not. We conclude that floral UV pigmentation shows a signature of EF, whereas flower size shows a signature of RCD. Discordant results between the traits suggest that the dominant selective agent differs between them.     

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.     

Differential divergence in autosomes and sex chromosomes is associated with intra‐island diversification at a very small spatial scale in a songbird lineage

Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome‐wide markers to investigate the diversification of the Reunion grey white‐eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping‐by‐sequencing and pooled RAD‐seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z‐linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process.     

The evolution of floral signals in relation to range overlap in a clade of California Jewelflowers (Streptanthus s.l.)

Because of their function as reproductive signals in plants, floral traits experience distinct selective pressures related to their role in speciation, reinforcement, and prolonged coexistence with close relatives. However, few studies have investigated whether population‐level processes translate into detectable signatures at the macroevolutionary scale. Here, we ask whether patterns of floral trait evolution and range overlap across a clade of California Jewelflowers reflect processes hypothesized to shape floral signal differentiation at the population level. We found a pattern of divergence in floral scent composition across the clade such that close relatives had highly disparate floral scents given their age. Accounting for range overlap with close relatives explained additional variation in floral scent over time, with sympatric species pairs having diverged more than allopatric species pairs given their age. However, three other floral traits (flower size, scent complexity and flower color) did not fit these patterns, failing to deviate from a null Brownian motion model of evolution. Together, our results suggest that selection for divergence among close relatives in the composition of floral scents may play a key, sustained role in mediating speciation and coexistence dynamics across this group, and that signatures of these dynamics may persist at the macroevolutionary scale.     

Lack of parallel genetic patterns underlying the repeated ecological divergence of beach and stream‐spawning kokanee salmon

Recent progress in methods for detecting adaptive population divergence in situ shows promise for elucidating the conditions under which selection acts to generate intraspecific diversity. Rapid ecological diversification is common in fishes; however, the role of phenotypic plasticity and adaptation to local environments is poorly understood. It is now possible to investigate genetic patterns to make inferences regarding phenotypic traits under selection and possible mechanisms underlying ecotype divergence, particularly where similar novel phenotypes have arisen in multiple independent populations. Here, we employed a bottom‐up approach to test for signatures of directional selection associated with divergence of beach‐ and stream‐spawning kokanee, the obligate freshwater form of sockeye salmon (Oncorhynchus nerka). Beach‐ and stream‐spawners co‐exist in many post‐glacial lakes and exhibit distinct reproductive behaviours, life‐history traits and spawning habitat preferences. Replicate ecotype pairs across five lakes in British Columbia, Canada were genotyped at 57 expressed sequence tag‐linked and anonymous microsatellite loci identified in a previous genome scan. Fifteen loci exhibited signatures of directional selection (high F_ST outliers), four of which were identified in multiple lakes. However, the lack of parallel genetic patterns across all lakes may be a result of: 1) an inability to detect loci truly under selection; 2) alternative genetic pathways underlying ecotype divergence in this system; and/or 3) phenotypic plasticity playing a formative role in driving kokanee spawning habitat differences. Gene annotations for detected outliers suggest pathogen resistance and energy metabolism as potential mechanisms contributing to the divergence of beach‐ and stream‐spawning kokanee, but further study is required.     

COMPETITION FOR HUMMINGBIRD POLLINATION SHAPES FLOWER COLOR VARIATION IN ANDEAN SOLANACEAE

One classic explanation for the remarkable diversity of flower colors across angiosperms involves evolutionary shifts among different types of pollinators with different color preferences. However, the pollinator shift model fails to account for the many examples of color variation within clades that share the same pollination system. An alternate explanation is the competition model, which suggests that color divergence evolves in response to interspecific competition for pollinators, as a means to decrease interspecific pollinator movements. This model predicts color overdispersion within communities relative to null assemblages. Here, we combine morphometric analyses, field surveys, and models of pollinator vision with a species‐level phylogeny to test the competition model in the primarily hummingbird‐pollinated clade Iochrominae (Solanaceae). Results show that flower color as perceived by pollinators is significantly overdispersed within sites. This pattern is not simply due to phylogenetic history: phylogenetic community structure does not deviate from random expectations, and flower color lacks phylogenetic signal. Moreover, taxa that occur in sympatry occupy a significantly larger volume of color space than those in allopatry, supporting the hypothesis that competition in sympatry drove the evolution of novel colors. We suggest that competition among close relatives may commonly underlie floral divergence, especially in species‐rich habitats where congeners frequently co‐occur.     

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.     

Toxicity and population structure of the Rough‐Skinned Newt (Taricha granulosa) outside the range of an arms race with resistant predators

Species interactions, and their fitness consequences, vary across the geographic range of a coevolutionary relationship. This spatial heterogeneity in reciprocal selection is predicted to generate a geographic mosaic of local adaptation, wherein coevolutionary traits are phenotypically variable from one location to the next. Under this framework, allopatric populations should lack variation in coevolutionary traits due to the absence of reciprocal selection. We examine phenotypic variation in tetrodotoxin (TTX) toxicity of the Rough‐Skinned Newt (Taricha granulosa) in regions of allopatry with its TTX‐resistant predator, the Common Garter Snake (Thamnophis sirtalis). In sympatry, geographic patterns of phenotypic exaggeration in toxicity and toxin‐resistance are closely correlated in prey and predator, implying that reciprocal selection drives phenotypic variation in coevolutionary traits. Therefore, in allopatry with TTX‐resistant predators, we expect to find uniformly low levels of newt toxicity. We characterized TTX toxicity in northwestern North America, including the Alaskan panhandle where Ta. granulosa occur in allopatry with Th. sirtalis. First, we used microsatellite markers to estimate population genetic structure and determine if any phenotypic variation in toxicity might be explained by historical divergence. We found northern populations of Ta. granulosa generally lacked population structure in a pattern consistent with northern range expansion after the Pleistocene. Next, we chose a cluster of sites in Alaska, which uniformly lacked genetic divergence, to test for phenotypic divergence in toxicity. As predicted, overall levels of newt toxicity were low; however, we also detected unexpected among‐ and within‐population variation in toxicity. Most notably, a small number of individuals contained large doses of TTX that rival means of toxic populations in sympatry with Th. sirtalis. Phenotypic variation in toxicity, despite limited neutral genetic divergence, suggests that factors other than reciprocal selection with Th. sirtalis likely contribute to geographic patterns of toxicity in Ta. granulosa.     

Geographic variation in hybridization across a reinforcement contact zone of chorus frogs (Pseudacris)

Reinforcement contact zones, which are secondary contact zones where species are diverging in reproductive behaviors due to selection against hybridization, represent natural laboratories for studying speciation‐in‐action. Here, we examined replicate localities across the entire reinforcement contact zone between North American chorus frogs Pseudacris feriarum and P. nigrita to investigate geographic variation in hybridization frequencies and to assess whether reinforcement may have contributed to increased genetic divergence within species. Previous work indicated these species have undergone reproductive character displacement (RCD) in male acoustic signals and female preferences due to reinforcement. We also examined acoustic signal variation across the contact zone to assess whether signal characteristics reliably predict hybrid index and to elucidate whether the degree of RCD predicts hybridization rate. Using microsatellites, mitochondrial sequences, and acoustic signal information from >1,000 individuals across >50 localities and ten sympatric focal regions, we demonstrate: (1) hybridization occurs and (2) varies substantially across the geographic range of the contact zone, (3) hybridization is asymmetric and in the direction predicted from observed patterns of asymmetric RCD, (4) in one species, genetic distance is higher between conspecific localities where one or both have been reinforced than between nonreinforced localities, after controlling for geographic distance, (5) acoustic signal characters strongly predict hybrid index, and (6) the degree of RCD does not strongly predict admixture levels. By showing that hybridization occurs in all sympatric localities, this study provides the fifth and final line of evidence that reproductive character displacement is due to reinforcement in the chorus frog contact zone. Furthermore, this work suggests that the dual action of cascade reinforcement and partial geographic isolation is promoting genetic diversification within one of the reinforced species.     

Sperm divergence in a passerine contact zone: Indication of reinforcement at the gametic level

Postcopulatory sexual selection may promote evolutionary diversification in sperm form, but the contribution of between‐species divergence in sperm morphology to the origin of reproductive isolation and speciation remains little understood. To assess the possible role of sperm diversification in reproductive isolation, we studied sperm morphology in two closely related bird species, the common nightingale (Luscinia megarhynchos) and the thrush nightingale (Luscinia luscinia), that hybridize in a secondary contact zone spanning Central and Eastern Europe. We found: (1) striking divergence between the species in total sperm length, accompanied by a difference in the length of the mitochondrial sperm component; (2) greater divergence between species in sperm morphology in sympatry than in allopatry, with evidence for character displacement in sperm head length detected in L. megarhynchos; (3) interspecific hybrids showing sperm with a length intermediate between the parental species, but no evidence for decreased sperm quality (the proportion of abnormal spermatozoa in ejaculates). Our results demonstrate that divergence in sperm morphology between the two nightingale species does not result in intrinsic postzygotic isolation, but may contribute to postcopulatory prezygotic isolation. This isolation could be strengthened in sympatry by reinforcement.     

Floral color variation and associations with fitness‐related traits in Malva moschata (Malvaceae)

Genetic polymorphisms for floral color are interesting phenomena to study because they are likely to be maintained by opposing selective forces. Pollinator preferences may exert direct selection on floral color; however, floral color might also be the indirect target of selection through genetic associations with other traits under selection. Malva moschata (Malvaceae) is a North American species that produces either red or white flowers. In the present study, we present reflectance spectrophotometry data that characterize the nature of floral color variation in this species and show that honey bees and bumble bees should be able to distinguish between the morphs through differential sensitivity at the green (long‐wavelength) photoreceptor. Second, we use a series of phenotypic measures to investigate whether the color morphs differ with respect to other floral traits, vegetative traits or female reproductive success, and use a series of correlation analyses to infer the relative independence of color from these other traits. We found that red‐flowered morphs produced more anthers per flower and had greater leaf area, and that white‐flowered morphs had greater percentage fruit set; however, there were no reproductive success differences between the morphs. The relationship between flower size and anther number was the only correlation that differed between the morphs. Finally, a series of pollinator‐choice experiments showed that bumble bees strongly prefer red morphs in terms of visit frequency and duration, but honey bees have no preference. Taken together, our results suggest that color is rather independent of other phenotypic traits, and that honey bee abundance is likely to play a role in maintaining color variation in this system.     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.     

Selection and demographic history shape the molecular evolution of the gamete compatibility protein bindin in Pisaster sea stars

Reproductive compatibility proteins have been shown to evolve rapidly under positive selection leading to reproductive isolation, despite the potential homogenizing effects of gene flow. This process has been implicated in both primary divergence among conspecific populations and reinforcement during secondary contact; however, these two selective regimes can be difficult to discriminate from each other. Here, we describe the gene that encodes the gamete compatibility protein bindin for three sea star species in the genus Pisaster. First, we compare the full‐length bindin‐coding sequence among all three species and analyze the evolutionary relationships between the repetitive domains of the variable second bindin exon. The comparison suggests that concerted evolution of repetitive domains has an effect on bindin divergence among species and bindin variation within species. Second, we characterize population variation in the second bindin exon of two species: We show that positive selection acts on bindin variation in Pisaster ochraceus but not in Pisaster brevispinus, which is consistent with higher polyspermy risk in P. ochraceus. Third, we show that there is no significant genetic differentiation among populations and no apparent effect of sympatry with congeners that would suggest selection based on reinforcement. Fourth, we combine bindin and cytochrome c oxidase 1 data in isolation‐with‐migration models to estimate gene flow parameter values and explore the historical demographic context of our positive selection results. Our findings suggest that positive selection on bindin divergence among P. ochraceus alleles can be accounted for in part by relatively recent northward population expansions that may be coupled with the potential homogenizing effects of concerted evolution.     

Reproductive character displacement of female,but not male song discrimination in an avian hybrid zone

Divergence of male sexual signals and female preferences for those signals often maintains reproductive boundaries between closely related, co‐occurring species. However, contrasting sources of selection, such as interspecific competition, can lead to weak divergence or even convergence of sexual signals in sympatry. When signals converge, assortative mating can be maintained if the mating preferences of females diverge in sympatry (reproductive character displacement; RCD), but there are few explicit examples. Pied flycatchers (Ficedula hypoleuca) are sympatric with collared flycatchers (F. albicollis) on the Baltic island of Öland, where males from both species compete over nestboxes, their songs converge, and the two species occasionally hybridize. We compare song discrimination of male and female pied flycatchers on Öland and in an allopatric population on the Swedish mainland. Using field choice trials, we show that male pied flycatchers respond similarly to the songs of both species in sympatry and allopatry, while female pied flycatchers express stronger discrimination against heterospecific songs in sympatry than in allopatry. These results are consistent with RCD of song discrimination of female pied flycatchers where they co‐occur with collared flycatchers, which should maintain species assortative mating despite convergence of male sexual signals.     

Maintenance of genetic and morphological identity in two sibling Syrrhopodon species (Calymperaceae,Bryopsida) despite extensive introgression

Bryophytes are a group of land plants in which the role of hybridization has long been challenged. Using genotyping by sequencing to circumvent the lack of molecular variation at selected loci previously used for phylogeny and morphology, we determine the level of genetic and morphological divergence and reproductive isolation between the sibling Syrrhopodon annotinus and S. simmondsii (Calymperaceae, Bryopsida) that occur in sympatry but in different habitats in lowland Amazonian rainforests. A clear morphological differentiation and a low (0.06), but significant Fst derived from the analysis of 183 single nucleotide polymorphisms were observed between the two species. Conspecific pairs of individuals consistently exhibited higher average kinship coefficients along a gradient of geographic isolation than interspecific pairs. The weak, but significant genetic divergence observed is consistent with growing evidence that ecological specialization can lead to genetic differentiation among bryophyte species. Nevertheless, the spatial genetic structures of the two species were significantly correlated, as evidenced by the significant slope of the Mantel test based on kinship coefficients between pairs of interspecific individuals and the geographic distance separating them. Interspecific pairs of individuals are thus more closely related when they are geographically closer, suggesting that isolation‐by‐distance is stronger than the interspecific reproductive barrier and pointing to interspecific gene flow. We conclude that interspecific introgression, whose role has long been questioned in bryophytes, may take place even in species wherein sporophyte production is scarce due to dioicy, raising the question as to what mechanisms maintain differentiation despite weak reproductive isolation.     

The extent and genetic basis of phenotypic divergence in life history traits in Mimulus guttatus

Differential natural selection acting on populations in contrasting environments often results in adaptive divergence in multivariate phenotypes. Multivariate trait divergence across populations could be caused by selection on pleiotropic alleles or through many independent loci with trait‐specific effects. Here, we assess patterns of association between a suite of traits contributing to life history divergence in the common monkey flower, Mimulus guttatus, and examine the genetic architecture underlying these correlations. A common garden survey of 74 populations representing annual and perennial strategies from across the native range revealed strong correlations between vegetative and reproductive traits. To determine whether these multitrait patterns arise from pleiotropic or independent loci, we mapped QTLs using an approach combining high‐throughput sequencing with bulk segregant analysis on a cross between populations with divergent life histories. We find extensive pleiotropy for QTLs related to flowering time and stolon production, a key feature of the perennial strategy. Candidate genes related to axillary meristem development colocalize with the QTLs in a manner consistent with either pleiotropic or independent QTL effects. Further, these results are analogous to previous work showing pleiotropy‐mediated genetic correlations within a single population of M. guttatus experiencing heterogeneous selection. Our findings of strong multivariate trait associations and pleiotropic QTLs suggest that patterns of genetic variation may determine the trajectory of adaptive divergence.     

Divergence and isolation of cryptic sympatric taxa within the annual legume Amphicarpaea bracteata

The amphicarpic annual legume Amphicarpaea bracteata is unusual in producing aerial and subterranean cleistogamous flowers that always self‐fertilize and, less commonly, aerial chasmogamous flowers that outcross. Although both morphologic and genetic variants are known in this highly selfing species, debate continues over whether this variation is continuous, reflecting the segregation of standing genetic variation, or discontinuous, reflecting distinct taxa that rarely intercross. We characterized SNP variation in 128 individuals in southern Wisconsin to assess within‐ and among‐population variation at 3928 SNPs. We also assessed genotype and leaf morphology in an additional 76 individuals to connect phenotypic variation with genetic variation. Genetic variation maps onto three strongly divergent and highly inbred genetic groups showing little relation to site location. Each group has a distinct phenotype, but the divergence of these groups differs from the varietal divisions previously identified based on morphological characters. Like previous authors, we argue that the taxonomy of this species should be revised. Despite extensive sympatry, estimates of among‐group migration rates are low, and hybrid individuals were at low frequency (<2%) in our dataset. Restricted gene flow likely results from high selfing rates and partial reproductive incompatibility as evidenced by the U‐shaped distribution of pairwise F_ST values reflecting “islands” of genomic divergence. These islands may be associated with hybrid incompatibility loci that arose in allopatry. The coexistence of lineages within sites may reflect density‐dependent attack by species‐specific strains of pathogenic fungi and/or root‐nodulating bacteria specializing on distinct genotypes.     

The birds, the bees, and the virtual flowers: can pollinator behavior drive ecological speciation in flowering plants?

Biologists have long assumed that pollinator behavior is an important force in angiosperm speciation, yet there is surprisingly little direct evidence that floral preferences in pollinators can drive floral divergence and the evolution of reproductive (ethological) isolation between incipient plant species. In this study, we expose computer-generated plant populations with a wide variation in flower color to selection by live and virtual hummingbirds and bumblebees and track evolutionary changes in flower color over multiple generations. Flower color, which was derived from the known genetic architecture and phenotypic variance of naturally occurring plant species pollinated by both groups, evolved in simulations through a genetic algorithm in which pollinator preference determined changes in flower color between generations. The observed preferences of live hummingbirds and bumblebees were strong enough to cause adaptive divergence in flower color between plant populations but did not lead to ethological isolation. However, stronger preferences assigned to virtual pollinators in sympatric and allopatric scenarios rapidly produced ethological isolation. Pollinators can thus drive ecological speciation in flowering plants, but more rigorous and comprehensive behavioral studies are required to specify conditions that produce sufficient preference levels in pollinators.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.