# Comparing genomic signatures of domestication in two Atlantic salmon (Salmo salar L.) populations with different geographical origins

Maria E López, Laura Benestan, Jean-Sébastien Moore, Charles Perrier, John Gilbey, Alex Di Genova, Alejandro Maass, Diego Diaz, Jean-Paul Lhorente, Katharina Correa, et al.

► To cite this version:

Maria E López, Laura Benestan, Jean-Sébastien Moore, Charles Perrier, John Gilbey, et al.. Comparing genomic signatures of domestication in two Atlantic salmon (Salmo salar L.) populations with different geographical origins. Evolutionary Applications, 2018, 12 (1), pp.137-156. 10.1111/eva.12689. hal-02921349

HAL Id: hal-02921349

https://hal.inrae.fr/hal-02921349

Submitted on 25 Aug 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Comparing genomic signatures of domestication in two Atlantic salmon (*Salmo salar* L.) populations with different geographical origins

Maria E. López<fnref n="1" />,<fnref n="2" />* | Laura Benestan<fnref n="3" />* | Jean-Sebastien Moore<fnref n="3" /> |

Charles Perrier<fnref n="4" /> | John Gilbey<fnref n="5" /> | Alex Di Genova<fnref n="6" /> | Alejandro Maass<fnref n="6" /> |

Diego Diaz<fnref n="6" /> | Jean-Paul Lhorente<fnref n="7" /> | Katharina Correa<fnref n="7" /> | Roberto Neira<fnref n="2" /> |

Louis Bernatchez<fnref n="3" /> | José M. Yáñez<fnref n="1" />,<fnref n="7" />,<fnref n="8" /> ID

## Correspondence

José M. Yáñez, Facultad de Ciencias
Agronómicas, Universidad de Chile,
Santiago, Chile.

Email: jmayanez@uchile.cl

## Funding information

Corporación de Fomento de la Producción,
Grant/Award Number: 11IEI-12843 and
12PIE-17669; CONICYT, Grant/Award
Number: 21120382; Canada-Chile
Leadership Exchange Scholarship; Ministerio
de Economía, Fomento y Turismo

## Abstract

Selective breeding and genetic improvement have left detectable signatures on the genomes of domestic species. The elucidation of such signatures is fundamental for detecting genomic regions of biological relevance to domestication and improving management practices. In aquaculture, domestication was carried out independently in different locations worldwide, which provides opportunities to study the parallel effects of domestication on the genome of individuals that have been selected for similar traits. In this study, we aimed to detect potential genomic signatures of domestication in two independent pairs of wild/domesticated Atlantic salmon populations of Canadian and Scottish origins, respectively. Putative genomic regions under divergent selection were investigated using a 200K SNP array by combining three different statistical methods based either on allele frequencies (LFMM, Bayesian) or haplotype differentiation (Rsb). We identified 337 and 270 SNPs potentially under divergent selection in wild and hatchery populations of Canadian and Scottish origins, respectively. We observed little overlap between results obtained from different statistical methods, highlighting the need to test complementary approaches for

<footnote n="1">¹Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Santiago, Chile</footnote>

<footnote n="2">²Facultad de Ciencias Agronómicas, Universidad de Chile, Santiago, Chile</footnote>

<footnote n="3">³IBIS, Institut de Biologie Intégrative et des Systèmes, Université Laval, Québec City, Québec, Canada</footnote>

<footnote n="4">⁴Centre d'Écologie Fonctionnelle et Évolutive, Unité Mixte de Recherche CNRS 5175, Montpellier, France</footnote>

<footnote n="5">⁵Marine Scotland Science, Freshwater Fisheries Laboratory, Faskally, Pitlochry, UK</footnote>

<footnote n="6">⁶Laboratory of Bioinformatics and Mathematics of the Genome, Center for Mathematical Modeling (UMI 2807 CNRS) and Center for Genome Regulation (Fondap 15090007), Universidad de Chile, Santiago, Chile</footnote>

<footnote n="7">⁷Aquainnovo, Puerto Montt, Chile</footnote>

<footnote n="8">⁸Núcleo Milenio INVASAL, Concepción, Chile</footnote>

<footnote>*Both authors contributed equally to the manuscript.</footnote>

<footnote>This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium,
provided the original work is properly cited.</footnote>

<footnote>© 2018 The Authors. *Evolutionary Applications* published by John Wiley & Sons Ltd</footnote>
detecting a broad range of genomic footprints of selection. The vast majority of the outliers detected were population-specific but we found four candidate genes that were shared between the populations. We propose that these candidate genes may play a role in the parallel process of domestication. Overall, our results suggest that genetic drift may have override the effect of artificial selection and/or point toward a different genetic basis underlying the expression of similar traits in different domesticated strains. Finally, it is likely that domestication may predominantly target polygenic traits (e.g., growth) such that its genomic impact might be more difficult to detect with methods assuming selective sweeps.

### KEYWORDS

*Salmo salar*, selective sweeps, single nucleotide polymorphisms

## 1 | INTRODUCTION

