Summary Marine stickleback fish have colonized and adapted to innumerable streams and lakes formed since the last ice age, providing an exceptional opportunity to characterize genomic mechanisms underlying repeated ecological adaptation in nature. Here we develop a high quality reference genome assembly for threespine sticklebacks. By sequencing the genomes of 20 additional individuals from a global set of marine and freshwater populations, we identify a genome-wide set of loci that are consistently associated with marine-freshwater divergence. Our results suggest that reuse of globally-shared standing genetic variation, including chromosomal inversions, plays an important role in repeated evolution of distinct marine and freshwater sticklebacks, and in the maintenance of divergent ecotypes during early stages of reproductive isolation. Both coding and regulatory changes occur in the set of loci underlying marine-freshwater evolution, with regulatory changes likely predominating in this classic example of repeated adaptive evolution in nature.

The molecular mechanisms underlying major phenotypic changes that have evolved repeatedly in nature are generally unknown. Pelvic loss in different natural populations of threespine stickleback fish has occurred by regulatory mutations deleting a tissue-specific enhancer of the Pituitary homeobox transcription factor 1 (Pitx1) gene. The high prevalence of deletion mutations at Pitx1 may be influenced by inherent structural features of the locus. Although Pitx1 null mutations are lethal in laboratory animals, Pitx1 regulatory mutations show molecular signatures of positive selection in pelvic-reduced populations. These studies illustrate how major expression and morphological changes can arise by single mutational leaps in natural populations, producing new adaptive alleles via recurrent regulatory alterations in a key developmental control gene.

Sexual antagonism, or conflict between the sexes, has been proposed as a driving force in both sex chromosome turnover and speciation. Although closely related species often have different sex chromosome systems, it is unknown whether sex chromosome turnover contributes to the evolution of reproductive isolation between species. In this study, we show that a newly evolved sex chromosome harbours genes that contribute to speciation in threespine stickleback fish (Gasterosteus aculeatus). We first identified a neo-sex chromosome system found only in one member of a sympatric species pair in Japan. We then performed genetic linkage mapping of male-specific traits important for reproductive isolation between the Japanese species pair. The neo-X chromosome harbours loci for male courtship display traits that contribute to behavioural isolation, while the ancestral X chromosome contains loci for both behavioural isolation and hybrid male sterility. Our work not only provides strong evidence for a large-X effect on reproductive isolation in a vertebrate system, but also provides direct evidence that a young neo-X chromosome contributes to reproductive isolation between closely related species. Our data suggest that sex chromosome turnover might play a greater role in speciation than previously appreciated.

Ecological differences often evolve early in speciation as divergent natural selection drives adaptation to distinct ecological niches, leading ultimately to reproductive isolation. Though this process is a major generator of biodiversity, its genetic basis remains poorly understood. Here we investigate the genetic architecture of niche differentiation in a sympatric species pair of threespine stickleback fish by mapping the environment-dependent effects of phenotypic traits on hybrid feeding and performance under semi-natural conditions. We show that multiple, unlinked loci act largely additively to determine position along the major niche axis separating these recently diverged species. We also find that functional mismatch between phenotypic traits reduces growth of some stickleback hybrids beyond that expected from an intermediate phenotype, suggesting a role for epistasis between the underlying genes. This functional mismatch might lead to hybrid incompatibilities that are analogous to those underlying intrinsic reproductive isolation but that depend on the ecological context.

Summary Genes underlying repeated adaptive evolution in natural populations are still largely unknown. Stickleback fish (Gasterosteus aculeatus) have undergone a recent dramatic evolutionary radiation, generating numerous examples of marine–freshwater species-pairs, and a small number of benthic–limnetic species-pairs found within single lakes [1]. We have developed a new genome-wide SNP genotyping array to study patterns of genetic variation in sticklebacks over a wide geographic range, and to scan the genome for regions that contribute to repeated evolution of marine–freshwater or benthic–limnetic species-pairs. Surveying 34 global populations with 1159 informative markers revealed substantial genetic variation, with predominant patterns reflecting demographic history and geographic structure. After correcting for geographic structure and filtering for neutral markers, we detect large repeated shifts in allele frequency at some loci, identifying both known and novel loci likely contributing to marine–freshwater and benthic–limnetic divergence. Several novel loci fall close to genes implicated in epithelial barrier or immune functions, which have likely changed as sticklebacks adapt to contrasting environments. Specific alleles differentiating sympatric benthic–limnetic species-pairs are shared in nearby solitary populations, suggesting an allopatric origin for adaptive variants, and selection pressures unrelated to sympatry in the initial formation of these classic vertebrate species-pairs.

