The Atlantic salmon genome provides insights into rediploidization.

The 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-quality genome assembly for Atlantic salmon (Salmo salar), and show that large genomic reorganizations, coinciding with bursts of transposon-mediated repeat expansions, were crucial for the post-Ss4R rediploidization process. Comparisons of duplicate gene expression patterns across a wide range of tissues with orthologous genes from a pre-Ss4R outgroup unexpectedly demonstrate far more instances of neofunctionalization than subfunctionalization. Surprisingly, we find that genes that were retained as duplicates after the teleost-specific whole-genome duplication 320 million years ago were not more likely to be retained after the Ss4R, and that the duplicate retention was not influenced to a great extent by the nature of the predicted protein interactions of the gene products. Finally, we demonstrate that the Atlantic salmon assembly can serve as a reference sequence for the study of other salmonids for a range of purposes.

Immersion challenge model for Flavobacterium psychrophilum infection of Atlantic salmon (Salmo salar L.) fry.

Flavobacterium psychrophilum is the causative agent of bacterial cold-water disease (CWBD) and rainbow trout fry syndrome (RTFS), which affect salmonids. To better understand this pathogen and its interaction with the host during infection, including to support the development of resistant breeds and new vaccines and treatments, there is a pressing need for reliable and reproducible immersion challenge models that more closely mimic natural routes of infection. The aim of this present study was to evaluate a challenge model developed previously for rainbow trout for use in Atlantic salmon. First, preliminary challenges were conducted in Atlantic salmon (n = 120) and rainbow trout (n = 80) fry using two F. psychrophilum isolates collected from each fish species, respectively; fish had been pretreated with 200 mg/L hydrogen peroxide for 1 h. Thereafter, the main challenge was performed for just one F. psychrophilum isolate for each species (at 2 × 107 CFU/mL) but using larger cohorts (Atlantic salmon: n = 1187; rainbow trout: n = 2701). Survival in the main challenge was 81.2% in Atlantic salmon (21 days post-challenge) and 45.3% in rainbow trout (31 days post-challenge). Mortalities progressed similarly during the preliminary and main challenges for both species, demonstrating the reproducibility of this model. This is the first immersion challenge model of F. psychrophilum to be developed successfully for Atlantic salmon.

Bayesian genomic models boost prediction accuracy for survival to Streptococcus agalactiae infection in Nile tilapia (Oreochromus nilioticus).

BACKGROUND: Streptococcosis is a major bacterial disease in Nile tilapia that is caused by Streptococcus agalactiae infection, and development of resistant strains of Nile tilapia represents a sustainable approach towards combating this disease. In this study, we performed a controlled disease trial on 120 full-sib families to (i) quantify and characterize the potential of genomic selection for survival to S. agalactiae infection in Nile tilapia, and (ii) identify the best genomic model and the optimal density of single nucleotide polymorphisms (SNPs) for this trait. METHODS: In total, 40 fish per family (15 fish intraperitoneally injected and 25 fish as cohabitants) were used in the challenge test. Mortalities were recorded every 3 h for 35 days. After quality control, genotypes (50,690 SNPs) and phenotypes (0 for dead and 1 for alive) for 2472 cohabitant fish were available. Genetic parameters were obtained using various genomic selection models (genomic best linear unbiased prediction (GBLUP), BayesB, BayesC, BayesR and BayesS) and a traditional pedigree-based model (PBLUP). The pedigree-based analysis used a deep 17-generation pedigree. Prediction accuracy and bias were evaluated using five replicates of tenfold cross-validation. The genomic models were further analyzed using 10 subsets of SNPs at different densities to explore the effect of pruning and SNP density on predictive accuracy. RESULTS: Moderate estimates of heritabilities ranging from 0.15 ± 0.03 to 0.26 ± 0.05 were obtained with the different models. Compared to a pedigree-based model, GBLUP (using all the SNPs) increased prediction accuracy by 15.4%. Furthermore, use of the most appropriate Bayesian genomic selection model and SNP density increased the prediction accuracy up to 71%. The 40 to 50 SNPs with non-zero effects were consistent for all BayesB, BayesC and BayesS models with respect to marker id and/or marker locations. CONCLUSIONS: These results demonstrate the potential of genomic selection for survival to S. agalactiae infection in Nile tilapia. Compared to the PBLUP and GBLUP models, Bayesian genomic models were found to boost the prediction accuracy significantly.

