Field of the invention

The present invention relates to a method for detection of a mutation, in particular a single polymorphism associated with organophosphate resistant Lepeophtheirus salmonis, and oligonucleotide sequences and kits useful in the method of the present invention. The present invention furthermore provides kits and reagents useful for the detection of organophosphate resistance associated mutations.

Background of the invention

Sea lice {Lepeophtheirus salmonis and Caligus spp.) are the major pathogens affecting global salmon farming industry and have a significant impact in many areas. The annual loss has recently been estimated to€300 million (Costello M. J. (2009)) and the aquaculture industry relays heavily on a few chemotherapeutants for lice control. Emerging resistance development to these drugs increase the necessity to develop new treatment methods (biological, prophylactic and drugs) and tools to avoid increased loss due to sea lice and to ensure a sustainable salmon farming industry in the future. Control measures have relied upon a limited number of chemotherapeutants since the 1970s. Parasite resistance and reduced efficacy have now been reported for the majority of these compounds (Sevatdal S., et al., (2005)).

Although various chemotherapeutants have been in use in the aquaculture industry for more than 30 years, it is only during the last 15 years that such use has been part of some kind of integrated pest management system. Management practices include coordinated salmon production within a defined area, use of single year class of fish, limited production period, fallowing, coordinated restocking, use of wrasse, synchronized treatments during the winter and targeting female lice to reduce the impact of settlement during the spring (Pike A., Wadsworth S. L., (2000)).

Organophosphates belong to a group of anti-parasitic agents that acts by inhibiting a vital enzyme, acetyl choline esterase, in cholinergic synapses. The inhibition blocks the cleavage of the transmitter substance acetyl choline and results in elevated levels of it in the synaptic cleft thereby causing exitation, paralysis and death of the organism.

Organophosphates have been used against salmon lice (Lepeophtheirus salmonis), a species within a group of ectoparasitic copepods on fish called sea lice, in Norwegian salmonid aquaculture since the late 1970s. The first agent used was metrifonate (Neguvon™), followed by dichlorvos in 1986 (Nuvan™) and azamethiphos (Salmosan™) in 1994 (Torrissen et al. 2013). In 1991 , the first cases of reduced efficacy of organophosphate treatments were noted in Mid-Norway (Denholm et al. 2002). When the use of

azamethiphos was terminated during 1999, the problem of reduced sensitivity in salmon lice against azamethiphos was wide-spread. At the time, the cause for resistance in salmon lice against azamethiphos was not determined. Azamethiphos was re-introduced as a treatment against salmon lice in 2008. In 2009, the first new reports of reduced efficacy of treatments with azamethiphos came from the field. A bioassay to test the sensitivity of the parasites was established by Veso (Sevatdal, unpublished). Since then several treatment failures have been verified as resistance using this bioassay (Sevatdal, unpublished).

Further controls are required to progress towards a true integrated pest management (IPM) system common to other forms of food production. One key to succeeding with an IPM is to develop tools for the management of resistance to the medicines in use (Brook K. (2009). Considerations in developing an integrated pest management program for control of sea lice on farmed salmon in Pacific Canada. Journal of Fish Diseases. 32. 59-73). So far, drug resistance in sea lice has been detected by various types of bioassays as a significant increase in

EC50/LC50. Although the bioassays can detect any type of resistance to a given drug, such methods are not very accurate or sensitive.

The genomes of all organisms undergo spontaneous mutations during their continuing evolution, forming variant forms of progenitor genetic sequences. A mutation may results in an evolutionary advantage or disadvantage relative to a progenitor form or may be neutral. A variant that result in an evolutionary advantage may eventually be incorporated in many members of the species and may thus effectively become the progenitor form.

Furthermore, often various variant forms survive and coexist in a species population. The coexistence of multiple forms of a genetic sequence gives rise to genetic polymorphism, including single-nucleotide polymorphisms (SNPs). A single-nucleotide polymorphism (SNP) is a DNA sequence variation occurring when a single nucleotide— A, T, C or G— in the genome (or other shared sequence) differs between members of a biological species or paired chromosomes in an organism. For example, two DNA fragments from different individuals, AAGCCTA to AAGCTTA, contain a difference in a single nucleotide, commonly referred to as two alleles. Almost all common SNPs have only two alleles. The genomic distribution of SNPs is not

homogenous; SNPs usually occur in non-coding regions more frequently than in coding regions or, in general, where natural selection is acting and fixating the allele of the SNP that constitutes the most favorable genetic adaptation. Up to date, no mutations or SNPs consequently associated with organophosphate resistance in sea lice have been reported.

Efficient and sensitive methods for diagnosing resistance are crucial in order to manage and control drug resistance. Early detection of reduced sensitivity to a chemical can enable effective countermeasures to be enforced at a time point when these have a greater probability of being effective. Therefore, accurate and speedy identification of

organophosphate resistant Lepeophtheirus salmonis is crucial. Detection of organoposphate resistance prior to treatment, and the use of such analyses after treatment to evaluate treatment efficacy constitutes an important determinant for the integrated pest management (IPM) in the aquaculture industry. Summary of invention

The present invention is based on the surprising finding that resistance towards organophosphates commonly used to combat sea lice infestation is linked to an identified mutation present in one of two novel acetylcholine esterase protein genes identified in sea lice (Lepeophtheirus salmonis) resulting in an amino acid substitution compared with the corresponding gene present in organophosphate sensistive sea lice. More particularly, the present invention is based on the identification of novel single-nucleotide polymorphisms (SNPs) in the gene encoding the acetylcholine esterase protein (LS-acel) of the sea lice {Lepeophtheirus salmonis) shown to be involved in the resistance towards

organophosphate-based chemotherapy. Specifically, the present inventors have identified one organophosphate resistance-associated SNPs located in SL-ace l resulting in the substitution of phenyl with tyrosine in position 362 (Phe362Tyr). According to one embodiment an in vitro method is provided for detection of organophosphate resistance in one or more crustaceans comprising the steps of detecting a mutation associated with organophosphate resistance in the

acetylcholine esterase gene 1 of the crustaceans to be analyzed. The crustacean may be one or more copepods, e.g. belonging to to the family Caligidae. According to one embodiment, the copepod is selected from the group consisting of

Lepeophtheirus salmonis, Caligus elongatus, and Caligus rogercresseyi.

According to another embodiment, said mutation is present in the codon of the

acetylcholine esterase gene 1 encoding the amino acid in position 362 of the ace 1 protein is determined. According to another embodiment, said codon encodes tyrosine.

According to another embodiment, wherein the mutation detected is a single nucleotid polymorphism.

According to another embodiment the present method comprise the steps of detecting a single nucleotide polymorphism (SNP) associated with organophosphate resistance in a crustacean to be analyzed, wherein said crustacean is resistant to organophosphate if comprising an acetylcholine esterase gene 1 encoding a acetylcholine esterase protein wherein the amino acid in position 362 is tyrosine.

According to another embodiment the present method comprise the steps of detecting a a single nucleotide polymorphism (SNP) associated with organophosphate resistance in a crustacean to be analyzed, wherein said crustacean is resistant to organophosphate if comprising a SL-ace l gene wherein the nucleotide in position 1085 of said SL-acel is A in the one or more crustacean to be analyzed, and wherein the numbering of said positions is in accordance with the sequence depicted in SEQ ID No. 2.

According to another embodiment the present method comprise the steps of detecting a the steps of:

a) collecting sea lice from infested fish or water samples; b) isolating genomic material from the any life stage of collected sea lice; c) determining the nucleotide polymorphic site at the positions 1085 of the isolated genomic material compared with of the nucleic acid sequence SEQ ID No. 2, wherein said sea lice is resistant to organophosphates if the nucleotide in said position is A, or a complementary oligonucleotide thereof.

According to another embodiment of the present method, said step c) is performed using a primer selected from the group consisting of SEQ ID No. 9- 12.

According to another embodiment of the present method, said step c) is performed using at least one probe selected from the group consisting of SEQ ID No. 13- 15.

According to another embodiment ofthe present method, said step c) comprises nucleic acid amplification.

According to another embodiment of the present method,the nucleic acid amplification is performed using polymerase chain reaction.

According to another embodiment of the present method, said step c) is performed by contacting the genomic material of the sea lice to be analyzed with a detection reagent, and determining which nucleotide is present in position 1085.

According to another embodiment of the present method, said step c) is performed using SNP specific probe hybridization, SNP specific primer extension, SPN specific

amplification, sequencing, 5' nuclease digestion, molecular beacon assay, oligonucleotide ligation assay, size analysis, single-stranded conformation polymorphism analysis, denaturing gradient gel electrophoresis.

According to another embodiment of the present method, the organophosphate is azamethiphos.

The present invention also provides an isolated oligonucleotide sequence comprising small nucleotide polymorphism (SNP) associated with organophosphate resistance in

crustaceans, wherein said isolated oligonucleotide sequence comprises nucleotides that distinguishes sea lice which are resistant to organophosphates from non-resistance to organophosphates, and wherein the SNP are present in the genomic DNA of the crustaceans.

According to another embodiment, the oliogonucleotide according to the present invention comprises a single-nucleotide polymorphism (SNP) associated with organophosphate resistance in sea lice, wherein said isolated oligonucleotide sequence comprises nucleotides that distinguishes sea lice which are resistant to organophosphates from non-resistant sea lice, and which enables the detection of the codon of the acel gene encoding amino acid 362 of the acetylcholine esterase protein.

According to yet antother embodiment, the the oligonucleotide sequence of the present invention is identical or has at least 80% sequence identity with a sequence selected from the group consisting of SEQ ID No. 9- 15 or a fragment thereof, and complementary sequences of SEQ ID No 9- 15 and fragments thereof. Furthermore, the present invention provides a kit for detection of organophosphate resistance in crustaceans comprising at least one oligonucleotide according to the present invention.

