Adaptive divergence between lake and stream populations of an East African cichlid fish

Divergent natural selection acting in different habitats may build up barriers to gene flow and initiate speciation. This speciation continuum can range from weak or no divergence to strong genetic differentiation between populations. Here, we focus on the early phases of adaptive divergence in the East African cichlid fish Astatotilapia burtoni, which occurs in both Lake Tanganyika (LT) and inflowing rivers. We first assessed the population structure and morphological differences in A. burtoni from southern LT. We then focused on four lake–stream systems and quantified body shape, ecologically relevant traits (gill raker and lower pharyngeal jaw) as well as stomach contents. Our study revealed the presence of several divergent lake–stream populations that rest at different stages of the speciation continuum, but show the same morphological and ecological trajectories along the lake–stream gradient. Lake fish have higher bodies, a more superior mouth position, longer gill rakers and more slender pharyngeal jaws, and they show a plant/algae and zooplankton‐biased diet, whereas stream fish feed more on snails, insects and plant seeds. A test for reproductive isolation between closely related lake and stream populations did not detect population‐assortative mating. Analyses of F1 offspring reared under common garden conditions indicate that the detected differences in body shape and gill raker length do not constitute pure plastic responses to different environmental conditions, but also have a genetic basis. Taken together, the A. burtoni lake–stream system constitutes a new model to study the factors that enhance and constrain progress towards speciation in cichlid fishes.     

VARIABLE PROGRESS TOWARD ECOLOGICAL SPECIATION IN PARAPATRY: STICKLEBACK ACROSS EIGHT LAKE-STREAM TRANSITIONS

Divergent selection between contrasting habitats can sometimes drive adaptive divergence and the evolution of reproductive isolation in the face of initially high gene flow. "Progress" along this ecological speciation pathway can range from minimal divergence to full speciation. We examine this variation for threespine stickleback fish that evolved independently across eight lake-stream habitat transitions. By quantifying stickleback diets, we show that lake-stream transitions usually coincide with limnetic-benthic ecotones. By measuring genetically based phenotypes, we show that these ecotones often generate adaptive divergence in foraging morphology. By analyzing neutral genetic markers (microsatellites), we show that adaptive divergence is often associated with the presence of two populations maintaining at least partial reproductive isolation in parapatry. Coalescent-based simulations further suggest that these populations have diverged with gene flow within a few thousand generations, although we cannot rule out the possibility of phases of allopatric divergence. Finally, we find striking variation among the eight lake-stream transitions in progress toward ecological speciation. This variation allows us to hypothesize that progress is generally promoted by strong divergent selection and limited dispersal across the habitat transitions. Our study thus makes a case for ecological speciation in a parapatric context, while also highlighting variation in the outcome.     

Diversifying selection drives parallel evolution of gill raker number and body size along the speciation continuum of European whitefish

Adaptive radiation is the evolution of ecological and phenotypical diversity. It arises via ecological opportunity that promotes the exploration of underutilized or novel niches mediating specialization and reproductive isolation. The assumed precondition for rapid local adaptation is diversifying natural selection, but random genetic drift could also be a major driver of this process. We used 27 populations of European whitefish (Coregonus lavaretus) from nine lakes distributed in three neighboring subarctic watercourses in northern Fennoscandia as a model to test the importance of random drift versus diversifying natural selection for parallel evolution of adaptive phenotypic traits. We contrasted variation for two key adaptive phenotypic traits correlated with resource utilization of polymorphic fish; the number of gill rakers and the total length of fish, with the posterior distribution of neutral genetic differentiation from 13 microsatellite loci, to test whether the observed phenotypic divergence could be achieved by random genetic drift alone. Our results show that both traits have been under diversifying selection and that the evolution of these morphs has been driven by isolation through habitat adaptations. We conclude that diversifying selection acting on gill raker number and body size has played a significant role in the ongoing adaptive radiation of European whitefish morphs in this region.     

Divergent Macroparasite Infections in Parapatric Swiss Lake-Stream Pairs of Threespine Stickleback (Gasterosteus aculeatus)