Sex-dependent dominance at a single locus maintains variation in age at maturity in salmon.

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.

Genomic regions and signaling pathways associated with indicator traits for feed efficiency in juvenile Atlantic salmon (Salmo salar).

BACKGROUND: One objective of this study was to identify putative quantitative trait loci (QTL) that affect indicator phenotypes for growth, nitrogen, and carbon metabolism in muscle, liver, and adipose tissue, and for feed efficiency. Another objective was to perform an RNAseq analysis (184 fish from all families), to identify genes that are associated with carbon and nitrogen metabolism in the liver. The material consisted of a family experiment that was performed in freshwater and included 2281 individuals from 23 full-sib families. During the 12-day feed conversion test, families were randomly allocated to family tanks (50 fish per tank and 2 tanks per family) and fed a fishmeal-based diet labeled with the stable isotopes 15N and 13C at inclusion levels of 2 and 1%, respectively. RESULTS: Using a linear mixed-model algorithm, a QTL for pre-smolt growth was identified on chromosome 9 and a QTL for carbon metabolism in the liver was identified on chromosome 12 that was closely related to feed conversion ratio on a tank level. For the indicators of feed efficiency traits that were derived from the stable isotope ratios (15N and 13C) of muscle tissue and growth, no convincing QTL was detected, which suggests that these traits are polygenic. The transcriptomic analysis showed that high carbon and nitrogen metabolism was associated with individuals that convert protein from the feed more efficiently, primarily due to higher expression of the proteasome, lipid, and carbon metabolic pathways in liver. In addition, we identified seven transcription factors that were associated with carbon and nitrogen metabolism and located in the identified QTL regions. CONCLUSIONS: Analyses revealed one QTL associated with pre-smolt growth and one QTL for carbon metabolism in the liver. Both of these traits are associated with feed efficiency. However, more accurate mapping of the putative QTL will require a more diverse family material. In this experiment, fish that have a high carbon and nitrogen metabolism in the liver converted protein from the feed more efficiently, potentially because of a higher expression of the proteasome, lipid, and carbon metabolic pathways in liver. Within the QTL regions, we detected seven transcription factors that were associated with carbon and nitrogen metabolism.

Beyond large-effect loci: large-scale GWAS reveals a mixed large-effect and polygenic architecture for age at maturity of Atlantic salmon.

BACKGROUND: Understanding genetic architecture is essential for determining how traits will change in response to evolutionary processes such as selection, genetic drift and/or gene flow. In Atlantic salmon, age at maturity is an important life history trait that affects factors such as survival, reproductive success, and growth. Furthermore, age at maturity can seriously impact aquaculture production. Therefore, characterizing the genetic architecture that underlies variation in age at maturity is of key interest. RESULTS: Here, we refine our understanding of the genetic architecture for age at maturity of male Atlantic salmon using a genome-wide association study of 11,166 males from a single aquaculture strain, using imputed genotypes at 512,397 single nucleotide polymorphisms (SNPs). All individuals were genotyped with a 50K SNP array and imputed to higher density using parents genotyped with a 930K SNP array and pedigree information. We found significant association signals on 28 of 29 chromosomes (P-values: 8.7 × 10-133-9.8 × 10-8), including two very strong signals spanning the six6 and vgll3 gene regions on chromosomes 9 and 25, respectively. Furthermore, we identified 116 independent signals that tagged 120 candidate genes with varying effect sizes. Five of the candidate genes found here were previously associated with age at maturity in other vertebrates, including humans. DISCUSSION: These results reveal a mixed architecture of large-effect loci and a polygenic component that consists of multiple smaller-effect loci, suggesting a more complex genetic architecture of Atlantic salmon age at maturity than previously thought. This more complex architecture will have implications for selection on this key trait in aquaculture and for management of wild salmon populations.

Transcriptomic analysis of spleen infected with infectious salmon anemia virus reveals distinct pattern of viral replication on resistant and susceptible Atlantic salmon (Salmo salar).