Finally, the present invention provides the use of an isolated oligonucleotide sequence comprising at least 8 contiguous nucleotides of the sequence SEQ ID No. 2 or a complementary oligonucleotide thereof, wherein said sequence is capable of detecting a codon encoding tyrosin in position 362 of the acetylcholine esterase protein (SEQ ID No. 6), in particular wherein the oligonucleotide sequence is capable of detecting the SNP T1085A.

The present invention and its various embodiments will be described in more detail in the following.

Figures

Figure 1 Alignment of the amino acid sequences for the acetyl cholin esterase from Torpedo californica (TC-acel) and Lepeophtheirus salmonis (LS-acel ). The number

(vertically displayed) over (Torpedo californica) or under (Lepeophtheirus salmonis) the amino acid gives the number in their respective amino acid sequence.

Figure 2 Alignment of LS-acel and LS-ace2 proteins with other typical AChE proteins from other species (Insects, Nematods, Arachnida and vertebrates). By convention, numbering is that of Torpedo californica. The three amino acids composing the catalytic triad (Ser200, Glu327 and His440) are indicated by arrows. The 14 conserved aromatic residues lining the active gorge are represented by circles. Out of these 14, 1 1 residues were present in both LS-ace l and LS-ace2 (shown by filled circles), whereas the other 3 non conserved residues (shown by open circles) were absent in both the proteins of salmon lice. The choline binding site (Trp at 84) is underlined. Three interchain disulphide bridges are drawn between conserved Cys residues (shown by arrows). The solid box represents the canonical *FGESAG* motif, characteristic of the active site of cholinesterases. The dotted box represents the typical sequence insertion/deletion domain that easily distinguishes ace l and ace2 proteins.

Figure 3 Phylogenetic relationship of LS-ace l and LS-ace2 with other

acetylcholinesterases from Insecta, Nematoda, Arachnidae and Veterbrata. The

phylogenetic tree was constructed using a ClustalW alignment and— tree (MrBayes). LS- ace l (SL acel) and LS-ace2 (SL_ace2) are boxed in the phylogenetic tree.

Figure 4 Nucleotide alignment of LS-acel from sensitive (LSAlta) and resistant (LSHitra) sea lice strain. The changes identified are boxed.

Figure 5 Nucleotide alignment of LS-ace2 from sensitive (LSAlta) and resistant (LSHitra) sea lice strain. The change identified is boxed.

Figure 6 High Resolution Melt (HRM) Analysis separated the samples based on differences in the sahpes of their melt curve that refelects the differences in their genotypes. The Green cluster represents samples homozygous for Phe362Tyr mutation, Blue cluster represents samples heterozygous for Phe362Tyr mutation and Red cluster represents the wild type samples without Phe362Tyr mutation.

Figure 7 3D model of the Lepeophtheirus salmonis protein LS-ace l using the ace protein from Drosophila melanogaster as template. The catalytic triade amino acid residues are indicated in black, and some of the "aromatic patch" amino acid residues are indicated in gray. The other amino acids are indicated as helixes, loops and sheets.

Figure 8 Positioning of some of the important amino acids in the active gorge (solid gray). The catalytic triade amino acids are displayed in black and some of the "aromatic patch" amino acids are displayed in grey. The position of the wild type Phe362 is displayed in the left panel (white) and the mutated type Tyr362 is displayed in the right panel (also white). The other amino acids are omitted. Figure 9 Plots of inhibition curves of acetyl cholin esterase from sea lice (L. salmonis) after inhibition with azamethiphos in a bi-molecular rate assay. Protein extracts from adult female sealice were incubated with one azamethiphos concentration, and aliquots were assayed for AChE activity at different time-points. The curves displayed either a mono- exponential (·) or bi-exponential (♦) decline in activity. The standard errors for the observations are indicated.

Figure 10. Plots of expected inhibition curves of acetyl cholin esterase from sea lice (L. salmonis) after inhibition with azamethiphos in a bi-molecular rate assay. SS =

homozygote sensitive parasites, SR=heterozygot resistant parasites, RR= homozygote resistant parasites.

Figure 11. SNP-analyses correlate well with known resistance status in field strains. This figure shows the average deviation values for populations of field strains of sea lice (L. salmonis) analyzed using the Real-Time PCR-assay using RNA as template. Each population is encoded by an AB-number, and had been characterized as sensitive or having reduced sensitivity to organophosphate treatment based on bioassays (as indicated). The study includes a total of 327 samples from individual sea lice, and for each population between 5 and 30 samples were analyzed. The figure shows that the Real-Time PCR-assay is suited for differentiation between sensitive and resistant sea lice populations.

Figure 12 Nucleotide sequences of the salmon louse (L. salmonis) ace l gene in the sensitive strain LSAlta. The positioning of primers and probes from Table 5 are underlined and marked in bold in both sequences, and the SNP T 1085A is marked in grey within the probe. Other nucleotide differences between the sensitive and resistant strain are marked in grey. The nucleotide numbering in this figure, and the positioning of the SNP T1085A is used for reference throughout the present specification.

Figure 13 Nucleotide sequences of the salmon louse (L. salmonis) ace l gene in the resistant strain LSHitra. The positioning of primers and probes from Table 5 are underlined and marked in bold in both sequences, and the SNP T 1085A is marked in grey within the probe. Other nucleotide differences between the sensitive and resistant strain are marked in grey. The nucleotide numbering in this figure, and the positioning of the SNP T1085A is used for reference throughout the present specification.

PAGE INTENTIONALLY LEFT BLANK

Detailed description of the invention

The present invention provides an in vitro method for determination of organophosphate resistance in crustaceans, including copepods, such as sea lice (Lepeophtheirus salmonis), based on the surprising findings that a mutation linked with resistance against

organophosphate in sea lice was found in the acethylcholine esterase 1 protein coding sequence (LS-Acel).

The present inventors have identified and sequenced two surprisingly similar genes coding for two variants of the AChE protein in Lepeophtheirus salmonis denoted LS-acel and LS- ace2. The screening for mutations revealed one mutation in the LS-ace l -gene resulting in an amino-acid change in the encoded protein which is linked to organophosphate resistance.

Several important challenges with the biological material, techniques and interpretation of the findings were encountered and overcome in the study:

1) Biological material and cultivation

Several different strains of salmon lice had to be sampled, cultivated and tested in search for the mechanism behind the azamethiphos resistance. The initial clue to reduced sensitivity was only that full-scale treatment with the product SalmosanTM comprising azamethiphos as the active ingredient had a poorer effect than expected. Reduced effect of a treatment can have several causes, as incorrect calculation of the treatment volume, incorrect calculation of the dosage, insufficient mixing of the active substance in the holding tank prior to dosing it to the cage, insufficient mixing of the stock solution in the treated cage, and/or too short holding time of the fish in the bath.The verification that the cause was insensitivity by the parasite towards azamethiphos therefore had to be done on salmon lice of that particular strain cultivated on live fish in a wet-lab. Cultivation of salmon lice in the lab can be quite challenging, as unexpected events easily can destroy the culture. As an example, oversaturation of the water with nitrogen gas destroyed some of our other salmon lice strains.

2) Enzymatic assays

The kinetic property of the mutated enzyme versus the wild type enzyme is an important clue to whether or not the enzyme is resistant towards organophosphates. An enzyme- kinetic inhibition assay technique was established in our lab some years ago (Fallang et al. 2004), but was difficult to reproduce in the current study. It is not known why this has been so difficult on these salmon lice strains. Without being bound by theory, one reason might be that the level of unspecific carboxylesterases in salmon lice is higher and more variable in salmon lice sampled during the last 2-3 years than it was when the Fallang et al. (2004) study was conducted. Unspecific esterases can cleave the substrate acetylcholine, causing interference with the assay (but not with the biological function of AChE). As the enzymatic studies on the strains in question have not given conclusive results so far, another approach to verify the correlation between phenotype and genotype had to be taken. 3) Verification of resistance

The parasites had to be tested by a bioassay method, i.e. exposure to increasing

concentrations of the agent to see how much they tolerated. As no bioassay method for salmon lice sensitivity towards azamethiphos was published, an unpublished protocol had to be used, with necessary adaptations. To be on the safe side, small scale laboratory treatments with fish infected with parasites containing only the wild-type enzyme versus fish infected with parasites containing the mutation were also successfully performed and confirmed the bioassay results. 4) Technical challenges with the identification and sequencing of ace genes

The cDNA sequences of the ace genes in salmon lice were not known prior to the present invention. In an attempt to identify homologous sequences in salmon lice genome, highly conserved acetylcholinesterase sequences from other published arthropod species was used. These sequences were further used to construct primers to amplify acetylcholinesterase in salmon lice. As small changes in the nucleotide sequences may result in non- working primers, the selection of good salmon lice-specific primers was challenging. In addition, the amplification of the S'prime and S'prime ends of the cDNA was challenging due to a high degree of homology (90%) between the two acetylcholinesterase cDNA sequences in salmon lice disclosed herein.

5) The two ace-genes

Similar to most of the other arthropods, salmon lice has two ace genes, but it came as a surprise that the two genes were very similar to each other with an amino acid homology of 84 %. Since most of the coding region of both these genes was very much homologous to each other, it was very challenging to find out whether there were one or two genes coding for ace in salmon lice. This was only discovered when the full-length cDNA was determined. This finding was unusual when comparing to ace-genes in other arthropods. Most arthropods have two ace-genes, but the amino-acid homology is normally only 30 - 50 %. The interaction between these two very similar enzymes in salmon lice is unknown. By using quantitative PCR (qPCR), no difference in the expression in neither of the ace genes between sensitive and resistant parasites have been found.

6) Interpretations

The discovery of two ace genes quite similar to each other, the influence of the Phe362Tyr mutation on the survival of the parasite during treatment and the implication of both homozygote resistant and heterozygote resistant parasites being present in the same azamethiphos-resistant salmon lice strains complicated the interpretation of the results. It was initially belived that the parasite carried two versions of the enzyme, one

azamethiphos-resistant and one sensitive version, and that resistance was caused by these being differentially regulated in sensitive and resistant parasites (Fallang et al. 2004). It was thus surpricing that the Phe362Tyr mutation in the ace l gene results in

