Cultural isolation is greater than genetic isolation across an avian hybrid zone

Elucidating the relationship between genetic and cultural evolution is important in understanding speciation, as learned premating barriers might be involved in maintaining species differences. Here, we test this relationship by examining a widely recognized premating barrier, bird song, in a hybrid zone between black‐throated green (Setophaga virens) and Townsend's warblers (S. townsendi). We use song analysis, genomic techniques and playback experiments to characterize the cultural and genetic backgrounds of individuals in this region, expecting that if song is an important reproductive barrier between these species, there should be a strong relationship between song and genotype. We show that songs in the hybrid zone correspond to the distinctly different songs found in allopatry but that song and genotype are not tightly coupled in sympatry. Allopatric individuals responded only to local songs, indicating that individuals may have learned to respond to songs they commonly hear. We observed discordance between song and genotype clines; a narrower cline suggests that cultural selection on song is stronger than natural selection on genotype. These findings indicate that song is unlikely to play a role in reproductive isolation between these species, and we suggest that spatial variation in song may nonetheless be maintained by frequency‐dependent cultural selection. This decoupling of genes and culture may contribute to hybridization in this region.     

Sexual selection and population divergence III: Interspecific and intraspecific variation in mating signals

A major challenge for studying the role of sexual selection in divergence and speciation is understanding the relative influence of different sexually selected signals on those processes in both intra‐ and interspecific contexts. Different signals may be more or less susceptible to co‐option for species identification depending on the balance of sexual and ecological selection acting upon them. To examine this, we tested three predictions to explain geographic variation in long‐ versus short‐range sexual signals across a 3,500 + km transect of two related Australian field cricket species (Teleogryllus spp.): (a) selection for species recognition, (b) environmental adaptation and (c) stochastic divergence. We measured male calling song and male and female cuticular hydrocarbons (CHCs) in offspring derived from wild populations, reared under common garden conditions. Song clearly differentiated the species, and no hybrids were observed suggesting that hybridization is rare or absent. Spatial variation in song was not predicted by geography, genetics or climatic factors in either species. In contrast, CHC divergence was strongly associated with an environmental gradient supporting the idea that the climatic environment selects more directly upon these chemical signals. In light of recently advocated models of diversification via ecological selection on secondary sexual traits, the different environmental associations we found for song and CHCs suggest that the impact of ecological selection on population divergence, and how that influences speciation, might be different for acoustic versus chemical signals.     

The genetics of speciation: genes of small effect underlie sexual isolation in the Hawaiian cricket Laupala

Sexual behaviours often evolve rapidly and are critical for sexual isolation. We suggest that coordinated sexual signals and preferences generate stabilizing selection, favouring the accumulation of many small‐effect mutations in sexual communication traits. Rapid radiation of a sexual behaviour used in signalling, song pulse rate, has been observed in the Hawaiian cricket genus Laupala. Using marker‐assisted introgression, we isolated five known quantitative trait loci (QTL) influencing species‐level differences in pulse rate from one species, L. paranigra, into a closely related species, L. kohalensis. All five QTL were found to have a significant effect on song and appear to be largely additive in backcross introgression lines. Furthermore, all effect sizes were small in magnitude. Our data provide support for the hypothesis that stabilizing selection on sexual signals in Laupala creates genetic conditions favourable to incremental divergence during speciation, through the evolution of alleles of minor rather than major phenotypic effects.     

Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S. aethnensis and S. chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome‐wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly‐differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs’ allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16–0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species – pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.     

A population genomic scan in Chorthippus grasshoppers unveils previously unknown phenotypic divergence

Understanding the genetics of speciation and the processes that drive it is a central goal of evolutionary biology. Grasshoppers of the Chorthippus species group differ strongly in calling song (and corresponding female preferences) but are exceedingly similar in other characteristics such as morphology. Here, we performed a population genomic scan on three Chorthippus species (Chorthippus biguttulus, C. mollis and C. brunneus) to gain insight into the genes and processes involved in divergence and speciation in this group. Using an RNA‐seq approach, we examined functional variation between the species by calling SNPs for each of the three species pairs and using F_ST‐based approaches to identify outliers. We found approximately 1% of SNPs in each comparison to be outliers. Between 37% and 40% of these outliers were nonsynonymous SNPs (as opposed to a global level of 17%) indicating that we recovered loci under selection. Among the outliers were several genes that may be involved in song production and hearing as well as genes involved in other traits such as food preferences and metabolism. Differences in food preferences between species were confirmed with a behavioural experiment. This indicates that multiple phenotypic differences implicating multiple evolutionary processes (sexual selection and natural selection) are present between the species.     

