Difference in skin immune responses to infection with salmon louse (Lepeophtheirus salmonis) in Atlantic salmon (Salmo salar L.) of families selected for resistance and susceptibility.

Atlantic salmon is susceptible to the salmon louse (Lepeophtheirus salmonis) and the variation in susceptibility within the species can be exploited in selective breeding programs for louse resistant fish. In this study, lice counts were completed on 3000 siblings from 150 families of Atlantic salmon identified as high resistant (HR) and low resistant (LR) families in two independent challenge trials. Skin samples behind the dorsal fin (nearby lice attachment) were collected from ten extreme families (HR or LR) and analyzed by qPCR for the expression of 32 selected genes, including a number of genes involved in T helper cell (Th) mediated immune responses, which have been previously implied to play important roles during salmon louse infections. Most genes showed lower expression patterns in the LR than in HR fish, suggesting an immunosuppressed state in LR families. The average number of lice (chalimi) was 9 in HR and 15 in LR fish. Large variation in lice counts was seen both within resistant and susceptible families, which enabled us to subdivide the groups into HR < 10 and HR > 10, and LR < 10 and LR > 10 to better understand the effect of lice burden per se. As expected, expression patterns were influenced both by genetic background and the number of attached parasites. Higher number of lice (>10) negatively affected gene expression in both HR and LR families. In general, strongest down-regulation was seen in LR > 10 and lesser down-regulation in HR < 10. HR in general and especially HR < 10 fish were better at resisting suppression of expression of both Th1 and Th2 genes. However, the best inverse correlation with infection level was seen for the prototypical Th1 genes, including several members from the interferon pathways. In addition, skin histomorphometry suggests that infected LR salmon had thicker epidermis in the area behind the dorsal fin and larger mucous cell size compared to infected HR fish, however marginally significant (p = 0.08). This histomorphometric finding was in line with the immune response being skewed in LR towards the Th2 rather than a Th1 profile. Our findings suggest that the ability to resist lice infection depends on the ability to avoid immunosuppression and not as much on the physical tissue barrier functions.

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.

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.

Epithelial Cadherin Determines Resistance to Infectious Pancreatic Necrosis Virus in Atlantic Salmon.

Infectious pancreatic necrosis virus (IPNV) is the cause of one of the most prevalent diseases in farmed Atlantic salmon (Salmo salar). A quantitative trait locus (QTL) has been found to be responsible for most of the genetic variation in resistance to the virus. Here we describe how a linkage disequilibrium-based test for deducing the QTL allele was developed, and how it was used to produce IPN-resistant salmon, leading to a 75% decrease in the number of IPN outbreaks in the salmon farming industry. Furthermore, we describe how whole-genome sequencing of individuals with deduced QTL genotypes was used to map the QTL down to a region containing an epithelial cadherin (cdh1) gene. In a coimmunoprecipitation assay, the Cdh1 protein was found to bind to IPNV virions, strongly indicating that the protein is part of the machinery used by the virus for internalization. Immunofluorescence revealed that the virus colocalizes with IPNV in the endosomes of homozygous susceptible individuals but not in the endosomes of homozygous resistant individuals. A putative causal single nucleotide polymorphism was found within the full-length cdh1 gene, in phase with the QTL in all observed haplotypes except one; the absence of a single, all-explaining DNA polymorphism indicates that an additional causative polymorphism may contribute to the observed QTL genotype patterns. Cdh1 has earlier been shown to be necessary for the internalization of certain bacteria and fungi, but this is the first time the protein is implicated in internalization of a virus.

Genomic prediction in an admixed population of Atlantic salmon (Salmo salar).

Reliability of genomic selection (GS) models was tested in an admixed population of Atlantic salmon, originating from crossing of several wild subpopulations. The models included ordinary genomic BLUP models (GBLUP), using genome-wide SNP markers of varying densities (1-220 k), a genomic identity-by-descent model (IBD-GS), using linkage analysis of sparse genome-wide markers, as well as a classical pedigree-based model. Reliabilities of the models were compared through 5-fold cross-validation. The traits studied were salmon lice (Lepeophtheirus salmonis) resistance (LR), measured as (log) density on the skin and fillet color (FC), with respective estimated heritabilities of 0.14 and 0.43. All genomic models outperformed the classical pedigree-based model, for both traits and at all marker densities. However, the relative improvement differed considerably between traits, models and marker densities. For the highly heritable FC, the IBD-GS had similar reliability as GBLUP at high marker densities (>22 k). In contrast, for the lowly heritable LR, IBD-GS was clearly inferior to GBLUP, irrespective of marker density. Hence, GBLUP was robust to marker density for the lowly heritable LR, but sensitive to marker density for the highly heritable FC. We hypothesize that this phenomenon may be explained by historical admixture of different founder populations, expected to reduce short-range lice density (LD) and induce long-range LD. The relative importance of LD/relationship information is expected to decrease/increase with increasing heritability of the trait. Still, using the ordinary GBLUP, the typical long-range LD of an admixed population may be effectively captured by sparse markers, while efficient utilization of relationship information may require denser markers (e.g., 22 k or more).

How specific MHC class I and class II combinations affect disease resistance against infectious salmon anaemia in Atlantic salmon (Salmo salar).

The aim was to evaluate the performance of selected individual MHC class I and class II alpha (A) alleles, and combinations of these on disease resistance against infectious salmon anaemia (ISA). The material consisting of 1966 fish from seven families, contained five MHC class I alleles and four MHC class II A alleles. Which representing given class II A and class II beta (B) haplotypes, totalling 19 MHC class I and class II A genotypes. The fish were challenged with infectious salmon anaemia virus (ISAV), the virus causing ISA. Dead fish were collected daily during the challenge experiment and the survivors were collected at termination. All fish were genotyped for MHC class I and class II A. The total mortality in the material was 85.14%. For MHC class I, UBA*0201 and UBA*0301 were significantly the most resistant alleles, while UBA*0601 for class I and DAA*0301 for class II A were the significantly most susceptible alleles. The analysis of combined MHC class I and class II A genotypes detected that fish with the genotype UBA*0201/*0301;DAA*0201/*0201 were the most resistant fish with a hazard ratio (HR) at 0.750, while the fish with the genotypes UBA*0601/*0801;DAA*0501/*0501 and UBA*0201/*0301;DAA*0301/*0501 were the most susceptible fish with HR of 1.334 and 1.425. In addition, Cox regression analysis within family detected combined MHC class I and class II A genotypes that contributed significantly to resistance and susceptibility. The study confirmed the expectation of performance of individual MHC class I and class II A alleles, and also detected an effect of MHC class I and class II A in combinations.