Strong selection pressures maintain divergence on genomic islands in Atlantic cod (Gadus morhua L.) populations.

BACKGROUND: Two distinct populations have been extensively studied in Atlantic cod (Gadus morhua L.): the Northeast Arctic cod (NEAC) population and the coastal cod (CC) population. The objectives of the current study were to identify genomic islands of divergence and to propose an approach to quantify the strength of selection pressures using whole-genome single nucleotide polymorphism (SNP) data. After applying filtering criteria, information on 93 animals (9 CC individuals, 50 NEAC animals and 34 CC × NEAC crossbred individuals) and 3,123,434 autosomal SNPs were used. RESULTS: Four genomic islands of divergence were identified on chromosomes 1, 2, 7 and 12, which were mapped accurately based on SNP data and which extended in size from 11 to 18 Mb. These regions differed considerably between the two populations although the differences in the rest of the genome were small due to considerable gene flow between the populations. The estimates of selection pressures showed that natural selection was substantially more important than genetic drift in shaping these genomic islands. Our data confirmed results from earlier publications that suggested that genomic islands are due to chromosomal rearrangements that are under strong selection and reduce recombination between rearranged and non-rearranged segments. CONCLUSIONS: Our findings further support the hypothesis that selection and reduced recombination in genomic islands may promote speciation between these two populations although their habitats overlap considerably and migrations occur between them.

QTLs Associated with Resistance to Cardiomyopathy Syndrome in Atlantic Salmon.

Cardiomyopathy syndrome (CMS) caused by piscine myocarditis virus is a major disease affecting the Norwegian Atlantic salmon industry. Three different populations of Atlantic salmon from the Mowi breeding program were used in this study. The first 2 populations (population 1 and 2) were naturally infected in a field outbreak, while the third population (population 3) went through a controlled challenged test. The aim of the study was to estimate the heritability, the genetic correlation between populations and perform genome-wide association analysis for resistance to this disease. Survival data from population 1 and 2 and heart atrium histology score data from population 3 was analyzed. A total of 571, 4312, and 901 fish from population 1, 2, and 3, respectively were genotyped with a noncommercial 55,735 Affymetrix marker panel. Genomic heritability ranged from 0.12 to 0.46 and the highest estimate was obtained from the challenge test dataset. The genetic correlation between populations was moderate (0.51-0.61). Two chromosomal regions (SSA27 and SSA12) contained single nucleotide polymorphisms associated with resistance to CMS. The highest association signal (P = 6.9751 × 10-27) was found on chromosome 27. Four genes with functional roles affecting viral resistance (magi1, pi4kb, bnip2, and ha1f) were found to map closely to the identified quantitative trait loci (QTLs). In conclusion, genetic variation for resistance to CMS was observed in all 3 populations. Two important quantitative trait loci were detected which together explain half of the total genetic variance, suggesting strong potential application for marker-assisted selection and genomic predictions to improve CMS resistance.

Training the salmon's genes: influence of aerobic exercise, swimming performance and selection on gene expression in Atlantic salmon.

BACKGROUND: Farmed and wild Atlantic salmon are exposed to many infectious and non-infectious challenges that can cause mortality when they enter the sea. Exercise before transfer promotes growth, health and survival in the sea. Swimming performance in juveniles at the freshwater parr stage is positively associated with resistance to some diseases. Genetic variation is likely to affect response to exercise. In this study we map genetic differences associated with aerobic exercise, swimming performance and genetic origin. Eggs from the selectively bred Bolaks salmon and wild Lærdal River salmon strains were reared until parr in a common environment. Swimming performance was assessed by subjecting the fish to either continuous hard exercise or control conditions for 18 days. Heart was sampled for examination of gene expression using RNA-seq (~60 fish/treatment). RESULTS: Lower expression of genes affecting immune function was found in domesticated than wild parr. Among wild parr under control exercise the expression of a large number of genes involved in general metabolism, stress and immune response was lower in superior swimmers suggesting that minimisation of energy expenditure during periods of low activity makes parr better able to sustain bursts of swimming for predator avoidance. A similar set of genes were down-regulated with training among wild parr with inferior swimming performance. These parr react to training in a way that their cardiac expression patterns become like the superior performing wild parr under control exercise conditions. Diversifying selection caused by breeding of domesticated stock, and adaptive pressures in wild stock, has affected the expression and frequency of single nucleotide polymorphisms (SNPs) for multiple functional groups of genes affecting diverse processes. SNPs associated with swimming performance in wild parr map to genes involved in energetic processes, coding for contractile filaments in the muscle and controlling cell proliferation. CONCLUSIONS: Domesticated parr have less phenotypic plasticity in response to training and lower expression of genes with functions affecting immune response. The genetic response to training is complex and depends on the background of parr and their swimming ability. Exercise should be tailored to the genetics and swimming performance of fish.