Mutation‐order divergence by sexual selection: diversification of sexual signals in similar environments as a first step in speciation

The origin of species remains a central question, and recent research focuses on the role of ecological differences in promoting speciation. Ecological differences create opportunities for divergent selection (i.e. ‘ecological’ speciation), a Darwinian hypothesis that hardly requires justification. In contrast, ‘mutation‐order’ speciation proposes that, instead of adapting to different environments, populations find different ways to adapt to similar environments, implying that speciation does not require ecological differences. This distinction is critical as it provides an alternative hypothesis to the prevailing view that ecological differences drive speciation. Speciation by sexual selection lies at the centre of debates about the importance of ecological differences in promoting speciation; here, we present verbal and mathematical models of mutation‐order divergence by sexual selection. We develop three general cases and provide a two‐locus population genetic model for each. Results indicate that alternative secondary sexual traits can fix in populations that initially experience similar natural and sexual selection and that divergent traits and preferences can remain stable in the face of low gene flow. This stable divergence can facilitate subsequent divergence that completes or reinforces speciation. We argue that a mutation‐order process could explain widespread diversity in secondary sexual traits among closely related, allopatric species.     

Interspecific competition promotes habitat and morphological divergence in a secondary contact zone between two hybridizing songbirds

Interspecific competition is assumed to play an important role in the ecological differentiation of species and speciation. However, empirical evidence for competition's role in speciation remains surprisingly scarce. Here, we studied the role of interspecific competition in the ecological differentiation and speciation of two closely related songbird species, the Common Nightingale (Luscinia megarhynchos) and the Thrush Nightingale (Luscinia luscinia). Both species are insectivorous and ecologically very similar. They hybridize in a secondary contact zone, which is a mosaic of sites where both species co‐occur (syntopy) and sites where only one species is present (allotopy). We analysed fine‐scale habitat data for both species in both syntopic and allotopic sites and looked for associations between habitat use and bill morphology, which have been previously shown to be more divergent in sympatry than in allopatry. We found that the two nightingale species differ in habitat use in allotopic sites, where L. megarhynchos occurred in drier habitats and at slightly higher elevations, but not in syntopic sites. Birds from allotopic sites also showed higher interspecific divergence in relative bill size compared to birds from syntopic sites. Finally, we found an association between bill morphology and elevation. Our results are consistent with the view that interspecific competition in nightingales has resulted in partial habitat segregation in sympatry and that the habitat‐specific food supply has in turn very likely led to bill size divergence. Such ecological divergence may enhance prezygotic as well as extrinsic postzygotic isolation and thus accelerate the completion of the speciation process.     

Bird song,ecology and speciation

The study of bird song dialects was once considered the most promising approach for investigating the role of behaviour in reproductive divergence and speciation. However, after a series of studies yielding conflicting results, research in the field slowed significantly. Recent findings, on how ecological factors may lead to divergence in both song and morphology, necessitate a re-examination. We focus primarily on species with learned song, examine conflicting results in the literature and propose some potential new directions for future studies. We believe an integrative approach, including an examination of the role of ecology in divergent selection, is essential for gaining insight into the role of song in the evolution of assortative mating. Habitat-dependent selection on both song and fitness-related characteristics can lead to parallel divergence in these traits. Song may, therefore, provide females with acoustic cues to find males that are most fit for a particular habitat. In analysing the role of song learning in reproductive divergence, we focus on post-dispersal plasticity in a conceptual framework. We argue that song learning may initially constrain reproductive divergence, while in the later stages of population divergence it may promote speciation.     

Isolation by distance,not incipient ecological speciation,explains genetic differentiation in an Andean songbird (Aves: Furnariidae: Cranioleuca antisiensis,Line‐cheeked Spinetail) despite near threefold body size change across an environmental gradient

