Assembly of the threespine stickleback Y chromosome reveals convergent signatures of sex chromosome evolution.

BACKGROUND: Heteromorphic sex chromosomes have evolved repeatedly across diverse species. Suppression of recombination between X and Y chromosomes leads to degeneration of the Y chromosome. The progression of degeneration is not well understood, as complete sequence assemblies of heteromorphic Y chromosomes have only been generated across a handful of taxa with highly degenerate sex chromosomes. Here, we describe the assembly of the threespine stickleback (Gasterosteus aculeatus) Y chromosome, which is less than 26 million years old and at an intermediate stage of degeneration. Our previous work identified that the non-recombining region between the X and the Y spans approximately 17.5 Mb on the X chromosome. RESULTS: We combine long-read sequencing with a Hi-C-based proximity guided assembly to generate a 15.87 Mb assembly of the Y chromosome. Our assembly is concordant with cytogenetic maps and Sanger sequences of over 90 Y chromosome BAC clones. We find three evolutionary strata on the Y chromosome, consistent with the three inversions identified by our previous cytogenetic analyses. The threespine stickleback Y shows convergence with more degenerate sex chromosomes in the retention of haploinsufficient genes and the accumulation of genes with testis-biased expression, many of which are recent duplicates. However, we find no evidence for large amplicons identified in other sex chromosome systems. We also report an excellent candidate for the master sex-determination gene: a translocated copy of Amh (Amhy). CONCLUSIONS: Together, our work shows that the evolutionary forces shaping sex chromosomes can cause relatively rapid changes in the overall genetic architecture of Y chromosomes.

Genetic Coupling of Female Mate Choice with Polygenic Ecological Divergence Facilitates Stickleback Speciation.

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.

Genetic mapping of natural variation in schooling tendency in the threespine stickleback.

Although there is a heritable basis for many animal behaviors, the genetic architecture of behavioral variation in natural populations remains mostly unknown, particularly in vertebrates. We sought to identify the genetic basis for social affiliation in two populations of threespine sticklebacks (Gasterosteus aculeatus) that differ in their propensity to school. Marine sticklebacks from Japan school strongly whereas benthic sticklebacks from a lake in Canada are more solitary. Here, we expanded on our previous efforts to identify quantitative trait loci (QTL) for differences in schooling tendency. We tested fish multiple times in two assays that test different aspects of schooling tendency: 1) the model school assay, which presents fish with a school of eight model sticklebacks; and 2) the choice assay, in which fish are given a choice between the model school and a stationary artificial plant. We found low-to-moderate levels of repeatability, ranging from 0.1 to 0.5, in schooling phenotypes. To identify the genomic regions that contribute to differences in schooling tendency, we used QTL mapping in two types of crosses: benthic × marine backcrosses and an F2 intercross. We found two QTL for time spent with the school in the model school assay, and one QTL for number of approaches to the school in the choice assay. These QTL were on three different linkage groups, not previously linked to behavioral differences in sticklebacks. Our results highlight the importance of using multiple crosses and robust behavioral assays to uncover the genetic basis of behavioral variation in natural populations.