The infectious salmon anemia virus (ISAv) produces a systemic infection in salmonids, causing large losses in salmon production. However, little is known regarding the mechanisms exerting disease resistance. In this paper, we perform an RNA-seq analysis in Atlantic salmon challenged with ISAv (using individuals coming from families that were highly susceptible or highly resistant to ISAv infection). We evaluated the differential expression of both host and ISAv genes in a target organ for the virus, i.e. the spleen. The results showed differential expression of host genes related to response to stress, immune response and protein folding (genes such as; atf3, mhc, mx1-3, cd276, cd2, cocs1, c7, il10, il10rb, il13ra2, ubl-1, ifng, ifngr1, hivep2, sigle14 and sigle5). An increased protein processing activity was found in susceptible fish, which generates a subsequent unfolded protein response. We observed extreme differences in the expression of viral segments between susceptible and resistant groups, demonstrating the capacity of resistant fish to overcome the virus replication, generating a very low viral load. This phenomenon and survival of this higher resistant fish seem to be related to differences in immune and translational process, as well as to the increase of HIV-EP2 (hivep2) transcript in resistant fish, although the causal mechanism is yet to be discovered. This study provides valuable information about disease resistance mechanisms in Atlantic salmon from a host-pathogen interaction point of view.

A dense SNP-based linkage map for Atlantic salmon (Salmo salar) reveals extended chromosome homeologies and striking differences in sex-specific recombination patterns.

BACKGROUND: The 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. RESULTS: A 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. CONCLUSION: The 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.

Identification of a novel allele of peroxisome proliferator-activated receptor gamma (PPARG) and its association with resistance to Aeromonas salmonicida in Atlantic salmon (Salmo salar).

Bacterial and viral diseases are major problems in Atlantic salmon aquaculture, but may be challenged through selection of brood stock with enhanced survival to diseases. Today's selection strategy is based on controlled challenge tests using siblings of the breeding candidates, and is thus indirect. Direct trait records on breeding candidates can potentially be provided through identification of genetic variation linked to the susceptibility to the disease. Peroxisome proliferator-activated receptor gamma (PPARG) is a lipid-sensing transcription factor primarily known for inducing fat-accumulation in adipocytes, but also in lipid-accumulating macrophages, in mammalian species. Here we report a novel allele of PPARG, pparg-2, in Atlantic salmon. pparg-2 has an insertion of sixty nucleotides that encodes two additional copies of the almost perfect decapeptide motif, (F/C/Y)NHSPDR(S/N)HS, compared to the previously described pparg-1. pparg-1 contains six copies of this repeat unit whereas eight copies are present in the novel pparg-2 allele. pparg-2 mRNA was detectable in kidney and spleen of random Atlantic salmon samples. Here, we studied the effect of pparg-1 and pparg-2 on survival upon challenge to a highly virulent bacterium, Aeromonas salmonicida, causing furunculosis, and the virus causing infectious salmon anaemia (ISA), respectively, in a Norwegian aquaculture population of Atlantic salmon. ppar alleles were found to be significantly associated with survival upon challenge to A. salmonicida, but not to ISA. pparg-2 was the better allele in terms of survival in the challenge test for furunculosis, survival rates being 0.32, 0.40 and 0.42 for animals with the pparg-1,-1, pparg-1, -2 and pparg-2, -2 genotypes, respectively. We conclude that pparg-2 is in linkage disequilibrium (LD) with, or identical to, a locus contributing to different susceptibility to furunculosis in Atlantic salmon. PPARG was mapped to linkage group eight (LG8) but could only be positioned on the male linkage map since all the informative parents in the mapping families were males. This is the first report showing an association between pparg alleles and an enhanced immune response in fish.

Mapping of quantitative trait loci for flesh colour and growth traits in Atlantic salmon (Salmo salar).

