- Research article
- Open Access
Atlantic salmon populations reveal adaptive divergence of immune related genes - a duplicated genome under selection
© The Author(s). 2016
- Received: 27 February 2016
- Accepted: 30 June 2016
- Published: 11 August 2016
Populations of Atlantic salmon display highly significant genetic differences with unresolved molecular basis. These differences may result from separate postglacial colonization patterns, diversifying natural selection and adaptation, or a combination. Adaptation could be influenced or even facilitated by the recent whole genome duplication in the salmonid lineage which resulted in a partly tetraploid species with duplicated genes and regions.
In order to elucidate the genes and genomic regions underlying the genetic differences, we conducted a genome wide association study using whole genome resequencing data from eight populations from Northern and Southern Norway. From a total of ~4.5 million sequencing-derived SNPs, more than 10 % showed significant differentiation between populations from these two regions and ten selective sweeps on chromosomes 5, 10, 11, 13–15, 21, 24 and 25 were identified. These comprised 59 genes, of which 15 had one or more differentiated missense mutation. Our analysis showed that most sweeps have paralogous regions in the partially tetraploid genome, each lacking the high number of significant SNPs found in the sweeps. The most significant sweep was found on Chr 25 and carried several missense mutations in the antiviral mx genes, suggesting that these populations have experienced differing viral pressures. Interestingly the second most significant sweep, found on Chr 5, contains two genes involved in the NF-KB pathway (nkap and nkrf), which is also a known pathogen target that controls a large number of processes in animals.
Our results show that natural selection acting on immune related genes has contributed to genetic divergence between salmon populations in Norway. The differences between populations may have been facilitated by the plasticity of the salmon genome. The observed signatures of selection in duplicated genomic regions suggest that the recently duplicated genome has provided raw material for evolutionary adaptation.
- Whole genome duplication
- Immune system
- Selective sweep
- Salmo salar
In addition to being one of the most highly prized freshwater fish for recreational fishing, the Atlantic salmon (Salmo salar L.) is one of the most economically important aquaculture species worldwide. Its natural distribution is throughout the North Atlantic, ranging from Long Island Sound to Ungava Bay in the west and from Northern Portugal to the Barents Sea in the east . This distribution is the result of postglacial colonization of ecosystems that became available when the glacial ice retreated about 10,000 years ago .
Atlantic salmon is characterised by highly significant, hierarchically structured population genetic divergence, with the largest differences observed between the European and North American lineages [3–5]. This divergence is also observed on a regional scale, presumably as a consequence of the colonization process associated with the retreat of the glacier [6, 7]. Moreover, local scale differentiation exists, for example between neighbouring rivers [8–10] and among tributaries within the same river which might be explained by restricted gene flow, genetic drift and adaptation [11–13].
Atlantic salmon exhibit a relatively complex life history that includes spawning and juvenile rearing in freshwater followed by extended ocean migrations to the feeding grounds . As a consequence, salmon go through several distinct transitions that are characterized by changes in behaviour and physiology . They are also able to adapt to varying local conditions throughout their range of environments , exemplified by their ability to inhabit rivers with a wide range of temperatures, from Spain to the colder Arctic latitudes . Previous studies have shown differences in temperature and climate to be associated with genetic differences between salmon populations [7, 18], and latitude also seems to be correlated with allele frequencies of markers relevant to immune response in American and European Atlantic salmon populations, possibly due to temperature induced differences in pathogen-driven selection or other environmental factors [19–21].
In the wild, Atlantic salmon are constantly confronted with a range of pathogens, and have consequently developed numerous innate and adaptive immune mechanisms to overcome infectious challenges . Recent studies suggest that the prevalence of parasites and infectious diseases is increasing in wild populations partly due to global warming [23, 24]. Given the commercial relevance of Atlantic salmon and the recent release of a reference genome , particular effort should be made to identify genes targeted by natural selection in wild Atlantic salmon populations that ultimately can lead to optimized aquaculture practices. The potential relevance of these findings for the Atlantic salmon farming industry is exemplified by the identification of Infectious Pancreatic Necrosis (IPN) Virus resistance  and age at maturity associated genes [27, 28]. A relatively recent whole genome duplication occurred in the salmonid lineage some 80 million years ago , resulting in a partly tetraploid genome undergoing rediploidization. Consequently the genome contains many paralogous regions that could provide raw material for evolution as paralogous genes and regions can diversify and acquire new functions .
Based upon the analysis of microsatellite and SNP markers, several studies have demonstrated that there are highly significant genetic differences between Atlantic salmon populations located in the north and south of Norway [31–33]. However, the genomic regions and genes behind the differences have not been investigated in detail, and consequently, the potential adaptive significance of this genetic divergence remains elusive.
Selective sweeps. FST and heterozygosity estimations in sliding windows were used to identify differentiated loci undergoing selection in either northern or southern populations of Atlantic salmon in Norway. Genes carrying differentiated missense mutations are shown in bold
Sweep length (bp)
Genes in selective sweeps
22,475,000 – 22,800,000
tnmd, elf1, nkrf , sept6, sowahc , rpl39, upf3b, nkap , zbtb33 , atp1b4, lamp2, cul4b, mcts1, c1galt1c1, clic2
76,250,000 – 76,450,000
sept7, anln , prmt7, coq9, polr2c, dgat1, nedd4, adat1, spire2
19,225,000 – 19,300,000
78,325,000 – 78,425,000
trpc2 , rrm1 , slc6a7, rb1 , lpar6
81,025,000 – 81,150,000
eda2r, ar, ophn1
64,700,000 – 65,275,000
pard6g, bloc1s4, nutf2, adnp2 , txnl4a, pqlc1, kcng2, ctdp1, adck5, cpsf1 , parp10 , st3gal1, khdrbs3
46,925,000 – 47,200,000
24,850,000 – 25,075,000
trim13, rnaseh2b , nr0b1, il1rapl1
34,225,000 – 34,525,000
47,075,000 – 47,225,000
stxbp5l, gtf2e1, mx1-1 , mx1-2 , mx2-1 , smcp-1a, mx1-3
Missense mutations in selective sweeps. Several missense mutations were discovered in the selective sweeps. The table lists resequencing derived reference alleles frequencies in the northern (N) and southern (S) populations of salmon in Norway. SNPs from Chr 5 and Chr 25 selected for genotyping are shown in bold
Reference allele frequency (S/N)
Ref/alt amino acid
Description of gene
NF-kB repressing factor.
NF-kB activating protein.
Actin binding protein.
Same as above.
Same as above.
Transient Receptor Cation Channel.
Regulator of entry into cell division.
Non catalytic subunit of RNase H2.
Same as above.
Same as above.
Same as above.
Same as above.
Same as above.
The second most significant selective sweep was found on Chr 5 (Fig. 2b) and included the stress and immune response transcription factor genes nkrf and nkap; zbtb33 encoding a transcriptional regulator binding to methylated CpG dinucleotides, and a gene with unknown function, sowahc. Both Nkrf and Nkap are transcription factors which regulate the NF-kB pathway in which Nkap activates many cell processes including inflammation, immunity, differentiation, cell growth and apoptosis, while Nkrf mediates transcriptional repression of certain Nkap responsive genes. Since NF-kB signaling pathways activate the immune system in the host, these proteins are key targets for proteases expressed by invading pathogens . Functional studies of the Nkap protein have revealed roles for this protein in T-cell maturation  and mRNA splicing . To our knowledge, no previous studies have identified functionally significant SNPs associated with any of the four genes located within this sweep, however one of the SNPs found in nkap is located in a highly conserved region necessary for transcriptional repression. Here the valine is conserved in other species representing the ancestral variant while in Northern Norway methionine is most common (Additional file 1: Figure S2). This finding may be related to differences in immune defense between salmon from these two regions, a suggestion supported by the fact that the NF-kB pathway is differently regulated in IPN resistant salmon . Further studies will reveal how these SNPs modulate the function of NF-kB and virus response or if other functional properties are associated with the selective sweep on Chr 5.
The most significant sweep was found on Chr 25 and contained a cluster of five mx (myxovirus resistance) genes known to be involved in defense against viruses. Three of these mx genes contained missense mutations; mx1-1, mx1-2 and mx2-1 (Fig. 2c). These proteins are dynamin-like GTPases induced upon virus infection through the innate interferon system. It has been shown that they can act broadly against both DNA and RNA viruses and specifically against certain viruses  and studies in mouse, human and chicken have shown that single missense mutations in Mx1 and Mx2 can confer such specific responses [45–48]. It is possible that the identified missense SNPs in the mx genes reflect specific adjustments to different viral disease pressures between northern and southern populations of salmon. We identified missense SNPs in all regions of the protein including a SNP in the antiviral specificity domain in exon 13 (Additional file 1: Figure S3). This SNP represents a structurally relevant amino acid substitution, where arginine seems to be the ancestral variant and cysteine the derived variant dominating in the northern population (Chr 25 position: 47,120,121). Likewise, SNPs in this domain have been associated with specific virus resistance in chicken [49, 50] and pig . SNPs in mx genes have also been investigated in another fish species, the turbot , however, properties related to protection against viruses were not investigated in this study. In rainbow trout (Oncorhynchus mykiss) genetic variation in mx between strains in exon 3–6, was correlated with susceptibility to infectious hematopoietic necrosis virus (IHNV) . This virus also infects Atlantic salmon and our discovery of a missense mutation in exon 6 suggests that salmon could have adapted to the IHNV (Additional file 1: Figure S3). In addition, different strains of rainbow trout display variable susceptibility to this virus . In this study we cannot elucidate the functional significance of the acquired SNPs in mx in Northern Norway, however, further studies will reveal whether any of these changes have been involved in host-virus adaptation .
Atlantic salmon aquaculture involves rearing domesticated fish that originate from commercial breeding programs. Forty wild populations from both the north and south of Norway were sampled when establishing the national breeding programs for salmon . However, analyses of genetic markers demonstrate that there is a dominance of salmon from Southern Norway in the domesticated lines currently in production . Genetic analyses of farmed salmon escapees in Norway have uncovered genetic introgression into native salmon populations in both Northern and Southern Norway, but the biological consequence remains unknown [32, 61, 62]. Consequently the results from the present study, where adaptive genetic divergence between wild salmon from populations located in the north and south of Norway was revealed, it is likely that the potential negative genetic impact of domesticated salmon introgression is greater in populations located in northern regions, since the farmed fish originate mostly from wild Southern Norway populations.
In this study we performed a GWAS by genome resequencing with the aim to screen the Atlantic salmon genome for genetic differentiation between the northern and southern populations in Norway. By investigating eight rivers we uncovered ten particularly striking sweeps including two clusters of immune related genes harboring missense mutations. A feasible interpretation is that different populations of Atlantic salmon have historically been exposed to different selection pressures in the form of pathogens. Some of these adapted alleles could be advantageous for aquaculture production which is currently hampered by a number of diseases, including virus infections . Future studies should include gene editing of immune genes found in these selective sweeps [64, 65] in combination with viral exposure experiments. Within these experiments, viruses relevant to salmon aquaculture should be the primary focus since finding specific resistance alleles can be of significant value to the industry and could also be used for protecting wild fish against high disease pressures posed by open cage aquaculture . Upon finding the protective alleles, selective breeding on individuals with beneficial haplotypes could lead to increased welfare for aquaculture salmon, decreased disease pressure on wild populations and could also be economically favorable for the industry. On the other hand, further studies should investigate the impact of genetic introgression from fertile aquaculture escapees on the adaptive genetic properties in wild populations. To reduce the risk of this unwanted loss of local adaptation and alteration of fitness-related traits, a sustainable solution would be the use of sterile fish in aquaculture, especially in Northern Norway. Future studies should also investigate whether paralogous regions of selective sweeps have undergone positive selection or not, as the latter scenario would suggest an evolutionary mechanism which provides higher adaptive possibilities when a genome is partially tetraploid.
Samples and sampling
Scales from 30 Atlantic salmon males per river were selected from a sample set of 26,000 samples collected in coastal fisheries in Northern Norway. In the Kolarctic Salmon project (http://prosjekt.fylkesmannen.no/Kolarcticsalmon), the multilocus genotypes of all individuals were compared to a genetic baseline consisting of over 180 rivers from Northern Russia and Norway and were assigned to river of origin. Samples that were assigned with high probability to four rivers in Northern Norway; Altaelva, Reisaelva, Lakselv and Tanaelva were generously made available to this study. 30 salmon males from each of four different rivers in Southern Norway, including Årdalselva, Eidselva, Gloppenelva and Suldalslågen were sampled and resequenced in a recent study . In addition to these, we also used male and female salmon DNA from 19 rivers along the Norwegian coast. These included 20 parr individuals from each of the rivers Grense Jakobselv, Neiden, Bergebyelva, Komagelva, Kongsfjordelva, Langfjordelva, Børselva, Stabburselva, Repparfjordselva, Målselv, Laukhelle, Alvsvågvassdraget, Årgårdsvassdraget, Stjørdalselva, Jølstra, Lyseelva, Bjerkreimselva, Storelva and Enningdalselva (represented by numbers in Fig. 1). With the exception of Enningdalselva where the sample was obtained from scales collected by recreational fisheries, these samples were obtained from fins collected by electrofishing of juvenile salmon from mulitple locations in the rivers.
DNA extraction and sequencing
DNA from the 19 rivers for genotyping was extracted from scales or fin samples using Qiagen DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany) according to manufacturer’s recommendations. From salmon belonging to the four populations in Northern Norway total DNA was extracted from scales using Qiagen DNeasy Blood and Tissue Kit. Equal amounts of DNA from ten individuals were pooled to make three pools per river, totaling 30 individuals from each river. Paired-end libraries were constructed using the Genomic DNA Sample Preparation Kit (Illumina, CA, USA) according to manufacturer’s instructions and sequenced on the Illumina HiSeq2000 platform (Illumina, CA, USA) at the Norwegian Sequencing center (https://www.sequencing.uio.no, Oslo, Norway) with each pool sequenced in separate lanes.
Sequence mapping and SNP calling
To ensure high quality sequences, sequenced reads were inspected with FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Adapter sequence removal and quality trimming was done with Cutadapt , resulting in 1,077,839,448 (SD 40,407,492) paired reads on average per river. Sequenced reads were mapped to the most recent release of the salmon genome (AGKD0000000.4) using Bowtie2 (v.2.1.0)  without soft clipping (end-to-end mode). To increase the sensitivity of the mapping, seed length (−L parameter) was set to 18 and the interval between extracted seeds (−i parameter) was set to S,1,1.5 corresponding to the function f (L) = 1 + 1.5*sqrt (L), where L is the read length. Additionally, the maximum number of mismatches per seed (−N parameter) was set to L,0,0.1, corresponding to the function f (L) = 0 + 0.1*L, where L is the read length, and minimum alignment score (−−score-min parameter) was set to L,-0.6,-0.4, corresponding to the function f (L) = −0.6 + −0.4*L, where L is the read length. To remove ambiguously mapped reads the mapping quality threshold was set to 20. To obtain higher sequence coverage, the three sequenced pools per river were merged to a single BAM file using SAMtools merge. SNPs were called using SAMtools mpileup  and the output was parsed using the PoPoolation2 package (mpileup2sync.jar)  with a minimum base quality threshold of 20. For a SNP to be included in the final set of high quality SNPs, minimum coverage of 10 and maximum coverage of 50 (99 % percentile) was required for each river. In addition, the total number of observed minor alleles was required to be at least 8. Recently published whole genome resequencing data from individuals  was downloaded and mapped to the reference genome. The data included three salmon from each of the rivers Tanaelva, Repparfjordelva, Altaelva, Namsenelva, Årgårdsvassdraget, Nausta and Jølstra, where the first three represent populations in Norhtern Norway and the last four represent Southern Norway. Accession numbers for the samples are shown in the caption of Additional file 1: Table S3.
Pairwise fixation index (FST) between all eight sequenced populations was calculated for all high quality SNPs using Hudson’s estimator for FST . FST values were averaged over all SNPs in each population to generate a distance matrix using FST as genetic distance. This matrix was converted to a newick tree using NEIGHBOR from the Phylip package  and a phylogenetic tree was created with NJplot . To find SNPs with significantly different allele frequencies (0.1 % FDR) between populations from Northern and Southern Norway the Cochran-Mantel-Haenszel test for repeated tests of independence from the PoPoolation2 package (cmh-test.pl)  was used. The FDR threshold was determined using the method described in . Allele counts for each river were merged to get the total allele count per SNP in Northern and Southern Norway, corresponding to 120 individuals per geographical region. From this, FST values between the northern and southern populations were estimated using the FST calculation from the PoPoolation2 package (fst-sliding.pl) for each SNP, with --pool-size parameter set to 120. Genomic regions with low values of heterozygosity may indicate SNPs under selection. Therefore heterozygosity values were estimated for north and south of Norway, separately, for each SNPs as 2 * (major allele frequency * minor allele frequency). Sliding windows of 50 kb with steps of 25 kb was used to find genomic regions with high FST values and with low heterozygosity values in either Northern or Southern Norway. This approach is similar to one used to discover genomic regions under selection in other animals . To identify putative selective sweeps it was required that the average FST value of the window was at least 0.17 (above 99.9 % percentile) and that average heterozygosity of the window in either Northern or Southern Norway was at most 0.15 (below 5 % percentile) (Additional file 1: Figure S1). The thresholds were chosen with focus on capturing the outliers in the FST and heterozygosity distributions. Putative sweeps were extended to the sides for as long as the neighboring windows had either average FST of at least 0.17 or heterozygosity of at most 0.15 in either Northern or Southern Norway. If identified sweeps were less than 50 kb apart these were joined to avoid fragmentation of the putative selective sweeps. Genomic windows containing more than 10 % ambiguous bases (Ns) in the reference assembly were discarded to exclude regions with high levels of uncertainty.
Genes in the sweep regions were obtained from the official genome annotation (NCBI Salmo salar Annotation Release 100). Missense mutations in selective sweep regions were identified by manual inspection of the coding sequences. Amino acid sequences of five mx genes found in a selective sweep on Chr 25 were aligned to the homologs Mx1 and Mx2 from human and MxD and MxG from Zebrafish using BLASTP (default parameters). Functional domains in the Mx proteins were assigned using domain information for human Mx1 from UniProt. Amino acid sequences from four genes containing missense mutations in a selective sweep on Chr 5 (nkrf, sowahc, nkap and zbtb33) were aligned to homologous zebrafish and Northern Pike genes using BLASTP with default parameters. Synteny between genes in the sweep on Chr 5 and other animals was found using the UCSC genome browser (https://genome.ucsc.edu) to inspect the syntenic regions of zebrafish, human and mouse. Paralogous regions of the sweeps were identified using TBLASTN (default parameters) with the genes in the sweeps against the salmon genome.
Twenty salmon from 19 rivers along the Norwegian coastline (n = 380) were genotyped using ten of the most significant missense mutations on a Sequenom MassARRAY iPLEX platform (San Diego, CA, USA). Primers and extension primers are listed in Additional file 1: Table S4. The genotyping primers were designed to not target any paralogous genes in the genome.
bp, base pai; Chr, chromosome; FDR, false discovery rate; FST, fixation index; GWAS, genome-wide association study; IHNV, infectious hematopoietic necrosis virus; IPN, infectious pancreatic necrosis; kb, kilo bases; SNP, single nucleotide polymorphism; WGD, whole genome duplication
Samples were generously supplied by a number of agencies and persons. We wish to thank coastal fishermen in Troms and Finnmark, Norwegian Institute of Nature Research, Statens Naturoppsyn, Rådgivende Biologer AS, Fylkesmannen i Nord-Trøndelag and local fishermen in Enningdalselva for providing samples. We would also like to thank Anne Grethe Sørvik for expert technical assistance.
This project was financed by the Norwegian research council (NFR) and their HAVBRUK-BIOTEK 2021 program (project number 226221- SALMAT). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
Genomic sequences from all sequenced pools used in this study have been made available on SRA with Bioproject number PRJNA305872. A list of high quality SNPs has been deposited at http://marineseq.imr.no/northsouth2016/.
AW, KAG, VW, CJR and RBE conceived and designed the experiments. FA, EKS and GD conducted laboratory experiments. FA, EKS, TF, CJR, AW and RBE analyzed the data. VW, EN, MO and JPV provided samples for analysis. EKS, FA, KAG, AW and RBE wrote the first draft of the paper. All authors read and approved to the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Scale samples from adult wild salmon in the rivers were collected by recreational anglers. Scale samples of salmon from four of the northern rivers were collected in commercial coastal fisheries for salmon. Thus no permits/licenses regarding the collection of these samples were required. Juvenile samples from other rivers were collected by own efforts or by several cooperating agencies with permits from the County Governor in the respective counties.
Open AccessThis 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.
- MacCrimmon HR, Gots BL. World distribution of Atlantic salmon, salmon salar. J Fish Res Board Can. 1979;36:422–57.View ArticleGoogle Scholar
- Verspoor E, McCarthy EM, Knox D, Bourke EA, Cross TF. The phylogeography of european Atlantic salmon (salmo salar L.) based on RFLP analysis of the ND1/16sRNA region of the mtDNA. Biol J Linn Soc. 1999;68:129–46.View ArticleGoogle Scholar
- King TL, Kalinowski ST, Schill WB, Spidle AP, Lubinski BA. Population structure of Atlantic salmon (Salmo salar L.): a range-wide perspective from microsatellite DNA variation. Mol Ecol. 2001;10:807–21.View ArticlePubMedGoogle Scholar
- Taggart JB, Verspoor E, Galvin PT, Moran P, Ferguson A. A minisatellite DNA marker for discriminating between European and North American Atlantic salmon (Salmo salar). Can J Fish Aquat Sci. 1995;52:2305–11.View ArticleGoogle Scholar
- Ståhl G. Genetic Population Structure of Atlantic Salmon, Population genetics & Fishery Management. Seattle: University of Washington Press; 1987.Google Scholar
- Tonteri A, Titov S, Veselov A, Zubchenko A, Koskinen MT, Lesbarreres D, Kaluzhin S, Bakhmet I, Lumme J, Primmer CR. Phylogeography of anadromous and non-anadromous Atlantic salmon (Salmo salar) from northern Europe. Ann Zoologici Fennici. 2005;42:1–22.Google Scholar
- Vincent B, Dionne M, Kent MP, Lien S, Bernatchez L. Landscape genomics in Atlantic salmon (salmo salar): searching for gene–environment interactions driving local adaptation. Evolution. 2013;67:3469–87.View ArticlePubMedGoogle Scholar
- Koljonen M-L, Tähtinen J, Säisä M, Koskiniemi J. Maintenance of genetic diversity of Atlantic salmon (Salmo salar) by captive breeding programmes and the geographic distribution of microsatellite variation. Aquaculture. 2002;212:69–92.View ArticleGoogle Scholar
- Ayllon F, Martinez JL, Garcia-Vazquez E. Loss of regional population structure in Atlantic salmon, Salmo salar L., following stocking. ICES J Mar Sci. 2006;63:1269–73.View ArticleGoogle Scholar
- Olafsson K, Pampoulie C, Hjorleifsdottir S, Gudjonsson S, Hreggvidsson GO. Present-Day genetic structure of Atlantic salmon (salmo salar) in Icelandic rivers and Ice-Cap retreat models. PLoS One. 2014;9(2):e86809.View ArticlePubMedPubMed CentralGoogle Scholar
- Primmer CR, Veselov AJ, Zubchenko A, Poututkin A, Bakhmet I, Koskinen MT. Isolation by distance within a river system: genetic population structuring of Atlantic salmon, Salmo salar, in tributaries of the Varzuga River in northwest Russia. Mol Ecol. 2006;15:653–66.View ArticlePubMedGoogle Scholar
- Vähä J-P, Erkinaro J, NiemelÄ E, Primmer CR. Life-history and habitat features influence the within-river genetic structure of Atlantic salmon. Mol Ecol. 2007;16:2638–54.View ArticlePubMedGoogle Scholar
- Dillane E, McGinnity P, Coughlan JP, Cross MC, de Eyto E, Kenchington E, Prodohl P, Cross TF. Demographics and landscape features determine intrariver population structure in Atlantic salmon (Salmo salar L.): the case of the River Moy in Ireland. Mol Ecol. 2008;17:4786–800.View ArticlePubMedGoogle Scholar
- Taylor EB. A review of local adaptation in Salmonidac, with particular reference to Pacific and Atlantic salmon. Aquaculture. 1991;98:185–207.View ArticleGoogle Scholar
- Ryman N, Utter F. Population genetics and fishery management. Seattle: University of Washington Press; 1987.Google Scholar
- de Leániz CG, Fleming IA, Einum S, Verspoor E, Consuegra S, Jordan WC, Aubin-Horth N, Lajus DL, Villanueva B, Ferguson A, et al.: Local Adaptation. In The Atlantic Salmon: Genetics, Conservation and Management. Blackwell Publishing Ltd; 2007: 195–235Google Scholar
- Thorstad EB, Whoriskey F, Rikardsen AH, Aarestrup K: Aquatic Nomads: The Life and Migrations of the Atlantic Salmon. In Atlantic Salmon Ecology. Wiley-Blackwell; 2010: 1–32Google Scholar
- Verspoor E, Jordan WC. Genetic variation at the Me-2 locus in the Atlantic salmon within and between rivers: evidence for its selective maintenance. J Fish Biol. 1989;35:205–13.View ArticleGoogle Scholar
- Tonteri A, Vasemägi A, Lumme J, Primmer CR. Beyond MHC: signals of elevated selection pressure on Atlantic salmon (Salmo salar) immune-relevant loci. Mol Ecol. 2010;19:1273–82.View ArticlePubMedGoogle Scholar
- Dionne M, Miller KM, Dodson JJ, Caron F, Bernatchez L. Clinal variation in MHC diversity with temperature: evidence for the role of host-pathogen interaction on local adaptation in Atlantic salmon. Evolution. 2007;61:2154–64.View ArticlePubMedGoogle Scholar
- Dionne M, Miller KM, Dodson JJ, Bernatchez L. MHC standing genetic variation and pathogen resistance in wild Atlantic salmon. Philos Trans R Soc Lon B. 2009;364:1555–65.View ArticleGoogle Scholar
- Dempsey PW, Vaidya SA, Cheng G. The Art of War: innate and adaptive immune responses. Cell Mol Life Sci. 2003;60:2604–21.View ArticlePubMedGoogle Scholar
- Burge CA, Mark Eakin C, Friedman CS, Froelich B, Hershberger PK, Hofmann EE, Petes LE, Prager KC, Weil E, Willis BL, et al. Climate change influences on marine infectious diseases: implications for management and society. Ann Rev Mar Sci. 2014;6:249–77.View ArticlePubMedGoogle Scholar
- Ficke AD, Myrick CA, Hansen LJ. Potential impacts of global climate change on freshwater fisheries. Rev Fish Biol Fisheries. 2007;17:581–613.View ArticleGoogle Scholar
- Davidson WS, Koop BF, Jones SJ, Iturra P, Vidal R, Maass A, Jonassen I, Lien S, Omholt SW. Sequencing the genome of the Atlantic salmon (Salmo salar). Genome Biol. 2010;11:403.PubMedPubMed CentralGoogle Scholar
- Moen T, Torgersen J, Santi N, Davidson WS, Baranski M, Ødegård J, Kjøglum S, Velle B, Kent M, Lubieniecki KP, et al. Epithelial cadherin determines resistance to infectious pancreatic necrosis virus in Atlantic salmon. Genetics. 2015;200:1313–26.View ArticlePubMedPubMed CentralGoogle Scholar
- Ayllon F, Kjærner-Semb E, Furmanek T, Wennevik V, Solberg M, Dahle G, Taranger GL, Glover KA, Almen MS, Rubin CJ, et al. The vgll3 locus controls Age at maturity in wild and domesticated Atlantic salmon (salmo salar L.) males. PLoS Genet. 2015;11(11):e1005628.View ArticlePubMedPubMed CentralGoogle Scholar
- Barson NJ, Aykanat T, Hindar K, Baranski M, Bolstad GH, Fiske P, Jacq C, Jensen AJ, Johnston SE, Karlsson S, et al.: Sex-dependent dominance at a single locus maintains variation in age at maturity in salmon. Nature. 2015.Google Scholar
- Near TJ, Eytan RI, Dornburg A, Kuhn KL, Moore JA, Davis MP, Wainwright PC, Friedman M, Smith WL. Resolution of ray-finned fish phylogeny and timing of diversification. Proc Natl Acad Sci U S A. 2012;109:13698–703.View ArticlePubMedPubMed CentralGoogle Scholar
- Qian W, Zhang J. Genomic evidence for adaptation by gene duplication. Genome Res. 2014;24:1356–62.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:532–51.View ArticlePubMedGoogle Scholar
- Glover KA, Quintela M, Wennevik V, Besnier F, Sorvik AG, Skaala O. Three decades of farmed escapees in the wild: a spatio-temporal analysis of Atlantic salmon population genetic structure throughout Norway. PLoS One. 2012;7:e43129.View ArticlePubMedPubMed CentralGoogle Scholar
- Ozerov M, Vasemagi A, Wennevik V, Diaz-Fernandez R, Kent M, Gilbey J, Prusov S, Niemela E, Vaha JP. Finding markers that make a difference: DNA pooling and SNP-arrays identify population informative markers for genetic stock identification. PLoS One. 2013;8:e82434.View ArticlePubMedPubMed CentralGoogle Scholar
- Bhatia G, Patterson N, Sankararaman S, Price AL. Estimating and interpreting FST: the impact of rare variants. Genome Res. 2013;23:1514–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith JM, Haigh J. The hitch-hiking effect of a favourable gene. Genet Res. 1974;23:23–35.View ArticlePubMedGoogle Scholar
- Carneiro M, Rubin CJ, Di Palma F, Albert FW, Alfoldi J, Barrio AM, Pielberg G, Rafati N, Sayyab S, Turner-Maier J, et al. Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication. Science. 2014;345:1074–9.View ArticlePubMedGoogle Scholar
- Omura M, Mombaerts P. Trpc2-expressing sensory neurons in the main olfactory epithelium of the mouse. Cell Rep. 2014;8:583–95.View ArticlePubMedGoogle Scholar
- Kaufmann M, Feijs KL, Luscher B. Function and regulation of the mono-ADP-ribosyltransferase ARTD10. Curr Top Microbiol Immunol. 2015;384:167–88.PubMedGoogle Scholar
- Kind B, Muster B, Staroske W, Herce HD, Sachse R, Rapp A, Schmidt F, Koss S, Cardoso MC, Lee-Kirsch MA. Altered spatio-temporal dynamics of RNase H2 complex assembly at replication and repair sites in Aicardi-Goutieres syndrome. Hum Mol Genet. 2014;23:5950–60.View ArticlePubMedGoogle Scholar
- Hodgson A, Wan F: Interference with nuclear factor kappaB signaling pathway by pathogen-encoded proteases: global and selective inhibition. Mol Microbiol. 2015.Google Scholar
- Thapa P, Das J, McWilliams D, Shapiro M, Sundsbak R, Nelson-Holte M, Tangen S, Anderson J, Desiderio S, Hiebert S, et al. The transcriptional repressor NKAP is required for the development of iNKT cells. Nat Commun. 2013;4:1582.View ArticlePubMedPubMed CentralGoogle Scholar
- Burgute BD, Peche VS, Steckelberg AL, Glockner G, Gassen B, Gehring NH, Noegel AA. NKAP is a novel RS-related protein that interacts with RNA and RNA binding proteins. Nucleic Acids Res. 2014;42:3177–93.View ArticlePubMedGoogle Scholar
- Cofre C, Gonzalez R, Moya J, Vidal R. Phenotype gene expression differences between resistant and susceptible salmon families to IPNV. Fish Physiol Biochem. 2014;40:887–96.View ArticlePubMedGoogle Scholar
- Mitchell PS, Emerman M, Malik HS. An evolutionary perspective on the broad antiviral specificity of MxA. Curr Opin Microbiol. 2013;16:493–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee SH, Vidal SM. Functional diversity of Mx proteins: variations on a theme of host resistance to infection. Genome Res. 2002;12:527–30.View ArticlePubMedGoogle Scholar
- Goujon C, Moncorge O, Bauby H, Doyle T, Ward CC, Schaller T, Hue S, Barclay WS, Schulz R, Malim MH. Human MX2 is an interferon-induced post-entry inhibitor of HIV-1 infection. Nature. 2013;502:559–62.View ArticlePubMedGoogle Scholar
- Sironi M, Biasin M, Cagliani R, Gnudi F, Saulle I, Ibba S, Filippi G, Yahyaei S, Tresoldi C, Riva S, et al. Evolutionary analysis identifies an MX2 haplotype associated with natural resistance to HIV-1 infection. Mol Biol Evol. 2014;31:2402–14.View ArticlePubMedGoogle Scholar
- Mitchell PS, Patzina C, Emerman M, Haller O, Malik HS, Kochs G. Evolution-guided identification of antiviral specificity determinants in the broadly acting interferon-induced innate immunity factor MxA. Cell Host Microbe. 2012;12:598–604.View ArticlePubMedPubMed CentralGoogle Scholar
- Sasaki K, Yoneda A, Ninomiya A, Kawahara M, Watanabe T. Both antiviral activity and intracellular localization of chicken Mx protein depend on a polymorphism at amino acid position 631. Biochem Biophys Res Commun. 2013;430:161–6.View ArticlePubMedGoogle Scholar
- Ko JH, Takada A, Mitsuhashi T, Agui T, Watanabe T. Native antiviral specificity of chicken Mx protein depends on amino acid variation at position 631. Anim Genet. 2004;35:119–22.View ArticlePubMedGoogle Scholar
- Nakajima E, Morozumi T, Tsukamoto K, Watanabe T, Plastow G, Mitsuhashi T. A naturally occurring variant of porcine Mx1 associated with increased susceptibility to influenza virus in vitro. Biochem Genet. 2007;45:11–24.View ArticlePubMedGoogle Scholar
- Abollo E, Ordas C, Dios S, Figueras A, Novoa B. Molecular characterisation of a turbot Mx cDNA. Fish Shellfish Immunol. 2005;19:185–90.View ArticlePubMedGoogle Scholar
- Trobridge GD, LaPatra SE, Kim CH, Leong JC. Mx mRNA expression and RFLP analysis of rainbow trout Oncorhynchus mykiss genetic crosses selected for susceptibility or resistance to IHNV. Dis Aquat Organ. 2000;40:1–7.View ArticlePubMedGoogle Scholar
- Purcell MK, Lapatra SE, Woodson JC, Kurath G, Winton JR. Early viral replication and induced or constitutive immunity in rainbow trout families with differential resistance to Infectious hematopoietic necrosis virus (IHNV). Fish Shellfish Immunol. 2010;28:98–105.View ArticlePubMedGoogle Scholar
- Daugherty MD, Malik HS. Rules of engagement: molecular insights from host-virus arms races. Annu Rev Genet. 2012;46:677–700.View ArticlePubMedGoogle Scholar
- Gidskehaug L, Kent M, Hayes BJ, Lien S. Genotype calling and mapping of multisite variants using an Atlantic salmon iSelect SNP array. Bioinformatics. 2011;27:303–10.View ArticlePubMedGoogle Scholar
- Ohno S, Wolf U, Atkin NB. Evolution from fish to mammals by gene duplication. Hereditas. 1968;59:169–87.View ArticlePubMedGoogle Scholar
- Warren IA, Ciborowski KL, Casadei E, Hazlerigg DG, Martin S, Jordan WC, Sumner S. Extensive local gene duplication and functional divergence among paralogs in Atlantic salmon. Genome Biol Evol. 2014;6:1790–805.View ArticlePubMedPubMed CentralGoogle Scholar
- Glasauer SM, Neuhauss SC. Whole-genome duplication in teleost fishes and its evolutionary consequences. Mol Genet Genomics. 2014;289:1045–60.View ArticlePubMedGoogle Scholar
- Gjedrem T, Gjoen HM, Gjerde B. Genetic-origin of Norwegian farmed Atlantic salmon. Aquaculture. 1991;98:41–50.View ArticleGoogle Scholar
- 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:4.Google Scholar
- Skaala O, Wennevik V, Glover KA. Evidence of temporal genetic change in wild Atlantic salmon, Salmo salar L., populations affected by farm escapees. Ices J Mar Sci. 2006;63:1224–33.View ArticleGoogle Scholar
- Collet B. Innate immune responses of salmonid fish to viral infections. Dev Comp Immunol. 2014;43:160–73.View ArticlePubMedGoogle Scholar
- Edvardsen RB, Leininger S, Kleppe L, Skaftnesmo KO, Wargelius A. Targeted mutagenesis in Atlantic salmon (Salmo salar L.) using the CRISPR/Cas9 system induces complete knockout individuals in the F0 generation. PLoS One. 2014;9:e108622.View ArticlePubMedPubMed CentralGoogle Scholar
- Wargelius A, Leininger S, Skaftnesmo KO, Kleppe L, Andersson E, Taranger GL, Schulz RW, Edvardsen RB. Dnd knockout ablates germ cells and demonstrates germ cell independent sex differentiation in Atlantic salmon. Sci Rep. 2016;6:21284.View ArticlePubMedPubMed CentralGoogle Scholar
- Taranger GL, Karlsen O, Bannister RJ, Glover KA, Husa V, Karlsbakk E, Kvamme BO, Boxaspen KK, Bjorn PA, Finstad B, et al. Risk assessment of the environmental impact of Norwegian Atlantic salmon farming. Ices J Mar Sci. 2015;72:997–1021.View ArticleGoogle Scholar
- Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–2.View ArticleGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. Genome project data processing S: the sequence alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Kofler R, Pandey RV, Schlotterer C. PoPoolation2: identifying differentiation between populations using sequencing of pooled DNA samples (Pool-Seq). Bioinformatics. 2011;27:3435–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Fink WL. Microcomputers and phylogenetic analysis. Science. 1986;234:1135–9.View ArticlePubMedGoogle Scholar
- Perriere G, Gouy M. WWW-query: an on-line retrieval system for biological sequence banks. Biochimie. 1996;78:364–9.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc B Met. 1995;57:289–300.Google Scholar