Understanding the genetics of adaptation is a central focus in evolutionary biology. Here, we use a population genomics approach to examine striking parallel morphological divergences of parapatric stream -lake ecotypes of threespine stickleback fish in three watersheds on the Haida Gwaii archipelago, western Canada. Genome-wide variation at greater than 1000 single nucleotide polymorphism loci indicate separate origin of giant lake and small-bodied stream fish within each watershed (mean F ST between watersheds ¼ 0.244 and within ¼ 0.114). Genome scans within watersheds identified a total of 21 genomic regions that are highly differentiated between ecotypes and are probably subject to directional selection. Most outliers were watershed-specific, but genomic regions undergoing parallel genetic changes in multiple watersheds were also identified. Interestingly, several of the stream -lake outlier regions match those previously identified in marine -freshwater and benthic -limnetic genome scans, indicating reuse of the same genetic loci in different adaptive scenarios. We also identified multiple new outlier loci, which may contribute to unique aspects of differentiation in stream-lake environments. Overall, our data emphasize the important role of ecological boundaries in driving both local and broadly occurring parallel genetic changes during adaptation.

How predictable is the genetic basis of phenotypic adaptation? Answering this question begins by estimating the repeatability of adaptation at the genetic level. Here, we provide a comprehensive estimate of the repeatability of the genetic basis of adaptive phenotypic evolution in a natural system. We used quantitative trait locus (QTL) mapping to discover genomic regions controlling a large number of morphological traits that have diverged in parallel between pairs of threespine stickleback (Gasterosteus aculeatus species complex) in Paxton and Priest lakes, British Columbia. We found that nearly half of QTL affected the same traits in the same direction in both species pairs. Another 40% influenced a parallel phenotypic trait in one lake but not the other. The remaining 10% of QTL had phenotypic effects in opposite directions in the two species pairs. Similarity in the proportional contributions of all QTL to parallel trait differences was about 0.4. Surprisingly, QTL reuse was unrelated to phenotypic effect size. Our results indicate that repeated use of the same genomic regions is a pervasive feature of parallel phenotypic adaptation, at least in sticklebacks. Identifying the causes of this pattern would aid prediction of the genetic basis of phenotypic evolution.KEYWORDS genetics of adaptation; genetic parallelism; parallel evolution; repeated evolution; QTL mapping W E still have a poor understanding of the predictability of the genetic basis of phenotypic adaptation (Stern and Orgogozo 2008;Conte et al. 2012;Martin and Orgogozo 2013;Stern 2013). One way to make progress is to quantify the repeatability of the genetic changes that underlie repeated phenotypic evolution. When organisms independently evolve similar phenotypes in response to similar selection pressures (a reliable signature of adaptive evolution (Endler 1986;Harvey and Pagel 1991;Schluter 2000; Losos 2011), we can ask: How similar are the genetic "solutions" underlying those phenotypes? Many features may influence similarity of genetic solutions, including availability of standing genetic variation, mutational biases, and functional constraints. The extent to which the genetic basis of repeatedly evolved phenotypes is shared indicates the extent to which the genetic basis of adaptation is predictable.Although there is some evidence indicating that gene reuse during repeated phenotypic evolution is common, current estimates of the frequency of gene reuse in adaptive evolution might not be accurate because of methodological limitations. A recent metaanalysis of studies of natural populations estimated that the average probability of gene reuse is 0.32-0.55 (depending on the type of data used to calculate it) across a diversity of taxa spanning divergence times from hundreds of years to hundreds of millions of years (Conte et al. 2012). The probability of gene reuse was highest among closely related species and it declined with increasing divergence time between the taxa being compared. This estimate of gene reuse was based on published...

Evolutionary studies are often limited by missing data that are critical to understanding the history of selection. Selection experiments, which reproduce rapid evolution under controlled conditions, are excellent tools to study how genomes evolve under selection. Here we present a genomic dissection of the Longshanks selection experiment, in which mice were selectively bred over 20 generations for longer tibiae relative to body mass, resulting in 13% longer tibiae in two replicates. We synthesized evolutionary theory, genome sequences and molecular genetics to understand the selection response and found that it involved both polygenic adaptation and discrete loci of major effect, with the strongest loci tending to be selected in parallel between replicates. We show that selection may favor de-repression of bone growth through inactivating two limb enhancers of an inhibitor, Nkx3-2. Our integrative genomic analyses thus show that it is possible to connect individual base-pair changes to the overall selection response.