BACKGROUND: Flesh colour and growth related traits in salmonids are both commercially important and of great interest from a physiological and evolutionary perspective. The aim of this study was to identify quantitative trait loci (QTL) affecting flesh colour and growth related traits in an F2 population derived from an isolated, landlocked wild population in Norway (Byglands Bleke) and a commercial production population. METHODS: One hundred and twenty-eight informative microsatellite loci distributed across all 29 linkage groups in Atlantic salmon were genotyped in individuals from four F2 families that were selected from the ends of the flesh colour distribution. Genotyping of 23 additional loci and two additional families was performed on a number of linkage groups harbouring putative QTL. QTL analysis was performed using a line-cross model assuming fixation of alternate QTL alleles and a half-sib model with no assumptions about the number and frequency of QTL alleles in the founder populations. RESULTS: A moderate to strong phenotypic correlation was found between colour, length and weight traits. In total, 13 genome-wide significant QTL were detected for all traits using the line-cross model, including three genome-wide significant QTL for flesh colour (Chr 6, Chr 26 and Chr 4). In addition, 32 suggestive QTL were detected (chromosome-wide P < 0.05). Using the half-sib model, six genome-wide significant QTL were detected for all traits, including two for flesh colour (Chr 26 and Chr 4) and 41 suggestive QTL were detected (chromosome-wide P < 0.05). Based on the half-sib analysis, these two genome-wide significant QTL for flesh colour explained 24% of the phenotypic variance for this trait. CONCLUSIONS: A large number of significant and suggestive QTL for flesh colour and growth traits were found in an F2 population of Atlantic salmon. Chr 26 and Chr 4 presented the strongest evidence for significant QTL affecting flesh colour, while Chr 10, Chr 5, and Chr 4 presented the strongest evidence for significant QTL affecting growth traits (length and weight). These QTL could be strong candidates for use in marker-assisted selection and provide a starting point for further characterisation of the genetic components underlying flesh colour and growth.

Publisher Correction: Sex-dependent dominance maintains migration supergene in rainbow trout.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genomic signatures of local directional selection in a high gene flow marine organism; the Atlantic cod (Gadus morhua).

BACKGROUND: Marine fishes have been shown to display low levels of genetic structuring and associated high levels of gene flow, suggesting shallow evolutionary trajectories and, possibly, limited or lacking adaptive divergence among local populations. We investigated variation in 98 gene-associated single nucleotide polymorphisms (SNPs) for evidence of selection in local populations of Atlantic cod (Gadus morhua L.) across the species distribution. RESULTS: Our global genome scan analysis identified eight outlier gene loci with very high statistical support, likely to be subject to directional selection in local demes, or closely linked to loci under selection. Likewise, on a regional south/north transect of central and eastern Atlantic populations, seven loci displayed strongly elevated levels of genetic differentiation. Selection patterns among populations appeared to be relatively widespread and complex, i.e. outlier loci were generally not only associated with one of a few divergent local populations. Even on a limited geographical scale between the proximate North Sea and Baltic Sea populations four loci displayed evidence of adaptive evolution. Temporal genome scan analysis applied to DNA from archived otoliths from a Faeroese population demonstrated stability of the intra-population variation over 24 years. An exploratory landscape genetic analysis was used to elucidate potential effects of the most likely environmental factors responsible for the signatures of local adaptation. We found that genetic variation at several of the outlier loci was better correlated with temperature and/or salinity conditions at spawning grounds at spawning time than with geographic distance per se. CONCLUSION: These findings illustrate that adaptive population divergence may indeed be prevalent despite seemingly high levels of gene flow, as found in most marine fishes. Thus, results have important implications for our understanding of the interplay of evolutionary forces in general, and for the conservation of marine biodiversity under rapidly increasing evolutionary pressure from climate and fisheries induced changes in local environments.

Confirmation and fine-mapping of a major QTL for resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar): population-level associations between markers and trait.