Spatial heterogeneity in diversity and intensity of parasitism is a typical feature of most host-parasite interactions, but understanding of the evolutionary implications of such variation is limited. One possible outcome of infection heterogeneities is parasite-mediated divergent selection between host populations, ecotypes or species which may facilitate the process of ecological speciation. However, very few studies have described infections in population-pairs along the speciation continuum from low to moderate or high degree of genetic differentiation that would address the possibility of parasite-mediated divergent selection in the early stages of the speciation process. Here we provide an example of divergent parasitism in freshwater fish ecotypes by examining macroparasite infections in threespine stickleback (Gasterosteus aculeatus) of four Swiss lake systems each harbouring parapatric lake-stream ecotype pairs. We demonstrate significant differences in infections within and between the pairs that are driven particularly by the parasite taxa transmitted to fish from benthic invertebrates. The magnitude of the differences tended to correlate positively with the extent of neutral genetic differentiation between the parapatric lake and stream populations of stickleback, whereas no such correlation was found among allopatric populations from similar or contrasting habitats. This suggests that genetic differentiation is unrelated to the magnitude of parasite infection contrasts when gene flow is constrained by geographical barriers while in the absence of physical barriers, genetic differentiation and the magnitude of differences in infections tend to be positively correlated.     

Genome-wide patterns of divergence during speciation: the lake whitefish case study

The nature, size and distribution of the genomic regions underlying divergence and promoting reproductive isolation remain largely unknown. Here, we summarize ongoing efforts using young (12 000 yr BP) species pairs of lake whitefish (Coregonus clupeaformis) to expand our understanding of the initial genomic patterns of divergence observed during speciation. Our results confirmed the predictions that: (i) on average, phenotypic quantitative trait loci (pQTL) show higher F(ST) values and are more likely to be outliers (and therefore candidates for being targets of divergent selection) than non-pQTL markers; (ii) large islands of divergence rather than small independent regions under selection characterize the early stages of adaptive divergence of lake whitefish; and (iii) there is a general trend towards an increase in terms of numbers and size of genomic regions of divergence from the least (East L.) to the most differentiated species pair (Cliff L.). This is consistent with previous estimates of reproductive isolation between these species pairs being driven by the same selective forces responsible for environment specialization. Altogether, dwarf and normal whitefish species pairs represent a continuum of both morphological and genomic differentiation contributing to ecological speciation. Admittedly, much progress is still required to more finely map and circumscribe genomic islands of speciation. This will be achieved through the use of next generation sequencing data but also through a better quantification of phenotypic traits moulded by selection as organisms adapt to new environmental conditions.     

Spatial phenotypic and genetic structure of threespine stickleback (Gasterosteus aculeatus) in a heterogeneous natural system,Lake Mývatn,Iceland

Eco‐evolutionary responses of natural populations to spatial environmental variation strongly depend on the relative strength of environmental differences/natural selection and dispersal/gene flow. In absence of geographic barriers, as often is the case in lake ecosystems, gene flow is expected to constrain adaptive divergence between environments – favoring phenotypic plasticity or high trait variability. However, if divergent natural selection is sufficiently strong, adaptive divergence can occur in face of gene flow. The extent of divergence is most often studied between two contrasting environments, whereas potential for multimodal divergence is little explored. We investigated phenotypic (body size, defensive structures, and feeding morphology) and genetic (microsatellites) structure in threespine stickleback (Gasterosteus aculeatus) across five habitat types and two basins (North and South) within the geologically young and highly heterogeneous Lake Mývatn, North East Iceland. We found that (1) North basin stickleback were, on average, larger and had relatively longer spines than South basin stickleback, whereas (2) feeding morphology (gill raker number and gill raker gap width) differed among three of five habitat types, and (3) there was only subtle genetic differentiation across the lake. Overall, our results indicate predator and prey mediated phenotypic divergence across multiple habitats in the lake, in face of gene flow.     

