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.

Populations occurring in similar habitats and displaying similar phenotypes are increasingly used to explore parallel evolution at the molecular level. This generally ignores the possibility that parallel evolution can be mimicked by the fragmentation of an ancestral population followed by genetic exchange with ecologically different populations. Here we demonstrate such an ecological vicariance scenario in multiple stream populations of threespine stickleback fish divergent from a single adjacent lake population. On the basis of demographic and population genomic analyses, we infer the initial spread of a stream-adapted ancestor followed by the emergence of a lake-adapted population, that selective sweeps have occurred mainly in the lake population, that adaptive lake–stream divergence is maintained in the face of gene flow from the lake into the streams, and that this divergence involves major inversion polymorphisms also important to marine-freshwater stickleback divergence. Overall, our study highlights the need for a robust understanding of the demographic and selective history in evolutionary investigations.

Life history divergence between populations inhabiting ecologically distinct habitats might be a potent source of reproductive isolation, but has received little attention in the context of speciation. We here test for life history divergence between threespine stickleback inhabiting Lake Constance (Central Europe) and multiple tributary streams. Otolith analysis shows that lake fish generally reproduce at two years of age, while their conspecifics in all streams have shifted to a primarily annual life cycle. This divergence is paralleled by a striking and consistent reduction in body size and fecundity in stream fish relative to lake fish. Stomach content analysis suggests that life history divergence might reflect a genetic or plastic response to pelagic versus benthic foraging modes in the lake and the streams. Microsatellite and mitochondrial markers further reveal that life history shifts in the different streams have occurred independently following the colonization by Lake Constance stickleback, and indicate the presence of strong barriers to gene flow across at least some of the lake-stream habitat transitions. Given that body size is known to strongly influence stickleback mating behavior, these barriers might well be related to life history divergence.

Genomic studies of parallel (or convergent) evolution often compare multiple populations diverged into two ecologically different habitats to search for loci repeatedly involved in adaptation. Because the shared ancestor of these populations is generally unavailable, the source of the alleles at adaptation loci, and the direction in which their frequencies were shifted during evolution, remain elusive. To shed light on these issues, we here use multiple populations of threespine stickleback fish adapted to two different types of derived freshwater habitats—basic and acidic lakes on the island of North Uist, Outer Hebrides, Scotland—and the present‐day proxy of their marine ancestor. In a first step, we combine genome‐wide pooled sequencing and targeted individual‐level sequencing to demonstrate that ecological and phenotypic parallelism in basic‐acidic divergence is reflected by genomic parallelism in dozens of genome regions. Exploiting data from the ancestor, we next show that the acidic populations, residing in ecologically more extreme derived habitats, have adapted by accumulating alleles rare in the ancestor, whereas the basic populations have retained alleles common in the ancestor. Genomic responses to selection are thus predictable from the ecological difference of each derived habitat type from the ancestral one. This asymmetric sorting of standing genetic variation at loci important to basic‐acidic divergence has further resulted in more numerous selective sweeps in the acidic populations. Finally, our data suggest that the maintenance in marine fish of standing variation important to adaptive basic‐acidic differentiation does not require extensive hybridization between the marine and freshwater populations. Overall, our study reveals striking genome‐wide determinism in both the loci involved in parallel divergence, and in the direction in which alleles at these loci have been selected.