organophosphate resistance and at the same time that both versions of ace were equally expressed in heterozygote parasites. Without being bound by theory, the dynamics of the resistance is believed to be as follows: a) Acel and ace2 are working together in the cholinergic synapses in the parasite's nervous system.

b) The wild-type acel is more effective in cleaving acetyl choline than the mutated ace l . The mutated acel has though a reasonable effect on cleaving acetyl choline.

c) During treatment with azamethiphos, wild-type ace l and ace2 are completely inhibited. This results in death of the parasite

d) When the parasite carries the Phe362Tyr mutation, either as homozygotes or as heterozygotes, a treatment with azamethiphos results in an almost complete inhibition of ace2. In heterozygote individuals, the wild-type ace l is also inhibited completely, but the mutated ace l is not significantly inhibited. The residual activity of mutated ace l is sufficient to maintain a basal activity of the enzyme in the nervous system of the parasite and thereby ensuring survival during the treatment. After treatment, the more effective wild-type acel and ace2 are regenerated. In homozygote resistant parasites, only ace2 is inhibited during treatment. Thus, these parasites will survive the treatment more easily. e) The fitness of heterozygote parasites is probably not influenced significantly, as the activity of the wild-type acel and ace2 are sufficient to maintain the enzyme activity at a normal level. The fitness of homozygote resistant parasites is probably impaired as the AChE activity is reduced.

f) Without a selection pressure of azamethiphos, the homozygote resistant parasites will over time be out-competed by the sensitive parasites. The heterozygote parasites will not be out-competed and the mutated allele(s) will therefore prevail along with the wild type allele in the population for a long time after termination of the use of the substance.

Based on the above findings, the present application provides for a method for determining organophosphate resistance in crustaceans comprising the steps of detecting a mutation associated with organophosphate resistance in the acetylcholine esterase gene 1 of the crustaceans to be analyzed.

Although organophosphate resistance linked SNP presented in the experimental data was identified in the sea lice species Lepeophtheirus salmonis, the skilled person will acknowledge, based on the teaching herein, that the present method and the present oligonucleotides may be used to determine organophosphate resistance is crustaceans, in particular copepods, in particular copepods belonging to the family Caligidae. In particular, it is to be understood that the present method and the present oligonucleotides may be used to determine organophosphate resistance in copepods affecting farmed fish, such as fish belonging to the family Salmonidae. According to one embodiment, the present method and present oligonucleotides are useful for detection of organophosphate resistance in copepods selected from the group consisting of Lepeophtheirus salmonis, Caligus elongatus, and Caligus rogercresseyi.

Throughout the application, the term "sea lice" is to be understood to mean any copepod belonging to the family Caligidae. However, whenever the term "sea lice" is used in connection with the experimental data in the present application, sea lice refers to the species Lepeophtheirus salmonis. Throughout the present specification, the numbering of amino acids in the protein is by convention from the Torpedo californica protein sequence since this was the first AChE sequence to be revealed. The T. californica and the corresponding L. salmonis amino acid number is shown in Figure 1. Furthermore, the L. salmonis amino acid number is in the present specification given in italic while the T. californica amino acid number is given as normal text. In tables and figures, the names of the respective amino acids are given using their one-letter abbrevations, while in the present specification, the amino acids are given using their full name or their three-letter abbrevations. Figure 2 presents an alignment of the amino acid sequences of the acel and ace2 proteins from a number of organisms. Figure 3 describes the phylogenetic relationship between the LS-ace l , LS-ace2 and acel and ace2 from a number of other organisms. Figure 4 shows an alignment of the nucleotide sequences of LS-acel from an organophosphate-sensitive strain (LSAlta) and an organophosphate-resistant strain

(LSHitra) with the mutation leading to the amino acid change indicated. Figure 5 is a similar alignment for LS-ace2. Table 1 lists the observed frequencies of the Phe362Tyr mutation in LS-acel in samples from the salmon lice strains not being under selection pressure by azamethiphos when collected for the study. Table 2 lists the correponding observed frequencies of the Phe362Tyr mutation in samples of the salmon lice strains being under selection pressure by azamethiphos. Table 3 lists the frequencies of the Iso433Thr mutation in LS-ace2 in the sea lice populations not being under selection pressure by azamethiphos. Figure 6 illustrates clustering of samples with different genotypes based on the differences in melt curve shapes using High Resolution Melt (HRM) analysis. Figure 7 illustrates the 3D model of the LS-ace l protein, modeled with ace from D. melanogaster as template, as helixes and coils with some of the important amino acids displayed. Figure 8 illustrates the active gorge of the LS-acel protein without and with the Phe362Tyr mutation, with some of the important amino acids displayed. Table 4 lists the association between the frequencies of the different alleles in LS-ace l , the sensitivity of the strain determined by bioassays and the treatment efficacy of

azamethiphos on two of these strains. Figure 9 is a plot of inhibition curves of acetyl cholin esterase from salmon lice after inhibition with azamethiphos in a bi-molecular rate assay. Table 5 show primers and probes useful in method according to one embodiment of the present application and table 6 shows the results of using olionucleotide sequences according to the present invention in a Real-Time PCR Assay, see example 5. Finally, table 7 list the amino acid sequences and nucleotid sequences contained in the sequence listing attached to the present application.

Based on the identification of the a SNP resulting in a mutation in the SL-Ace l gene linked with organophosphate resistance in sea lice, a method for determination of

organophosphate resistance in crustaceans, such as copepods, in particular in sea lice, e.g. Lepeophtheirus salmonis have been provided. More particular, a method for detection of organophosphate resistance is provided comprising the steps of detecting a mutation, such as a single nucleotide polymorphism (SNP) in the genomic DNA of a sea lice to be analyzed, in particular a mutation present in the acetylchiline esterase 1 gene (SL-ace l), wherein said sea lice is resistant to organophosphates if the mutation results in an amino acid substitution in posision 362.

According to one embodiment, the present invention provides to a method for detection of a single nucleotide polymorphism in the SL-ace 1 gene of a crustaceans, wherein said crustaceans is resistant to organophosphate if comprising a SL-ace l gene, wherein the nucleotide in position 1085 of SL-acel gene is A, wherein the numbering of said position is in accordance with the sequence depicted in SEQ ID No. 1.

It is to be understood that the detection of any mutation in the SL-ace l gene resulting in organophosphate resistance in crustaceans, such as copepods, such as sea lice, is covered by the present application. The skilled person will e.g. acknowledge that due to the degeneration of the genetic code, different mutations may result in the substitution of the phenylalanine 362 with a tyrosin. For example, as shown in figure 4, the triplet encoding phenylalanine 362 have changed from TTC to TAC. Due to the degeneration of the genetic code, a mutation resulting in a change from TTC to TAC would likewise result in the substitution of phenylalanine with tyrosine.

In addition to the qualitative determining a mutation involved in organophosphate resistance according to the present invention, the present invention also provides for a method for gradation of crustaceans population being susceptible of developing resistance, i.e. taking into account the composition of haplotypes that are determined within a crustaceans population.

According to the present invention, "single-nucleotide polymorphisms (SNP)" is to be understood to refer to a nucleotide sequence variation occurring when a singe nucleotide, A, T (U), C or G in the genome, or other shared sequences (e.g. RNA) or fragments thereof, differs between members of a biological species, such as between variants of Lepeophtheirus salmonis. SNPs may fall within coding sequences of genes, or non-coding regions of genes, or in the regions between the genes in a genome. The SNP identified according to the present invention fall within the genes encoding the acetylcholine esterase 1 gene.

Furthermore, as used herein, an "oligonucleotide sequence" or "nucleic acid sequence" is generally an oligonucleotide sequence or a nucleic acid sequence containing a SNP described herein, or one that hybridizes to such molecule such as a nucleic acid sequence with a complementary sequence.

The skilled person is well aware of the fact that nucleic acid molecules may be double- stranded or single-stranded, and that reference to a particular site of one strand refers, as well, to the corresponding site on a complementary strand. Thus, when defining a SNP position, reference to an adenine (A), a thymine (T) (uridine (U)), a cytosine (C) or a guanine (G) at a particular site on one strand of a nucleic acid is also to be understood to define a thymine (uridine), adenine, guanine, or cytosine, respectively, at the

corresponding site on a complementary strand of the nucleic acid molecule. Thus, reference may be made to either strand in order to refer to a particular mutation or SNP position. The oligonucleotide probes and oligonucleotide primers according to the present invention may be designed to hybridize to either strand, and SNP detection methods disclosed herein may thus also in general target either strand.

An "isolated nucleic acid" as used herein is generally one that contains at least one of the SNPs described herein or one that hybridizes to such molecule, e.g. a nucleic acid with a complementary sequence, and is separated from most other nucleic acids present in the natural source of the nucleic acid, and is thus substantially free of other cellular material.

Oligonucleotide probes and oligonucleotide primers

The present invention provides oligonucleotide probes and oligonucleotide primers that may be used for detection of the presence of the SNPs according to the present invention in DNA or RNA of a crustacean to be tested, and thus for determination of organophosphate resistance. The detection of nucleic acids present in a biological sample is widely applied in both human and veterinary diagnosis, wherein nucleic acids from e.g. pathogens present in biological samples are isolated and hybridized to one or more hybridizing probes or primers are used in order to amplify a target sequence.

One or more oligonucleotide probes may be constructed based on the teaching herein and used in hybridization based detection methods where upon the binding of the

oligonucleotides to the target sequence enables detection of the presence of at least one of the SNPs described herein if present in the sample to be tested.

The skilled person will acknowledge that an oligonucleotide probe according to the present invention may be a fragment of DNA or RNA of variable length used herein in order to detect an SNP in a target sequence, e.g. single-stranded DNA or RNA, upon hybridization of the oligonucleotide probe to complementary sequence(s) of the said target sequence to be analyzed. The oligonucleotide probe according to the present invention may furthermore be labeled with a molecular marker in order to easily visualize that hybridization, and thus detection of the SNPs disclosed herein, have been achieved. Molecular markers commonly known to the skilled person may be used, e.g. a radiolabel, and more preferably, a luminescent molecule or a fluorescent molecule enabling the visualisation of the binding of the probe(s) to a target sequence.