Along the speciation continuum: Quantifying intrinsic and extrinsic isolating barriers across five million years of evolutionary divergence in California jewelflowers

Understanding the relative roles of intrinsic and extrinsic reproductive barriers, and their interplay within the geographic context of diverging taxa, remains an outstanding challenge in the study of speciation. We conducted a comparative analysis of reproductive isolation in California Jewelflowers (Streptanthus, s.l., Brassicaceae) by quantifying potential barriers to gene flow at multiple life history stages in 39 species pairs spanning five million years of evolutionary divergence. We quantified nine potential pre‐ and postzygotic barriers and explored patterns of reproductive isolation in relation to genetic distance. Intrinsic postzygotic isolation was initially weak, increased at intermediate genetic distances, and reached a threshold characterized by complete genetic incompatibility. Climatic niche differences were strong at shallow genetic distances, and species pairs with overlapping ranges showed slight but appreciable phenological isolation, highlighting the potential for ecological barriers to contribute to speciation. Geographic analyses suggest that speciation is not regionally allopatric in the California Jewelflowers, as recently diverged taxa occur in relatively close proximity and display substantial range overlap. Young pairs are characterized by incomplete intrinsic postzygotic isolation, suggesting that extrinsic barriers or fine‐scale spatial segregation are more important early in the divergence process than genetic incompatibilities.     

Phenotypic divergence and reproductive isolation between sympatric forms of Japanese threespine sticklebacks

The threespine stickleback ( Gasterosteus aculeatus ) species complex is well suited for identifying the types of phenotypic divergence and isolating barriers that contribute to reproductive isolation at early stages of speciation. In the present study, we characterize the patterns of genetic and phenotypic divergence as well as the types of isolating barriers that are present between two sympatric pairs of threespine sticklebacks in Hokkaido, Japan. One sympatric pair consists of an anadromous and a resident freshwater form and shows divergence in body size between the forms, despite the lack of genetic differentiation between them. The second sympatric pair consists of two anadromous forms, which originated before the last glacial period and are currently reproductively isolated. These two anadromous forms have diverged in many morphological traits as well as in their reproductive behaviours. Both sexual isolation and hybrid male sterility contribute to reproductive isolation between the anadromous species pair. We discuss the shared and unique aspects of phenotypic divergence and reproductive isolation in the Japanese sympatric pairs compared with postglacial stickleback species pairs. Further studies of these divergent species pairs will provide a deeper understanding of the mechanisms of speciation in sticklebacks.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 671–685.     

Coexistence and origin of trophic ecotypes of pygmy whitefish,Prosopium coulterii,in a south‐western Alaskan lake

Ecologically, morphologically and genetically distinct populations within single taxa often coexist in postglacial lakes and have provided important model systems with which to investigate ecological and evolutionary processes such as niche partitioning and ecological speciation. Within the Salmonidae, these species complexes have been well studied, particularly within the Coregonus clupeaformis–C. laveratus (lake and European whitefish, respectively) group, but the phenomenon has been less well documented in the other whitefish genera, Prosopium and Stenodus. Here, we examined the morphology, feeding biology and genetic structure of three putative forms of the pygmy whitefish, Prosopium coulterii (Eigenmann & Eigenmann, 1892), first reported from Chignik Lake, south‐western Alaska, over 40 years ago. Field collections and morphological analyses resolved a shallow water (< 5 m depth) low gill raker count form (< 15 first arch gill rakers), a deepwater (> 30 m), low gill raker form and a deepwater, high gill raker count (> 15 gill rakers) form. The two low gill raker count forms fed almost exclusively on benthic invertebrates (mostly chironomids), while the deepwater, high gill raker count form fed almost exclusively on zooplankton; differences in diet were also reflected in differences both in δ¹³C and δ¹⁵N stable isotopes. All three forms were characterized by the same major mitochondrial DNA clade that has been associated with persistence in, and postglacial dispersal from, a Beringian glacial refugium. Analysis of variation at nine microsatellite DNA loci indicated low, but significant differentiation among forms, especially between the two low gill raker count forms and the high gill raker count form. The extent of differentiation along phenotypic (considerable) and genetic (subtle) axes among the Chignik Lake forms is similar to that found among distinct taxa of Prosopium found in pre‐glacial Bear Lake (Utah–Idaho, USA) which is probably at least ten times older than Chignik Lake. Our analyses illustrate the potential for the postglacial differentiation in traits subject to divergent natural selection across variable environments.     

