Identification and Genomic Localization of Autosomal sdY Locus in a Population of Atlantic Salmon (Salmo salar).

The determination of sex in salmonid fishes is controlled by genetic mechanisms, with males being the heterogametic sex. The master sex-determining gene, the sexually dimorphic gene on the Y chromosome (sdY), is a conserved gene across various salmonid species. Nevertheless, variations in the genomic location of sdY have been observed both within and between species. Furthermore, different studies have reported discordances in the association between the sdY and the phenotypic gender. While some males seem to lack this locus, there have been reports of females carrying sdY. Although the exact reasons behind this discordance remain under investigation, some recent studies have proposed the existence of an autosomal, non-functional copy of sdY as a potential cause. In this study, we confirmed the presence of this autosomal sdY in the SalmoBreed strain of Atlantic salmon using a genotyping platform through a novel approach that allows for high-throughput screening of a large number of individuals. We further characterized the segregation profile of this locus across families and found the ratio of genetically assigned female-to-male progeny to be in accordance with the expected profile of a single autosomal sdY locus. Additionally, our mapping efforts localized this locus to chromosome 3 and suggested a putative copy on chromosome 6.

Transcriptomic response to ISAV infection in the gills, head kidney and spleen of resistant and susceptible Atlantic salmon.

BACKGROUND: Infectious Salmon Anaemia virus (ISAV) is an orthomyxovirus responsible for large losses in Atlantic salmon (Salmo salar) aquaculture. Current available treatments and vaccines are not fully effective, and therefore selective breeding to produce ISAV-resistant strains of Atlantic salmon is a high priority for the industry. Genomic selection and potentially genome editing can be applied to enhance the disease resistance of aquaculture stocks, and both approaches can benefit from increased knowledge on the genomic mechanisms of resistance to ISAV. To improve our understanding of the mechanisms underlying resistance to ISAV in Atlantic salmon we performed a transcriptomic study in ISAV-infected salmon with contrasting levels of resistance to this virus. RESULTS: Three different tissues (gills, head kidney and spleen) were collected on 12 resistant and 12 susceptible fish at three timepoints (pre-challenge, 7 and 14 days post challenge) and RNA sequenced. The transcriptomes of infected and non-infected fish and of resistant and susceptible fish were compared at each timepoint. The results show that the responses to ISAV are organ-specific; an important response to the infection was observed in the head kidney, with up-regulation of immune processes such as interferon and NLR pathways, while in gills and spleen the response was more moderate. In addition to immune related genes, our results suggest that other processes such as ubiquitination and ribosomal processing are important during early infection with ISAV. Moreover, the comparison between resistant and susceptible fish has also highlighted some interesting genes related to ubiquitination, intracellular transport and the inflammasome. CONCLUSIONS: Atlantic salmon infection by ISAV revealed an organ-specific response, implying differential function during the infection. An immune response was observed in the head kidney in these early timepoints, while gills and spleen showed modest responses in comparison. Comparison between resistance and susceptible samples have highlighted genes of interest for further studies, for instance those related to ubiquitination or the inflammasome.

Use of DNA pools of a reference population for genomic selection of a binary trait in Atlantic salmon.