BACKGROUND: Infectious pancreatic necrosis (IPN) is one of the most prevalent and economically devastating diseases in Atlantic salmon (Salmo salar) farming worldwide. The disease causes large mortalities at both the fry- and post-smolt stages. Family selection for increased IPN resistance is performed through the use of controlled challenge tests, where survival rates of sib-groups are recorded. However, since challenge-tested animals cannot be used as breeding candidates, within-family selection is not performed and only half of the genetic variation for IPN resistance is being exploited. DNA markers linked to quantitative trait loci (QTL) affecting IPN resistance would therefore be a powerful selection tool. The aim of this study was to identify and fine-map QTL for IPN-resistance in Atlantic salmon, for use in marker-assisted selection to increase the rate of genetic improvement for this trait. RESULTS: A genome scan was carried out using 10 large full-sib families of challenge-tested Atlantic salmon post-smolts and microsatellite markers distributed across the genome. One major QTL for IPN-resistance was detected, explaining 29% and 83% of the phenotypic and genetic variances, respectively. This QTL mapped to the same location as a QTL recently detected in a Scottish Atlantic salmon population. The QTL was found to be segregating in 10 out of 20 mapping parents, and subsequent fine-mapping with additional markers narrowed the QTL peak to a 4 cM region on linkage group 21. Challenge-tested fry were used to show that the QTL had the same effect on fry as on post-smolt, with the confidence interval for QTL position in fry overlapping the confidence interval found in post-smolts. A total of 178 parents were tested for segregation of the QTL, identifying 72 QTL-heterozygous parents. Genotypes at QTL-heterozygous parents were used to determine linkage phases between alleles at the underlying DNA polymorphism and alleles at single markers or multi-marker haplotypes. One four-marker haplotype was found to be the best predictor of QTL alleles, and was successfully used to deduce genotypes of the underlying polymorphism in 72% of the parents of the next generation within a breeding nucleus. A highly significant population-level correlation was found between deduced alleles at the underlying polymorphism and survival of offspring groups in the fry challenge test, parents with the three deduced genotypes (QQ, Qq, qq) having mean offspring mortality rates of 0.13, 0.32, and 0.49, respectively. The frequency of the high-resistance allele (Q) in the population was estimated to be 0.30. Apart from this major QTL, one other experiment-wise significant QTL for IPN-resistance was detected, located on linkage group 4. CONCLUSION: The QTL confirmed in this study represents a case of a major gene explaining the bulk of genetic variation for a presumed complex trait. QTL genotypes were deduced within most parents of the 2005 generation of a major breeding company, providing a solid framework for linkage-based MAS within the whole population in subsequent generations. Since haplotype-trait associations valid at the population level were found, there is also a potential for MAS based on linkage disequilibrium (LD). However, in order to use MAS across many generations without reassessment of linkage phases between markers and the underlying polymorphism, the QTL needs to be positioned with even greater accuracy. This will require higher marker densities than are currently available.

Genomic and functional gene studies suggest a key role of beta-carotene oxygenase 1 like (bco1l) gene in salmon flesh color.

Red coloration of muscle tissue (flesh) is a unique trait in several salmonid genera, including Atlantic salmon. The color results from dietary carotenoids deposited in the flesh, whereas the color intensity is affected both by diet and genetic components. Herein we report on a genome-wide association study (GWAS) to identify genetic variation underlying this trait. Two SNPs on ssa26 showed strong associations to the flesh color in salmon. Two genes known to be involved in carotenoid metabolism were located in this QTL- region: beta-carotene oxygenase 1 (bco1) and beta-carotene oxygenase 1 like (bco1l). To determine whether flesh color variation is caused by one, or both, of these genes, functional studies were carried out including mRNA and protein expression in fish with red and pale flesh color. The catalytic abilities of these two genes were also tested with different carotenoids. Our results suggest bco1l to be the most likely gene to explain the flesh color variation observed in this population.

Sex-dependent dominance maintains migration supergene in rainbow trout.

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.

A linkage map of the Atlantic salmon (Salmo salar) based on EST-derived SNP markers.

BACKGROUND: The Atlantic salmon is a species of commercial and ecological significance. Like other salmonids, the species displays residual tetrasomy and a large difference in recombination rate between sexes. Linkage maps with full genome coverage, containing both type I and type II markers, are needed for progress in genomics. Furthermore, it is important to estimate levels of linkage disequilibrium (LD) in the species. In this study, we developed several hundred single nucleotide polymorphism (SNP) markers for the Atlantic salmon, and constructed male and female linkage maps containing SNP and microsatellite markers. We also investigated further the distribution of male and female recombination events across the genome, and estimated levels of LD between pairs of markers. RESULTS: The male map had 29 linkage groups and was 390 cM long. The female map had 30 linkage groups as was 1983 cM long. In total, the maps contained 138 microsatellite markers and 304 SNPs located within genes, most of which were successfully annotated. The ratio of male to female recombination events was either close to zero or very large, indicating that there is little overlap between regions in which male and female crossovers occur. The female map is likely to have close to full genome coverage, while the majority of male linkage groups probably lack markers in telomeric regions where male recombination events occur. Levels of r2 increased with decreasing inter-marker distance in a bimodal fashion; increasing slowly from approximately 60 cM, and more rapidly more from approximately 12 cM. Long-ranging LD may be consequence of recent admixture in the population, the population being a 'synthetic' breeding population with contributions from several distinct rivers. Levels of r2 dropped to half its maximum value (above baseline) within 15 cM, and were higher than 0.2 above baseline for unlinked markers ('useful LD') at inter-marker distances less than 5 cM. CONCLUSION: The linkage map presented here is an important resource for genetic, comparative, and physical mapping of the Atlantic salmon. The female map is likely to have a map coverage that is not far from complete, whereas the male map length is likely to be significantly shorter than the true map, due to suboptimal marker coverage in the apparently small physical regions where male crossovers occur. 'Useful LD' was found at inter-marker distances less than 5 cM.

