The Atlantic salmon genome provides insights into rediploidization

In Atlantic salmon, as in most salmonids, males can mature early in the life cycle, as small freshwater fish, termed parr, and/or undergo a sea migration before maturing as full-size adults. The alternative life histories are contingent on environmental and social circumstances, such as growth rate, territory quality or any other factor that affects the individual's state. In order to model the choice of life history in this group of commercially valuable species, it is necessary to understand not only the relative contribution of the different male types to subsequent generations, but also to know the factors that affect reproductive success in each type. In this paper we present the results of a study designed to investigate the factors that affect the reproductive success of mature parr. We used highly polymorphic minisatellite DNA markers to analyse paternity in a series of mating experiments where the number and body size of parr were manipulated. The fraction of eggs fertilized by mature parr ranged from 26 to 40 per cent, with individual parr fertilizing up to 26 per cent of the eggs. A strong positive correlation was found between parr size and reproductive success. The relative success of parr decreased with increasing parr number. Data from this and other studies on variation in the timing and degree of parr reproductive success are discussed in relation to the evolution of male mating strategies and life history in salmonids.

Characterization of the Atlantic salmon (Salmo salar) brain-type fatty acid binding protein (fabp7) genes reveals the fates of teleost fabp7 genes following whole genome duplications

BackgroundQuantitative trait locus (QTL) studies show that variation in salinity tolerance in Arctic charr and rainbow trout has a genetic basis, even though both these species have low to moderate salinity tolerance capacities. QTL were observed to localize to homologous linkage group segments within putative chromosomal regions possessing multiple candidate genes. We compared salinity tolerance QTL in rainbow trout and Arctic charr to those detected in a higher salinity tolerant species, Atlantic salmon. The highly derived karyotype of Atlantic salmon allows for the assessment of whether disparity in salinity tolerance in salmonids is associated with differences in genetic architecture. To facilitate these comparisons, we examined the genomic synteny patterns of key candidate genes in the other model teleost fishes that have experienced three whole-genome duplication (3R) events which preceded a fourth (4R) whole genome duplication event common to all salmonid species.ResultsNine linkage groups contained chromosome-wide significant QTL (AS-2, -4p, -4q, -5, -9, -12p, -12q, -14q -17q, -22, and −23), while a single genome-wide significant QTL was located on AS-4q. Salmonid genomes shared the greatest marker homology with the genome of three-spined stickleback. All linkage group arms in Atlantic salmon were syntenic with at least one stickleback chromosome, while 18 arms had multiple affinities. Arm fusions in Atlantic salmon were often between multiple regions bearing salinity tolerance QTL. Nine linkage groups in Arctic charr and six linkage group arms in rainbow trout currently have no synteny alignments with stickleback chromosomes, while eight rainbow trout linkage group arms were syntenic with multiple stickleback chromosomes. Rearrangements in the stickleback lineage involving fusions of ancestral arm segments could account for the 21 chromosome pairs observed in the stickleback karyotype.ConclusionsSalinity tolerance in salmonids from three genera is to some extent controlled by the same loci. Synteny between QTL in salmonids and candidate genes in stickleback suggests genetic variation at candidate gene loci could affect salinity tolerance in all three salmonids investigated. Candidate genes often occur in pairs on chromosomes, and synteny patterns indicate these pairs are generally conserved in 2R, 3R, and 4R genomes. Synteny maps also suggest that the Atlantic salmon genome contains three larger syntenic combinations of candidate genes that are not evident in any of the other 2R, 3R, or 4R genomes examined. These larger synteny tracts appear to have resulted from ancestral arm fusions that occurred in the Atlantic salmon ancestor. We hypothesize that the superior hypo-osmoregulatory efficiency that is characteristic of Atlantic salmon may be related to these clusters.

Genomic organization and characterization of two vomeronasal 1 receptor-like genes ( ora1 and ora2) in Atlantic salmon Salmo salar