Genomic selection has a great potential in aquaculture breeding since many important traits are not directly measured on the candidates themselves. However, its implementation has been hindered by staggering genotyping costs because of many individual genotypes. In this study, we explored the potential of DNA pooling for creating a reference population as a tool for genomic selection of a binary trait. Two datasets from the SalmoBreed population challenged with salmonid alphavirus, which causes pancreas disease, were used. Dataset-1, that includes 855 individuals (478 survivors and 377 dead), was used to develop four DNA pool samples (i.e., 2 pools each for dead and survival). Dataset-2 includes 914 individuals (435 survivors and 479 dead) belonging to 65 full-sibling families and was used to develop in-silico DNA pools. SNP effects from the pool data were calculated based on allele frequencies estimated from the pools and used to calculate genomic breeding values (GEBVs). The correlation between SNP effects estimated based on individual genotypes and pooled data increased from 0.3 to 0.912 when the number of pools increased from 1 to 200. A similar trend was also observed for the correlation between GEBVs, which increased from 0.84 to 0.976, as the number of pools per phenotype increased from 1 to 200. For dataset-1, the accuracy of prediction was 0.71 and 0.70 when the DNA pools were sequenced in 40× and 20×, respectively, compared to an accuracy of 0.73 for the SNP chip genotypes. For dataset-2, the accuracy of prediction increased from 0.574 to 0.691 when the number of in-silico DNA pools increased from 1 to 200. For this dataset, the accuracy of prediction using individual genotypes was 0.712. A limited effect of sequencing depth on the correlation of GEBVs and prediction accuracy was observed. Results showed that a large number of pools are required to achieve as good prediction as individual genotypes; however, alternative effective pooling strategies should be studied to reduce the number of pools without reducing the prediction power. Nevertheless, it is demonstrated that pooling of a reference population can be used as a tool to optimize between cost and accuracy of selection.

Raman and near Infrared Spectroscopy for Quantification of Fatty Acids in Muscle Tissue-A Salmon Case Study.

The aim of the present study was to critically evaluate the potential of using NIR and Raman spectroscopy for prediction of fatty acid features and single fatty acids in salmon muscle. The study was based on 618 homogenized salmon muscle samples acquired from Atlantic salmon representing a one year-class nucleus, fed the same high fish oil feed. NIR and Raman spectra were used to make regression models for fatty acid features and single fatty acids measured by gas chromatography. The predictive performance of both NIR and Raman was good for most fatty acids, with R2 above 0.6. Overall, Raman performed marginally better than NIR, and since the Raman models generally required fewer components than respective NIR models to reach high and optimal performance, Raman is likely more robust for measuring fatty acids compared to NIR. The fatty acids of the salmon samples co-varied to a large extent, a feature that was exacerbated by the overlapping peaks in NIR and Raman spectra. Thus, the fatty acid related variation of the spectroscopic data of the present study can be explained by only a few independent principal components. For the Raman spectra, this variation was dominated by functional groups originating from long-chain polyunsaturated FAs like EPA and DHA. By exploring the independent EPA and DHA Raman models, spectral signatures similar to the respective pure fatty acids could be seen. This proves the potential of Raman spectroscopy for single fatty acid prediction in muscle tissue.

Identification of a New Infectious Pancreatic Necrosis Virus (IPNV) Variant in Atlantic Salmon (Salmo salar L.) that can Cause High Mortality Even in Genetically Resistant Fish.

Infectious pancreatic necrosis (IPN) is an important viral disease of salmonids that can affect fish during various life cycles. In Atlantic salmon, selecting for genetically resistant fish against IPN has been one of the most highly praised success stories in the history of fish breeding. During the late 2000s, the findings that resistance against this disease has a significant genetic component, which is mainly controlled by variations in a single gene, have helped to reduce the IPN outbreaks to a great extent. In this paper, we present the identification of a new variant of the IPN virus from a field outbreak in Western Norway that had caused mortality, even in genetically resistant salmon. We recovered and assembled the full-length genome of this virus, following the deep-sequencing of the head-kidney transcriptome. The comparative sequence analysis revealed that for the critical amino acid motifs, previously found to be associated with the degree of virulence, the newly identified variant is similar to the virus's avirulent form. However, we detected a set of deduced amino acid residues, particularly in the hypervariable domain of the VP2, that collectively are unique to this variant compared to all other reference sequences assessed in this study. We suggest that these mutations have likely equipped the virus with the capacity to escape the host defence mechanism more efficiently, even in the genetically deemed IPN resistant fish.

The heritable landscape of near-infrared and Raman spectroscopic measurements to improve lipid content in Atlantic salmon fillets.

