Highly pleiotropic genes are predicted to be used less often during adaptation, as mutations in these loci are more likely to have negative fitness consequences. Following this logic, we tested whether pleiotropy impacts the probability that a locus will be used repeatedly in adaptation. We used two proxies to estimate pleiotropy: number of phenotypic traits affected by a given genomic region and gene connectivity. We first surveyed 16 independent stream-lake and three independent benthic-limnetic ecotype pairs of threespine stickleback to estimate genome-wide patterns in parallel genomic differentiation. Our analysis revealed parallel divergence across the genome; 30%-37% of outlier regions were shared between at least two independent pairs in either the stream-lake or benthic-limnetic comparisons. We then tested whether parallel genomic regions are less pleiotropic than nonparallel regions. Counter to our a priori prediction, parallel genomic regions contained genes with significantly more pleiotropy; that is, influencing a greater number of traits and more highly connected. The increased pleiotropy of parallel regions could not be explained by other genomic factors, as there was no significant difference in mean gene count, mutation or recombination rates between parallel and nonparallel regions. Interestingly, although nonparallel regions contained genes that were less connected and influenced fewer mapped traits on average than parallel regions, they also tended to contain the genes that were predicted to be the most pleiotropic. Taken together, our findings are consistent with the idea that pleiotropy only becomes constraining at high levels and that low or intermediate levels of pleiotropy may be beneficial for adaptation.

The repeated evolution of similar phenotypes in independent populations (i.e. parallel or convergent evolution) provides an opportunity to identify genetic and ecological factors that influence the process of adaptation. Threespine stickleback fish ( Gasterosteus aculeatus) are an excellent model for such studies, as they have repeatedly adapted to divergent habitats across the Northern hemisphere. Here, we use genomic, ecological and morphological data from 16 independent pairs of stickleback populations adapted to divergent lake and stream habitats. We combine a population genomic approach to identify regions of the genome that are likely under selection in these divergent habitats with an association mapping approach to identify regions of the genome that underlie variation in ecological factors and morphological traits. Over 37% of genomic windows are repeatedly differentiated across lake-stream pairs. Similarly, many genomic windows are associated with variation in abiotic factors, diet items and morphological phenotypes. Both the highly differentiated windows and candidate trait windows are non-randomly distributed across the genome and show some overlap. However, the overlap is not significant on a genome-wide scale. Together, our data suggest that adaptation to divergent food resources and predation regimes are drivers of differentiation in lake-stream stickleback, but that additional ecological factors are also important. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Recombination often differs markedly between males and females. Here we present the first analysis of sex-specific recombination in Gasterosteus sticklebacks. Using whole-genome sequencing of 15 crosses between G. aculeatus and G. nipponicus, we localized 698 crossovers with a median resolution of 2.3 kb. We also used a bioinformatic approach to infer historical sex-averaged recombination patterns for both species. Recombination is greater in females than males on all chromosomes, and overall map length is 1.64 times longer in females. The locations of crossovers differ strikingly between sexes. Crossovers cluster toward chromosome ends in males, but are distributed more evenly across chromosomes in females. Suppression of recombination near the centromeres in males causes crossovers to cluster at the ends of long arms in acrocentric chromosomes, and greatly reduces crossing over on short arms. The effect of centromeres on recombination is much weaker in females. Genomic differentiation between G. aculeatus and G. nipponicus is strongly correlated with recombination rate, and patterns of differentiation along chromosomes are strongly influenced by male-specific telomere and centromere effects. We found no evidence for fine-scale correlations between recombination and local gene content in either sex. We discuss hypotheses for the origin of sexual dimorphism in recombination and its consequences for sexually antagonistic selection and sex chromosome evolution.

Ecological speciation with gene flow is widespread in nature [1], but it presents a conundrum: how are associations between traits under divergent natural selection and traits that contribute to assortative mating maintained? Theoretical models suggest that genetic mechanisms inhibiting free recombination between loci underlying these two types of traits (hereafter, "genetic coupling") can facilitate speciation [2-4]. Here, we perform a direct test for genetic coupling by mapping both divergent traits and female mate choice in a classic model of ecological speciation: sympatric benthic and limnetic threespine stickleback (Gasterosteus aculeatus). By measuring mate choice in F2 hybrid females, we allowed for recombination between loci underlying assortative mating and those under divergent ecological selection. In semi-natural mating arenas in which females had access to both benthic and limnetic males, we found that F2 females mated with males similar to themselves in body size and shape. In addition, we found two quantitative trait loci (QTLs) associated with female mate choice that also predicted female morphology along the benthic-limnetic trait axis. Furthermore, a polygenic genetic model that explains adaptation to contrasting benthic and limnetic feeding niches [5] also predicted F2 female mate choice. Together, these results provide empirical evidence that genetic coupling of assortative mating with traits under divergent ecological selection helps maintain species in the face of gene flow, despite a polygenic basis for adaptation to divergent environments.

Hybrid incompatibilities occur when interactions between opposite ancestry alleles at different loci reduce the fitness of hybrids. Most work on incompatibilities has focused on those that are "intrinsic," meaning they affect viability and sterility in the laboratory. Theory predicts that ecological selection can also underlie hybrid incompatibilities, but tests of this hypothesis using sequence data are scarce. In this article, we compiled genetic data for F2 hybrid crosses between divergent populations of threespine stickleback fish (Gasterosteus aculeatus L.) that were born and raised in either the field (seminatural experimental ponds) or the laboratory (aquaria). Because selection against incompatibilities results in elevated ancestry heterozygosity, we tested the prediction that ancestry heterozygosity will be higher in pond-raised fish compared to those raised in aquaria. We found that ancestry heterozygosity was elevated by approximately 3% in crosses raised in ponds compared to those raised in aquaria. Additional analyses support a phenotypic basis for incompatibility and suggest that environment-specific single-locus heterozygote advantage is not the cause of selection on ancestry heterozygosity. Our study provides evidence that, in stickleback, a coarse-albeit indirect-signal of environment-dependent hybrid incompatibility is reliably detectable and suggests that extrinsic incompatibilities can evolve before intrinsic incompatibilities.

Understanding the genetic architecture of evolutionary change remains a long-standing goal in biology. In vertebrates, skeletal evolution has contributed greatly to adaptation in body form and function in response to changing ecological variables like diet and predation. Here we use genome-wide linkage mapping in threespine stickleback fish to investigate the genetic architecture of evolved changes in many armor and trophic traits. We identify >100 quantitative trait loci (QTL) controlling the pattern of serially repeating skeletal elements, including gill rakers, teeth, branchial bones, jaws, median fin spines, and vertebrae. We use this large collection of QTL to address long-standing questions about the anatomical specificity, genetic dominance, and genomic clustering of loci controlling skeletal differences in evolving populations. We find that most QTL (76%) that influence serially repeating skeletal elements have anatomically regional effects. In addition, most QTL (71%) have at least partially additive effects, regardless of whether the QTL controls evolved loss or gain of skeletal elements. Finally, many QTL with high LOD scores cluster on chromosomes 4, 20, and 21. These results identify a modular system that can control highly specific aspects of skeletal form. Because of the general additivity and genomic clustering of major QTL, concerted changes in both protective armor and trophic traits may occur when sticklebacks inherit either marine or freshwater alleles at linked or possible "supergene" regions of the stickleback genome. Further study of these regions will help identify the molecular basis of both modular and coordinated changes in the vertebrate skeleton.

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.