Salmonids represent an intriguing taxonomical group for investigating genome evolution in vertebrates due to their relatively recent last common whole genome duplication event, which occurred between 80 and 100 million years ago. Here, we report on the chromosome-level genome assembly of European grayling (Thymallus thymallus), which represents one of the earliest diverged salmonid subfamilies. To achieve this, we first generated relatively long genomic scaffolds by using a previously published draft genome assembly along with long-read sequencing data and a linkage map. We then merged those scaffolds by applying synteny evidence from the Atlantic salmon (Salmo salar) genome. Comparisons of the European grayling genome assembly to the genomes of Atlantic salmon and Northern pike (Esox lucius), the latter used as a nonduplicated outgroup, detailed aspects of the characteristic chromosome evolution process that has taken place in European grayling. While Atlantic salmon and other salmonid genomes are portrayed by the typical occurrence of numerous chromosomal fusions, European grayling chromosomes were confirmed to be fusion-free and were characterized by a relatively large proportion of paracentric and pericentric inversions. We further reported on transposable elements specific to either the European grayling or Atlantic salmon genome, on the male-specific sdY gene in the European grayling chromosome 11A, and on regions under residual tetrasomy in the homeologous European grayling chromosome pairs 9A-9B and 25A-25B. The same chromosome pairs have been observed under residual tetrasomy in Atlantic salmon and in other salmonids, suggesting that this feature has been conserved since the subfamily split.

BackgroundThe Atlantic salmon genome is in the process of returning to a diploid state after undergoing a whole genome duplication (WGD) event between 25 and100 million years ago. Existing data on the proportion of paralogous sequence variants (PSVs), multisite variants (MSVs) and other types of complex sequence variation suggest that the rediplodization phase is far from over. The aims of this study were to construct a high density linkage map for Atlantic salmon, to characterize the extent of rediploidization and to improve our understanding of genetic differences between sexes in this species.ResultsA linkage map for Atlantic salmon comprising 29 chromosomes and 5650 single nucleotide polymorphisms (SNPs) was constructed using genotyping data from 3297 fish belonging to 143 families. Of these, 2696 SNPs were generated from ESTs or other gene associated sequences. Homeologous chromosomal regions were identified through the mapping of duplicated SNPs and through the investigation of syntenic relationships between Atlantic salmon and the reference genome sequence of the threespine stickleback (Gasterosteus aculeatus). The sex-specific linkage maps spanned a total of 2402.3 cM in females and 1746.2 cM in males, highlighting a difference in sex specific recombination rate (1.38:1) which is much lower than previously reported in Atlantic salmon. The sexes, however, displayed striking differences in the distribution of recombination sites within linkage groups, with males showing recombination strongly localized to telomeres.ConclusionThe map presented here represents a valuable resource for addressing important questions of interest to evolution (the process of re-diploidization), aquaculture and salmonid life history biology and not least as a resource to aid the assembly of the forthcoming Atlantic salmon reference genome sequence.

The ability of Atlantic salmon ( Salmo salar L.) parr to discriminate between the urine from sibling and non-sibling fish was studied by using electrophysiological and behavioural techniques. Urine was collected from four sibling groups of Atlantic salmon parr which had been reared separately. Urine from all four groups were potent odorants in the Atlantic salmon parr, eliciting responses that were recordable from the olfactory epithelium. However, the mean recorded responses were significantly greater in Atlantic salmon parr that were stimulated with the urine from sibling than from non-sibling fish. Urine from all four groups elicited display behaviour in the Atlantic salmon parr which often occurs during aggressive interactions between parr when defending territories. Atlantic salmon parr also moved towards the source of sibling urine but moved away from non-sibling urine. The results are discussed in relation to the role of urine and the adaptive significance of sibling recognition in the Atlantic salmon.

The melanocortin system is a key neuroendocrine network involved in the control of food intake and energy homeostasis in vertebrates. Within the hypothalamus, the system comprises two main distinct neuronal cell populations that express the neuropeptides proopiomelanocortin (POMC; anorexigenic) or agouti-related protein (AGRP; orexigenic). Both bind to the melanocortin-4 receptor (MC4R) in higher order neurons that control both food intake and energy expenditure. This system is relatively well-conserved among vertebrates. However, in Atlantic salmon (Salmo salar L.), the salmonid-specific fourth round whole-genome duplication led to the presence of several paralog genes which might result in divergent functions of the duplicated genes. In the current study, we report the first comprehensive comparative identification and characterization of Mc4r and extend the knowledge of Pomc and Agrp in appetite control in Atlantic salmon. In silico analysis revealed multiple paralogs for mc4r (a1, a2, b1, and b2) in the Atlantic salmon genome and confirmed the paralogs previously described for pomc (a1, a2, and b) and agrp (1 and 2). All Mc4r paralogs are relatively well-conserved with the human homolog, sharing at least 63% amino acid sequence identity. We analyzed the mRNA expression of mc4r, pomc, and agrp genes in eight brain regions of Atlantic salmon post-smolt under two feeding states: normally fed and fasted for 4 days. The mc4ra2 and b1 mRNAs were predominantly and equally abundant in the hypothalamus and telencephalon, the mc4rb2 in the hypothalamus, and a1 in the telencephalon. All pomc genes were highly expressed in the pituitary, followed by the hypothalamus and saccus vasculosus. The agrp genes showed a completely different expression pattern from each other, with prevalent expression of the agrp1 in the hypothalamus and agrp2 in the telencephalon. Fasting did not induce any significant changes in the mRNA level of mc4r, agrp, or pomc paralogs in the hypothalamus or in other highly expressed regions between fed and fasted states. The identification and wide distribution of multiple paralogs of mc4r, pomc, and agrp in Atlantic salmon brain provide new insights and give rise to new questions of the melanocortin system in the appetite regulation in Atlantic salmon.