BACKGROUND: Product quality and production efficiency of Atlantic salmon are, to a large extent, influenced by the deposition and depletion of lipid reserves. Fillet lipid content is a heritable trait and is unfavourably correlated with growth, thus genetic management of fillet lipid content is needed for sustained genetic progress in these two traits. The laboratory-based reference method for recording fillet lipid content is highly accurate and precise but, at the same time, expensive, time-consuming, and destructive. Here, we test the use of rapid and cheaper vibrational spectroscopy methods, namely near-infrared (NIR) and Raman spectroscopy both as individual phenotypes and phenotypic predictors of lipid content in Atlantic salmon. RESULTS: Remarkably, 827 of the 1500 individual Raman variables (i.e. Raman shifts) of the Raman spectrum were significantly heritable (heritability (h2) ranging from 0.15 to 0.65). Similarly, 407 of the 2696 NIR spectral landscape variables (i.e. wavelengths) were significantly heritable (h2 = 0.27-0.40). Both Raman and NIR spectral landscapes had significantly heritable regions, which are also informative in spectroscopic predictions of lipid content. Partial least square predicted lipid content using Raman and NIR spectra were highly concordant and highly genetically correlated with the lipid content values ([Formula: see text] = 0.91-0.98) obtained with the reference method using Lin's concordance correlation coefficient (CCC = 0.63-0.90), and were significantly heritable ([Formula: see text] = 0.52-0.67). CONCLUSIONS: Both NIR and Raman spectral landscapes show substantial additive genetic variation and are highly genetically correlated with the reference method. These findings lay down the foundation for rapid spectroscopic measurement of lipid content in salmonid breeding programmes.

Identification of genetic loci associated with higher resistance to pancreas disease (PD) in Atlantic salmon (Salmo salar L.).

BACKGROUND: Pancreas disease (PD) is a contagious disease caused by salmonid alphavirus (SAV) with significant economic and welfare impacts on salmon farming. Previous work has shown that higher resistance against PD has underlying additive genetic components and can potentially be improved through selective breeding. To better understand the genetic basis of PD resistance in Atlantic salmon, we challenged 4506 smolts from 296 families of the SalmoBreed strain. Fish were challenged through intraperitoneal injection with the most virulent form of the virus found in Norway (i.e., SAV3). Mortalities were recorded, and more than 900 fish were further genotyped on a 55 K SNP array. RESULTS: The estimated heritability for PD resistance was 0.41 ± 0.017. The genetic markers on two chromosomes, ssa03 and ssa07, showed significant associations with higher disease resistance. Collectively, markers on these two QTL regions explained about 60% of the additive genetic variance. We also sequenced and compared the cardiac transcriptomics of moribund fish and animals that survived the challenge with a focus on candidate genes within the chromosomal segments harbouring QTL. Approximately 200 genes, within the QTL regions, were found to be differentially expressed. Of particular interest, we identified various components of immunoglobulin-heavy-chain locus B (IGH-B) on ssa03 and immunoglobulin-light-chain on ssa07 with markedly higher levels of transcription in the resistant animals. These genes are closely linked to the most strongly QTL associated SNPs, making them likely candidates for further investigation. CONCLUSIONS: The findings presented here provide supporting evidence that breeding is an efficient tool for increasing PD resistance in Atlantic salmon populations. The estimated heritability is one of the largest reported for any disease resistance in this species, where the majority of the genetic variation is explained by two major QTL. The transcriptomic analysis has revealed the activation of essential components of the innate and the adaptive immune responses following infection with SAV3. Furthermore, the complementation of the genomic with the transcriptomic data has highlighted the possible critical role of the immunoglobulin loci in combating PD virus.

Genome-Wide Association Study Confirms Previous Findings of Major Loci Affecting Resistance to Piscine myocarditis virus in Atlantic Salmon (Salmo salar L.).

