- Research article
- Open Access
Human-induced evolution caught in action: SNP-array reveals rapid amphi-atlantic spread of pesticide resistance in the salmon ecotoparasite Lepeophtheirus salmonis
© Besnier et al.; licensee BioMed Central Ltd. 2014
- Received: 10 June 2014
- Accepted: 16 October 2014
- Published: 26 October 2014
The salmon louse, Lepeophtheirus salmonis, is an ectoparasite of salmonids that causes huge economic losses in salmon farming, and has also been causatively linked with declines of wild salmonid populations. Lice control on farms is reliant upon a few groups of pesticides that have all shown time-limited efficiency due to resistance development. However, to date, this example of human-induced evolution is poorly documented at the population level due to the lack of molecular tools. As such, important evolutionary and management questions, linked to the development and dispersal of pesticide resistance in this parasite, remain unanswered. Here, we introduce the first Single Nucleotide Polymorphism (SNP) array for the salmon louse, which includes 6000 markers, and present a population genomic scan using this array on 576 lice from twelve farms distributed across the North Atlantic.
Our results support the hypothesis of a single panmictic population of lice in the Atlantic, and importantly, revealed very strong selective sweeps on linkage groups 1 and 5. These sweeps included candidate genes potentially connected to pesticide resistance. After genotyping a further 576 lice from 12 full sibling families, a genome-wide association analysis established a highly significant association between the major sweep on linkage group 5 and resistance to emamectin benzoate, the most widely used pesticide in salmonid aquaculture for more than a decade.
The analysis of conserved haplotypes across samples from the Atlantic strongly suggests that emamectin benzoate resistance developed at a single source, and rapidly spread across the Atlantic within the period 1999 when the chemical was first introduced, to 2010 when samples for the present study were obtained. These results provide unique insights into the development and spread of pesticide resistance in the marine environment, and identify a small genomic region strongly linked to emamectin benzoate resistance. Finally, these results have highly significant implications for the way pesticide resistance is considered and managed within the aquaculture industry.
- Linkage Group
- Selective Sweep
- Single Nucleotide Polymorphism Array
- Pesticide Resistance
- Louse Infestation
The salmon louse Lepeophtheirus salmonis (Krøyer, 1838) is an ectoparasitic copepod that feeds on the mucus, skin and blood of salmonid fishes in the marine environment. Through its feeding action, L. salmonis leads to stress [1, 2] and in severe untreated cases, open wounds and death of the host fish . Within the Atlantic salmon (Salmo salar L.) and rainbow trout (Onchorhynchus mykiss) aquaculture industry that is primarily conducted in marine cages, L. salmonis infections cause highly significant economic losses [4, 5]. In addition, due to the rapid expansion of aquaculture, the evolutionary dynamics of the relationship between L. salmonis and its hosts has changed . Furthermore, there is evidence that L. salmonis is associated with declines of wild salmonid populations [5, 7, 8], mortality in the marine environment [9, 10], and modifications of life history traits such as the age of return from the sea . As a consequence, lice control regimes have been put into place to reduce the level of sea lice infestation in marine farms.
A range of integrated methods are currently employed or being developed to control infestations of sea lice on salmonids in marine farms [12–16]. However, the industry continues to rely heavily upon the use of a few certified pesticides to remove lice from infested fish. This is not considered as a sustainable approach [17, 18] due to a high risk of resistance development. Currently, over 500 arthropod pest-species are documented as having developed resistance to pesticides , and L. salmonis, as for arthropods with short generation times and high reproduction rates, has a strong capacity to adapt to new environmental constraints. Indeed, L. salmonis is known for developing resistance to chemical treatments [19–22], and loss of efficiency has been reported for a number of chemicals, for example organophosphates  and pyrethroids , and no new chemicals have been introduced since 1999 in Norway when SLICE® came on the market. During the past ten years, the main pesticide used for treating lice infestations in the North-Atlantic has been SLICE® emamectin benzoate (EMB) . Emamectin benzoate is an avermectin where the target is gamma-aminobutyrate (GABA) and glutamate-gated chloride channels (GABA-Cls and Glu-Cls) . However, the target site of this chemical in L. salmonis is not confirmed, and over reliance on the product has also led to the development of reduced sensitivity and resistance to EMB in sea lice populations [27–31].
Despite thorough efforts to unravel the mechanisms behind the development of pesticide resistance in L. salmonis[20, 21, 32, 33], little is known about the genetics and genomics of this species at the population level. Whether resistance to pesticides has developed in L. salmonis in multiple farms and regions in parallel, or has developed in one location and spread thereafter, is completely unknown. Development of resistance in a pest population submitted to intense control is however expected to significantly alter the pattern of allelic frequencies around the positively selected genomic region [34, 35]. In the case of the salmon louse, addressing this issue is not only essential for the development of management guidelines, it also has a high level of relevance to questions regarding the evolutionary relationships among lice throughout the Atlantic. This is because the question of population division remains a subject that has not been resolved, with previous studies being hampered by either a limited number of samples [36–38], or limited numbers of polymorphic loci [39, 40].
A 6 k SNP array was produced for L. salmonis, based on an initial set of 640,000 polymorphic sites detected in five pooled samples of L. salmonis originating from five regions of the North Atlantic. Among the 6000 SNPs selected for inclusion on the array, 5080 were selected to cover the genome with an average interval of 100 kb (“genome distributed”), 174 were selected to cover a few larger scaffolds with average inter-SNP interval of 10 kb (“LD”), and 190 were selected based upon being located in or close to putative genes, based on matches to Expressed Sequence Tag (“ESTs”) regions. In addition, 556 SNPs were selected to provide a potential diagnostic tool to differentiate between populations throughout the Atlantic (“diagnostic”).
After production, 5540 of the original 6000 SNPs (92.3%) were successfully converted into on-chip assays including 508 diagnostic (91.2%), 170 EST (89.5%), 159 LD (91.4%), and 4703 (92.6%) genome distributed.
The SNP array was then used to genotype DNA collected from 576 individuals sampled from 12 commercial farms situated in the North Atlantic (Figure 1). After automatic clustering and manual validation, the average sample call rate was 99.5%, as was the average SNP call rate. Of the 5540 from 6000 SNPs, 121 (2.2%) were monomorphic, 209 (3.8%) displayed evidence of secondary SNPs within the flanking sequence making genotype calls unreliable, and 110 (2.0%) assays failed. Data from these three SNP categories, totaling 440 assays, was excluded from the present study. The remaining 5100 SNPs were further categorized as showing typical biallelic distribution with good cluster separation (SNP; 4666, 84.2%), showing a tight but distinct cluster distribution between allelotypes and requiring manual validation (n = 68, 1.2%), and 366 (6.6%) presenting some evidence of atypical clustering and also demanding manual checking. In a parallel project to the present study, a genetic linkage map was constructed (unpublished data) which included 5091 of the successfully genotyped SNPs. Data from these 5091 SNPs are included in the present study (Additional file 1).
In addition to F statistics, the genotypes of all 576 individuals on the neutral markers were analyzed with Bayesian clustering implemented in STRUCTURE , and with the R package ADEgenet  using the “find.clusters” function. With those two approaches, the parameter K, which stands for the number of populations or clusters in the dataset, was tested for values from K = 1 to 12. Both log likelihood of the data calculated by STRUCTURE and the Bayesian Information Criterion (BIC) calculated from ADEgenet reached a minimum value (-351000 and 3790 respectively) for K = 1 as the number of clusters in the data. The two approaches gave consistent results, and indicated that L. salmonis consists of one single population throughout the North Atlantic.
Correlation (r 2 ) between five genome scans for outlier markers
Arlequin FCT 12pop
Estimated effective population size from each sampling site
Lower 95% CI
Upper 95% CI
Genomic regions affected by directional selection
Detection of selective sweeps from the two summary statistics, CLR and Omega, on the twelve sampling sites
Chr1 (92–94 cM)
Chr5 (63–67 cM)
P < 0.05
P < 0.05
P < 0.05
P < 0.05
P < 0.05
P < 0.05
P < 0.05
P < 0.05
P < 0.05
P < 0.05
P < 0.05
P < 0.05
P < 0.05
P < 0.05
P < 0.05
Matches for coding DNA were found in the genomic regions under positive selection, in 3 and 15 scaffolds on LG1 and LG5 respectively. Those scaffolds were 440 Kb length on average, and covered together 1.7 Mb on LG1 and 6.2 Mb on LG5. The number of predicted genes was respectively 13 and 110 in the two regions (Figure 6E and F, Additional file 2: Table S1). Functional annotation from BLAST hits revealed that some of these genes could play a role in drug resistance mechanisms. This was the case for Q96J66 an ATP-binding cassette (ABC) transporter protein predicted near marker c8501.30679 on LG1, Q28FI8, a cytochrome b5 domain-containing protein predicted near marker c8125.67718 on LG 5, or P17970 and P08510, two Potassium voltage-gated channel proteins predicted on LG 5 near markers c10059.23896, and c15581.45069 respectively (Additional file 2: Table S1, and Additional file 3: Table S2).
Despite showing a weaker signals for selective sweep than the reported regions on LG1 and LG5, one genomic region on LG14 did gather an important number of SNPs with high FST values (Figure 2B). Moreover, a portion of LG14 between 28 and 31 cM showed significant signal for selective sweeps in the Faeroe Island and Shetland samples (Additional file 4: Figure S1). This region was annotated following the same protocol as for the two previous regions. Matches for coding DNA were found for 6 scaffolds within which 19 genes were predicted. Among those, 1 gene (Q86GC8) coded for acetycholinesterase (Additional file 2: Table S1). This is the target for organophosphates and could be involved in resistance as a knock down mutation.
Association between selective sweep regions and drug resistance
Extrapolation of selection time frame
In LG5, we detected conserved haplotypes on a 5 cM region, which corresponds to a genomic region where chromosome recombination would occur with 5% frequency every generation i.e. once every 20 meiosis. Assuming the most drastic selection scenario where the frequency of the advantageous variant is nearly 100% in the population after a punctual selection event (exposure to chemical treatment of lice infestation), the frequency of a 5 cM haplotype associated to this mutation would decrease by 5% every generation following selection. A lapse of 15 generations would thus be needed to observe a decrease of haplotype frequency to 46% (as in sample S856) and, similarly, 43 generations to decrease to 10%.
The population genomic structure of L. salmonis throughout the Atlantic ocean was investigated using F statistics , principal component , and Bayesian clustering analysis  on genetic data from >5000 SNPs and 12 sets of samples from six regions. These analyses did not reveal any geographical genetic structure of L. salmonis throughout the North Atlantic. The overall global FST value was very low, with 98% of the markers displaying FST lower than 0.1. Furthermore, no significant correlation between genetic distance, measured as FST, and geographical distance, were detected. Both Bayesian clustering and principal component analyses suggested K = 1 as the most likely number of clusters in the samples. Previous studies of population genetic structure of L. salmonis throughout the Atlantic have revealed contrasting results, ranging from highly significant genetic differentiation among nearby farms, to no or little structure throughout the entire North Atlantic. However, many of these earlier studies have been hindered by technical challenges, low sample sizes, low numbers of polymorphic markers, or combinations of these challenges [36–40]. Nevertheless, the most rigorous of the earlier studies, using data from four mitochondrial DNA (mtDNA) genes , and microsatellite markers on large numbers of samples [39, 40] have all revealed statistically weak or non-significant genetic differentiation throughout the Atlantic. The results of these former studies together with the present population genomic analysis strongly suggest that L. salmonis is characterized by a single panmictic population throughout the North Atlantic. The most likely mechanism for the absence of isolation by distance in such a large geographical area is transport by the salmon host. Wild Atlantic salmon migrating from rivers located on both sides of the Atlantic are known to mix on distant feeding grounds in the high seas of the Atlantic ocean . Seaward migrating smolts originating from these American and European rivers may often be infected by salmon lice in their respective coastal zones, and thus transport these to the offshore feeding grounds where cross-infection of salmon from other geographic regions can occur. Such a mechanism would explain the lack of isolation by distance reported here. This hypothesis is further supported by the short time frame needed for an advantageous mutation to spread in the whole North Atlantic as reported here.
In several sites, the estimate (Nê) of the effective population size was infinite. This is usually interpreted as a signal of a large population where genetic drift is unlikely to be of major influence on the evolutionary properties of the population. The estimate of infinity occurs because the method implemented in LDNe , as other methods for estimating Ne, depends on a signal that is a function of 1/Ne . As a consequence, the precision in estimating Ne is good for small populations, but the method has difficulties in distinguishing between large and infinite populations.
Among the 12 sampling sites studied here, Nê varied between 333 and infinity. This estimation remains imprecise, but is sufficient to suggest that a high proportion of the population inheriting large haplotype segments from a common ancestor, as observed in this study, is unlikely to happen by chance.
Genomic regions affected by directional selection
Selective sweeps are the signature of a strong selection event where the incremented frequency of a selected locus in a population is accompanied by hitchhiking of the flanking genomic regions [34, 35]. The size of a genomic region affected by hitchhiking depends upon the recombination rate and the strength of the selection. In cases of high recombination rate and weak selection, the area may be small and thus, difficult to detect with a moderate density of molecular markers. In the literature, occurrence of selective sweeps have been reported in extensively studied organisms for which molecular resources are highly available, such as human [49, 50], or Drosophila[51–53], and in domestic species for which directional selection was very strong in the recent history [54–57]. Another example of reported selective sweeps is the development of chemical resistance in pest species or human disease vectors that are subject to intense population control [58–62]. In this latter example, the selection pressure is often very strong and led the selected variant to colonize the population in a short period of time. Such examples of human-induced evolution appear to be easier to detect as they happen on a contemporary time scale where recombination has not had the chance to erode the haplotype segments in the hitchhiking regions. In such cases, selective sweeps may be detectable with a medium density of molecular markers. This last case scenario has remarkable similarities with the results presented here, which re-enforces the main hypothesis of selection driven by drug resistance.
Combining two independent methods increases the power of sweep detection . Here, we focused on genomic regions where selective sweeps were reported using both the CLR and Omega methods. Using these criteria, two selective sweeps were identified on linkage groups 1 and 5 respectively. These two genomic regions were already identified to contain a large number of the outlier markers identified as under positive selection, and displayed local disturbance in LD. It is thus highly likely that these two regions contain one or several advantageous mutations that were positively selected in L. salmonis.
The two selective sweep regions contained large haplotypes with frequencies higher that 10% within several sampling sites across the entire North Atlantic. Because the selection of an advantageous mutation can be accompanied by genetic hitchhiking , which creates LD patterns around the target of positive selection [64, 65], it is not surprising that haplotypes under the regions where selective sweep were detected are present in high frequencies in each sample. This is an expected consequence of the stronger LD observed in those regions. Interestingly, identical haplotypes were also frequent in geographically distant regions. The presence of selective sweeps coinciding with large haplotype segments conserved across distant localities strongly suggests that a de-novo mutation appeared between 10 to 40 generations before the sample used in the present study were collected. The causative mutation is likely to have occurred in one single geographical site, before being spread to the other regions rather than originating from multiple and independent selection events. This consolidates the conclusion above that L. salmonis is represented by a single panmictic population throughout the North Atlantic, and importantly, demonstrates that alleles conveying resistance to pesticides may be quickly spread over very large areas in the marine environment. This unique documentation of genetic resistance dispersal in the ocean reveals a major challenge for the management of resistance development within the aquaculture industry in the entire Atlantic. Because resistance can quickly spread across the ocean, management of pesticide resistance in this parasite needs to be addressed on an ocean-wide level rather than on a regional level. Moreover, due to the heterogeneous structure of the landscape where L. salmonis is evolving, with patches of high host density in salmon farms and coastal areas, and large areas of low host density in the offshore regions, standard population genetic models are likely to produce biased estimates of the evolutionary dynamics of this organism .
In order to better identify the possible targets of positive selection detected on linkage groups 1 and 5, the two genomic regions were annotated to identify sequences of coding DNA together with the possible function of such sequences. The number of predicted genes on the scaffolds situated in the genomic regions under selection was 13 and 110 on linkage groups 1 and 5 respectively. Functional annotation was obtained from BLAST hits . In most cases, BLAST from a given sequence returned several times the same gene in several organisms. In such cases, we reported the version of the gene that gave the best hit (Figure 5E and F). One of those genes, Q96J66, reported on LG1, codes for an ATP binding cassette protein that has previously been described as multi-drug resistance associated protein . In addition, the P-glycoprotein (Pgp), a member of the ATP-binding cassette (ABC) transporter protein superfamily was recently reported for playing an important role in EMB resistance [33, 69]. However, other studies could not establish significant correlation between Pgp polymorphisms and EMB susceptibility [70, 71]. A second candidate gene for drug resistance is Q28FI8 on linkage group 5, coding for Cytochrome b5 domain-containing protein. Cytochrome b5 is a known activator of cytochrome P450 , which has documented insecticide resistance effects in several species [73, 74]. This is concordant with the fact that EMB is one of the most extensively used chemical pesticides to control lice infestation in sea cages , and that intensive use of the substance led to development of resistance [27–30].
Because exposure to pesticides used for treating lice infestation is probably the strongest directional selection pressure applied to L. salmonis at the temporal scale of the present study, it is likely that the selective sweeps reported here are connected to positive selection on variants that provide better resistance to pesticides. The observations reported in the present paper would fit with this hypothesis: the region of the sweep on LG1 contains one candidate gene for a family of protein that has already been identified for playing a role in resistance to the most commonly used delousing chemical in salmon farms. Moreover, the use of the same chemical in all regions, as it is the case for in-feed lice infestation treatment, is expected to induce a strong selection pressure on the same resistance locus in all regions. This second fact also coincides with our observations of haplotypes under selective sweeps that are highly conserved across distant regions of the North Atlantic. Together, those facts point at pesticide resistance as the most likely cause for the observed selective sweeps.
One genomic region on LG14 displayed a number of SNPs with high FST values (Figure 2B), and was reported with significant signal for selective sweeps in Faeroe Island and Shetland samples (Additional file 4: Figure S1). Within the 19 genes predicted in this region, one coded for acetycholinesterase, which is the well-known target for organophosphates. Because it has been established that a variant type of acetylcholinesterase provides resistance to organophosphates in lice , it is possible that the high values of FST among sampling sites and the weak signal the of selective sweep reported on LG14 are due to an older event of selection for lice resistant to organophosphates. Organophosphates were the only chemical used to control sea lice up to around 1995 and major resistance problems occurred and left the chemical largely ineffective. However, due to emerging resistance against both EMB and phyretroids, Azametiphos (an organophosphate) was reintroduced and have been used alone and in combination with other medicines during the last 6–7 years.
Link between EMB resistance and selective sweeps
The possible link between the detected selective sweeps and pesticide resistance was specifically investigated for EMB, which is the most commonly used pesticide to control L. salmonis in the Atlantic in the period 2000–2010. A genome-wide association analysis revealed that a broad region around the selective sweep of LG5 had a strong and significant contribution to the variation in EMB resistance (Figure 7). This last analysis demonstrates a clear link between genetic variation on LG5 around the selective sweep region, and EMB resistance. The selective sweep region on LG1 did however not show any link to EMB resistance. This could indicate that despite coinciding with a coding sequence for a protein that has previously been described as multi-drug resistance associated , the selective sweep on LG1 is caused by another factor than resistance to EMB. Alternatively, we point out that the lice that were included in the genome wide association dataset were all sampled in South Norway, whereas in the population genetic dataset only the samples from Ireland, Shetland and Faeroe Island displayed highly frequent haplotypes under the selective weep region on LG1. It is thus possible that the variant under positive selection on LG1 were not present in the association dataset.
We also point out that the power to detect association between genomic regions and EMB resistance was better in the F2 dataset than from the F1 cross between sensitive and resistant forms. In the F2 dataset, 5 genomic regions were significantly associated with the resistance while only two in the F1 data. When the same region was significantly associated in both datasets, the association score was larger in the F2 data.
An F2 intercross represents a notoriously better experimental design to maximize the genotypic variability in a population, while F1 crosses are not suitable for populations with fixed allelic frequencies. Assuming that the resistant and sensitive lines had fixed allelic frequencies at loci associated with EMB resistance, the resulting F1 hybrids would all be heterozygous at those loci, thus reducing to zero the genotypic variability of that population together with the power to associate genotypes to phenotypes.
We produced the first SNP array for L. salmonis, and used it to genotype a set of geographically distinct samples in order to assess whether it was possible to identify signatures of pesticide resistance and dispersal, as well as to evaluate population genomic structure of this economically and ecologically important parasite. All results supported the hypothesis that this parasite is characterized by a single panmictic population through its distribution in the Atlantic. Importantly, the search for signatures of positive selection revealed two genomic regions under strong positive selection. Both of the selective sweeps identified coincided with the coding sequence of one or several possible candidate genes for pesticide resistance. Additionally, one genomic region on LG4 displayed a weak signal for selection, which may correspond to an older selective sweep caused by the development of resistance to organophosphates which have been used for a longer period. Genome-wide association analysis revealed a strong link between the selective sweep on linkage group 5 and resistance to EMB which was the most commonly used pesticide between years 2000 to 2010, when the samples upon which the study is based were collected. The strong selective sweeps identified on linkage groups 1 and 5 were both characterized by haplotypes that were present in high frequency in samples from distant regions of the North Atlantic. This last observation strongly suggests that L. salmonis has a high capacity to spread new advantageous mutations across continents in the time scale of just a few generations or years (at most 11 years), thus corroborating concerns that pesticide resistance can develop and rapidly spread over large areas on an ecological time-scale. These results have very significant implications for the way in which pesticide resistance for L salmonis, and potentially other parasitic copepods, is managed in marine aquaculture.
Generating sequences from L. salmonis
DNA sequences were obtained from L. salmonis samples collected in 2009 from five different regions of the north Atlantic. Four of these sets of samples have been described previously ; C858 (Canada), S856 (Shetland), I852 (Ireland), N849 (Northern Norway) (Figure 1). The fifth sample was collected in September 2008 from an emamectin-benzoate-desensitized population in Austevoll, Western Norway. Eight adult female L. salmonis were sampled from each of the five regions . For all samples, DNA was isolated in a 96-well format using the DNeasy kit according to the manufacturer’s instructions (Qiagen, Hilden, Germany). Equal amounts of DNA from each of the eight individuals from each station were pooled to meet concentration demands and were sequenced by Fasteris SA using the Illumina HiSeq 2000 platform following their standard protocols . In addition to the five regional samples, DNA was also sequenced from Expressed Sequence Tags (EST). In total, 715.106 reads were generated from the five regions and EST samples, with an average length of 100 bp/read.
SNP detection and chip design
The largest fraction of the SNPs were selected to provide a regular coverage over the 700 Mb genome assembly. To optimize the genome coverage, SNPs were chosen with intervals of 100 Kb. SNPs were thus chosen from the largest scaffolds available, with 100 Kb bases interval on scaffolds larger that 100 Kb, or only one SNP per scaffold, for scaffolds smaller than 100 Kb. SNPs also had to match the following quality criterion: having a minimum of two reads depth per region in at least four regions out of five, and to be polymorphic (displaying both reference and alternative nucleotide) in at least three regions.
A second set of 186 SNPs was selected in genomic regions aligned with EST sequences. These SNPs had a minimum of two reads depth in at least three regions out of five.
A third set of SNPs was selected to help discriminate individuals from the five geographical regions, based on genotype information. In this group, SNPs were selected for having coverage deeper than two reads in at least three regions. There are ten pair-wise combinations of the five geographical regions. For each pair, a set of SNPs was chosen for being monomorphic for one allele in the first region, monomorphic for the alternative allele in the second region, and polymorphic for the remaining region(s), providing thus a possible tool for determining the origin of individuals, based on genotype information. Fifty of such diagnostic SNPs were selected for each regional pair, with the exception of two pairs: Ireland-Shetland and Ireland-Canada where only 35 diagnostic SNPs could be found for each pair. Based on the same principle of selection, fifty SNPs were selected for discriminating West Atlantic (Canada) from the East Atlantic (Ireland Shetland Austevoll Norland) regions, and fifty SNPs were also selected for discriminating Austevoll (showing resistance to chemical treatment) from the other regions.
Finally, a set of 174 SNPs was selected at narrow intervals within the same scaffolds in order to provide information about linkage disequilibrium (LD) in salmon lice populations. These SNPs had a minimum of two reads depth in at least two regions out of five, and were polymorphic in at east one region. These 174 “LD” markers were distributed on 24 scaffolds among 11 linkage groups, and separated by an average of 10Kb when situated on the same scaffold.
A total of 576 samples were collected from twelve sampling sites in six regions in the North Atlantic (Figure 1). Five of the sampling sites were the same as those used to produce the DNA libraries previous to SNP detection: C858 (Canada), S856 (Shetland), I852 (Ireland), N849 (Northern Norway), and N813 (South Norway) . In addition, genotypes were also obtained from seven sites also described in previous reports ; N854 (South Norway), N837 (North Norway), F850 and F851 (Faeroe Islands), S855 (Shetland), I853 (Ireland), C857 (Canada), (Figure 1). The samples are thus divided in six geographical regions and two sampling sites per region. A total of 48 individuals were genotyped in each sampling site.
All 576 samples were quantified using picogreen fluorescent stain (Invitrogen, USA) and an aliquot from each examined on 1% agarose gel to subjectively assess DNA quality based on the presence of high molecular weight DNA. Where necessary, samples were concentrated using a Speed-vac to fall within the target concentration range of 25-75 ng/ul required for genotyping. Genotyping was performed according to the infinium HD assay protocol (Illumina, San Diego). SNP genotyping results were quality checked to eliminate unreliable markers. Clusters were inspected using the Genotyping Module within Genome Studio. Genotype data from 5090 polymorphic SNPs were obtained for the 576 samples.
Markers were tested for Hardy-Weinberg equilibrium in each population. Markers with allele frequencies that deviated significantly (p <0.01 with Bonferroni correction) from HW equilibrium in at least one of the populations were considered outliers. Scans for loci under selection were conducted with three different software: i) The FST outlier based approach Fdist  implemented in LOSITAN , ii) The Bayesian approach implemented in BayeScan , and iii) the hierarchical F statistic FCT approach  implemented in Arlequin 3.5 . Both approaches use the FST measurement of genetic distance as basis for detecting directional selection. In i) and iii), directional selection is revealed by the distribution of the ratio of the genetic distance by heterozygicity (FST/He), whereas in ii), the genetic distance among populations is decomposed in two components: a population-specific component (Beta) shared by all loci and a locus-specific component (Alfa). In the later, directional selection is revealed by a significant value of Alfa. In iii) the data is stratified by population in a hierarchical manner, from local demes to the largest population. The three approaches handle the detection of markers under selection in a different manner and are expected to produce slightly different results.
Two alternative datasets were generated: one dataset (noted 6P) where the two sampling sites within each region are assumed to come from the same population, and one dataset (noted 12P) where each sampling site is considered as a different population. In total, five genome scans were performed: Both 6P and 12P datasets were scanned with LOSITAN and BayeScan, while the FCT approach was only used on the 12P dataset because it already takes into account the hierarchical structure of the data and uncertainty regarding the number of demes in the population structure [81, 82]. After detection of outlier markers, the data were split into two sets of markers: markers under selection if reported as outlier by three or more of the above-mentioned genome scans, and neutral markers if reported as neutral by both genome scans. This criterion produced one set of 4810 neutral markers and 167 markers under selection.
Genomic regions affected by directional selection
In addition to scan for individual outlier marker, the genome was scanned for regions of reduced genetic variation linked to a recently fixed beneficial mutation, or selective sweeps. At the molecular level, selective sweeps are characterized by a reduction of the genetic variability in the flanking sequences from both sides of the selected locus  and a local increment of linkage disequilibrium, also from both sides of the selected locus , Based on these two features, two summary statistics were used to assess the presence of selective sweeps. The first approach is a composite likelihood ratio (CLR) based on site frequency spectrum [44, 83], and implemented in SweeD software . This method consists in scanning the genome for regions of reduced genetic diversity. The second approach is the Omega statistic  implemented in the OmegaPlus software , which scans the genome for local patterns of high linkage disequilibrium. Both methods were applied to scan the 5090 SNPs genotypes from 576 individuals. For genomic regions reporting significant evidence of selective sweeps, haplotypes were reconstructed independently within each sampling site with Phase2.1 software [86, 87].
The significance thresholds of the scans for selective sweeps were obtained by bootstrapping after simulating 1000 data sets under neutral model. Simulations were performed with MSMS software .
600,000 polymorphic sites were simulated on 576 individuals.
The 576 individuals were divided in 12 groups of 48 individual each, as in the sampling sites of the present data.
Eight individuals were randomly picked from five groups in order to replicate the pooled DNA sequences that were used for the SNP selection process.
From the sites that matched the SNP selection criteria (being polymorphic in at least three groups out of five), 6000 were randomly picked to generate a neutral dataset of 6000 SNP genotypes on 576 individuals.
SNPs determined to be under selection were mapped to preliminary genome scaffolds and compared to a draft gene annotation (http://sealouse.imr.no/). The identified scaffolds had an average length of 440 Kb, and possible genes present on these scaffolds were predicted by i- ab initio gene prediction [89, 90] on the scaffold sequence, and ii- functional annotation of genes based on BLAST hits . The predicted peptide sequences were compared to SwissProt using BLASTP with E-value threshold of 10-6.
Genome-wide association mapping for drug resistance
In order to investigate the possible link between selective sweeps and pesticide resistance, a second set of lice were also genotyped on the SNP array. This second dataset consisted of 576 individuals divided in 12 full sibling families that originated from three L. salmonis source strains with different tolerances of EMB . Strains had been maintained in the salmon louse rearing facility using protocols previously described . Five families were produced as hybrids between two laboratory strains, expected to be ‘tolerant’ and ‘susceptible’ to EMB , while the remaining seven full sibling families were produced as the F2 generation of the third source strain, that had uncertain level of tolerance to EMB. The resistant strain was collected in September 2008 near site N813 (Figure 1) after reports of Slice® treatment failures, while the susceptible strain was collected in October 2009 in Oslofjorden, east of Norway where no commercial salmon farming is taking place and exposure to EMB is low. The third strain was collected in April 2010 near site N813 (Figure 1), in a salmon farm where no reduced efficiency of Slice® had been reported despite being in a dense farming area.
The offspring of each family were exposed to EMB according to the following protocol: all individuals were placed in 5 L glass beakers and exposed to a one-dose (50 ppb) EMB overnight trial, performed by the method for bioassay testing described in Handbook in resistance management . After 20 hours of exposure, each individual was evaluated in accordance with the response criteria, defined as ‘Living’ or ‘Moribund and/or dead’, and stored on 95% ethanol for subsequent DNA analysis . All individuals were genotyped on the SNP array, and the software FastPhase  was used to reconstruct haplotypes in all linkage groups. SNPs were then pooled in groups of three flanking markers haplotypes, in sliding windows along the genome. At each genomic position, a hierarchical generalized mixed model with binomial family was fitted with the R package hglm . The model included the status after exposure (dead or alive) as binary response to one fixed full sib family effect and one random haplotype effect. The dispersion coefficient of the random haplotype effect was kept as indicator of the genetic contribution to the variance in survival after EMB exposure.
The genome scan was repeated 1000 times with a randomized response vector to obtain an empirical distribution of the random genetic effect under null hypothesis. The upper 95 percentile of this empirical distribution was kept as 5% genome-wide significance threshold.
The genetic distance between samples was measured for each SNP with an R implementation of Weir and Cockerham formula [46, 94]. In order to test for possible isolation by distance among regions, pair wise genetic distance and geographical distance were compared with standard Mantel test  implemented in the R package vegan [94, 96]. Statistical significance of Mantel test was corrected a posteriori for multiple testing by applying Bonferroni correction. To investigate for possible population structure on the North Atlantic, the genotype data of the 4810 neutral SNPs was analyzed in the population genetic software STRUCTURE 2.3  without population information, using 50,000 burnin, 200,000 iterations. The number of cluster (K) was tested for K = 1 to K = 12, with ten replicates for each value of K. Due to the long computation it requires, this analysis was performed on an 8-core workstation using the R package ParallelStructure  to distribute parallel runs of STRUCTURE on multiple cores. The number of cluster in the populations, K was then assessed by Evanno’s method  using the STRUCTURE HARVESTER implementation . In addition, the genotype data was also analyzed with the “find.clusters” function from the R package ADEgenet [42, 94], also in order to infer the number of populations in the data.
Linkage disequilibrium was tested for each pair of markers by calculating the coefficient of correlation between genotype pairs r2. In addition, the effective size of each population (Ne) was estimated from population linkage disequilibrium with the LDNe software . This program was implemented to estimate Ne from unlinked markers; therefore, we use a subset of 100 unlinked SNPs from the total set of 5090 (Additional file 5).
The data sets supporting the results of this article are included within the article and its additional files.
This study was funded by the Norwegian Research Council through the PrevenT platform project. We would also like to acknowledge Ole Torrissen and fish farmers for providing samples of lice from their farms, as well as two anonymous reviewers for their constructive comments that helped to improve the manuscript.
- Skugor S, Glover KA, Nilsen F, Krasnov A: Local and systemic gene expression responses of Atlantic salmon (Salmo salar L.) to infection with the salmon louse (Lepeophtheirus salmonis). BMC Genomics. 2008, 9: 498-10.1186/1471-2164-9-498.PubMed CentralPubMedView ArticleGoogle Scholar
- Krasnov A, Skugor S, Todorcevic M, Glover K, Nilsen F: Gene expression in Atlantic salmon skin in response to infection with the parasitic copepod Lepeophtheirus salmonis, cortisol implant, and their combination. BMC Genomics. 2012, 13: 130-10.1186/1471-2164-13-130.PubMed CentralPubMedView ArticleGoogle Scholar
- Bjørn PA, Finstad B: The development of salmon lice (Lepeophtheirus salmonis) on artificially infected post smolts of sea trout (Salmo trutta). Can J Zool. 1998, 76: 970-977. 10.1139/cjz-76-5-970.View ArticleGoogle Scholar
- Costello MJ: The global economic cost of sea lice to the salmonid farming industry. J Fish Dis. 2009, 32: 115-118. 10.1111/j.1365-2761.2008.01011.x.PubMedView ArticleGoogle Scholar
- Torrissen O, Jones S, Asche F, Guttormsen A, Skilbrei OT, Nilsen F, Horsberg TE, Jackson D: Salmon lice –impact on wild salmonids and salmon aquaculture. J Fish Dis. 2013, 36: 171-194. 10.1111/jfd.12061.PubMed CentralPubMedView ArticleGoogle Scholar
- Mennerat A, Hamre L, Ebert D, Nilsen F, Davidova M, Skorping A: Life history and virulence are linked in the ectoparasitic salmon louse Lepeophtheirus salmonis. J Evol Biol. 2012, 25: 856-861. 10.1111/j.1420-9101.2012.02474.x.PubMedView ArticleGoogle Scholar
- Krkosek M, Ford JS, Morton A, Lele S, Myers RA, Lewis MA: Declining wild salmon populations in relation to parasites from farm salmon. Science. 2007, 318: 1772-1775. 10.1126/science.1148744.PubMedView ArticleGoogle Scholar
- Krkosek M, Connors BM, Morton A, Lewis MA, Dill LM, Hilborn R: Effects of parasites from salmon farms on productivity of wild salmon. Proc Natl Acad Sci U S A. 2011, 108: 14700-14704. 10.1073/pnas.1101845108.PubMed CentralPubMedView ArticleGoogle Scholar
- Krkosek M, Revie CW, Gargan PG, Skilbrei OT, Finstad B, Todd CD: Impact of parasites on salmon recruitment in the Northeast Atlantic Ocean. Proc Biol Sci. 2013, 280: 20122359-PubMed CentralPubMedView ArticleGoogle Scholar
- Costello MJ: How sea lice from salmon farms may cause wild salmonid declines in Europe and North America and be a threat to fishes elsewhere. Proc R Soc B Biol Sci. 2009, 276: 3385-3394. 10.1098/rspb.2009.0771.View ArticleGoogle Scholar
- Vollset KW, Barlaup BT, Skoglund H, Normann ES, Skilbrei OT: Salmon lice increase the age of returning Atlantic salmon. Biol Lett. 2014, 10: 20130896-10.1098/rsbl.2013.0896.PubMed CentralPubMedView ArticleGoogle Scholar
- Glover KA, Aasmundstad T, Nilsen F, Storset A, Skaala Ø: Variation of Atlantic salmon families (Salmo salar L.) in susceptibility to the sea lice Lepeophtheirus salmonis and Caligus elongatus. Aquaculture. 2005, 245: 19-30. 10.1016/j.aquaculture.2004.11.047.View ArticleGoogle Scholar
- Gharbi K, Glover K, Stone L, MacDonald E, Matthews L, Grimholt U, Stear M: Genetic dissection of MHC-associated susceptibility to Lepeophtheirus salmonis in Atlantic salmon. BMC Genet. 2009, 10: 20-PubMed CentralPubMedView ArticleGoogle Scholar
- Dempster T, Kristiansen TS, Korsoen OJ, Fosseidengen JE, Oppedal F: Technical note: Modifying Atlantic salmon (Salmo salar) jumping behavior to facilitate innovation of parasitic sea lice control techniques. J Anim Sci. 2011, 89: 4281-4285. 10.2527/jas.2011-3894.PubMedView ArticleGoogle Scholar
- Bui S, Oppedal F, Korsøen ØJ, Dempster T: Modifying Atlantic salmon behaviour with light or feed stimuli may improve parasite control techniques. Aquac Environ Interact. 2013, 3: 125-133. 10.3354/aei00055.View ArticleGoogle Scholar
- Skiftesvik AB, Bjelland RM, Durif CMF, Johansen IS, Browman HI: Delousing of Atlantic salmon (Salmo salar) by cultured vs. wild ballan wrasse (Labrus bergylta). Aquaculture. 2013, 402–403: 113-118.View ArticleGoogle Scholar
- Boxaspen K: A review of the biology and genetics of sea lice. ICES J Mar Sci J Cons. 2006, 63: 1304-1316. 10.1016/j.icesjms.2006.04.017.View ArticleGoogle Scholar
- Igboeli OO, Burka JF, Fast MD: Lepeophtheirus salmonis: a persisting challenge for salmon aquaculture. Anim Front. 2014, 4: 22-32. 10.2527/af.2014-0004.View ArticleGoogle Scholar
- Denholm I, Devine GJ, Horsberg TE, Sevatdal S, Fallang A, Nolan DV, Powell R: Analysis and management of resistance to chemotherapeutants in salmon lice, Lepeophtheirus salmonis (Copepoda: Caligidae). Pest Manag Sci. 2002, 58: 528-536. 10.1002/ps.482.PubMedView ArticleGoogle Scholar
- Fallang A, Ramsay JM, Sevatdal S, Burka JF, Jewess P, Hammell KL, Horsberg TE: Evidence for occurrence of an organophosphate-resistant type of acetylcholinesterase in strains of sea lice (Lepeophtheirus salmonis Kroyer). Pest Manag Sci. 2004, 60: 1163-1170. 10.1002/ps.932.PubMedView ArticleGoogle Scholar
- Fallang A, Denholm I, Horsberg TE, Williamson MS: Novel point mutation in the sodium channel gene of pyrethroid-resistant sea lice Lepeophtheirus salmonis (Crustacea: Copepoda). Aquat Organ. 2005, 65: 129-136.View ArticleGoogle Scholar
- Sevatdal S, Copley L, Wallace C, Jackson D, Horsberg TE: Monitoring of the sensitivity of sea lice (Lepeophtheirus salmonis) to pyrethroids in Norway, Ireland and Scotland using bioassays and probit modelling. Aquaculture. 2005, 244: 19-27. 10.1016/j.aquaculture.2004.11.009.View ArticleGoogle Scholar
- Jones MW, Sommerville C, Wootten R: Reduced sensitivity of the salmon louse, Lepeophtheirus salmonis, to the organophosphate dichlorvos. J Fish Dis. 1992, 15: 197-202. 10.1111/j.1365-2761.1992.tb00654.x.View ArticleGoogle Scholar
- Sevatdal S, Horsberg TE: Determination of reduced sensitivity in sea lice (Lepeophtheirus salmonis Krøyer) against the pyrethroid deltamethrin using bioassays and probit modelling. Aquaculture. 2003, 218: 21-31. 10.1016/S0044-8486(02)00339-3.View ArticleGoogle Scholar
- Stone J, Sutherland IH, Sommerville CS, Richards RH, Varma KJ: The efficacy of emamectin benzoate as an oral treatment of sea lice, Lepeophtheirus salmonis (Krøyer), infestations in Atlantic salmon, Salmo salar L. J Fish Dis. 1999, 22: 261-270.View ArticleGoogle Scholar
- Ozoe Y: gamma-Aminobutyrate- and Glutamate-gated Chloride Channels as Targets of Insecticides. Target Recept Control Insect Pests Pt I, Volume 44. Edited by: Cohen E. 2013, 24–28 Oval Road, London Nw1 7dx, England: Academic Press Ltd-Elsevier Science Ltd, 211-286. Advances in Insect PhysiologyGoogle Scholar
- Westcott JD, Stryhn H, Burka JF, Hammell KL: Optimization and field use of a bioassay to monitor sea lice Lepeophtheirus salmonis sensitivity to emamectin benzoate. Aquat Organ. 2008, 79: 119-131.View ArticleGoogle Scholar
- Westcott J, Revie C, Giffin B, Hamm K: Evidence of Sea Lice Lepeophtheirus Salmonis Tolerance to Emamectin Benzoate in New Brunswick, Canada. 2010Google Scholar
- Jones PG, Hammell KL, Gettinby G, Revie CW: Detection of emamectin benzoate tolerance emergence in different life stages of sea lice, Lepeophtheirus salmonis, on farmed Atlantic salmon, Salmo salar L. J Fish Dis. 2013, 36: 209-220. 10.1111/jfd.12022.PubMedView ArticleGoogle Scholar
- Espedal PG, Glover KA, Horsberg TE, Nilsen F: Emamectin benzoate resistance and fitness in laboratory reared salmon lice (Lepeophtheirus salmonis). Aquaculture. 2013, 416–417: 111-118.View ArticleGoogle Scholar
- Ljungfeldt LE, Espedal PG, Nilsen F, Skern-Mauritzen M, Glover KA: A common-garden experiment to quantify evolutionary processes in copepods: the case of emamectin benzoate resistance in the parasitic sea louse Lepeophtheirus salmonis. BMC Evol Biol. 2014, 14: 108-10.1186/1471-2148-14-108.PubMed CentralPubMedView ArticleGoogle Scholar
- Carmichael SN, Bron JE, Taggart JB, Ireland JH, Bekaert M, Burgess ST, Skuce PJ, Nisbet AJ, Gharbi K, Sturm A: Salmon lice (Lepeophtheirus salmonis) showing varying emamectin benzoate susceptibilities differ in neuronal acetylcholine receptor and GABA-gated chloride channel mRNA expression. BMC Genomics. 2013, 14: 408-10.1186/1471-2164-14-408.PubMed CentralPubMedView ArticleGoogle Scholar
- Igboeli OO, Fast MD, Heumann J, Burka JF: Role of P-glycoprotein in emamectin benzoate (SLICE®) resistance in sea lice, Lepeophtheirus salmonis. Aquaculture. 2012, 344–349: 40-47.View ArticleGoogle Scholar
- Smith JM, Haigh J: The hitch-hiking effect of a favourable gene. Genet Res. 1974, 23: 23-35. 10.1017/S0016672300014634.PubMedView ArticleGoogle Scholar
- Kaplan NL, Hudson RR, Langley CH: The “hitchhiking effect” revisited. Genetics. 1989, 123: 887-899.PubMed CentralPubMedGoogle Scholar
- Nolan DV, Martin SAM, Kelly Y, Glennon K, Palmer R, Smith T, McCormack GP, Powell R: Development of microsatellite PCR typing methodology for the sea louse, Lepeophtheirus salmonis (Krøyer). Aquac Res. 2000, 31: 815-822. 10.1046/j.1365-2109.2000.00514.x.View ArticleGoogle Scholar
- Tjensvoll K, Glover KA, Nylund A: Sequence variation in four mitochondrial genes of the salmon louse Lepeophtheirus salmonis. Aquat Organ. 2006, 68: 251-259.View ArticleGoogle Scholar
- Nolan D, Powell R: Geographic and temporal genetic structure in Lepeophtheirus salmonis from four salmon farms along the northwest and west coasts of Ireland: results from a microsatellite analysis. Hydrobiologia. 2009, 617: 55-63. 10.1007/s10750-008-9525-7.View ArticleGoogle Scholar
- Glover KA, Stølen Åsta B, Messmer A, Koop BF, Torrissen O, Nilsen F: Population genetic structure of the parasitic copepod Lepeophtheirus salmonis throughout the Atlantic. Mar Ecol Prog Ser. 2011, 427: 161-172.View ArticleGoogle Scholar
- Todd CD, Walker AM, Ritchie MG, Graves JA, Walker AF: Population genetic differentiation of sea lice (Lepeophtheirus salmonis) parasitic on Atlantic and Pacific salmonids: analyses of microsatellite DNA variation among wild and farmed hosts. Can J Fish Aquat Sci. 2004, 61: 1176-1190. 10.1139/f04-069.View ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMed CentralPubMedGoogle Scholar
- Jombart T: Adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008, 24: 1403-1405. 10.1093/bioinformatics/btn129.PubMedView ArticleGoogle Scholar
- Waples RS, Do C: Linkage disequilibrium estimates of contemporary Ne using highly variable genetic markers: a largely untapped resource for applied conservation and evolution. Evol Appl. 2010, 3: 244-262. 10.1111/j.1752-4571.2009.00104.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim Y, Stephan W: Detecting a local signature of genetic hitchhiking along a recombining chromosome. Genetics. 2002, 160: 765-777.PubMed CentralPubMedGoogle Scholar
- Ronnegard L, Shen X, Alam M: hglm: A Package for Fitting Hierarchical Generalized Linear Models. R J. 2010, 2: 20-28.Google Scholar
- Weir BS, Cockerham CC: Estimating F-Statistics for the Analysis of Population Structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Jacobsen JA, Hansen LP, Bakkestuen V, Halvorsen R, Reddin DG, White J, O Maoiléidigh N, Russell IC, Potter ECE (T), Fowler M, Smith GW, Mork KA, Isaksson A, Oskarsson S, Karlsson L, Pedersen S: Distribution by origin and sea age of Atlantic salmon (Salmo salar) in the sea around the Faroe Islands based on analysis of historical tag recoveries. ICES J Mar Sci J Cons. 2012, 69: 1598-1608. 10.1093/icesjms/fss115.View ArticleGoogle Scholar
- Waples RS, DO C: ldne: a program for estimating effective population size from data on linkage disequilibrium. Mol Ecol Resour. 2008, 8: 753-756. 10.1111/j.1755-0998.2007.02061.x.PubMedView ArticleGoogle Scholar
- Enard W, Przeworski M, Fisher SE, Lai CSL, Wiebe V, Kitano T, Monaco AP, Paabo S: Molecular evolution of FOXP2, a gene involved in speech and language. Nature. 2002, 418: 869-872. 10.1038/nature01025.PubMedView ArticleGoogle Scholar
- Pickrell JK, Coop G, Novembre J, Kudaravalli S, Li JZ, Absher D, Srinivasan BS, Barsh GS, Myers RM, Feldman MW, Pritchard JK: Signals of recent positive selection in a worldwide sample of human populations. Genome Res. 2009, 19: 826-837. 10.1101/gr.087577.108.PubMed CentralPubMedView ArticleGoogle Scholar
- Glinka S, De Lorenzo D, Stephan W: Evidence of gene conversion associated with a selective sweep in Drosophila melanogaster. Mol Biol Evol. 2006, 23: 1869-1878. 10.1093/molbev/msl069.PubMedView ArticleGoogle Scholar
- Jensen JD, Thornton KR, Andolfatto P: An Approximate Bayesian Estimator Suggests Strong. Recurrent Selective Sweeps in Drosophila. PLoS Genet. 2008, 4: e1000198-10.1371/journal.pgen.1000198.PubMed CentralPubMedView ArticleGoogle Scholar
- Singh ND, Jensen JD, Clark AG, Aquadro CF: Inferences of Demography and Selection in an African Population of Drosophila melanogaster. Genetics. 2013, 193: 215-228. 10.1534/genetics.112.145318.PubMed CentralPubMedView ArticleGoogle Scholar
- Ramey H, Decker J, McKay S, Rolf M, Schnabel R, Taylor J: Detection of selective sweeps in cattle using genome-wide SNP data. BMC Genomics. 2013, 14: 382-10.1186/1471-2164-14-382.PubMed CentralPubMedView ArticleGoogle Scholar
- Petersen JL, Mickelson JR, Rendahl AK, Valberg SJ, Andersson LS, Axelsson J, Bailey E, Bannasch D, Binns MM, Borges AS, Brama P, da Camara MA, Capomaccio S, Cappelli K, Cothran EG, Distl O, Fox-Clipsham L, Graves KT, Guerin G, Haase B, Hasegawa T, Hemmann K, Hill EW, Leeb T, Lindgren G, Lohi H, Lopes MS, McGivney BA, Mikko S, Orr N, et al: Genome-wide analysis reveals selection for important traits in domestic horse breeds. PLoS Genet. 2013, 9: e1003211-10.1371/journal.pgen.1003211.PubMed CentralPubMedView ArticleGoogle Scholar
- Quilez J, Short A, Martinez V, Kennedy L, Ollier W, Sanchez A, Altet L, Francino O: A selective sweep of >8 Mb on chromosome 26 in the Boxer genome. BMC Genomics. 2011, 12: 339-10.1186/1471-2164-12-339.PubMed CentralPubMedView ArticleGoogle Scholar
- Johansson AM, Pettersson ME, Siegel PB, Carlborg O: Genome-wide effects of long-term divergent selection. PLoS Genet. 2010, 6: e1001188-10.1371/journal.pgen.1001188.PubMed CentralPubMedView ArticleGoogle Scholar
- McCollum AM, Basco LK, Tahar R, Udhayakumar V, Escalante AA: Hitchhiking and selective sweeps of Plasmodium falciparum sulfadoxine and pyrimethamine resistance alleles in a population from central Africa. Antimicrob Agents Chemother. 2008, 52: 4089-4097. 10.1128/AAC.00623-08.PubMed CentralPubMedView ArticleGoogle Scholar
- Powles SB, Yu Q: Evolution in Action: Plants Resistant to Herbicides. Annu Rev Plant Biol. 2010, 61: 317-347. 10.1146/annurev-arplant-042809-112119.PubMedView ArticleGoogle Scholar
- Cheeseman IH, Miller BA, Nair S, Nkhoma S, Tan A, Tan JC, Al Saai S, Phyo AP, Moo CL, Lwin KM, McGready R, Ashley E, Imwong M, Stepniewska K, Yi P, Dondorp AM, Mayxay M, Newton PN, White NJ, Nosten F, Ferdig MT, Anderson TJC: A Major Genome Region Underlying Artemisinin Resistance in Malaria. Science. 2012, 336: 79-82. 10.1126/science.1215966.PubMed CentralPubMedView ArticleGoogle Scholar
- Jones CM, Liyanapathirana M, Agossa FR, Weetman D, Ranson H, Donnelly MJ, Wilding CS: Footprints of positive selection associated with a mutation (N1575Y) in the voltage-gated sodium channel of Anopheles gambiae. Proc Natl Acad Sci. 2012, 109: 6614-6619. 10.1073/pnas.1201475109.PubMed CentralPubMedView ArticleGoogle Scholar
- Lynd A, Weetman D, Barbosa S, Egyir Yawson A, Mitchell S, Pinto J, Hastings I, Donnelly MJ: Field, Genetic, and Modeling Approaches Show Strong Positive Selection Acting upon an Insecticide Resistance Mutation in Anopheles gambiae s.s. Mol Biol Evol. 2010, 27: 1117-1125. 10.1093/molbev/msq002.PubMed CentralPubMedView ArticleGoogle Scholar
- Pavlidis P, Jensen JD, Stephan W: Searching for footprints of positive selection in whole-genome SNP data from nonequilibrium populations. Genetics. 2010, 185: 907-922. 10.1534/genetics.110.116459.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim Y, Nielsen R: Linkage disequilibrium as a signature of selective sweeps. Genetics. 2004, 167: 1513-1524. 10.1534/genetics.103.025387.PubMed CentralPubMedView ArticleGoogle Scholar
- Stephan W, Song YS, Langley CH: The hitchhiking effect on linkage disequilibrium between linked neutral loci. Genetics. 2006, 172: 2647-2663.PubMed CentralPubMedView ArticleGoogle Scholar
- Schneider KA, Kim Y: Genetic Hitchhiking under Heterogeneous Spatial Selection Pressures. PLoS One. 2013, 8: e61742-10.1371/journal.pone.0061742.PubMed CentralPubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1016/S0022-2836(05)80360-2.PubMedView ArticleGoogle Scholar
- Xu M, Molento M, Blackhall W, Ribeiro P, Beech R, Prichard R: Ivermectin resistance in nematodes may be caused by alteration of P-glycoprotein homolog. Mol Biochem Parasitol. 1998, 91: 327-335. 10.1016/S0166-6851(97)00215-6.PubMedView ArticleGoogle Scholar
- Tribble ND, Burka JF, Kibenge FSB: Evidence for changes in the transcription levels of two putative P-glycoprotein genes in sea lice (Lepeophtheirus salmonis) in response to emamectin benzoate exposure. Mol Biochem Parasitol. 2007, 153: 59-65. 10.1016/j.molbiopara.2007.02.002.PubMedView ArticleGoogle Scholar
- Heumann J, Carmichael S, Bron JE, Tildesley A, Sturm A: Molecular cloning and characterisation of a novel P-glycoprotein in the salmon louse Lepeophtheirus salmonis. Comp Biochem Physiol C Toxicol Pharmacol. 2012, 155: 198-205. 10.1016/j.cbpc.2011.08.004.PubMedView ArticleGoogle Scholar
- Heumann J, Carmichael SN, Bron JE, Sturm A: Isolation and characterisation of four partial cDNA sequences encoding multidrug resistance-associated proteins (MRPs) in the salmon louse Lepeophtheirus salmonis (Krøyer, 1837). Aquaculture. 2014, 424–425: 207-214.View ArticleGoogle Scholar
- Schenkman JB, Jansson I: The many roles of cytochrome b5. Pharmacol Ther. 2003, 97: 139-152. 10.1016/S0163-7258(02)00327-3.PubMedView ArticleGoogle Scholar
- Scott JG: Cytochromes P450 and insecticide resistance. Insect Biochem Mol Biol. 1999, 29: 757-777. 10.1016/S0965-1748(99)00038-7.PubMedView ArticleGoogle Scholar
- Scott JG, Wen Z: Cytochromes P450 of insects: the tip of the iceberg. Pest Manag Sci. 2001, 57: 958-967. 10.1002/ps.354.PubMedView ArticleGoogle Scholar
- Skern-Mauritzen R, Malde K, Besnier F, Nilsen F, Jonassen I, Reinhardt R, Koop B, Dalvin S, Mæhle S, Kongshaug H, Glover KA: How does sequence variability affect de novo assembly quality?. J Nat Hist. 2013, 47: 901-910. 10.1080/00222933.2012.738833.View ArticleGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralPubMedView ArticleGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079. 10.1093/bioinformatics/btp352.PubMed CentralPubMedView ArticleGoogle Scholar
- Beaumont MA, Nichols RA: Evaluating Loci for Use in the Genetic Analysis of Population Structure. Proc Biol Sci. 1996, 263: 1619-1626. 10.1098/rspb.1996.0237.View ArticleGoogle Scholar
- Antao T, Lopes A, Lopes RJ, Beja-Pereira A, Luikart G: LOSITAN: a workbench to detect molecular adaptation based on a Fst-outlier method. BMC Bioinformatics. 2008, 9: 323-10.1186/1471-2105-9-323.PubMed CentralPubMedView ArticleGoogle Scholar
- Foll M, Gaggiotti O: A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics. 2008, 180: 977-993. 10.1534/genetics.108.092221.PubMed CentralPubMedView ArticleGoogle Scholar
- Excoffier L, Hofer T, Foll M: Detecting loci under selection in a hierarchically structured population. Heredity. 2009, 103: 285-298. 10.1038/hdy.2009.74.PubMedView ArticleGoogle Scholar
- Excoffier L, Lischer HEL: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010, 10: 564-567. 10.1111/j.1755-0998.2010.02847.x.PubMedView ArticleGoogle Scholar
- Nielsen R, Williamson S, Kim Y, Hubisz MJ, Clark AG, Bustamante C: Genomic scans for selective sweeps using SNP data. Genome Res. 2005, 15: 1566-1575. 10.1101/gr.4252305.PubMed CentralPubMedView ArticleGoogle Scholar
- Pavlidis P, Zivkovic D, Stamatakis A, Alachiotis N: SweeD: likelihood-based detection of selective sweeps in thousands of genomes. Mol Biol Evol. 2013, 30: 2224-2234. 10.1093/molbev/mst112.PubMed CentralPubMedView ArticleGoogle Scholar
- Alachiotis N, Stamatakis A, Pavlidis P: OmegaPlus: a scalable tool for rapid detection of selective sweeps in whole-genome datasets. Bioinformatics. 2012, 28: 2274-2275. 10.1093/bioinformatics/bts419.PubMedView ArticleGoogle Scholar
- Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68: 978-989. 10.1086/319501.PubMed CentralPubMedView ArticleGoogle Scholar
- Stephens M, Donnelly P: A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003, 73: 1162-1169. 10.1086/379378.PubMed CentralPubMedView ArticleGoogle Scholar
- Ewing G, Hermisson J: MSMS: a coalescent simulation program including recombination, demographic structure and selection at a single locus. Bioinformatics. 2010, 26: 2064-2065. 10.1093/bioinformatics/btq322.PubMed CentralPubMedView ArticleGoogle Scholar
- Stanke M, Waack S: Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics. 2003, 19 (Suppl 2): ii215-ii225.PubMedView ArticleGoogle Scholar
- Stanke M, Diekhans M, Baertsch R, Haussler D: Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008, 24: 637-644. 10.1093/bioinformatics/btn013.PubMedView ArticleGoogle Scholar
- Hamre LA, Glover KA, Nilsen F: Establishment and characterisation of salmon louse (Lepeophtheirus salmonis (Krøyer 1837)) laboratory strains. Parasitol Int. 2009, 58: 451-460. 10.1016/j.parint.2009.08.009.PubMedView ArticleGoogle Scholar
- Consortium S: Sea Lice Resistance to Chemotherapeutants: A Handbook in Resistance Management. 2006Google Scholar
- Scheet P, Stephens M: A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet. 2006, 78: 629-644. 10.1086/502802.PubMed CentralPubMedView ArticleGoogle Scholar
- Team RC: R: A Language and Environment for Statistical Computing. 2013, Vienna: AustriaGoogle Scholar
- Mantel N: The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27: 209-220.PubMedGoogle Scholar
- Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H: Vegan: Community Ecology Package. 2013Google Scholar
- Besnier F, Glover KA: ParallelStructure: a R package to distribute parallel runs of the population genetics program STRUCTURE on multi-core computers. PLoS One. 2013, 8: e70651-10.1371/journal.pone.0070651.PubMed CentralPubMedView ArticleGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.PubMedView ArticleGoogle Scholar
- Earl D, von Holdt B: STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012, 4: 359-361. 10.1007/s12686-011-9548-7.View ArticleGoogle Scholar
- Hill WG, Robertson A: Linkage disequilibrium in finite populations. Theor Appl Genet. 1968, 38: 226-231. 10.1007/BF01245622.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.