The functional and evolutionary outcomes of whole genome duplication (WGD) events are driven by global remodeling of gene expression. Most investigations of gene expression changes following WGD have applied bulk transcriptomics using tissue samples, thus failing to resolve affected cell types. Here, we leverage single-cell transcriptomics of liver tissue in Atlantic salmon (Salmo salar L.) to quantify cell-specific expression and transcriptional responses to a bacterial infection with Aeromonas salmonicida for thousands of duplicate gene pairs (ohnologs) retained from WGD ancestral to all salmonids. The major liver cell types showed hundreds of differentially expressed ohnolog pairs, with hepatocytes showing the greatest number and immune cells the least number of uniquely differentially expressed pairs. Many more differentially expressed ohnolog pairs were identified after accounting for cell type heterogeneity within a cell lineage, despite a reduction in statistical power. The degree of conservation in ohnolog expression responses to bacterial infection also varied significantly among cell types, both in terms of the number of differentially expressed pairs and the direction of responses. Overall, this study highlights the importance of resolving cell-specific gene expression to understand the functional and evolutionary outcomes of WGD events.

Transposable elements (TEs) are hypothesized to play important roles in shaping genome evolution following whole-genome duplications (WGDs), including rewiring of gene regulation. In a recent analysis, duplicate gene copies that had evolved higher expression in liver following the salmonid WGD ∼100 million years ago were associated with higher numbers of predicted TE-derived cis-regulatory elements (TE-CREs). Yet, the ability of these TE-CREs to recruit transcription factors (TFs) in vivo and impact gene expression remains unknown. Here, we evaluated the gene-regulatory functions of 11 TEs using luciferase promoter reporter assays in Atlantic salmon (Salmo salar) primary liver cells. Canonical Tc1-Mariner elements from intronic regions showed no or small repressive effects on transcription. However, other TE-CREs upstream of transcriptional start sites increased expression significantly. Our results question the hypothesis that TEs in the Tc1-Mariner superfamily, which were extremely active following WGD in salmonids, had a major impact on regulatory rewiring of gene duplicates, but highlights the potential of other TEs in post-WGD rewiring of gene regulation in the Atlantic salmon genome.

Physical localization and characterization of the BglI element in the genomes of Atlantic salmon ( Salmo salar L.) and brown trout ( S. trutta L.)

BackgroundGenetic linkage maps are useful tools for mapping quantitative trait loci (QTL) influencing variation in traits of interest in a population. Genotyping-by-sequencing approaches such as Restriction-site Associated DNA sequencing (RAD-Seq) now enable the rapid discovery and genotyping of genome-wide SNP markers suitable for the development of dense SNP linkage maps, including in non-model organisms such as Atlantic salmon (Salmo salar). This paper describes the development and characterisation of a high density SNP linkage map based on SbfI RAD-Seq SNP markers from two Atlantic salmon reference families.ResultsApproximately 6,000 SNPs were assigned to 29 linkage groups, utilising markers from known genomic locations as anchors. Linkage maps were then constructed for the four mapping parents separately. Overall map lengths were comparable between male and female parents, but the distribution of the SNPs showed sex-specific patterns with a greater degree of clustering of sire-segregating SNPs to single chromosome regions. The maps were integrated with the Atlantic salmon draft reference genome contigs, allowing the unique assignment of ~4,000 contigs to a linkage group. 112 genome contigs mapped to two or more linkage groups, highlighting regions of putative homeology within the salmon genome. A comparative genomics analysis with the stickleback reference genome identified putative genes closely linked to approximately half of the ordered SNPs and demonstrated blocks of orthology between the Atlantic salmon and stickleback genomes. A subset of 47 RAD-Seq SNPs were successfully validated using a high-throughput genotyping assay, with a correspondence of 97% between the two assays.ConclusionsThis Atlantic salmon RAD-Seq linkage map is a resource for salmonid genomics research as genotyping-by-sequencing becomes increasingly common. This is aided by the integration of the SbfI RAD-Seq SNPs with existing reference maps and the draft reference genome, as well as the identification of putative genes proximal to the SNPs. Differences in the distribution of recombination events between the sexes is evident, and regions of homeology have been identified which are reflective of the recent salmonid whole genome duplication.