During the process of ecological speciation, reproductive isolation results from divergent natural selection and leads to a positive correlation between genetic divergence and adaptive phenotypic divergence, that is, isolation by adaptation (IBA). In natural populations, phenotypic differentiation is often autocorrelated with geographic distance, making IBA difficult to distinguish from the neutral expectation of isolation by distance (IBD). We examined these two alternatives in a dramatic case of clinal phenotypic variation in an Andean songbird, the Line‐cheeked Spinetail (Cranioleuca antisiensis). At its geographic extremes, this species shows a near threefold difference in body mass (11.5 to 31.0 g) with marked plumage differences. We analysed phenotypic, environmental and genetic data (5,154 SNPs) from 172 individuals and 19 populations sampled along its linear distribution in the Andes. We found that body mass was tightly correlated with environmental temperature, consistent with local adaptation as per Bergmann's rule. Using a P_ST–F_ST analysis, we found additional support for natural selection driving body mass differentiation, but these results could also be explained by environment‐mediated phenotypic plasticity. When we assessed the relative support for patterns of IBA and IBD using variance partitioning, we found that IBD was the best explanation for genetic differentiation along the cline. Adaptive phenotypic or environmental divergence can reduce gene flow, a pattern interpreted as evidence of ecological speciation's role in diversification. Our results provide a counterexample to this interpretation. Despite conditions conducive to ecological speciation, our results suggest that dramatic size and environmental differentiation within C. antisiensis are not limiting gene flow.     

Model‐based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi‐permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome‐wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post‐glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.     

THE IMPORTANCE OF FEMALE CHOICE,MALE–MALE COMPETITION,AND SIGNAL TRANSMISSION AS CAUSES OF SELECTION ON MALE MATING SIGNALS

Selection on advertisement signals arises from interacting sources including female choice, male–male competition, and the communication channel (i.e., the signaling environment). To identify the contribution of individual sources of selection, we used previously quantified relationships between signal traits and each putative source to predict relationships between signal variation and fitness in Enchenopa binotata treehoppers (Hemiptera: Membracidae). We then measured phenotypic selection on signals and compared predicted and realized relationships between signal traits and mating success. We recorded male signals, then measured lifetime mating success at two population densities in a realistic environment in which sources of selection could interact. We identified which sources best predicted the relationship between signal variation and mating success using a multiple regression approach. All signal traits were under selection in at least one of the two breeding seasons measured, and in some cases selection was variable between years. Female preference was the strongest source of selection shaping male signals. The E. binotata species complex is a model of ecological speciation initiated by host shifts. Signal and preference divergence contribute to behavioral isolation within the complex, and the finding that female mate preferences drive signal evolution suggests that speciation in this group results from both ecological divergence and sexual selection.     

Panmixia supports divergence with gene flow in Darwin's small ground finch, Geospiza fuliginosa, on Santa Cruz, Galápagos Islands

The divergence‐with‐gene‐flow model of speciation has a strong theoretical basis with a growing number of plausible examples in nature, but remains hotly debated. Darwin’s finches of the Galápagos Archipelago have played an important role in our understanding of speciation processes. Recent studies suggest that this group may also provide insights into speciation via divergence with gene flow. On the island of Santa Cruz, recent studies found evidence for adaptive divergence in Darwin’s small ground finch, Geospiza fuliginosa, between ecologically contrasting arid and humid zones. Despite the short geographical distance between these zones, strong disruptive selection during low rainfall periods is expected to generate and maintain adaptive divergence. Conversely, during high rainfall periods, when disruptive selection is predicted to be weakened, population divergence in adaptive traits is expected to break down. Because periods of low and high rainfall irregularly alternate, the geographical pattern of adaptive divergence can be assumed to break down and, importantly, regenerate in situ. Here, we use microsatellite allele frequency data to assess the genetic population structure of G. fuliginosa on Santa Cruz. We sample 21 sites and four ecological zones across the island. We reject hypotheses of population substructure linked to ecological and geographical differences among sites in favour of a single panmictic population. Panmixia implies high levels of gene flow within Santa Cruz, which favours selection over genetic drift as a valid process generating phenotypic divergence in G. fuliginosa on Santa Cruz. We discuss how our findings may support classic adaptation, phenotypic plasticity, matching habitat choice or any combination of these three processes.     

