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 nonrecombining 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 testisbiased 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.

No abstract

The stickleback family (Gasterosteidae) of fish is less than 40 million years old, yet stickleback species have diverged in both diploid chromosome number (2n) and morphology. We used comparative fluorescence in situ hybridization (FISH) on 2 stickleback species, Gasterosteus aculeatus (2n = 42) and Apeltes quadracus (2n = 46), to ascertain the types of chromosome rearrangements that differentiate these species. The A. quadracus karyotype contains more acrocentric and telocentric chromosomes than the G. aculeatus karyotype. By using bacterial artificial chromosome probes from G. aculeatus in our FISH screen, we found that 6 pericentric inversions and 2 chromosome fusions/fissions are responsible for the greater number of acrocentric and telocentric chromosomes in A. quadracus. While most populations of G. aculeatus have an XX/XY sex chromosome system, A. quadracus has a ZZ/ZW sex chromosome system, as previously reported. However, we discovered that a population of A. quadracus from Connecticut lacks heteromorphic sex chromosomes, providing evidence for unexpected sex chromosome diversity in this species.

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.