Cardiomyopathy syndrome is a viral disease of Atlantic salmon, mostly affecting fish during the late stages of production, resulting in significant losses to the industry. It has been shown that resistance to this disease has a strong genetic component, with quantitative trait loci (QTL) on chromosomes 27 (Ssa27) and Ssa12 to explain most of the additive genetic variance. Here, by analysing animals from a different year-class and a different population, we further aimed to confirm and narrow down the locations of these QTL. The data support the existence of the two QTL and suggest that the causative mutation on Ssa27 is most likely within the 10-10.5 Mbp segment of this chromosome. This region contains a cluster of major histocompatibility complex class I (MHC I) genes with the most strongly associated marker mapped to one of these loci. On Ssa12, the data confirmed the previous finding that the location of the causative mutation is within the 61.3 to 61.7 Mbp region. This segment contains several immune-related genes, but of particular interest are genes related to MHC II. Together, these findings highlight the likely key role of MHC genes in Atlantic salmon following infection with Piscine myocarditis virus (PMCV) and their potential impact on influencing the trajectory of this disease.

Genome-Wide Association Study of Piscine Myocarditis Virus (PMCV) Resistance in Atlantic Salmon (Salmo salar).

Cardiomyopathy syndrome is a severe, viral disease of Atlantic salmon that mostly affects farmed animals during their late production stage at sea. Caused by piscine myocarditis virus (PMCV), over the past few years outbreaks due to this disease have resulted in significant losses to the aquaculture industry. However, there is currently no vaccine that has proven effective against this virus. In this study, using a challenge model, we investigated the genetic variation for resistance to PMCV, by screening a large number of animals using a 55 K SNP array. In particular, we aimed to identify genetic markers that are tightly linked to higher disease resistance and can potentially be used in breeding programs. Using genomic information, we estimated a heritability of 0.51 ± 0.06, suggesting that resistance against this virus, to a great extent, is controlled by genetic factors. Through association analysis, we identified a significant quantitative trait locus (QTL) on chromosome 27, explaining approximately 57% of the total additive genetic variation. The region harboring this QTL contains various immune-related candidate genes, many of which have previously been shown to have a different expression profile between the naïve and infected animals. We also identified a suggestive association on chromosome 12, with the QTL linked markers located in 2 putatively immune-related genes. These results are of particular interest, as they can readily be implemented into breeding programs, can further assist in fine-mapping the causative mutations, and help in better understanding the biology of the disease and the immunological mechanisms underlying resistance against PMCV.

Genomic and Transcriptomic Analysis of Amoebic Gill Disease Resistance in Atlantic Salmon (Salmo salar L.).

Amoebic gill disease (AGD) is one of the most important parasitic diseases of farmed Atlantic salmon. It is a source of major economic loss to the industry and poses significant threats to animal welfare. Previous studies have shown that resistance against this disease has a moderate, heritable genetic component, although the genes and the genetic pathways that contribute to this process have yet to be elucidated. In this study, to identify the genetic mechanisms of AGD resistance, we first investigated the molecular signatures of AGD infection in Atlantic salmon through a challenge model, where we compared the transcriptome profiles of the naïve and infected animals. We then conducted a genome-wide association analysis with 1,333 challenged tested fish to map the AGD resistance genomic regions, supported by the results from the transcriptomic data. Further, we investigated the potential of incorporating gene expression analysis results in genomic prediction to improve prediction accuracy. Our data suggest thousands of genes have modified their expression following infection, with a significant increase in the transcription of genes with functional properties in cell adhesion and a sharp decline in the abundance of various components of the immune system genes. From the genome-wide association analysis, QTL regions on chromosomes ssa04, ssa09, and ssa13 were detected to be linked with AGD resistance. In particular, we found that QTL region on ssa04 harbors members of the cadherin gene family. These genes play a critical role in target recognition and cell adhesion. The QTL region on ssa09 also is associated with another member of the cadherin gene family, protocadherin Fat 4. The associated genetic markers on ssa13 span a large genomic region that includes interleukin-18-binding protein, a gene with function essential in inhibiting the proinflammatory effect of cytokine IL18. Incorporating gene expression information through a weighted genomic relationship matrix approach decreased genomic prediction accuracy and increased bias of prediction. Together, these findings help to improve our breeding programs and animal welfare against AGD and advance our knowledge of the genetic basis of host-pathogen interactions.