The genomic bases of morphological divergence and reproductive isolation driven by ecological speciation in Senecio (Asteraceae)

Ecological speciation, driven by adaptation to contrasting environments, provides an attractive opportunity to study the formation of distinct species, and the role of selection and genomic divergence in this process. Here, we focus on a particularly clear‐cut case of ecological speciation to reveal the genomic bases of reproductive isolation and morphological differences between closely related Senecio species, whose recent divergence within the last ~200 000 years was likely driven by the uplift of Mt. Etna (Sicily). These species form a hybrid zone, yet remain morphologically and ecologically distinct, despite active gene exchange. Here, we report a high‐density genetic map of the Senecio genome and map hybrid breakdown to one large and several small quantitative trait loci (QTL). Loci under diversifying selection cluster in three 5 cM regions which are characterized by a significant increase in relative (F_ST), but not absolute (d_XY), interspecific differentiation. They also correspond to some of the regions of greatest marker density, possibly corresponding to ‘cold‐spots’ of recombination, such as centromeres or chromosomal inversions. Morphological QTL for leaf and floral traits overlap these clusters. We also detected three genomic regions with significant transmission ratio distortion (TRD), possibly indicating accumulation of intrinsic genetic incompatibilities between these recently diverged species. One of the TRD regions overlapped with a cluster of high species differentiation, and another overlaps the large QTL for hybrid breakdown, indicating that divergence of these species may have occurred due to a complex interplay of ecological divergence and accumulation of intrinsic genetic incompatibilities.     

Phenotypic variation and covariation indicate high evolvability of acoustic communication in crickets

Studying the genetic architecture of sexual traits provides insight into the rate and direction at which traits can respond to selection. Traits associated with few loci and limited genetic and phenotypic constraints tend to evolve at high rates typically observed for secondary sexual characters. Here, we examined the genetic architecture of song traits and female song preferences in the field crickets Gryllus rubens and Gryllus texensis. Song and preference data were collected from both species and interspecific F1 and F2 hybrids. We first analysed phenotypic variation to examine interspecific differentiation and trait distributions in parental and hybrid generations. Then, the relative contribution of additive and additive‐dominance variation was estimated. Finally, phenotypic variance–covariance ( P ) matrices were estimated to evaluate the multivariate phenotype available for selection. Song traits and preferences had unimodal trait distributions, and hybrid offspring were intermediate with respect to the parents. We uncovered additive and dominance variation in song traits and preferences. For two song traits, we found evidence for X‐linked inheritance. On the one hand, the observed genetic architecture does not suggest rapid divergence, although sex linkage may have allowed for somewhat higher evolutionary rates. On the other hand, P matrices revealed that multivariate variation in song traits aligned with major dimensions in song preferences, suggesting a strong selection response. We also found strong covariance between the main traits that are sexually selected and traits that are not directly selected by females, providing an explanation for the striking multivariate divergence in male calling songs despite limited divergence in female preferences.     

Divergent selection and drift shape the genomes of two avian sister species spanning a saline–freshwater ecotone

The role of species divergence due to ecologically based divergent selection—or ecological speciation—in generating and maintaining biodiversity is a central question in evolutionary biology. Comparison of the genomes of phylogenetically related taxa spanning a selective habitat gradient enables discovery of divergent signatures of selection and thereby provides valuable insight into the role of divergent ecological selection in speciation. Tidal marsh ecosystems provide tractable opportunities for studying organisms' adaptations to selective pressures that underlie ecological divergence. Sharp environmental gradients across the saline–freshwater ecotone within tidal marshes present extreme adaptive challenges to terrestrial vertebrates. Here, we sequence 20 whole genomes of two avian sister species endemic to tidal marshes—the saltmarsh sparrow (Ammospiza caudacutus) and Nelson's sparrow (A. nelsoni)—to evaluate the influence of selective and demographic processes in shaping genome‐wide patterns of divergence. Genome‐wide divergence between these two recently diverged sister species was notably high (genome‐wide F_ST = 0.32). Against a background of high genome‐wide divergence, regions of elevated divergence were widespread throughout the genome, as opposed to focused within islands of differentiation. These patterns may be the result of genetic drift resulting from past tidal march colonization events in conjunction with divergent selection to different environments. We identified several candidate genes that exhibited elevated divergence between saltmarsh and Nelson's sparrows, including genes linked to osmotic regulation, circadian rhythm, and plumage melanism—all putative candidates linked to adaptation to tidal marsh environments. These findings provide new insights into the roles of divergent selection and genetic drift in generating and maintaining biodiversity.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION DURING SPECIATION‐WITH‐GENE‐FLOW IN LAKE WHITEFISH SPECIES PAIRS ASSESSED BY RAD SEQUENCING