A oligonucleotide probe according to the present invention is able to hybridize to another nucleic acid molecule, such as the single strand of DNA or RNA originating from a crustacean to be analysed, under appropriate conditions of temperature and solution ionic strength, cf. e.g. Sambrook et al., Molecular Cloning: A laboratory Manual (third edition), 2001 , CSHL Press, (ISBN 978-087969577-4). The condition of temperature and ionic strength determine what the skilled person will recognise as the "stringency" of the hybridization. The suitable stringency for hybridisation of a probe to target nucleic acids depends on inter alia the length of the probe and the degree of complementation, variables well known to the skilled person. A oligonucleotide probe according to the present invention typically comprises a nucleotide sequence which under stringent conditions hybridize to at least 8, 10, 12, 16, 20, 22, 25, 30, 40, 50 (or any other number in-between) or more consecutive nucleotides in a target nucleic acid molecule, e.g. single- stranded DNA or RNA isolated from the crustacean to be analyzed according to the present invention. According to one embodiment, the oligonucleotide probe according to the present invention comprises about 13 to 25 consecutive nucleotides. It is to be understood that the oligonucleotide probe according one embodiment comprise one of the SNPs described herein or the complement thereof. New technology like specific Locked Nucleic Acid (LNA) hybridization probes allows for the use of extremely short oligonucleotide probes (You Y.; Moreira B.G.; Behlke M .A. and Owczarzy R. (2006). "Design of LNA. probes that improve mismatch discrimination". Nucleic Acids Res. 34 (8): e60. doi: 10.1093/nar/gkl 175. PMC 1456327. PM.ID 16670427.) According to one embodiment, probes are provided which are selected from the group consisting of SEQ ID No. 13- 15.

The present invention furthermore provides oligonucleotide primers useful for

amplification of any given region of a nucleotide sequence, in particular a region containing one of the SNPs described herein. An oligonucleotide primer according to the present invention typically comprises a nucleotide sequence at least 8, 10, 12, 16, 20, 22, 25, 30, 40, 50 (or any other number in-between) or more consecutive nucleotides.

According to one embodiment, the oligonucleotide primer according to the present invention comprises about 18 - 25 consecutive nucleotides, more preferably about 20 nucleotides.

As used herein, the term "oligonucleotide primer" is to be understood to refer to a nucleic acid sequence suitable for directing an activity to a region of a nucleic acid, e.g. for amplification of a target nucleic acid sequence by polymerase chain reaction (PCR).

According to one embodiment of the present invention, "oligonucleotide primer pairs" is provided suitable for amplification of a region of genome material comprising the SNPs according to the present invention.

The skilled person will acknowledge that an oligonucleotide primer according to the present invention may be a fragment of DNA or RNA of variable length used herein in order to detect an SNP in a target sequence, e.g. single-stranded DNA or RNA, upon alignment of the oligonucleotide probe to complementary sequence(s) of the said target sequence to be analyzed. An oligonucleotide primer according to the present invention may furthermore be labeled with a molecular marker in order to enable visualization of the results obtained. Various molecular markers or labels are available, dependent on the SNP detection method used.

An oligonucleotide primer according to the present invention typically comprises the appropriate number of nucleotides allowing that said primer align with the target sequence to be analyzed. It is to be understood that the oligonucleotide primer according to the present invention according to one embodiment comprises the SNP described herein or the complement thereof. According to one embodiment, the primers useful in order to determine organophosphate resistant sea lice is selected from the group consisting of SEQ ID No. 9- 12.

Oligonucleotide probes and oligonucleotide primers according to the present invention may be synthesizes according to methods well known to the skilled person.