BackgroundThe genomes of salmonids are considered pseudo-tetraploid undergoing reversion to a stable diploid state. Given the genome duplication and extensive biological data available for salmonids, they are excellent model organisms for studying comparative genomics, evolutionary processes, fates of duplicated genes and the genetic and physiological processes associated with complex behavioral phenotypes. The evolution of the tetrapod hemoglobin genes is well studied; however, little is known about the genomic organization and evolution of teleost hemoglobin genes, particularly those of salmonids. The Atlantic salmon serves as a representative salmonid species for genomics studies. Given the well documented role of hemoglobin in adaptation to varied environmental conditions as well as its use as a model protein for evolutionary analyses, an understanding of the genomic structure and organization of the Atlantic salmon α and β hemoglobin genes is of great interest.ResultsWe identified four bacterial artificial chromosomes (BACs) comprising two hemoglobin gene clusters spanning the entire α and β hemoglobin gene repertoire of the Atlantic salmon genome. Their chromosomal locations were established using fluorescence in situ hybridization (FISH) analysis and linkage mapping, demonstrating that the two clusters are located on separate chromosomes. The BACs were sequenced and assembled into scaffolds, which were annotated for putatively functional and pseudogenized hemoglobin-like genes. This revealed that the tail-to-tail organization and alternating pattern of the α and β hemoglobin genes are well conserved in both clusters, as well as that the Atlantic salmon genome houses substantially more hemoglobin genes, including non-Bohr β globin genes, than the genomes of other teleosts that have been sequenced.ConclusionsWe suggest that the most parsimonious evolutionary path leading to the present organization of the Atlantic salmon hemoglobin genes involves the loss of a single hemoglobin gene cluster after the whole genome duplication (WGD) at the base of the teleost radiation but prior to the salmonid-specific WGD, which then produced the duplicated copies seen today. We also propose that the relatively high number of hemoglobin genes as well as the presence of non-Bohr β hemoglobin genes may be due to the dynamic life history of salmon and the diverse environmental conditions that the species encounters.Data deposition: BACs S0155C07 and S0079J05 (fps135): GenBank GQ898924; BACs S0055H05 and S0014B03 (fps1046): GenBank GQ898925

Domestication and artificial selection for desirable traits have driven significant phenotypic changes and left detectable genomic footprints in farmed animals. Since the 1960s, intensive breeding has led to the rapid domestication of Atlantic salmon (Salmo salar), with multiple independent events that make it a valuable model for studying early domestication stages and the parallel evolution of populations of different origins subjected to similar selection pressures. Some aquatic species, including Atlantic salmon, have undergone whole-genome duplication (WGD), raising the possibility that genetic redundancy resulting from WGD has contributed to adaptation in captive environments, as seen in plants. Here, we examined the genomic responses to domestication in Atlantic salmon, focusing on potential signatures of parallel selection, including those associated with WGD. Candidate genomic regions under selection were identified by comparing whole-genome sequences from aquaculture and wild populations across 2 independently domesticated lineages (Western Norway and North America) using a genome-wide scan that combined 3 statistical methods: allele frequencies (FST), site frequency (Tajima's D), and haplotype differentiation (XP-EHH). These analyses revealed shared selective sweeps on identical SNPs in major histocompatibility complex (MHC) genes across aquaculture populations. This suggests that a combination of long-term balancing selection and recent human-induced selection has shaped MHC gene evolution in domesticated salmon. Additionally, we observed selective sweeps on a small number of gene pairs in homeologous regions originating from WGD, offering insights into how historical genome duplication events may intersect with recent selection pressures in aquaculture species.

Basal and immune-responsive transcript expression of two Atlantic salmon interferon regulatory factor 2 (irf2) paralogues.