- Research article
- Open Access
Effects of the vertically transmitted microsporidian Facilispora margolisi and the parasiticide emamectin benzoate on salmon lice (Lepeophtheirus salmonis)
BMC Genomicsvolume 18, Article number: 630 (2017)
Microsporidia are highly specialized, parasitic fungi that infect a wide range of eukaryotic hosts from all major taxa. Infections cause a variety of damaging effects on host physiology from increased stress to death. The microsporidian Facilispora margolisi infects the Pacific salmon louse (Lepeophtheirus salmonis oncorhynchi), an economically and ecologically important ectoparasitic copepod that can impact wild and cultured salmonids.
Vertical transmission of F. margolisi was demonstrated by using PCR and in situ hybridization to identify and localize microsporidia in female L. salmonis and their offspring. Spores and developmental structures of F. margolisi were identified in 77% of F1 generation copepods derived from infected females while offspring from uninfected females all tested negative for the microsporidia. The transcriptomic response of the salmon louse to F. margolisi was profiled at both the copepodid larval stage and the pre-adult stage using microarray technology. Infected copepodids differentially expressed 577 transcripts related to stress, ATP generation and structural components of muscle and cuticle. The infection also impacted the response of the copepodid to the parasiticide emamectin benzoate (EMB) at a low dose of 1.0 ppb for 24 h. A set of 48 transcripts putatively involved in feeding and host immunomodulation were up to 8-fold underexpressed in the F. margolisi infected copepodids treated with EMB compared with controls or either stressor alone. Additionally, these infected lice treated with EMB also overexpressed 101 transcripts involved in stress resistance and signalling compared to the other groups. In contrast, infected pre-adult lice did not display a stress response, suggesting a decrease in microsporidian virulence associated with lice maturity. Furthermore, copepodid infectivity and moulting was not affected by the microsporidian infection.
This study demonstrated that F. margolisi is transmitted vertically between salmon louse generations and that biological impacts of infection differ depending on the stage of the copepod host. The infection caused significant perturbations of larval transcriptomes and therefore must be considered in future studies in which impacts to host development and environmental factors are assessed. Fitness impacts are probably minor, although the interaction between pesticide exposure and microsporidian infection merits further study.
Microsporidia are a diverse group of obligate intracellular, spore-forming parasites of invertebrates and vertebrates with over 150 genera and 1200 recognized species . Molecular phylogenetic evidence indicates that microsporidia are highly specialized fungi that parasitize a wide range of hosts . Many microsporidia have damaging impacts in agriculture, apiculture, and aquaculture [3,4,5], as well as in human health  particularly among the immunocompromised . Microsporidia are also of interest due to their use as biological control agents against insect pests [7, 8]. To date, microsporidia infecting humans and important insects have received the most research attention despite nearly half of all known microsporidia having aquatic hosts . For example, microsporidian infections of crustaceans, a diverse subphylum with approximately 67,000 species [10, 11], are only beginning to be uncovered. Given the ubiquity of crustaceans in the aquatic environment, the lack of information on interacting stressors within these organisms, and the economic importance of crustacean culture (i.e. 6.9 M tonnes worth $36 billion USD; ), studies examining the consequences of microsporidian infections are needed.
The microsporidian spore is the only stage capable of surviving outside of a host cell  and is the infective stage. Spores contain sporoplasm (i.e. infectious cytoplasm) and a polar filament that erupts under pressure and penetrates the host cell, providing a route through which the sporoplasm and nucleolar material can enter. Merogonic development of the parasite enlarges the host cell and this is followed by sporogonic development of the parasite. Spores can be released to infect adjacent cells (i.e. autoinfection), or to infect other individuals (e.g. through urine, faeces, decomposition) . Passage of infective spores among individual hosts in water or food may be the most common mode of horizontal transmission. Spores can also transmit vertically through eggs to infect offspring. Given the intimate association of microsporidia within host cells, it is not surprising that these parasites can have major effects on their hosts . Drastic host transcriptomic reductions in various functions have been observed in vivo . Pathological and physiological consequences of microsporidia infections have been characterized for a few terrestrial species, but those of infections in the marine environment remain poorly understood.
Microsporidia have evolved reduced genomes and other biological components (e.g. the known genome sizes are in the range of bacteria, 2.3–19.5 Mb) . The genome of Encephalitozoon cuniculi is well characterised and consists of ~2.9 Mb across 11 chromosomes, with approximately 2000 potential protein coding genes . Most studied microsporidia have a conserved set of microsporidia-specific genes, suggesting that genome reduction may have occurred prior to the diversification of the lineage [16, 17], but lineage-specific gene expansion and novelty is expected. The molecular basis and mechanisms of host manipulation are still being uncovered (recently reviewed in ).
Salmon lice are ectoparasitic copepods belonging to the family Caligidae (suborder Siphonostomatoida) that infect both wild and farmed salmonids, causing losses of more than $500 M USD/annum globally to the aquaculture industry . The salmon louse (Lepeophtheirus salmonis) is the most well studied species and has a direct life cycle with two free-living naupliar stages, an infective copepodid stage, two sessile parasitic chalimus stages, and three motile parasitic preadult and adult stages during which the lice are larger and more damaging to the host . Genetically distinct subspecies of L. salmonis occur in the Atlantic and Pacific Oceans [21, 22]. Infections on farmed Atlantic salmon (Salmo salar) are controlled by treatment with in-feed emamectin benzoate (EMB), although lice in the Atlantic Ocean have become resistant to this drug and to many others (reviewed by ). Alternative, non-chemical management options are required, and given the relevance of microsporidia to pest control combined with recent descriptions of microsporidia infections in sea lice, this area requires investigation.
Infections with two microsporidian species have been described in sea lice. Desmozoon lepeophtherii was originally identified in L. salmonis infecting farmed S. salar in Scotland and Norway [24,25,26] and a genetic variant of D. lepeophtherii was described from L. salmonis on farmed S. salar in the Pacific Ocean [27, 28]. The second microsporidian, Facilispora margolisi, was identified in L. salmonis infecting S. salar and Pacific salmon (Oncorhynchus spp.), and in two other Lepeophtheirus species from salmonid and non-salmonid fishes in the northeast Pacific Ocean . The prevalence of D. lepeophtherii in Pacific L. salmonis ranges from 5 to 15% whereas that of F. margolisi ranges from 50 to 90% .
Histological evidence of F. margolisi spores within ovarian tissue  suggests the possibility of vertical transmission. Furthermore, nothing is known regarding the effects of the infection on the survival, infectivity and fecundity of L. salmonis. Here we provide further evidence of vertical transmission of F. margolisi and determine whether the microsporidian affects copepodid survival and infectivity. We then characterize the impacts of the microsporidian infection on the copepod transcriptome, compare this response to known stress genes for L. salmonis, and determine whether the response is affected by the addition of a stressor (low dose EMB). Finally, we use this dataset and identify probes that are probably of microsporidian origin on the commonly used L. salmonis microarray. We then use these probes to diagnose infection status in an existing published transcriptome dataset of pre-adult L. salmonis, validate the findings using RT-qPCR, and compare the responses of pre-adult and copepodid lice to the microsporidian infection.
Vertical transmission of F. margolisi
Four experiments (Exp.) were conducted to evaluate different aspects of the biology of F. margolisi, including vertical transmission in L. salmonis (Table 1). Exp. 1 was used to evaluate vertical transmission of the microsporidian. Exp. 2 was also for this purpose, but included an experimental infection to test the effect on infectivity and development. Exp. 3 was used for in situ hybridization in order to determine locations of F. margolisi within the louse. Exp. 4 was used to perform transcriptome profiling of the interaction of F. margolisi infection and EMB exposure. Further details regarding each experiment are provided in Table 1. F. margolisi was detected in the cephalothorax and egg string samples of approximately 60% of the screened F0 individuals within these four experiments (n = 124). The microsporidian was detected in both the cephalothorax and matched egg string from 96% (n = 47) of F margolisi positive (MS+) copepods in which both tissues were tested. In the first experiment (Exp. 1), when individual egg strings were used to cultivate F1 individuals, 92% of the MS+ egg strings gave MS+ copepodid pools and none of the MS- egg strings produced MS+ offspring.
At 15 days post infection, Chum Salmon (Oncorhynchus keta) were infected with L. salmonis copepodids of known MS infection status (Exp. 2, Table 1) and there was no significant difference in the intensity of infection for MS+ or MS- lice (MS+: 6.5 ± 0.56; MS-: 7.0 ± 0.52 lice / fish; Chi2-test, p = 0.53). Most copepods in both infection groups (97.4% and 95.2%, respectively) had developed to the late chalimus II stage at the time of sampling. F. margolisi was detected in 72% of the chalimus II lice derived from MS+ egg strings but not in any from MS- egg strings (see Exp. 2 in Table 1).
Histological identification of F. margolisi in L. salmonis
In-situ hybridization (ISH) identified F. margolisi DNA in histological sections from adult female L. salmonis shown previously to be PCR-positive for F. margolisi (see Exp. 3). In the adult louse, ISH-positive spores and non-sporogonic structures were associated with cells immediately below the cuticle, in striated muscle, glandular tissues, ovary and egg strings (Fig. 1). ISH-positive reactions were also visualised within unidentified tissues of larval copepods derived from PCR-positive females (Fig. 1). No ISH reactions were observed in PCR-positive samples incubated with unlabelled probe or in PCR-negative samples incubated with the labelled probe.
Transcriptomic response of copepodid salmon lice to F. margolisi
Copepodid pools (F1; see Exp. 4), each originating from a single adult female (MS+ or MS-), were incubated either in seawater alone (unexposed control) or in seawater containing 1 ppb emamectin benzoate (EMB) to evaluate the effect of MS infection and the effect of combined MS and EMB on the L. salmonis transcriptome (7 pools per condition, 28 pools in total; Fig. 2). Principal components analysis (PCA) of all quality filtered transcripts indicated a larger effect of MS than of EMB (Fig. 3). The first axis differentiated the MS+ and MS- pools and explained 29.5% of the total variation. Clustering associated with EMB exposure was less obvious. Treatment groupings most distant from one another in the PCA were MS−/EMB- and MS+/EMB+. The number of differentially expressed transcripts also supported a greater effect of MS than 1 ppb EMB (Fig. 4). MS infection affected 577 transcripts whereas EMB affected 228 transcripts (main effects with no interaction effect; p ≤ 0.01 and fold change (FC) ≥ 1.5; Additional file 1). The range of FC was also greater for transcripts with a main effect of MS (FC = −9.7 to 4.5) compared with those showing a main effect of EMB (FC = −2.8 to 2.4; Additional file 1). Similarities between the effects MS and EMB were indicated from 150 transcripts concordantly differentially expressed by EMB and MS (Fig. 4; Additional file 1).
To better characterize response genes, transcripts were classified into the following categories: i) MS-specific response genes; ii) general stress response genes, responding concordantly to MS and EMB; and iii) interacting stressor genes, whose response depended on the presence of the second stressor (see Fig. 4). As the majority of transcripts affected by EMB exposure also responded to MS infection (Fig. 4), an EMB-specific response category was not included. All transcripts with an EMB-specific response can be found in Additional file 1.
MS-specific response genes
Of the 577 transcripts responding to F. margolisi infection (main effect, no MS/EMB interaction effect), 427 (74%) were MS-specific (Fig. 4). Of these, 212 were overexpressed and 215 were underexpressed in MS+ L. salmonis.
Overexpressed MS-specific transcripts included 12 cuticle-associated transcripts (e.g. cuticle protein 6, cuticle protein 7, cuticle protein C14.6) and 18 transcripts associated with muscle (e.g. myosin heavy chain muscle, myosin-3, tropomyosin, tropomyosin-2, tropomyosin alpha-1 chain, troponin I and troponin T; Additional file 1). Cytoskeleton-related transcripts were overexpressed and enriched (GO: cytoskeleton, 14 transcripts; p = 0.003) in MS+ lice including restin homolog, actin clone 403, actin-related protein 2/3 complex subunit 4, actin cytoplasmic 1, tubulin alpha-1A chain, tubulin alpha-3 chain, and tubulin beta-1A chain.
Underexpressed MS-specific transcripts were related to the mitochondria (e.g. GO: mitochondrial part, 16 transcripts; p < 0.01; Additional file 2) and ribosomes. Mitochondria-related transcripts included NADH dehydrogenase 1 alpha subcomplex subunits 2, 10 and 12, cytochrome c oxidase subunit 1, 5B and 6B, CYP450 2 J2 and CYP450 2 L1 (Additional file 1). Ribosomal proteins (rp) underexpressed in MS+ lice included those present in the mitochondria (e.g. 28S rp s32, 39S rp L22, 39S rp L23, 39S rp L32, 40S rp s12 and 60S rp L12) and those in the cytoplasm (e.g. 40S rp s16, 60S rp L3, 60S rp L29, 60S rp L30, 60S rp L44, rp s6 kinase beta 2), suggesting an overall reduction in ribosome constituents in the MS+ lice.
General stress response genes
A total of 150 transcripts were overexpressed in both the MS+ and EMB+ treatment groups (main effects, no interaction effect). These were enriched for protein folding and ATP binding functions (Table 2; Additional file 2). Overexpressed transcripts included several previously associated with stress in L. salmonis such as 60 kDa heat shock protein mitochondrial, DnaJ homolog subfamily A member 1, heat shock 70 kDa protein cognate 4, heat shock protein 90-alpha, heat shock cognate protein 90-beta and T-complex proteins (TCP) TCP-1-alpha, TCP-1-beta, TCP-1-epsilon, TCP-1-theta, and TCP-1-zeta (Additional file 1; ). Fold changes of these transcripts were more modest than those of MS-specific response genes, and were typically in the range of 1.5- to 2.1-fold (Additional file 1).
Underexpressed general stress response transcripts included seven proteases (e.g. trypsin-1, neprilysin-11, carboxypeptidase B, and hypodermin-B; Table 2; Additional file 1). All underexpressed proteases were annotated as secreted by SP_PIR_Keyword, and this list was enriched for peptidase activity acting on L-amino acid peptides (Table 2). Furthermore, multiple solute carrier (SLC) family members were underexpressed including high affinity copper uptake protein 1 (slc31A1), solute carrier family 15 member 1 (slc15A1) and sodium-dependent phosphate transport protein 2B (slc34A2). These transcripts, and others, enriched the GO category substrate-specific transmembrane transporter activity (Table 2). Transporters were underexpressed in MS+ lice and in EMB+ lice, and five of these were underexpressed in both groups.
Interacting stressor response genes
The expression of 290 transcripts showed a significant interaction between EMB exposure and MS (p < 0.01 and FC ≥ 1.5; Fig. 4). The expression similarities among these transcripts were characterized by k-means clustering to separate them into six general patterns (clusters A-F; Fig. 5). Transcripts in clusters A, C, and E were overexpressed and B, D, and F were under-expressed in at least one of the treatment groups relative to the controls.
The majority of the transcripts with significant interaction effects belong to two classes. Genes in cluster A (32/290 transcripts; Fig. 5) had equal overexpression in both single stressors (MS or EMB) but were not overexpressed additively in the double stressor condition (MS+/EMB+). These transcripts were similar to the general stress response genes described above, including calreticulin, Dnaj homolog subfamily B member 11, and stress-induced-phosphoprotein 1 (Additional file 3). The second main class included transcripts in clusters E (101/290) and F (49/290), that were only affected in the double stressor group (EMB+/MS+) (Fig. 5). Cluster E included overexpressed catalase, DnaJ homolog subfamily A member 1, heat shock protein 81–1, or structural transcripts such as myosin heavy chain, myosin light chain alkili, muscle M-line assembly protein unc-89, tropomyosin, and many others (Additional file 3). Cluster F included transcripts only underexpressed in MS+/EMB+ lice such as the transporters slc25A36, protein spinster, and RhGB. RhGB was previously shown to be down-regulated in L. salmonis by EMB, cypermethrin, and hyposalinity .
Cluster B had the largest fold changes, and contained 48 transcripts 1.8- to 8-fold underexpressed in the double stressor condition (Fig. 5). Some of these transcripts were also slightly downregulated (FC < 2) in MS+ lice. Cluster B largely contained serine-type endopeptidases including those annotated as trypsin-1 (9 different contigs), anionic trypsin-1, chymotrypsin A chain C, chymotrypsin BI (2 different contigs), hypodermin-B, ovochymase-1, and collagenase (Fig. 6; Additional file 3). Also included were other degradative enzymes such as cathepsin-L light chain, cathepsin-D, and carboxylesterase (2 different contigs), and detoxification-related transcripts such as cytochrome p450 2 J2 and carboxypeptidase B (4 different contigs). Of the 34 annotated transcripts in Cluster B, 22 were annotated with the SP_PIR_Keyword secreted, 16 of which contained at least one trypsin-like protease domain (CDD: smart00020, NCBI).
Transcripts within cluster B were further compared with other published L. salmonis microarray datasets due to their putative secreted degradative enzyme function. Of the 48 unique transcripts (based on unique contig identifiers), 36 were overexpressed in actively feeding lice relative to starved lice (Fig. 6; ), 20 transcripts were overexpressed in EMB-resistant relative to EMB-sensitive Atlantic lice, and 32 transcripts were down-regulated by 50 ppb EMB exposure in pre-adult Pacific lice . Furthermore, 20 transcripts were affected in Pacific copepodids by hyposaline conditions (Fig. 6 and Additional file 4) . These genes may therefore be important for a range of processes, including feeding, stress response and drug resistance.
Microsporidian genes on the salmon louse microarray
Archived individual transcriptome profiles of pre-adult Pacific salmon lice exposed to very low concentrations of EMB that did not affect the transcriptome (GEO accession: GSE73734; ) were tested for the presence of F. margolisi using PCR probes (see Methods). Of the 40 samples, 22 (55%) were MS+ and 18 were MS-, permitting the analysis of genes affected by MS within these pre-adult individuals.
In the pre-adult transcriptomes, few differences were observed between MS+ and MS- individuals (Additional file 5). However, the expression of 20 transcripts was detected in all 22 MS+ individuals, but not in any of the 18 MS- samples (i.e. 100% of MS- samples had a Cy5 value below the QC filter background). Interestingly, these transcripts did not pass background fluorescence thresholds in any of the copepodid samples from the previously described MS/EMB study, regardless of MS infection status. Of the 20 transcripts expressed exclusively in MS+ pre-adults, 14 (70%) were annotated to genes from other microsporidian parasites such as Spraguea lophii (4 transcripts) and Nosema ceranae (3 transcripts; Table 3). Two of the remaining six transcripts were annotated to the slime mold Dictyostelium discoideum and the causative agent of malaria Plasmodium falciparum, and four had no annotation. The presence and absence of expression in the MS+ and MS- samples identified on the microarray was validated using RT-qPCR for subunit alpha of phenylalanyl-tRNA synthetase, heat shock protein 90, and 40S ribosomal protein S4 (Table 3; Additional file 6).
The exclusive presence of these transcripts in MS+ individuals, combined with their annotation to genes from other microsporidian species suggests that they are not from L. salmonis but rather of microsporidia origin in the original samples used to create the microarray. The concordance in profiles between the diagnostic primers , microarray probes, and RT-qPCR primers suggests that these sequences may be useable to screen for the presence of F. margolisi in individual adult L. salmonis samples, but not in copepodid pools. These probes will be flagged as microsporidian origin in future updates of the 38 k microarray annotation file.
Microsporidia can have large effects on their hosts and have potential as biological control agents in pest management. In this study, we demonstrate that Facilispora margolisi, a parasite of L. salmonis in the Eastern Pacific Ocean , is vertically transmitted and has stage-specific impacts on host energetic and structural gene expression. Infection was associated with transcriptomic evidence of a stress response in copepodids, which was absent in the pre-adults. However, infection with F. margolisi did not affect the number of L. salmonis copepodids able to attach and moult on Chum Salmon. Combinations of drug and MS infection impacted transcripts previously identified as related to feeding and drug response and resistance. These results are discussed below in terms of the evolution and biology of key traits in salmon lice.
Microsporidia transmission and stage-specific effects
Vertical transmission of microsporidia is common in other crustacean hosts and can be the sole means of transmission  or as part of a mixed transmission strategy . Vertical transmission of F. margolisi was identified using PCR and in situ hybridization to locate the parasite within the ovary and developing embryos of infected L. salmonis females and their offspring. Spore dimorphism in F. margolisi, specifically the presence of larger spores with longer polar filaments, indicates the possibility of horizontal transmission , but this was not yet confirmed experimentally. Vertically transmitted microsporidia generally are less virulent than horizontally transmitted species as they depend on host survival for replication . Nonetheless, vertically transmitted MS can affect host growth , and induce male feminization , thereby impacting sex ratios .
Despite the absence of measurable effects on L. salmonis infectivity or development, infection with F. margolisi impacted host transcriptomes particularly in larval stages. MS+ copepodids differentially expressed 577 transcripts that were enriched for energetic and stress functions whereas in adults only 123 transcripts were differentially expressed. Of these transcripts, only five were differentially expressed in both developmental stages, including a mitochondrial chaperone DnaJ homolog subfamily A member 1, calreticulin, and three transcripts without annotation (Additional files 1 and 5). The reason for the reduced host impact at later lice life stages is not known, but it may be related to the reduced virulence required for vertical transmission. Furthermore, copepodids are generally more sensitive to stress than adults. For example, pre-adult L. salmonis show a transcriptional response to 50 ppb EMB but not 25 or 10 ppb EMB  while copepodids respond at 1 ppb EMB, as observed here. Improved understanding of the sites of microsporidian infections during early copepod development may inform our understanding of the physiological impacts of these infections.
Lice response to F. margolisi
MS+ copepodids underexpressed mitochondrial component genes, including ribosomal subunits and those involved in mitochondrial organization. Close proximity between early stages of F. margolisi and mitochondria of L. salmonis has been observed . This is a common strategy among microsporidia that results from reduced mitosomes and an inability to generate ATP . Intricate manipulations of the host mitochondria structure, and function permit energy acquisition for microsporidian growth and development [39, 40]. In the host this can cause prolonged energetic stress by up-regulating important metabolic pathways [41,42,43]. Increased metabolic activity was suggested in this study as MS+ copepodids overexpressed genes involved in cellular respiration and ATP binding while transcripts related to mitochondrial structure were underexpressed. It is therefore possible that F. margolisi is associated with manipulation of mitochondrial structure and function in L. salmonis.
The overexpression of muscle- and cuticle-related transcripts in F. margolisi infected copepodids may also be a result of manipulation by the microsporidia. Mature spores and other developmental stages of F. margolisi were observed within striated muscle and subcuticular connective tissue of L. salmonis using in situ hybridization, supporting previous findings from electron microscopy . Although there was no indication of an innate immune response against F. margolisi in the lice, this may be due to the generally poor annotation of immunity genes in crustaceans . Microsporidian infections are known to impact a variety of immune-related genes in other arthropods [45, 46] and nematodes . Transcripts without annotation accounted for 47% of those affected by MS and may be useful to consider for future studies on biotic stress in L. salmonis.
Non-additive impacts on lice hosts from parasiticide and microsporidia
Our observed interactions between drug treatment and microsporidian infection may be important in understanding the development of drug resistance in L. salmonis. In this study, differences from controls in the transcriptome of MS+ L. salmonis exposed to EMB were greater than differences from either stressor alone, and a total of 290 transcripts showed differential expression in a non-additive manner. In the western honey bee Apis mellifera, additive impacts were observed between MS infection and chemical exposure that resulted in changes to fitness costs depending on the mechanism of drug resistance . Pesticide exposure may also affect microsporidian virulence and host susceptibility, as in Flabelliforma magnivora and Nosema ceranae infecting Daphnia magna  and A. mellifera , respectively.
The largest fold changes between MS+/EMB+ lice and controls were for transcripts in cluster B (Fig. 5), which were underexpressed in only the MS+/EMB+ lice. These specific transcripts are also involved in feeding , stress response , and EMB resistance . Braden et al.  found 36 of the 48 cluster B transcripts to be underexpressed after 24 h and 48 h of starvation (Additional file 4). Based on annotation, these transcripts may be involved in feeding and digestion, as well as host immunosuppression (Fig. 6) and included 16 trypsin-like serine proteases. Trypsin-like serine proteases are present in lice secretions [51,52,53] and many of these transcripts were most highly expressed in the gut relative to other tissues . Another transcript from this cluster, hypodermin B, inhibits the activation of complement component C3 in vitro  and may therefore play a role in immunomodulating salmonid hosts. Although F. margolisi infection or EMB exposure alone had minimal effects on these transcripts, in combination they resulted in underexpression (FC = 1.5–8.0). Future studies should address the localized effects of F. margolisi infection based on large gene expression changes associated with infected lice tissues (i.e. muscle, cuticle, mitochondria, and glandular tissue) (Fig. 1 and ) in order to determine possible interactive effects of microsporidia infection and louse feeding, digestion, mobility, moulting, energy expenditure, and overall fitness.
In addition, 32 of the 48 cluster B transcripts were also found to be down-regulated from an exposure to 50 ppb EMB in Pacific pre-adult lice, and overexpressed in an EMB resistant Atlantic L. salmonis population (Additional file 4; ). Resistance to EMB occurs in L. salmonis in the Northeast and Northwest Atlantic Ocean while lice in the Pacific subspecies remain sensitive [23, 56]. Although there is no evidence that the presence of F. margolisi contributes to the maintenance of EMB sensitivity in the Pacific lice, the correlation and the non-additive impacts of F. margolisi and EMB on genes involved in stress, feeding, and host attachment indicate the need for further research to better understand a role of microsporidian infections on drug tolerance, selection, and fitness.
Non-host material in next-generation sequencing data
We identified a set of 20 transcripts within the L. salmonis microarray that most likely belong to the microsporidian F. margolisi. The presence of non-host material in a sequence database highlights the importance of identifying and/or removing non-target eukaryotic sequences from assemblies, as wild individuals or lines with persisting, vertically transmitted intracellular eukaryotic parasites represent a mixture of organisms . These microsporidian-like transcripts were only detectable in cDNA generated from MS+ pre-adults, but not in that from MS+ copepodids, potentially indicating stage-specific expression of F. margolisi genes or a higher effective microsporidian load in pre-adult lice. Future RNA-sequencing studies can potentially use these genes to quantify microsporidian infections in L. salmonis as was done with microsporidia of humans .
Although F. margolisi induced a stress-like transcriptomic response in copepodids of L. salmonis, signatures of stress were absent in pre-adults and the infection did not cause changes in salmon lice development or infection potential towards fish hosts. Collectively our observations of limited impacts within the later developmental stages are consistent with vertical transmission of the microsporidian and its high prevalence within L. salmonis populations. However, the non-additive effects from combinations of microsporidia infection and parasiticide treatment merit further study. In addition to microsporidia, the salmon louse also hosts viruses  and bacteria  and the present observations indicate a need to better understand the possible influence of hyperparasitism with respect to copepod sensitivity to parasiticides as well as the confounding effects these infections may have on transcriptomic studies.
Ovigerous Lepeophtheirus salmonis were collected from wild adult Chum Salmon (Oncorhynchus keta) captured in a regulatory test fishery or from farmed Atlantic Salmon (Salmo salar) during harvest and maintained in aerated seawater at 10 °C for no longer than 48 h before processing.
Chum Salmon from the Nanaimo River Hatchery were reared from swim-up fry in a mixture of freshwater and seawater (9.2 °C – 10.0 °C) and provided a daily ration of commercial pellets.
Culture of L. salmonis larvae
Four experiments (Exp. 1–4) were conducted in which the F. margolisi infection status of individual adult female L. salmonis and their nauplius or copepodid offspring were determined and compared (Table 1). For each experiment between 20 and 40 ovigerous copepods bearing two intact and pigmented egg strings were selected. From each, one egg string was placed into a sterile flask containing 300 mL of aerated seawater and incubated at 10 °C for 6 to 8 days. The second egg string and the anterior third of the cephalothorax were preserved separately in 95% ethanol and the remainder of the specimen was preserved in neutral buffered 10% formalin (NBF). All dissections were conducted using tools rinsed sequentially with 4% sodium hypochlorite, water and 95% ethanol between specimens. The infection status of each specimen was determined from separate polymerase chain reactions (PCR, see below) of the cephalothorax and egg-string.
In Exp. 1, 20 larval cultures were separately filtered through 47 mm cellulose acetate/cellulose nitrate membranes with a pore size of 8.0 μm (EMD Millipore). The membranes were flash-frozen in liquid nitrogen and stored at −80 °C for subsequent PCR analysis. In Exp. 2, 24 larval cultures were pooled according to infection status of the source copepod. The numbers of copepodids and nauplii in each pool were determined by microscopic examination and used to infect naïve salmon (see below). In Exp. 3, 40 larval cultures were separately pooled according to infection status of the source copepod and fixed in NBF for 24 h followed by storage in 70% ethanol for in situ hybridisation assays (see below). In Exp. 4, 40 larval cultures were pooled according to the F. margolisi infection status of the source copepod. Each pool was divided into 16 sub-pools, 14 of which were maintained in aerated seawater for use in the emamectin benzoate (EMB) exposure study (see below) while the remaining two pools from each of infected and uninfected females were filtered and stored as above for PCR confirmation.
Polymerase chain reaction (PCR) detection
DNA was extracted from ethanol-preserved specimens using the DNeasy Animal Tissue protocol (Qiagen), as per manufacturer’s instructions. A region of F. margolisi SSU rDNA was amplified with PCR as described by . Samples were scored F. margolisi positive (MS+) or negative (MS-) using 1.5% agarose gel electrophoresis.
Exposure of salmon to L. salmonis
Juvenile Chum Salmon (11.0 ± 0.3 g; n = 6 per tank) were acclimated for 6 days in duplicate 30-L tanks containing aerated seawater at 30 ppt and 9 °C, flowing at 1 L min−1. Fish in each tank were exposed to 1570 and 1680 7-day post hatch (dph) copepodids derived from infected and uninfected female copepods (Exp. 2, above), respectively, using the method described earlier . At 15 dph, the salmon were euthanized by immersion in 200 mg L−1 tricaine methane sulphonate and all copepods were removed and stored in 95% ethanol. Each copepod was assessed microscopically for stage of development and for the presence of F. margolisi by using PCR.
NBF-preserved specimens from Exp. 3 were processed routinely for histology and subsequent microscopic examination as previously described . During processing the larvae were aggregated prior to embedding to facilitate detection. Sections of 5 μm were applied to uncoated or silane coated glass microscope slides and stained routinely (Gram or Giemsa stains), or processed for in-situ hybridisation, respectively.
In-situ hybridization probe design, synthesis and assay
An 84 bp region of the F. margolisi SSU rRNA gene was amplified (see primers in Additional file 6) to serve as template for subsequent digoxygenin (DIG) labelling reactions using a probe synthesis kit (Roche Applied Science). Each reaction contained 2 μl of genomic DNA, 1× PCR reaction buffer (Invitrogen), 1.5 mM MgCl2 (Invitrogen), 0.2 mM of each dNTP, 0.025 U/μl of Platinum® Taq DNA polymerase (Invitrogen), and sterilized water. Positive (F. margolisi DNA) and negative (water) controls were included alongside all reactions. The PCR profile consisted of 95 °C for 5 min, 35 cycles at 95 °C for 30 s, 55 °C for 30 s, and extension at 72 °C for 2 min, followed by a final extension at 72 °C for 10 min. Amplified PCR products were visualized on ethidium bromide stained, 1% agarose gels. PCR products were purified using QIAquick® PCR Purification Kit (Qiagen) and quantified with NanoDrop-1000 Spectrophotometer. DIG-labelling reactions were carried out according to the manufacturer’s instructions using 30 ng of purified PCR product with the same primers and thermal profile as the conventional PCR and DIG labelled products were purified as described above.
The ISH protocol was adapted from earlier work (Antonio et al. 1998, Jones et al. 2003). Tissue sections on silane-coated slides were routinely deparaffinized and permeabilized with 20 ng/μl of proteinase K (Qiagen) for 10 min. Sections were incubated with 50 μl of hybridization buffer (5.1 ml deionised formamide, 2.0 ml 20× standard saline citrate (SSC; 3 M NaCl, 0.3 M sodium citrate, pH 7.0), 2.0 ml 5% dextran sulphate, 0.5 ml denatured sperm DNA (10 mg/ml), 0.2 ml 50× Denhardt’s solution (1% acetylated BSA, 1% polyvinylpyrrolidone, and 1% Ficoll 400 in molecular biology grade water) and 25 μl 10% SDS per 10 ml of solution) for 1 h at 37 °C. A volume of 50 μl of digoxygenin (DIG) - labelled probe (0.1 ng/μl in hybridization buffer) was placed on each tissue section. Hybridization, formazan precipitation, and probe visualization were performed at 37 °C following an established protocol . Controls included sections from PCR-negative copepods stained with DIG - labelled probe and sections from PCR-positive copepods stained with unlabelled probe.
Exposure to EMB and extraction of RNA
Four experimental conditions were tested (Exp. 4) based on the presence of the microsporidian F. margolisi (MS+ or MS-) and EMB (EMB+ or EMB-). These conditions are labeled as MS+/EMB+, MS+/EMB-, MS−/EMB+, MS−/EMB-, and each with seven biological replicates. Each replicate was a flask containing between 50 and 75 copepodids. EMB (Sigma-Aldrich) was dissolved at a concentration of 1.0 ppb and aerated cultures were maintained for 24 h at 10 °C before filtration, flash-freezing and storage as above. The frozen filters were homogenized (mixer mill; Retsch® MM 301), and total RNA was extracted using TRIzol® (Invitrogen), as per manufacturers’ instructions. Total RNA was purified through RNeasy spin columns with an on-column DNase I treatment (Qiagen). Total RNA was quantified by spectrophotometry (NanoDrop-1000), and quality checked by electrophoresis on a 1% agarose gel. Samples were then randomized for all downstream nucleic acid manipulations.
Labeled cRNA was generated from total purified RNA using Low Input Quick Amp Labeling kits (Agilent) as per manufactures’ instructions and as reported previously . A Cy3-cRNA pool was generated for a reference design by synthesizing Cy3-cRNA from three randomly selected samples from each of the four experimental conditions and combining equimolar amounts from each into a common pool. Samples were hybridized as per manufacturer’s instructions as previously described  to a 38 K oligo microarray designed using previously annotated ESTs from both Pacific and Atlantic L. salmonis  using eArray (Agilent). Slides were scanned on a ScanArray® Express (Perkin Elmer) at 5 μm resolution using PMT settings optimized to have the median signal of ~1–2% of array spots saturated (Cy5: 70; Cy3: 75).
Images were quantified in Imagene 8.1 (Biodiscovery), poor spots flagged, and background corrected as reported previously . Sample files were loaded into GeneSpring 13.0 (Agilent). Raw sample files have been uploaded to GEO (GSE94692). Samples were normalized as follows: raw value threshold of 1.0; intensity-dependent Lowess normalization; and baseline transformation to the median of all samples. Control spots, and any probes not passing the following filter were removed from the analysis: raw values ≥500 in both Cy3 and Cy5 channels and no poor quality flags in at least 65% of samples in any one condition. A principle component analysis was performed within GeneSpring (Agilent) on the samples to investigate for any clustering of technical or biological variables.
Probes were tested using a two-way ANOVA (p < 0.01; FC ≥ 1.5) for main effects of microsporidia, main effects of EMB, and a significant interaction effect of microsporidia and EMB. Probes with a significant interaction effect were subtracted from the main effect lists to be considered with the interaction effect list. Probes with a significant interaction effect were separated by k-means clustering using the least number of clusters until cluster redundancy was visible (Euclidean distance metric; 6 clusters; 50 iterations; Agilent). Gene Ontology and pathway analysis was performed on significant entity lists using UniProt accession numbers in the DAVID bioinformatics tool [63,64,65] using a modified Fisher’s exact test (p < 0.05; genes/enrichment category ≥4) with a background list of all entities passing quality control filters. This background list was used for a principal components analysis that clustered samples according to expression of genes. Samples were colored or shaped differently in order to identify trends associated with experimental conditions.
Pre-adult L. salmonis exposed to EMB
A fifth experiment, which was previously described , included exposure of F1 generation pre-adult II L. salmonis (21 females and 19 males) to either seawater, 0.01 ppb EMB (Sigma-Aldrich), or 0.1 ppb EMB. These lice were collected from farmed Atlantic salmon in the Broughton Archipelago, British Columbia (BC). Using the previously extracted RNA, as described , cDNA for detection of the microsporidian and qPCR of selected genes in these archived samples was synthesized for individual lice from 1 μg of purified RNA using iScript™ Reverse Transcription Supermix kit (BioRad) as per manufacturer’s instructions (see ). Detection of F. margolisi was completed as described earlier  with the exception that cDNA derived from purified RNA was used as template and a touch-down PCR assay with decreasing annealing temperatures from 63 °C to 58 °C was completed before standard PCR for an additional 25 cycles to amplify the target sequence.
The microarray results were validated using RT-qPCR of five targets to confirm expression patterns. Primer efficiencies were confirmed by creating standard curves with 5–6 points and a 5-fold series dilution with a r2 > 0.95 and effeciencies between 0.90 and 1.05. For RT-qPCR amplification, SsoAdvanced SYBR Green Supermix (BioRad) was used in 11 μl reactions with 1 μl template and 0.1 μM of each primer using the following program: 95 °C for 2 min, followed by 40 cycles of 95 °C for 5 s and 60 °C for 15 s. Melt curve analysis was performed at the end of each qPCR assay by increasing the temperature in 0.5 °C increments (65 °C to 95 °C) every 5 s. No template controls (NTC) and no RT controls were negative for all genes assayed. Normalization for genes of interest was completed using qbase-PLUS (Biogazelle; Gent, Belgium) with an output of log2 ratios relative to the reference genes elongation factor 1α (ef1α) and ribosomal protein subunit 20 (rps20). Reference gene stability was assessed using geNorm , which showed an M value of 0.35 and a coefficient of variation (CV) of 0.13.
Normalized log2 expression measured by RT-qPCR was compared to Cy5/Cy3 log2 expression ratios produced by the microarray. For genes specifically expressed in MS+ lice on the microarray, this was validated by RT-qPCR. For genes expressed in all conditions, a Pearson’s correlation (R v. 3; ) used to validate the microarray (Additional file 6).
Analysis of variance
Days post hatch
- F1 :
In situ hybridization
Microsporidian F. margolisi
Neutral buffered formalin
Principal components analysis
Polymerase chain reaction
- SSU rRNA:
Small subunit ribosomal RNA
Keeling PJ, Fast NM. Microsporidia: biology and evolution of highly reduced intracellular parasites. AnnuRevMicrobiol. 2002;56:93–116.
Keeling P. Five questions about microsporidia. PLoS Pathog. 2009;5:e1000489.
Cali A, Takvorian PM. Developmental morphology and life cycles of the microsporidia. In: Weiss LM, Becenel JJ, editors. Microsporidia Pathog. Oppor. First ed. Ames: Wiley; 2014. p. 71–133.
Higes M, García-Palenciab P, Martín-Hernándeza R, Meanac A. Experimental infection of Apis mellifera honeybees with Nosema ceranae (Microsporidia). J Invertebr Pathol. 2007;94:211–7.
Nylund S, Andersen L, Sævareid I, Plarre H, Watanabe K, Arnesen CE, et al. Diseases of farmed Atlantic salmon Salmo salar associated with infections by the microsporidian Paranucleospora theridion. Dis Aquat Org. 2011;94:41–57.
Didier ES, Weiss LM. Microsporidiosis: current status. Curr Opin Infect Dis. 2006;19:485–92.
Hajek AE, Delalibera I. Fungal pathogens as classical biological control agents against arthropods. BioControl. 2010;55:147–58.
Lomer C, Bateman R, Johnson D, Langewald J, Thomas M. Biological control of locusts and grasshoppers. Annu Rev Entomol. 2001;46:667–702.
Stentiford GD, Feist SW, Stone DM, Bateman KS, Dunn AM. Microsporidia: diverse, dynamic, and emergent pathogens in aquatic systems. Trends Parasitol. 2013;29:567–78.
Ahyong ST, Lowry JK, Alonso M, Bamber RN, Boxshall GA, Castro P, et al. Animal biodiversity: An outline of higher-level classification and survey of taxonomic richness. Zootaxa. 2011;3148:165–91.
Martin JW, Davis GE. An updated classification of the recent crustacea. Natural History Museum of Los Angeles County. 2001.
FAO. Global aquaculture production statistics. Glob Crustac Cult. 2014; Available from: http://www.fao.org/fishery/statistics/en. Accessed 4 Nov 2016.
Vávra J, Ronny Larsson JI. Structure of microsporidia. In: Weiss LM, Becenel JJ, editors. Microsporidia Pathog. Oppor. First ed. Ames: Wiley; 2014. p. 1–70.
Watson AK, Williams TA, Williams BAP, Moore KA, Hirt RP, Embley TM. Transcriptomic profiling of host-parasite interactions in the microsporidian Trachipleistophora hominis. BMC Genomics. 2015;16:983.
Katinka MD, Duprat S, Cornillot E, Méténier G, Thomarat F, Prensier G, et al. Genome sequence and gene compaction of the eukaryote parasite Encephalitozoon cuniculi. Nature. 2001;414:450–3.
Nakjang S, Williams TA, Heinz E, Watson AK, Foster PG, Sendra KM, et al. Reduction and expansion in microsporidian genome evolution: new insights from comparative genomics. Genome Biol Evol. 2013;5:2285–303.
Heinz E, Williams TA, Nakjang S, Noël CJ, Swan DC, Goldberg AV, et al. The genome of the obligate intracellular parasite Trachipleistophora hominis: new insights into microsporidian genome dynamics and reductive evolution. PLoS Pathog. 2012;8(10): e1002979.
Dean P, Hirt RP, Embley TM. Microsporidia: Why Make Nucleotides if You Can Steal Them? PLoS Pathog. 2016;12:e1005870.
Costello MJ. The global economic cost of sea lice to the salmonid farming industry. J Fish Dis. 2009;32:115–8.
Hamre LA, Eichner C, Caipang CMA, Dalvin ST, Bron JE, Nilsen F, et al. The salmon louse Lepeophtheirus salmonis (Copepoda: Caligidae) life cycle has only two chalimus stages. PLoS One. 2013;8(9):e73539.
Yazawa R, Yasuike M, Leong J, von Schalburg KR, Cooper GA, Beetz-Sargent M, et al. EST and mitochondrial DNA sequences support a distinct Pacific form of salmon louse, Lepeophtheirus salmonis. Mar Biotechnol. 2008;10:741–9.
Skern-Mauritzen R, Torrissen O, Glover KA. Pacific and Atlantic Lepeophtheirus salmonis (Krøyer, 1838) are allopatric subspecies: Lepeophtheirus salmonis salmonis and L . salmonis oncorhynchi subspecies novo. BMC Genet. 2014;15:32–40.
Aaen SM, Helgesen KO, Bakke MJ, Kaur K, Horsberg TE. Drug resistance in sea lice: a threat to salmonid aquaculture. Trends Parasitol. 2015;31:72–81.
Freeman M. A hyperparasitic microsporidian infecting the salmon louse, Lepeophtheirus salmonis: an rDNA-based molecular phylogentic study. J Fish Dis. 2003;26:667–76.
Freeman M, Sommerville C. Desmozoon lepeophtherii n. gen., n. sp., (Microsporidia: Enterocytozoonidae) infecting the salmon louse Lepeophtheirus salmonis (Copepoda: Caligidae). Parasit Vectors. 2009;2:58.
Nylund S, Nylund A, Watanabe K, Arnesen CE, Karlsbakk E. Paranucleospora theridion n. gen., n. sp. (Microsporidia, Enterocytozoonidae) with a life cycle in the salmon louse (Lepeophtheirus salmonis, Copepoda) and Atlantic Salmon (Salmo salar). J Eukaryot Microbiol. 2010;57:95–114.
Freeman M, Sommerville C. Original observations of Desmozoon lepeophtherii, a microsporidian hyperparasite infecting the salmon louse Lepeophtheirus salmonis, and its subsequent detection by other researchers. Parasit Vectors. 2011;4:231.
Jones SRM, Prosperi-Porta G, Kim E. The diversity of microsporidia in parasitic copepods (Caligidae: Siphonostomatoida) in the northeast pacific ocean with description of Facilispora margolisi n. g., n. sp. and a new family Facilisporidae n. fam. J Eukaryot Microbiol. 2012;59:206–17.
Sutherland BJG, Jantzen SG, Yasuike M, Sanderson DS, Koop BF, Jones SRM. Transcriptomics of coping strategies in free-swimming Lepeophtheirus salmonis (Copepoda) larvae responding to abiotic stress. Mol Ecol. 2012;21:6000–14.
Poley JD, Braden LM, Messmer AM, Whyte SK, Koop BF, Fast MD. Cypermethrin exposure induces metabolic and stress-related gene expression in copepodid salmon lice (Lepeophtheirus salmonis). Comp Biochem Physiol - Part D Genomics Proteomics. 2016;20:74–84.
Braden LM, Sutherland BJG, Koop BF, Jones SRM. Enhanced transcriptomic responses in the Pacific salmon louse Lepeophtheirus salmonis oncorhynchi to the non-native Atlantic Salmon Salmo salar suggests increased parasite fitness. BMC Genomics. 2017;18:110.
Sutherland BJG, Poley JD, Igboeli OO, Jantzen JR, Fast MD, Koop BF, et al. Transcriptomic responses to emamectin benzoate in Pacific and Atlantic Canada salmon lice Lepeophtheirus salmonis with differing levels of drug resistance. Evol Appl. 2015;8:133–48.
Poley JD, Sutherland BJG, Jones SRM, Koop BF, Fast MD. Sex-biased gene expression and sequence conservation in Atlantic and Pacific salmon lice (Lepeophtheirus salmonis). BMC Genomics. 2016;17:483.
Kelly A, Hatcher MJ, Dunn AM. The impact of a vertically transmitted microsporidian, Nosema granulosis on the fitness of its Gammarus duebeni host under stressful environmental conditions. Parasitology. 2003;126:119–24.
Andreadis TG. Evolutionary strategies and adaptations for survival between mosquito-parasitic microsporidia and their intermediate copepod hosts: A comparative examination of Amblyospora connecticut and Hyalinocysta chapmani (Microsporidia: Amblyosporidae). Folia Parasitol. 2005;52:23–35.
Dunn AM, Smith JE. Microsporidian life cycles and diversity: The relationship between virulence and transmission. Microbes Infect. 2001;3:381–8.
Terry RS, Smith JE, Bouchon D, Rigaud T, Duncanson P, Sharpe RG, et al. Ultrastructural characterisation and molecular taxonomic identification of Nosema granulosis n. sp., a transovarially transmitted feminising (TTF) microsporidium. J Eukaryot Microbiol. 1999;46:492–9.
Terry RS, Smith JE, Sharpe RG, Rigaud T, Littlewood DTJ, Ironside JE, et al. Widespread vertical transmission and associated host sex-ratio distortion within the eukaryotic phylum Microspora. Proc Biol Sci. 2004;271:1783–9.
Hacker C, Howell M, Bhella D, Lucocq J. Strategies for maximizing ATP supply in the microsporidian Encephalitozoon cuniculi: direct binding of mitochondria to the parasitophorous vacuole and clustering of the mitochondrial porin VDAC. Cell Microbiol. 2014;16:565–79.
Van Der Giezen M, Tovar J, Clark CG. Mitochondrion-derived organelles in protists and fungi. Int Rev Cytol. 2005;244:175–225.
Dussaubat C, Brunet J-L, Higes M, Colbourne JK, Lopez J, Choi J-H, et al. Gut pathology and responses to the microsporidium Nosema ceranae in the honey bee Apis mellifera. PLoS One. 2012;7:e37017.
Holt HL, Aronstein KA, Grozinger CM. Chronic parasitization by Nosema microsporidia causes global expression changes in core nutritional, metabolic and behavioral pathways in honey bee workers (Apis mellifera). BMC Genomics. 2013;14:799.
Mayack C, Naug D. Energetic stress in the honeybee Apis mellifera from Nosema ceranae infection. J Invertebr Pathol. 2009;100:185–8.
Clark KF, Greenwood SJ. Next-generation sequencing and the crustacean immune system: The need for alternatives in immune gene annotation. Integr Comp Biol. 2016;56:1113–30.
Aufauvre J, Misme-Aucouturier B, Viguès B, Texier C, Delbac F, Blot N. Transcriptome analyses of the honeybee response to Nosema ceranae and insecticides. PLoS One. 2014;9:e91686.
Biron DG, Agnew P, Marché L, Renault L, Sidobre C, Michalakis Y. Proteome of Aedes aegypti larvae in response to infection by the intracellular parasite Vavraia culicis. Int J Parasitol. 2005;35:1385–97.
Bakowski MA, Desjardins CA, Smelkinson MG, Dunbar TA, Lopez-Moyado IF, Rifkin SA, et al. Ubiquitin-mediated response to microsporidia and virus Infection in C. elegans. PLoS Pathog. 2014;10:e1004200.
Agnew P, Berticat C, Bedhomme S, Sidobre C, Michalakis Y. Parasitism increases and decreases the costs of insecticide resistance in mosquitoes. Evolution. 2004;58:579–86.
Coors A, Decaestecker E, Jansen M, De Meester L. Pesticide exposure strongly enhances parasite virulence in an invertebrate host model. Oikos. 2008;117:1840–6.
Wu JY, Smart MD, Anelli CM, Sheppard WS. Honey bees (Apis mellifera) reared in brood combs containing high levels of pesticide residues exhibit increased susceptibility to Nosema (Microsporidia) infection. J Invertebr Pathol. 2012;109:326–9.
Firth KJ, Johnson SC, Ross NW. Characterization of proteases in the skin mucus of Atlantic salmon (Salmo salar) infected with the salmon louse (Lepeophtheirus salmonis) and in whole-body louse homogenate. J Parasitol. 2000;86:1199–205.
Fast AMD, Burka JF, Johnson SC, Ross NW. Enzymes released from Lepeophtheirus salmonis in response to mucus from different salmonids. Am Soc Parasitol. 2003;89:7–13.
Fast MD, Johnson SC, Eddy TD, Pinto D, Ross NW. Lepeophtheirus salmonis secretory/excretory products and their effects on Atlantic salmon immune gene regulation. Parasite Immunol. 2007;29:179–89.
Edvardsen RB, Dalvin S, Furmanek T, Malde K, Mæhle S, Kvamme BO, et al. Gene expression in five salmon louse (Lepeophtheirus salmonis, Krøyer 1837) tissues. Mar Genomics. 2014;18:39–44.
Hu A-K, Chen R-J, Zhu X-R, Gu J-H, Liu Z-P. Immunosuppressive effect of hypodermin C on complement component 3 in vitro. Cell Biochem Biophys. 2015;72:93–8.
Saksida SM, Morrison D, McKenzie P, Milligan B, Downey E, Boyce B, et al. Use of Atlantic salmon, Salmo salar L., farm treatment data and bioassays to assess for resistance of sea lice, Lepeophtheirus salmonis, to emamectin benzoate (SLICE®) in British Columbia, Canada. J Fish Dis. 2013;36:515–20.
Kumar S, Jones M, Koutsovoulos G, Clarke M, Blaxter M. Blobology: Exploring raw genome data for contaminants, symbionts and parasites using taxon-annotated GC-coverage plots. Front Genet. 2013;4:237.
Økland AL, Nylund A, Øvergård AC, Blindheim S, Watanabe K, Grotmol S, et al. Genomic characterization and phylogenetic position of two new species in Rhabdoviridae infecting the parasitic copepod, salmon louse (Lepeophtheirus salmonis). PLoS One. 2014;9:e112517.
Barker DE, Braden LM, Coombs MP, Boyce B. Preliminary studies on the isolation of bacteria from sea lice, Lepeophtheirus salmonis, infecting farmed salmon in British Columbia, Canada. Parasitol Res. 2009;105:1173–7.
Jones SRM, Kim E, Dawe S. Experimental infections with Lepeophtheirus salmonis (Kroyer) on threespine sticklebacks, Gasterosteus aculeatus L., and juvenile Pacific salmon, Oncorhynchus spp. J Fish Dis. 2006;29:489–95.
Jones SRM, Goh B, Prosperi-Porta G. Duration and method of fixation affects the sensitivity of a digoxygenin-labelled DNA probe in detecting Kudoa thyrsites in Atlantic salmon skeletal muscle. Aquaculture. 2003;220:157–64.
Yasuike M, Leong J, Jantzen SG, von Schalburg KR, Nilsen F, Jones SRM, et al. Genomic resources for sea lice: analysis of ESTs and mitochondrial genomes. Mar Biotechnol. 2012;14:155–66.
Huang D, Sherman B, Lempick R. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008;4:44–57.
Huang DW, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, et al. DAVID Bioinformatics Resources: Expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007;35:169–75.
Huang D, Sherman B, Lempick R. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:RESEARCH0034.
Team RC. A language and environment for statistical computing. R Found. Stat. Comput. Vienna, Austria. ISBN 3–900051–07-0, URL http://www.R-project.org/. 2016.
UniProt: A hub for protein information. Nucleic Acids Res. 2015;43:204–12.
We thank Sara Purcell for help in designing and executing PCR assays to detect F. margolisi from total RNA, Eliah Kim for collecting ovigerous lice and for culturing and filtering salmon lice larvae, and Lisa Zhang for assisting in the design and implementation of the in situ hybridization assay. Thanks to Stuart Jantzen and Niko Rebenich for contributing to the study design figure (Fig. 2) and for discussions on the analysis.
Gene expression data has been uploaded to GEO (NCBI) under the accession GSE94692.
This work was supported by Genome BC, Fisheries and Oceans Canada (DFO), Novartis/Elanco Animal Health (MDF), ACOA (AIF-199308; MDF) and NSERC discovery (RGPIN/402288–12; MDF). During this work, JDP was supported by an NSERC PGS D3 and BJGS was supported by an NSERC CGS graduate fellowship. In no case was a funding body involved in the design of the study, in the collection, analysis, and interpretation of data, or in the writing of the manuscript.
Availability of data and materials
The datasets supporting the results of this article are available in the NCBI Gene Expression Omnibus (GEO) repository (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi) GSE94692 for the copepodid data and GSE73734  for the pre-adult data.
This study was carried out in accordance with the recommendations of the Canadian Council for Animal Care. All procedures were approved by the Pacific Regional Animal Care Committee at the Pacific Biological Station under the Animal Use Protocol 11–004. Chum salmon used in experiments were obtained from the Nanaimo River Hatchery.
Consent for publication
The authors declare that they have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Differentially expressed transcripts by microsporidia and EMB. (XLSX 244 kb)
Gene Ontology. (XLSX 30 kb)
Interaction of microsporidia and EMB within k-means clustering. (XLSX 95 kb)
Consensus transcripts within Cluster B across L. salmonis studies. (XLSX 107 kb)
Pre-adult L. salmonis responses to F. margolisi infection and microsporidian genes on the microarray. (XLSX 85 kb)
Primers for RT-qPCR validation and in situ hybridization. (XLSX 10 kb)