Identification and characterisation of novel SNP markers in Atlantic cod: evidence for directional selection.

BACKGROUND: The Atlantic cod (Gadus morhua) is a groundfish of great economic value in fisheries and an emerging species in aquaculture. Genetic markers are needed to identify wild stocks in order to ensure sustainable management, and for marker-assisted selection and pedigree determination in aquaculture. Here, we report on the development and evaluation of a large number of Single Nucleotide Polymorphism (SNP) markers from the alignment of Expressed Sequence Tag (EST) sequences in Atlantic cod. We also present basic population parameters of the SNPs in samples of North-East Arctic cod and Norwegian coastal cod obtained from three different localities, and test for SNPs that may have been targeted by natural selection. RESULTS: A total of 17,056 EST sequences were used to find 724 putative SNPs, from which 318 segregating SNPs were isolated. The SNPs were tested on Atlantic cod from four different sites, comprising both North-East Arctic cod (NEAC) and Norwegian coastal cod (NCC). The average heterozygosity of the SNPs was 0.25 and the average minor allele frequency was 0.18. FST values were highly variable, with the majority of SNPs displaying very little differentiation while others had FST values as high as 0.83. The FST values of 29 SNPs were found to be larger than expected under a strictly neutral model, suggesting that these loci are, or have been, influenced by natural selection. For the majority of these outlier SNPs, allele frequencies in a northern sample of NCC were intermediate between allele frequencies in a southern sample of NCC and a sample of NEAC, indicating a cline in allele frequencies similar to that found at the Pantophysin I locus. CONCLUSION: The SNP markers presented here are powerful tools for future genetics work related to management and aquaculture. In particular, some SNPs exhibiting high levels of population divergence have potential to significantly enhance studies on the population structure of Atlantic cod.

A New Single Nucleotide Polymorphism Database for Rainbow Trout Generated Through Whole Genome Resequencing.

Single-nucleotide polymorphisms (SNPs) are highly abundant markers, which are broadly distributed in animal genomes. For rainbow trout (Oncorhynchus mykiss), SNP discovery has been previously done through sequencing of restriction-site associated DNA (RAD) libraries, reduced representation libraries (RRL) and RNA sequencing. Recently we have performed high coverage whole genome resequencing with 61 unrelated samples, representing a wide range of rainbow trout and steelhead populations, with 49 new samples added to 12 aquaculture samples from AquaGen (Norway) that we previously used for SNP discovery. Of the 49 new samples, 11 were double-haploid lines from Washington State University (WSU) and 38 represented wild and hatchery populations from a wide range of geographic distribution and with divergent migratory phenotypes. We then mapped the sequences to the new rainbow trout reference genome assembly (GCA_002163495.1) which is based on the Swanson YY doubled haploid line. Variant calling was conducted with FreeBayes and SAMtools mpileup, followed by filtering of SNPs based on quality score, sequence complexity, read depth on the locus, and number of genotyped samples. Results from the two variant calling programs were compared and genotypes of the double haploid samples were used for detecting and filtering putative paralogous sequence variants (PSVs) and multi-sequence variants (MSVs). Overall, 30,302,087 SNPs were identified on the rainbow trout genome 29 chromosomes and 1,139,018 on unplaced scaffolds, with 4,042,723 SNPs having high minor allele frequency (MAF > 0.25). The average SNP density on the chromosomes was one SNP per 64 bp, or 15.6 SNPs per 1 kb. Results from the phylogenetic analysis that we conducted indicate that the SNP markers contain enough population-specific polymorphisms for recovering population relationships despite the small sample size used. Intra-Population polymorphism assessment revealed high level of polymorphism and heterozygosity within each population. We also provide functional annotation based on the genome position of each SNP and evaluate the use of clonal lines for filtering of PSVs and MSVs. These SNPs form a new database, which provides an important resource for a new high density SNP array design and for other SNP genotyping platforms used for genetic and genomics studies of this iconic salmonid fish species.

