Skip to main content
  • Research article
  • Open access
  • Published:

Epistatic regulation of growth in Atlantic salmon revealed: a QTL study performed on the domesticated-wild interface

Abstract

Background

Quantitative traits are typically considered to be under additive genetic control. Although there are indications that non-additive factors have the potential to contribute to trait variation, experimental demonstration remains scarce. Here, we investigated the genetic basis of growth in Atlantic salmon by exploiting the high level of genetic diversity and trait expression among domesticated, hybrid and wild populations.

Results

After rearing fish in common-garden experiments under aquaculture conditions, we performed a variance component analysis in four mapping populations totaling ~ 7000 individuals from six wild, two domesticated and three F1 wild/domesticated hybrid strains. Across the four independent datasets, genome-wide significant quantitative trait loci (QTLs) associated with weight and length were detected on a total of 18 chromosomes, reflecting the polygenic nature of growth. Significant QTLs correlated with both length and weight were detected on chromosomes 2, 6 and 9 in multiple datasets. Significantly, epistatic QTLs were detected in all datasets.

Discussion

The observed interactions demonstrated that the phenotypic effect of inheriting an allele deviated between half-sib families. Gene-by-gene interactions were also suggested, where the combined effect of two loci resulted in a genetic effect upon phenotypic variance, while no genetic effect was detected when the two loci were considered separately. To our knowledge, this is the first documentation of epistasis in a quantitative trait in Atlantic salmon. These novel results are of relevance for breeding programs, and for predicting the evolutionary consequences of domestication-introgression in wild populations.

Background

The process of domestication results in a set of genetic changes as a population is taken from the wild and bred over multiple generations in captivity [see 1]. It typically consists of a mixture of selective breeding for desired traits, inadvertent selection, relaxation of natural selection, and the stochastic process of genetic drift. Traditionally, directional selection in breeding programs was practised by phenotypic selection of individuals displaying a greater than average magnitude or frequency of the trait(s) of interest, working on the premise that at least part of the trait’s variance is heritable. As genomic resources have become more accessible, the most recent developments in selective breeding have utilized the statistical correlation between genotypes and phenotypes to predict phenotypic gain in the framework of an additive genetic model. This approach is commonly referred to as genomic selection (GS), and involves analysis of genome-distributed single nucleotide polymorphisms (SNPs) on the population under selection [2, 3].

By making use of large-scale Genome Wide Association (GWAS) studies, GS is often capable of detecting loci that have low or medium contribution to the trait (< 1%) [3, 4]. GS is also more precise than traditional phenotypic selection alone due to its higher ability at quantifying Mendelian sampling across siblings. In fact, GS can result in rapid improvements of livestock and crops [5,6,7]. However, due to its focus on heritability (i.e., additive genetic variance) GS does not necessarily use the full complexity of the genetic architecture of polygenic traits. A possible way to improve the performance of GS would therefore be to account for non-additive genetic variation [8, 9]. Polygenic traits are often simultaneously influenced by non-additive genetic mechanisms (dominance/epistasis) [10], and while heritability is the main target of focus to improve the performance of a given breed, non-additive genetic effects may play a key role in explaining phenotypic diversity. This is especially important when looking at fitness related traits (e.g., growth, shyness, foraging, predator awareness) in populations that are subjected to domestication, but still occasionally interbreed with their wild conspecifics (e.g., Atlantic salmon, Salmo salar L.). Deviating selection pressures in wild and captive environments are expected to favour very different genotypes leading to distinct phenotypes in the two environments [11]. In order to understand the genetic architecture shaping the fitness of wild and domesticated fish as well as their hybrids, it is important to consider all possible types of genetic mechanisms, including dominance and epistasis.

Quantitative trait loci (QTL) mapping in experimental intercrosses represents a useful tool to provide insights into the underlying genetic basis of the variability of the investigated trait(s). By focusing on an experimental pedigree with a high degree of relatedness between individuals, QTL studies focus on loci with large to medium effects and may not only reveal the number and distribution of loci contributing to the trait of interest [12], but also elucidate potential interactions between genes affecting the phenotype [13]. Epistasis, is often ignored in quantitative trait studies [13] and breeding programs, yet has the potential to significantly contribute to the phenotypic expression of the trait [10]. The benefit of including non-additive genetic effects in genomic prediction of complex traits is however not a resolved question [14, 15].

Atlantic salmon is one of the world’s most domesticated fishes [see 11]. Commercial aquaculture, including selective breeding programs, was first initiated in Norway in the early 1970’s [16]. Increasing growth rate through both phenotypic and pedigree-based selection has been the major target for all breeding programs [17, 18]. Atlantic salmon display a high heritability in growth rate, h2 > 30% [18], and the genetic gain per generation has been estimated to be ~ 10–15% of the breeding values [18, 19]. Furthermore, after ~ 12 generations of directional selection, genetic gains in growth rate are still being achieved each generation [see Fig. 4 in 11]. Therefore, after nearly half a century of domestication and selective breeding, farmed salmon now out-grow wild salmon several-fold when reared together under commercial conditions [20, 21]. Nevertheless, despite large and well-documented gains in growth rate through selection, the underlying genetic mechanisms remain largely elusive.

In salmon, as in most organisms, growth is regarded as a polygenic trait [22]. The importance of additive genetic inheritance on this trait is revealed by the large heritability. Nevertheless, there are indications that growth may also be influenced by non-additive factors [23]. Several growth-linked QTLs have been documented in Atlantic salmon, and QTLs on some of the same chromosomes have been reported across studies [24,25,26,27,28,29,30,31]. However, previous QTL-studies related to growth of Atlantic salmon have been restricted to domesticated salmon, with exceptions such as the study by Baranski et al. [24] which also included a landlocked population and the study by Besnier et al. [31] that focused on domesticated/wild interactions in a natural environment.

Atlantic salmon displays several key features making it ideal to investigate additive and non-additive genetic factors on growth [32]. Firstly, fast growing domesticated salmon originate from slow-growing wild salmon [16]. The latter clearly hold the genetic potential for elevated growth rates, although this is not selected for in the wild. Second, wild Atlantic salmon display substantial population genetic structure throughout their native range [33,34,35,36], including genetic-based phenotypic and life-history variation [37]. Third, genomic resources for the Atlantic salmon are now widely available [38]. Fourth, well-established rearing systems combined with the ability to cross genetically diverse wild and domesticated salmon provides a good experimental framework in which to investigate genetic inheritance of growth. By crossing fish from unrelated multiple wild populations and domesticated strains, the chances of including allelic combinations that may reveal the genetic architecture of a complex trait such as growth, and the potential for epistasis, is increased.

In this study, we conducted genome-wide QTL mapping on multiple domesticated, hybrid and wild populations of Atlantic salmon originating from geographically distinct regions (Fig. 1). Using a variance component mapping analysis, we studied ~ 7000 individuals originating from six wild, two domesticated and three F1 wild/domesticated hybrid strains, and identified genome-wide significant QTLs for freshwater growth. By creating a mapping population of diverse genetic background, we utilized the population genetic variance observed in Norwegian Atlantic salmon. This allowed us to search for replicable QTLs in four distinct datasets representing multiple strains, indicating their relevance for several of the Norwegian salmon populations, as well as detecting evidence of non-additive genetic architecture of the quantitative trait investigated.

Fig. 1
figure 1