Domestication, "the process by which captive animals adapted to man and the environment he provides," has led to several genetic changes over generations in various animal and plant species (Price, 1984). This evolutionary process may induce similar phenotypic changes in populations of independent origins within the same species when the same phenotypic traits are subjected to similar selective pressures. Phenotypes of these animals have progressively evolved due to the combined influence of domestication selection through reproduction in captivity and human directional selective breeding (Andersson, 2012). This may result in a remarkable phenotypic diversity within domestic species as well as a wide variety of genetic adaptations to both environmental conditions and production systems (Andersson, 2012, 2013). Domestication has thus shaped the genetic diversity of these species throughout history, and their present genomes may contain traceable signatures of selection (Utsunomiya, Pérez, O'Brien, Sonstegard, & Garcia, 2015).

In most fishes, domestication is a recent process compared with any other livestock species (Gjedrem, 2005; López, Neira, & Yáñez, 2015). Therefore, important knowledge gaps still remain about the consequences of domestication on behavior, physiology, and morphology of fishes (Teletchea & Fontaine, 2014) compared to domesticated mammals and birds, for which more data are available (Kelley, Brown, Therkildsen, & Foote, 2016; Lorenzen, Beveridge, & Mangel, 2012). Nevertheless, there is no reason to consider domestication of terrestrial and aquatic animals distinctively (Balon, 2004; Teletchea & Fontaine, 2014). Hence, we might expect that, similar to birds and mammals, selection for the phenotypes contributing to domestication goals has also impacted the extent and distribution of variability within the genomes of fishes. Indeed, domestication has greatly impacted the phenotypes of domesticated aquatic species (e.g., Atlantic salmon [*Salmo salar* L.], rainbow trout [*Oncorhynchus mykiss* W.], tilapia [*Oreochromis niloticus* L.], common carp [*Cyprinus carpio* L.]) in which selective breeding programs have been implemented, focusing on specific economic objectives (Gjedrem, 2012; Rye, Gjerde, & Gjedrem, 2010).

Detecting genomic signatures of selection is a major goal of modern population genetics (Fariello et al., 2014; Nielsen, Hellmann, Hubisz, Bustamante, & Clark, 2007) as it enhances our knowledge of the molecular mechanisms shaping the genome as well as providing functional information on specific genes/genomic regions that would be of interest for breeding programs (Qanbari & Simianer, 2014). Recently, the progress of high-throughput and cost-effective genotyping techniques has offered a unique opportunity to analyze large datasets of domesticated species to study genome changes in response to domestication events (Druet, Pérez-Pardal, Charlier, & Gautier, 2013; Ma et al., 2015). In livestock species, genomewide analyses have already shown promising results in mapping traits of economic interest, such as genomic regions related to milk production in cattle (*Bos taurus*) (Qanbari et al., 2011), muscle development in pig (*Sus scrofa*) (Amaral et al., 2011), coat pigmentation and skeletal morphology in sheep (*Ovis aries*) (Kijas et al., 2012), or gait and size in horses (*Equus caballus*) (Petersen et al., 2013). The practical aspect of these population genomics studies lies on the possibility of detecting selected genes associated with traits of economic interest and acts as complement to gene mapping approaches (e.g., genomewide association studies [GWAS]) that may help to further genetically improve these traits on domestic species (Qanbari & Simianer, 2014). Furthermore, this knowledge is important from an evolutionary perspective by highlighting traits that have been exposed to natural and artificial selection as well as using this information to design and/or update breeding programs for conservation purposes worldwide (Brito et al., 2017; Cesconeto et al., 2017; Zhao, McParland, Kearney, Lixin, & Berry, 2015).

Recently, methods such as whole genome sequencing or genomewide SNP arrays have enabled the screening of a large part of the genome to detect signatures of selection in domestic and natural populations (Druet et al., 2013; Fuentes-Pardo & Ruzzante, 2017; Ma et al., 2015). Moreover, new genome scan approaches have been developed with the goal of efficiently and accurately identifying genomic footprints of selection out of the thousands of markers screened (Jensen, Foll, & Bernatchez, 2016). Several analytical methods are now available, ranging from population
differentiation analyses, based on $F_{ST}$ calculation (Porto-Neto, Lee, Lee, & Gondro, 2013), to environmental association methods (Cesconeto et al., 2017; Frichot, Schoville, Bouchard, & François, 2013), which aim to detect genetic variants associated with specific environmental factors. Population differentiation methods are expected to detect only strong signatures of selection (e.g., hard sweep), whereas alternative methods such as environmental association would be able to detect subtle signal of selection (e.g., soft sweep), which makes the combination of the two a particularly promising tool for uncovering any genomic footprints of selection (Benestan et al., 2016; Rellstab, Gugerli, Eckert, Hancock, & Holderegger, 2015).

Atlantic salmon aquaculture has a high economic value, and breeding programs have been implemented worldwide to improve genetic performance of domesticated individuals (Gjedrem, Robinson, & Rye, 2012; Rye et al., 2010). In this perspective, several studies have already detected evidences for genomic signatures of domestication in Atlantic salmon (Bourret, O'Reilly, Carr, Berg, & Bernatchez, 2011; Mäkinen, Vasemägi, McGinnity, Cross, & Primmer, 2014; Martinez et al., 2013; Roberge, Einum, Guderley, & Bernatchez, 2006; Vasemägi et al., 2012) and showed that these signatures were related to important biological functions including immunity, lipid transport, fatty acids, steroid metabolism, and cellular cycle/growth (Devlin, Sakhrani, Tymchuk, Rise, & Goh, 2009; Sauvage et al., 2010; Wringe et al., 2010). Yet, little is known about how the process of domestication altered the ancestral wild Atlantic salmon genome, partly due to the fact that the demographic history of wild Atlantic salmon populations is complex. Indeed, wild populations of Atlantic salmon exhibit complex patterns of hierarchical genetic differentiation at continental, regional, and local scales (Dillane et al., 2008; Dionne, Caron, Dodson, & Bernatchez, 2008; Moore et al., 2014; Perrier, Guyomard, Bagliniere, & Evanno, 2011; Spidle et al., 2003; Verspoor et al., 2005). These complex patterns are the product of postglacial colonization (King, 2007; Rougemont & Bernatchez, 2018) combined with contemporary ecological processes such as homing behavior, which maintains genetic differentiation among rivers, as well as spatially varying environmental conditions (Bourret et al., 2013; Moore et al., 2014; Rougemont & Bernatchez, 2018).

Across the entire Atlantic salmon natural range, European and North American populations exhibit a deep genetic divergence and belong to two divergent glacial lineages which likely separated more than 1,000,000 years ago (Rougemont & Bernatchez, 2018). While these two lineages did not share the same demographic history, domestic populations derived from them have been selected for similar economically important phenotypic traits such as growth, age at sexual maturity, disease resistance, and flesh quality (Gjedrem, 2012), which suggests that domestication may involve parallel evolution. Finding parallel genetic differences between these two populations would reinforce the hypothesis that candidate genomic regions may be linked to molecular processes underlying domestication. Indeed, signals of selection detected among farmed fish belonging to different origins are likely to result from selection rather

than having been inherited by chance or from migration of selected allele across breeding strains, which limits false-positive detection. In Atlantic salmon, previous transcriptomic studies showed that European and North American populations exhibit parallel differential gene expression profiles between hatchery strains and wild populations in the same set of genes (Roberge et al., 2006; Sauvage et al., 2010). Yet, Vasemägi et al. (2012) using 261 SNPs and 70 microsatellites loci found little support for parallel evolution as they found none of the candidate genomic regions potentially associated with domesticated and improved traits overlapped across three different lineages (Ireland, Sweden, and Canada). Using a genomewide analysis with a 6.5K SNP array, Mäkinen et al. (2014) revealed that candidate loci associated with domestication were located on different chromosomes for the North American and European lineages. In this study, we aimed at further investigating the genomic footprints of domestication as well as parallel genomic changes between two pairs of wild/domestic populations of Atlantic salmon with European and North American origins using a 200K SNP chip that may offer higher resolution than achieved in previous studies. We applied three complementary statistical approaches to identify genomic regions putatively under selection: (a) a genome scan based on population differentiation indices ($F_{ST}$) (BAYESCAN), (b) genome-wide environmental association (LFMM), and (c) haplotype extended patterns (Rsb). From the set of candidate genes detected by these approaches, we identified their biological function to delineate their potential importance in the process of domestication of Atlantic salmon populations.

## 2 | MATERIALS AND METHODS

### 2.1 | Studied populations

This study was conducted using two wild/domestic population pairs of Atlantic salmon from two different geographical origins: Canada and Scotland, for a total of 183 individuals (see Table 1 for details). The purpose of this design was to perform two independent comparisons between domestic strains and their wild counterpart populations. Samples from Canadian domestic population (Can-D) were collected from a domestic strain of Atlantic salmon cultured in Chile. This strain was established from eggs obtained from the Gaspé Bay (where the Saint-Jean River flows), Québec, in the 1950s (Withler, Supernault, & Miller, 2005). Using 11 microsatellites Withler et al. (2005) showed no evidence of introgression from other strains. We presume that fish of this strain were brought to the United States in the late 80s and early 90s and kept in a hatchery located in the state of Washington for about two generations before being introduced to Chile between 1996 and 1998 (J. P. Lhorente, personal communication). This strain has been subjected to intensive selection for rapid growth rate as the main objective of breeding, and also for low incidence of early sexual maturation (J. P. Lhorente, personal communication). This Chilean farmed population has been reared in the environmental conditions of the Los Lagos Region (42°S 72°O,
<table><caption>TABLE 1 Atlantic salmon populations analyzed in this study with geographical origin and sample size (n)</caption><thead><tr><th>Population code</th><th>Geographical origin</th><th>n<sup>a</sup></th><th>n<sub>QC</sub><sup>b</sup></th><th>%<sup>c</sup></th><th>Observations</th></tr></thead><tbody><tr><td>Can-W</td><td>Canada</td><td>46</td><td>44</td><td>95.7</td><td>Wild population from Saint-Jean River, Québec, Canada</td></tr><tr><td>Can-D</td><td>Canada</td><td>40</td><td>37</td><td>92.5</td><td>Farmed population selected for fast growth rate and characterized for late sexual maturity</td></tr><tr><td>Sct-W</td><td>Scotland</td><td>47</td><td>41</td><td>87.2</td><td>Wild population from the East Coast of Scotland.</td></tr><tr><td>Sct-D</td><td>Scotland</td><td>50</td><td>43</td><td>86</td><td>Farmed population selected for fast growth rate and characterized for early sexual maturity</td></tr></tbody></table>

Chile) for about four generations. The wild population from Canada (Can-W) included adult individuals sampled through angling in Saint-Jean River on the same year (Dionne et al., 2008). The Saint-Jean River flows into Gaspe Bay, along with York and Dartmouth rivers, and salmon population from Saint-Jean River is the largest one of the three. Furthermore, it has been shown that genetic differentiation among individuals from these three rivers is weak ($\theta_{ST}$ (Weir & Cockerham, 1984) = 0.011) (Dionne et al., 2008), enough to consider our samples as a good representation of the possible wild ancestors of domestic Canadian population used in this study. On the other hand, the Scottish farmed population (Sct-D) is a strain cultivated in Chile and it was established with fish from Loch Lochy, located on the West Coast of Scotland. This strain was maintained in a commercial genetic improvement program based in Scotland (Johnston, Alderson, Sandham, Mitchell, et al., 2000) and is characterized by a high proportion of early maturing fish and faster growth rate than strains with high proportion of late maturing fish (Johnston, Alderson, Sandham, Dingwall, et al., 2000; Johnston, Alderson, Sandham, Mitchell, et al., 2000). In 1986, this Scottish farmed population was transferred to Los Lagos Region (42°S 72°O, Chile), presenting a high incidence of fish with early sexual maturity and rapid growth as well, based on field records from 2004. Since then, this population has been adapted to environmental, geographical, and climatic conditions in Chile and selected for rapid growth rate using mass selection for at least six generations (C. Soto, personal communication). The wild population from Scotland (Sct-W) comprised juvenile individuals from South Esk River on the East Coast of Scotland collected in 2011, which is a region with sparse aquaculture activity; therefore, there is a low probability that escaped farmed fish had affected the genetic diversity of wild individuals. We used samples from this location, instead of fish from the West Coast of Scotland, as the latter has an important aquaculture activity and has also been stocked from different sources. Therefore, wild fish from this region does not represent a pure Scottish wild Atlantic salmon. Individuals were collected using electrofishing followed by anesthesia and partial fin-clipping. All fish were released back to the capture location (Yáñez et al., 2016). Fish collecting work has been reviewed both by Marine Scotland Review Committee and the United Kingdom Home Office (Project License 60/4251).

## 2.2 | Genotyping and quality control

Genomic DNA was extracted from fin clips using a DNeasy Blood & Tissue Kit (QIAGEN). For genotyping, we used an Affymetrix Axiom® myDesign Custom Array of 200K SNPs, which was developed using thirteen European fish (Scottish and Norwegian origin) and seven North American fish, as described by Yáñez et al. (2016). SNP quality control (QC) was carried out using Axiom Genotyping Console (GTC, Affymetrix) and SNPolisher (an R package developed by Affymetrix), based on SNP clustering metrics and call rate status. Thus, the QC was performed by (a) removing SNP that did not match with high quality clustering patterns, which are defined according to the best practices recommended by Affymetrix through using SNPs with (i) PolyHighResolution (good cluster resolutions for homozygote and heterozygote samples and at least two occurrences of minor allele) and (ii) NoMinorHom (two distinctive cluster with no-minor homozygous genotypes), (b) removing SNP loci with call rate lower than 0.95, and (c) discarding individuals with call rate lower than 0.90. Further filtering steps included removing loci deviating from Hardy-Weinberg equilibrium calculated for each population separately, and subsequently, we removed SNPs shared among all populations, and this stage was part of the validation of the SNP chip as described by Yáñez et al. (2016). Loci with minor allele frequency (MAF) lower than 5% within each pair of populations joined (i.e., Can-W and Can-D comprised one group; Sct-W and Sct-D comprised another group) using PLINK v1.07 (Purcell et al., 2007).

## 2.3 | Basic statistics and structure analysis

To investigate population differentiation, we calculated the pairwise Weir and Cockerham's $F_{ST}$ (1984) estimator across all loci among populations, using VCFTools (Danecek et al., 2011). Observed and expected heterozygosity ($H_O$ and $H_E$) were estimated using PLINK v1.07 (Purcell et al., 2007), and confidence intervals for these statistics were calculated using the R package boot with 1,000 bootstrap replicates (Canty & Ripley, 2016). Genomic distribution of $F_{ST}$ was further analyzed by a kernel-smoothing approach using the R package Lokern. A local bandwidth of ~400 was used for fitting the kernel-smoothed regression line. Two clustering methods were used

<footnote><sup>a</sup>Initial number of individuals. <sup>b</sup>Number of individuals passing sample call rate of 90%. <sup>c</sup>Percentage of individuals passing sample call rate of 90%.</footnote>
to assess the extent of genetic structure among the populations. First, we performed a principal component analysis (PCA) implemented in the R package *adegenet* (Jombart, 2008). Second, we used the maximum likelihood estimation of individual ancestries through ADMIXTURE (Alexander, Novembre, & Lange, 2009) software. ADMIXTURE analysis was run using 2,000 bootstraps, and the number of ancestral populations was set from 1 to 10 (K). The optimal K was selected based on the lowest cross-validation error and a visual inspection of the co-ancestry values for each individual.

## 2.4 | Identification of selection signatures

Three methods were implemented to detect signatures of selection: (a) population differentiation, (b) environmental association, and (c) haplotype extended patterns. First, we used the population differentiation approach with a Bayesian likelihood method, implemented in BAYESCAN v.2.1, which identifies candidate loci under selection using differences in allele frequencies between populations. The algorithm uses a reversible-jump Markov Chain Monte Carlo to explore models with or without selection. To estimate the probability that a locus is under selection, the program uses a Bayes factor for two models: one assuming selection and another assuming neutrality given the data (Foll & Gaggiotti, 2008). Bayesian was run using a pairwise differentiation test (i.e., wild vs. domestic strains). For both independent runs on Canadian and Scottish pair of populations, we used BAYESCAN v.2.1 default parameters. A Bayes factor between 10 and 32 ($\log_{10} = 1-1.5$) for one locus is considered putatively under divergent selection, while between 32 and 100 ($\log_{10} = 1.5-2$) and Bayes factors above 100 ($\log_{10} > 2$) are considered putatively under very strong and decisive divergent selection, respectively. We used a Bayes factor of 10 ($\log_{10} = 1.30$) as threshold to select candidate loci under selection.

Second, we performed an environmental association approach using latent factor mixed models implemented in the software LFMM (Frichot et al., 2013). This method calculates the correlations between allele frequencies and environmental variables. We defined *Domesticated* and *Wild* conditions as an environmental dichotomous variable (0 = wild; 1 = domesticated). This method can efficiently estimate random effects due to population history and isolation-by-distance and can reduce false-positive associations compared to genome scans (Frichot et al., 2013). We applied the *latent factor mixed models* algorithm and calculated the $|z|$ scores for all of the SNPs in both independent comparisons for Canadian and Scottish populations, using 20,000 iterations, 10,000 burn-in iterations in the Gibbs Sampling algorithm and a seed of 1,000. The number of latent factors, which is the best number of clusters describing population structure of the original data, was chosen according to the results of the clustering analyses (i.e., PCA and ADMIXTURE) and as recommended by Frichot et al. (2013).

Finally, we investigated haplotype extended patterns using the *Standardized log-ratio of integrated EHHS (iES) between pairs of populations* (Rsb) analysis (Tang, Thornton, & Stoneking, 2007). This test is based on the extended haplotype homozygosity (EHH) statistic and contrasts EHH patterns of the same haplotype between populations. The SNPs

position on the genome is needed estimating EHH, which was possible thanks to the availability of the chromosomal-level genome assembly for Atlantic salmon (Lien et al., 2016). Rsb is defined as the natural log of the ratio between $iES_{POP1}$ and $iES_{POP2}$, where iES is the integrated EHHS (site-specific EHH) for each SNP within each population. As the Rsb values are normally distributed, a Z-test was applied to identify significant SNPs under selection between wild and domestic strains of both Canadian and Scottish origins. A positive value of Rsb indicate $iES_{POP1}$ is larger than $iES_{POP2}$; therefore, Pop1 has longer haplotypes than Pop2, and hence, positive Rsb scores suggest that selection occurred in the alternative population (domestic population), whereas the negative scores suggest that selection occurred in the reference population (wild population). One-sided *p*-values were obtained as $-\log_{10}(1-\Phi(Rsb))$, where $\Phi(Rsb)$ represents the Gaussian distribution function for the Rsb statistic. We used $-\log_{10}(p\text{-value}) = 3$ as a threshold to define significant Rsb values *Standardized log-ratio of integrated EHHS (iES) between pairs of populations* using rehh R package (Gautier & Vitalis, 2012).

Bayesian and LFMM methods can directly handle the SNP genotypes, whereas Rsb requires haplotypes. We imputed the missing genotypes and inferred haplotypes from genotypes using BEAGLE v.3 (Browning & Browning, 2009) with default parameters. Phasing was carried out separately for each population; the SNP position on chromosome was obtained from the Atlantic salmon genome reference ICSAG_v2 (Lien et al., 2016).

## 2.5 | SNP annotation and functional enrichment analysis

Genomic regions harboring each locus putatively under selection, detected by at least one method (BAYESCAN, LFMM or Rsb), were interrogated for genes annotated to the Atlantic salmon genome reference ICSAG_v2 (GenBank: GCA_000233375.4) using SnpEff v4.3 (Cingolani et al., 2012). This functional annotation is a complete representation of the genome containing 37,206 high-confidence protein-coding genes that have been assigned to a putative functional annotation based on homology within the SwissProt database (Lien et al., 2016).

Using the salmon transcripts, we performed a blastx on the zebrafish (*Danio rerio*) peptide reference database (downloaded from http://www.ensembl.org). For each unique peptide in this species, the corresponding gene id and associated gene name were retrieved from ensembl biomart database (retrieved from http://www.ensembl.org/biomart). Enrichment analysis was conducted using the online tool Panther Classification System (http://www.pantherdb.org) with default parameters.

# 3 | RESULTS

## 3.1 | Genotyping and SNP filtering

After genotyping QC, following the conditions recommended by Affymetrix described in Materials and Methods, 159,099 SNPs were categorized in two good quality metrics (a) poly high resolution (good
cluster resolution of homozygote and heterozygote samples and at least two occurrences of minor allele) and (b) no-minor allele homozygous (good cluster resolution with no-minor homozygous genotypes). Thus, a total of 151,509 SNPs were anchored to a unique location on the last version of Atlantic salmon genome assembly (GenBank: GCA_000233375.4) (Yáñez et al., 2016). We removed 149 markers deviating from Hardy-Weinberg equilibrium. After filtering for minor allele frequency, we obtained 61,199 loci with MAF > 5% for both Canadian populations joined (i.e., Can-W and Can-D), and 130,586 loci with MAF > 5% for Scottish populations joined (i.e., Sct-W and Sct-D). The 55,406 common markers among four populations that passed MAF filtering were used for the following analyses. We also discarded individuals with call rate less than 90%, and 165 individuals were retained for further analyses. The detailed number of samples for each population is shown in Table 1.

## 3.2 | Population differentiation and structure analyses

Using the set of markers with MAF > 5% common to all four populations revealed that measures of observed heterozygosity ($H_O$) and expected heterozygosity ($H_E$) were higher in Scottish than Canadian stocks. Overall, $H_O$ was slightly higher than $H_E$ in every population evaluated and heterozygosity levels were higher in both wild populations in comparison with their domesticated counterpart. Indeed, $H_O$ was 0.404 (CI = 0.400-0.407) versus 0.353 (CI = 0.350-0.356) and 0.297 (CI = 0.294-0.299) versus 0.266 (CI = 0.264-0.269) in wild versus domesticated Scottish and Canadian populations, respectively. A similar trend was observed for $H_E$ with wild Scottish exhibited the highest level ($H_E$ = 0.397, CI = 0.394-0.400), followed by domestic Scottish ($H_E$ = 0.343, CI = 0.340-0.346), while wild and domestic Canadian $H_E$ was lower ($H_E$ = 0.294, CI = 0.292: 0.297 and 0.244, CI = 0.242-0.247 respectively). All these genetic diversity measures were statistically significant from each other ($p$ < 0.05, Kruskal-Wallis test).

The highest level of differentiation was found between both domestic populations originating from Canada and Scotland ($F_{ST}$ = 0.393, CI = 0.390-0.396), although the $F_{ST}$ between both wild populations was 0.286 (CI = 0.284-0.288). When combining wild and domestic individuals (of each geographical origin), Canadian and Scottish populations still showed a high level of genetic differentiation ($F_{ST}$ = 0.293, CI = 0.290-0.295). A lower $F_{ST}$ was observed between Canadian wild and domestic populations ($F_{ST}$ = 0.171, CI = 0.169-0.173) which was twice as large as the $F_{ST}$ estimated between wild and domestic Scottish populations ($F_{ST}$ = 0.080, CI = 0.078-0.081).

The PCA using 55,406 polymorphic SNPs revealed four distinct clusters corresponding to the wild Canadian (Can-W), wild Scottish (Sct-W), domesticated Canadian (Can-D), and domesticated Scottish populations (Sct-D; Figure 1a). Principal components (PCs) 1 and 2 jointly accounted for 29.6% of total variance, with PC 1 (23.8%) separating Canadian and Scottish populations. Then, PC 2 (5.8%) discriminated Sct-W, Sct-D, and Can-W populations as well as some individual clusters existing within the

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** Population genetic structure. (a) Principal components analysis (PCA) of genetic differentiation among the 165 individuals, based on 55,406 polymorphic markers. Each point represents one individual, and different colors represent populations: Canadian domestic = dark green, Canadian wild = light green, Scottish domestic = dark blue, and Scottish wild = light blue. (b) Admixture analysis of the Atlantic salmon populations. The analysis was performed using the program ADMIXTURE with K = 4, based on 55,406 polymorphic markers. Each color represents a different theoretical ancestral population and each vertical line represent a single individual</figcaption></figure>

Can-W population. Regarding this axis, Can-D individuals formed a more heterogeneous group as some individuals are geometrically spread along PC 2, with even three individuals being distant from the main cluster. In agreement with PCA results, ADMIXTURE revealed the same four clusters (Figure 1b), based on the lowest cross-validation error. As ADMIXTURE analysis does not account for the bias generated by MAF filtering, we conducted this analysis using all SNPs, finding four distinctive clusters corresponding to the four original populations, which is congruent with the results including MAF filtering.

## 3.3 | Genomic markers putatively under selection distributed across the genome

### 3.3.1 | Bayesian $F_{ST}$ method

In Canadian populations, Bayesian analysis detected 160 SNPs putatively under divergent selection (Figure 2a). These markers were distributed among almost all the chromosomes (Figure 3a). The largest
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** Venn diagrams showing shared SNPs with evidence of selection among three independent tests in Canada populations (a) and Scottish populations (b). Also shown is the overlap of identified loci between both Canada and Scotland across all three tests (BayeScan, LFMM, Rsb [c])</figcaption></figure>

chromosome, Ssa01, showed 42 loci above the significance threshold, of which 22 markers had $-\log_{10} \ge 2$, meaning that these markers are potentially decisive evidence of selection. This chromosome presented the highest number of outliers concentrated within two regions, one of 370 Kb (13 SNPs) and the second one of 2,672 Kb (26 SNPs), respectively. Chromosomes Ssa03 and Ssa04 also contained, respectively, 15 and 27 loci potentially under divergent selection, the majority with a $-\log_{10} \ge 2$. Other chromosomes with several outlier markers were Ssa13, Ssa15, and Ssa21, with 13, 10, and 15 loci, respectively.

In Scottish populations, Bayesian scan detected 197 outliers. SNP distributed along the vast majority of the chromosomes, similarly to what we observed for Canadian populations (Figure 2b). Yet, in contrast to the analysis in Canadian populations, outliers were distributed in larger regions and were more evenly spaced within a given chromosome (Figure 3d). For instance, chromosome Ssa01 presented 30 outliers located in spanning regions of only 11 to 6,689 Kb long. The chromosome Ssa25 showed the highest number of outliers with 38 in total. Up to 20 outliers were also found in Ssa01, Ssa14, and Ssa18 chromosomes. A summary of Bayesian scan results for both populations is shown in Table 2.

### 3.3.2 | LFMM method

We considered a $|z|$ scores higher than 10 as the threshold to select SNPs associated with domestication using this approach. This cutoff value indicated significant SNP effect at the level of $p$-value < 10e-7 after applying a standard Bonferroni correction for $\alpha = 0.05$. We found 115 outliers in Canadian populations (Figure 2a) and nine outliers in Scottish populations (Figure 2b). In Canadian populations, these outliers were distributed in all chromosomes, except for Ssa14 chromosome (Figure 3b). The most representative chromosomes were Ssa01, Ssa03, and Ssa04 with 15, 11, and 14 outlier markers, respectively. In Scottish populations, the nine outliers found were distributed in eight chromosomes with Ssa16 being the only chromosomes harboring two outliers (Figure 3e). Details of these results are shown in Table 3.

### 3.3.3 | Rsb method

The Rsb analysis between Canadian wild and domestic populations revealed 79 SNPs under selection (Figure 2a), distributed on 13 chromosomes (Figure 3c). The highest scores were found in Ssa03 where 21 SNPs passed the threshold of $-\log_{10}(p\text{-value})$ equal to three within a
<figure><img src="image_5.png" /></figure>

**FIGURE 3** Manhattan plots showing the results of the BayeScan, LFMM, and Rsb outlier analyses applied across the domesticated strain/wild population pairs for Canada (Can; a-c) and Scotland (Sct; d-f). Salmon linkage groups were mapped are indicated on the x-axis. The dashed red line indicates the significant threshold applied to detect the outlier SNPs from the entire dataset. In BAYESCAN analysis this line represents the *q*-value threshold at $log_{10}$(0.05) = 1.32, which corresponds to Bayes factor of 10. In LFMM, the threshold was set at *z*-score = 10, while *Rsb* significant threshold was set at -$log_{10}$ one-tailed *p*-value = 3. Green dots in Canadian populations are the SNP common between BAYESCAN and LFMM analysis. In the case of Scottish population, blue dots are loci shared between BAYESCAN and Rsb
<table><caption>TABLE 2 Results of Bayesian analysis showing SNPs identified by chromosome</caption><thead><tr><th rowspan="2">Chromosome</th><th colspan="5">Can-D/Can-W</th><th colspan="5">Sct-D/Sct-W</th></tr><tr><th>Initial SNPs number</th><th>Strong</th><th>Very strong</th><th>Decisive</th><th>Total</th><th>Initial SNPs number</th><th>Strong</th><th>Very strong</th><th>Decisive</th><th>Total</th></tr></thead><tbody><tr><td>Ssa01</td><td>4,649</td><td>7</td><td>13</td><td>22</td><td>42</td><td>10,510</td><td>13</td><td>8</td><td>9</td><td>30</td></tr><tr><td>Ssa02</td><td>1,687</td><td></td><td></td><td></td><td></td><td>3,341</td><td></td><td>2</td><td></td><td>2</td></tr><tr><td>Ssa03</td><td>1,821</td><td>4</td><td>3</td><td>8</td><td>15</td><td>5,571</td><td></td><td></td><td>2</td><td>2</td></tr><tr><td>Ssa04</td><td>1,839</td><td>2</td><td>5</td><td>20</td><td>27</td><td>3,954</td><td>4</td><td>5</td><td>1</td><td>10</td></tr><tr><td>Ssa05</td><td>2,658</td><td></td><td>1</td><td>1</td><td>2</td><td>5,237</td><td></td><td></td><td>1</td><td>1</td></tr><tr><td>Ssa06</td><td>1,759</td><td></td><td>1</td><td></td><td>1</td><td>4,672</td><td></td><td></td><td></td><td></td></tr><tr><td>Ssa07</td><td>1,448</td><td>1</td><td></td><td></td><td>1</td><td>3,836</td><td></td><td></td><td></td><td></td></tr><tr><td>Ssa08</td><td>525</td><td></td><td></td><td>2</td><td>2</td><td>1,243</td><td></td><td></td><td></td><td></td></tr><tr><td>Ssa09</td><td>4,166</td><td>2</td><td>1</td><td></td><td>3</td><td>7,981</td><td>3</td><td>2</td><td>2</td><td>7</td></tr><tr><td>Ssa10</td><td>3,626</td><td>1</td><td></td><td></td><td>1</td><td>6,882</td><td>2</td><td>5</td><td>2</td><td>9</td></tr><tr><td>Ssa11</td><td>2,590</td><td></td><td>2</td><td>1</td><td>3</td><td>5,400</td><td>4</td><td>4</td><td>1</td><td>9</td></tr><tr><td>Ssa12</td><td>3,027</td><td></td><td></td><td></td><td></td><td>5,440</td><td>2</td><td>1</td><td>2</td><td>5</td></tr><tr><td>Ssa13</td><td>2,902</td><td>2</td><td>6</td><td>5</td><td>13</td><td>5,811</td><td>1</td><td>2</td><td>1</td><td>4</td></tr><tr><td>Ssa14</td><td>2,683</td><td>1</td><td>1</td><td></td><td>2</td><td>5,583</td><td>8</td><td>5</td><td>2</td><td>15</td></tr><tr><td>Ssa15</td><td>2,674</td><td>1</td><td>4</td><td>5</td><td>10</td><td>6,150</td><td></td><td></td><td>1</td><td>1</td></tr><tr><td>Ssa16</td><td>2,547</td><td>2</td><td></td><td></td><td>2</td><td>4,669</td><td>4</td><td>5</td><td>1</td><td>10</td></tr><tr><td>Ssa17</td><td>1,501</td><td></td><td></td><td>1</td><td>1</td><td>2,696</td><td></td><td>3</td><td></td><td>3</td></tr><tr><td>Ssa18</td><td>1,954</td><td>4</td><td>3</td><td></td><td>7</td><td>4,037</td><td>5</td><td>8</td><td>3</td><td>16</td></tr><tr><td>Ssa19</td><td>2,269</td><td></td><td></td><td></td><td></td><td>4,625</td><td>1</td><td></td><td></td><td>1</td></tr><tr><td>Ssa20</td><td>2,581</td><td></td><td>3</td><td>1</td><td>4</td><td>5,022</td><td>2</td><td>3</td><td></td><td>5</td></tr><tr><td>Ssa21</td><td>1,172</td><td>2</td><td>9</td><td>4</td><td>15</td><td>2,894</td><td></td><td>2</td><td>4</td><td>6</td></tr><tr><td>Ssa22</td><td>2,076</td><td></td><td></td><td></td><td></td><td>3,921</td><td></td><td></td><td></td><td></td></tr><tr><td>Ssa23</td><td>1,369</td><td>1</td><td>1</td><td></td><td>2</td><td>2,995</td><td>1</td><td>1</td><td>2</td><td>4</td></tr><tr><td>Ssa24</td><td>1,384</td><td></td><td>5</td><td></td><td>5</td><td>2,982</td><td>1</td><td>1</td><td>1</td><td>3</td></tr><tr><td>Ssa25</td><td>1,569</td><td>1</td><td></td><td></td><td>1</td><td>3,316</td><td>9</td><td>13</td><td>16</td><td>38</td></tr><tr><td>Ssa26</td><td>765</td><td></td><td></td><td></td><td></td><td>2,003</td><td>2</td><td></td><td></td><td>2</td></tr><tr><td>Ssa27</td><td>1,014</td><td></td><td></td><td></td><td></td><td>2,873</td><td>2</td><td>1</td><td>1</td><td>4</td></tr><tr><td>Ssa28</td><td>744</td><td></td><td></td><td>1</td><td>1</td><td>2,042</td><td></td><td>5</td><td></td><td>5</td></tr><tr><td>Ssa29</td><td>1,307</td><td></td><td></td><td></td><td></td><td>2,982</td><td></td><td>1</td><td></td><td>1</td></tr><tr><td>Unknown</td><td>893</td><td></td><td></td><td></td><td></td><td>1,918</td><td></td><td></td><td></td><td>4</td></tr><tr><td>Total</td><td>61,199</td><td>31</td><td>58</td><td>71</td><td>160</td><td>130,586</td><td>64</td><td>79</td><td>54</td><td>197</td></tr></tbody></table>

region of 1,269 Kb long. The region with the second highest score was on Ssa27 with 14 significant markers inside a 280 Kb region. For the comparison between wild and domestic in Scottish populations, Rsb detected 67 SNPs under selection (Figure 3f). The most representative chromosome was Ssa04 with 53 markers located within a region of 26,190 Kb long. Other significant SNPs are detailed in Table 3.

## 3.4 | Genomic regions putatively under selection distributed across the chromosomes

Overall, kernel-smoothing $F_{ST}$ results did not reveal any obvious island of divergence in common between wild versus farmed salmon

in both part of the Atlantic (i.e., Canada and Scotland; Figure 4). Accordingly, outliers found in each population were distributed differently along the chromosomes (Figure 5). For instance, more than 24 outliers were found as belonging to Ssa03 chromosome in Canadian populations, whereas only one was located in Ssa03 for the Scottish population. Conversely, 25 outliers were detected on Ssa25 in the Scottish populations and only two for the Canadian ones. Ssa01 and Ssa04 were the only chromosomes showing a high number of outliers detected in both systems (Figure 5). Yet, the chromosomes that harbored the highest number of outliers were different for the Canadian populations and Scottish populations; Ssa01 and Ssa04 respectively.
<table><caption>**TABLE 3** Results of LFMM and Rsb
Analyses SNPs identified by chromosome</caption><thead><tr><th rowspan="2">Chromosome</th><th colspan="3">Can-D/Can-W</th><th colspan="3">Sct-D/Sct-W</th></tr><tr><th>Initial SNPs number</th><th>LFMM</th><th>Rsb</th><th>Initial SNPs number</th><th>LFMM</th><th>Rsb</th></tr></thead><tbody><tr><td>Ssa01</td><td>4,649</td><td>15</td><td>16</td><td>10,510</td><td></td><td>6</td></tr><tr><td>Ssa02</td><td>1,687</td><td>1</td><td></td><td>3,341</td><td></td><td></td></tr><tr><td>Ssa03</td><td>1,821</td><td>11</td><td>21</td><td>5,571</td><td></td><td></td></tr><tr><td>Ssa04</td><td>1,839</td><td>14</td><td>1</td><td>3,954</td><td>1</td><td>53</td></tr><tr><td>Ssa05</td><td>2,658</td><td>2</td><td>4</td><td>5,237</td><td></td><td></td></tr><tr><td>Ssa06</td><td>1,759</td><td>6</td><td>1</td><td>4,672</td><td></td><td>3</td></tr><tr><td>Ssa07</td><td>1,448</td><td>2</td><td></td><td>3,836</td><td></td><td></td></tr><tr><td>Ssa08</td><td>525</td><td>4</td><td></td><td>1,243</td><td></td><td></td></tr><tr><td>Ssa09</td><td>4,166</td><td>7</td><td>12</td><td>7,981</td><td>1</td><td></td></tr><tr><td>Ssa10</td><td>3,626</td><td>2</td><td></td><td>6,882</td><td></td><td></td></tr><tr><td>Ssa11</td><td>2,590</td><td>3</td><td></td><td>5,400</td><td></td><td></td></tr><tr><td>Ssa12</td><td>3,027</td><td>3</td><td>1</td><td>5,440</td><td></td><td>2</td></tr><tr><td>Ssa13</td><td>2,902</td><td>4</td><td></td><td>5,811</td><td>1</td><td></td></tr><tr><td>Ssa14</td><td>2,683</td><td></td><td>4</td><td>5,583</td><td>1</td><td></td></tr><tr><td>Ssa15</td><td>2,674</td><td>4</td><td></td><td>6,150</td><td></td><td></td></tr><tr><td>Ssa16</td><td>2,547</td><td>6</td><td></td><td>4,669</td><td>2</td><td></td></tr><tr><td>Ssa17</td><td>1,501</td><td>1</td><td></td><td>2,696</td><td></td><td></td></tr><tr><td>Ssa18</td><td>1,954</td><td>1</td><td>1</td><td>4,037</td><td></td><td></td></tr><tr><td>Ssa19</td><td>2,269</td><td>3</td><td></td><td>4,625</td><td>1</td><td></td></tr><tr><td>Ssa20</td><td>2,581</td><td>2</td><td></td><td>5,022</td><td>1</td><td></td></tr><tr><td>Ssa21</td><td>1,172</td><td>2</td><td></td><td>2,894</td><td></td><td></td></tr><tr><td>Ssa22</td><td>2,076</td><td>7</td><td>1</td><td>3,921</td><td></td><td></td></tr><tr><td>Ssa23</td><td>1,369</td><td>3</td><td>2</td><td>2,995</td><td></td><td></td></tr><tr><td>Ssa24</td><td>1,384</td><td>1</td><td>1</td><td>2,982</td><td></td><td>2</td></tr><tr><td>Ssa25</td><td>1,569</td><td>3</td><td></td><td>3,316</td><td></td><td>1</td></tr><tr><td>Ssa26</td><td>765</td><td>2</td><td></td><td>2,003</td><td></td><td></td></tr><tr><td>Ssa27</td><td>1,014</td><td>2</td><td>14</td><td>2,873</td><td>1</td><td></td></tr><tr><td>Ssa28</td><td>744</td><td>1</td><td></td><td>2,042</td><td></td><td></td></tr><tr><td>Ssa29</td><td>1,307</td><td>1</td><td></td><td>2,982</td><td></td><td></td></tr><tr><td>Unknown</td><td>893</td><td>2</td><td></td><td>1,918</td><td></td><td></td></tr><tr><td>Total</td><td>61,199</td><td>115</td><td>79</td><td>130,586</td><td>9</td><td>67</td></tr></tbody></table>

## 3.5 | Genomic footprints of selection: Inconsistencies among methods

We uncovered a total of 337 and 270 SNPs putatively under divergent selection for Canadian and Scottish populations, respectively, according to the three complementary methods (Figure 2c). Overlaps are summarized in Venn diagrams in Figure 2c. In Canadian populations, only Bayesian and LFMM approaches jointly identified 17 outlier SNPs (Figures 2a and 3), whereas for Scottish population, only Rsb and Bayesian approaches found three outlier SNPs in common (Figures 2b and 3). Among these 17 shared outliers between Bayesian and LFMM in Canadian populations, seven regions were found to belong to seven genes (*gdnf*; *opn4x1b2*; *svil-like*; *tbck-like*; *plxnb2-like* and two uncharacterized proteins) in Atlantic salmon

genome annotation. In Scottish populations, three shared SNPs belong to two predicted genes (*brwd3-like* and *sh3bgrl3-like*). Only one common outlier on chromosome 25 was shared between the two geographical origins (Figure 2c), even when we considered all the 337 and 270 outliers detected. This common outlier (Affx-87919237) is located in the gene coding for ubiquitin-conjugating enzyme E2F putative (*ube2f*).

## 3.6 | Investigating for parallel signal of selection among the North American and European lineages

Overall, we found a total of 138 and 121 genes among the 337 and 270 outlier sequences detected as being potential targets of selection between wild versus domesticated Atlantic
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

**FIGURE 4** Kernel-smoothing *F*<sub>ST</sub>
calculated in 500-kb with steps of
50 kb. (a) Pattern of heterozygosity in
domestic (blue) and wild (green) Canadian
populations. (b) Pattern of *F*<sub>ST</sub> between
domestic/wild Canadian populations. (c)
Pattern of heterozygosity in domestic
(blue) and wild (green) Scottish
populations. (d) Pattern of *F*<sub>ST</sub> between
domestic/wild Scottish populations

**FIGURE 5** Distribution of the number
total of outliers along the Atlantic salmon
chromosome considering Canadian (CAN
in black) and Scottish population (SCT)

salmon, for Canada and Scotland, respectively. Among these
sets of genes, we discovered four genes overlapping continents:
ubiquitin-conjugating enzyme E2F putative (ube2f), collagen
alpha-1XIII chain (coda1), autism susceptibility 2 protein-like
(LOC106585083), and transient receptor potential cation chan-
nel subfamily M member 3-like (LOC106580056). A list of the

most important genes is described in Tables 4 and 5 for Canadian
and Scottish populations, respectively, while the complete list of
genes is available in Supporting Information Tables S1 and S2. A
gene ontology (GO) enrichment analysis on the total of outlier se-
quences found in both populations revealed that 11 and 12 GO
categories were significantly over-represented in Canadian and
<table><caption>**TABLE 4** Summary of genomic regions under selection and associated genes in Canadian populations</caption><thead><tr><th>Chromosome</th><th>Position SNPs</th><th>Test</th><th>Gene</th><th>Protein name</th></tr></thead><tbody><tr><td>Ssa01</td><td>118333490</td><td>RSB</td><td>LOC106563782</td><td>PREDICTED: unconventional myosin-Vb-like</td></tr><tr><td>Ssa01</td><td>119321920</td><td>LFMM_BYS</td><td>opn4x1b2</td><td>Melanopsin</td></tr><tr><td>Ssa01</td><td>119472753</td><td>LFMM_BYS</td><td>gdnf</td><td>Glial cell-derived neurotrophic factor</td></tr><tr><td>Ssa01</td><td>120337869<br>120350430</td><td>BYS</td><td>LOC106564357</td><td>PREDICTED: polypeptide N-acetylgalactosaminyltransferase 9-like</td></tr><tr><td>Ssa03</td><td>35341493<br>35353272</td><td>RSB</td><td>LOC106600280</td><td>PREDICTED: matrix metalloproteinase-14-like</td></tr><tr><td>Ssa03</td><td>58772401</td><td>LFMM_BYS</td><td>LOC106600794</td><td>PREDICTED: supervillin-like</td></tr><tr><td>Ssa03</td><td>58814947</td><td>LFMM_BYS</td><td>LOC106600792</td><td>PREDICTED: uncharacterized protein LOC106600792</td></tr><tr><td>Ssa04</td><td>1988842</td><td>LFMM_BYS</td><td>LOC106602233</td><td>PREDICTED: uncharacterized protein LOC106602233</td></tr><tr><td>Ssa04</td><td>2491374<br>2503265</td><td>BYS</td><td>LOC106602228</td><td>PREDICTED: dual specificity testis-specific protein kinase 2-like</td></tr><tr><td>Ssa04</td><td>41701375</td><td>RSB</td><td>LOC106603279</td><td>PREDICTED: unconventional myosin-XVIIIa-like</td></tr><tr><td>Ssa04</td><td>54112595</td><td>LFMM</td><td>LOC106603471</td><td>PREDICTED: collagen alpha-1(XXVI) chain-like</td></tr><tr><td>Ssa06</td><td>37003521</td><td>LFMM</td><td>LOC106607315</td><td>PREDICTED: unconventional myosin-X-like</td></tr><tr><td>Ssa08</td><td>19722947</td><td>LFMM_BYS</td><td>LOC106610360</td><td>PREDICTED: TBC domain-containing protein kinase-like protein</td></tr><tr><td>Ssa17</td><td>54473681</td><td>LFMM_BYS</td><td>LOC106576598</td><td>PREDICTED: plexin-B2-like</td></tr><tr><td>Ssa18</td><td>13286009</td><td>RSB</td><td>coda1</td><td>Collagen alpha-1 XIII chain</td></tr><tr><td>Ssa18</td><td>62851227</td><td>BYS</td><td>LOC106577990</td><td>PREDICTED: E3 ubiquitin-protein ligase TRIM39-like</td></tr><tr><td>Ssa20</td><td>15597033</td><td>BYS</td><td>LOC106580056</td><td>PREDICTED: transient receptor potential cation channel subfamily M member 3-like</td></tr><tr><td>Ssa22</td><td>29365659<br>29365735</td><td>LFMM</td><td>LOC106583145</td><td>PREDICTED: collagen alpha-1(VII) chain-like</td></tr><tr><td>Ssa24</td><td>12213523</td><td>BYS</td><td>LOC106585083</td><td>PREDICTED: autism susceptibility gene 2 protein-like</td></tr><tr><td>Ssa25</td><td>19817836</td><td>LFMM</td><td>ube2f</td><td>Ubiquitin-conjugating enzyme E2F (putative)</td></tr><tr><td>Ssa27</td><td>10232307<br>10238834</td><td>RSB</td><td>brd2</td><td>Bromodomain containing 2</td></tr><tr><td>Ssa27</td><td>10251922</td><td>RSB</td><td>LOC106588396</td><td>PREDICTED: collagen alpha-2(XI) chain-like</td></tr></tbody></table>

Scottish populations, respectively. The most representative categories in both populations were cellular and metabolic processes (GO:0009987 and GO:0008152) including 38.8%–51.4% and 33.9%–34.7% of the genes, respectively. Other categories were biological regulation, developmental processes, localization, multicellular organismal processes, response to stimulus, biological adhesion, immune system processes, and reproduction. All these categories were shared in both populations, although locomotion category was only present in Scottish population (this includes 0.9% of the genes).

4 | DISCUSSION

Understanding how the process of domestication may shape the genome of wild to domesticated animals is particularly useful for enhancing our knowledge on how human-driven selection may induce genetic changes as it may provide a practical framework to guide genetic improvement practices (Wang, Xie, Peng, Irwin, & Zhang, 2014). Here, we performed a large genomewide scan, using a dense SNP array and three outlier methods, to investigate genomic footprints of domestication in two geographically isolated Atlantic salmon
<table><caption>TABLE 5 Summary of genomic regions under selection and associated genes in Scottish populations</caption><thead><tr><th>Chromosome</th><th>Position SNPs</th><th>Test</th><th>Gene</th><th>Protein name</th></tr></thead><tbody><tr><td>Ssa01</td><td>48495280<br>48506349</td><td>BYS</td><td>LOC106604517</td><td>PREDICTED: myopalladin-like</td></tr><tr><td>Ssa03</td><td>50571649</td><td>BYS</td><td>LOC106600700</td><td>PREDICTED: transient receptor potential<br>cation channel subfamily M member 4-like</td></tr><tr><td>Ssa04</td><td>44415367</td><td>BYS_RSB</td><td>LOC106603326</td><td>PREDICTED: SH3 domain-binding glutamic<br>acid-rich-like protein</td></tr><tr><td>Ssa04</td><td>44432033<br>44436726</td><td>BYS_RSB</td><td>LOC106603328</td><td>PREDICTED: bromodomain and WD<br>repeat-containing protein 3-like</td></tr><tr><td>Ssa04</td><td>45658450<br>45667448</td><td>RSB</td><td>LOC106603377</td><td>PREDICTED: long-chain fatty acid-CoA<br>ligase 4-like</td></tr><tr><td>Ssa04</td><td>47256953<br>47265274</td><td>RSB</td><td>LOC106603586</td><td>PREDICTED: E3 ubiquitin-protein ligase<br>Siah2</td></tr><tr><td>Ssa04</td><td>51300396<br>51426545</td><td>RSB</td><td>LOC106603511</td><td>PREDICTED: active breakpoint cluster<br>region-related protein-like</td></tr><tr><td>Ssa10</td><td>16255371<br>16284803</td><td>BYS</td><td>LOC106613666</td><td>PREDICTED: unconventional<br>myosin-IXb-like</td></tr><tr><td>Ssa14</td><td>47278390</td><td>BYS</td><td>tekt2</td><td>Tektin 2 (testicular)</td></tr><tr><td>Ssa14</td><td>65526017</td><td>BYS</td><td>LOC106570217</td><td>PREDICTED: collagen alpha-1(IX) chain-like</td></tr><tr><td>Ssa14</td><td>65541964<br>65552172</td><td>BYS</td><td>LOC106570216</td><td>PREDICTED: collagen alpha-1(XXII)<br>chain-like</td></tr><tr><td>Ssa14</td><td>79347457<br>79378912</td><td>BYS</td><td>LOC106570512</td><td>PREDICTED: collagen alpha-2(IX) chain-like</td></tr><tr><td>Ssa18</td><td>13282526<br>13283308</td><td>BYS</td><td>coda1</td><td>Collagen alpha-1XIII chain</td></tr><tr><td>Ssa18</td><td>47314912</td><td>BYS</td><td>LOC106577539</td><td>PREDICTED: neurexin-2-like</td></tr><tr><td>Ssa20</td><td>15586930</td><td>BYS</td><td>LOC106580056</td><td>PREDICTED: transient receptor potential<br>cation channel subfamily M member 3-like</td></tr><tr><td>Ssa24</td><td>12125145<br>12127325</td><td>RSB</td><td>LOC106585083</td><td>PREDICTED: autism susceptibility gene 2<br>protein-like</td></tr><tr><td>Ssa25</td><td>19817836</td><td>BYS</td><td>ube2f</td><td>Ubiquitin-conjugating enzyme E2F<br>(putative)</td></tr><tr><td>Ssa25</td><td>32313916<br>32325227</td><td>BYS</td><td>mspd2</td><td>Motile sperm domain-containing protein 2</td></tr><tr><td>Ssa27</td><td>18583143</td><td>BYS</td><td>LOC106588644</td><td>PREDICTED: zona pellucida sperm-binding<br>protein 3-like</td></tr></tbody></table>

strains. Overall, our results revealed only a few footprints of domestication that overlap among methods and populations, suggesting that detection of selection depends on the method used and parallel genomic changes resulting from domestication events are rare in this system. Indeed, domestication is a recent and ongoing process in Atlantic salmon populations, which may explain this low number of genomic changes. Furthermore, as Canadian and Scottish populations have diverged independently for possibly as much as 1 million years (Rougemont & Bernatchez, 2018), it is likely that this had led to different standing genetic variation upon which selection may be acting, limiting the potential for parallel evolution at the genome level in this study system. This is in agreement with previous population genomic work on Atlantic salmon (Mäkinen et al., 2014; Vasemägi et al., 2012) and, more broadly, on several studies showing that selection may act faster on standing variation compared to new mutations (Barrett &

Schluter, 2008). Limited evidence for parallel impact of artificial selection at the genome level contrasts with the pronounced pattern of parallelism previously documented at the transcriptome level (Roberge et al., 2006), suggesting the different genomic architecture may result in similar pattern of gene expression as well as parallel phenotypic changes occurring during domestication. Clearly, shedding light on the complete genomic basis of domestication of Atlantic salmon would require using alternative methods rather than genome scan and collecting a larger dataset with known phenotypes. Future studies testing for parallel evolution in Atlantic salmon should consider this avenue.

## 4.1 | Limitations of genome scan methods

Traits of interest selected in aquaculture programs, such as growth and fat content, may be polygenic, that is, controlled by many genes
with high variance of effect size (Gagnaire & Gaggiotti, 2016). Polygenic selection may imply a subtle change of allele frequencies at several loci that would be hard to detect using genome scan methods that mainly focus on moderate to high shift in allele frequency at few independent loci. For instance, Bourret, Dionne, and Bernatchez (2014) genotyped 5,568 SNPs to test for differential allelic and genotypic frequencies between juveniles (smolts) migrating to sea and adults (grilses) returning to freshwater after 1 year at sea. Although numerous outliers were identified by the single-locus analysis, no evidence for parallel, temporally repeated selection was found. In contrast, the multilocus approach detected repeated patterns of selection for a multilocus group of 34 co-varying SNPs. Overall, the large difference in the sampling design of wild populations between Canada and Scotland may partly explain why so few common outliers were highlighted by genome scan methods when comparing the two locations. The minor allele frequency and Hardy-Weinberg filtering could also have affected the final result of the analysis, leading to lower number of SNPs detected. However, Bayesian analysis conducted with no MAF and HWE filters still revealed none loci putatively selected in Canadian populations, although 104 markers (Supporting Information Table S3) were detected in Scottish population, which is a lower number than those detected using filtering options, revealing the suitableness of applying these filters to discard low-quality genotypes.

Alternative models to common genome scan approaches, such as a quantitative genetic framework, provide powerful tools to investigate polygenic selection. Earlier Atlantic salmon studies of quantitative trait loci (QTL) mapping have already reported a large number of different QTLs for growth (Baranski, Moen, & Våge, 2010; Gutierrez et al., 2012; Houston et al., 2009; Reid, Szanto, Glebe, Danzmann, & Ferguson, 2005), although genomewide studies revealed very low levels of association between markers and growth (Gutierrez, Yáñez, & Davidson, 2015; Gutierrez, Yáñez, Fukui, Swift, & Davidson, 2015). This indicates that growth-related traits are likely to be controlled by a number of population-specific loci of low to moderate effect with an important polygenic component (Tsai et al., 2015), making it difficult to compare our results with those obtained from GWAS or QTL mapping. Future studies combining population genomic and quantitative genetic framework would then provide a more inclusive approach for uncovering the molecular footprints of domestication in Atlantic salmon. Furthermore, the complete functional annotation of the Atlantic salmon genome will undoubtedly enhance opportunities for a comprehensive understanding of the effects of domestication in Atlantic salmon populations (Macqueen et al., 2017).

## 4.2 | Standing genetic variation and loss of genetic diversity

Of the 151,509 SNPs successfully genotyped, 61,199 and 130,586 SNPs were polymorphic for Canadian and Scottish populations, respectively. This striking difference on the number of polymorphic SNPs is in accordance with earlier studies based on various types of markers (mtDNA, microsatellites, AFLPs, SNPs) showing that North

American Atlantic salmon populations have lower genetic diversity than European populations (Bourret et al., 2013; Mäkinen et al., 2014). This may also be explained to some extent by ascertainment bias as the 200K SNP Array used in the present study was created based on genomic information mostly from European salmon (Yáñez et al., 2016). Both lineages only shared 55,406 SNPs in common and the genetic diversity analysis of these SNPs revealed that the lowest and highest levels of heterozygosity were found in Canadian domestic and Scottish wild population, respectively. Indeed, lower levels of genetic diversity in both domestic populations compared to their wild homologs were expected as most of domestic populations experience an inevitable loss of genetic diversity due to both selective breeding as well as the absence of genetic connectivity with other populations as found in nature (Baumung, Simianer, & Hoffmann, 2004; Johansson & Rendel, 1968). Concomitantly, and as expected from reduced diversity within population, domesticated populations showed a higher proportion of genetic differentiation ($F_{ST}$ = 0.393) than wild populations ($F_{ST}$ = 0.286). This outcome is in agreement with the work of Mäkinen et al. (2014) who reported that the mean value of genetic differentiation between two populations of hatcheries strain is higher than between wild populations. On the other hand, the low level of genetic differentiation ($F_{ST}$ = 0.171) between domestic and wild Scottish populations showed that using a wild population from a nearby geographical region as a genetic proxy was appropriate, in the absence of a true ancestral population.

## 4.3 | Detecting the impacts of domestication with genome scans

Using a genomewide SNP array, we found limited evidence of potential signatures of selection that were jointly detected between domestic and natural populations, in Canada and Scotland, considering the outcomes of three different outlier tests. This lack of congruence among methods is likely to result from the statistics underlying each approach: LFMM is based on linear associations, Bayesian is an $F_{ST}$ differentiation method, and Rsb is a haplotypes comparison method. As LFMM is based on linear model and it detects a few numbers of outliers in our study system, this approach may not be the most appropriate approach for this kind of experimental design. More particularly, Rsb test is likely to detect more recent selection signatures in comparison with $F_{ST}$-based methods (Oleksyk, Smith, & O'Brien, 2009). On the other hand, Bayesian was also the method that resulted in the highest number of outliers identified but it should be noted that we used a permissive threshold (PO = 10), as used by Gutierrez, Yáñez, and Davidson (2015); Chen et al. (2016); Lin et al. (2017) that may have led to a high rate of false positives (Lotterhos & Whitlock, 2015). This nonoverlap in outlier results tests has already been observed and discussed by Mäkinen et al. (2014) who also investigated genomic signal of domestication in Atlantic salmon using Bayesian and LnRH methods. As the three different tests can detect different patterns of molecular variation (de Simoni Gouveia, da Silva, Paiva, & de Oliveira, 2014) and at different time scales (Hohenlohe et al., 2010; Oleksyk et al.,
2009) even if SNPs are only detected by one test, they still can be
real targets of selection. Otherwise, it is important to mention that
chromosomal differences between European and North American
Atlantic salmon lineages can have an impact on analyses based on
haplotype measures (Rsb). According to the linkage map of 3K SNPs,
gross chromosomal re-arrangements between European and North
American Atlantic salmon have been identified involving Ssa01 (that
presents a translocation of its p arm to Ssa23) and two chromosomal
fusions: Ssa26 with Ssa28 and Ssa08 and Ssa29 (Brenna-Hansen
et al., 2012). Nevertheless, while these major structural differ-
ences were important, marker order within these chromosomal
regions was remarkably conserved (Brenna-Hansen et al., 2012).
Chromosomal re-arrangements may lead to spurious results when
using European Atlantic salmon as the reference genome sequence
in North American populations. In this study, the only chromosome
presenting re-arrangements in Canadian populations that harbors
loci putatively under selection was Ssa01 with 16 markers; conse-
quently, these results should be interpreted with caution.

## 4.4 | Investigating the existence of parallel genomic footprints of Atlantic salmon domestication

Over the past 50 years, intense artificial selection in Atlantic salmon
populations has led to the development of strains specialized in cer-
tain phenotypic traits which has caused large phenotypic and geno-
typic changes between wild and domestic populations (Glover et al.,
2009). Despite Atlantic salmon populations having been intensely
selected for growth-related traits on both sides of the Atlantic ocean,
we found little overlap between the outlier genomic regions identi-
fied in previous studies of other Atlantic salmon strains selected for
high growth rate (Gutierrez, Yáñez, & Davidson, 2015; Gutierrez,
Yáñez, Fukui, et al., 2015). These population-specific signals sug-
gest that selection may have acted upon different genes, which was
already shown by one study documenting genomewide footprints
of pig domestication (Amaral et al., 2011). More particularly, several
population genomic studies of Atlantic salmon populations have sug-
gested that the same phenotype may arise from different genetic
pathways among geographically isolated populations (Elmer et al.,
2014; Mäkinen et al., 2014; Perrier, Bourret, Kent, & Bernatchez,
2013; Pujolar, Ferchaud, Bekkevold, & Hansen, 2017; Vasemägi
et al., 2012). Evolution of complex parallel phenotypes can indeed
arise from different evolutionary routes and this is likely to happen
when inbreeding and genetic drift play a greater role than selection.

## 4.5 | Parallel genomic regions detected as being putatively under artificial selection

Despite limited evidence for parallelism at the genome level, evi-
dence for parallel evolution was observed at a few potentially im-
portant genes. Here, genomic regions harboring collagen alpha-1XIII
chain gene (coda1) were identified in both comparisons. There is
evidence that collagen genes may be involved in the response of
Atlantic salmon to sea lice as initiators of inflammatory cytokine

signaling (Castillo-Briceño et al., 2009; Correa et al., 2016; Krasnov,
Skugor, Todorcevic, Glover, & Nilsen, 2012). Therefore, this gene
might be involved in the immune response to cope with specific
diseases present in aquaculture environment and may be a relevant
target of selection. Despite none of the farmed populations having
been directly selected for traits associated with disease resistance,
farmed populations are often subjected to high levels of patho-
gens, such as parasite infections, which are known to be among the
strongest selective forces driving the evolution of host populations
(Roberge et al., 2006; Zueva et al., 2014). A second common gene
was ubiquitin-conjugating enzyme E2F putative (ube2f), coding for
ubiquitin-conjugating enzyme E2-F. Growth and development of
skeletal muscle undergo breakdown and replacement of proteins
during different periods, and in teleost fish, this process involves E2-
ubiquitin-conjugating enzymes and E3-ubiquitin ligases (Johnston,
Bower, & Macqueen, 2011). Therefore, this gene could be involved
in muscle development and growth in Atlantic salmon, which is an
important trait selected for enhancing salmon aquaculture pro-
duction. This gene family was also suggested to be potentially in-
volved in immune response of marine vertebrates (Núñez-Acuña,
Aguilar-Espinoza, Chávez-Mardones, & Gallardo-Escárate, 2012),
and more particularly, ubiquitin genes have already been detected
to be potential targets of selection in the domestication process
of shrimp (Rotllant et al., 2015), auroch (Braud et al., 2017) and
rice (Li & Li, 2014). Finally, the autism susceptibility 2 protein-like
(LOC106585083) was also found to overlap between continents.
This gene (Auts2-like) has been associated with autism and mental
retardation in humans (Bedogni et al., 2010; Oksenberg & Ahituv,
2013). Interestingly, Auts2 has been found to be under selection in
domestic cattle breeds (Consortium, Bovine HapMap, 2009), sug-
gesting that domestication may act on behavioral traits in salmon as
in other domestic animals (Clutton-Brock, 1999). Finally, a potential
parallel signal of selection was located in the transient receptor po-
tential cation channel subfamily M member 3-like (LOC106580056),
which belongs to Trpm3, a gene associated with the reception of nox-
ious temperature stimuli in mammals (Vriens et al., 2011). Transient
receptor potential cation channel subfamily genes are involved in
melanocyte pigmentation (Cieslak, Reissmann, Hofreiter, & Ludwig,
2011). Therefore, this gene may play a role in the immune system
of domesticated Atlantic salmon and Kittilsen, Johansen, Braastad,
and Øverli (2012) showed that pigmentation seems to be correlated
with the development of the ectoparasitic lice in Atlantic salmon.
Interestingly, Yang, Li, Li, Fan, and Tang (2014) also detected selec-
tion signature linked to domestication process in this gene by com-
paring Chinese indigenous and commercial pig breeds.

## 4.6 | Biological function of nonparallel genomic regions identified as putatively being under artificial selection

Some of the genes that overlap between statistical methods in
Canadian populations (Bayescan and LFMM) included supervillin-like
gene, which is a protein involved in actin and myosin II assembly that
promotes cell growth (Fang & Luna, 2013). This gene has been shown to be involved in muscle fiber type determination in large white pigs (Zhu et al., 2016). We suggest that this gene may be associated with muscle development in Atlantic salmon and be under selection for improving growth-related traits. We also identified *plexin-b2-like*, which in mammals is involved in several processes during development of the nervous, cardiovascular, renal, and skeletal system (Worzfeld et al., 2014). Another plexin from the same subfamily has been associated with behavior-related traits, such as tameness and temperament in rat and cattle (Friedrich, Brand, & Schwerin, 2015; Heyne et al., 2014). This finding is in line with the identification of the gene autism susceptibility 2 protein-like and both suggest that domestication could be acting on behavior-related traits. On the other hand, *opn4x1b2* detected in Canadian populations is a melanopsin that can be involved in signaling of photoperiodic information in Atlantic salmon (Sandbakken, Ebbesson, Stefansson, & Helvik, 2012), indicating some role in the adaptation to the aquaculture environment.

We detected a potential gene candidate in Scottish populations, *brwd3*, that has been associated to mental retardation in humans (Field et al., 2007) and which was also identified as a target of selection during the domestication process in cattle (Consortium, Bovine HapMap, 2009). Other genes detected in this study included *myopal-ladin-like* gene in Scottish populations, which has been associated with late maturation in Atlantic salmon (Gutierrez et al., 2014). The long-chain fatty acid-CoA ligase 4-like gene, located in Ssa04, which belongs to Acsl family of genes, was found as a potential gene influenced by domestication in Scottish populations. This gene is involved in lipid metabolism, although the activation of long-chain fatty acids for synthesis and degradations of cellular lipids (Golej et al., 2011). In rainbow trout, lipid metabolism has been reported to be associated with growth (Xu et al., 2011). Similarly, another gene from the same family, the long-chain fatty acid-CoA ligase 1 has been shown to be associated with growth in clam (*Meretrix meretrix*, Dai, Huan, Wang, Xia, & Liu, 2015).

No other overlaps were found with previous studies related to sexual maturity in Atlantic salmon, for both Canadian and Scottish populations. Nevertheless, we found putative outliers in the bromodomain-containing two genes (*br2*) and the zona pellucida (ZP) sperm-binding protein 3-like gene (*Zp3*) in Canadian and Scottish populations, respectively, both localized in chromosome Ssa27. The bromodomain-containing 2 gene may be involved in spermatogenesis or folliculogenesis (Rhee, Brunori, Besset, Trousdale, & Wolgemuth, 1998), while ZPs are responsible for the initial sperm binding and the subsequent induction of the acrosome reaction that allows sperm penetration in mammals (Lin, Roy, Yan, Burns, & Matzuk, 2007). Thus, these genes are interesting functional candidates for reproduction-related traits under selection on these populations. To better understand the molecular functions of these genes, we investigated their GO classification. Many genes were categorized into cellular process and metabolic process. Both categories can involve cell growth or anabolic/catabolic process resulting in cell growth, which suggest that these genes could be associated with growth in Atlantic salmon. Admittedly, these observations need to be explored and verified in further studies.

## 5 | CONCLUSION

Domestication events in Atlantic salmon have produced pronounced phenotypic changes in domesticated fish compared to their wild counterpart, and the selection process has left important signals throughout the genome of these populations. This study highlighted the existence of nonparallel and parallel signatures of selection between domesticated versus farmed fish, detected by different statistical approaches. A small number of parallel genomic regions putatively under selection belong to genes with molecular functions that might be associated with traits under domestication, such as behavior, immune response, and reproduction. However, the vast majority of genomic regions detected as being putatively under selection are not shared either among method or lineages, suggesting that domestication may involve a different set of genes, perhaps because of the distinct evolutionary history and very different standing genetic variation of Canadian and Scottish populations. Another possibility is that alternative methods to genome scan that test for genomic signal of polygenic adaptation are required for revealing genetic architecture underlying trait of interest. Indeed, the most important trait for which Atlantic salmon has been artificially selected is growth, whereas few selection signatures possibly associated with this trait have been detected, perhaps due to the polygenic nature of this trait. Nevertheless, our results contribute to the detection of variants that may underlie important traits of practical interest for aquaculture and to fully understanding the early patterns and effect of domestication in Atlantic salmon.

## ACKNOWLEDGEMENTS

Financial support was provided by CORFO (11IEI-12843 and 12PIE17669). We are grateful to Doctoral scholarship from CONICYT (21120382), Government of Chile; Canada-Chile Leadership Exchange Scholarship program, Government of Canada. We acknowledge to Carlos Soto for his contribution with samples from Scottish domestic population. We also thank the Associate Editor and two anonymous reviewers for or their constructive comments and suggestions on the manuscript. JMY is supported by Núcleo Milenio INVASAL funded by Chile's government program, Iniciativa Científica Milenio from Ministerio de Economía, Fomento y Turismo.

## CONFLICT OF INTEREST

None declared.

## DATA ARCHIVING

Data available from the Dryad Digital Repository: https://doi.org/10.5061/dryad.60b9p56.

## ORCID

José M. Yáñez http://orcid.org/0000-0002-6612-4087
REFERENCES

Alexander, D. H., Novembre, J., & Lange, K. (2009). Fast model-based estimation of ancestry in unrelated individuals. *Genome Research*, 19(9), 1655-1664. https://doi.org/10.1101/gr.094052.109

Amaral, A. J., Ferretti, L., Megens, H.-J., Crooijmans, R. P. M. A., Nie, H., Ramos-Onsins, S. E., ... Groenen, M. A. M. (2011). Genome-wide footprints of pig domestication and selection revealed through massive parallel sequencing of pooled DNA. *PLoS ONE*, 6(4), e14782. https://doi.org/10.1371/journal.pone.0014782

Andersson, L. (2012). How selective sweeps in domestic animals provide new insight into biological mechanisms. *Journal of Internal Medicine*, 271(1), 1-14. https://doi.org/10.1111/j.1365-2796.2011.02450.x

Andersson, L. (2013). Molecular consequences of animal breeding. *Current Opinion in Genetics & Development*, 23(3), 295-301. https://doi.org/10.1016/j.gde.2013.02.014

Balon, E. K. (2004). About the oldest domesticates among fishes. *Journal of Fish Biology*, 65, 1-27. https://doi.org/10.1111/j.0022-1112.2004.00563.x

Baranski, M., Moen, T., & Våge, D. I. (2010). Mapping of quantitative trait loci for flesh colour and growth traits in Atlantic salmon (*Salmo salar*). *Genetics Selection Evolution*, 42(1), 1-14.

Barrett, R. D. H., & Schluter, D. (2008). Adaptation from standing genetic variation. *Trends in Ecology & Evolution*, 23(1), 38-44. https://doi.org/10.1016/j.tree.2007.09.008

Baumung, R., Simianer, H., & Hoffmann, I. (2004). Genetic diversity studies in farm animals—A survey. *Journal of Animal Breeding and Genetics*, 121(6), 361-373. https://doi.org/10.1111/j.1439-0388.2004.00479.x

Bedogni, F., Hodge, R. D., Nelson, B. R., Frederick, E. A., Shiba, N., Daza, R. A., & Hevner, R. F. (2010). Autism susceptibility candidate 2 (Auts2) encodes a nuclear protein expressed in developing brain regions implicated in autism neuropathology. *Gene Expression Patterns: GEP*, 10(1), 9. https://doi.org/10.1016/j.gep.2009.11.005

Benestan, L., Quinn Brady, K., Maaroufi, H., Martin Laporte, K., Clark Fraser, J., Spencer, G., ... Bernatchez, L. (2016). Seascape genomics provides evidence for thermal adaptation and current-mediated population structure in American lobster (*Homarus americanus*). *Molecular Ecology*, 25(20), 5073-5092. https://doi.org/10.1111/mec.13811

Bourret, V., Dionne, M., & Bernatchez, L. (2014). Detecting genotypic changes associated with selective mortality at sea in Atlantic salmon: Polygenic multilocus analysis surpasses genome scan. *Molecular Ecology*, 23(18), 4444-4457. https://doi.org/10.1111/mec.12798

Bourret, V., Kent, M. P., Primmer, C. R., Vasemagi, A., Karlsson, S., Hindar, K., ... Lien, S. (2013). SNP-array reveals genome wide patterns of geographical and potential adaptive divergence across the natural range of Atlantic salmon (*Salmo salar*). *Molecular Ecology*, 22(3), 532-551. https://doi.org/10.1111/mec.12003

Bourret, V., O'Reilly, P. T., Carr, J. W., Berg, P. R., & Bernatchez, L. (2011). Temporal change in genetic integrity suggests loss of local adaptation in a wild Atlantic salmon (*Salmo salar*) population following introduction by farmed escapees. *Heredity*, 106(3), 500-510. https://doi.org/10.1038/hdy.2010.165

Braud, M., Magee, D. A., Park, S. D. E., Sonstegard, T. S., Waters, S. M., MacHugh, D. E., & Spillane, C. (2017). Genome-wide microRNA binding site variation between extinct wild aurochs and modern cattle identifies candidate microRNA-regulated domestication genes. *Frontiers in Genetics*, 8, 3.

Brenna-Hansen, S., Li, J., Kent, M. P., Boulding, E. G., Dominik, S., Davidson, W. S., & Lien, S. (2012). Chromosomal differences between European and North American Atlantic salmon discovered by linkage mapping and supported by fluorescence in situ hybridization analysis. *BMC Genomics*, 13(1), 1-13.

Brito, L. F., Kijas, J. W., Ventura, R. V., Sargolzaei, M., Porto-Neto, L. R., Cánovas, A., ... Schenkel, F. S. (2017). Genetic

diversity and signatures of selection in various goat breeds revealed by genome-wide SNP markers. *BMC Genomics*, 18(1), 229. https://doi.org/10.1186/s12864-017-3610-0

Browning, B. L., & Browning, S. R. (2009). A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. *American Journal of Human Genetics*, 84(2), 210-223. https://doi.org/10.1016/j.ajhg.2009.01.005

Canty, A., & Ripley, B. (2016). *boot: Bootstrap R (S-Plus) Functions*. Paquete de R, versión 1.3-18. Retrieved from https://CRAN.R-project.org/package=boot

Castillo-Briceño, P., Sepulcre, M. P., Chaves-Pozo, E., Meseguer, J., García-Ayala, A., & Mulero, V. (2009). Collagen regulates the activation of professional phagocytes of the teleost fish gilthead seabream. *Molecular Immunology*, 46(7), 1409-1415. https://doi.org/10.1016/j.molimm.2008.12.005

Cesconeto, R. J., Joost, S., McManus, C. M., Paiva, S. R., Cobuci, J. A., & Braccini, J. (2017). Landscape genomic approach to detect selection signatures in locally adapted Brazilian swine genetic groups. *Ecology and Evolution*, 7(22), 9544-9556. https://doi.org/10.1002/eece.3323

Chen, J., Källman, T., Ma, X.-F., Zaina, G., Morgante, M., & Lascoux, M. (2016). Identifying genetic signatures of natural selection using pooled population sequencing in Picea abies. *G3: Genes|Genomes|Genetics*, 6(7), 1979-1989. https://doi.org/10.1534/g3.116.028753

Cieslak, M., Reissmann, M., Hofreiter, M., & Ludwig, A. (2011). Colours of domestication. *Biological Reviews*, 86(4), 885-899. https://doi.org/10.1111/j.1469-185X.2011.00177.x

Cingolani, P., Platts, A., Wang, L., Coon, M., Nguyen, T., Land, S. J., ... Ruden, D. M. (2012). A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of *Drosophila melanogaster* strain w1118; iso-2; iso-3. Fly, 6(2), 80-92. https://doi.org/10.4161/fly.19695

Clutton-Brock, J. (1999). A natural history of domesticated mammals. Cambridge, UK: Cambridge University Press.

Consortium, Bovine HapMap (2009). Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds. *Science*, 324(5926), 528-532.

Correa, K., Lhorente, J. P., Bassini, L., López, M. E., Di Genova, A., Maass, A., ... Yáñez, J. M. (2016). Genome wide association study for resistance to Caligus rogercresseyi in Atlantic salmon (*Salmo salar* L.) using a 50K SNP genotyping array. *Aquaculture*, 472, 61-65.

Dai, P., Huan, P., Wang, H., Xia, L., & Liu, B. (2015). Characterization of a long-chain fatty acid-CoA ligase 1 gene and association between its SNPs and growth traits in the clam *Meretrix meretrix*. *Gene*, 566(2), 194-200. https://doi.org/10.1016/j.gene.2015.04.047

Danecek, P., Auton, A., Abecasis, G., Albers, C. A., Banks, E., DePristo, M. A., ... 1000 Genomes Project Analysis Group (2011). The variant call format and VCFtools. *Bioinformatics*, 27, 2156-2158. https://doi.org/10.1093/bioinformatics/btr330

de Simoni Gouveia, J. J., da Silva, M. V. G. B., Paiva, S. R., & de Oliveira, S. M. P. (2014). Identification of selection signatures in livestock species. *Genetics and Molecular Biology*, 37(2), 330-342. https://doi.org/10.1590/S1415-47572014000300004

Devlin, R. H., Sakhrani, D., Tymchuk, W. E., Rise, M. L., & Goh, B. (2009). Domestication and growth hormone transgenesis cause similar changes in gene expression in coho salmon (*Oncorhynchus kisutch*). *Proceedings of the National Academy of Sciences of the United States of America*, 106(9), 3047-3052. https://doi.org/10.1073/pnas.0809798106

Dillane, E., McGinnity, P., Coughlan, J. P., Cross, M. C., De Eyto, E., Kenchington, E., ... Cross, T. F. (2008). Demographics and landscape features determine intrariver population structure in Atlantic salmon (*Salmo salar* L.): The case of the River Moy in
Ireland. *Molecular Ecology*, 17(22), 4786-4800. https://doi.org/10.1111/j.1365-294X.2008.03939.x

Dionne, M., Caron, F., Dodson, J. J., & Bernatchez, L. (2008). Landscape genetics and hierarchical genetic structure in Atlantic salmon: The interaction of gene flow and local adaptation. *Molecular Ecology*, 17(10), 2382-2396. https://doi.org/10.1111/j.1365-294X.2008.03771.x

Druet, T., Pérez-Pardal, L., Charlier, C., & Gautier, M. (2013). Identification of large selective sweeps associated with major genes in cattle. *Animal Genetics*, 44(6), 758-762. https://doi.org/10.1111/age.12073

Elmer, K. R., Fan, S., Kusche, H., Spreitzer, M. L., Kautt, A. F., Franchini, P., & Meyer, A. (2014). Parallel evolution of Nicaraguan crater lake cichlid fishes via non-parallel routes. *Nature Communications*, 5, 5168. https://doi.org/10.1038/ncomms6168

Fang, Z., & Luna, E. J. (2013). Supervillin-mediated suppression of p53 protein enhances cell survival. *Journal of Biological Chemistry*, 288(11), 7918-7929. https://doi.org/10.1074/jbc.M112.416842

Fariello, M. I., Servin, B., Tosser-Klopp, G., Rupp, R., Moreno, C.; International Sheep Genomics Consortium, ... Boitard, S. (2014). Selection signatures in worldwide sheep populations. *PLoS ONE*, 9(8), e103813. https://doi.org/10.1371/journal.pone.0103813

Field, M., Tarpey, P. S., Smith, R., Edkins, S., O' Meara, S., Stevens, C., ... Raymond, F. L. (2007). Mutations in the BRWD3 gene cause X-linked mental retardation associated with macrocephaly. *American Journal of Human Genetics*, 81(2), 367-374. https://doi.org/10.1086/520677

Foll, M., & Gaggiotti, O. (2008). A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: A Bayesian perspective. *Genetics*, 180, 977-993. https://doi.org/10.1534/genetics.108.092221

Frichot, E., Schoville, S. D., Bouchard, G., & François, O. (2013). Testing for associations between loci and environmental gradients using latent factor mixed models. *Molecular Biology and Evolution*, 30(7), 1687-1699. https://doi.org/10.1093/molbev/mst063

Friedrich, J., Brand, B., & Schwerin, M. (2015). Genetics of cattle temperament and its impact on livestock production and breeding—A review. *Archiv fuer Tierzucht*, 58(1), 13.

Fuentes-Pardo, A. P., & Ruzzante, D. E. (2017). Whole-genome sequencing approaches for conservation biology: Advantages, limitations, and practical recommendations. *Molecular Ecology*, 26(20), 5369-5406. https://doi.org/10.1111/mec.14264

Gagnaire, P.-A., & Gaggiotti, O. E. (2016). Detecting polygenic selection in marine populations by combining population genomics and quantitative genetics approaches. *Current Zoology*, 62(6), 603-616. https://doi.org/10.1093/cz/zow088

Gautier, M., & Vitalis, R. (2012). rehh: An R package to detect footprints of selection in genome-wide SNP data from haplotype structure. *Bioinformatics*, 28(8), 1176-1177. https://doi.org/10.1093/bioinformatics/bts115

Gjedrem, T. (2005). *Selection and breeding programs in aquaculture*. Dordrecht, The Netherlands: Springer. https://doi.org/10.1007/1-4020-3342-7

Gjedrem, T. (2012). Genetic improvement for the development of efficient global aquaculture: A personal opinion review. *Aquaculture*, 344-349, 12-22. https://doi.org/10.1016/j.aquaculture.2012.03.003

Gjedrem, T., Robinson, N., & Rye, M. (2012). The importance of selective breeding in aquaculture to meet future demands for animal protein: A review. *Aquaculture*, 350-353, 117-129. https://doi.org/10.1016/j.aquaculture.2012.04.008

Glover, K., Ottera, H., Olsen, R., Slinde, E., Taranger, G., & Skaala, O. (2009). A comparison of farmed, wild and hybrid Atlantic salmon (*Salmo salar*) reared under farming conditions. *Aquaculture*, 286, 203-210. https://doi.org/10.1016/j.aquaculture.2008.09.023

Golej, D. L., Askari, B., Kramer, F., Barnhart, S., Vivekanandan-Giri, A., Pennathur, S., & Bornfeldt, K. E. (2011). Long-chain acyl-CoA synthetase 4 modulates prostaglandin E(2) release from human arterial

smooth muscle cells. *Journal of Lipid Research*, 52(4), 782-793. https://doi.org/10.1194/jlr.M013292

Gutierrez, A. P., Lubieniecki, K. P., Davidson, E. A., Lien, S., Kent, M. P., Fukui, S., ... Davidson, W. S. (2012). Genetic mapping of quantitative trait loci (QTL) for body-weight in Atlantic salmon (*Salmo salar*) using a 6.5K SNP array. *Aquaculture*, 358-359, 61-70. https://doi.org/10.1016/j.aquaculture.2012.06.017

Gutierrez, A. P., Lubieniecki, K. P., Fukui, S., Withler, R. E., Swift, B., & Davidson, W. S. (2014). Detection of quantitative trait loci (QTL) related to grilsing and late sexual maturation in Atlantic salmon (*Salmo salar*). *Marine Biotechnology* (New York, N.Y.), 16(1), 103-110. https://doi.org/10.1007/s10126-013-9530-3

Gutierrez, A. P., Yáñez, J. M., & Davidson, W. S. (2015). Evidence of recent signatures of selection during domestication in an Atlantic salmon population. *Marine Genomics*, 26, 41-50.

Gutierrez, A. P., Yáñez, J. M., Fukui, S., Swift, B., & Davidson, W. S. (2015). Genome-wide association study (GWAS) for growth rate and age at sexual maturation in Atlantic salmon (*Salmo salar*). *PLoS ONE*, 10(3), e0119730. https://doi.org/10.1371/journal.pone.0119730

Heyne, H. O., Lautenschläger, S., Nelson, R., Besnier, F., Rotival, M., Cagan, A., ... Albert, F. W. (2014). Genetic influences on brain gene expression in rats selected for tameness and aggression. *Genetics*, 198(3), 1277. https://doi.org/10.1534/genetics.114.168948

Hohenlohe, P. A., Bassham, S., Etter, P. D., Stiffler, N., Johnson, E. A., & Cresko, W. A. (2010). Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags. *PLoS Genetics*, 6(2), e1000862.

Houston, R. D., Bishop, S. C., Hamilton, A., Guy, D. R., Tinch, A. E., Taggart, J. B., ... Haley, C. S. (2009). Detection of QTL affecting harvest traits in a commercial Atlantic salmon population. *Animal Genetics*, 40, 753-755. https://doi.org/10.1111/j.1365-2052.2009.01883.x

Jensen, J. D., Foll, M., & Bernatchez, L. (2016). The past, present and future of genomic scans for selection. *Molecular Ecology*, 25(1), 1-4. https://doi.org/10.1111/mec.13493

Johansson, I., & Rendel, J. (1968). *Genetics and animal breeding*. San Francisco, CA: W. H. Freeman.

Johnston, I. A., Alderson, R., Sandham, C., Dingwall, A., Mitchell, D., Selkirk, C., ... Springate, J. (2000). Muscle fibre density in relation to the colour and texture of smoked Atlantic salmon (*Salmo salar* L.). *Aquaculture*, 189(3-4), 335-349. https://doi.org/10.1016/S0044-8486(00)00373-2

Johnston, I. A., Alderson, R., Sandham, C., Mitchell, D., Selkirk, C., Dingwall, A., ... Springate, J. (2000). Patterns of muscle growth in early and late maturing populations of Atlantic salmon (*Salmo salar* L.). *Aquaculture*, 189(3-4), 307-333. https://doi.org/10.1016/S0044-8486(00)00372-0

Johnston, I. A., Bower, N. I., & Macqueen, D. J. (2011). Growth and the regulation of myotomal muscle mass in teleost fish. *Journal of Experimental Biology*, 214(10), 1617-1628. https://doi.org/10.1242/jeb.038620

Jombart, T. (2008). adegenet: An R package for the multivariate analysis of genetic markers. *Bioinformatics*, 24(11), 1403-1405. https://doi.org/10.1093/bioinformatics/btn129

Kelley, J. L., Brown, A. P., Therkildsen, N. O., & Foote, A. D. (2016). The life aquatic: Advances in marine vertebrate genomics. *Nature Reviews Genetics*, 17(9), 523-534. https://doi.org/10.1038/nrg.2016.66

Kijas, J. W., Lenstra, J. A., Hayes, B., Boitard, S., Porto Neto, L. R., San Cristobal, M., ... International Sheep Genomics Consortium Members (2012). Genome-wide analysis of the world's sheep breeds reveals high levels of historic mixture and strong recent selection. *PLoS Biology*, 10(2), e1001258. https://doi.org/10.1371/journal.pbio.1001258

King, T. L. (2007). *Biodiversity and population structure. The Atlantic salmon: Genetics, conservation and management*. E. Verspoor, L. Stradmeyer, & J. L. Nielsen (Eds.). Oxford: Blackwell Publishing.
Kittilsen, S., Johansen, I. B., Braastad, B. O., & Øverli, Ø. (2012). Pigments, parasites and personality: Towards a unifying role for steroid hormones? *PLoS ONE*, 7(4), e34281. https://doi.org/10.1371/journal.pone.0034281

Krasnov, A., Skugor, S., Todorcevic, M., Glover, K. A., & Nilsen, F. (2012). Gene expression in Atlantic salmon skin in response to infection with the parasitic copepod *Lepeophtheirus salmonis*, cortisol implant, and their combination. *BMC Genomics*, 13(1), 1-15.

Li, N., & Li, Y. (2014). Ubiquitin-mediated control of seed size in plants. *Frontiers in Plant Science*, 5, 332.

Lien, S., Koop, B. F., Sandve, S. R., Miller, J. R., Kent, M. P., Nome, T., ... Davidson, W. S. (2016). The Atlantic salmon genome provides insights into rediploidization. *Nature*, 533(7602), 200-205. https://doi.org/10.1038/nature17164

Lin, Y., Chen, Y., Yi, C., Fong, J. J., Kim, W., Rius, M., & Zhan, A. (2017). Genetic signatures of natural selection in a model invasive ascidian. *Scientific Reports*, 7, 44080. https://doi.org/10.1038/srep44080

Lin, Y.-N., Roy, A., Yan, W., Burns, K. H., & Matzuk, M. M. (2007). Loss of zona pellucida binding proteins in the acrosomal matrix disrupts acrosome biogenesis and sperm morphogenesis. *Molecular and Cellular Biology*, 27(19), 6794-6805. https://doi.org/10.1128/MCB.01029-07

López, M. E., Neira, R., & Yáñez, J. M. (2015). Applications in the search for genomic selection signatures in fish. *Frontiers in Genetics*, 5, 458.

Lorenzen, K., Beveridge, M. C. M., & Mangel, M. (2012). Cultured fish: Integrative biology and management of domestication and interactions with wild fish. *Biological Reviews*, 87(3), 639-660. https://doi.org/10.1111/j.1469-185X.2011.00215.x

Lotterhos, K. E., & Whitlock, M. C. (2015). The relative power of genome scans to detect local adaptation depends on sampling design and statistical method. *Molecular Ecology*, 24(5), 1031-1046. https://doi.org/10.1111/mec.13100

Ma, Y., Wei, J., Zhang, Q., Chen, L., Wang, J., Liu, J., & Ding, X. (2015). A genome scan for selection signatures in pigs. *PLoS ONE*, 10(3), e0116850. https://doi.org/10.1371/journal.pone.0116850

Macqueen, D. J., Primmer, C. R., Houston, R. D., Nowak, B. F., Bernatchez, L., Bergseth, S., ... Yáñez, J. M. (2017). Functional Annotation of All Salmonid Genomes (FAASG): An international initiative supporting future salmonid research, conservation and aquaculture. *BMC Genomics*, 18(1), 484. https://doi.org/10.1186/s12864-017-3862-8

Mäkinen, H., Vasemägi, A., McGinnity, P., Cross, T. F., & Primmer, C. R. (2014). Population genomic analyses of early-phase Atlantic Salmon (*Salmo salar*) domestication/captive breeding. *Evolutionary Applications*, 8(1), 93-107.

Martinez, V., Dettleff, P., Lopez, P., Fernandez, G., Jedlicki, A., Yáñez, J. M., & Davidson, W. S. (2013). Assessing footprints of selection in commercial Atlantic salmon populations using microsatellite data. *Animal Genetics*, 44(2), 223-226. https://doi.org/10.1111/j.1365-2052.2012.02387.x

Moore, J.-S., Bourret, V., Dionne, M., Bradbury, I., O'Reilly, P., Kent, M., ... Bernatchez, L. (2014). Conservation genomics of anadromous Atlantic salmon across its North American range: Outlier loci identify the same patterns of population structure as neutral loci. *Molecular Ecology*, 23(23), 5680-5697. https://doi.org/10.1111/mec.12972

Nielsen, R., Hellmann, I., Hubisz, M., Bustamante, C., & Clark, A. G. (2007). Recent and ongoing selection in the human genome. *Nature Reviews Genetics*, 8, 857-868. https://doi.org/10.1038/nrg2187

Núñez-Acuña, G., Aguilar-Espinoza, A., Chávez-Mardones, J., & Gallardo-Escárate, C. (2012). Ubiquitin-conjugating enzyme E2-like gene associated to pathogen response in *Concholepas concholepas*: SNP identification and transcription expression. *Fish & Shellfish Immunology*, 33(4), 1065-1068. https://doi.org/10.1016/j.fsi.2012.08.024

Oksenberg, N., & Ahituv, N. (2013). The role of AUTS2 in neurodevelopment and human evolution. *Trends in Genetics*, 29(10), 600-608. https://doi.org/10.1016/j.tig.2013.08.001

Oleksyk, T. K., Smith, M. W., & O'Brien, S. J. (2009). Genome-wide scans for footprints of natural selection. *Philosophical Transactions of the Royal Society of London B: Biological Sciences*, 365(1537), 185-205.

Perrier, C., Bourret, V., Kent, M. P., & Bernatchez, L. (2013). Parallel and nonparallel genome-wide divergence among replicate population pairs of freshwater and anadromous Atlantic salmon. *Molecular Ecology*, 22(22), 5577-5593. https://doi.org/10.1111/mec.12500

Perrier, C., Guyomard, R., Bagliniere, J. L., & Evanno, G. (2011). Determinants of hierarchical genetic structure in Atlantic salmon populations: Environmental factors vs. anthropogenic influences. *Molecular Ecology*, 20(20), 4231-4245. https://doi.org/10.1111/j.1365-294X.2011.05266.x

Petersen, J. L., Mickelson, J. R., Rendahl, A. K., Valberg, S. J., Andersson, L. S., Axelsson, J., ... McCue, M. E. (2013). Genome-wide analysis reveals selection for important traits in domestic horse breeds. *PLOS Genetics*, 9(1), e1003211. https://doi.org/10.1371/journal.pgen.1003211

Porto-Neto, L. R., Lee, S. H., Lee, H. K., & Gondro, C. (2013). Detection of signatures of selection using $F_{ST}$. In *Genome-wide association studies and genomic prediction* (pp. 423-436). Totowa, NJ: Humana Press.

Price, E. O. (1984). Behavioral aspects of animal domestication. *The Quarterly Review of Biology*, 59(1), 1-32. https://doi.org/10.1086/413673

Pujolar, J. M., Ferchaud, A. L., Bekkevold, D., & Hansen, M. M. (2017). Non-parallel divergence across freshwater and marine three-spined stickleback *Gasterosteus aculeatus* populations. *Journal of Fish Biology*, 91(1), 175-194. https://doi.org/10.1111/jfb.13336

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., ... Sham, P. C. (2007). PLINK: A tool set for whole-genome association and population-based linkage analyses. *American Journal of Human Genetics*, 81(3), 559-575. https://doi.org/10.1086/519795

Qanbari, S., Gianola, D., Hayes, B., Schenkel, F., Miller, S., Moore, S., ... Simianer, H. (2011). Application of site and haplotype-frequency based approaches for detecting selection signatures in cattle. *BMC Genomics*, 12(1), 318. https://doi.org/10.1186/1471-2164-12-318

Qanbari, S., & Simianer, H. (2014). Mapping signatures of positive selection in the genome of livestock. *Livestock Science*, 166, 133-143. https://doi.org/10.1016/j.livsci.2014.05.003

Reid, D. P., Szanto, A., Glebe, B., Danzmann, R. G., & Ferguson, M. M. (2005). 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*, 94(2), 166-172. https://doi.org/10.1038/sj.hdy.6800590

Rellstab, C., Gugerli, F., Eckert, A. J., Hancock, A. M., & Holderegger, R. (2015). A practical guide to environmental association analysis in landscape genomics. *Molecular Ecology*, 24(17), 4348-4370. https://doi.org/10.1111/mec.13322

Rhee, K., Brunori, M., Besset, V., Trousdale, R., & Wolgemuth, D. J. (1998). Expression and potential role of Fsg1, a murine bromodomain-containing homologue of the *Drosophila* gene female sterile homeotic. *Journal of Cell Science*, 111(23), 3541-3550.

Roberge, C., Einum, S., Guderley, H., & Bernatchez, L. (2006). Rapid parallel evolutionary changes of gene transcription profiles in farmed Atlantic salmon. *Molecular Ecology*, 15(1), 9-20.

Rotllant, G., Wade, N. M., Arnold, S. J., Coman, G. J., Preston, N. P., & Glencross, B. D. (2015). Identification of genes involved in reproduction and lipid pathway metabolism in wild and domesticated shrimps. *Marine Genomics*, 22, 55-61. https://doi.org/10.1016/j.margen.2015.04.001

Rougemont, Q., & Bernatchez, L. (2018). Reconstructing the demographic history of Atlantic salmon (*Salmo salar*) across its distribution range using approximate Bayesian computations. *BioRxiv*, 142372.

Rye, M., Gjerde, B., & Gjedrem, T. (2010). Genetic improvement programs for aquaculture species in developed countries. Paper read at 9th World
congress on genetics applied to livestock production, Liepzing, Germany.

Sandbakken, M., Ebbesson, L., Stefansson, S., & Helvik, J. V. (2012). Isolation and characterization of melanopsin photoreceptors of Atlantic salmon (*Salmo salar*). *Journal of Comparative Neurology*, **520**(16), 3727-3744. https://doi.org/10.1002/cne.23125

Sauvage, C., Derôme, N., Normandeau, E., St-Cyr, J., Audet, C., & Bernatchez, L. (2010). Fast transcriptional responses to domestication in the brook charr *Salvelinus fontinalis*. *Genetics*, **185**(1), 105-112. https://doi.org/10.1534/genetics.110.115071

Spidle, A. P., Kalinowski, S. T., Lubinski, B. A., Perkins, D. L., Beland, K. F., Kocik, J. F., & King, T. L. (2003). Population structure of Atlantic salmon in maine with reference to populations from Atlantic Canada. *Transactions of the American Fisheries Society*, **132**(2), 196-209. https://doi.org/10.1577/1548-8659(2003)132&lt;0196:PSOASI&gt;2.0.CO;2

Tang, K., Thornton, K. R., & Stoneking, M. (2007). A new approach for using genome scans to detect recent positive selection in the human genome. *PLoS Biology*, **5**(7), e171. https://doi.org/10.1371/journal.pbio.0050171

Teletchea, F., & Fontaine, P. (2014). Levels of domestication in fish: Implications for the sustainable future of aquaculture. *Fish and Fisheries*, **15**(2), 181-195. https://doi.org/10.1111/faf.12006

Tsai, H. Y., Hamilton, A., Guy, D. R., Tinch, A. E., Bishop, S. C., & Houston, R. D. (2015). The genetic architecture of growth and fillet traits in farmed Atlantic salmon (*Salmo salar*). *BMC Genetics*, **16**(1), 1-11.

Utsunomiya, Y. T., Pérez, A. M., O'Brien, T. S., Sonstegard, J. S., & Garcia, J. F. (2015). Genomic data as the “hitchhiker’s guide” to cattle adaptation: Tracking the milestones of past selection in the bovine genome. *Frontiers in Genetics*, **6**, 36.

Vasemägi, A., Nilsson, J., McGinnity, P., Cross, T., O'Reilly, P., Glebe, B., ... Primmer, C. R. (2012). Screen for footprints of selection during domestication/captive breeding of Atlantic salmon. *Comparative and Functional Genomics*, **2012**(1), 1-14. https://doi.org/10.1155/2012/628204

Verspoor, E., Beardmore, J. A., Consuegra, S., Garcia de Leaniz, C., Hindar, K., Jordan, W. C., ... Sanchez, J. A. (2005). Population structure in the Atlantic salmon: Insights from 40 years of research into genetic protein variation. *Journal of Fish Biology*, **67**, 3-54. https://doi.org/10.1111/j.0022-1112.2005.00838.x

Vriens, J., Owsianik, G., Hofmann, T., Philipp, S. E., Stab, J., Chen, X., ... Voets, T. (2011). TRPM3 is a nociceptor channel involved in the detection of noxious heat. *Neuron*, **70**(3), 482-494. https://doi.org/10.1016/j.neuron.2011.02.051

Wang, G.-D., Xie, H.-B., Peng, M.-S., Irwin, D., & Zhang, Y.-P. (2014). Domestication genomics: Evidence from animals. *Annual Review of Animal Biosciences*, **2**(1), 65-84. https://doi.org/10.1146/annurev-animal-022513-114129

Weir, B. S., & Cockerham, C. C. (1984). Estimating *F*-statistics for the analysis of population structure. *Evolution*, **38**(6), 1358-1370.

Withler, R. E., Supernault, K. J., & Miller, K. M. (2005). Genetic variation within and among domesticated Atlantic salmon broodstocks in British Columbia, Canada. *Animal Genetics*, **36**(1), 43-50. https://doi.org/10.1111/j.1365-2052.2004.01220.x

Worzfeld, T., Swiercz, J. M., Sentürk, A., Genz, B., Korostylev, A., Deng, S., ... Offermanns, S. (2014). Genetic dissection of plexin signaling in vivo. *Proceedings of the National Academy of Sciences of the United States of America*, **111**(6), 2194-2199. https://doi.org/10.1073/pnas.1308418111

Wringe, B., Devlin, R., Ferguson, M., Moghadam, H., Sakhrani, D., & Danzmann, R. (2010). Growth-related quantitative trait loci in domestic and wild rainbow trout (*Oncorhynchus mykiss*). *BMC Genetics*, **11**(1), 63. https://doi.org/10.1186/1471-2156-11-63

Xu, P., McIntyre, L. M., Scardina, J., Wheeler, P. A., Thorgaard, G. H., & Nichols, K. M. (2011). Transcriptome profiling of embryonic development rate in rainbow trout advanced backcross introgression lines. *Marine Biotechnology*, **13**(2), 215-231. https://doi.org/10.1007/s10126-010-9283-1

Yáñez, J. M., Naswa, S., López, M. E., Bassini, L., Correa, K., Gilbey, J., ... Maass, A. (2016). Genomewide single nucleotide polymorphism discovery in Atlantic salmon (*Salmo salar*): Validation in wild and farmed American and European populations. *Molecular Ecology Resources*, **16**(4), 1002-1011.

Yang, S., Li, X., Li, K., Fan, B., & Tang, Z. (2014). A genome-wide scan for signatures of selection in Chinese indigenous and commercial pig breeds. *BMC Genetics*, **15**, 7. https://doi.org/10.1186/1471-2156-15-7

Zhao, F., McParland, S., Kearney, F., Lixin, D., & Berry, D. P. (2015). Detection of selection signatures in dairy and beef cattle using high-density genomic information. *Genetics, Selection, Evolution: GSE*, **47**(1), 49. https://doi.org/10.1186/s12711-015-0127-3

Zhu, J., Shi, X., Lu, H., Xia, B., Li, Y., Li, X., ... Yang, G. (2016). RNA-seq transcriptome analysis of extensor digitorum longus and soleus muscles in large white pigs. *Molecular Genetics and Genomics*, **291**(2), 687-701. https://doi.org/10.1007/s00438-015-1138-z

Zueva, K. J., Lumme, J., Veselov, A. E., Kent, M. P., Lien, S., & Primmer, C. R. (2014). Footprints of directional selection in wild Atlantic salmon populations: Evidence for parasite-driven evolution? *PLoS ONE*, **9**(3), e91672. https://doi.org/10.1371/journal.pone.0091672

## SUPPORTING INFORMATION

Additional supporting information may be found online in the Supporting Information section at the end of the article.

**How to cite this article:** López ME, Benestan L, Moore J-S, et al. Comparing genomic signatures of domestication in two Atlantic salmon (*Salmo salar* L.) populations with different geographical origins. *Evol Appl.* 2018;00:1-20. https://doi.org/10.1111/eva.12689