Myostatin (MSTN) gene duplications in Atlantic salmon (Salmo salar): evidence for different selective pressure on teleost MSTN-1 and -2.

Whereas the negative muscle regulator myostatin (MSTN) in mammals is almost exclusively expressed in the muscle by a single encoding gene, teleost fish possess at least two MSTN genes which are differentially expressed in both muscular and non-muscular tissues. Duplicated MSTN-1 genes have previously been identified in the tetraploid salmonid genome. From Atlantic salmon we succeeded in isolating the paralogous genes of MSTN-2, which shared about 70% identity with MSTN-1a and -1b. The salmon MSTN-2a cDNA encoded a predicted protein of 363 residues and included the conserved C-terminal bioactive domain. MSTN-2a seemed to be primarily expressed in the brain, and a functional role of teleost MSTN-2 in the neurogenesis similar to the inhibitory action of the closely related GDF-11 in the mammalian brain was proposed. In contrast, a frame-shift mutation in exon 1 of salmon MSTN-2b would lead to the synthesis of a putatively non-functional truncated protein. The absence of processed MSTN-2b mRNA in the examined tissues indicated that this gene has become a non-functional pseudogene. The differential, but partially overlapping, expression patterns of salmon MSTN-2a, -1a and -1b in muscular and non-muscular tissues are probably due to the different arrangement of the potential cis-acting regulatory elements identified in their putative promoter regions. Single and paired E-boxes in the MSTN-1b promoter were shown to bind both homo-and hetero-dimers of the myogenic regulatory factor MyoD and E47 in vitro of importance for initiating the myogenic program. Analyses of nucleotide substitution patterns indicated that the teleost MSTNs essentially have evolved under purifying selection, but a subset of amino acid sites under positive selective pressure were identified within the MSTN1 branch. The results may reflect the evolutionary forces related to adoption of the different functional roles proposed for the teleost MSTN isoforms. The phylogenetic analysis of multiple vertebrate MSTNs suggested at least two separate gene duplication events in the fish lineage. Linkage analysis of polymorphic microsatellites within intron 2 of salmon MSTN-1a and -1b mapped the two genes to different linkage groups in agreement with the tetraploid origin of the duplicated salmonid MSTN-1 and MSTN-2 genes.

Mapping of a quantitative trait locus for resistance against infectious salmon anaemia in Atlantic salmon (Salmo salar): comparing survival analysis with analysis on affected/resistant data.

BACKGROUND: Infectious Salmon Anaemia (ISA) is a viral disease affecting farmed Atlantic salmon (Salmo salar) worldwide. The identification of Quantitative Trait Loci (QTL) affecting resistance to the disease could improve our understanding of the genetics underlying the trait and provide a means for Marker-Assisted Selection. We previously performed a genome scan on commercial Atlantic salmon families challenge tested for ISA resistance, identifying several putative QTL. In the present study, we set out to validate the strongest of these QTL in a larger family material coming from the same challenge test, and to determine the position of the QTL by interval mapping. We also wanted to explore different ways of performing QTL analysis within a survival analysis framework (i.e. using time-to-event data), and to compare results using survival analysis with results from analysis on the dichotomous trait 'affected/resistant'. RESULTS: The QTL, located on Atlantic salmon linkage group 8 (following SALMAP notation), was confirmed in the new data set. Its most likely position was at a marker cluster containing markers BHMS130, BHMS170 and BHMS553. Significant segregation distortion was observed in the same region, but was shown to be unrelated to the QTL. A maximum likelihood procedure for identifying QTL, based on the Cox proportional hazard model, was developed. QTL mapping was also done using the Haley-Knott method (affected/resistant data), and within a variance-component framework (affected/resistant data and time-to-event data). In all cases, analysis using affected/resistant data gave stronger evidence for a QTL than did analysis using time-to-event data. CONCLUSION: A QTL for resistance to Infectious Salmon Anaemia in Atlantic salmon was validated in this study, and its more precise location on linkage group eight was determined. The QTL explained 6% of the phenotypic variation in resistance to the disease. The linkage group also displayed significant segregation distortion. Survival models proved in this case not to be more suitable than models based on the dichotomous trait 'affected/resistant' for analysing the data.