Map of wild populations and location of the Matre research station. Wild Atlantic salmon, Salmo salar L., populations originating from six Norwegian rivers were included in the study. Parental salmon were collected directly from the rivers, except for salmon of the River Skibotn and Driva strain that had been conserved and reared in the Norwegian Gene Bank for Atlantic salmon. Parental salmon of the River Vosso strain had been reared by the gene bank until the smolt stage, and then released into the wild. All six wild populations, two domesticated strains and three F1 wild/domesticated strains were produced and reared at the Matre research station. The genome-wide quantitative trait locus (QTL) scan for traits related to freshwater growth were performed on a total of four datasets/experiment, including 134 families and ~ 7000 individuals. The map was produced by using the software QGIS 2.8 (https://www.qgis.org/en)

Results

Identification of QTLs within datasets

Genome scans for QTL detection were performed on four datasets consisting of different crosses between five wild populations and two domestic strains. The variation of fish weight and length between crosses is presented (Fig. 2).

Fig. 2
figure 2

Phenotypic growth measurements of Atlantic salmon of domesticated, hybrid and wild origin. Freshwater growth of Atlantic salmon, Salmo salar L., of all origins, in dataset 1–4, reared communally under hatchery conditions. Weight in grams on the left y-axis, length in cm on the right y-axis. The solid line illustrates the average weight of all fish in the respective datasets, while the dotted line illustrates the average length. Error bars show standard errors. Salmon in dataset 1 and 2 were sampled after their first winter as 1+ parr/smolt, while salmon in dataset 3 and 4 were sampled after their first summer as 0+ young of the year

By using similar statistical models in all datasets, i.e., not accounting for sex in dataset 1 and dataset 2, seven QTLs for weight were reported in dataset 1 on chromosomes 1, 2, 3, 7, 8, 9 and 19 (Fig. 3a), 12 QTLs in dataset 2 on chromosomes 3, 6, 7, 9, 10, 11, 12, 15, 20, 21, 22 and 28 (Fig. 3b), two QTLs in dataset 3 on chromosomes 2 and 17 (Fig. 3c), and five QTLs in dataset 4 on chromosomes 2, 6, 10, 21 and 23 (Fig. 3d) (Table 1). Each QTL accounted for 2.5–20% of the genetic variance for weight within each dataset (Table 1). By fitting all QTL effects in the same model, we estimated that the cumulated variance explained by all seven QTLs in datasets 1 to 4 accounted for 33.0, 36.2, 11.3 and 43.1%, respectively. When sex was accounted for, as a fixed covariate, the number of QTLs detected reduced to three in dataset 1 on chromosome 3, 9, and 19, and four in dataset 2, on chromosome 12, 15, 20, 22.

Fig. 3
figure 3

QTL scan for weight dataset 1 (3.a), dataset 2 (3.b), dataset 3 (3.c) and dataset 4 (3.d). Genome-wide scan at regular intervals (20–30 CM) of the Atlantic salmon linkage map for QTLs affecting freshwater weight, performed on the wild/domesticated interface. Horizontal lines indicate the 5 and 1% genome-wide significance threshold, based upon the likelihood ratio between a HGLM fitted at each genomic position with and without a QTL effect. Vertical lines separate chromosomes

Table 1 Weight QTLs

Also using similar models in all datasets, six QTLs for length were reported in dataset 1 on chromosomes 1, 7, 8, 9, 19 and 20 (Fig. 4a), eight QTLs in dataset 2 on chromosomes 2, 6, 9, 11, 15, 20, 22 and 28 (Fig. 4b), two QTLs in dataset 3 on chromosomes 2 and 17 (Fig. 4c), and five QTLs in dataset 4 on chromosomes 2, 6, 10, 21 and 23 (Fig. 4d) (Table 2). Each QTL accounted for 2.8–19.3% of the genetic variance for length, within each dataset (Table 2). By fitting all QTL effects in the same model, we estimated that the cumulated variance explained by all six QTLs in datasets 1 to 4 accounted for 23.0, 33.1, 12.5 and 40.6%, respectively. When sex was accounted for as a fix covariate, the number of QTLs detected was reduced to two in dataset 1 on chromosome 9 and 19, and one in dataset 2, on chromosome 2.

Fig. 4
figure 4

QTL scan for length dataset 1 (3.a), dataset 2 (3.b), dataset 3 (3.c) and dataset 4 (3.d). Genome-wide scan at regular intervals (20–30 CM) of the Atlantic salmon linkage map for QTLs affecting freshwater length, performed on the wild/domesticated interface. Horizontal lines indicate the 5 and 1% genome-wide significance threshold, based upon the likelihood ratio between a HGLM fitted at each genomic position with and without a QTL effect. Vertical lines separate chromosomes

Table 2 Length QTLs

In datasets 1 and 2, QTL scans were performed to detect genomic regions associated with phenotypic sex. In dataset 1, chromosome 2, 3 and 6 were significantly associated with sex, whereas in dataset 2, only chromosome 2 and 6 were associated with the phenotype.

Identification of QTLs across datasets

Several loci were consistently correlated with weight or length across datasets. QTLs located on chromosome 2 were correlated with weight in datasets 1, 3 and 4, and with length on datasets 2, 3 and 4. Similarly, QTLs on chromosome 6 were correlated with weight and length in datasets 2 and 4, while QTLs on chromosome 9 were correlated with weight and length in datasets 1 and 2. Detection of significant QTLs on the same chromosome in multiple datasets does not automatically demonstrate that the same polymorphic site is associated with the phenotypes(s) across datasets. However, between datasets, QTLs were detected at overlapping, but not always, identical positions (Tables 1 and 2).

Family contribution to QTL

The allelic substitution effect of parental alleles was estimated separately in each family, and the relative contribution of a wild versus domestic parent was calculated for each QTL. In all datasets, both wild and domesticated parents contributed to the genome-wide significant QTLs (Tables 1 and 2). A total of 107 wild and 73 domesticated parents were used in the four datasets. On average for the four datasets, the number of domesticated parental salmon was smaller than of the wild parents. This is reflected in the ratio of domestic versus wild alleles transmitted to the offspring generation. Domestic parents accounted for 40.6% of the allelic contribution to the F1 generation, while they accounted for 73.6 and 55.0% of the relative contribution towards the genome-wide significant QTLs for weight and length respectively (Tables 1 and 2). The proportion of domestic parents contributing to the F1 generation varied between the datasets: the domesticated parents accounted for 54.1, 29.1, 52.6, and 34.0% of the F1 alleles in datasets 1–4 respectively. In comparison, the domestic parents accounted for 76.3, 89.5, 51.0 and 40.8% of the relative contribution towards the genome-wide significant QTLs for weight and 60.5, 73.5, 36.0, and 26.6% of the relative contribution towards the genome-wide significant QTLs for length.

The parental contribution to phenotypic sex was also estimated in each family separately. This per-family scan for sex determination loci indicated that a locus associated with sex determination was present on chromosome 2 for 62 and 80% of the families from datasets 1 and 2 respectively. Loci associated with sex determination was also present on chromosome 3 for 19% of the families from dataset 1 and on chromosome 6 for 19 and 20% of the families from dataset 1 and dataset 2 respectively.

Epistasis

Two approaches were implemented to look for evidence of epistasis. First, the four datasets were screened for gene-by-gene interactions. Here, only suggestive evidence of interaction was detected which did not reach the genome-wide significant threshold. For example, in dataset 2, interaction was detected between the alleles on chromosome 20 and 24. The interaction occurred in the full sib families from parent 43 (domesticated) crossed with parent 15 (wild) and parent 36 (domesticated) where no genetic effect could be detected when considering either loci separately, while the combined effect of the two loci explained 10% of the phenotype variance in the family. Despite a nominal p value of 2.2.e− 06, the interaction between chromosomes 20 and 24 was not genome-wide significant (p ≈ 0.09) after conservative bootstrapping.

The second approach implemented to investigate epistasis was to estimate the interaction between allelic substitution effect and polygenic effect in a half-sib family context. This could only be performed when a given male was used to fertilize the egg of two distinct females (or a female having her eggs fertilized by two distinct males). In such case, genetic interaction is detectable when the allelic substitution effect changes significantly depending on which genetic background it is measured in.

From this approach, significant epistatic QTLs for weight and length were detected in all four datasets on a total of 9 chromosomes. For example, in dataset 1, female parent 7 (domesticated) was crossed with two male parents, 35 (wild) and 17 (domesticated), to produce two half-sib families. On chromosome 17, the allelic substitution effect was significantly linked with variation in the offspring length. Notably however, the one allele inherited from female parent 7 was associated with smaller offspring in the female 7 X male 35 family, yet the same allele inherited from female parent 7 was associated with larger offspring in the female 7 X male 17 family (Fig. 5a). Therefore, the direction of the allelic substitution effect changed between the offspring of male 35 and male 17.

Fig. 5
figure 5

Epistatic QTLs. Phenotypic response of epistatic QTLs affecting freshwater growth in half-sibling families of Atlantic salmon. Significant gene-by-parent interactions, detected in all four datasets, demonstrate that the effect at a given allele was affected by the genetic background of the parents in the half-sib family. Selected examples from all four datasets; a, dataset 1; b, dataset 2; c, dataset 3 and; d, dataset 4

A similar pattern to the above example was observed in dataset 2 on chromosome 9 for offspring of female parent 30 (wild), when crossed with male parent 39 (domesticated) and 18 (wild) (Fig. 5b); in dataset 3, on chromosome 4 for offspring of male parent 21 (wild), crossed with female parent 41 (domesticated), and 17 (wild) (Fig. 5c); and in dataset 4, on chromosome 2 for male parent 10 (wild), crossed with female parent 12 (wild) and 9 (wild) (Fig. 5d).

In total, significant gene-by-parent interactions were observed in a total of 9 half-sib families, on 9 different chromosomes, across all four datasets. These gene-by-parent interactions were responsible for between 2.5 to 16.4% of the phenotype variance within these families (Table 3), thus demonstrating a significant influence of non-additive variation on the trait.

Table 3 Epistatic QTLs

Discussion

In the present study, the genetic architecture underlying freshwater growth of Atlantic salmon was investigated through a genome-wide QTL scan performed on four independent datasets on the genetically divergent wild/domesticated interface. In addition to multiple QTLs, we identified multiple epistatic QTLs where the allelic substitution effect of a given locus changed depending on the genetic background it was measured in. For all epistatic QTLs, both maternally and paternally inherited, the phenotypic growth effect of inheriting a specific allele deviated between half-sib families, demonstrating that the effect of the allele was dependent on the parental genetic background. Epistasis has been defined as “the situation where the phenotype of a given genotype cannot be predicted by the sum of its component single-locus effects” [13]. Therefore, the present study provides the first experimental evidence of epistasis affecting phenotypic trait expression in one of the world’s most studied and socio-economically important fishes, the Atlantic salmon.

Detection of QTLs and parental contribution

This study was designed to exploit the large degree of genetic divergence among domesticated and wild salmon populations [11] in order to provide knowledge regarding the genetic basis of inheritance of the investigated quantitative trait, and to identify genomic regions affecting freshwater growth of Atlantic salmon populations in general. Significant QTLs on chromosome 2, affecting freshwater growth, were detected in all four datasets. QTLs affecting growth have previously been documented on chromosome 2 in domesticated salmon of varying age, both during freshwater [27,28,29] and marine rearing [24,25,26,27,28,29]. The latter studies included several domesticated strains originating from both European Atlantic salmon [24,25,26], North American Atlantic salmon [28] as well as trans-Atlantic backcrosses [27, 29] and a landlocked population [24]. Significant QTLs on chromosome 2 affecting growth have also been documented in domesticated and wild European Atlantic salmon studied in the wild [31]. Although these results strongly suggest the presence of QTLs on chromosome 2 linked to both juvenile and adult growth in Atlantic salmon, it is worth mentioning that a QTL on chromosome 2 is not reported in all studies [30], nor at all sampling points over time, when the same mapping material is sampled repeatedly [26]. The genetic background of the mapping population, the life stage investigated or the size (growth rate) of the fish could influence deviations in the detection of significant QTLs on chromosome 2 between these studies. Deviating results among studies could also reflect gene-by-environment interactions. In fact, deviating QTLs for juvenile growth in the same mapping material of salmon families reared in two separate environments, i.e., siblings reared in a hatchery and in the wild, has been documented [39]. We also detected genome-wide significant QTLs related to weight and/or length in two or more datasets on chromosomes 3, 6, 7, 9, 10, 20 and 21. Significant or suggestive QTLs linked to growth on these chromosomes have all previously been documented across the above-mentioned QTLs studies of Atlantic salmon. Growth is a highly polygenic trait, thus the detection of significant QTLs on several chromosomes here was expected. Likewise, the four independent datasets were obtained from the crosses of different genetically divergent salmon populations. It is therefore not surprising that we observed different QTLs across datasets.

Scanning the genome also showed that salmon chromosomes 2, 3 and 6 were linked with phenotypic sex, which concurs with previous publications on sex determination in salmon [40, 41]. Phenotypic sex was included as covariate in the search for growth QTLs in datasets 1 and 2, where fish were kept long enough after termination of the experiment to reach sexual maturity and thus be phenotypically determined. In both datasets 1 and 2, the number of significant QTLs dropped strongly after including sex as covariate in the model. However, that the number of fish included in the analysis dropped from 2000 to 1170 in dataset 1 and from 2400 to 1564 in dataset 2, when including phenotypic sex, is likely to have contributed to the reduction in the number of QTLs detected (i.e., reduced statistical power). For the QTL on chromosome 6, that was both canceled by the inclusion of phenotypic sex in the model, and linked with sex determination, it is not possible to exclude the possibility that the initially detected QTL reflects the effect of sex on growth rather than an independent locus affecting growth.

In the present study, parental salmon from the two domesticated strains and the six wild populations all contributed to the significant QTLs for freshwater weight and length. In total, domesticated parents accounted for 73.6 and 55.0% of the relative contribution towards the genome-wide significant QTLs linked to body weight and length respectively, while they contributed to 40.6% of the allelic composition of the F1. However, variations in parental contribution towards significant QTLs between parents of domesticated and wild origins were detected in the four datasets. Relative to the allelic contribution, domesticated parents contributed both more and less than expected in the different datasets. Differences in their contribution towards QTLs for weight and length were also detected. Parents contributing to the epistatic QTLs were also of both wild (60%) and domesticated origin (40%). The overall result suggests low allelic fixation for the loci contributing to the growth in the domesticated parental salmon included here. Lower allelic variation at highly polymorphic microsatellites has been reported on domesticated salmon strains [42, 43]. Furthermore, in dataset 3, reduced genetic variation for growth in the domesticated salmon strain had previously been suggested due to the documentation of lower heritability estimates for this trait in the domesticated relative to wild salmon [see 44]. Here, the domesticated parents contributed to the significant QTLs for length in a lower number than they contributed with gametes, but as expected for QTLs linked to growth.

The detection of similar, even higher levels of domesticated parental contribution to the overall QTL phenotypic and genotypic variation from growth, as revealed in this study, demonstrates that selection for growth, even after more than 12 generations of directional selection for this trait, has not reached saturation. This is in agreement with results from a suite of comparative growth studies between salmon of domesticated and wild parentage across the different generation of domestication [see Fig. 4 in 11]. Thus, it appears that despite the several-fold increase of growth rate in domesticated salmon, they still contain considerable growth potential that can be utilised for further genetic gains via directional selection.

Evidence of epistasis and its practical implications

Quantitative genetics and QTL mapping mainly focus on detecting loci that contribute additively to the phenotypic trait variation [45]. However, genetic interactions, such as epistasis and dominance, may also influence the phenotypic trait variation [13]. Here, we detected a non-additive genetic architecture of the trait investigated, i.e., growth. Notably, epistatic QTLs were observed in all four datasets included in the study. As the phenotypic growth effect of inheriting a specific allele was different between half-sib families, the parental genetic background affected the phenotypic expression of these alleles. Due to the opposing effect of inheriting a specific allele, the overall effect across the families might be evened out, and therefore a genome-wide significant QTL may or may not have been detected on that particular chromosome. We also found suggestions of gene-by-gene interactions, where the combining effect of two loci resulted in a genetic effect upon phenotypic variance, while no genetic effect was detected when the two loci were considered separately.

Non-additive inheritance in gene expressions has previously been documented [46,47,48,49,50]. However, and to the best of our knowledge, this is the first empirical study to reveal epistatic regulation of a quantitative trait in Atlantic salmon. Epistasis-influences on growth have been documented in other animals such as mice [10], chicken [51] and pigs [52]. The documentation of epistasis as a genetic basis for quantitative trait variation in Atlantic salmon could have implications for selective breeding programs; a topic that has been debated in both MAS [53] and GS [54]. The inclusion of non-additive effects have the potential to improve the accuracy of the predicted genetic values [54], however, the actual benefit of including non-additive genetic effects in genomic prediction is not a resolved question [14, 15].

Epistasis in domesticated, wild and hybrid half-sibling families may have implications for the long term consequences of genetic interactions between domesticated escapees and wild conspecifics [11]. Introgression of domesticated salmon has been detected in multiple wild populations in Norway and elsewhere [55,56,57,58]. In turn, this has also influenced life history traits, such as age and size at maturation [59]. Gene-flow from domesticated escapees is thought to occur mainly through domesticated females spawning with wild males [60], similar to the mating design used in dataset 1 (Fig. 5a) where we also documented epistasis. Individual-based eco-genetic models developed to study and quantify responses in wild populations to challenges such as genetic introgression, e.g., IBSEM: an individual-based Atlantic salmon population model [61, 62], are based on additive inheritance of traits. Although this approach is defendable given that additive genetic variation explains most of the genetic differences in survival in the wild between domesticated and wild salmon [31, 63,64,65,66,67], as well as growth under controlled conditions [20, 21, 44], the demonstrated possibility of epistasis, as revealed here, suggests that non-additive variation should not be overlooked.

Potential limitations

The present study is based on four datasets which all consist in two generations of intercrosses between wild and domestic Salmon. While this type of experimental design is powerful to detect QTLs, it also has a number of limitations. Due to the small recombination rate in Salmon, the parents and offspring share haplotype blocks consisting of large chromosome segments. This makes it possible to follow the parent to offspring allelic transmission with only a small number of genetic markers. It also conveys very little precision regarding the actual location the causative polymorphism associated with the trait. In short, the detected QTLs often cover half a chromosome or more, and it is therefore impossible to determine whether the association between QTL and phenotype is due to one or several causative polymorphic sites. In the case of epistatic QTL, the nature of the observed interaction is not accessible either. The change of allelic substitution effect between two half-sib families can be in fact due to the change of effect of a single site, or the average change of effect between several loci.

Conclusion

QTLs associated with weight and length were detected on 18 chromosomes, while three of these were consistent across multiple datasets, this indicates that these results are relevant for a wide range of salmon populations. Significantly, we had multiple observations that the effect of several QTL alleles changed between half-sib families, indicating epistatic regulation of growth (Table 3, Fig. 5a-d). To our knowledge, this is the first documentation of epistasis in a quantitative trait in Atlantic salmon. These novel results are of relevance for breeding programs, and for predicting the evolutionary consequences of domestication-introgression in wild populations.

Methods

Overall experimental design and phenotypic measurements

A genome-wide QTL scan for traits related to freshwater growth (weight and length) was performed on wild, domesticated, and F1 wild-domesticated hybrid Atlantic salmon from four independent experiments (datasets 1–4). In total, ~ 7000 salmon from 134 families, representing six wild populations, two domesticated strains and three wild-domesticated hybrid strains were included (Fig. 1). Fish within each dataset were communally-reared under standard fish farming conditions from the eyed-egg stage onwards. Growth measurements were collected from all offspring during the freshwater stage, either after their first summer as 0+ young of the year or after their first winter as 1+ parr/smolt. Prior to sampling, all individuals were anesthetised or euthanised with an overdose of metacain (Finquel® Vet, ScanVacc, Årnes, Norway). Fish were then wet weighed, fork length measured, and adipose or caudal fin clipped. Fins were preserved in 95% ethanol. All ~ 7000 individuals were originally genotyped with microsatellites markers [68,69,70,71,72] in order to perform parental assignment, and later with SNPs in order to link genetic variation with phenotypic variation. In addition, fish from datasets 1 and 2 were kept in the experimental facility until the adult stage, and phenotypic sex was recorded for those individuals surviving until sexual maturity. This was not performed for fish in datasets 3 or 4.

Experimental populations

The six wild populations included in this study originate from rivers spread along the coastline of Norway (Fig. 1), encompassing both of the identified main phylogenetic groups revealed in Norway [73, 74].

For the rivers Figgjo, Arna and Etne, wild parental salmon were caught in the river, and scale samples were analysed to verify that broodfish to be used in the experiments were wild salmon and not escaped domesticated salmon from fish-farms [75]. The salmon populations in the rivers Driva and Skibotn are conserved by the Norwegian Gene bank for Atlantic salmon, and parental salmon had been reared in freshwater at Haukvik, central Norway, for between one and three generations. In the gene bank, maintenance of wild salmon and their offspring are performed without any form of directional selection, although inadvertent selection may occur. The Vosso strain is also conserved by the Norwegian Gene bank; however, at the smolt stage fish are released in the wild and only returning salmon have been used as parents in this study. For more details on the Norwegian Gene Bank program for Atlantic salmon, see [76].

The commercial Mowi strain owned by Marine Harvest AS (recently renamed Mowi AS) is the oldest Norwegian domesticated strain [16]. This strain was established in 1969 when large multi-sea winter fish were collected from the River Bolstad in the Vosso watercourse and the River Åroy, in addition to wild salmon caught at sea off Western Norway near Osterfjord and Sotra [77, 78]. Individuals from the 10th and 11th generation, i.e., both overlapping and non-overlapping year classes, were used asparents to generate the biological material (offspring) in this study. The SalmoBreed domesticated strain was commercially established in 1999 but is based upon genetic material from several Norwegian domesticated strains that have been under commercial selection since the early 1970’s. Individuals from the approximately 11th generation were used as selected parents in this study. Fin tissue samples were collected from all parental salmon to extract DNA for parental assignment and QTL-mapping.

Dataset 1

In 2010, adult salmon from Figgjo and Mowi were used to generate 9 families of wild, 10 families of domesticated and 10 families of F1 hybrid (domesticated ♀ x wild ♂) origin. Thus, the F1 hybrids were maternal and paternal half-siblings to the domesticated and wild salmon respectively. Two replicates of 100 individuals per family (2900 individuals/replicate) were communally reared and later split into four replicates due to increasing biomass (c. 1450 individuals /replicate). In March 2012, 500 smolt/replicate (2000 smolt in total) were sampled for growth measurements (mean ± sd; weight (g): 62.2 ± 32.3, length (cm): 16.6. ±3.0, Fig. 2). One thousand one hundred seventy of these individuals reached the adult stage within the experimental facility, and had their phenotypic sex accurately recorded. All sampled individuals were thereafter genotyped and included in this study. More information about production and rearing of these groups can be found elsewhere [21]. (Additional file 1: Table S1).

Dataset 2

In 2011, adult salmon from Figgjo, Arna, Vosso and Mowi were used to generate 20 families of wild, 6 families of domesticated and 15 families of F1 hybrid (reciprocal Figgjo x Mowi) origin. Thus, the F1 hybrids represented both maternal and paternal half-siblings to the domesticated and wild salmon respectively. Two replicates of 50 individuals per family (2050 individuals/replicate) were communally reared and later split into four replicates due to increasing biomass (c. 1025 individuals /replicate). In March 2013, 600 smolts per replicate (2400 smolt in total), and 71 parr were sampled for growth measurements (mean ± sd; weight (g): 103.1 ± 72.9, length (cm): 19.4 ± 4.3, Fig. 2). One thousand five hundred sixty-four of these individuals reached the adult stage within the experimental facility, and had their phenotypic sex accurately recorded. All sampled individuals were thereafter genotyped and included in this study. More information about production of these experimental populations can be found elsewhere [79]. (Additional file 1: Table S2).

Dataset 3

In 2009, adult salmon from Etne and Mowi were used to generate 10 families of wild, 10 families of domesticated and 9 families of hybrid (domesticated ♀ x wild ♂) origin. Thus, the F1 hybrids were maternal and paternal half-siblings to the domesticated and wild salmon respectively. Two replicates of 50 individuals per family were communally reared until September 2010, when 750 individuals/replicate were sampled for growth measurements (mean ± sd; weight (g): 22.6 ± 12.0, length (cm): 11.6. ±2.2, Fig. 2). One thousand one hundred twenty-eight of the sampled individuals (564 individuals/replicate) were genotyped and included in this study. All individuals were terminated after the experiment, therefore, phenotypic sex was not recorded. More information about production and rearing of these groups can be found elsewhere [44]. (Additional file 1: Table S3).

Dataset 4

In 2013, adult salmon from Figgjo, Arna, Vosso, Driva, Skibotn, Mowi and SalmoBreed were used to generate 19 families of wild, 8 families of domesticated and 8 families of F1 hybrid origin. Both domesticated strains will from here on be referred to as Farm1 and Farm2 (random order). In total, two wild-domesticated F1 hybrid strains were produced by crossing the wild Figgjo and Vosso population with the two domesticated strains. Thus, the F1 hybrids were paternal and maternal half-siblings to one of the domesticated strains and the Figgjo strain, or maternal and paternal half-siblings to one of the domesticated strains and the wild Vosso strain respectively. Two replicates of 30 individuals per family (1040 individuals/replicate) were communally reared until September 2014 when 700 fry per replicate were sampled for growth measurements (mean ± sd; weight (g): 32.4 ± 13.5, length (cm): 13.5 ± 1.9, Fig. 2). All sampled individuals were genotyped and included in this study. All individuals were terminated after the experiment, therefore, phenotypic sex was not recorded. More information about production and rearing of these groups can be found elsewhere [20]. (Additional file 1: Table S4).

Genotyping and parentage assignment

For the QTL analysis of fish from all four datasets described above, an initial set of 116 genome-wide SNP markers were selected. These were selected from the genome to optimize genomic information content for QTL mapping and was based on both the marker position on the salmon genome and allelic frequencies in the parental generation. SNPs were selected for providing genotype information at regular intervals of 20–30 cm in the female recombination map, and for being polymorphic within full sib families. Where possible, SNPs displaying heterozygous parental genotypes in each family were prioratised. Previously, these SNPs have been shown to provide the best information content using some of the same strains and populations as in the present study [31]. SNP genotyping was performed on a MassARRAY Analyzer 4 from Agena Bioscience™, according to the manufacturer’s instructions. From the initial 116 SNPs, we produced 4 multiplexes containing a total of 114 SNPs (32, 32, 29 and 21 SNPs/ multiplex). A limited number of the SNPs did not amplify, leaving the final mapping dataset to include 109 genome-wide distributed SNPs.

Parentage testing of offspring from all four datasets was initially performed using FAP v3.6 [80], and six microsatellites. Following the exclusion-based approach implemented in FAP, 97 to 99% of the offspring were unambiguously assigned to their family of origin. After SNP genotyping, the microsatellite parental assignment was double-checked using the 109 SNP markers and a custom R script [R Core Team 81]. No discrepancies in parentage assignment were detected between the marker classes, and all fish not unambiguously assigned using microsatellites were subsequently assigned using SNPs. Therefore, all individuals were used in the analysis.

QTL mapping

Variance components

In each of the four datasets, the mapping population consisted of two generations (parents and offspring) from a total of 29 to 41 full and half-sib families. In each family, each offspring allele originated from one of the four parental haplotypes: two maternal and two paternal haplotypes. In order to perform QTL mapping, we first reconstructed the haplotypes of both parents and offspring based on pedigree and genotype data [82]. Compared to raw genotype data, the haplotyped data contains additional information regarding the parental (maternal or paternal) origin of each offspring allele, and identifies parental alleles linked within the same gametic haplotype. The next step consisted of estimating the Identity By Descent (IBD) coefficient between each pair of individuals at each locus along the genome. IBD coefficients were obtained from a recursive approach [83] implemented to account for haplotype information as input. The IBD coefficients are therefore estimated from the combined information at several markers. The QTL scan was thereafter performed by fitting a Mixed Linear Model at each genomic location as:

$$ y= X\beta + Ga+ Zq+e $$
(1)

where y is the phenotype vector, X the design matrix for fixed effects, β the vector of fixed effects, G the kinship matrix, a the vector of normal-distributed random polygenic effects, Z the design matrix for allelic effects, q the vector of normal-distributed random QTL effects, and e the normal-distributed random residuals. Note that GG’ is equivalent to the square kinship matrix, and covariance structure for the random polygenic effects, and ZZ’ is equivalent to the square IBD matrix and covariance structure of the random QTL effects. The fixed part of the model () consisted of the effect of strain and replicate. When phenotypic sex was available, sex was also incorporated in the model as fixed covariate. The Mixed model was fitted with the R package HGLM [84].

At each tested genomic position, the likelihood of model 1 is compared to the likelihood of the model without QTL effect:

$$ y= X\beta + Ga+e $$
(2)

In both models, we consider the adjusted profile log-likelihood profiled over random effects as provided by HGLM [84]. The likelihood ratio between model 1 and model 2 is then considered as the indicator for QTL i.e., correlation between genotype and phenotype variance. To account for multiple testing along the genome, the genome wide significance threshold for likelihood ratio was obtain through permutation test as in GA Churchill and RW Doerge [85].

After the initial genome scan, all significant QTL were fitted into the same model (model3) in order to estimate the proportion of genetic variance explained by all QTL simultaneously.

$$ y= X\beta + Ga+{\sum}_{i=1}^n{Z}_i{q}_i+e $$
(3)

Where Zi is the design matrix of QTL effect at QTL i, qi the vector of random QTL effects at QTL i, and n is the number of significant QTL detected in the dataset.

Family-based fixed effect model

To identify the parental alleles associated with phenotypic variation, a simpler linear model was fitted with fixed genetic effects. The model was applied on the offspring of each parent separately:

$$ y= Zq+e $$
(4)

where z is a two-column design matrix indicating for each offspring, the probability of having inherited either of the two parental alleles at a given locus, q is the allelic substitution effect at each locus i.e., the average difference between the phenotype of the offspring that inherit allele 1 or allele 2 from the same parent. When a given parent had offspring with two different mates or more, a fixed mate affect was added to the model, as well as an interaction term between the allelic substitution effect and the mate effect, as follows:

$$ y= X\beta + Zq+ ZXp+e $$
(5)

where X is the design matrix for the mate effect a, β the mate effect, and p the interaction term.

Epistasis

Investigation of non-additive genetic effects was also performed in order to detect possible gene-by-gene or gene-by-parent interactions. In the case of gene-by-gene interaction, a two loci model was fitted for each locus pair across the genome:

$$ y= X\beta +{Z}_1q+{Z}_2p+e $$
(6)
$$ y= X\beta +{Z}_1q+{Z}_2p+{Z}_1{Z}_2r+e $$
(7)

where q and p are the allelic effects at locus 1 and locus 2 respectively, and Z1 and Z2 are the corresponding design matrices, and r is the vector of the interaction terms between effects of locus 1 and locus 2. The magnitude of the interaction effect between locus 1 and locus 2 was tested by likelihood ratio between model 5 and model 6. In a similar fashion, interactions between parents and genotypes were tested by comparing the likelihood of model 5 and model 8 that include an interaction term between allelic effect and mate:

$$ y= X\beta + Zq+ XZr+e $$
(8)

Availability of data and materials

All data analysed during this study are included in this published article as supplementary files (File S1 to S4).

Abbreviations

GS:

Genomic Selection

GWAS:

Genome Wide Association Study

HGLM:

Hierarchical Generalized Linear Model

MAS:

Marker Assisted Selection

QTL:

Quantitative Trait Locus (loci)

SNP:

Single Nucleotide Polymorphism

References

  1. Diamond J. Evolution, consequences and future of plant and animal domestication. Nature. 2002;418(6898):700–7.

    Article  PubMed  Google Scholar 

  2. Goddard ME, Hayes BJ. Genomic selection. J Anim Breed Genet. 2007;124(6):323–30.

    Article  PubMed  Google Scholar 

  3. Meuwissen THE, Hayes BJ, Goddard ME. Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001;157(4):1819–29.

    PubMed  PubMed Central  CAS  Google Scholar 

  4. Powell JE, Kranis A, Floyd J, Dekkers JCM, Knott S, Haley CS. Optimal use of regression models in genome-wide association studies. Anim Genet. 2012;43(2):133–43.

    Article  PubMed  Google Scholar 

  5. Heffner EL, Jannink JL, Sorrells ME. Genomic selection accuracy using multifamily prediction models in a wheat breeding program. Plant Genome. 2011;4(1):65–75.

    Article  Google Scholar 

  6. Schaeffer LR. Strategy for applying genome-wide selection in dairy cattle. J Anim Breed Genet. 2006;123(4):218–23.

    Article  PubMed  Google Scholar 

  7. Nielsen HM, Sonesson AK, Yazdi H, Meuwissen THE. Comparison of accuracy of genome-wide and BLUP breeding value estimates in sib based aquaculture breeding schemes. Aquaculture. 2009;289(3–4):259–64.

    Article  Google Scholar 

  8. Hemani G, Knott S, Haley C. An evolutionary perspective on epistasis and the missing heritability. PLoS Genet. 2013;9(2):e1003295.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Robinson MR, Wray NR, Visscher PM. Explaining additional genetic variation in complex traits. Trends Genet. 2014;30(4):124–32.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Routman EJ, Cheverud JM. Gene effects on a quantitative trait: two-locus epistatic effects measured at microsatellite markers and at estimated QTL. Evolution. 1997;51(5):1654–62.

    Article  PubMed  Google Scholar 

  11. Glover K, Solber MF, McGinnity P, Hindar K, Verspoor E, Coulson MW, Hansen MM, Araki H, Skaala Ø, Svåsand T. Half a century of genetic interaction between farmed and wild Atlantic salmon: status of knowledge and unanswered questions. Fish Fish. 2017;18(5).

    Article  Google Scholar 

  12. Yue GH. Recent advances of genome mapping and marker-assisted selection in aquaculture. Fish Fish. 2014;15(3):376–96.

    Article  Google Scholar 

  13. Carlborg O, Haley CS. Epistasis: too often neglected in complex trait studies? Nat Rev Genet. 2004;5(8):618–U614.

    Article  PubMed  Google Scholar 

  14. Morgante F, Huang W, Maltecca C, Mackay TFC. Effect of genetic architecture on the prediction accuracy of quantitative traits in samples of unrelated individuals. Heredity. 2018;120(6):500–14.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Mäki-Tanila A, Hill WG. Influence of gene interaction on complex trait variation with multilocus models. Genetics. 2014;198(1):355–67.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Gjedrem T, Gjøen HM, Gjerde B. Genetic-origin of Norwegian farmed Atlantic salmon. Aquaculture. 1991;98(1–3):41–50.

    Article  Google Scholar 

  17. Gjedrem T. The first family-based breeding program in aquaculture. Rev Aquac. 2010;2(1):2–15.

    Article  Google Scholar 

  18. Gjedrem T. Genetic improvement of cold-water fish species. Aquac Res. 2000;31(1):25–33.

    Article  Google Scholar 

  19. Thodesen J, Gjedrem T. Breeding programs on Atlantic salmon in Norway: lessons learned. In: Development of aquatic animal genetic improvement and dissemination programs: current status and action plans. Edited by Ponzoni RW, Acosta BO, Ponniah AG; 2006:22–26.

  20. Harvey CA, Glover KA, Taylor MI, Creer S, Carvalho GR. A common garden design reveals population-specific variability in potential impacts of hybridization between populations of farmed and wild Atlantic salmon, Salmo salar L. Evol Appl. 2016;9:435–49.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Solberg MF, Zhang Z, Nilsen F, Glover KA. Growth reaction norms of domesticated, wild and hybrid Atlantic salmon families in response to differing social and physical environments. BMC Evol Biol. 2013;13:234.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Bentsen HB. Genetic effects of selection on polygenic traits with examples from Atlantic salmon, Salmo salar L. Aquac Fish Manag. 1994;25(1):89–102.

    Google Scholar 

  23. Debes PV, Fraser DJ, Yates M, Hutchings JA. The between-population genetic architecture of growth, maturation, and plasticity in atlantic salmon. Genetics. 2014;196(4):1277.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Baranski M, Moen T, Vage DI. Mapping of quantitative trait loci for flesh colour and growth traits in Atlantic salmon (Salmo salar). Genet Sel Evol. 2010;42:17.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Houston RD, Bishop SC, Hamilton A, Guy DR, Tinch AE, Taggart JB, Derayat A, McAndrew BJ, Haley CS. Detection of QTL affecting harvest traits in a commercial Atlantic salmon population. Anim Genet. 2009;40(5):753–5.

    Article  PubMed  Google Scholar 

  26. Gutierrez AP, Lubieniecki KP, Davidson EA, Lien S, Kent MP, Fukui S, Withler RE, Swift B, Davidson WS. Genetic mapping of quantitative trait loci (QTL) for body-weight in Atlantic salmon (Salmo salar) using a 6.5 K SNP array. Aquaculture. 2012;358:61–70.

    Article  Google Scholar 

  27. Boulding EG, Culling M, Glebe B, Berg PR, Lien S, Moen T. Conservation genomics of Atlantic salmon: SNPs associated with QTLs for adaptive traits in parr from four trans-Atlantic backcrosses. Heredity. 2008;101(4):381–91.

    Article  PubMed  Google Scholar 

  28. Reid DP, Szanto A, Glebe B, Danzmann RG, Ferguson MM. QTL for body weight and condition factor in Atlantic salmon (Salmo salar): comparative analysis with rainbow trout (Oncorhynchus mykiss) and Arctic charr (Salvelinus alpinus). Heredity. 2005;94(2):166–72.

    Article  PubMed  Google Scholar 

  29. Pedersen S, Berg PR, Culling M, Danzmann RG, Glebe B, Leadbeater S, Lien S, Moen T, Vandersteen W, Boulding EG. Quantitative trait loci for precocious parr maturation, early smoltification, and adult maturation in double-backcrossed trans-Atlantic salmon (Salmo salar). Aquaculture. 2013;410:164–71.

    Article  Google Scholar 

  30. Tsai HY, Hamilton A, Guy DR, Tinch AE, Bishop SC, Houston RD. The genetic architecture of growth and fillet traits in farmed Atlantic salmon (Salmo salar). BMC Genet. 2015;16.

  31. Besnier F, Glover KA, Lien S, Kent M, Hansen MM, Shen X, Skaala O. Identification of quantitative genetic components of fitness variation in farmed, hybrid and native salmon in the wild. Heredity. 2015;115(1):47–55.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Mackay TFC. Epistasis and quantitative traits: using model organisms to study gene-gene interactions. Nat Rev Genet. 2014;15(1):22–33.

    Article  PubMed  Google Scholar 

  33. Dionne M, Caron F, Dodson JJ, Bernatchez L. Landscape genetics and hierarchical genetic structure in Atlantic salmon: the interaction of gene flow and local adaptation. Mol Ecol. 2008;17(10):2382–96.

    Article  PubMed  Google Scholar 

  34. Ståhl G. Genetic population structure of Atlantic salmon. In: Ryman N, editor. Population Genetics and fishery management. Seattle: University of Washington Press; 1987. p. 121–36.

    Google Scholar 

  35. Verspoor E, Beardmore JA, Consuegra S, García de Leániz C, Hindar K, Jordan WC, Koljonen M-L, Mahkrov AA, Paaver T, Sánchez JA, et al. Population structure in the Atlantic salmon: insights from 40 years of research into genetic protein variation. J Fish Biol. 2005;67:3–54.

    Article  Google Scholar 

  36. King TL, Verspoor E, Spidle AP, Gross R, Philips RB, Koljonen M-L, Sánchez JA, Morrison CL. Biodiversity and population structure. In: Verspoor E, Stradmeyer L, Nielsen J, editors. The Atlantic salmon: genetics, conservation and management. Oxford: Blackwell Publishing Ltd; 2007. p. 115–66.

    Google Scholar 

  37. García de Leániz C, Fleming IA, Einum S, Verspoor E, Jordan WC, Consuegra S, Aubin-Horth N, Lajus D, Letcher BH, Youngson AF, et al. A critical review of adaptive genetic variation in Atlantic salmon: implications for conservation. Biol Rev. 2007;82(2):173–211.

    Article  PubMed  Google Scholar 

  38. Lien S, Gidskehaug L, Moen T, Hayes BJ, Berg PR, Davidson WS, Omholt SW, Kent MP. A dense SNP-based linkage map for Atlantic salmon (Salmo salar) reveals extended chromosome homeologies and striking differences in sex-specific recombination patterns. BMC Genomics. 2011;12:615.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Vasemagi A, Kahar S, Ozerov MY. Genes that affect Atlantic salmon growth in hatchery do not have the same effect in the wild. Funct Ecol. 2016;30(10):1687–95.

    Article  Google Scholar 

  40. Kijas J, McWilliam S, Naval Sanchez M, Kube P, King H, Evans B, Nome T, Lien S, Verbyla K. Evolution of sex determination loci in Atlantic Salmon. Sci Rep. 2018;8(1):5664.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Eisbrenner WD, Botwright N, Cook M, Davidson EA, Dominik S, Elliott NG, Henshall J, Jones SL, Kube PD, Lubieniecki KP, et al. Evidence for multiple sex-determining loci in Tasmanian Atlantic salmon (Salmo salar). Heredity (Edinb). 2014;113(1):86–92.

    Article  Google Scholar 

  42. Norris AT, Bradley DG, Cunningham EP. Microsatellite genetic variation between and within farmed and wild Atlantic salmon (Salmo salar) populations. Aquaculture. 1999;180(3–4):247–64.

    Article  Google Scholar 

  43. Skaala O, Hoyheim B, Glover K, Dahle G. Microsatellite analysis in domesticated and wild Atlantic salmon (Salmo salar L.): allelic diversity and identification of individuals. Aquaculture. 2004;240(1–4):131–43.

    Article  Google Scholar 

  44. Solberg MF, Skaala Ø, Nilsen F, Glover KA. Does domestication cause changes in growth reaction norms? A study of farmed, wild and hybrid Atlantic salmon families exposed to environmental stress. PLoS One. 2013;8(1):e54469.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Barton NH, Keightley PD. Understanding quantitative genetic variation. Nat Rev Genet. 2002;3(1):11–21.

    Article  PubMed  Google Scholar 

  46. Bicskei B, Taggart JB, Glover K, Bron JE. Comparing the transcriptomes of embryos from domesticated and wild Atlantic salmon (Salmo salar L.) stocks and examining factors that influence heritability of gene expression. Genet Sel Evol. 2016;48:20.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Bicskei B, Bron JE, Glover KA, Taggart JB. A comparison of gene transcription profiles of domesticated and wild Atlantic salmon (Salmo salar L.) at early life stages, reared under controlled conditions. BMC Genomics. 2014;15:884.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Debes PV, Normandeau E, Fraser DJ, Bernatchez L, Hutchings JA. Differences in transcription levels among wild, domesticated, and hybrid Atlantic salmon (Salmo salar) from two environments. Mol Ecol. 2012;21:2574–87.

    Article  PubMed  Google Scholar 

  49. Normandeau E, Hutchings JA, Fraser DJ, Bernatchez L. Population-specific gene expression responses to hybridization between farm and wild Atlantic salmon. Evol Appl. 2009;2(4):489–503.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Roberge C, Normandeau E, Einum S, Guderley H, Bernatchez L. Genetic consequences of interbreeding between farmed and wild Atlantic salmon: insights from the transcriptome. Mol Ecol. 2008;17(1):314–24.

    Article  PubMed  Google Scholar 

  51. Carlborg O, Kerje S, Schutz K, Jacobsson L, Jensen P, Andersson L. A global search reveals epistatic interaction between QTL for early growth in the chicken. Genome Res. 2003;13(3):413–21.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Crooks L, Guo Y: Consequences of epistasis on growth in an erhualian x white duroc pig cross. PLoS One 2017, 12(1):e0162045.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Liu PY, Zhu J, Lou XY, Lu Y. A method for marker-assisted selection based on QTLs with epistatic effects. Genetica. 2003;119(1):75–86.

    Article  PubMed  Google Scholar 

  54. Wang D, El-Basyoni IS, Baenziger PS, Crossa J, Eskridge KM, Dweikat I. Prediction of genetic values of quantitative traits with epistatic effects in plant breeding populations. Heredity. 2012;109(5):313–9.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Glover KA, Pertoldi C, Besnier F, Wennevik V, Kent M, Skaala Ø. Atlantic salmon populations invaded by farmed escapees: quantifying genetic introgression with a Bayesian approach and SNPs. BMC Genet. 2013;14:74.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Karlsson S, Diserud OH, Fiske P, Hindar K. Widespread genetic introgression of escaped farmed Atlantic salmon in wild salmon populations. ICES J Mar Sci. 2016;73(10):2488–98.

    Article  Google Scholar 

  57. Clifford SL, McGinnity P, Ferguson A. Genetic changes in an Atlantic salmon population resulting from escaped juvenile farm salmon. J Fish Biol. 1998;52(1):118–27.

    Article  Google Scholar 

  58. Keyser F, Wringe BF, Jeffery NW, Dempson JB, Duffy S, Bradbury IR. Predicting the impacts of escaped farmed Atlantic salmon on wild salmon populations. Can J Fish Aquat Sci. 2018;75:1–7.

    Article  Google Scholar 

  59. Bolstad GH, Hindar K, Robertsen G, Jonsson B, Sægrov H, Diserud OH, Fiske P, Jensen AJ, Urdal K, Næsje TF, et al. Gene flow from domesticated escapes alters the life history of wild Atlantic salmon. 2017;1:0124.

  60. Fleming IA, Jonsson B, Gross MR, Lamberg A. An experimental study of the reproductive behaviour and success of farmed and wild Atlantic salmon (Salmo salar). J Appl Ecol. 1996;33(4):893–905.

    Article  Google Scholar 

  61. Castellani M, Heino M, Gilbey J, Araki H, Svasand T, Glover KA. IBSEM: an individual-based atlantic salmon population model. PLoS One. 2015;10(9):e0138444.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Castellani M, Heino M, Gilbey J, Araki H, Svåsand T, Glover K. Modeling fitness changes in wild Atlantic salmon populations faced by spawning intrusion of domesticated escapees. Evol Appl. 2018;11:1–16.

    Article  Google Scholar 

  63. McGinnity P, Stone C, Taggart JB, Cooke D, Cotter D, Hynes R, McCamley C, Cross TF, Ferguson A. Genetic impact of escaped farmed Atlantic salmon (Salmo salar L.) on native populations: use of DNA profiling to assess freshwater performance of wild, farmed, and hybrid progeny in a natural river environment. ICES J Mar Sci. 1997;54(6):998–1008.

    Google Scholar 

  64. McGinnity P, Prodohl P, Ferguson K, Hynes R, O’Maoileidigh N, Baker N, Cotter D, O'Hea B, Cooke D, Rogan G, et al. Fitness reduction and potential extinction of wild populations of Atlantic salmon, Salmo salar, as a result of interactions with escaped farm salmon. Proc R Soc Lond B. 2003;270(1532):2443–50.

    Article  Google Scholar 

  65. Reed TE, Prodöhl P, Hynes R, Cross T, Ferguson A, McGinnity P. Quantifying heritable variation infitness-related traits of wild, farmed and hybrid Atlantic salmon families in a wild river environment. Heredity. 2015;115(2):173–84.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Skaala Ø, Glover KA, Barlaup BT, Svåsand T, Besnier F, Hansen MM, Borgstrøm R. Performance of farmed, hybrid and wild Atlantic salmon (Salmo salar) families in a natural river environment. Can J Fish Aquat Sci. 2012;69:1–13.

    Article  Google Scholar 

  67. Fleming IA, Hindar K, Mjølnerød IB, Jonsson B, Balstad T, Lamberg A. Lifetime success and interactions of farm salmon invading a native population. Proc R Soc Lond B. 2000;267(1452):1517–23.

    Article  Google Scholar 

  68. Grimholt U, Drablos F, Jørgensen SM, Hoyheim B, Stet RJM. The major histocompatibility class I locus in Atlantic salmon (Salmo salar L.): polymorphism, linkage analysis and protein modelling. Immunogenetics. 2002;54(8):570–81.

    Article  PubMed  Google Scholar 

  69. Sánchez JA, Clabby C, Ramos D, Blanco G, Flavin F, Vazquez E, Powell R. Protein and microsatellite single locus variability in Salmo salar L. (Atlantic salmon). Heredity. 1996;77:423–32.

    Article  PubMed  Google Scholar 

  70. O'Reilly PT, Hamilton LC, McConnell SK, Wright JM. Rapid analysis of genetic variation in Atlantic salmon (Salmo salar) by PCR multiplexing of dinucleotide and tetranucleotide microsatellites. Can J Fish Aquat Sci. 1996;53(10):2292–8.

    CAS  Google Scholar 

  71. Stet RJM, de Vries B, Mudde K, Hermsen T, van Heerwaarden J, Shum BP, Grimholt U. Unique haplotypes of co-segregating major histocompatibility class II a and class II B alleles in Atlantic salmon (Salmo salar) give rise to diverse class II genotypes. Immunogenetics. 2002;54(5):320–31.

    Article  PubMed  Google Scholar 

  72. Slettan A, Olsaker I, Lie O. Atlantic salmon, Salmo salar, microsattelites at the SsOSL25, SsOSL85, SsOSL311, SsOSL417 loci. Anim Genet. 1995;26(4):281–2.

    Article  PubMed  Google Scholar 

  73. Bourret V, Kent MP, Primmer CR, Vasemagi A, Karlsson S, Hindar K, McGinnity P, Verspoor E, Bernatchez L, Lien S. SNP-array reveals genome-wide patterns of geographical and potential adaptive divergence across the natural range of Atlantic salmon (Salmo salar). Mol Ecol. 2013;22(3):532–51.

    Article  PubMed  Google Scholar 

  74. Glover KA, Quintela M, Wennevik V, Besnier F, Sørvik AGE, Skaala Ø. Three decades of farmed escapees in the wild: a spatio-temporal analysis of Atlantic salmon population genetic structure throughout Norway. PLoS One. 2012;7(8):e43129.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Lund RA, Hansen LP, Jarvi T. Identifisering av oppdrettslaks og vill-laks ved ytre morfologi, finnestørrelse og skjellkarakterer. Res Rep. 1989;1:1–54 (in Norwegian) from the Norwegian Institute for Nature Research.

    Google Scholar 

  76. Gausen D. The Norwegian gene bank programme for Atlantic salmon (Salmo salar). In: Cloud JG, Thorgaard GH, editors. Genetic conservation of salmonid fishes, vol. 248. Boston: Springer; 1993.

    Google Scholar 

  77. Glover KA, Otterå H, Olsen RE, Slinde E, Taranger GL, Skaala Ø. A comparison of farmed, wild and hybrid Atlantic salmon (Salmo salar L.) reared under farming conditions. Aquaculture. 2009;286(3–4):203–10.

    Article  Google Scholar 

  78. Skaala Ø, Taggart JB, Gunnes K. Genetic differences between five major domesticated strains of Atlantic salmon and wild salmon. J Fish Biol. 2005;67:118–28.

    Article  Google Scholar 

  79. Solberg MF, Fjelldal PG, Nilsen F, Glover KA. Hatching time and alevin growth prior to the onset of exogenous feeding in farmed, wild and hybrid Norwegian Atlantic salmon. PLoS One. 2014;9(12):e113697.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Taggart JB. FAP: an exclusion-based parental assignment program with enhanced predictive functions. Mol Ecol Notes. 2007;7(3):412–5.

    Article  Google Scholar 

  81. R: a language and environment for statistical computing [http://www.r-project.org/].

  82. Besnier F, Carlborg O. A genetic algorithm based method for stringent haplotyping of family data. BMC Genet. 2009;10:57.

    Article  PubMed  PubMed Central  Google Scholar 

  83. Pong-Wong R, George AW, Woolliams JA, Haley CS. A simple and rapid method for calculating identity-by-descent matrices using multiple markers. Genet Sel Evol. 2001;33(5):453–71.

    Article  PubMed  PubMed Central  Google Scholar 

  84. Alam M, Roennegard L, Shen X. Fitting conditional and simultaneous autoregressive spatial models in hglm. R Journal. 2015;7(2):5–18.

    Article  Google Scholar 

  85. Churchill GA, Doerge RW. Empirical threshold values for quantitative trait mapping. Genetics. 1994;138(3):963–71.

    PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We greatly acknowledge Lise Dyrhovden, Ivar Helge Matre, Kåre Storsæter and Jan Olav Fosse at the Matre research station (IMR) for rearing of fish included in this study, Anne Grete E. Sørvik and Zhiwei Zhang for microsatellite and SNP genotyping, and Marine Harvest (Mowi) and SalmoBreed for generously and unconditionally providing the domesticated gametes for the experimental crosses.

Funding

This study was funded by resources from the Research Council of Norway via the projects INTERACT and QUANTESCAPE (respective grant numbers 200510 and 216105), and via the European Community’s Seventh Framework Programme for the project AQUATRACE (grant number KBBE-311920). The funding bodies had no influence on the design, acquisition, analysis, interpretation, nor presentation of the data and results of this study.

Author information

Authors and Affiliations

Authors

Contributions

Conceived and designed the study: All authors. Conducted the experiments (fish rearing and sampling): MFS, ACH, ØS, KAG. Conducted genotyping and QC of genetic data: GD, FB. Analyzed the data: FB, MFS. Wrote the first draft of the paper: MFS, FB, KAG. Contributed to final draft of paper: All authors. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Francois Besnier.

Ethics declarations

Ethics and consent to participate

The experimental protocols (permit numbers 4268, 5296, 2648) were approved by the Norwegian Animal Research Authority (NARA). Welfare and use of experimental animals were performed in strict accordance with the Norwegian Animal Welfare Act. This included anesthesia or euthanasia of fish using metacain (Finquel® Vet, ScanVacc, Årnes, Norway), during all samples referred to. In addition, all personnel involved in this experiment had undergone training approved by the Norwegian Food Safety Authority, which is mandatory for all personnel running experiments involving animals included in the Animal Welfare Act.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Table S1-S4. Four datasets analyzed in this work. Supplementary tables are formatted as semi-column separated text files. For each dataset, the table contains the individual ID in column1, the parental IDs in column 2 and 3, the family ID in column 4, the type of cross in column 5, length and weight in column 6 and 7. The following columns contain the SNP genotype for each individual.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Besnier, F., Solberg, M.F., Harvey, A.C. et al. Epistatic regulation of growth in Atlantic salmon revealed: a QTL study performed on the domesticated-wild interface. BMC Genet 21, 13 (2020). https://doi.org/10.1186/s12863-020-0816-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12863-020-0816-y

Keywords