Individual differences in EPA and DHA content of Atlantic salmon are associated with gene expression of key metabolic processes.

The aim of this study was to explore how individual differences in content of the omega-3 fatty acids EPA and DHA in skeletal muscle of slaughter-sized Atlantic salmon, are associated with expression of genes involved in key metabolic processes. All experimental fish were fed the same diet throughout life and fasted for 14 days prior to slaughter. Still, there were relatively large individual variations in EPA and DHA content of skeletal muscle. Higher DHA content was concurrent with increased expression of genes of the glycolytic pathway and the production of pyruvate and lactate, whereas EPA was associated with increased expression of pentose phosphate pathway and glycogen breakdown genes. Furthermore, EPA, but not DHA, was associated with expression of genes involved in insulin signaling. Expression of genes specific for skeletal muscle function were positively associated with both EPA and DHA. EPA and DHA were also associated with expression of genes related to eicosanoid and resolvin production. EPA was negatively associated with expression of genes involved in lipid catabolism. Thus, a possible reason why some individuals have a higher level of EPA in the skeletal muscle is that they deposit - rather than oxidize - EPA for energy.

Genetic effects of fatty acid composition in muscle of Atlantic salmon.

BACKGROUND: The replacement of fish oil (FO) and fishmeal with plant ingredients in the diet of farmed Atlantic salmon has resulted in reduced levels of the health-promoting long-chain polyunsaturated omega-3 fatty acids (n-3 LC-PUFA) eicosapentaenoic (EPA; 20:5n-3) and docosahexaenoic acid (DHA; 22:6n-3) in their filets. Previous studies showed the potential of selective breeding to increase n-3 LC-PUFA levels in salmon tissues, but knowledge on the genetic parameters for individual muscle fatty acids (FA) and their relationships with other traits is still lacking. Thus, we estimated genetic parameters for muscle content of individual FA, and their relationships with lipid deposition traits, muscle pigmentation, sea lice and pancreas disease in slaughter-sized Atlantic salmon. Our aim was to evaluate the selection potential for increased n-3 LC-PUFA content and provide insight into FA metabolism in Atlantic salmon muscle. RESULTS: Among the n-3 PUFA, proportional contents of alpha-linolenic acid (ALA; 18:3n-3) and DHA had the highest heritability (0.26) and EPA the lowest (0.09). Genetic correlations of EPA and DHA proportions with muscle fat differed considerably, 0.60 and 0.01, respectively. The genetic correlation of DHA proportion with visceral fat was positive and high (0.61), whereas that of EPA proportion with lice density was negative. FA that are in close proximity along the bioconversion pathway showed positive correlations with each other, whereas the start (ALA) and end-point (DHA) of the pathway were negatively correlated (- 0.28), indicating active bioconversion of ALA to DHA in the muscle of fish fed high FO-diet. CONCLUSIONS: Since contents of individual FA in salmon muscle show additive genetic variation, changing FA composition by selective breeding is possible. Taken together, our results show that the heritabilities of individual n-3 LC-PUFA and their genetic correlations with other traits vary, which indicates that they play different roles in muscle lipid metabolism, and that proportional muscle contents of EPA and DHA are linked to body fat deposition. Thus, different selection strategies can be applied in order to increase the content of healthy omega-3 FAin the salmon muscle. We recommend selection for the proportion of EPA + DHA in the muscle because they are both essential FA and because such selection has no clear detrimental effects on other traits.