Evolutionary diversification is often initiated by adaptive divergence between populations occupying ecologically distinct environments while still exchanging genes. The genetic foundations of this divergence process are largely unknown and are here explored through genome scans in multiple independent lake-stream population pairs of threespine stickleback. We find that across the pairs, overall genomic divergence is associated with the magnitude of divergence in phenotypes known to be under divergent selection. Along this same axis of increasing diversification, genomic divergence becomes increasingly biased towards the centre of chromosomes as opposed to the peripheries. We explain this pattern by within-chromosome variation in the physical extent of hitchhiking, as recombination is greatly reduced in chromosome centres. Correcting for this effect suggests that a great number of genes distributed widely across the genome are involved in the divergence into lake vs. stream habitats. Analyzing additional allopatric population pairs, however, reveals that strong divergence in some genomic regions has been driven by selection unrelated to lake-stream ecology. Our study highlights a major contribution of large-scale variation in recombination rate to generating heterogeneous genomic divergence and indicates that elucidating the genetic basis of adaptive divergence might be more challenging than currently recognized.

Heterogeneity in recombination rate may strongly influence genome evolution and entail methodological challenges to genomic investigation. Nevertheless, a solid understanding of these issues awaits detailed information across a broad range of taxa. Based on 282 F(2) individuals and 1872 single nucleotide polymorphisms, we characterize recombination in the threespine stickleback fish genome. We find an average genome-wide recombination rate of 3.11 cm/Mb. Crossover frequencies are dramatically elevated in the chromosome peripheries as compared to the centres, and are consistent with one obligate crossover per chromosome (but not chromosome arm). Along the sex chromosome, we show that recombination is restricted to a small pseudoautosomal domain of c. 2 Mb, spanning only 10% of that chromosome. Comparing female to male RAD sequence coverage allows us to identify two discrete levels of degeneration on the Y chromosome, one of these 'evolutionary strata' coinciding with a previously inferred inverted region. Using polymorphism data from two young (<10 000 years old) ecologically diverged lake-stream population pairs, we demonstrate that recombination rate correlates with both the magnitude of allele frequency shifts between populations and levels of genetic diversity within populations. These associations reflect genome-wide heterogeneity in the influence of selection on linked sites. We further find a strong relationship between recombination rate and GC content, possibly driven by GC-biased gene conversion. Overall, we highlight that heterogeneity in recombination rate has profound consequences on genome evolution and deserves wider recognition in marker-based genomic analyses.

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. Coalescentbased 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.

To what extent are patterns of biological diversification determined by natural selection? We addressed this question by exploring divergence in foraging morphology of threespine stickleback fish inhabiting lake and stream habitats within eight independent watersheds. We found that lake fish generally displayed more developed gill structures and had more streamlined bodies than did stream fish. Diet analysis revealed that these morphological differences were associated with limnetic vs. benthic foraging modes, and that the extent of morphological divergence within watersheds reflected differences in prey resources utilized by lake and stream fish. We also found that patterns of divergence were unrelated to patterns of phenotypic trait (co)variance within populations (i.e. the ‘line of least resistance’). Instead, phenotypic (co)variances were more likely to have been shaped by adaptation to lake vs. stream habitats. Our study thus implicates natural selection as a strong deterministic force driving morphological diversification in lake–stream stickleback. The strength of this inference was obtained by complementing a standard analysis of parallel divergence in means between discrete habitat categories (lake vs. stream) with quantitative estimates of selective forces and information on trait (co)variances.

Speciation can be viewed as a continuum, potentially divisible into several states: (1) continuous variation within panmictic populations, (2) partially discontinuous variation with minor reproductive isolation, (3) strongly discontinuous variation with strong but reversible reproductive isolation and (4) complete and irreversible reproductive isolation. Research on sticklebacks (Gasterosteidae) reveals factors that influence progress back and forth along this continuum, as well as transitions between the states. Most populations exist in state 1, even though some of these show evidence of disruptive selection and positive assortative mating. Transitions to state 2 seem to usually involve strong divergent selection coupled with at least a bit of geographic separation, such as parapatry (e.g. lake and stream pairs and mud and lava pairs) or allopatry (e.g. different lakes). Transitions to state 3 can occur when allopatric or parapatric populations that evolved under strong divergent selection come into secondary contact (most obviously the sympatric benthic and limnetic pairs), but might also occur between populations that remained in parapatry or allopatry. Transitions to state 4 might be decoupled from these selective processes, because the known situations of complete, or nearly complete, reproductive isolation (Japan Sea and Pacific Ocean pair and the recognized gasterosteid species) are always associated with chromosomal rearrangements and environment-independent genetic incompatibilities. Research on sticklebacks has thus revealed complex and shifting interactions between selection, adaptation, mutation and geography during the course of speciation.