Salmon lice (Lepeophtheirus salmonis) showing varying emamectin benzoate susceptibilities differ in neuronal acetylcholine receptor and GABA-gated chloride channel mRNA expression
© Carmichael et al.; licensee BioMed Central Ltd. 2013
Received: 12 March 2013
Accepted: 13 June 2013
Published: 18 June 2013
Caligid copepods, also called sea lice, are fish ectoparasites, some species of which cause significant problems in the mariculture of salmon, where the annual cost of infection is in excess of €300 million globally. At present, caligid control on farms is mainly achieved using medicinal treatments. However, the continued use of a restricted number of medicine actives potentially favours the development of drug resistance. Here, we report transcriptional changes in a laboratory strain of the caligid Lepeophtheirus salmonis (Krøyer, 1837) that is moderately (~7-fold) resistant to the avermectin compound emamectin benzoate (EMB), a component of the anti-salmon louse agent SLICE® (Merck Animal Health).
Suppression subtractive hybridisation (SSH) was used to enrich transcripts differentially expressed between EMB-resistant (PT) and drug-susceptible (S) laboratory strains of L. salmonis. SSH libraries were subjected to 454 sequencing. Further L. salmonis transcript sequences were available as expressed sequence tags (EST) from GenBank. Contiguous sequences were generated from both SSH and EST sequences and annotated. Transcriptional responses in PT and S salmon lice were investigated using custom 15 K oligonucleotide microarrays designed using the above sequence resources. In the absence of EMB exposure, 359 targets differed in transcript abundance between the two strains, these genes being enriched for functions such as calcium ion binding, chitin metabolism and muscle structure. γ-aminobutyric acid (GABA)-gated chloride channel (GABA-Cl) and neuronal acetylcholine receptor (nAChR) subunits showed significantly lower transcript levels in PT lice compared to S lice. Using RT-qPCR, the decrease in mRNA levels was estimated at ~1.4-fold for GABA-Cl and ~2.8-fold for nAChR. Salmon lice from the PT strain showed few transcriptional responses following acute exposure (1 or 3 h) to 200 μg L-1 of EMB, a drug concentration tolerated by PT lice, but toxic for S lice.
Avermectins are believed to exert their toxicity to invertebrates through interaction with glutamate-gated and GABA-gated chloride channels. Further potential drug targets include other Cys-loop ion channels such as nAChR. The present study demonstrates decreased transcript abundances of GABA-Cl and nAChR subunits in EMB-resistant salmon lice, suggesting their involvement in avermectin toxicity in caligids.
KeywordsDrug resistance Sea lice Avermectin Ligand-gated chloride channel Cys-loop receptor
The development of synthetic insecticides and parasiticides has revolutionised the control of arthropod pests and parasites of agricultural, medical and veterinary importance. However, the overuse of chemical control agents tends to promote the development of heritable insecticide or drug resistance , with life-cycle traits of targeted organisms often favouring the ability to adapt genetically, such as high reproductive potential and short generation span. Drug resistance can follow from pre-existing mutations in resistance genes, termed resistance alleles, which under normal conditions are rare and have a negative or neutral effect on fitness but which, under conditions of exposure to control agents, afford fitness benefits and can become enriched in the target species’ gene pool given persisting selection pressure from the control agent . By the time resistance becomes apparent as treatment failure, resistance alleles have usually already reached high frequencies in the gene pool .
Caligid copepods, also called sea lice, are common ectoparasites of marine fish . One species, the salmon louse (Lepeophtheirus salmonis (Krøyer, 1837)), has emerged as a serious problem in mariculture of Atlantic salmon (Salmo salar Linnaeus, 1758) in the Northern hemisphere. The annual cost of sea louse infection to the global salmon farming industry has been estimated at €300 million, with the majority of this accounted for through expenses accrued from treatments with veterinary medicines . Only a limited range of anti-sea louse drugs are available and licensed for the treatment of fish, and the continued use of a relatively small number of compounds creates a situation potentially favouring the development of drug resistance . In the salmon louse, losses of efficacy have been reported for a number of control agents including organophosphates , pyrethroids , hydrogen peroxide  and avermectins (AVMs) [9, 10].
The commonly used anti-sea louse treatment SLICE® (Merck Animal Health) contains the avermectin compound emamectin benzoate (EMB) (Stone et al., 1999). SLICE® is administered orally and a one-week treatment provides prolonged protection against all host-attached life stages of sea lice . Avermectins are also used against external and internal parasites of humans and livestock, including parasitic nematodes causing the human diseases onchocerciasis (River blindness) and lymphatic filariasis, as well as gastrointestinal parasites of sheep, cows and horses . The selective toxicity of avermectins against ecdysozoan invertebrates is believed to be based on the binding and blockage of glutamate-gated (GluCl) and γ-aminobutyric acid (GABA)-gated (GABA-Cl) chloride channels in the invertebrate nervous system . Several molecular mechanisms have been suggested as contributing factors to the resistance of parasitic nematodes to the AVM compound ivermectin (IVM) . Functional studies revealed that resistant nematodes can have single amino acid mutations in subunits of GluCl and GABA-Cl that decrease the channels’ sensitivities to the drug [15, 16]. Furthermore, resistant nematodes may show increased expression of ABC (ATP-binding cassette) transporters, a group of membrane proteins with members capable of mediating the cellular efflux of drugs [17, 18]. Finally, avermectin resistance in insects has been connected to alterations in drug metabolism .
Previous studies on potential molecular mechanisms of EMB resistance in salmon lice have used the candidate gene approach, i.e. the study of genes that have previously been linked to drug resistance in other organisms. In particular, such studies have investigated salmon louse ABC transporters [9, 20] and GABA-Cl and GluCl subunits . However, potential outcomes resulting from such candidate approaches are limited by the current knowledge of the biological process being studied. To avoid this limitation, a transcriptomic approach was followed in the present study, in which microarray analysis was used to compare mRNA responses between drug susceptible and moderately (~7-fold) EMB resistant laboratory strains of salmon lice.
Custom microarray design
Composition of features on custom L. salmonis oligo microarrays
Experiment 1 (AMADID # 039612)
Experiment 2 (AMADID # 033382)
Analysis of strain differences in constitutive gene expression
Genes showing differential expression between salmon louse strains differing in EMB susceptibility
Fold change (PT/S)
Fold change (PT/S)
Ligand gated ion channel (<1%)
Neuronal acetylcholine receptor subunit alpha-3
Metabolism of xenobiotics (<1%)
Cytochrome P450 3A24
Regulation of synapse development (<1%)
Eye lens proteins (<1%)
Solute carrier family 8 (sodium/calcium exchanger)
Alpha subunit of Na+/K+ ATPase
Excitatory amino acid transporter 3
Cytoskeleton proteins (26%)
Muscle LIM protein
Fast myosin heavy chain
Myosin light chain alkali
Myosin light chain alkali
Collagen alpha-1 chain
Myosin light chain
Troponin C, isoform 1
Troponin C, isoform 1
Troponin C, isoform 1
Troponin C, isoform 1
Torso-like protein precursor
Troponin C, isoform 1
Regulation of actin cytoskeleton (<1%)
Calcium transport (<1%)
Sarco/endoplasmic reticulum calcium transporting ATPase
Calcium binding (1%)
Sarcoplasmic calcium-binding protein, beta chain
Cuticle proteins (4%)
Cuticle protein 6
Cuticle protein CP14.6 precursor
Serine carboxypeptidase CPVL precursor
Intestinal trypsin 2 precursor
Intestinal trypsin 4 precursor
Gene expression measured by RT-qPCR in salmon lice from two strains differing in EMB susceptibility
Fold change (PT/S)
Fold change (PT/S)
Fold change (PT/S)
Fold change (PT/S)
Neuronal acetylcholine receptor subunit α3
GABA receptor subunit alpha
Cytochrome p450 18a1
Cytochrome p450 15a1
Effects of short-term exposure to EMB on transcript profiles in salmon lice
Number of differentially expressed features identified in microarray analysis of samples from experiment 2
Number of features
Strain x Treatment
Strain x Time
Time x Treatment
Strain x Treatment x Time
The most significantly differentially regulated genes from cluster 1, and the responses to EMB exposure in both salmon louse strains are summarised in Additional file 2. An analogous selection of genes from cluster 2 is provided in Additional file 3. Genes in cluster one include the GABA-Cl and nAChR α-3 subunits which are potential targets for EMB, with nAChR α-3 also having been identified as constitutively differentially expressed between salmon louse strains. Genes in cluster two included a glutathione-S-transferase isoform and a nAChR α-3 precursor. Enrichment analysis of the features in cluster one is detailed in Additional file 4, showing significant over representation of twelve and under representation of four GO attributes. Chitin binding, calcium ion binding and hydrolase activity were the most significantly over-represented attributes and nucleic acid binding was identified as the most significant under represented attribute.
Using transcriptional profiling in comparative studies of a drug-susceptible and an EMB-resistant salmon louse strain, this study demonstrated the reduced constitutive mRNA expression of subunits of certain ligand-gated ion channels (LGIC) in the EMB resistant strain, namely a GABA-gated chloride channel subunit (GABA-Cl, ~1.4-fold decreased) and a neuronal acetylcholine receptor subunit (nAChR α-3; ~2.8-fold decreased). The toxicity of EMB and other AVMs against ecdysozoan invertebrates is reported to be based mainly on their interaction with another class of LGICs, the glutamate-gated chloride channels (GluCl) [23–25] although GABA-Cls are also believed to be pharmacological targets of AVMs . While nAChRs are traditionally not considered to be implicated in the toxic action of AVMs in ecdysozoans, they can be allosterically modulated by IVM . This study’s finding that mRNA levels of GABA-Cl and nAChR subunits are decreased in EMB-resistant salmon lice suggests that these LGICs may represent potentially additional target sites for AVMs in sea lice. When transcriptomic profiles of EMB-resistant and -susceptible laboratory salmon louse strains were investigated in response to short term (1–3 h) aqueous EMB exposure, a number of transcriptional responses to the treatment were observed in the drug-susceptible S lice, but few responses were found in the EMB-resistant PT strain. While we cannot exclude the possibility that EMB exposure might provoke more pronounced transcriptomic responses in PT lice at later time points, differential toxicity of EMB between the louse strains is apparent as early as 5 hours post-exposure (data not shown). This suggests that the mechanism of resistance, whether constitutive, responsive or a combination of both, must have occurred before this time point. The present study only considered levels of mRNA. Molecular mechanisms underlying differential susceptibilities between the studied strains could also include post-transcriptional regulatory mechanisms, such as mRNA processing and degradation, translation and protein degradation , however, these lay outside the scope of the present study.
AVMs are used against parasitic and pest species of ecdysozoan invertebrates including nematodes, insects and mites, and more recently against crustaceans. While the molecular target sites of AVMs in crustaceans are unknown, GluCls are generally considered to be the main pharmacological targets of IVMs in nematodes and insects [29, 30]. The GluCls form an invertebrate-specific subgroup of the large ‘Cys-loop’ subfamily of LGIC. ‘Cys-loop’ LGICs have a pentameric structure, and are composed of either the same type of subunits or two to three different subunit isoforms. Consistent with the role of GluCl as the main target of AVMs, IVM-resistant strains of invertebrates can show mutations changing the expression levels or the peptide sequence of channel subunits [23, 24]. A GluClα subunit has been cloned in L. salmonis and while GluClα was represented amongst the microarray targets used for this study, no difference in mRNA expression was observed between salmon lice of the two studied strains, or amongst those subjected to control and sublethal EMB treatments (data not shown).
Apart from GluCl, further LGICs are known to interact with AVMs. For instance, IVM modulates the activity of nematode GABA-Cl , and can exert directly activating or potentiating effects on vertebrate glycine-gated chloride channels . Moreover, AVMs can modulate the activity of cation-LGICs such as the α-7 nAChR  and the ATP-gated P2X4 receptors . A number of observations involving drug-resistant insects and nematodes support the hypothesis that LGICs other than GluCl constitute further toxicologically relevant targets of AVMs in invertebrates. Cyclodiene-resistant fruit flies having a single amino acid mutation in a GABA-Cl showed a moderate degree of cross-resistance to IVM . A null mutation in a histamine-gated chloride channel also conferred moderate IVM resistance in Drosophila melanogaster Meigen, 1830 , and a novel dopamine-gated ion channel (HcGGR3) was significantly down-regulated in an AVM-selected strain of the nematode Haemonchus contortus (Rudolphi, 1803) . The observation in this study that EMB-resistant salmon lice show decreased mRNA levels of nAChR and GABA-Cl is consistent with findings in the literature cited above, and suggests a role for nAChR and GABA-Cl as additional pharmacological targets of EMB in salmon lice. It is worth noting in this context that observed changes in nAChR expression could also relate to previous exposure of PT lice to compounds interfering with cholinergic neurotransmission such as the organophosphate (OP) anti-sea louse drug azamethiphos . However, decreases in nAChR expression are not among typical molecular responses associated with OP resistance in insects [1, 36]. While decreased expression of nAChR has been observed in nematodes resistant to imidazothiazoles and other drugs that directly target nAChR channels , these are classes of compound that have not been used against sea lice.
Apart from modifications of the molecular targets, biocide resistance in pests and parasites can result from increased elimination of the chemical as a result of up-regulation of biotransformation enzymes and/or drug transporters. The superfamily of cytochrome P450s (CYPs) contains heme-thiolate proteins that function as monooxygenases, many of which are involved in drug metabolism [37, 38]. CYPs play important roles as chemical resistance genes in insects [1, 39], but their roles in the biochemical defence against toxicants in crustaceans are less well understood. The microarrays used in the present study included probes representing a number of target sequences that are annotated as CYPs (experiment 1: 18 probes; experiment 2: 14 probes). These partial L. salmonis CYP sequences could not be unequivocally attributed to specific CYP families, as this would require establishing the overall degree of amino acid similarity to CYP family members . However, based on BLASTx annotation (e-values of ≤ 10-7), most of these sequences could be provisionally allocated to CYP clans, which are higher-order groupings that combine phylogenetically related CYP families  (Clan 2: 6 sequences; clan 3: 5 sequences; mitochondrial clan: 3 sequences). Two targets showing similarity to CYP15A1 and CYP18A1 (both clan 2) differed in mRNA expression in L. salmonis from experiment 1, but not experiment 2. Moreover, compared to the EMB susceptible strain a target showing similarity to CYP3A24 (clan 3) had lower expression levels in the EMB resistant strain, with a moderate (1.48-fold) difference in transcript abundance found between the strains. In insects, CYP15A1 is involved in juvenile hormone synthesis , whereas CYP18A1 functions to inactivate ecdysteroids . Crustaceans possess homologues to both these clan 2 CYPs , but little is known of their function. In the green shore crab (Carcinus maenas) expression levels of two CYPs from clan 2 were affected by both the moulting cycle and previous exposure to xenobiotics . The differences in CYP mRNA levels found in this study were relatively small and/or variable between experiments and therefore do not provide clear evidence for an involvement of CYPs in the differential EMB susceptibility found in the salmon louse strains that were studied.
Carboxylesterases are another class of enzymes that can confer insecticide resistance . In this study, the expression of one carboxylesterase was moderately enhanced in EMB-resistant salmon lice in experiment 1, but no significant differences in expression were observed between strains in experiment 2. Accordingly, the data provides no evidence for a role of carboxylesterase in EMB resistance of salmon lice.
ABC (ATP-binding cassette) proteins are a family of membrane-bound transporters mediating the transport of a diverse array of substrates across biological membranes . Certain ABC proteins are drug efflux transporters located in the cell membrane, and have roles in the biochemical defence against toxicants . The ABC transporter P-glycoprotein transports IVM  and has relevance as a biochemical factor limiting the drug’s toxicity in mice and nematodes [18, 48]. It has been suggested that P-glycoprotein could be implicated in the resistance of pests and parasites to AVMs [49–51], and a salmon louse homologue of P-glycoprotein called SL-PGY1 has been cloned . In the present study, mRNA expression of SL-PGY1 was unaffected by EMB treatment and did not differ between S and PT strains (data not shown), confirming similar data previously reported for these strains .
In addition to comparing transcriptional profiles between louse strains in the absence of drug exposure, the present study also addresses transcriptomic responses of both louse strains to short-term exposures (1 to 3 hours) to 200 μg L-1 of EMB. The selected EMB concentration is not toxic to the resistant salmon louse strain, but in the drug-susceptible strain it results in close to 100% immotility after 24 hours exposure. After 1 to 3 hours of exposure, few transcriptomic responses were observed in the resistant strain, and therefore these data do not provide evidence for the presence of specific resistance mechanisms requiring transcriptional regulation. In the drug-susceptible strain, a complex array of transcriptional responses to EMB was observed. This finding was expected, as this level of EMB exposure was known to result in significant toxic responses after 24 hours. Early transcriptional responses are therefore likely to comprise both adaptive and general stress responses, as well as transcriptional changes reflecting the toxic action of the drug.
The transcriptional profiling of a drug susceptible and an EMB resistant laboratory strain of salmon lice revealed reduced mRNA expression of a number of LGIC subunits in the EMB resistant lice. GluCl, considered to be the major target site of AVM drugs in invertebrates, displayed similar mRNA expression levels in EMB-resistant and reference strains. In contrast, subunits of GABA-Cl and nAChR showed decreased mRNA abundances in the EMB resistant compared to the reference strain. While GABA-Cl is considered a secondary target of AVMs in invertebrates, nAChR is not traditionally considered a target site for AVMs, even though it has been shown to interact with AVMs in vertebrates. It is possible that nAChR and GABA-Cl represent additional EMB target sites in salmon lice, and that the down-regulation of these channel subunits in this EMB-resistant strain could be related to the resistance phenotype. In the present study no changes were seen in the expression levels of biotransformation enzymes and drug transporters, both of which classes have been suggested to contribute to AVM resistance in other species. Further studies are needed to investigate potential relationships between the transcriptional changes observed and the susceptibility phenotype.
Salmon louse strains
Two laboratory-maintained salmon louse (L. salmonis) strains differing in susceptibility to EMB , were used in this study. The susceptible strain S was established in 2003 using salmon lice from a Scottish farm site where no chemical control agents other than hydrogen peroxide had been used. The moderately EMB-resistant salmon louse strain PT was established in December 2008 using salmon lice from another Scottish production site with reports of variable EMB treatment efficacies. The strains have since been cultured under identical laboratory conditions, as described in detail elsewhere . In brief, salmon lice were maintained on Atlantic salmon (S. salar) with an initial weight of 500–1000 g in circular tanks supplied with fresh seawater at ambient temperature, using a photoperiod corresponding to natural day length. To propagate salmon louse cultures, egg strings were allowed to hatch and develop to copepodids, which were used to inoculate a tank containing fresh host fish. Prior to the collection of salmon lice from hosts, fish were anaesthetised with 100 mg L-1 2-phenoxyethanol. Infection rates were maintained at levels that were unlikely to compromise fish welfare. All laboratory infections were carried out under UK Home Office licence and appropriate veterinary supervision.
Salmon louse exposure experiments
Adult male salmon lice were collected from anaesthetised host fish as described above and allowed to recover for 2 hours in aerated filtered seawater at ambient sea temperature. To analyse transcript expression in salmon louse strains S and PT in the absence of drug exposure (microarray experiment 1), adult males were collected and preserved in an RNA stabilisation solution (4.54 M ammonium sulphate, 25 mM trisodium citrate, 20 mM EDTA, pH 5.4) prior to storage at -70°C. For both strains six pooled samples comprising four individuals each were taken. The effects of acute seawater-borne exposure to EMB on transcript expression were investigated in microarray experiment 2. In order to establish the EMB concentration for this experiment, 24 hour bioassays  were carried out on salmon louse strains S and PT. EMB (technical grade, a gift from Merck Animal Health) was solubilised in seawater with PEG300 (final concentration 0.01% (v/v)). EMB concentrations used in bioassays were 32.5, 75, 150, 300 and 600 μg L-1 with S lice, and 200, 400, 800, 1,600 and 3,200 μg L-1 with PT lice. Duplicate glass dishes containing 10 salmon lice and 200 mL of exposure solution were used per EMB concentration, control (seawater) or solvent control (seawater with 0.01% (v/v) PEG300). At the end of 24 hours of exposure, salmon lice were recorded as normally motile or immotile upon visual examination and stimulation with a fine brush. To generate samples for microarray experiment 2, salmon lice of either strain were subjected to short term (1 and 3 hour) exposures to 200 μg L-1 EMB, a concentration that would result in >95% immotility in S lice after 24 hours but have no effects in PT lice. In addition, the experiment comprised seawater and solvent (0.01% (v/v) PEG300) controls. For each combination of strain, exposure period and treatment, three pooled samples consisting of four salmon lice each were collected for later RNA extraction as above. None of the treatments had effects on louse motility. At the end of the experiment (3 hours), water samples were taken and sent to a commercial laboratory (Eclipse Scientific, Chatteris, UK) for EMB residue analysis (liquid chromatography with detection by MSMS). The measured EMB concentration in the nominal 200 μg L-1 EMB treatment was 99.5 ± 5.2 μg L-1 EMB. This depletion of solubilised active ingredient may be attributed to EMB adsorption to the glass containers used for exposure assays .
RNA extraction and purification
In microarray and RT-qPCR experiments, samples were pools of four adult male salmon lice. Frozen samples were ground in liquid nitrogen using a pestle and mortar, and total RNA was immediately extracted from the homogenised sample using TRI Reagent® (Sigma-Aldrich, UK), following the manufacturer’s protocol. After phase separation, RNA was precipitated from the aqueous phase by addition of 0.25 volumes isopropanol and 0.25 volumes of a high salt buffer (0.8 M trisodium citrate; 1.2 M sodium chloride), as recommended for samples with high polysaccharide content . The total RNA was resuspended in nuclease-free water and further purified using RNeasy columns (Qiagen, UK). For the construction of subtracted cDNA libraries, total RNA from 60 untreated adult males from either strain (S or PT) were pooled and subjected to poly (A)+ RNA isolation using the Poly (A) Purist™ kit (Ambion®, UK). UV spectroscopy (NanoDrop ND-1000, Thermo Scientific, USA) was used to confirm purity of the RNA samples and establish concentrations, whereas RNA integrity was assessed by agarose gel electrophoresis and ethidium bromide staining.
Subtracted cDNA library construction and sequencing
Suppression subtractive hybridisation (SSH) was used to prepare cDNA libraries enriched in transcripts differentially expressed between strains S and PT using commercial methods (PCR-Select™, Clontech, Takara Bio Inc., USA) . Subtractions were performed in both directions, i.e. using cDNA derived from each strain (S or PT) either as the tester or the driver. A pool of cDNA from each subtraction, containing an equal amount of both subtracted cDNA libraries, was used for generating a 454 sequencing library using the GS FLX Titanium Rapid Library Preparation kit (Roche Applied Science, UK), following manufacturer’s instructions. Adaptive Focus Acoustics™ (AFA™) using the S220 High Performance Ultrasonicator (Covaris® Inc., KBiosciences, UK) was employed to randomly shear the cDNA, blunt ends were repaired and MID adapters ligated to the DNA fragments prior to sequencing using the Genome Sequencer™ (GS) Titanium FLX instrument (Roche Applied Science, UK) (EBI Sequence Read Archive (SRA) study ERP002190). GS FLX Titanium library preparation and sequencing was performed by The GenePool Genomics Facility (University of Edinburgh, UK). Sequence reads were assembled using the GS De Novo Assembler (Newbler) v2.5.3 software (Roche Applied Science, UK) using default parameters after trimming of MID adapter and primer sequences.
Sequence assembly from existing L. salmonisEST resource
A total of 129,225 sequences (> 100 bp) were downloaded during December 2010 for L. salmonis from the GenBank EST database (as described by Yasuike et al.) , and assembled into contigs using default assembly settings of The Gene Indices Clustering Tools (TGICL), obtained from the Computational Biology and Functional Genomics Laboratory (The Gene Index Project, USA). Prior to sequence assembly, vector sequences were removed using SeqMan II 6.1 (DNAStar Inc., USA).
Salmon louse microarray design
The assembled contig sequences were annotated using BLASTx (Basic Local Alignment Search Tool) searches against the non-redundant proteins (nr), UniprotKB/Swiss-Prot (Swissprot) and Reference Proteins (Refseq_Proteins) GenBank databases at the National Centre for Biotechnology Information (NCBI), with an annotation hit having an expectation value (e-value) of <1 × 10-4 being considered significant. All sequences were further annotated with GO identifiers using Blast2Go software for Windows® using Java Webstart (Centro de Investigación Príncipe Felipe, Spain). Oligonucleotide probes (60mers) were designed to target contig sequences using the eArray Gene Expression (GE) probe design tool (Agilent Technologies, UK), employing the base composition and best probe methodologies, and designed in sense orientation with 3’ bias. For each sequence without a significant BLASTx based annotation two probes were selected; designed to both forward and reverse complement sequences. Standard expression microarrays were designed using the eArray custom microarray design wizard (Agilent Technologies, UK) for an 8 × 15 K design format. Each microarray comprised 15,744 features including 536 obligatory controls (Table 1). An initial design was used for experiment 1 interrogations (Agilent AMADID No 039612; EBI ArrayExpress design A-MEXP-2285). This design incorporated probes designed to target 2,699 sequences that were identified when sequencing the subtracted cDNA libraries enriched for transcripts differentially expressed between the EMB-resistant (PT) and drug susceptible (S) salmon louse strains. Two probes were designed for each of the SSH targets. Experiment 2 employed a modified design (Agilent AMADID No 033382; EBI ArrayExpress design AMEXP-2284). The array designs shared 10,251 identical features.
Labelling protocols are described in detail elsewhere . Briefly, for each test sample 250 ng total RNA was used as template for the amplification of antisense RNA with the incorporation of the modified nucleotide 5-(3-aminoallyl)-UTP (aaUTP) into the amplified RNA (aRNA) during the in vitro transcription step (TargetAmp™ Aminoallyl-aRNA Amplification Kit 101; Epicentre®, Cambio Ltd. UK). A common reference pool was created through pooling equal amounts of all aRNA test samples to be used in the experiment. The individual test samples were labelled with cyanine 3 (Cy3) and the common reference pool labelled with Cy5 mono-reactive dye (GE Healthcare, UK) in dye coupling reactions. Unincorporated dye was removed by column purification (Illustra Autoseq™ G-50 spin columns; GE healthcare, UK), and then dye incorporation was assessed by spectrophotometry (NanoDrop ND-1000, Thermo Scientific, USA) and fluorescent gel electrophoresis. Three hundred nanograms of each Cy3-labelled test sample was competitively hybridised with 300 ng Cy5-labelled common reference pool on a 15 K feature custom microarray, following the manufacturer’s instructions (Agilent Technologies, UK). The hybridisation reactions were incubated at 65°C with 10 rpm rotation for 17 hours in an Agilent rotary hybridisation oven and then washed with Gene Expression Wash Buffers 1 and 2, with a final wash using Stabilisation and Drying solution, again following the manufacturer’s instructions (Agilent Technologies, UK). The hybridised microarrays were scanned using an Axon Genepix 4200A scanner with Genepix Pro 6.1 image acquisition software (Molecular Devices, UK) using 40% laser power, 5 μm pixel size resolution and auto photo-multiplier tube (auto-PMT) function with 0.05 saturation tolerances. The raw microarray images were processed using Agilent Feature Extraction (FE) software version 18.104.22.168 that performed feature grid alignment, extraction and quantification. The fluorescence intensity results files from the FE software were imported into the GeneSpring GX version 12 software (Agilent Technologies, UK) for differential gene expression analysis. Data were normalised using Lowess normalisation of log2-expression ratios without baseline transformation. Features showing low quality according to Agilent quality control metrics were discarded to provide a final feature set for analysis. Details of microarray experiments 1 and 2 have been submitted to ArrayExpress and assigned accession numbers E-MTAB-1484 (Experiment 1) and E-MTAB-1478 (Experiment 2). The recording of the microarray experimental metadata complies with Minimum Information About a Microarray Experiment (MIAME) guidelines.
To validate gene expression results from microarray experiments, the abundance of selected differentially expressed transcripts was determined by reverse transcription quantitative PCR (RT-qPCR). Per experiment, three targets that showed stable expression levels in microarray hybridisations were selected as reference genes. (Experiment 1: 60S ribosomal protein S20, 40S ribosomal protein L44 and RMD-5 homolog; experiment 2: Hypoxanthine-guanine phosphoribosyltransferase (HGPRT), RMD-5 homolog and Elongation factor 1α). For each target sequence, primers were designed with a melting temperature (Tm) of ~60°C using Primer 3 software (Additional file 5). Aliquots (1 μg) of total RNA samples previously used in microarray analyses were reverse transcribed (Superscript III, Invitrogen, UK) using random hexamers and anchored oligo-dT in a 3:1 molar ratio. No-template controls and controls omitting RT enzyme were included on each assay plate to detect potential DNA contamination. A cDNA pool containing equal amounts of all samples was made and included on each assay plate, serving as a calibration sample (20-fold dilution) and for derivation of a standard curve from serial dilutions. RT-qPCR reactions were performed in duplicate in a total volume of 20 μL containing 5 μL sample cDNA (20-fold dilution), 0.3 μM of each primer and 10 μL Absolute SYBR Green I mix (ThermoFisher Scientific, UK), using the Mastercycler ep realplex2 (Eppendorf, UK) with the following amplification conditions: 95°C for 15 minutes, followed by 40 cycles of 94°C for 30 seconds, 15 seconds at the specific primer pair annealing temperature (Ta; Additional file 5) and 72°C for 30 seconds. After amplification a melt curve from 55°C to 95°C at 0.5°C increments for 15 seconds each was performed to ensure that a single product was amplified in each reaction. Threshold cycles were analysed using the PCR cycler software. Standard curves were derived from plots of the threshold cycle against the logarithm of the relative concentration of cDNA pool. Primer efficiency (E) was derived from linear fits to the standard curve according to the equation E = 10(-1/slope). The BestKeeper tool  was employed to analyse expression stability of three reference genes and determine a robust BestKeeper expression index as a geometric mean for the three reference genes, which was in turn used to establish relative gene expression ratios using the ΔΔCt method using the Relative Expression Software Tool (REST) Multiple Condition Solver (MCS) .
Microarray gene expression data were analysed using GeneSpring GX version 12 (Agilent Technologies). The analysis of constitutive differential gene expression in experiments 1 and 2 used Student’s t-test adapted for samples with unequal variance (Welch) using a fold change threshold of 1.3. The analysis of differential gene expression induced in microarray experiment 2 employed two-way ANOVA to compare exposure of both salmon louse strains at two time points against control conditions. Multiple testing corrections were not applied to any statistical analysis of this gene expression study as this can often be over-conservative when studying potentially subtle gene expression responses to stimuli [56, 59]. This decision is supported by confirmation of differential expression by RT-qPCR in the current study. Network analysis of microarray experiment 2 expression data was performed using the BioLayout Express3D application . A network graph was constructed using the Pearson correlation coefficient (threshold of 0.94) to determine similarities between expression profiles, which were then arranged into groups of features with similar profiles using the Markov clustering algorithm (MCL) with the default inflation setting (2.2) for optimal clustering. Gene enrichment analysis was performed on lists of features chosen based on differential gene expression patterns using default settings of the FuncAssociate 2.0 web application . Gene enrichment was calculated according to the significance (p< 0.05) of the association between the list of features and the GO attributes represented on the microarray. Relative expression ratios from RT-qPCR experiments were tested for normality and equal variance and log transformed to allow assumptions to be satisfied before being subjected to one way ANOVA using Minitab 16.1 software (Minitab Inc., UK). The significance level was set at p<0.05 in all tests.
This project was financed by the Scottish Salmon Producers Organisation (SSPO), The University of Stirling and the Moredun Foundation. We are grateful to all the staff at the Marine Environmental Research Laboratory (MERL) for provision of the salmon lice used in this study. We would also like to thank Anna Montazam and Stephen Bridget from the GenePool Genomics Facility for assistance with Roche 454 sequencing and assembly.
- ffrench-Constant RH, Daborn PJ, Le Goff G: The genetics and genomics of insecticide resistance. Trends Genet. 2004, 20: 163-170. 10.1016/j.tig.2004.01.003.View ArticlePubMedGoogle Scholar
- Wolstenholme AJ, Fairweather I, Prichard R, von Samson-Himmelstjerna G, Sangster NC: Drug resistance in veterinary helminths. Trends Parasitol. 2004, 20: 469-476. 10.1016/j.pt.2004.07.010.View ArticlePubMedGoogle Scholar
- Kabata Z: Parasitic copepoda of British fishes. 1979, London: British Museum (Natural History) for the Ray SocietyGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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 Kroyer) against the pyrethroid deltamethrin using bioassays and probit modelling. Aquaculture. 2003, 218: 21-31. 10.1016/S0044-8486(02)00339-3.View ArticleGoogle Scholar
- Treasurer JW, Wadsworth S, Grant A: Resistance of sea lice, Lepeophtheirus salmonis (Kroyer), to hydrogen peroxide on farmed Atlantic salmon, Salmo salar L. Aquacult Res. 2000, 31: 855-860. 10.1046/j.1365-2109.2000.00517.x.View 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. Comparative Biochemistry and Physiology - C Toxicology and Pharmacology. 2012, 155: 198-205. 10.1016/j.cbpc.2011.08.004.View ArticlePubMedGoogle 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
- Stone J, Sutherland IH, Sommerville CS, Richards RH, Varma KJ: The efficacy of emamectin benzoate as an oral treatment of sea lice, Lepeophtheirus salmonis (Kroyer), infestations in Atlantic salmon, Salmo salar L. J Fish Dis. 1999, 22: 261-270.View ArticleGoogle Scholar
- Geary TG: Ivermectin 20 years on: maturation of a wonder drug. Trends Parasitol. 2005, 21: 530-532. 10.1016/j.pt.2005.08.014.View ArticlePubMedGoogle Scholar
- McCavera S, Walsh TK, Wolstenholme AJ: Nematode ligand-gated chloride channels: an appraisal of their involvement in macrocyclic lactone resistance and prospects for developing molecular markers. Parasitology. 2007, 134: 1111-1121. 10.1017/S0031182007000042.View ArticlePubMedGoogle Scholar
- Beech RN, Skuce P, Bartley DJ, Martin RJ, Prichard RK, Gilleard JS: Anthelmintic resistance: markers for resistance, or susceptibility?. Parasitology. 2011, 138: 160-174. 10.1017/S0031182010001198.PubMed CentralView ArticlePubMedGoogle Scholar
- Njue AI, Prichard RK: Genetic variability of glutamate-gated chloride channel genes in ivermectin-susceptible and -resistant strains of Cooperia oncophora. Parasitology. 2004, 129: 741-751. 10.1017/S0031182004006183.View ArticlePubMedGoogle Scholar
- McCavera S, Rogers AT, Yates DM, Woods DJ, Wolstenholme AJ: An ivermectin-sensitive glutamate-gated chloride channel from the parasitic nematode Haemonchus contortus. Mol Pharmacol. 2009, 75: 1347-1355. 10.1124/mol.108.053363.PubMed CentralView ArticlePubMedGoogle Scholar
- Ardelli BF, Prichard RK: Identification of variant ABC-transporter genes among onchocerca volvulus collected from ivermectin-treated and untreated patients in Ghana, West Africa. Ann Trop Med Parasitol. 2004, 98: 371-384. 10.1179/000349804225003415.View ArticlePubMedGoogle Scholar
- James CE, Davey MW: Increased expression of ABC transport proteins is associated with ivermectin resistance in the model nematode Caenorhabditis elegans. Int J Parasitol. 2009, 39: 213-220. 10.1016/j.ijpara.2008.06.009.View ArticlePubMedGoogle Scholar
- Chen X, Yuan L, Du Y, Zhang Y, Wang J: Cross-resistance and biochemical mechanisms of abamectin resistance in the western flower thrips, Frankliniella occidentalis. Pesticide Biochemistry and Physiology. 2011, 101: 34-38. 10.1016/j.pestbp.2011.07.001.View 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.View ArticlePubMedGoogle Scholar
- Tribble ND, Burka JF, Kibenge FSB: Identification of the genes encoding for putative gamma aminobutyric acid (GABA) and glutamate-gated chloride channel (GluCl) alpha receptor subunits in sea lice (Lepeophtheirus salmonis). J Vet Pharmacol Ther. 2007, 30: 163-167. 10.1111/j.1365-2885.2007.00823.x.View ArticlePubMedGoogle Scholar
- Theocharidis A, van Dongen S, Enright AJ, Freeman TC: Network visualization and analysis of gene expression data using BioLayout express (3D). Nat Protoc. 2009, 4: 1535-1550. 10.1038/nprot.2009.177.View ArticlePubMedGoogle Scholar
- Dent JA, Smith MM, Vassilatis DK, Avery L: The genetics of ivermectin resistance in Caenorhabditis elegans. Proc Natl Acad Sci USA. 2000, 97: 2674-2679. 10.1073/pnas.97.6.2674.PubMed CentralView ArticlePubMedGoogle Scholar
- Kane NS, Hirschberg B, Qian S, Hunt D, Thomas B, Brochu R, Ludmerer SW, Zheng YC, Smith M, Arena JP, Cohen CJ, Schmatz D, Warmke J, Cully DF: Drug-resistant Drosophila indicate glutamate-gated chloride channels are targets for the antiparasitics nodulisporic acid and ivermectin. Proc Natl Acad Sci USA. 2000, 97: 13949-13954. 10.1073/pnas.240464697.PubMed CentralView ArticlePubMedGoogle Scholar
- Bloomquist JR: Chloride channels as tools for developing selective insecticides. Arch Insect Biochem Physiol. 2003, 54: 145-156. 10.1002/arch.10112.View ArticlePubMedGoogle Scholar
- Feng XP, Hayashi J, Beech RN, Prichard RK: Study of the nematode putative GABA type-A receptor subunits: evidence for modulation by ivermectin. J Neurochem. 2002, 83: 870-878. 10.1046/j.1471-4159.2002.01199.x.View ArticlePubMedGoogle Scholar
- Krause RM, Buisson B, Bertrand S, Corringer PJ, Galzi JL, Changeux JP, Bertrand D: Ivermectin: A positive allosteric effector of the α7 neuronal nicotinic acetylcholine receptor. Mol Pharmacol. 1998, 53: 283-294.PubMedGoogle Scholar
- Vogel C, Marcotte EM: Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat Rev Genet. 2012, 13: 227-232.PubMed CentralPubMedGoogle Scholar
- Cully DF, Vassilatis DK, Liu KK, Paress PS, VanderPloeg LHT, Schaeffer JM, Arena JP: Cloning of an avermectin sensitive glutamate-gated chloride channel from Caenorhabditis elegans. Nature. 1994, 371: 707-711. 10.1038/371707a0.View ArticlePubMedGoogle Scholar
- Cully DF, Paress PS, Liu KK, Schaeffer JM, Arena JP: Identification of a drosophila melanogaster glutamate-gated chloride channel sensitive to the antiparasitic agent avermectin. J Biol Chem. 1996, 271: 20187-20191. 10.1074/jbc.271.33.20187.View ArticlePubMedGoogle Scholar
- Shan Q, Haddrill JL, Lynch JW: Ivermectin, an unconventional agonist of the glycine receptor chloride channel. J Biol Chem. 2001, 276: 12556-12564. 10.1074/jbc.M011264200.View ArticlePubMedGoogle Scholar
- Silberberg SD, Li M, Swartz KJ: Ivermectin interaction with transmembrane helices reveals widespread rearrangements during opening of P2X receptor channels. Neuron. 2007, 54: 263-274. 10.1016/j.neuron.2007.03.020.View ArticlePubMedGoogle Scholar
- Yusein S, Velikova N, Kupenova P, Hardie R, Wolstenholme A, Semenov E: Altered ivermectin pharmacology and defective visual system in Drosophila mutants for histamine receptor HCLB. Invert Neurosci. 2008, 8: 211-222. 10.1007/s10158-008-0078-2.View ArticlePubMedGoogle Scholar
- Rao VTS, Siddiqui SZ, Prichard RK, Forrester SG: A dopamine-gated ion channel (HcGGR3*) from Haemonchus contortus is expressed in the cervical papillae and is associated with macrocyclic lactone resistance. Mol Biochem Parasitol. 2009, 166: 54-61. 10.1016/j.molbiopara.2009.02.011.View ArticlePubMedGoogle Scholar
- Burridge L, Weis JS, Cabello F, Pizarro J, Bostick K: Chemical use in salmon aquaculture: A review of current practices and possible environmental effects. Aquaculture. 2010, 306: 7-23. 10.1016/j.aquaculture.2010.05.020.View ArticleGoogle Scholar
- Labbé P, Berticat C, Berthomieu A, Unal S, Bernard C, Weill M, Lenormand T: Forty years of erratic insecticide resistance evolution in the mosquito Culex pipiens. PLoS Genet. 2007, 3: e205-10.1371/journal.pgen.0030205.PubMed CentralView ArticlePubMedGoogle Scholar
- Nebert DW, Gonzalez FJ: P450 genes: structure, evolution, and regulation. Annu Rev Biochem. 1987, 56: 945-993. 10.1146/annurev.bi.56.070187.004501.View ArticlePubMedGoogle Scholar
- Bernhardt R: Cytochrome P450: Structure, function, and generation of reactive oxygen species. Rev Physiol Biochem Pharmacol. 1995, 127: 137-221.Google Scholar
- Heckel DG: Insecticide resistance after silent spring. Science. 2012, 337: 1612-1614. 10.1126/science.1226994.View ArticlePubMedGoogle Scholar
- Nelson DR: Metazoan cytochrome P450 evolution. Comp Biochem Physiol C Pharmacol Toxicol Endocrinol. 1998, 121: 15-22. 10.1016/S0742-8413(98)10027-0.View ArticlePubMedGoogle Scholar
- Helvig C, Koener JF, Unnithan GC, Feyereisen R: CYP15A1, the cytochrome P450 that catalyzes epoxidation of methyl farnesoate to juvenile hormone III in cockroach corpora allata. Proc Natl Acad Sci USA. 2004, 101: 4024-4029. 10.1073/pnas.0306980101.PubMed CentralView ArticlePubMedGoogle Scholar
- Guittard E, Blais C, Maria A, Parvy JP, Pasricha S, Lumb C, Lafont R, Daborn PJ, Dauphin-Villemant C: CYP18A1, a key enzyme of Drosophila steroid hormone inactivation, is essential for metamorphosis. Dev Biol. 2011, 349: 35-45. 10.1016/j.ydbio.2010.09.023.View ArticlePubMedGoogle Scholar
- Baldwin WS, Marko PB, Nelson DR: The cytochrome P450 (CYP) gene superfamily in Daphnia pulex. BMC Genomics. 2009, 10: 169-10.1186/1471-2164-10-169.PubMed CentralView ArticlePubMedGoogle Scholar
- Dam E, Rewitz KF, Styrishave B, Andersen O: Cytochrome P450 expression is moult stage specific and regulated by ecdysteroids and xenobiotics in the crab Carcinus maenas. Biochem Biophys Res Commun. 2008, 377: 1135-1140. 10.1016/j.bbrc.2008.10.125.View ArticlePubMedGoogle Scholar
- Jones PM, O'Mara ML, George AM: ABC transporters: a riddle wrapped in a mystery inside an enigma. Trends Biochem Sci. 2009, 34: 520-531. 10.1016/j.tibs.2009.06.004.View ArticlePubMedGoogle Scholar
- Leslie EM, Deeley RG, Cole SPC: Multidrug resistance proteins: role of P-glycoprotein, MRP1, MRP2, and BCRP (ABCG2) in tissue defense. Toxicol Appl Pharmacol. 2005, 204: 216-237. 10.1016/j.taap.2004.10.012.View ArticlePubMedGoogle Scholar
- Bain LJ, LeBlanc GA: Interaction of structurally diverse pesticides with the human MDR1 gene product P-glycoprotein. Toxicol Appl Pharmacol. 1996, 141: 288-298.View ArticlePubMedGoogle Scholar
- Schinkel AH, Smit JJM, Van Tellingen O, Beijnen JH, Wagenaar E, Van Deemter L, Mol CAAM, Van Der Valk MA, Robanus-Maandag EC, Te Riele HPJ, Berns AJM, Borst P: Disruption of the mouse mdr1a P-glycoprotein gene leads to a deficiency in the blood–brain barrier and to increased sensitivity to drugs. Cell. 1994, 77: 491-502. 10.1016/0092-8674(94)90212-7.View ArticlePubMedGoogle Scholar
- Lanning CL, Ayad HM, Abou-Donia MB: P-glycoprotein involvement in cuticular penetration of [14C]thiodicarb in resistant tobacco budworms. Toxicol Lett. 1996, 85: 127-133. 10.1016/0378-4274(96)03654-5.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle Scholar
- Ardelli BF, Guerriero SB, Prichard RK: Ivermectin imposes selection pressure on P-glycoprotein from Onchocerca volvulus: Linkage disequilibrium and genotype diversity. Parasitology. 2006, 132: 375-386.View ArticlePubMedGoogle 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
- Helgesen KO, Horsberg TE: Influence of different materials on the concentration of delousing agents in sea water during bioassays. J Fish Dis. 2013, 36: 529-532. 10.1111/jfd.12046.View ArticlePubMedGoogle Scholar
- Chomczynski P, Mackey K: Modification of the TRI-reagent procedure for isolation of RNA from polysaccharide-rich and proteoglycan-rich sources. Biotechniques. 1995, 19: 942-945.PubMedGoogle Scholar
- Yasuike M, Leong J, Jantzen SG, von Schalburg KR, Nilsen F, Jones SRM, Koop BF: Genomic resources for sea lice: analysis of ESTs and mitochondrial genomes. Marine Biotechnol. 2012, 14: 155-166. 10.1007/s10126-011-9398-z.View ArticleGoogle Scholar
- Morais S, Pratoomyot J, Taggart JB, Bron JE, Guy DR, Bell JG, Tocher DR: Genotype-specific responses in Atlantic salmon (Salmo salar) subject to dietary fish oil replacement by vegetable oil: a liver transcriptomic analysis. BMC Genomics. 2011, 12: 255-10.1186/1471-2164-12-255.PubMed CentralView ArticlePubMedGoogle Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP: Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper - excel-based tool using pair-wise correlations. Biotechnol Lett. 2004, 26: 509-515.View ArticlePubMedGoogle Scholar
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.PubMed CentralView ArticlePubMedGoogle Scholar
- Leaver MJ, Villeneuve LAN, Obach A, Jensen L, Bron JE, Tocher DR, Taggart JB: Functional genomics reveals increases in cholesterol biosynthetic genes and highly unsaturated fatty acid biosynthesis after dietary substitution of fish oil with vegetable oils in Atlantic salmon (Salmo salar). BMC Genomics. 2008, 9: 299-10.1186/1471-2164-9-299.PubMed CentralView ArticlePubMedGoogle Scholar
- Berriz GF, Beaver JE, Cenik C, Tasan M, Roth FP: Next generation software for functional trend analysis. Bioinformatics. 2009, 25: 3043-3044. 10.1093/bioinformatics/btp498.PubMed CentralView ArticlePubMedGoogle 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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.