Genome-Wide Association Study for Identifying Loci that Affect Fillet Yield, Carcass, and Body Weight Traits in Rainbow Trout (Oncorhynchus mykiss).

Fillet yield (FY, %) is an economically-important trait in rainbow trout aquaculture that affects production efficiency. Despite that, FY has received little attention in breeding programs because it is difficult to measure on a large number of fish and cannot be directly measured on breeding candidates. The recent development of a high-density SNP array for rainbow trout has provided the needed tool for studying the underlying genetic architecture of this trait. A genome-wide association study (GWAS) was conducted for FY, body weight at 10 (BW10) and 13 (BW13) months post-hatching, head-off carcass weight (CAR), and fillet weight (FW) in a pedigreed rainbow trout population selectively bred for improved growth performance. The GWAS analysis was performed using the weighted single-step GBLUP method (wssGWAS). Phenotypic records of 1447 fish (1.5 kg at harvest) from 299 full-sib families in three successive generations, of which 875 fish from 196 full-sib families were genotyped, were used in the GWAS analysis. A total of 38,107 polymorphic SNPs were analyzed in a univariate model with hatch year and harvest group as fixed effects, harvest weight as a continuous covariate, and animal and common environment as random effects. A new linkage map was developed to create windows of 20 adjacent SNPs for use in the GWAS. The two windows with largest effect for FY and FW were located on chromosome Omy9 and explained only 1.0-1.5% of genetic variance, thus suggesting a polygenic architecture affected by multiple loci with small effects in this population. One window on Omy5 explained 1.4 and 1.0% of the genetic variance for BW10 and BW13, respectively. Three windows located on Omy27, Omy17, and Omy9 (same window detected for FY) explained 1.7, 1.7, and 1.0%, respectively, of genetic variance for CAR. Among the detected 100 SNPs, 55% were located directly in genes (intron and exons). Nucleotide sequences of intragenic SNPs were blasted to the Mus musculus genome to create a putative gene network. The network suggests that differences in the ability to maintain a proliferative and renewable population of myogenic precursor cells may affect variation in growth and fillet yield in rainbow trout.

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.

Two adjacent inversions maintain genomic differentiation between migratory and stationary ecotypes of Atlantic cod.

Atlantic cod is composed of multiple migratory and stationary populations widely distributed in the North Atlantic Ocean. The Northeast Arctic cod (NEAC) population in the Barents Sea undertakes annual spawning migrations to the northern Norwegian coast. Although spawning occurs sympatrically with the stationary Norwegian coastal cod (NCC), phenotypic and genetic differences between NEAC and NCC are maintained. In this study, we resolve the enigma by revealing the mechanisms underlying these differences. Extended linkage disequilibrium (LD) and population divergence were demonstrated in a 17.4-Mb region on linkage group 1 (LG1) based on genotypes of 494 SNPs from 192 parents of farmed families of NEAC, NCC or NEACxNCC crosses. Linkage analyses revealed two adjacent inversions within this region that repress meiotic recombination in NEACxNCC crosses. We identified a NEAC-specific haplotype consisting of 186 SNPs that was fixed in NEAC sampled from the Barents Sea, but segregating under Hardy-Weinberg equilibrium in eight NCC stocks. Comparative genomic analyses determine the NEAC configuration of the inversions to be the derived state and date it to ~1.6-2.0 Mya. The haplotype block harbours 763 genes, including candidates regulating swim bladder pressure, haem synthesis and skeletal muscle organization conferring adaptation to long-distance migrations and vertical movements down to large depths. Our results suggest that the migratory ecotype experiences strong directional selection for the two adjacent inversions on LG1. Despite interbreeding between NEAC and NCC, the inversions are maintaining genetic differentiation, and we hypothesize the co-occurrence of multiple adaptive alleles forming a 'supergene' in the NEAC population.

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.

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.

Analysis of immune-related ESTs and differential expression analysis of few important genes in lines of rohu (Labeo rohita) selected for resistance and susceptibility to Aeromonas hydrophila infection.