The present invention furthermore relates to isolated nucleic acid sequences and variants or fragments thereof having at least 70% identity with the nucleic acid sequences depicted in SEQ ID NO. 2, or fragments thereof. The term "% identity" is to be understood to refer to the percentage of nucleotides that two or more sequences or fragments thereof contains, that are the same. A specified percentage of nucleotides can be referred to as e.g. 70% identity, 80% identity, 85% identity, 90% identity, 95% identity, 99% identity or more (or any number in between) over a specified region when compared and aligned for maximum correspondence. The skilled person will acknowledge that various means for comparing sequences are available. For example, one non-limiting example of a useful computer homology or identity program useful for determining the percent homology between sequences includes the Basic Local Alignment Search Tool (BLAST) (Altschul et al., 1990, J. of Molec. Biol., 215 :403-410, "The BLAST Algorithm; Altschul et al., 1997, Nuc. Acids Res. 25:3389-3402, , Karlin and Altschul 1990, Proc. Nat'l Acad. Sci. USA, 87:2264-68; 1993, Proc. Nat'l Acad. Sci. USA 90:5873-77).

SNP identification methods

Upon the identification of the SNPs associated with organophosphate resistance in crustacean according to the present invention, the skilled person will acknowledge that various methods commonly used in order to detect polymorphisms within a population of crustacean may be used. Both methods based on genome sequencing, hybridization and enzyme based methods are applicable for determining whether a sea louse is

organophosphate resistant in accordance with the present inventions.

Various enzyme based methods are available for the skilled person, of which a number of polymerase chain reaction (PCR) based methods are available.

For example, Perkin Elmer Life Sciences provides SNP detection kit that may be used in order to determine whether a sea louse is organophosphate resistant (AcycloPrime™-FP SNP Detection). In said method, a thermostable polymerase is used which extends an oligonucleotide primer according to the present invention by one base, then ending the oligonucleotide primer one nucleotide immediately upstream of the relevant SNP position by the incorporation of fluorescent dye-labeled terminators. The identity of the base added is then determined by the increase fluorescence polarization of its linked dye

(http://shop.perkinelmer.com/Content/Manuals/MA AcvcloPrimeSNPKit.pdf ).

Oligonucleotide primers according to the present invention useful in such a method would thus be constructed in order to facilitate the extension of the primer by one base in the position selected from the group 8605, 9035, 13957, 14017, or 14064, relative to SEQ ID No. 1.

Another enzyme based method that may be used is restriction fragment length

polymorphism (RLFP), utilizing that various, highly specific endonuclease upon digestion of the target sample, i.e. genome, results in different fragments that may be separated by gel electrophoresis.

Yet another enzyme based method that may be utilized in accordance with the present invention is the flap endocuclease (FEN) method. In said method, a structure-specific endonuclease is used to cleave a three-dimensional complex formed by hybridization with the target DNA, e.g. genomic material from sea lice, and where annealing with a target sequence comprising the SNP of interest triggers cleavage by the endonuclease (Oliver, 2005, Mutat. Res., vol. 573 (1 -2), pp. 103- 1 10). Yet another method applicable in respect of the present invention is based on the use of TaqMan® Assays (Invitrogen). In said assay the oligonucleotide primers used in order to detect an SNP is labeled in both the 5'- and the 3 ' end, i.e. with a fluorophore at the 5 'end of the oligonucleotide primer, and a quencher at the 3 '-end of the oligonucleotide primer. Upon annealing of the oligonucleotide primer with a target sequence, the Taq polymerase will extend the oligonucleotide primer and form a nascent strand, followed by degradation of the oligonucleotide primer being annealed to the target, said degradation eventually resulting in the release of the fluorophore and provide a cleavage close to the quencher. The fluorescence signal produced is proportional to the fluorophore released. Various fluorophore labels may be used, such as e.g. 6-carboxyfluorescein, tetrafluorofluorescein. As quenchers, tetramethylrhodamine or dihydrocyclopyrroloindol may be used.

Several hybridization methods for detection of SNPs are available to the skilled person, and which may be utilized in accordance with the method of the present invention. For example, the SNPs according to the present invention may be detected utilizing molecular beacon technology. According to this aspect of the present invention, oligonucleotide primers may be synthesized comprising complementary regions at each end allowing the formation of a hairpin loop, and wherein a fluorophore is attached at one end of the oligonucleotide primer, and a quenching agent is attached to the other end, and wherein fluorescence signal is produced upon binding to a DNA target of interest, i.e. genomic material isolated from the sea louse to be analyzed.

Yet another method applicable in respect of the present invention is based on DNA or RNA sequencing, which is the process of determining the precise order of nucleotides within a molecule. It includes any method or technology that is used to determine the order of the four bases (adenine, guanine, cytosine, and thymine) in a strand of DNA. The skilled person is well known with the various commonly known DNA and RNA sequencing methods that may be used according to the present invention, such as e.g. shotgun sequencing or bridge PCR sequencing.

Still another method applicable for detecting SNP is High Resolution Melting Analysis (HRM) enabling rapid detection of SNPs and determination of genetic variation within a population. The first step of a HRM protocol is often amplification of the region of interest, suing standard nucleotide sequence amplification techiches well known to the skilled person, and wherein the amplification is performed in the presense of a specialized double-stranded DNA binding dye being highly fluorescent when bound to dsDNA and poorly fluorescent in unbound state. This difference provides for the monitoring of the DNA amplification. After amplification, the target is gradually denatured by increasing the temperature in small increments, resulting in a characteristic melting profile. As the amplified DNA is denatured gradually, dye is released, thus resulting in a drop in fluorescence. As seen in example 3, HRM may be used to detect a single base change between otherwise identical nucleotide sequence regions. According to one embodiment, probes useful in high resolution melting analysis for the detection of organophosphate resistance in crustaceans, such as copepods, are provided, selected from the group consisting of SEQ ID No. 9- 10. Isolation of sea lice genomic material

The method according to the present invention may according to one embodiment involve the isolation of a biological sample from a sea lice and testing for the presence of a SNP associated with organophosphate resistance in the genome. The skilled person will acknowledge that the mutation /SNP identified according to the present invention may be detected by analyzing genomic DNA as well as genomic RNA, dependent upon the detection method used.

In order to determine whether a sea louse is organophosphate resistant in accordance with the present invention, genomic material may be isolated. Various methods for obtaining genomic material well known to the skilled person are available. The skilled person will acknowledge that any tissue (i.e. any part of the sea lice) may be used in order to extract genomic material. Furthermore, the genomic material to be analyzed according to the present invention may be obtained from sea lice of any life stages, e.g. the free swimming stages (nauplius stage I and II), the copepod stage, the pre-adult (chalimus stages 1 -4), or the adult stage (adult male or adult female). According to one embodiment, tissue removed from sea lice to be tested is maintained in 70% ethanol or other conservation liquid prior to further isolation of genomic material. DNA may be extracted from the obtained tissue using commonly available DNA extraction/isolation methods, such as e.g. DNeasy DNA Tissue Kit according to the protocol of the manufacturer

(http://lycofs01 . lvcoming.edu/~gcat- seek/protocols/DNeasy Blood & Tissue Handbook.pdf ).

SNP detection Kits

Based on the teaching herein, the skilled person will acknowledge that , based on the identified mutation/SNP and associated sequence information disclosed herein, detection reagents can be developed and used to determine any SNP described herein individually or in combination, and that such detection reagents can be readily incorporated into kits used for SNP detection known in the art. The term "kit" as used herein in the context of SNP detection reagents are intended to cover e.g. combinations of multiple SNP detection reagents, or one or more SNP detection reagents, such as oligonucleotide probe(s) and oligonucleotide primer(s) or primer sets, arrays/micro arrays of nucleic acid molecules, and beads that contain one ore more oligonucleotide probe(s), oligonucleotide primer(s) or other detection reagents useful in the method of the present invention. It is furthermore to be understood that the SNP detection reagents in a kit according to the present invention may furthermore include other components commonly included in such kits, e.g. such as various types of biochemical reagents (buffers, DNA polymerase, ligase, deoxynucleotide triphosphates for chain extension/amplification, etc.), containers, packages, substrates to which SNP detection reagents are attached., etc. necessary to carry the method according to the present invention. According to one embodiment of the present invention, a kit is provided which comprises the necessary reagents to carry out one or more assays in order to detect the SNP disclosed herein according to the method of the present invention. A kit according to the present invention may preferably comprise one or more oligonucleotide probes that hybridize to a nucleic acid target molecule (i.e. genetic material) enabling detection of each target SNP position if present in the material analyzed. Multiple pairs of probes may be included in the kit to simultaneously analyze for the presence of the SNP disclosed herein at the same time. The probes contained in the kit according to the present invention may according to one embodiment be immobilized on a carrier, such as e.g. an array or a bead.

According to one embodiment, the oligonucleotide probes are suitable for the detection of the SNP T1085A.

According to one embodiment, a kit according to the present invention comprises oligonucleotide primer(s) and optionally further SNP detection reagents useful in SNP detection methods utilizing oligonucleotide primers or primer pair(s). According to one embodiment, the kit according to the present invention comprises a forward primer and a reverse primer for amplifying a region containing a a mutation in the triplet encoding amino acid 362 of the SL-ace l gene. Said kit may furthermore optionally comprise further SNP detection reagents (enzymes and nucleotide triphosphates) necessary for conducting PCR or real time PCR. According to one embodiment, the primer pairs are suitable for the detection of the SNP T1085A. According to yet another embodiment, the kit according to the present invention comprises primer pairs suitable for detection of all the SNPs described herein.

Organophosphates as used herein refers to a group of insecticides acting on the enzyme acetylcholine esterase. According to one embodiment of the present invention, a method is provided for detection of the organophosphate azamethiphos (S-[(6-Chloro-2- oxo[l ,3]oxazolo[4,5-b]pyridin-3(2H)-yl)methyl] Ο,Ο-dimethyl phosphorothioate, CAS No. 35575-96-3). It is well known that resistance against one type of organophosphate in general provides the possibility of resistance also against all type of organophosphates; see e.g. Sevatdal S, Fallang A, Ingebrigtsen K, Horsberg TE. 2005. Monooxygenase mediated organophosphate detoxification in sea lice (Lepeophtheirus salmonis). Pest Manag Sci. 2005 Aug;61 (8):772-8, Du Y, Nomura Y, Satar G, Hu Z, Nauen R, He SY, Zhorov BS, Dong K. 2013. Molecular evidence for dual organophosphate-receptor sites on a mosquito sodium channel. Proc Natl Acad Sci U S A. 2013 Jul 2. [Epub ahead of print], and Zhu F, Gujar H, Gordon JR, Haynes KF, Potter MF, Palli SR. 2013. Bed bugs evolved unique adaptive strategy to resist organophosphate insecticides. Sci Rep. 2013;3 : 1456. doi:

10.1038/srep01456.

In health, agriculture, and government, the word "organophosphates" refers to a group of insecticides or nerve agents acting on the enzyme acetylcholinesterase. The skilled person will thus acknowledge, based on the teaching of the present invention, that the method according to the invention may be used to determine resistance towards various types of organophosphates in crustaceans, such as e.g. azametiphos, parathion, malathion, methyl parathion, chlorpyrifos, diazinon, dichlorvos, phosmet, fenitrothion, tetrachlorvinphos, and azinphos methyl. Examples

Example 1 : Identification and sequencing of acetylcholine esterase genes in

Lepeophteirus salmonis Samples of salmon lice

Salmon lice samples with different histories of sensitivity / resistance towards the organophosphate azamethiphos were collected in the field and cultivated in the laboratory at the NIVA Marine Research Station at So lb erg strand, Drobak or at the Institute of Biology, University of Bergen.

The strains were:

LSAlta: Sampled 2010, inbred for approx. 10 generations. No history of insensivity towards azamethiphos.

LSBjugn: Sampled 2008, inbred for approx. 15 generations. History of insensitivity towards azamethiphos.

LSGulen: Sampled in 2004, inbred for approx. 25 generations. No history of

insensivity towards azamethiphos.

LSHitra: Sampled in 201 1 , inbred for approx. 7 generations. History of insensitivity towards azamethiphos.

LSHitra-s: Subset of LSHitra where fish with parasites were treated with 100 ppb azamethiphos (Salmosan™) for 30 minutes in a 200 liter fibre-glass tank.

The subset is the surviving parasites after treatment.

LSVeso: Sampled in 2012 in the field (Otteroya) shortly after a full-scale treatment with azamethiphos in the farm. The site had been treated with azamethiphos at several occations prior to sampling. Only surviving parasites were sampled.

Total RNA extraction and cDNA synthesis

Total RNA was extracted using RNeasy Mini kit (Qiagen, CA, USA), from female adult individuals, as per manufactures 's protocol. The RNA was quantified and qualified on ND-100 Spectrophotometer (Thermo Fisher Scientific, DE, USA). First strand cDNA was synthesized from total RNA (lug) using qScript reverse transcriptase (Quanta Biosciences, MD, USA).

Partial cDNA fragments of LS-acel and LS-ace2

Conserved sequences of AChEs were selected from other species using GenBank database, which were then compared against the sealice genome database (Viroblast; sealouse.imr.no) to detect the homologous sequences in sea lice.

Rapid amplification of cDNA ends

5' and 3' ends of partial cDNAs, obtained by homology blast, were amplified using Rapid amplification of cDNA ends (RACE) with sequence specific primers (mentioned below) and SMART RACE kit (Clontech, Palo Alto, CA, USA) as per manufacturer's instructions. RACE PCR was performed under conditions: 5 cycles at 94°C for 30 s, 72°C for 3 min followed by 5 cycles at 94°C for 30 s, 70°C for 30 s, 72°C for 3 min followed by 25 cycles at 94°C for 30 s, 68°C for 30 s and 72°C for 3 min.

Primers used for 5 'RACE

LS-acel primer AATCCCCTCCCAGGCATCAATGTGTCTA

LS-ace 2 primer CATCCCATCCCAGGCATTAATGTCTTTG Primers used for 3 'RACE

LS-acel primer CCTATTGTGGATGGAAGTTTCTTGGATGAG

LS-ace 2 primer CAACCAGGAGGGGAGTGGTACCGTCTT

Both 5'RACE and 3'RACE PCR products were cloned using TOPO TA Cloning Kit for sequencing (Invitrogen, CA, USA) followed by isolation of plasmid DNA from positive colonies using PureLink Quick Plasmid Miniprep kit (Invitrogen, CA, USA) under manufactures 's instructions. Amplicons were obtained using the plasmid DNA and TOPO vector specific primers (mentioned below) under PCR conditions: 94°C for 4 min followed by 30 cycles at 94°C for 1 min, 52°C for 1 min, 72°C for 1 min and followed by final extension at 72°C for 5 min.

Amplicons were then subjected to direct sequencing using BIG Dye Terminator v3.1 cycle sequencing kit (Life technologies, Invitrogen, CA, USA) on 3130x1 Genetic Analyzer (ABI Prism, Life technologies, Invitrogen, CA, USA) to obtain the full length sequence of cDNAs.

TOPO vector specific primers used

M13Forward primer: GTAAAACGACGGCCAG

M13Reverse primer: CAGGAAACAGCTATGAC Screening of full length LS-acel and LS-ace2

Full length cDNAs (LS-acel and LS-ace2), from 3 sensitive (LSAlta) and 3 resistant (LSHitra) adult female sea lice samples, were amplified using gene specific primers (mentioned below). PCR reactions were performed using Phusion high- fidelity DNA polymerase (New England BioLabs, MA, USA) under the conditions: 98°C for 30 s, followed by 35 cycles at 98°C for 10 s, 55°C for 15 s, 72°C for 2 min followed by a final extension at 72°C for 10 min. Amplicons were then subjected to direct sequencing using BIG Dye Terminator v3.1 cycle sequencing kit (Life technologies, Invitrogen, CA, USA) on 3130x1 Genetic Analyzer (ABI Prism, Life technologies, Invitrogen, CA, USA). Primers used to amplify the whole cDNA

LS-acel forward primer: CTCTGCTGCTACACCGACTCCTGTT

LS-acel reverse primer: TCGAGGATGTTTGACACTGATGGTC

LS-ace2 forward primer: TGTTTTAGATGTGGATTCAAGTCCGAA

LS-ace2 reverse primer: CGATGGATGGTACGTACGTATGAACATA

The deduced amino acid sequence of both LS-acel and LS-ace2 was aligned with previously published acetylcholinesterases (AChE) from other insects, nematodes, arachnida and vertebrates (Insects: Liposcelis entomophila acel (Lip_ent_acel),Liposcelis entomophila ace2

(Lip_ent_ace2), Bemisia tabaci acel (Bern Jab _acel), Bemisia tabaci ace2 (Bern Jab _ace2), Blattella germanica ace 1 (Bla_ger_acel), Blattella germanica ace 2 (Bla_ger_ace2),

Nephotettix cincticeps acel (Nep in icel), Nephotettix cincticeps ace2(Nep in ice2), Ctenocephalides felis acel (Cte Jel icel), Ctenocephalides felis ace2 (Cte Jel ice2), Culex pipiens ace(Cul j)ip), Chilo suppressalis acel (Chi_sup icel), Chilo suppressalis ace2

(Chi_sup ice2), Apis mellifera acel (Api nel icel), Apis mellifera ace2 (Api_mel ice2), Cimex lectularius acel (Cimjec icel), Cimex lectularius ace2 (Cimjec ice2), Bombyx mandarina acel (Bom man acel), Bombyx mandarina ace2 (Bom nan ice2), Bombyx mori acel (Bom Mor acel), Bombyx mori ace2(Bom_Mor_ace2), Leptinotarsa decemlineata

(Lep_dec) Drosophila (Dros), Musca domestica (Mus_dom), Anopheles gambiae (Ano gam), Aedes albopictus (Aed ilb), Culex quinquefasciatus (Cul_qui), Nematods: Caenorhabditis elegans acel (Celeg_acel), Caenorhabditis elegans ace2 (Celeg_ace2), Caenorhabditis briggsae acel(Cae_bri_acel), Caenorhabditis briggsae ace2 (Cae_bri_ace2), Arachnida: Tetranychus urticae acel (Tet_urt_acel), Rhipicephalus decoloratus (Rhi_dec), Vertebrates: Torpedo californica (Tor_cal), Homo sapiens (Homosap)), using ClustalW alignment.

Surprisingly, both the LS-ace proteins exhibited high degree of homology with each other (84%). This is very much different from other arthropods, where only a moderate or low level of homology has been reported between acel and ace2 genes. For example, 35.6%, 39% and 53% similarity has been reported between Cimex lectularius, Chilo suppressalis and Anopheles gambiae acel and ace2, respectively (Seong et al, 2012; Jiang et al, 2009; Weill et al, 2002).

Besides, both the LS-ace proteins exhibited the typical feature of lacking 32 amino acids at the sequence insertion/deletion domain, which has been reported to be a common feature of all the insect acel type genes (Kim et al, 2006). All these observations indicate the possibility that two different forms of acel gene are present in sea lice unlike most of the other arthropods that have two different ace (acel and ace2) genes. This could possibly be due to duplication of a single ace gene in sea lice.

Moreover, both of these proteins showed highest homology to acetylcholinesterase- 1 proteins from Blattella germanica (63%) followed by Liposcelis entomophila (62%), Nephotettix cincticeps (62%), Cimex lectularius (60%), Bombyx mori (60%), Leptinotarsa decemlineata (60%), Bemisia tabaci (59%), Chilo suppressalis (59%) and Bombyx mandarina (58%).

However, only moderate level of cross homology was observed between LS-ace 1 and LS-ace2 with acetylcholinesterase-2 from different species, ranging from 45% (Chilo suppressalis) to 33% (Blattella germanica). In addition, both LS-acel and LS-ace2 showed only 44% homology to Drosophila melanogaster AChE and 46% to Torpedo californica AChE.

The protein alignment also revealed that both LS-acel and LS-ace2 have the characteristic features of AChE, including the anionic choline binding site at Trp84, the three residues of the catalytic triad (Ser200, Glu327 and His440), the six cysteines potentially involved in three conserved disulphide bonds (67-94; 254-265; 402-521), the characteristic *FGESAG* motif surrounding the active serine and the 14 aromatic residues lining the active site gorge, 11 of which were well conserved and present in both LS-acel and LS-ace2 (Figure 2). This includes the acyl pocket residues (Trp233, Phe290 and Phe331) that accommodate the acyl moiety of the active site. In addition, the residues (Glyl 18, Glyl 19 and Ala201) that form the oxyanion hole, that helps to stabilize the tetrahedral molecule during catalysis, was also present in both the proteins. The 3 non conserved amino acids (70, 121 and 442) were absent in both the salmon lice ace proteins.

Interesting, most of the other species (including LS-acel and LS-ace2) have aspartic acid at 442 instead of tyrosine that is present in Torpedo californica (Figure 2). Moreover, Phe290 is present and Phe288 is absent in both the sequences, a characteristic property of all invertebrate AChE sequences, explaining a wider substrate specificity than vertebrate AChE (Vellom et al. 1993).

To construct the phylogenetic tree, the ClustalW alignment of LS-ace 1 and LS-ace2 along with AChEs from 35 different species (already deposited in GenBank) was subjected to MrBayes software. The phylogenetic tree suggested that both LS-acel and LS-ace2 were clustered with other AChE-typel proteins, and were clearly separated from AChE-type 2 proteins that form a separate clad in the phylogenetic tree (Figure 3). Example 2: Identification of missense changes in LS-acel and LS-ace2

Missense changes in LS-acel and LS-ace2 were screened in 50 samples each from 2 sensitive (LSAlta and LSGulen) and 2 resistant populations (LSHitra, LSBjugn), by direct sequencing. These samples were not under any organophosphate treatment pressure when enrolled. In addition, 2 L. salmonis populations (24 samples of LSHitra-s and 20 samples of LSVeso) that survived the organophosphate treatment were also screened for these missense changes by direct sequencing. The sensitivity status of all the L. salmonis populations has been mentoned in detail above (Samples of salmon lice).

Primers to amplify missense change in LS-ace 1

LS-acel forward primer: GTGGATGGAAGTTTCTTGGATGAGAG

LS-acel reverse primer: CTCAAAAGTTATTGCCTCTCTTCCCATT Primers to amplify missense change in LS-ace2

LS-ace2 forward primer: ACGAGCAAAGTCAGCAGTTG

LS-ace2 reverse primer: TTTCATCCGCAGTGTTTCAG

To determine the role of LS-acel and LS-ace2 with resistance, whole cDNA sequence of both the genes was screened in 3 sensitive (LSAlta) and 3 resistant (LSHitra) sea lice organisms.

Sequencing revealed three single nucleotide substitutions in LS-acel . Of these 3 changes, only one led to an amino acid change (Phenylalanine to Tyrosine at codon 362, which corresponds to codon 331 in the Torpedo californica amino acid sequence), whereas the other 2 substitutions were silent changes (Figure 4). All the three changes were identified only in the resistant salmon lice organisms.

A single change was identified in LS-ace2 that led to Isoleucine to Threonine substitution at codon 433 (Figure 5).

Both the non-synonymous changes (Phe362Tyr in LS-acel and Iso433Thr in LS-ace2) were screened, by direct sequencing, in lab cultured sea lice populations, including 2 sensitive strains (LSAlta, LSGulen) and 2 resistant strains (LSBjugn and LSHitra) to determine their association with resistance against organophosphates (azamethiphos). 50 samples from each population were enrolled for screening. None of these populations were under any treatment when enrolled.

Besides, 20 highly resistant samples that survived higher doses of azametiphos along with 24 first generation samples from a resistant (LSHitra) population that survived the azametiphos treatment were also screened for the two changes. Table 1: Frequencies of Phe362Tyr in LS-acel in the salmon lice populations without any treatment.

Table 2: Frequencies of Iso433Thr change in LS-ace2 in the salmon lice populations without any treatment.

Screening revealed a higher frequency of Phe362Tyr in LS-acel in the resistant population indicating its association with resistance against azemetiphos in sea lice (Table 1). Moreover, the frequency observed was very much higher in the samples that survived the azametiphos treatment (Table 2). This observation clearly indicates an association between Phe362Tyr and survival of the sea lice under azametiphos treatment. Moreover, Phe362 (Phe331 in T. californica) is very much conserved throughout the species as evident from the multiple sequence alignment of LS- acel and LS-ace2 with acetylcholinesterases from other species (Figure 2).

However, no such association was observed with Iso433Thr change in LS-ace2 (Table 3). Thus, this change was considered to be a polymorphism. Table 3: Frequencies of the Iso433Thr change in LS-ace2 in the salmon lice populations without any treatment.

Example 3: High Resolution Melt analysis (HRM)

The Phe362Tyr mutation in LS-ace l was validated by HRM, which is a simple rapid tool to screen single base changes (mutations/polymorphisms) with high sensitivity and accuracy. The methodology included the generation of specific PCR product using gene specific primers (mentioned below) and Precision Melt supermix (Bio Rad, CA, USA), as per manufactures 's instructions, with a sensitive fluorescent dye (EvaGreen) that binds specifically only to double stranded DNA, followed by subjecting the amplicon to gradual increase in temperature (65°C to 95°C), which led to the denaturation of double stranded amplicon and decrease in fluorescence. This change in florescence was recorded by the C I 000 Touch thermal cycler (Bio-Rad, CA, USA) as a melt curve (fluorescence versus temperature). The samples were assembled into different groups based on difference in the shapes of their melt curves. HRM primers for the Phe362Tyr mutation

Forward primer: TTTTAATTG G AG CG AATAAG G A (SEQ ID No. 9)

Reverse Primer: TCTGTTCGATCAACATAGACG (SEQ ID No. 10)

High Resolution Melt analysis (HRM) was performed to validate the sequencing results and in an attempt to develop a rapid diagnostic tool for the detection of Phe362Tyr mutation in LS-acel . After standardising the technique with samples of known genotypes determined by direct sequencing (wild type, heterozygous and homozygous for Phe362Tyr mutation), samples with unknown genotypes were run to confirm the results obtained. These were also confirmed by direct sequencing. HRM analysis could distinguish between the samples of different genotypes with high accuracy. As shown in Figure 6, the samples were very well separated based on their genotypes. Example 4: 3D modeling of the enzymes

The three-dimentional structure of the enzyme was modeled using SWISS MODEL in the automated mode (Arnold et al. 2006), http://swissmodel.expasv.org/. The models were generated using native acetyl choline esterase (AChE or ace) from Torpedo californica and Drosophila melanogaster as templates. The templates were generated on basis of the crystaline structure of the T. californica or D. melanogaster ace proteins, determined by X- ray diffraction (Harel et al. 2000, PDB-ID lqo9; Harel & Sussman (unpublished), PDB-ID 2j 3 d) . Through the process, 3D files (pdb-files) of the salmon louse protein showing the best fit with the templates from D. melanogaster and T. californica were generated.

The 3D modeling was done using SWISS MODEL (Arnold et al. 2006),

http://swissmodel.expasy.org/. As templates, the 3D structure of the native enzyme from Drosophila melanogaster and Torpedo californica were tested. These 3D structures were based on X-ray diffraction. The protein from L. salmonis could be fitted to both templates, but the fits were not optimal. The QMEAN4 score (a parameter between 0 and 1 where a higher number indicates a better fit) was 0.541 when fitted against the AChE protein from Drosophila melanogaster and 0.614 when fitted against the protein from Torpedo californica. Even though the fit was slightly better with the T. californica protein, the results here describe the fit against the Drosphila protein. This because Drosophila melanogaster is another arthropod (insect) while Torpedo californica is a vertebrate (fish). The esteratic subsite of the active gorge of the acetylcholine esterase enzyme consists of a triade of three amino acids, serine (Ser), glutamic acid (Glu) and histidine (His). These are in T. californica located in positions Ser200, Glu327 and His440, in the LS-ace l protein in positions Ser230, Glu358 and His472. They are responsible for cleaving the acetylcholine molecule.

The catalytic triade is located almost in the center of the protein. A groove leads from the surface to the triade. Several aromatic amino acid residues (the "aromatic patch") are lining the gorge, binding the choline part of the AChE molecule and providing an exit route for the degradation products after hydrolysis. The residues considered to be of highest importance are Trp84 (Trpll5), Tyrl 30 (Tyr 161), Phe330 (Tyr 361) and Phe331 (Phe362) (Pezzementi et al. 2003). One of these aromatic residues, Phe331 (Phe362) is also considered a part of the "histidine trap", fixing the histidine at position 440 (His472) in its correct position during hydrolysis of acetyl choline. Several other amino acids are also of importance for binding the substrate in the active gorge. The amino acids at the esteratic subsite, the "aromatic patch" and the anionic subsites of the enzyme (not listed) are highly conserved across species (Ordentlich et al. 1993).

In the figure, the positions of the three amino acids forming the catalytic triade of the enzyme, Ser230, Glu358 and His472 within the whole protein are indicated in black. The "aromatic patch" amino acid residues (Trpl 15, Tyrl61, Tyr361 and Phe362) are indicated in gray. The other amino acids are indicated as coils, sheets and helixes.

The mutation described here leads to a substitution of the amino acid Phe362 to Tyr362. The aromatic ring of Tyr is turned approximately 50 degrees compared to the aromatic ring of Phe. Tyrosine has a hydroxyl group in the para-position, which enters the groove, decreasing the volume of the pocket from 743 A ³ (wild type) to 714 A ³ (mutated type) according to the model. The substitution of the nonpolar Phe with the more polar Tyr changes the polarity of the active gorge and thereby also the binding site for

organophosphates, probably affecting binding of these molecules in the enzyme.

In addition to being a part of the "aromatic patch", the amino acid Phe in position 331 (Phe362) participates in aromatic trapping of histidine in position 440 (His472) in the catalytic triade. Phe331 is important for maintaining histidine in its correct position for catalysis (Nabeshima et al. 2004). The model estimated the shortest distance between Phe362 and His472 (wild type) to be 3.52 A, and the shortest distance between Tyr362 and His472 (mutated type) to be 3.20 A. So possibly, the Phe362Tyr mutation interferes directly with the catalytic activity of histidine in the catalytic triade by leaving less room for the acetyl choline molecule. This may affect the normal catalytic function of the enzyme, thus posing a fitness cost associated with the mutation.

The Phe362Tyr mutation has not been described as naturally occurring in any other species. It has though been investigated by site-directed mutagenesis in Drosophila melanogaster (Phe371Tyr, Drosophila amino acid sequence, Boublik et al. 2002). The mutation was demonstrated to decrease the sensitivity of the enzyme to the

organophosphate paraoxon by 200-fold. Table 4: Sensitivity tests of different strains of sea lice (L. salmonis) and the frequency of alleles, given as S/S (homozygote Phe362), R/R. (homozygote Tyr 362) and S/R (heterozygote Phe362 / Tyr 362).

The sensitivity of the strains analyzed was tested by traditional bioassays and small scale treatments (Helgesen et al., 2013). The history of the strains were

LSAlta: Sampled 2010, inbred for approx. 10 generations. Fully sensitive to all agents

LSBjugn: Sampled 2008, inbred for approx. 15 generations. Moderately resistant to azamethiphos, pyrethroids and emamectin benzoate.

LSGulen: Sampled in 2004, inbred for approx. 25 generations. Fully sensitive to all agents LSHitra: Sampled in 201 1. Resistant towards azamethiphos (not tested for other agents). Inbred for approx. 7 generations.

LSHitra-s: Subset of LSHitra where fish with parasites were treated with 100 ppb azamethiphos (Salmosan™) for 30 minutes in a 200 liter fibre-glass tank. The subset is the surviving parasites after treatment.

LSVeso: Sampled in 2012 in the field (Otteroya) shortly after a full-scale treatment with azamethiphos in the farm. Only surviving parasites were sampled.

The table demonstrates a high negative correlation between the sensitivity (EC ₅ o values) and the frequency of S/S alleles in the samples. The correlation between EC ₅₀ s and the S/S frequency was

Person 's r (with 95 % confidence interval) = -0.9687 (- 0.9980 - -0.5947).

Since the confidence interval does not contain the value zero, the correlation is statistically significant (p<0.05). As the relationship between dose and effect is non-linear, the data were also correlated using the non-paramethric Spearman's p. This was also statistically significant.

Spearman 's p = -0.9000 (p<0.0374)

The small scale treatments demonstrated 100 % efficacy of 100 ppb azamethiphos towards the fully sensitive LSAlta strain, and 50 % efficacy towards the resistant LSHitra strain. This demonstrated that when the R allele was present in some samples, the treatment efficacy was lower than when it was not present. Since the mean value from one group was not included in the 95 % confidence interval of the other group, the treatment effects were statistically significantly different (p<0.05). In samples collected in 2000 - 2002 from various parts of Norway, the inhibition of the salmon louse acetylcholin esterase by azamethiphos was determined (Fallang et al. 2004). The study has not yet been repeated for samples collected in 2010 - 2013.

In the previous study, the sea lice were homogenized, centrifuged and filtered. Each well in a column was prepared with aliquots of 30 μΐ of the sample supernatant, 0.12 mg DTNB and buffer up to 140 μΐ. Azamethiphos was added to the first well giving a final concentration of 1.6 μΜ. After 30 min of equilibration at 25°C, azamethiphos giving a final concentration of 0.01 μΜ, was added to wells 2-7. The same amount of buffer was added to well 8. Aliquots of 0.4 mg acetylthiocholine (substrate) were then added to the first two wells, after which the reading (405 nm, each min for 30 min) was started. The same amount of substrate was added to wells 3-6, respectively, after 1 , 3, 5 and 10 min.

After 20 min, substrate was added to wells 7 and 8. The final volume in each well was 200 μΐ. The maximum slope for each well, recorded within 10 min after addition of azamethiphos and substrate, corrected for spontaneous hydrolysis of the substrate, was used in the calculation.

The enzyme activity was calculated for each well. The enzyme activity was followed down to approx. 10% of the activity in the uninhibited sample, after which the readings became too unreliable. The percent remaining AChE activity in each single sample was plotted against the incubation time using a semi-logarithmic plot. When the curves in the semi- logarithmic plot formed a straight line, they were considered being best described by a mono-exponential function. The line for the average value was fitted using the following equation:

Curves not forming a straight line in the plot were considered being best described by a bi- exponential function, and the line for the average value was fitted using the following equation:

where E is the enzyme activity in the well, Ei and E ₂ are the zero-time intercepts for the exponential terms, k _! ^(obs) and k ₂ ^(obs) are the observed inhibition rate constants, e is the base of the natural logarithm and t is the time after addition of the inhibitor azamethiphos. The concentration-independent bimolecular rate constants, k;i and k _i2 were calculated by dividing ki ^(obs) and k ₂ ^(obs) with the concentration of azamethiphos used in the assay.

The figure demonstrates that the enzyme activity in some samples was inhibited completely within 2 minutes (mono-exponential decline rate), most likely as a result of the

acetylcholin esterase being of the wild type. In other samples, the acetyl cholin esterase was inhibited approx. 80 % by azamethiphos within 2 minutes, but the activity remained at this level (bi-exponential decline rate), most likely as a result of heterozygot resistant acetylcholin esterase. In homozygot resistant parasites, a very slow decline rate corresponding to the last part of the curve for the heterozygot parasites would have been expected. This was not seen in the samples collected in 2000 - 2002.

Example 5: Testing of organophosphate resistance

Based on the identification of the SNP responsible for organophosphate resistance in L. salmonis, we developed a sensitive Real-Time PCR (TaqMan) 5 '-nuclease assays for single nucleotide polymorphism (SNP) detection (Real-Time PCR SNP-assay) (Table 5), and successfully applied it to differentiate between organophosphate resistant and sensitive sea lice. Fluorogenic PCR probes, chemically modified with a minor groove binding agent to increase duplex stability, were used in single and multiplex probe closed tube formats. The probe were tested in a commercially available thermocycling fluorimeter (Applied

Biosystems 7500 Real-Time PCR System) at PatoGen Analyse AS laboratory in Alesund. The assay were used qualitatively to determine the relative amount of the SNPs in individual samples, and samples from different life stages of sea lice were use. Comparison of results obtained using this assay showed no discrepancies from results obtained by genome sequencing. The reported SNPs and the Real-Time PCR SNP-assay is ideal for the differentiating between organophosphate resistant and sensitive sea lice in the aquaculture industry.

Sampling of sea lice

Sea lice samples were collected from sea lice with known sensitivity-status to

organophosphate. The sea lice collected in fish farms had been tested with respect to sensitivity to organophosphate by bioassays as previously described. Sea lice used for comparison between genome sequencing and the Real-Time PCR-assay had been tested using bioassays and the genome sequence were known from previous genomic sequencing performed as described in example 1 above. In addition, sea lice were sampled from fish farms on the west coast of Norway. Sea lice were collected using forceps, and approximately 10-50 lice per site were conserved in 70% ethanol and kept at 4 ^° C. Samples were sent refrigerated to PatoGen by express mail carrier.

Nucleic acid purification

In PatoGens laboratory, RNA and/or DNA were extracted from samples by methods well known to the skilled person. In short, tissue samples were transferred to Micro Collection Tubes and lysed and homogenized using QIAzol Lysis Reagent, steel beads and vigorous shaking using a TissueLyser system, followed by nucleic acid extraction using an RNAeasy kit (Qiagen) or DNAeasy kit (Qiagen), all according to the manufacturer's instructions, and by methods well known to the skilled person. Chloroform were added to the samples and shaken vigorously. After resting and centrifuging, the relevant liquid phase were collected for further extraction of either RNA or DNA by vacuum technology using a Qiagen robot system, all according to the manufacturer's instructions, and by methods well known to the skilled person. Finally, nucleic acids were eluted in 25 ml of elution buffer and used for PCR by methods well known to the skilled person.

SNP -detection by Real-Time PCR

The primers and probes listed in table 5 are TaqMan® MGB Probe SNP Genotyping Assays using TaqMan® 5 ' nuclease assay chemistry for amplifying and detecting specific SNP alleles in purified genomic DNA or RNA. The primers and probe used to identify the SNPs is listed in table 5, and were ordered from Life Technologies Corporation. One-step amplification (45 cycles) was performed on an Applied Biosystems 7500 Real-Time PCR System performed at PatoGen Analyse AS laboratory in Alesund, all according to the manufacturer's instructions, and by methods well known to the skilled person. The assay was used qualitatively to determine the relative amount of the relevant SNPs in individual samples.

Table 5: Primers and probes The primers and probe used to identify the SNPs in salmon lice. Probes are both listed with IUPAC nucleotide codes for SNPs, and with the alternative nucleotides for each position as indicated. IUPAC codes used are in accordance with Nomenclature for Incompletely Specified Bases in Nucleic Acid Sequences, of which only one is relevant here: W=A or T (Nomenclature Committee of the International Union of Biochemistry (NC-IUB) (1984) http://www.chem.qmul.ac.uk/iubmb/misc/naseq.html).

Results

The result from Real-Time PCR-assay is interpreted by looking at the deviation in Ct- values between the probes detecting the two variants of the SNP. Sensitive sea lice are the type homozygote SS, and have a deviation higher than zero. Genetically resistant sea lice have a deviation lower than zero. Genetically resistant sea lice can be either heterozygote (RS) or homozygote (RR), and the heterozygote (RS) has a deviation between zero and -3,4, while the homozygote resistant (RR) has a deviation lower than -3,4. In addition, there is a quantitative aspect where the highly resistant strains have a lower deviation value than moderately resistant sea lice. Most likely, genetically resistant strains that have not been exposed to organophosphates for some time has a lower deviation value than genetically resistant strains that been repeatedly and newly exposed to

organophosphates.

The results from the Real-Time PCR SNP-analyses showed good correlation with results obtained by genome sequencing, and correlated with the known resistance status in sea lice populations (Table 6).

Also, analyses performed on a higher number of samples field strains showed good correlation with the known resistance status and as shown in Figure 1 1.

The Real-Time PCR-SNP-assay showed surprisingly good results, and we believe that it will serve as a practical tool in the differentiation between organophosphate sensitive and resistant sea lice.

Thus the reported SNP and the relevant Real-Time PCR SNP-assay represent an ideal tool for differentiating between organophosphate resistant and sensitive sea lice in the aquaculture industry. Also, the prevalence of sensitive versus resistant sea lice in a population, and the deviation value for individual lice or average deviation values for populations, can be used to predict the best possible outcome of a treatment using organophosphates in the population. Also, the technique can become an important tool for optimizing sea lice treatments using organophosphates, combination treatments using organophosphate in combination with other insecticides or other measures to reduce sea lice infection pressure, and to monitor the resistant status of populations of sea lice before treatment.

Table 6. Conclusions regarding genotype based on the Real-Time PCR-assay showed good correlation with genotyping. The figure shows that the Real-Time PCR-assay is suited for differentiation between genotypes related to resistance towards

organophosphates.

Lice population Parallel Deviation Conclusion

3 SS

SS I 2 3 SS

SS2 liis!iiiiiiii 12,8 SS

SS2 2 11,2 SS

SS4 4,6 SS

SS4 2 4,1 SS

SS5 »i»lil ftJI§W 8,8 SS

SS5 2 4,7 SS

RS2 -0,9 RS

RS3 1 -0,8 RS

RS3 l ls llJPiw!i -0,6 RS

RS4 1 -1,2 RS

RS4 ftiiSiiilillfSil -1,4 RS

RS5 1 -1,2 RS

RS5 l#llllSfcl¾5 -1,3 RS

RS6 1 -1,7 RS

RS6 -1,9 RS

Average -1,09

RR1 1 -4,7 RR

RR1 2 -4,5 RR

RR2 ifSf!lililllfi _* ' -4,4 RR

RR2 2 -4,5 RR

RR4 -4,7 RR

RR4 2 -4,2 RR

RR5 IfliijJlllllil -4,9 RR

RR5 2 -4,6 RR

Average -4,56

Table 7: Sequences related to the invention

SEQ ID No. Description of sequence

1. SEQ ID Acetylcholinesterase gene sequence (LS-ace l) isolated from

No. 1 organophosphate sensitive sea lice strain (LS-Alta)

2. SEQ ID Acetylcholinesterase gene sequence (LS-ace l) isolated from

No. 2 organophosphate restistant sea lice strain (LS-Hitra)

3. SEQ ID Acetylcholinesterase gene sequence (LS-ace2) isolated from

No. 3 organophosphate sensitive sea lice strain (LS-Alta)

4. SE ID Acetylcholinesterase gene sequence (LS-ace2) isolated from

No. 4 organophosphate restistant sea lice strain (LS-Hitra)

5. SEQ ID Acetylcholinesterase amino acid sequence encoded by the

No. 5 acetylcholinesterase gene sequence (LS-acel , SEQ ID No. 1) isolated from organophosphate sensitive sea lice strain (LS-Alta)

6. SEQ ID Acetylcholinesterase amino acid sequence encoded by the

No. 6 acetylcholinesterase gene sequence (LS-acel , SEQ ID No. 2) isolated from organophosphate sensitive sea lice strain (LS-Hitra)

7. SEQ ID Acetylcholinesterase amino acid sequence encoded by the

No. 7 acetylcholinesterase gene sequence (LS-ace2, SEQ ID No. 3) isolated from organophosphate sensitive sea lice strain (LS-Alta)

8. SEQ ID Acetylcholinesterase amino acid sequence encoded by the

No. 8 acetylcholinesterase gene sequence (LS-ace2, SEQ ID No. 4) isolated from organophosphate sensitive sea lice strain (LS-Hitra)

9. SEQ ID forward primer used in example 3

No.9

10. SEQ ID reverse primer used in example 3

No.10

1 1. SEQ ID Forward primer used in example 5

No.1 1

12. SEQ ID Forward primer used in example 5

No.12

13. SEQ ID Probe used in example 5 in accordance with IUPAC nucleotide code No. 13

14. SEQ ID Probe used in example 5 (sensitive)

No. 14

15. SEQ ID Probe used in example 5 (resistant)

No. 15 References

Arnold et al. (2006). Bioinformatics 22, 195-201

Boublik et al. (2002). Protein engeneering 15, 43-50

Costello M. J. (2009). Journal of Fish Diseases. 32. 115-118

Denholm et al. (2002). Pest Manag Sci 58, 528-536

Fallang et al. (2004). Pest Manag Sci, 60, 1163-1170

Harel et al. (2000). Protein Science 9, 1063-1072

Helgesen et al. (2013). J Fish Dis 36, 261-272

Jiang et al. (2009). Biochem Biophys Res Comm 378, 269-272

Kim et al. (2006). Insect Mol Biol 15, 513-522

Nabeshima et al. (2004). Biochem Biophys Res Comm 313, 794-801

Ordentlich et al. (1993). J Biol Chem 268, 17083-95

Pezzementi et al. (2003). Comp Biochem Physiol B 136, 813-832 Pike A., Wadsworth S. L., (2000), Academic Press. 44. 232-337).

Seong et al. (2012). Insect Mol Biol 21, 149-159

Sevatdal S., et al. (2005), Aquaculture 244. 19-27

Torrissen et al. (2013). J Fish Dis 36, 171-194

Vellom et al. (1993). Biochemistry 32, 12-17

Weill et al. (2002). Proc Roy Soc B Biol Sci 269, 2007-2016