During speciation‐with‐gene‐flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome‐wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome‐wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation‐with‐gene‐flow mechanism.     

Phylogeographic and phenotypic outcomes of brown anole colonization across the Caribbean provide insight into the beginning stages of an adaptive radiation

Some of the most important insights into the ecological and evolutionary processes of diversification and speciation have come from studies of island adaptive radiations, yet relatively little research has examined how these radiations initiate. We suggest that Anolis sagrei is a candidate for understanding the origins of the Caribbean Anolis adaptive radiation and how a colonizing anole species begins to undergo allopatric diversification, phenotypic divergence and, potentially, speciation. We undertook a genomic and morphological analysis of representative populations across the entire native range of A. sagrei, finding that the species originated in the early Pliocene, with the deepest divergence occurring between western and eastern Cuba. Lineages from these two regions subsequently colonized the northern Caribbean. We find that at the broadest scale, populations colonizing areas with fewer closely related competitors tend to evolve larger body size and more lamellae on their toepads. This trend follows expectations for post‐colonization divergence from progenitors and convergence in allopatry, whereby populations freed from competition with close relatives evolve towards common morphological and ecological optima. Taken together, our results show a complex history of ancient and recent Cuban diaspora with populations on competitor‐poor islands evolving away from their ancestral Cuban populations regardless of their phylogenetic relationships, thus providing insight into the original diversification of colonist anoles at the beginning of the radiation. Our research also supplies an evolutionary framework for the many studies of this increasingly important species in ecological and evolutionary research.     

The shape of things to come in the study of the origin of species?

Perhaps Darwin would agree that speciation is no longer the mystery of mysteries that it used to be. It is now generally accepted that evolution by natural selection can contribute to ecological adaptation, resulting in the evolution of reproductive barriers and, hence, to the evolution of new species (Schluter & Conte 2009 ; Meyer 2011 ; Nosil 2012 ). From genes that encode silencing proteins that cause infertility in hybrid mice (Mihola et al. 2009 ), to segregation distorters linked to speciation in fruit flies (Phadnis & Orr 2009 ), or pollinator‐mediated selection on flower colour alleles driving reinforcement in Texan wildflowers (Hopkins & Rausher 2012 ), characterization of the genes that drive speciation is providing clues to the origin of species (Nosil & Schluter 2011 ). It is becoming apparent that, while recent work continues to overturn historical ideas about sympatric speciation (e.g. Barluenga et al. 2006 ), ecological circumstances strongly influence patterns of genomic divergence, and ultimately the establishment of reproductive isolation when gene flow is present (Elmer & Meyer 2011 ). Less clear, however, are the genetic mechanisms that cause speciation, particularly when ongoing gene flow is occurring. Now, in this issue, Franchini et al. ( 2014 ) employ a classic genetic mapping approach augmented with new genomic tools to elucidate the genomic architecture of ecologically divergent body shapes in a pair of sympatric crater lake cichlid fishes. From over 450 segregating SNPs in an F2 cross, 72 SNPs were linked to 11 QTL associated with external morphology measured by means of traditional and geometric morphometrics. Annotation of two highly supported QTL further pointed to genes that might contribute to ecological divergence in body shape in Midas cichlids, overall supporting the hypothesis that genomic regions of large phenotypic effect may be contributing to early‐stage divergence in Midas cichlids.     

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.