The 22,000-year-old cave painting of an Atlantic salmon (Salmo salar) near the Vézère River in France is a reminder of our fascination with, and dependence on, Atlantic salmon throughout human history. Atlantic salmon belongs to the salmonid lineage which comprises 11 genera, with at least 70 species that exhibit a wide range of ecological adaptations and use a variety of marine and freshwater life history strategies 1 . Salmonids hold important positions as socially iconic species and economic resources within aquaculture, wild fisheries and recreational sport fisheries. Moreover, they serve as key indicator species of the health of North Atlantic and Pacific coastal and river ecosystems.All teleosts share at least three rounds of whole-genome duplication (WGD), 1R and 2R before the divergence of lamprey from the jawed vertebrates 2 , and a third teleost-specific WGD (Ts3R) at the base of the teleosts ~320 million years ago (Mya) [3][4][5] . Very little is known about the mechanisms of genomic and chromosomal reorganization after WGD in vertebrates because the 1R, 2R and Ts3R occurred so long ago that few clear signatures of post-WGD reorganization events remain. In contrast, a fourth WGD (the Ss4R salmonid-specific autotetraploidization event) occurred in the common ancestor of salmonids ~80 Mya after their divergence from Esociformes ~125 Mya 6-8 (Fig. 1), and the continued presence of multivalent pairing at meiosis and evidence of tetrasomic inheritance in salmonid species suggests that diploidy is not yet fully re-established 6,9,10 . Salmonids thus appear to provide an unprecedented opportunity for studying vertebrate genome evolution after an autotetraploid WGD 11,12 over a time period that is long enough to reveal long-term evolutionary patterns, but short enough to give a high-resolution picture of the process. In addition, they provide an excellent setting for contextualizing genome evolution with a dramatic post-WGD species radiation and intricate adaptations to a whole range of life history regimes.Here we present a high-quality reference genome assembly of the Atlantic salmon, and use it to describe major patterns characterizing the post-Ss4R salmonid genome evolution over the past 80 million years (Myr). Our results challenge the recent claim that rediploidization in salmonids has been a gradual process unlinked to significant genome rearrangements 13 . They also challenge current views about the relative importance of sub-and neofunctionalization in vertebrate genomes (reviewed in ref. 14), and the importance of dosage balance as a gene duplicate retention mechanism 15 . Genome characterizationThe Atlantic salmon reference genome assembly (GenBank: GCA_000233375.4) adds up to 2.97 gigabases (Gb) with aThe whole-genome duplication 80 million years ago of the common ancestor of salmonids (salmonid-specific fourth vertebrate whole-genome duplication, Ss4R) provides unique opportunities to learn about the evolutionary fate of a duplicated vertebrate genome in 70 extant lineages. Here we present a high...

Males and females share many traits that have a common genetic basis; however, selection on these traits often differs between the sexes, leading to sexual conflict. Under such sexual antagonism, theory predicts the evolution of genetic architectures that resolve this sexual conflict. Yet, despite intense theoretical and empirical interest, the specific loci underlying sexually antagonistic phenotypes have rarely been identified, limiting our understanding of how sexual conflict impacts genome evolution and the maintenance of genetic diversity. Here we identify a large effect locus controlling age at maturity in Atlantic salmon (Salmo salar), an important fitness trait in which selection favours earlier maturation in males than females, and show it is a clear example of sex-dependent dominance that reduces intralocus sexual conflict and maintains adaptive variation in wild populations. Using high-density single nucleotide polymorphism data across 57 wild populations and whole genome re-sequencing, we find that the vestigial-like family member 3 gene (VGLL3) exhibits sex-dependent dominance in salmon, promoting earlier and later maturation in males and females, respectively. VGLL3, an adiposity regulator associated with size and age at maturity in humans, explained 39% of phenotypic variation, an unexpectedly large proportion for what is usually considered a highly polygenic trait. Such large effects are predicted under balancing selection from either sexually antagonistic or spatially varying selection. Our results provide the first empirical example of dominance reversal allowing greater optimization of phenotypes within each sex, contributing to the resolution of sexual conflict in a major and widespread evolutionary trade-off between age and size at maturity. They also provide key empirical evidence for how variation in reproductive strategies can be maintained over large geographical scales. We anticipate these findings will have a substantial impact on population management in a range of harvested species where trends towards earlier maturation have been observed.