Selection on outlier loci and their association with adaptive phenotypes in Littorina saxatilis contact zones

A fundamental issue in speciation research is to evaluate phenotypic variation and the genomics driving the evolution of reproductive isolation between sister taxa. Above all, hybrid zones are excellent study systems for researchers to examine the association of genetic differentiation, phenotypic variation and the strength of selection. We investigated two contact zones in the marine gastropod Littorina saxatilis and utilized landmark‐based geometric morphometric analysis together with amplified fragment length polymorphism (AFLP) markers to assess phenotypic and genomic divergence between ecotypes under divergent selection. From genetic markers, we calculated the cline width, linkage disequilibrium and the average effective selection on a locus. Additionally, we conducted an association analysis linking the outlier loci and phenotypic variation between ecotypes and show that a proportion of outlier loci are associated with key adaptive phenotypic traits.     

Generation of a neutral FST baseline for testing local adaptation on gill raker number within and between European whitefish ecotypes in the Baltic Sea basin

Divergent selection at ecologically important traits is thought to be a major factor driving phenotypic differentiation between populations. To elucidate the role of different evolutionary processes shaping the variation in gill raker number of European whitefish (Coregonus lavaretus sensu lato) in the Baltic Sea basin, we assessed the relationships between genetic and phenotypic variation among and within three whitefish ecotypes (sea spawners, river spawners and lake spawners). To generate expected neutral distribution of F_ST and to evaluate whether highly variable microsatellite loci resulted in deflated F_ST estimates compared to less variable markers, we performed population genetic simulations under finite island and hierarchical island models. The genetic divergence observed among (F_CT = 0.010) and within (F_ST = 0.014–0.041) ecotypes was rather low. The divergence in gill raker number, however, was substantially higher between sea and river spawners compared to observed microsatellite data and simulated neutral baseline (P_CT > F_CT). This suggests that the differences in gill raker number between sea and river spawners are likely driven by divergent natural selection. We also found strong support for divergent selection on gill raker number among different populations of sea spawners (P_ST > F_ST), most likely caused by highly variable habitat use and diverse diet. The putative role of divergent selection within lake spawners initially inferred from empirical microsatellite data was not supported by simulated F_ST distributions. This work provides a first formal test of divergent selection on gill raker number in Baltic whitefish, and demonstrates the usefulness of population genetic simulations to generate informative neutral baselines for P_ST–F_ST analyses helping to disentangle the effects of stochastic evolutionary processes from natural selection.     

Heterogeneous genomic differentiation between walking-stick ecotypes: "isolation by adaptation" and multiple roles for divergent selection

Genetic differentiation can be highly variable across the genome. For example, loci under divergent selection and those tightly linked to them may exhibit elevated differentiation compared to neutral regions. These represent "outlier loci" whose differentiation exceeds neutral expectations. Adaptive divergence can also increase genome-wide differentiation by promoting general barriers to neutral gene flow, thereby facilitating genomic divergence via genetic drift. This latter process can yield a positive correlation between adaptive phenotypic divergence and neutral genetic differentiation (described here as "isolation-by-adaptation"). Here, we examine both these processes by combining an AFLP genome scan of two host plant ecotypes of Timema cristinae walking-sticks with existing data on adaptive phenotypic divergence and ecological speciation in these insects. We found that about 8% of loci are outliers in multiple population comparisons. Replicated comparisons between population-pairs using the same versus different host species revealed that 1-2% of loci are subject to host-related selection specifically. Locus-specific analyses revealed that up to 10% of putatively neutral (nonoutlier) AFLP loci exhibit significant isolation-by-adaptation. Our results suggest that selection may affect differentiation directly, via linkage, or by facilitating genetic drift. They thus illustrate the varied and sometimes nonintuitive contributions of selection to heterogeneous genomic differentiation.     

