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.

Structural variants (SVs) are a major source of genetic and phenotypic variation, but remain challenging to accurately type and are hence poorly characterized in most species. We present an approach for reliable SV discovery in non-model species using whole genome sequencing and report 15,483 high-confidence SVs in 492 Atlantic salmon (Salmo salar L.) sampled from a broad phylogeographic distribution. These SVs recover population genetic structure with high resolution, include an active DNA transposon, widely affect functional features, and overlap more duplicated genes retained from an ancestral salmonid autotetraploidization event than expected. Changes in SV allele frequency between wild and farmed fish indicate polygenic selection on behavioural traits during domestication, targeting brain-expressed synaptic networks linked to neurological disorders in humans. This study offers novel insights into the role of SVs in genome evolution and the genetic architecture of domestication traits, along with resources supporting reliable SV discovery in non-model species.

BackgroundClinical mastitis is an inflammation of the mammary gland and causes significant costs to dairy production. It is unfavourably genetically correlated to milk production, and, thus, knowledge of the mechanisms that underlie these traits would be valuable to improve both of them simultaneously through breeding. A quantitative trait locus (QTL) that affects both clinical mastitis and milk production has recently been fine-mapped to around 89 Mb on bovine chromosome 6 (BTA6), but identification of the gene that underlies this QTL was not possible due to the strong linkage disequilibrium between single nucleotide polymorphisms (SNPs) within this region. Our aim was to identify the gene and, if possible, the causal polymorphism(s) responsible for this QTL through association analysis of high-density SNPs and imputed full sequence data in combination with analyses of transcript and protein levels of the identified candidate gene.ResultsAssociations between SNPs and the studied traits were strongest for SNPs that were located within and immediately upstream of the group-specific component (GC) gene. This gene encodes the vitamin D-binding protein (DBP) and has multiple roles in immune defense and milk production. A 12-kb duplication that was identified downstream of this gene covered its last exon and segregated with the QTL allele that is associated with increased mastitis susceptibility and milk production. However, analyses of GC mRNA levels on the available samples revealed no differences in expression between animals having or lacking this duplication. Moreover, we detected no differences in the concentrations of DBP and its ligand vitamin D between the animals with different GC genotypes that were available for this study.ConclusionsOur results suggest GC as the gene that underlies the QTL for clinical mastitis and milk production. However, since only healthy animals were sampled for transcription and expression analyses, we could not draw any final conclusion on the absence of quantitative differences between animals with different genotypes. Future studies should investigate GC RNA expression and protein levels in cows with different genotypes during an infection.Electronic supplementary materialThe online version of this article (doi:10.1186/s12711-016-0257-2) contains supplementary material, which is available to authorized users.

Boar taint is characterized by an unpleasant taste or odor in intact male pigs and is primarily attributed to increased concentrations of androstenone and skatole and to a lesser extent by increased indole. The boar taint compounds skatole and indole are produced by gut bacteria, metabolized in the liver, and stored in the fat tissue. Androstenone, on the other hand, is synthesized in the testis along with testosterone and estrogens, which are known to be important factors affecting fertility. The main goal of this study was to investigate the relationship between genetic factors involved in the primary boar taint compounds in an attempt to discover ways to reduce boar taint without decreasing fertility-related compounds. Heritabilities and genetic correlations between traits were estimated for compounds related to boar taint (androstenone, skatole, indole) and reproduction (testosterone, 17β-estradiol, and estrone sulfate). Heritabilities in the range of 0.47 to 0.67 were detected for androstenone concentrations in both fat and plasma, whereas those for skatole and indole were slightly less (0.27 to 0.41). The genetic correlations between androstenone in plasma and fat were extremely high (0.91 to 0.98) in Duroc and Landrace. In addition, genetic correlations between androstenone (both plasma and fat) and the other sex steroids (estrone sulfate, 17β-estradiol, and testosterone) were very high, in the range of 0.80 to 0.95. Furthermore, a genome-wide association study (GWA) and a combined linkage disequilibrium and linkage analysis (LDLA) were conducted on 1,533 purebred Landrace and 1,027 purebred Duroc to find genome regions involved in genetic control of the boar taint compounds androstenone, skatole, and indole, and sex hormones related to fertility traits. Up to 3,297 informative SNP markers were included for both breeds, including SNP from several boar taint candidate genes. From the GWA study, we found that altogether 27 regions were significant at a genome-wide level (P < 0.05) and an additional 7 regions were significant at a chromosomal level. From the LDLA study, 7 regions were significant on a genome-wide level and an additional 7 regions were significant at a chromosomal level. The most convincing associations were obtained in 6 regions affecting skatole and indole in fat on chromosomes 1, 2, 3, 7, 13, and 14, 1 region on chromosome 6 affecting androstenone in plasma only, and 5 regions on chromosomes 3, 4, 13, and 15 affecting androstenone, testosterone, and estrogens.