A total of 137,629 contigs generated via de novo transcriptome assembly from resistant and susceptible lines of rohu (first generation) raised against aeromoniasis were further analyzed in terms of defence-related genes. Out of 1,939 contigs showing homology to genes involved in immune processes, 1,866 were further categorised into different functional subgroups. Comparative analysis revealed five genes for the first time in any carp species out of which apolipoprotein h, septin 4 isoform 3 and septin isoform cra_c were identified for the first time in fish. Differential expression analysis of ten genes viz., heat shock proteins (Hsps) (Hsp30, Hsp70 and Hsp90), serum lectin isoform 1 (SLI1), linker histone H1M (LHH1M), NAD(P)H quinone 1 (NQO1), zona pellucida 2 (ZP2) and three unknown genes that were highly up-expressed in first generation resistant line fish from mRNA-seq coverage data, was carried out using susceptible and resistant individuals of the second generation selected populations in eight different tissues viz. liver, kidney, intestine, gill, brain, spleen, skin and muscle using qPCR. Significant up-regulation in Hsp90, NQO1, C_116914 and C_22454 in specific tissues of resistant line and variable expression in Hsp30 and LHH1M genes in different tissues of both lines were noticed. The expression of Hsp70 was lower in many tissues of the resistant line than in susceptible line rohu. The expression of ZP2, SLI1 and C_94589 genes was not significantly different in terms of fold difference between the two lines. Differentially expressed genes need further characterisation to explore their role in resistance to Aeromonas hydrophila infection in rohu.

QTL for white spot syndrome virus resistance and the sex-determining locus in the Indian black tiger shrimp (Penaeus monodon).

BACKGROUND: Shrimp culture is a fast growing aquaculture sector, but in recent years there has been a shift away from tiger shrimp Penaeus monodon to other species. This is largely due to the susceptibility of P. monodon to white spot syndrome virus disease (Whispovirus sp.) which has impacted production around the world. As female penaeid shrimp grow more rapidly than males, mono-sex production would be advantageous, however little is known about genes controlling or markers associated with sex determination in shrimp. In this study, a mapped set of 3959 transcribed single nucleotide polymorphisms were used to scan the P. monodon genome for loci associated with resistance to white-spot syndrome virus and sex in seven full-sibling tiger shrimp families challenged with white spot syndrome virus. RESULTS: Linkage groups 2, 3, 5, 6, 17, 18, 19, 22, 27 and 43 were found to contain quantitative trait loci significantly associated with hours of survival after white spot syndrome virus infection (P < 0.05 after Bonferroni correction). Nine QTL were significantly associated with hours of survival. Of the SNPs mapping to these and other regions with suggestive associations, many were found to occur in transcripts showing homology to genes with putative immune functions of interest, including genes affecting the action of the ubiquitin-proteasome pathway, lymphocyte-cell function, heat shock proteins, the TOLL pathway, protein kinase signal transduction pathways, mRNA binding proteins, lectins and genes affecting the development and differentiation of the immune system (eg. RUNT protein 1A). Several SNPs significantly associated with sex were mapped to linkage group 30, the strongest associations (P < 0.001 after Bonferroni correction) for 3 SNPs located in a 0.8 cM stretch between positions 43.5 and 44.3 cM where the feminisation gene (FEM-1, affecting sexual differentiation in Caenorhabditis elegans) mapped. CONCLUSIONS: The markers for disease resistance and sexual differentiation identified by this study could be useful for marker assisted selection to improve resistance to WSSV and for identifying homogametic female individuals for mono-sex (all female) production. The genes with putative functions affecting immunity and sexual differentiation that were found to closely map to these loci provide leads about the mechanisms affecting these important economic traits in shrimp.

A linkage map of transcribed single nucleotide polymorphisms in rohu (Labeo rohita) and QTL associated with resistance to Aeromonas hydrophila.