Mate preference across the speciation continuum in a clade of mimetic butterflies

Premating behavioral isolation is increasingly recognized as an important part of ecological speciation, where divergent natural selection causes the evolution of reproductive barriers. A number of studies have now demonstrated that traits under divergent natural selection also affect mate preferences. However, studies of single species pairs only capture a snapshot of the speciation process, making it difficult to assess the role of mate preferences throughout the entire process. Heliconius butterflies are well known for their brightly colored mimetic warning patterns, and previous studies have shown that these patterns are also used as mate recognition cues. Here, we present mate preference data for four pairs of sister taxa, representing different stages of divergence, which together allow us to compare diverging mate preferences across the continuum of Heliconius speciation. Using a novel Bayesian approach, our results support a model of ecological speciation in which strong premating isolation arises early, but continues to increase throughout the continuum from polymorphic populations through to "good," sympatric ecologically divergent species.     

Molecular differentiation at nuclear loci in French host races of the European corn borer (Ostrinia nubilalis)

French populations of the European corn borer consist of two sympatric and genetically differentiated host races. As such, they are well suited to study processes that could be involved in sympatric speciation, but the initial conditions of host-race divergence need to be elucidated. Gene genealogies can provide insight into the processes involved in speciation. We used DNA sequences of four nuclear genes to (1) document the genetic structure of the two French host races previously delineated with allozyme markers, (2) find genes directly or indirectly involved in reproductive isolation between host races, and (3) estimate the time since divergence of the two taxa and see whether this estimate is compatible with this divergence being the result of a host shift onto maize after its introduction into Europe approximately 500 years ago. Gene genealogies revealed extensive shared polymorphism, but confirmed the previously observed genetic differentiation between the two host races. Significant departures from the predictions of neutral molecular evolution models were detected at three loci but were apparently unrelated to reproductive isolation between host races. Estimates of time since divergence between French host races varied from approximately 75,000 to approximately 150,000 years, suggesting that the two taxa diverged recently but probably long before the introduction of maize into Europe.     

Ecological speciation in postglacial European whitefish: rapid adaptive radiations into the littoral,pelagic, and profundal lake habitats

Understanding how a monophyletic lineage of a species diverges into several adaptive forms has received increased attention in recent years, but the underlying mechanisms in this process are still under debate. Postglacial fishes are excellent model organisms for exploring this process, especially the initial stages of ecological speciation, as postglacial lakes represent replicated discrete environments with variation in available niches. Here, we combine data of niche utilization, trophic morphology, and 17 microsatellite loci to investigate the diversification process of three sympatric European whitefish morphs from three northern Fennoscandian lakes. The morphological divergence in the gill raker number among the whitefish morphs was related to the utilization of different trophic niches and was associated with reproductive isolation within and across lakes. The intralacustrine comparison of whitefish morphs showed that these systems represent two levels of adaptive divergence: (1) a consistent littoral–pelagic resource axis; and (2) a more variable littoral–profundal resource axis. The results also indicate that the profundal whitefish morph has diverged repeatedly from the ancestral littoral whitefish morph in sympatry in two different watercourses. In contrast, all the analyses performed revealed clustering of the pelagic whitefish morphs across lakes suggesting parallel postglacial immigration with the littoral whitefish morph into each lake. Finally, the analyses strongly suggested that the trophic adaptive trait, number of gill rakers, was under diversifying selection in the different whitefish morphs. Together, the results support a complex evolutionary scenario where ecological speciation acts, but where both allopatric (colonization history) and sympatric (within watercourse divergence) processes are involved.     

Effects of assortative mate choice on the genomic and morphological structure of a hybrid zone between two bird subspecies

Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between the alba and personata subspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437 SNP loci on a 3,000 km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome‐wide divergence. Variation in only one trait—head plumage patterning—was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100 km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome‐wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.     