Atlantic cod is composed of multiple migratory and stationary populations widely distributed in the North Atlantic Ocean. The Northeast Arctic cod (NEAC) population in the Barents Sea undertakes annual spawning migrations to the northern Norwegian coast. Although spawning occurs sympatrically with the stationary Norwegian coastal cod (NCC), phenotypic and genetic differences between NEAC and NCC are maintained. In this study, we resolve the enigma by revealing the mechanisms underlying these differences. Extended linkage disequilibrium (LD) and population divergence were demonstrated in a 17.4-Mb region on linkage group 1 (LG1) based on genotypes of 494 SNPs from 192 parents of farmed families of NEAC, NCC or NEACxNCC crosses. Linkage analyses revealed two adjacent inversions within this region that repress meiotic recombination in NEACxNCC crosses. We identified a NEAC-specific haplotype consisting of 186 SNPs that was fixed in NEAC sampled from the Barents Sea, but segregating under Hardy-Weinberg equilibrium in eight NCC stocks. Comparative genomic analyses determine the NEAC configuration of the inversions to be the derived state and date it to ~1.6-2.0 Mya. The haplotype block harbours 763 genes, including candidates regulating swim bladder pressure, haem synthesis and skeletal muscle organization conferring adaptation to long-distance migrations and vertical movements down to large depths. Our results suggest that the migratory ecotype experiences strong directional selection for the two adjacent inversions on LG1. Despite interbreeding between NEAC and NCC, the inversions are maintaining genetic differentiation, and we hypothesize the co-occurrence of multiple adaptive alleles forming a 'supergene' in the NEAC population.

Males and females often differ in their fitness optima for shared traits that have a shared genetic basis, leading to sexual conflict. Morphologically differentiated sex chromosomes can resolve this conflict and protect sexually antagonistic variation, but they accumulate deleterious mutations. However, how sexual conflict is resolved in species that lack differentiated sex chromosomes is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 55-Mb double-inversion supergene that mediates sex-specific migratory tendency through sex-dependent dominance reversal, an alternative mechanism for resolving sexual conflict. The double inversion contains key photosensory, circadian rhythm, adiposity and sex-related genes and displays a latitudinal frequency cline, indicating environmentally dependent selection. Our results show sex-dependent dominance reversal across a large autosomal supergene, a mechanism for sexual conflict resolution capable of protecting sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutations associated with typical heteromorphic sex chromosomes. Methodology ReplicatesDescribe the experimental replicates, specifying number, type and replicate agreement. Sequencing depthDescribe the sequencing depth for each experiment, providing the total number of reads, uniquely mapped reads, length of reads and whether they were paired-or single-end. AntibodiesDescribe the antibodies used for the ChIP-seq experiments; as applicable, provide supplier name, catalog number, clone name, and lot number. Peak calling parametersSpecify the command line program and parameters used for read mapping and peak calling, including the ChIP, control and index files used. Data qualityDescribe the methods used to ensure data quality in full detail, including how many peaks are at FDR 5% and above 5-fold enrichment. SoftwareDescribe the software used to collect and analyze the ChIP-seq data. For custom code that has been deposited into a community repository, provide accession details. Flow Cytometry PlotsConfirm that:The axis labels state the marker and fluorochrome used (e.g. CD4-FITC).The axis scales are clearly visible. Include numbers along axes only for bottom left plot of group (a 'group' is an analysis of identical markers).All plots are contour plots with outliers or pseudocolor plots.A numerical value for number of cells or percentage (with statistics) is provided. Methodology Sample preparationDescribe the sample preparation, detailing the biological source of the cells and any tissue processing steps used. InstrumentIdentify the instrument used for data collection, specifying make and model number. SoftwareDescribe the software used to collect and analyze the flow cytometry data. For custom code that has been deposited into a community repository, provide accession details.Cell population abundance Describe the abundance of the relevant cell populations within post-sort fractions, providing details on the...

Traits with different fitness optima in males and females cause sexual conflict when they have a shared genetic basis. Heteromorphic sex chromosomes can resolve this conflict and protect sexually antagonistic polymorphisms but accumulate deleterious mutations. However, many taxa lack differentiated sex chromosomes, and how sexual conflict is resolved in these species is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 56 Mb double-inversion supergene that mediates sex-specific migration through sex-dependent dominance, a mechanism that reduces sexual conflict. The double-inversion contains key photosensory, circadian rhythm, adiposity, and sexual differentiation genes and displays frequency clines associated with latitude and temperature, revealing environmental dependence. Our results constitute the first example of sex-dependent dominance across a large autosomal supergene, a novel mechanism for sexual conflict resolution capable of protecting polygenic sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutation load of heteromorphic sex chromosomes.