BACKGROUND: Production of carp dominates world aquaculture. More than 1.1 million tonnes of rohu carp, Labeo rohita (Hamilton), were produced in 2010. Aeromonas hydrophila is a bacterial pathogen causing aeromoniasis in rohu, and is a major problem for carp production worldwide. There is a need to better understand the genetic mechanisms affecting resistance to this disease, and to develop tools that can be used with selective breeding to improve resistance. Here we use a 6 K SNP array to genotype 21 full-sibling families of L. rohita that were experimentally challenged intra-peritoneally with a virulent strain of A. hydrophila to scan the genome for quantitative trait loci associated with disease resistance. RESULTS: In all, 3193 SNPs were found to be informative and were used to create a linkage map and to scan for QTL affecting resistance to A. hydrophila. The linkage map consisted of 25 linkage groups, corresponding to the number of haploid chromosomes in L. rohita. Male and female linkage maps were similar in terms of order, coverage (1384 and 1393 cM, respectively) and average interval distances (1.32 and 1.35 cM, respectively). Forty-one percent of the SNPs were annotated with gene identity using BLAST (cut off E-score of 0.001). Twenty-one SNPs mapping to ten linkage groups showed significant associations with the traits hours of survival and dead or alive (P <0.05 after Bonferroni correction). Of the SNPs showing significant or suggestive associations with the traits, several were homologous to genes of known immune function or were in close linkage to such genes. Genes of interest included heat shock proteins (70, 60, 105 and "small heat shock proteins"), mucin (5b precursor and 2), lectin (receptor and CD22), tributyltin-binding protein, major histocompatibility loci (I and II), complement protein component c7-1, perforin 1, ubiquitin (ligase, factor e4b isoform 2 and conjugation enzyme e2 c), proteasome subunit, T-cell antigen receptor and lymphocyte specific protein tyrosine kinase. CONCLUSIONS: A panel of markers has been identified that will be validated for use with both genomic and marker-assisted selection to improve resistance of L. rohita to A. hydrophila.

The development of a high density linkage map for black tiger shrimp (Penaeus monodon) based on cSNPs.

Transcriptome sequencing using Illumina RNA-seq was performed on populations of black tiger shrimp from India. Samples were collected from (i) four landing centres around the east coastline (EC) of India, (ii) survivors of a severe WSSV infection during pond culture (SUR) and (iii) the Andaman Islands (AI) in the Bay of Bengal. Equal quantities of purified total RNA from homogenates of hepatopancreas, muscle, nervous tissue, intestinal tract, heart, gonad, gills, pleopod and lymphoid organs were combined to create AI, EC and SUR pools for RNA sequencing. De novo transcriptome assembly resulted in 136,223 contigs (minimum size 100 base pairs, bp) with a total length 61 Mb, an average length of 446 bp and an average coverage of 163× across all pools. Approximately 16% of contigs were annotated with BLAST hit information and gene ontology annotations. A total of 473,620 putative SNPs/indels were identified. An Illumina iSelect genotyping array containing 6,000 SNPs was developed and used to genotype 1024 offspring belonging to seven full-sibling families. A total of 3959 SNPs were mapped to 44 linkage groups. The linkage groups consisted of between 16-129 and 13-130 markers, of length between 139-10.8 and 109.1-10.5 cM and with intervals averaging between 1.2 and 0.9 cM for the female and male maps respectively. The female map was 28% longer than the male map (4060 and 2917 cM respectively) with a 1.6 higher recombination rate observed for female compared to male meioses. This approach has substantially increased expressed sequence and DNA marker resources for tiger shrimp and is a useful resource for QTL mapping and association studies for evolutionarily and commercially important traits.

Characterization of the ceruloplasmin gene and its potential role as an indirect marker for selection to Aeromonas hydrophila resistance in rohu, Labeo rohita.

Ceruloplasmin is an acute phase protein found to be activated by the host immune system during stress conditions. The ceruloplasmin gene has been reported in several teleosts and here we characterize the gene and test its association with resistance to Aeromonas hydrophila in rohu, Labeo rohita. A ceruloplasmin mRNA sequence of 3355 base pairs (bp) was derived (GenBank ID: JX010736). The coding sequence (CDS) comprised of 3276 bp that coded for 1092 amino acids. Alignment results showed the greatest similarity with zebrafish followed by channel catfish sequence, and a phylogenetic tree constructed on the basis of amino acid sequences showed that rohu shares a common clade with these two species. In the ontogeny study, the expression of ceruloplasmin was detected at 9 h post-fertilization onwards, and a strong level of expression was detected at 24 h (38-fold) and 15 days (34-fold) post-fertilization. The ceruloplasmin transcripts were evident in liver, spleen, stomach and heart. Expression was undetectable in gill, brain, eye, skin, muscle, intestine, anterior and posterior kidney tissues. Expression of ceruloplasmin after A. hydrophila infection was up-regulated 6 h post-challenge and was modulated until 15 days post-challenge. The level of ceruloplasmin was also compared in rohu selectively bred for higher growth and disease resistance. The gene showed a 4.58-fold higher level of expression in resistant line over susceptible line rohu selected based on family challenge test survival to A. hydrophila. Serum ceruloplasmin levels in three year classes of rohu selected for higher growth showed a positive correlation (0.49 ± 1.11) with survival against challenge with A. hydrophila. The estimated heritability was also found to be quite high (0.50 ± 0.22) for this parameter. Thus, ceruloplasmin could be one of the useful marker traits for selection against A. hydrophila resistance in fish.