Discrete foraging niches promote ecological, phenotypic, and genetic divergence in sympatric whitefish (Coregonus lavaretus)

Natural populations often vary in their degree of ecological, morphological and genetic divergence. This variation can be arranged along an ecological speciation continuum of increasingly discrete variation, with high inter-individual variation at one end and well defined species in the other. In postglacial fishes, evolutionary divergence has commonly resulted in the co-occurrence of a pelagic and a benthic specialist. We studied three replicate lakes supporting sympatric pelagic and benthic European whitefish (Coregonus lavaretus (L.)) morphs in search for early signs of possible further divergence into more specialized niches. Using stomach content data (recent diet) and stable isotope analyses (time-integrated measure of trophic niche use), we observed a split in the trophic niche within the benthic whitefish morph, with individuals specializing on either littoral or profundal resources. This divergence in resource use was accompanied by small but significant differences in an adaptive morphological trait (gill raker number) and significant genetic differences between fish exploiting littoral and profundal habitats and foraging resources. The same pattern of parallel divergence was found in all three lakes, suggesting similar natural selection pressures driving and/or maintaining the divergence. The two levels of divergence (a clear and robust benthic – pelagic and a more subtle littoral – profundal divergence) observed in this study apparently represent different stages in the process of ecological speciation.     

Quantitative genetic inheritance of morphological divergence in a lake-stream stickleback ecotype pair: implications for reproductive isolation

Ecological selection against hybrids between populations occupying different habitats might be an important component of reproductive isolation during the initial stages of speciation. The strength and directionality of this barrier to gene flow depends on the genetic architecture underlying divergence in ecologically relevant phenotypes. We here present line cross analyses of inheritance for two key foraging-related morphological traits involved in adaptive divergence between stickleback ecotypes residing parapatrically in lake and stream habitats within the Misty Lake watershed (Vancouver Island, Canada). One main finding is the striking genetic dominance of the lake phenotype for body depth. Selection associated with this phenotype against first- and later-generation hybrids should therefore be asymmetric, hindering introgression from the lake to the stream population but not vice versa. Another main finding is that divergence in gill raker number is inherited additively and should therefore contribute symmetrically to reproductive isolation. Our study suggests that traits involved in adaptation might contribute to reproductive isolation qualitatively differently, depending on their mode of inheritance.     

Divergent selection and heterogeneous genomic divergence

Levels of genetic differentiation between populations can be highly variable across the genome, with divergent selection contributing to such heterogeneous genomic divergence. For example, loci under divergent selection and those tightly physically linked to them may exhibit stronger differentiation than neutral regions with weak or no linkage to such loci. Divergent selection can also increase genome‐wide neutral differentiation by reducing gene flow (e.g. by causing ecological speciation), thus promoting divergence via the stochastic effects of genetic drift. These consequences of divergent selection are being reported in recently accumulating studies that identify: (i) ‘outlier loci’ with higher levels of divergence than expected under neutrality, and (ii) a positive association between the degree of adaptive phenotypic divergence and levels of molecular genetic differentiation across population pairs [‘isolation by adaptation’ (IBA)]. The latter pattern arises because as adaptive divergence increases, gene flow is reduced (thereby promoting drift) and genetic hitchhiking increased. Here, we review and integrate these previously disconnected concepts and literatures. We find that studies generally report 5–10% of loci to be outliers. These selected regions were often dispersed across the genome, commonly exhibited replicated divergence across different population pairs, and could sometimes be associated with specific ecological variables. IBA was not infrequently observed, even at neutral loci putatively unlinked to those under divergent selection. Overall, we conclude that divergent selection makes diverse contributions to heterogeneous genomic divergence. Nonetheless, the number, size, and distribution of genomic regions affected by selection varied substantially among studies, leading us to discuss the potential role of divergent selection in the growth of regions of differentiation (i.e. genomic islands of divergence), a topic in need of future investigation.     

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

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