Expressed sequences and polymorphisms in rohu carp (Labeo rohita, Hamilton) revealed by mRNA-seq.

Expressed genes and polymorphisms were identified in lines of rohu Labeo rohita selected for resistance or susceptibility to Aeromonas hydrophila, an important bacterial pathogen causing aeromoniasis. All animals were grown in a common environment and RNA from ten individuals from each line pooled for Illumina mRNA-seq. De novo transcriptome assembly produced 137,629 contigs with 40× average coverage.Forty-four percent of the assembled sequences were annotated with gene names and ontology terms. Of these, 3,419 were assigned biological process terms related to "stress response" and 1,939 "immune system". Twenty-six contigs containing 38 single nucleotide polymorphisms (SNPs) were found to map to the Cyprinus carpio mitochondrial genome and over 26,000 putative SNPs and 1,700 microsatellite loci were detected. Seventeen percent of the 100 transcripts with coverage data most indicative of higher-fold expression(>5.6 fold) in the resistant line pool showed homology to major histocompatibility (MH), heat shock proteins (HSP)30, 70 and 90, glycoproteins or serum lectin genes with putative functions affecting immune response. Forty-one percent of these 100 transcripts showed no or low homology to known genes. Of the SNPs identified, 96 showing the highest allele frequency differences between susceptible and resistant line fish included transcripts with homology to MH class I and galactoside-binding soluble lectin, also with putative functions affecting innate and acquired immune response. A comprehensive sequence resource for L. rohita, including annotated microsatellites and SNPs from a mixture of A. hydrophila-susceptible and -resistant individuals, was created for subsequent experiments aiming to identify genes associated with A. hydrophila resistance.

Characterisation of a novel paralog of scavenger receptor class B member I (SCARB1) in Atlantic salmon (Salmo salar).

BACKGROUND: Red flesh colour is a unique trait found in some salmonid genera. Carotenoid pigments are not synthesized de novo in the fish, but are provided by dietary uptake. A better understanding of the molecular mechanisms underlying the cellular uptake and deposition of carotenoids could potentially be used to improve the low muscle deposition rate that is typically found in farmed Atlantic salmon. In addition, from an evolutionary point of view, the establishment and maintenance of this trait is still poorly understood. It has been demonstrated in several species that scavenger receptor class B, member 1 (SCARB1) is involved in intestinal absorption of carotenoids, which makes this gene a possible source of genetic variation in salmonid flesh pigmentation. RESULTS: In this study, a novel paralog of SCARB1 (SCARB1-2) was detected through screening for genetic variation in Atlantic salmon SCARB1. Full length SCARB1-2 encodes a protein with 89% identity to Atlantic salmon SCARB1, except for the C-terminal cytoplasmic tail that shows only 12% identity. The most prominent site of SCARB1 mRNA expression was in the mid gut, while a five-fold lower level was detected in Atlantic salmon skeletal muscle and liver. The SCARB1-2 mRNA was equally expressed in liver, muscle and mid gut, and at a lower level than SCARB1 mRNA. A total of seven different SCARB1-2 alleles comprising repetitive enhancer of zeste motifs (EZH2) were identified in the founding parents of a resource Atlantic salmon population. We mapped the SCARB1-2 paralog to a region on Atlantic salmon chromosome 1, containing a putative QTL for flesh colour. Addition of the SCARB1-2 marker increased the significance of this QTL, however the large confidence interval surrounding the QTL precludes confirmation of SCARB1-2 as a causative gene underlying variation in this trait. CONCLUSION: We have characterised a novel paralog of SCARB1 (SCARB1-2), have mapped it to Atlantic salmon chromosome 1 and have described its expression in various tissues. Mapping with SCARB1-2 alleles added further evidence for a QTL affecting flesh colour on this chromosome, however further studies are needed to confirm a functional role for this gene in flesh colour pigmentation.

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.

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.

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.