Skip to main content

Advertisement

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 volume 18, Article number: 110 (2017) Cite this article

Abstract

Background

Outcomes of infections with the salmon louse Lepeophtheirus salmonis vary considerably among its natural hosts (Salmo, Oncorhynchus spp.). Host-parasite interactions range from weak to strong host responses accompanied by high to low parasite abundances, respectively. Parasite behavioral studies indicate that the louse prefers the host Atlantic Salmon (Salmo salar), which is characterized by a weak immune response, and that this results in enhanced parasite reproduction and growth rates. Furthermore, parasite-derived immunosuppressive molecules (e.g., proteases) have been detected at higher amounts in response to the mucus of Atlantic Salmon relative to Coho Salmon (Oncorhynchus kisutch). However, the host-specific responses of the salmon louse have not been well characterized in either of the genetically distinct sub-species that occur in the Atlantic and Pacific Oceans.

Results

We assessed and compared the transcriptomic feeding response of the Pacific salmon louse (L. salmonis oncorhynchi,) while parasitizing the highly susceptible Atlantic Salmon and Sockeye Salmon (Oncorhynchus nerka) or the more resistant Coho Salmon (Oncorhynchus kisutch) using a 38 K oligonucleotide microarray. The response of the louse was enhanced both in the number of overexpressed genes and in the magnitude of expression while feeding on the non-native Atlantic Salmon, compared to either Coho or Sockeye Salmon. For example, putative virulence factors (e.g., cathepsin L, trypsin, carboxypeptidase B), metabolic enzymes (e.g., cytochrome B, cytochrome C), protein synthesis enzymes (e.g., ribosomal protein P2, 60S ribosomal protein L7), and reproduction-related genes (e.g., estrogen sulfotransferase) were overexpressed in Atlantic-fed lice, indicating heightened parasite fitness with this host species. In contrast, responses in Coho- or Sockeye-fed lice were more similar to those of parasites deprived of a host. To test for host acclimation by the parasite, we performed a reciprocal host transfer experiment and determined that the exaggerated response to Atlantic Salmon was independent of the initial host species, confirming our conclusion that the Pacific salmon louse exhibits an enhanced response to Atlantic Salmon.

Conclusions

This study characterized global transcriptomic responses of Pacific salmon lice during infection of susceptible and resistant hosts. Similar parasite responses during infection of Coho or Sockeye Salmon, despite differences in natural immunity to infection between these host species, indicate that host susceptibility status alone does not drive the parasite response. We identified an enhanced louse response after feeding on Atlantic Salmon, characterized by up-regulation of virulence factors, energy metabolism and reproductive-associated transcripts. In contrast, the responses of lice infecting Coho or Sockeye Salmon were weaker, with reduced expression of virulence factors. These observations indicate that the response of the louse is independent of host susceptibility and suggest that co-evolutionary host-parasite relationships may influence contemporary host-parasite interactions. This research improves our understanding of the susceptibility of Atlantic Salmon and may assist in the development of novel control measures against the salmon louse.

Background

In a host–parasite relationship at equilibrium, parasite fitness is optimized and costs to the host are minimized. Disequilibrium of this relationship can cause harm to the host in the case of heightened virulence [1], or conversely, rejection of the parasite caused by heightened host immunity [2]. Aggressive host environments (e.g. a robust immune response) have been shown to decrease parasite fitness. For example, reduced reproductive output in ticks (Rhipicephalus microplus) is associated with feeding on more resistant hosts [3]. Increased parasite fitness while parasitizing more susceptible hosts is also observed in the differential developmental rates of the parasitic copepods Caligus rogercresseyi and Lepeophtheirus salmonis while infecting susceptible species (Oncorhynchus mykiss and Salmo salar, respectively), compared with that of resistant species [4, 5]. In addition, fewer low molecular weight proteases (virulence factors) are secreted by L. salmonis in response to mucus from resistant salmon (O. kisutch) relative to that from more susceptible hosts such as S. salar or O. mykiss [6] indicating that in addition to reducing reproductive output and development rate, an aggressive host response may also interfere with physiological responses associated with parasite feeding.

The salmon louse, L. salmonis, is a naturally occurring parasitic copepod that parasitizes anadromous salmonids belonging to the genera Oncorhynchus and Salmo spp, and is consequently an important pest of salmonid mariculture [7–9] throughout the Northern Hemisphere with sub-species in the Pacific (L. salmonis oncorhynchi) and Atlantic (L. salmonis salmonis) Oceans [10]. During heavy infestations, degradation of the epidermis and mucosal layer leads to osmoregulatory distress, anaemia, lethargy, secondary infections and a general stress response [11–13]. Resistance to the parasite varies among juvenile salmon such that Coho Salmon (Oncorhynchus kisutch) [14, 15] and Pink Salmon (Oncorhynchus gorbuscha) [16–18] display a resistant phenotype characterized by well-developed inflammation at the attachment site, rapid parasite rejection and limited additional pathology. In contrast, Atlantic Salmon (S. salar) [15, 17–21], Chum Salmon (Oncorhynchus keta) [16–18], and Sockeye Salmon (Oncorhynchus nerka) [15, 22] display a susceptible phenotype characterized by weak or absent local inflammation, higher parasite burden and associated pathology.

The differing host responses to the salmon louse suggest that there are also host-specific parasite responses. Behavioural studies indicate a preference of L. salmonis for salmonid over non-salmonid hosts [23, 24], and in particular a preference for Atlantic Salmon [25, 26]. Secretions from L. salmonis elicited by Atlantic Salmon mucus or dopamine extraction contain prostaglandin E2 (PGE2) and trypsin proteases, which are potent immune-modulators [27]. However, there are many questions that remain including the mechanisms involved in the differential host responses, and whether the secretion of virulence factors is influenced by the host species. Furthermore, most of the foundational studies on the differences in susceptibility in hosts (described above) have not considered the co-evolutionary history of the host-parasite relationship.

There is a significant knowledge gap regarding the response of the salmon louse during feeding. To this end, we applied a transcriptomic approach to assess the effect of host resistance on louse feeding responses. We hypothesized that the transcriptomic feeding response on more susceptible species (Atlantic, Sockeye Salmon) would be enhanced in virulence factors and fitness-related pathways (e.g. feeding, reproduction, energy metabolism) relative to that elicited by a resistant species (Coho Salmon). However, our results indicate a specific enhanced response to Atlantic Salmon that was not explained by host acclimation but that may be due to differences in co-evolutionary history of the parasite and the hosts.

Results

Using a L. salmonis 38 K oligonucleotide microarray (eArray, Agilent) designed with expressed sequence tags (ESTs) from L. salmonis salmonis and L. salmonis oncorhynchi [28], we detected 15,718 probes that passed quality control filters. Excluding duplicate probes, this included 8,776 unique transcripts. Differences in the expression of these genes were investigated L. salmonis infecting Atlantic, Coho, or Sockeye Salmon and in lice withheld from a host (i.e., starved; Fig. 1a).

Fig. 1
figure1

Experimental design. In the host-effect hypothesis experiment (a) Atlantic (AT), Coho (CO), and Sockeye (SK) Salmon were infected with L. salmonis sourced during commercial Atlantic Salmon harvest. A sub-set of unattached lice served as the starved (STV) control. At 24 and 48 h lice (n = 10) were removed from every species and processed for down-stream microarray and RT-qPCR analysis. In the acclimation hypothesis experiment (b) L. salmonis were sourced from Atlantic Salmon (AT) during harvest as well as from Sockeye Salmon (SK) from a test fishery. Lice from each species (AT, SK) were used to infect both Atlantic and Sockeye Salmon. At 24 and 48 h, lice (n = 8) were removed and processed for downstream RT-qPCR analysis

Full size image

Profiling the feeding response of the salmon louse

The total number of differentially expressed genes (DEGs) at 24 and 48 hpi in lice parasitizing each host (Atlantic, Coho, and Sockeye Salmon) was compared to a group of lice withheld from hosts (i.e., ‘starved’; Additional file 1: Table S1). The feeding response of L. salmonis was determined by identifying DEGs common to lice parasitizing salmon regardless of species (n = 10 individual lice per host species per time point; total = 60 feeding lice and 20 starved). Gene Ontology analysis of overexpressed feeding DEGs revealed enrichment in biological processes such as proteolysis (GO:0006508), hatching (GO:0035188), blood coagulation (GO:0007596), oxidation reduction (GO:0055114), and collagen metabolic process (GO:0032963) (Fig. 2).

Fig. 2
figure2

Profiling the feeding response of L. salmonis. Overexpressed transcripts in lice feeding on Atlantic, Coho or Sockeye Salmon relative to starved lice were compared to produce a list of unique genes involved in the feeding response of L. salmonis. These genes were analyzed using DAVID to produce GO enrichment results after a 24 and b 48 hpi. aFold Enrichment

Full size image

Using RPS-BLAST and the Conserved Domain Database [29] as previously described [30], conserved domains were identified within the feeding response genes. This analysis revealed 67 distinct domains overexpressed at 24 hpi with frequently identified domains including “trypsin-like” (29 unique contigs, smart00020), “saposin-like” (7 unique contigs, smart00741), and “cysteine-like” (4 unique contigs, cd03860). At 48 hpi there were 71 distinct domains, with the most frequently identified domains including “trypsin-like” (19 unique contigs, smart00020) and “glutathione peroxidase-like” (7 unique contigs, cd00340). As expected, these domains are associated with feeding and digestion. The “trypsin-like serine protease” domain (smart00020, cd00190) was present in 14 genes including trypsin-1, collagenase, chymotrypsin B1, coagulation factor IX and hypodermin B. Also well represented were “peptidase” domains, including “peptidase C1” (e.g., cathepsin L; pfam00112), “peptidase C13” (e.g., legumain; pfam01650), “peptidase S28” (e.g., putative serine protease K12H4.7; pfam05577), “peptidase M13” (e.g., neprilysin-2; pfam01431), and “peptidase M14” (e.g., carboxypeptidase B; cd03860) (Additional file 2).

Although expression of these transcripts was evident during feeding on salmon regardless of species, Atlantic-fed lice exhibited the highest expression of proteolytic-, metabolic- (oxidative-reduction), and reproductive-associated (hatching) genes (Fig. 3).

Fig. 3
figure3

Enhanced transcriptomic response in Atlantic-fed lice. Hierarchical clustering of transcripts at 48 hpi enriched in GO biological process categories of a proteolysis, b oxidative reduction, and c hatching revealed an enhanced response specific to L. salmonis feeding on Atlantic Salmon (AT). Expression levels are mean log2 normalized intensities of each transcript, with red and green representing low and high expression, respectively. These genes were significantly differentially expressed in response to feeding on Atlantic Salmon compared to starved (p < 0.01; FC ≥ 1.5)

Full size image

Some specific transcripts were identified in Atlantic-fed lice that may be related to virulence in the salmon louse. For example, we detected a “phospholipase A2 bee venom-like” domain (cd04704) that was overexpressed at 48 hpi in Atlantic-fed lice compared to either Coho- or Sockeye-fed lice (FC = 2.03, 2.38, respectively). Additionally, “saposin B-like” domains (smart00741), known to be important in lipid-interacting proteins such as NK-lysins, were well represented on the array (e.g., 12 contigs passing fold-change filters at 48 hpi), and expression of this gene was significantly higher in Atlantic-fed lice (FC = 2.57–5.31, p < 0.001) compared to Coho-fed lice (FC = 1.5–1.94, p < 0.001), while Sockeye-fed lice expression of these domains did not pass the fold-change filter.

Profiling the starvation response of the salmon louse

To determine the response to starvation, we evaluated overexpressed transcripts in starved lice compared to lice feeding on salmon, regardless of species. DEGs overexpressed in starved relative to salmon-fed lice represented the starvation response (Additional file 3: Figure S1). There were 47 and 143 transcripts overexpressed in starved L. salmonis at 24 and 48 h, respectively. Gene Ontology analysis of these unique “starved” DEGs revealed one significantly enriched category at 24 h (GO:0007165, signal transduction), and nine categories at 48 h including negative regulation of transcription (GO:0045892), cell differentiation (GO:0030154), and sarcomere organization (GO:0045214) (Additional file 3: Figure S1).

Host-specific feeding responses of the salmon louse

We then assessed the host-specific responses of L. salmonis by focusing on the DEGs with FC ≥ 1.5 in Atlantic- compared to either Coho- or Sockeye-fed lice. In Atlantic-fed lice, an enhanced response was identified with an enrichment for digestion, reproduction and energy metabolism. Conversely in either Coho- or Sockeye-fed lice, the response was enriched for skeletal and sensory system development (Sockeye Salmon), or cytoskeletal organization, regulation of growth and nucleosome organization (Coho Salmon) (Additional file 4: Table S2).

At 24 hpi the most significantly enriched GO category in Atlantic-fed lice (p = 6.5 × 10−4 compared to Coho, p = 3.6 × 10−9 compared to Sockeye) was proteolysis (GO:0006508), and included proteases (cathepsin L, trypsin-1, chymotrypsin A chain C, aspartic proteinase oryzasin-1, anionic trypsin-1), carboxypeptidases (carboxypeptidase B, zinc carboxypeptidase A1), and matrix metalloproteases (matrix metalloprotease 2, zinc metalloprotease nas-4) (Table 1). In Coho-fed lice, the most significantly enriched GO category (p = 0.013 compared to Atlantic) was cytoskeletal organization (GO:0007010) and included genes like troponin C isoform 1, and myosin heavy chain. In Sockeye-fed lice, the most significantly enriched GO category (p = 0.006 compared to Atlantic) was skeletal system development (GO:0001501) and included genes like zinc finger protein 16, and bone morphogenetic protein 2-B (Additional file 4: Table S2).

Table 1 Proteolytic genes specific to Atlantic-fed lice
Full size table

At 48 hpi, there was no significant enrichment of biological processes in Coho- or Sockeye-fed lice, whereas enriched pathways remained abundant and highly populated in Atlantic-fed lice. Some of the most pronounced overexpression specific to Atlantic-fed lice was observed in the genes involved in energy metabolism and protein synthesis, including genes enriched in the GO category oxidative reduction (GO:0055114), including cytochrome C oxidase subunit 2 (FC = 4.1–8.4 compared to Coho-fed) and 3 (FC = 6.7–7.3 and 24.1–33.2, compared to Coho- and Sockeye-fed, respectively), and cytochrome B (FC = 19.8–32.3, 37.7–41.8, compared to Coho- and Sockeye-fed, respectively). Additionally, genes associated with protein synthesis such as ribosomal protein P2 and 60S ribosomal protein L7, were significantly overexpressed in Atlantic-fed lice (FC = 7.0–8.1 and 5.7, compared to Coho-fed lice, respectively; FC = 6.7–8.6 and 3.6, compared to Sockeye-fed lice, respectively) (Additional file 2).

Enrichment of genes associated with reproductive-type processes was also specific to Atlantic-fed lice and included the biological process GO categories of hatching (GO:0035188), female pregnancy (GO:0007565), blastocyst development (GO:0001824) and blastocyst hatching (GO:0001835). These categories were enriched in the response to Atlantic Salmon but not to either Coho or Sockeye Salmon, and included genes such as placental protein 11, neutral ceramidase, granulin-7, and estrogen sulfotransferase (Additional file 4: Table S2).

To characterize genes responding differently over time during infection with each species, k-means clustering of DEGs was performed and this further indicated the responses of Sockeye- and Coho-fed lice were similar to the starved lice and distinct from the Atlantic-fed lice. In one of the five clusters, 45 transcripts associated with stress, including heat shock protein 90, T-complex protein 1 subunit alpha, heat shock protein SSA1, and heat shock protein homolog ECU03_0520 were up-regulated over time in Coho- and Sockeye-fed lice while down-regulated in Atlantic-fed lice (Fig. 4). Transcripts involved in feeding, energy metabolism and reproduction were most highly expressed while feeding on Atlantic Salmon (Figs. 3 and 5). Thus the feeding response of L. salmonis indicated that Atlantic Salmon is the most desirable host compared Coho or Sockeye Salmon.

Fig. 4
figure4

Similar expression profiles over time of stress-related transcripts in Coho-fed, Sockeye-fed and starved L. salmonis. Using k-means clustering analysis, genes associated with stress were shown to be up-regulated over time in lice feeding on Coho and Sockeye Salmon, and by lice withheld from hosts (Starved). In contrast, these genes were down-regulated over time in lice feeding on Atlantic Salmon

Full size image
Fig. 5
figure5

Enhanced expression of feeding and energy in Atlantic-fed lice. Differentially expressed transcripts identified by the microarray were profiled using RT-qPCR. Expression of genes involved in proteolysis/digestion (a) and energy metabolism (b) are shown as log2 calibrated normalized relative quantities (CNRQ). Expression in Atlantic-fed lice (AT) increased over time (24 → 48 hpi) and was significantly higher than in Coho-fed (CO) or Sockeye-fed (SK) lice. Expression of SK-fed and CO-fed lice was not significantly different from starved lice (ST). Significance was identified by two-way ANOVA (p < 0.05) followed by post-hoc Tukey HSD test to determine pairwise significance. Within a time point, lower-case letters denote differences between groups where groups that do not share a letter are significantly different. Asterisks denote differences within a group between time points

Full size image

Temporal activation of the feeding response

We hypothesized that while parasitizing a more optimal host, responses associated with increased fitness (i.e., digestion, energy, reproduction) would increase over time (i.e. 24 → 48 hpi). The responses of Atlantic-fed lice increased over time, measured by the number of DEGs and the magnitude of expression of the DEGs. In contrast, the response in Coho-fed or Sockeye-fed lice either did not change, or was reduced. There was significant enrichment of up-regulated genes in the biological process categories such as proteolysis, digestion, and oxidative reduction in Atlantic-fed lice while enrichment of down-regulated transcripts was only observed in Coho-fed lice. There was no enrichment for up- or down-regulation over time in Sockeye-fed lice (Table 2).

Table 2 Enrichment in genes over time (24 → 48 hpi) in feeding L. salmonis
Full size table

Conserved domain analysis was used to characterize the types of genes that were significantly different over time. Of the 435 genes affected by time (main effect time, no interaction effect), 211 were unannotated, 198 were annotated with domains, while 26 were annotated but did not contain a conserved domain. In total, 99 domains were present in the genes significantly overexpressed by L. salmonis over time (Additional file 2). Most domains were only represented by one gene, although 7 genes contained “trypsin-like serine protease” (smart00020, cd00190) domains, and 2 contained “peptidase M14 carboxypeptidase subfamily A/B-like” (cd03860), “FKBP-type peptidyl-prolyl cis-trans isomerase” (pfam00254), and “fasciclin” (pfam02469) domains. Genes with “trypsin-like serine protease” domains included trypsin-1, ovochymase-2, neurotrypsin, chymotrypsin A, and trypsin-like serine protease. Importantly, although these genes changed over time and contained similar protein domains, they did not change in the same way. For example, expression of trypsin-1 increased over time in Atlantic-fed lice, and either declined or remained constant in Coho- or Sockeye-fed lice. In contrast, expression of ovochymase-2 increased in Sockeye-fed lice, but remained constant in Atlantic- or Coho-fed lice.

Positively correlated expression profiles were observed for genes containing other domains. For example, genes with the “metallocarboxypeptidase” domain (cd03860; carboxypeptidase B and zinc carboxypeptidase A1) were both down-regulated over time in Coho-fed lice, and although expression in Sockeye-fed lice did not pass fold-change filters, there was a trend towards negative regulation for both these genes (FC = −1.4, −1.3, respectively).

Several other domains associated with proteolytic-enzymes were up-regulated over time only in Atlantic-fed lice including “peptidase M14 carboxypeptidase subfamily N/E-like” (cd03868), “zinc-dependent metalloprotease astacin-like subfamily” (cd04280), “papain family cysteine protease” (pfam00112), and “serine carboxypeptidase S28” (pfam05577).

One gene, tyrosine aminotransferase (TIGR01265), was concordantly down-regulated over time in Atlantic- (FC = −4.0), Coho- (FC = −1.6) or Sockeye-fed lice (FC = −2.0), indicating this gene is involved in the starvation stress response. Similarly, the “alpha crystallin” domain (cd06464; programmed cell death 4) was specifically associated with starvation.

Enhanced transcriptomic response to Atlantic Salmon is not explained by acclimation

We considered the possibility that enhanced transcriptomic responses in Atlantic-fed lice were explained by acclimation to the host species, as the parasites had been collected from Atlantic Salmon. In a follow-up experiment, we assessed the expression of a subset of genes identified from the microarray experiment described above by exposing Atlantic and Sockeye Salmon to lice collected from either Atlantic or Sockeye Salmon (Fig. 1b). In the event of acclimation, we predicted enhanced responses in lice attached to salmon belonging to the same species as the source host. Contrary to this, our results showed enhanced responses of Atlantic-fed lice regardless of the source host species, for genes associated with protein synthesis (ribosomal protein L2, 60S ribosomal protein L7; Additional file 5: Figure S2A), metabolism (cytochrome oxidase subunit 2, cytochrome B; Additional file 5 : Figure S2B), and proteolysis (cathepsin L; Additional file 5: Figure S2C). This indicates that the observed trends in expression were not due to acclimation, but rather to host-specific factors.

Validation of the microarray

To confirm findings from microarray analyses, transcript abundance was analyzed for a sub-set of DEGs using RT-qPCR. Genes with potential relation to feeding and energy production in L. salmonis were chosen for validation. There was high correlation between the RT-qPCR and microarray data (p < 0.05, n = 14 gene comparisons, Additional file 6: Table S3). Additionally, temporal trends observed in the genes from Atlantic-fed lice by microarray analysis were also detected by RT-qPCR analysis which showed overexpression of proteases (cathepsin L, trypsin-1, neprilysin-1; Fig. 5a), mitochondrial enzymes (cytochrome B, cytochrome C oxidase subunit 2; Fig. 5b), ribosomal proteins (ribosomal protein P2, ribosomal protein 60S; Fig. 6a), and oxidative stress-associated genes (ferritin, high affinity copper uptake protein 11; Fig. 6b). In contrast, in the RT-qPCR results, the overexpression of stress-associated genes (programmed cell death 4, T-complex protein 1; Fig. 6c) was specific to starved lice.

Fig. 6
figure6

Enhanced expression of protein synthesis and oxidative-stress genes in Atlantic-fed lice. Differentially expressed transcripts identified by the microarray were profiled using RT-qPCR, and are shown as log2 calibrated normalized relative quantities (CNRQ). Expression of genes involved in protein synthesis (a) and oxidative stress (b) were the highest in AT-fed lice and increased over time (24 → 48 hpi). (c) Expression of two genes associated with stress were highest in starved lice compared to either AT-, CO- or SK-fed lice. Significance was identified by two-way ANOVA (p < 0.05) followed by post-hoc Tukey test to determine pairwise significance. Within a time point, lower-case letters denote differences between groups where groups that do not share a letter are significantly different. Asterisks denote differences within a group between time points

Full size image

Discussion

This study tested the hypothesis that transcriptomic responses of adult female L. salmonis on susceptible salmon are enhanced relative to the responses measured on resistant salmon. Furthermore, we predicted that characteristics of the enhanced response would be consistent with parasite fitness. Compared to Atlantic Salmon-fed lice, differential expression of genes in lice feeding on either the susceptible or resistant Pacific salmon species was relatively weak with low magnitude fold-changes and sparsely populated Gene Ontology categories. Moreover, in lice feeding on the Pacific salmonids, there was no increase in the response over time and often gene expression profiles were more similar to those observed in starved L. salmonis (e.g., heat shock protein 90, tubulin alpha chain, and T-complex protein 1). The similarity of transcriptional responses in the Pacific salmon-fed lice, despite differences in natural resistance of Coho and Sockeye Salmon to L. salmonis, indicates that host resistance status does not explain the differential parasite response. We considered the possibility that the enhanced transcriptomic response in lice feeding on Atlantic Salmon was because of host acclimation as the lice were originally collected from Atlantic Salmon [31]. This possibility was addressed by conducting a reciprocal host exposure study using lice collected from Sockeye Salmon or Atlantic Salmon. In the event of acclimation, we expected an enhanced parasite response when lice were allowed to feed on the species from which they had been collected. However, we consistently observed that lice responded strongly to Atlantic Salmon, irrespective of the source host species, thus corroborating the host-effect hypothesis and further illustrating the desirability of Atlantic Salmon to the Pacific salmon louse.

The characteristics of the transcriptomic response to Atlantic Salmon suggest increased parasite fitness. Host blood is a main dietary component of the adult female salmon louse [32] suggesting a need for haemolytic enzymes and anti-coagulants to maintain a free flow of blood to the site of feeding as observed for other hematophagous parasites [33–37]. Our data show that proteolytic and other digestive-associated enzymes were a major component of the salmon louse response to Atlantic Salmon (Table 1). The overexpression of cathepsin L, trypsin-1, neprilysin-1, carboxypeptidase B, zinc carboxypeptidase A1, and legumain was more pronounced on Atlantic-fed L. salmonis, consistent with earlier reports of the secretion of proteases by this parasite [38–40]. Many of these genes are known virulence factors in other ectoparasites, and modulate the host immune response during feeding: cathepsin L is a virulence factor found in numerous parasites [41–43] that suppresses the host immune response while aiding in tissue digestion [42, 44]; legumain is important in blood digestion [41]; neprilysin is associated with regulating inflammation [45]; and carboxypeptidase-B prevents clotting and inhibits inflammation [46–48]. These virulence factors likely provide similar functions for L. salmonis.

We detected the expression of several other feeding-associated genes in L. salmonis that may represent virulence factors based on functions observed in other organisms. For example phospholipase A2 (PLA2) associated with feeding on Atlantic Salmon, is an important constituent of bee and snake venom [49, 50] and is also found in the secretions of hematophagous ectoparasites [49]. Interestingly, PLA2 induces a type-2 immune response in mice [50] and possesses potent hemolytic activity [37]. Another potential virulence factor detected in feeding lice was L-amino acid oxidase, a toxin found in snake venoms [51]. At 48 hpi, this gene was overexpressed in feeding L. salmonis irrespective of the host species. Finally, a saposin-B like protein was significantly overexpressed in the feeding salmon louse transcriptome, and most highly by parasites feeding on Atlantic Salmon. Saposin-like proteins (SAPLIPs) have been described from Fasciola spp. [52], Schistosoma mansoni [53], Entamoeba histolytica [54], and Amblyomma americanum [55] and are involved in cytolysis and lipid metabolism [56]. The high abundance of unique contigs containing the saposin-like protein domain indicates a need to further characterize the function of SAPLIPs in L. salmonis. Similarly, the potential for hemorrhagic, fibrinolytic, cytolytic, and apoptotic effects, among others, warrants further investigation of the role of L-amino acid oxidase and PLA2 in the salmon louse.

Genes related to oxidative stress and iron homeostasis (e.g., ferritin, high-affinity copper-uptake protein-1) were up-regulated in feeding L. salmonis, and expression was highest on Atlantic Salmon. In crustaceans, copper is critical as a cofactor for enzymes involved in many physiological processes including oxidative phosphorylation and mobilization of iron [57, 58]. The mitochondrial enzyme cytochrome c oxidase (cox) is a particularly abundant cuproprotein [59]. We observed prominent overexpression of metabolism-associated genes including cytochrome c oxidase (subunit 2 and 3) in lice feeding on Atlantic Salmon. Enhanced mitochondrial activity, combined with significant overexpression of protein synthesis-related genes (e.g., ribosomal protein P2, 60S ribosomal subunit), indicate that feeding on Atlantic Salmon was associated with pronounced metabolic activity.

Enhanced reproductive output is another proxy for higher parasite fitness, which influences host-choice and parasite virulence [2]. The salmon louse has been shown to prefer Atlantic Salmon, and parasite reproductive output and growth rates are increased while parasitizing Atlantic Salmon compared to Coho Salmon [6, 60]. We provide evidence for enrichment of reproductive-associated transcripts (e.g., placental protein-11, neutral ceramidase, granulin-7) in L. salmonis feeding on Atlantic Salmon. Further to this, enrichment for reproduction-like pathways was only present in Atlantic-fed lice compared to either Coho- or Sockeye-fed lice.

We show distinct responses to Atlantic and Pacific salmonids, possibly reflecting alternative mechanisms by which the parasite contributes to the outcome of infection. The enhanced exploitation of Atlantic Salmon may be due to differences in host-specific factors such as skin structure (e.g., low mucous cell density [61]) or physiology (e.g., low immune response [62]), or due to reduced genetic diversity in farmed populations [63, 64]. Furthermore, coevolution of the Pacific louse subspecies with Pacific Salmon for between 4.6 and 11 Ma [65] may have led to adaptations in the host-parasite relationships including variations in natural host resistance among salmon species, similar to what is observed in other host-parasite systems [66], and possibly related to host life history strategies [67–69] and life stages [70]. For example, juvenile Chum and Sockeye Salmon support high infections with L. salmonis and exhibit weakened cellular and humoral inflammatory responses at the louse attachment site, compared to those of juvenile Coho or Pink Salmon [15–18, 60]. Coho Salmon also exhibit heightened resistance towards other ectoparasitic copepods [4]. Similarly, juvenile Pink Salmon display a resistant phenotype that results in rapid rejection of the parasite [17, 18, 71, 72]. However mature Pink Salmon lose much of this natural resistance [70], and large abundances of the parasite are observed on mature Pink Salmon [73].

Our data support the hypothesis that Atlantic Salmon provides a host environment more permissive for L. salmonis fitness as shown by overexpression of transcripts related to virulence factors, energy metabolism, and reproduction. Functional enrichment for reproduction was observed during louse infection of Atlantic relative to either Coho or Sockeye. High energy metabolism is correlated with high reproductive output, and as such parasites have evolved to prefer hosts in which they are able to extract the highest available energy and maximize reproductive output [74]. Energy metabolism and reproductive output may be used as a proxy for the relative nutritional value of the host and our data suggest the nutritional gain from Atlantic Salmon exceeded that from Coho or Sockeye Salmon. Balancing host immunity with nutritive value is a driver of parasite host-specificity observed in many host-parasite relationships [75, 76], and may explain the observed host preference of L. salmonis for Atlantic Salmon.

Conclusions

Understanding host-specific feeding responses of L. salmonis may help explain the variable outcomes of infection among host species that have previously been associated with diverse host responses, such as the delayed or muted inflammation in Atlantic, Chum and Sockeye Salmon [13, 15–19, 77]. Contrary to our original expectations, the responses were similar for parasites feeding on Coho or Sockeye Salmon, despite differences in their susceptibility, but differed strongly with those of lice feeding on Atlantic salmon. We suggest that the prolonged evolutionary relationship between L. salmonis oncorhynchi and Oncorhynchus spp. explains the more limited capacity for parasite response. In contrast, the non-native Atlantic Salmon elicits an enhanced parasite feeding response, which may dampen local host response mechanisms and facilitate an environment more conducive to parasite fitness. The identification of host-specific factors involved in enhanced functional and transcriptomic lice responses (e.g. constituents of mucous) will further improve our knowledge on this system. Furthermore, a comparison of the responses of the Atlantic salmon louse (L. salmonis salmonis) on Pacific and Atlantic salmonids will be valuable to confirm the role of co-evolutionary host-parasite interactions in enhanced parasite fitness.

Methods

Experimental fish and infection procedures

Host-effect hypothesis

Atlantic Salmon parr were obtained from a commercial salmonid hatchery, Coho Salmon parr were obtained from the Chase River hatchery on Vancouver Island, British Columbia (B.C.), Canada and Sockeye Salmon parr were obtained from the Inch Creek hatchery, Chilliwack, B.C., as previously described [15]. All fish were reared in brackish water (~15 ppt) until smoltification, after which they were maintained on ultraviolet-treated salt water (~33 ppt) in single-pass flow-through tanks and on a 12:12 h light:dark cycle. Fish were fed 1% total biomass daily with a commercially available diet (EWOS). Fish of each species were randomly divided among 12 tanks (330 L), with four tanks used for each species (2X infection tanks, 2X control tanks). Fish were acclimated for approximately 7 days and starved at least 24 h prior to any manipulation as previously described [15].

All experiments utilized the Pacific salmon louse, Lepeophtheirus salmonis oncorhynchi, referred to here as L. salmonis. Adult female L. salmonis were collected during harvest of Atlantic Salmon at a commercial aquaculture site on Vancouver Island, B.C. After collection, the lice were rinsed in fresh seawater and transported to the Pacific Biological Station, Nanaimo, B.C. in 8 °C aerated previously-sterilised sea water. Only lice that were firmly attached to the collection container were included in the study. The time between collection of the salmon lice and initial infection was < 24 h [15].

For experimental infections, the water level of each tank was reduced by 50%, and fish (n = 25/species) were sedated in seawater containing 0.2 mg/L metomidate hydrochloride (M-HCl; Aquacalm, Syndel Laboratories). Sedated fish were transferred to a temporary tank containing M-HCl (0.2 mg/L) to which 5 lice/fish were added and allowed to settle and attach [15]. Once infected with 5 lice, fish were gently removed from the infection tank and returned to their original tank. Another group of lice were maintained at 8 °C in aerated seawater and thus served as non-attached controls.

At both 24 and 48 h post-infection (hpi), one louse was sampled from each of ten Atlantic (mean weight 218 ± 29 g), Coho (mean weight 192 ± 35 g) and Sockeye Salmon (mean weight 167 ± 17 g). These 60 lice comprised the feeding lice and a further 20 lice (10 from each sample time) comprised the non-feeding controls (i.e., starved lice, Fig. 1a). Upon sampling, each louse was individually snap-frozen in liquid nitrogen for gene expression profiling.

Acclimation hypothesis experiment

To test the importance of the source host species on subsequent parasite responses (i.e. parasite acclimation), a second experiment was conducted using lice collected from either Atlantic or Sockeye Salmon (Fig. 1b). Adult female L. salmonis were collected during harvest at a commercial aquaculture site (Atlantic-acclimated; L. salmonis-AT) or during a test fishery (Sockeye-acclimated; L. salmonis-SK) and transported to the Pacific Biological Station. Salmon of each species were divided into two tanks (4 tanks; n = 15 per tank) and allowed to acclimate for 7 days. Infections were performed as above except that for each species, fish in one tank were infected with L. salmonis-AT (n = 5 lice per fish) and those in the second tank infected with L. salmonis-SK (n = 5 lice per fish). At 24 and 48 hpi, 8 lice were collected from each tank (total number of lice = 32 per time) and snap-frozen in liquid nitrogen.

RNA extraction

Frozen lice were homogenized using 5 mm stainless-steel beads and a Tissue-lyser (Qiagen). RNA was extracted using TRIzol (Invitrogen) following manufacturer’s instructions with modifications. Specifically, following the organic phase extraction, the supernatant was removed and RNA was then purified using RNeasy spin columns (Qiagen) with an on-column DNase I digestion to remove genomic DNA as per manufacturers’ instruction. Total RNA was eluted in 30 μL ultra-pure water and quantified by spectrophotometry (Nanodrop-1000, Thermo Fisher). RNA quality was determined using Experion Automated Electrophoresis (Bio-Rad) with all samples having an RQI < 9.

cRNA synthesis and reference pool generation

Purified total RNA (200 ng) was reverse-transcribed to cDNA and then amplified to Cy5-labeled cRNA with Cy5-CTP-labeled nucleotides (Perkin Elmer) as previously described [78] using Low Input Quick Amp Labeling kits (Agilent), as per manufacturer’s instructions for hybridization to a 4-pack oligo gene expression microarray. Labelled cRNA was purified through RNeasy columns as per manufacturer’s instructions (QIAGEN) and quantified using spectrophotometry (NanoDrop-1000), ensuring specific activity of all samples > 6 pmol dye per microgram cRNA (Agilent). Samples were kept at −80 °C until hybridization. A reference pool of Cy3-cRNA was synthesized by amplifying experimental samples as described above, but with Cy3-CTP-labeled nucleotides (Perkin Elmer). In each experiment, a reference pool of equimolar cRNA was generated from each experimental condition (n = 10 individuals).

Microarray hybridization, quantification, normalization and filtering

A 38 K oligo microarray was designed using previously annotated ESTs from both Pacific and Atlantic L. salmonis [28] using eArray (Agilent). Probes were preferentially selected at 3′ untranslated regions and each EST was represented by duplicate probes (19 k ESTs represented on the array). Each individual louse was hybridized to a single array (i.e. total of 20, 4-pack microarrays hybridized in the host-effect experiment). Sample and reference combinations (825 ng cRNA each) were fragmented and then hybridized at 65 °C for 17 h at 10 rpm as per manufacturer’s instructions (Agilent) using SureHyb chambers (Agilent). Washing was performed as per manufacturers’ instructions, using the optional protocol to prevent ozone degradation. All slides were transferred to a dark box and kept at low ozone until scanned on a Perkin Elmer ScanArray® Express at 5 μm resolution using PMT settings optimized to have the median signal of ~1–2% of array spots saturated (Cy5: 65; Cy3: 68). Spot intensities were quantified in Imagene 8.1 (Biodiscovery) using an eArray GAL file (Design ID: 024389; Agilent). Poor spots and control spots were flagged by the software for downstream filtering. The background of each spot was subtracted from the foreground median, and samples were imported into GeneSpring 11.5.1 (Build 138755; Agilent). Each experiment was normalized and filtered separately 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 at least 65% of samples in any one condition and no flags in at least 65% of samples in any one condition as described elsewhere [78]. Raw quantified microarray files have been submitted to NCBI Gene Expression Omnibus (GEO) under the accession GSE80220.

Differential expression and functional analysis of microarray data

Probes were tested for differential expression using a two-way ANOVA without equal variance assumption using host-species and time as explanatory variables, followed by a post hoc Tukey’s HSD (p < 0.01) and using an FDR (Bonferroni’s test for multiple test correction) and fold change ≤ 1.5 from starved controls. All differentially expressed probes were used as an input for k-means clustering to identify co-expressed gene clusters (Euclidean distance metric; 5 clusters; 50 iterations; GeneSpring 11.5.1 Agilent). Gene Ontology (GO) and pathway enrichment were performed on annotated probes using the DAVID online bioinformatics tool (modified Fisher’s exact test) [79], with UniProt accessions of clustered probes compared to a background list of all probes passing quality control filters (n = 15,718 probes). Overlap between lists of differentially expressed genes was evaluated using VENNY [80]. After Gene Ontology enrichment analysis, GO Trimming was used to reduce redundancy of enriched Gene Ontology categories with a soft trim threshold of 40% [81].

Reverse-transcriptase quantitative polymerase chain reaction (RT-qPCR)

The same RNA samples analyzed with microarrays in the host effect experiment were used for RT–qPCR. This included a total of ten individual lice from each of four conditions (feeding on Atlantic Salmon, Coho Salmon, Sockeye Salmon, or not-feeding; see above for more details) at 24 and 48 h, to produce a total of 80 samples. Synthesis of cDNA was performed with 2 μg of total RNA in 20 μl reactions using oligo (dT) primers and AffinityScript cDNA Synthesis kits (Agilent), as per manufacturer’s instructions. Each cDNA sample was diluted 10-fold. To generate a standard curve, one sample from each of the four conditions (starved, Atlantic-fed, Coho-fed, Sockeye-fed) was randomly selected and synthesized as described previously [15, 78]. These samples were pooled and diluted 10-fold and this pool was then used for a serial dilution (5-point, 10-fold each point) for efficiency calculation. RT- qPCR amplification was performed using Brilliant UltraFast SYBR III® (Agilent) in 20 μl reactions with 0.1 μM of each primer using the following thermal regime: 95 °C for 3 min, followed by a combined annealing and extension step of 60 °C for 40 cycles. For qPCR technical replication, each sample/gene combination was run in triplicate. Genes of interest were selected from the microarray results based on biological relevance, high significance level or presence in enriched GO categories [78]. Reference gene candidates were selected from microarray results based on stable expression across conditions, consistency across replicate spots and moderate levels of expression, as well as from previous literature [78]. Primers were designed in Primer3 [82] selecting amplicon sizes of 80–150 base pairs (Additional file 7: Table S4). Amplicons were checked for single products by melt curve analysis and were sequenced on an ABI 3130 (Applied Biosystems) to confirm identity. RT–qPCR data analysis was performed using qBase-PLUS (Biogazelle). Stability of reference genes was tested using geNorm [83]. Selected reference genes included the previously identified structural ribosomal protein S20 and tubulin beta chain, with a collective M value of 0.382 and CV of 0.146, which is within the range typically observed for stably expressed reference genes in heterogeneous sample [83]. Other tested reference genes that were not used to normalize due to higher variability included elongation factor 1-α and HPGRT. NTC and RT controls showed no amplification. Statistical significance was identified by two-way ANOVA (p < 0.05) with pairwise significance determined by post-hoc Tukey test (SigmaPlot V11.1). Correlation between methods (RT–qPCR and array) was tested using by the correlation of log2 normalized expression values for RT–qPCR samples against microarray normalized log2 expression ratios (Cy5/Cy3) for the probe corresponding to the contig used for primer design, as previously described (Additional file 6: Table S3; [18]).

Abbreviations

CDD:

Conserved domain database

FC:

Fold change

GEO:

Gene expression omnibus

GO:

Gene ontology

References

  1. 1.

    Poulin R, Combes C. The concept of virulence: Interpretations and implications. Parasitol Today. 1999;15:474–5.

  2. 2.

    Tschirren B, Bischoff LL, Saladin V, Richner H. Host condition and host immunity affect parasite fitness in a bird-ectoparasite system. Funct Ecol. 2007;21:372–8.

  3. 3.

    Popara M, Villar M, Mateos-Hernández L, Fernández De Mera IG, De La Fuente J. Proteomics approach to the study of cattle tick adaptation to white tailed deer. Biomed Res Int. 2013;2013:319812.

  4. 4.

    González L, Carvajal J, George-Nascimento M. Differential infectivity of Caligus flexispina (Copepoda, Caligidae) in three farmed salmonids in Chile. Aquaculture. 2000;183:13–23.

  5. 5.

    Johnson SC. A comparison of development and growth rates of Lepeophtheirus salmonis (Copepoda: Caligidae) on naïve Atlantic (Salmo salar) and Chinook (Oncorhynchus tshawytscha) Salmon. In: Boxshall GA, Defaye D, editors. A comparison of development and growth rates of Lepeophtheirus salmonis (Copepoda: Caligidae) on naïve Atlantic (Salmo salar) and Chinook (Oncorhynchus tshawytscha) Salmon. Pathog: Wild Farmed Fish Sea Lice. Ellis Horwood; 1993. p. 68–82.

  6. 6.

    Fast MD, Burka JF, Johnson SC, Ross NW. Enzymes released from Lepeophtheirus salmonis in response to mucus from different salmonids. J Parasitol. 2003;89:7–13.

  7. 7.

    Torrissen O, Jones S, Asche F, Guttormsen A, Skilbrei OT, Nilsen F, et al. Salmon lice - impact on wild salmonids and salmon aquaculture. J Fish Dis. 2013;36:171–94.

  8. 8.

    Costello MJ. Ecology of sea lice parasitic on farmed and wild fish. Trends Parasitol. 2006;22:475–83.

  9. 9.

    Boxaspen K. A review of the biology and genetics of sea lice. ICES J Mar Sci. 2006;63:1304–16.

  10. 10.

    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.

  11. 11.

    Wagner GN, Fast MD, Johnson SC. Physiology and immunology of Lepeophtheirus salmonis infections of salmonids. Trends Parasitol. 2008;24:176–83.

  12. 12.

    Grimnes A, Jakobsen PJ. The physiological effects of salmon lice infection on post-smolt of Atlantic Salmon. Journal Fish Biol. 1996;48:1179–94.

  13. 13.

    Fast MD, Muise DM, Easy RE, Ross NW, Johnson SC. The effects of Lepeophtheirus salmonis infections on the stress response and immunological status of Atlantic Salmon (Salmo salar). Fish Shellfish Immunol. 2006;21:228–41.

  14. 14.

    Johnson SC, Albright LJ. Comparative susceptibility and histopathology of the response of naive Atlantic, Chinook and Coho Salmon to experimental infection with Lepeophtheirus salmonis (Copepoda : Caligidae). Dis Aquat Organ. 1992;14:179–93.

  15. 15.

    Braden LM, Koop BF, Jones SRM. Signatures of resistance to Lepeophtheirus salmonis include a Th2-type response at the louse-salmon interface. Dev Comp Immunol. 2015;48:178–91.

  16. 16.

    Braden LM, Barker DE, Koop BF, Jones SRM. Comparative defense-associated responses in salmon skin elicited by the ectoparasite Lepeophtheirus salmonis. Comp Biochem Physiol Part D Genomics Proteomics. 2012;7:100–9.

  17. 17.

    Sutherland BJG, Jantzen SG, Sanderson DS, Koop BF, Jones SRM. Differentiating size-dependent responses of juvenile Pink Salmon (Oncorhynchus gorbuscha) to sea lice (Lepeophtheirus salmonis) infections. Comp Biochem Physiol Part D Genomics Proteomics. 2011;6:213–23.

  18. 18.

    Sutherland BJG, Koczka KW, Yasuike M, Jantzen SG, Yazawa R, Koop BF, et al. Comparative transcriptomics of Atlantic Salmo salar, Chum Oncorhynchus keta and Pink Salmon O. gorbuscha during infections with salmon lice Lepeophtheirus salmonis. BMC Genomics. 2014;15:200.

  19. 19.

    Skugor S, Glover KA, Nilsen F, Krasnov A. Local and systemic gene expression responses of Atlantic Salmon (Salmo salar L.) to infection with the salmon louse (Lepeophtheirus salmonis). BMC Genomics. 2008;9:498.

  20. 20.

    Krasnov A, Skugor S, Todorcevic M, Glover KA, Nilsen F. Gene expression in Atlantic Salmon skin in response to infection with the parasitic copepod Lepeophtheirus salmonis, cortisol implant, and their combination. BMC Genomics. 2012;13:130.

  21. 21.

    Tadiso TM, Krasnov A, Skugor S, Afanasyev S, Hordvik I, Nilsen F. Gene expression analyses of immune responses in Atlantic Salmon during early stages of infection by salmon louse (Lepeophtheirus salmonis) revealed bi-phasic responses coinciding with the copepod-chalimus transition. BMC Genomics. 2011;12:141.

  22. 22.

    Jakob E, Sweeten T, Bennett W, Jones SRM. Development of the salmon louse Lepeophtheirus salmonis and its effects on juvenile Sockeye Salmon Oncorhynchus nerka. Dis Aquat Organ. 2013;106:217–27.

  23. 23.

    Devine GJ, Ingvarsdóttir A, Mordue W, Pike AW, Pickett J, Duce I, et al. Salmon lice, Lepeophtheirus salmonis, exhibit specific chemotactic responses to semiochemicals originating from the salmonid, Salmo salar. J Chem Ecol. 2000;26:1833–47.

  24. 24.

    Ingvarsdóttir A, Birkett M, Duce I, Genna RL, Mordue W, Pickett J, et al. Semiochemical strategies for sea louse control: Host location cues. Pest Manag Sci. 2002;58:537–45.

  25. 25.

    Mordue AJ, Birkett MA. A review of host finding behaviour in the parasitic sea louse, Lepeophtheirus salmonis (Caligidae: Copepoda). J Fish Dis. 2009;32:3–13.

  26. 26.

    Bailey RJE, Birkett M, Ingvarsdóttir A, Luntz JM, Mordue W, Shea BO, et al. The role of semiochemicals in host location and non-host avoidance by salmon louse (Lepeophtheirus salmonis) copepodids. Can J Fish Aquat Sci. 2006;63:448–56.

  27. 27.

    Fast MD, Ross NW, Craft C, Locke SJ, MacKinnon SL, Johnson SC. Lepeophtheirus salmonis: characterization of prostaglandin E(2) in secretory products of the salmon louse by RP-HPLC and mass spectrometry. Exp Parasitol. 2004;107:5–13.

  28. 28.

    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.

  29. 29.

    Marchler-Bauer A, Derbyshire MK, Gonzales NR, Lu S, Chitsaz F, Geer LY, et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2015;43:D222–6.

  30. 30.

    Poley JD, Sutherland BJ, Jones SR, Koop BF, Fast MD. Sex-biased gene expression and sequence conservation in Pacific and Atlantic salmon Lice (Lepeophtheirus salmonis). BMC Genomics. 2016;17:1–16.

  31. 31.

    Raeymaekers JM, Wegner KM, Huyse T, Volckaert FM. Infection dynamics of the monogenean parasite Gyrodactylus gasterostei on sympatric and allopatric populations of the Three-spined Stickleback Gasterosteus aculeatus. Folia Parasitol. 2011;58:27–34.

  32. 32.

    Brandal PO, Egidius E, Romslo I. Host blood: a major food component for the parasitic copepod Lepeophtheirus salmonis Kroyeri, 1838 (Crustacea: Caligidae). Norwegian J Zool. 1976;24:341–3.

  33. 33.

    Abdul Alim M, Tsuji N, Miyoshi T, Khyrul Islam M, Huang X, Motobu M, et al. Characterization of asparaginyl endopeptidase, legumain induced by blood feeding in the ixodid tick Haemaphysalis longicornis. Insect Biochem Mol Biol. 2007;37:911–22.

  34. 34.

    Williamson AL, Brindley PJ, Knox DP, Hotez PJ, Loukas A. Digestive proteases of blood-feeding nematodes. Trends Parasitol. 2003;19:417–23.

  35. 35.

    Motobu M, Tsuji N, Miyoshi T, Huang X, Islam MK, Alim M, et al. Molecular characterization of a blood-induced serine carboxypeptidase from the ixodid tick Haemaphysalis longicornis. FEBS J. 2007;274:3299–312.

  36. 36.

    Manship BM, Walker AJ, Jones LA, Davies AJ. Blood feeding in juvenile Paragnathia formica (Isopoda: Gnathiidae): biochemical characterization of trypsin inhibitors, detection of anticoagulants, and molecular identification of fish hosts. Parasitology. 2012;139:744–54.

  37. 37.

    Zhu K, Dillwith JW, Bowman AS, Sauer JR. Identification of hemolytic activity in saliva of the lone star tick (Acari:Ixodidae). J Med Entomol. 1997;34:160–6.

  38. 38.

    Cunningham E, McCarthy E, Copley L, Jackson D, Johnson D, Dalton JP, et al. Characterisation of cathepsin B-like cysteine protease of Lepeophtheirus salmonis. Aquaculture. 2010;310:38–42.

  39. 39.

    McCarthy E, Cunningham E, Copley L, Jackson D, Johnston D, Dalton JP, et al. Cathepsin L proteases of the parasitic copepod, Lepeophtheirus salmonis. Aquaculture. 2012;356–357:264–71.

  40. 40.

    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.

  41. 41.

    Ahammed Shareef PA, Abidi SMA. Cysteine protease is a major component in the excretory/secretory products of Euclinostomum heterostomum (Digenea: Clinostomidae). Parasitol Res. 2013;113:65–71.

  42. 42.

    Robinson MW, Dalton JP, Donnelly S. Helminth pathogen cathepsin proteases: it’s a family affair. Trends Biochem Sci. 2008;33:601–8.

  43. 43.

    Maldonado-Aguayo W, Chávez-Mardones J, Gonçalves AT, Gallardo-Escárate C. Cathepsin gene family reveals transcriptome patterns related to the infective stages of the salmon louse Caligus rogercresseyi. PLoS One. 2015;10, e0123954.

  44. 44.

    Turk V, Stoka V, Vasiljeva O, Renko M, Sun T, Turk B, et al. Cysteine cathepsins: From structure, function and regulation to new frontiers. Biochim Biophys Acta Proteins Proteomics. 1824;2012:68–88.

  45. 45.

    Barbara G, Giorgio De R, Stanghellini V, Corinaldesi R, Cremon C, Gerard N, et al. Neutral endopeptidase (EC 3.4.24.11) downregulates the onset of intestinal inflammation in the nematode infected mouse. Gut. 2003;52:1457–64.

  46. 46.

    Sakharov DV, Plow EF, Rijken DC. On the mechanism of the antifibrinolytic activity of plasma carboxypeptidase B. J Biol Chem. 1997;272:14477–82.

  47. 47.

    Yan J, Cheng Q, Li CB, Aksoy S. Molecular characterization of three gut genes from Glossina morsitans morsitans: cathepsin B, zinc-metalloprotease and zinc-carboxypeptidase. Insect Mol Biol. 2002;11:57–65.

  48. 48.

    Song JJ, Hwang I, Cho KH, Garcia M, Kim AJ, Wang TH, et al. Plasma carboxypeptidase B downregulates inflammatory responses in autoimmune arthritis. J Clin Invest. 2011;121:3517–27.

  49. 49.

    Fry BG, Roelants K, Champagne DE, Scheib H, Tyndall JD, King GF, et al. The toxicogenomic multiverse: convergent recruitment of proteins into animal venoms. Annu Rev Genomics Hum Genet. 2009;10:483–511.

  50. 50.

    Palm NW, Rosenstein RK, Yu S, Schenten DD, Florsheim E, Medzhitov R. Bee venom phospholipase A2 induces a primary type 2 response that is dependent on the receptor ST2 and confers protective immunity. Immunity. 2013;39:976–85.

  51. 51.

    Izidoro LFM, Sobrinho JC, Mendes MM, Costa TR, Grabner AN, Rodrigues VM, et al. Snake venom L-amino acid oxidases: Trends in pharmacology and biochemistry. Biomed Res Int. 2014;2014.

  52. 52.

    Espino AM, Hillyer GV. Molecular cloning of a member of the Fasciola hepatica saposin-like protein family. J. Parasitol. 2003;89:545–52.

  53. 53.

    Don TA, Bethony JM, Loukas A. Saposin-like proteins are expressed in the gastrodermis of Schistosoma mansoni and are immunogenic in natural infections. Int J Infect Dis. 2008;12:39–47.

  54. 54.

    Winkelmann J, Leippe M, Bruhn H. A novel saposin-like protein of Entamoeba histolytica with membrane-fusogenic activity. Mol Biochem Parasitol. 2006;147:85–94.

  55. 55.

    Radulović ZM, Kim TK, Porter LM, Sze S-H, Lewis L, Mulenga A. A 24–48 h fed Amblyomma americanum tick saliva immuno-proteome. BMC Genomics. 2014;15:518.

  56. 56.

    Bruhn H. A short guided tour through functional and structural features of saposin-like proteins. Biochem J. 2005;389:249–57.

  57. 57.

    Robinson NJ, Winge DR. Copper metallochaperones. Annu Rev Biochem. 2010;79:537–62.

  58. 58.

    Burmester T. Evolutionary history and diversity of arthropod hemocyanins. Micron. 2004;35:121–2.

  59. 59.

    Ahearn GA, Mandal PK, Mandal A. Mechanisms of heavy-metal sequestration and detoxification in crustaceans: A review. J Comp Physiol B Biochem Syst Environ Physiol. 2004;174:439–52.

  60. 60.

    Fast MD, Ross NW, Mustafa A, Sims DE, Johnson SC, Conboy GA, et al. Susceptibility of Rainbow Trout Oncorhynchus mykiss, Atlantic Salmon Salmo salar and Coho Salmon Oncorhynchus kisutch to experimental infection with sea lice Lepeophtheirus salmonis. Dis Aquat Organ. 2002;52:57–68.

  61. 61.

    Fast MD, Sims DE, Burka JF, Mustafa A, Ross NW. Skin morphology and humoral non-specific defence parameters of mucus and plasma in Rainbow Trout, Coho and Atlantic Salmon. Comp Biochem Physiol A Mol Integr Physiol. 2002;132:645–57.

  62. 62.

    Muona M, Soivio A. Changes in plasma lysozyme and blood leucocyte levels of hatchery-reared Atlantic Salmon (Salmo salar L.) and Sea Trout (Salmo trutta L.) during parr-smolt transformation. Aquaculture. 1992;106:75–87.

  63. 63.

    Yu JN, Azuma N, Abe S. Genetic differentiation between collections of hatchery and wild Masu Salmon (Oncorhynchus masou) inferred from mitochondrial and microsatellite DNA analyses. Environ Biol Fishes. 2012;94:259–71.

  64. 64.

    Wang L, Shi X, Su Y, Meng Z, Lin H. Loss of genetic diversity in the cultured stocks of the Large Yellow Croaker, Larimichthys crocea, revealed by microsatellites. Int J Mol Sci. 2012;13:5584–97.

  65. 65.

    Yazawa R, Yasuike M, Leong J, Von Schalburg KR, Cooper G, 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.

  66. 66.

    Price PW. Evolutionary Biology of Parasites. Princeton: Princeton University Press; 1980.

  67. 67.

    Pecquerie L, Johnson LR, Kooijman SALM, Nisbet RM. Analyzing variations in life-history traits of Pacific salmon in the context of Dynamic Energy Budget (DEB) theory. J Sea Res. 2011;66:424–33.

  68. 68.

    Norris K, Evans MR. Ecological immunology : life history trade-offs and immune defense in birds. Behav Ecol. 2000;11:19–26.

  69. 69.

    Zuk M, Stoehr AM. Immune defense and host life history. Am Nat. 2002;160(Suppl):S9–22.

  70. 70.

    Braden LM, Barker DE, Koop BF, Jones SRM. Differential modulation of resistance biomarkers in skin of juvenile and mature Pink Salmon, Oncorhynchus gorbuscha by the salmon louse, Lepeophtheirus salmonis. Fish Shellfish Immunol. 2015;47:7–14.

  71. 71.

    Jones SRM, Fast MD, Johnson SC, Groman DB. Differential rejection of salmon lice by Pink and Chum Salmon: Disease consequences and expression of proinflammatory genes. Dis Aquat Organ. 2007;75:229–38.

  72. 72.

    Jones S, Kim E, Bennett W. Early development of resistance to the salmon louse, Lepeophtheirus salmonis (Krøyer), in juvenile Pink Salmon, Oncorhynchus gorbuscha (Walbaum). J Fish Dis. 2008;31:591–600.

  73. 73.

    Nagasawa K. Annual changes in the population size of the salmon louse Lepeophtheirus salmonis (Copepoda: Caligidae) on high-seas Pacific salmon (Oncorhynchus spp.), and relationship to host abundance. Hydrobiologia. 2001;453–454:411–6.

  74. 74.

    Bize P, Jeanneret C, Klopfenstein A, Roulin A. What makes a host profitable? Parasites balance host nutritive resources against immunity. Am Nat. 2008;171:107–18.

  75. 75.

    Heylen DJ, Matthysen E. Experimental evidence for host preference in a tick parasitizing songbird nestlings. Oikos. 2011;120:1209–16.

  76. 76.

    Heylen DJ, White J, Elst J, Jacobs I, Van De Sande C, Matthysen E. Nestling development and the timing of tick attachments. Parasitology. 2012;139:766–73.

  77. 77.

    Johnson SC, Blaylock RB, Elphick J, Hyatt KD. Disease induced by the sea louse (Lepeophtheirus salmonis) (Copepoda: Caligidae) in wild Sockeye Salmon (Oncorhynchus nerka) stocks of Alberni Inlet, British Columbia. Can J Fish Aquat Sci. 1996;53:2888–97.

  78. 78.

    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.

  79. 79.

    Huang DW, Sherman BT, Lempicki R. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

  80. 80.

    Oliveros JC. VENNY. An interactive tool for comparing lists with Venn Diagrams. 2007.

  81. 81.

    Jantzen SG, Sutherland BJG, Minkley DR, Koop BF. GO Trimming: Systematically reducing redundancy in large Gene Ontology datasets. BMC Res Notes. 2011;4:267.

  82. 82.

    Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000;132:365–86.

  83. 83.

    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.

Download references

Acknowledgements

The authors acknowledge Holly Hicklin and Robert Kennedy for assistance with fish culture. Coho Salmon were a generous gift from Dr. Chris Foote and the Chase River Hatchery, and Sockeye Salmon were a generous gift from the Inch Creek Hatchery. The authors thank Marine Harvest Canada for assisting with parasite collection. The authors thank D. Minkley and J. Poley for assistance with bioinformatics analysis. T. Hori and J. Poley are gratefully acknowledged for comments on earlier versions of this manuscript. The authors are thankful for valuable comments from two anonymous reviewers.

Funding

This research was funded by Genome BC, the Province of BC, the Department of Fisheries and Oceans Canada (DFO), and the University of Victoria. LMB and BJGS were funded by NSERC-CGS fellowships during this work.

Availability of data and materials

The datasets supporting the results of this article are available in the Gene Expression Omnibus repository GSE80220.

The temporary reviewer access link for Pacific lice (GSE80220): http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE80220

Authors’ contributions

Conceived and designed the experiments: LMB, BJGS, SRMJ. Performed the experiments: LMB. Analyzed the data: LMB, BJGS, BFK, SRMJ. Contributed reagents/materials/analysis tools: LMB, BJGS, BFK, SRMJ. Wrote the manuscript: LMB, BJGS, BFK, SRMJ. All authors have read and approved the manuscript.

Authors’ information

Not applicable.

Competing interests

The authors declare they have no competing interests.

Consent for publication

Not applicable.

Ethics approval

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 14–001.

Author information

Affiliations

  1. Centre for Biomedical Research, University of Victoria, Victoria, British Columbia, Canada
    • Laura M. Braden
    • , Ben J. G. Sutherland
    •  & Ben F. Koop
  2. Pacific Biological Station, Fisheries & Oceans Canada, Nanaimo, British Columbia, Canada
    • Simon R. M. Jones
  3. Present Address: Department of Pathology and Microbiology, Atlantic Veterinary College, Charlottetown, Prince Edward Island, Canada
    • Laura M. Braden
  4. Present Address: Département de biologie, Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Québec, Québec, Canada
    • Ben J. G. Sutherland
Authors
  1. Search for Laura M. Braden in:
  2. Search for Ben J. G. Sutherland in:
  3. Search for Ben F. Koop in:
  4. Search for Simon R. M. Jones in:

Corresponding author

Correspondence to Simon R. M. Jones.

Additional files

Additional file 1: Table S1.

Numbers of differentially expressed genes detected in L. salmonis while feeding on different host species (i) compared to starved lice, or (ii) compared to each other. (XLSX 25 kb)

Additional file 2:

Transcriptomic response of L. salmonis while feeding on Atlantic, Coho or Sockeye Salmon. (XLSX 429 kb)

Additional file 3: Figure S1.

Profiling the starvation response of L. salmonis. Overexpressed transcripts in lice withheld from hosts were compared to Atlantic-, Coho- or Sockeye-fed lice to produce a list of unique genes involved in the “starvation response” of L. salmonis. These genes were analyzed using DAVID to produce enriched gene lists after (A) 24 and (B) 48 hpi. aFold Enrichment. (TIF 14475 kb)

Additional file 4: Table S2.

Gene Ontology enrichment of the host-specific feeding response in L. salmonis. Functional enrichment of DEGs from Atlantic-fed compared to either Coho-fed or Sockeye-fed lice at 24 and 48 hpi, and from either Coho-fed or Sockeye-fed lice compared to Atlantic-fed lice at 24 hpi. (XLSX 11 kb)

Additional file 5: Figure S2.

The salmon louse response cannot be explained by acclimation to the host. A reciprocal experiment was conducted by placing lice originally collected from Atlantic (at) or Sockeye Salmon (sk) on either Atlantic Salmon (AT) or Sockeye Salmon (SK) hosts. Expression of protein synthesis (A), energy metabolism (B) and digestion (C) was significantly higher in lice feeding on Atlantic Salmon (blue boxplots) irrespective of their original host, thus negating the acclimation hypothesis. Differentially expressed transcripts identified by the microarray were profiled using RT-qPCR, and are shown as log2 calibrated normalized relative quantities (CNRQ). Significance was identified by two-way ANOVA (p < 0.05) followed by post-hoc Tukey test to determine pairwise significance. Differences between groups are denoted by lower case letters, while differences over time within a group is denoted by an asterisk (*). (TIF 17046 kb)

Additional file 6: Table S3.

Validation of microarray experiments. The expression profiles of a select subset of genes were compared between RT-qPCR and microarray results. Expression levels were validated by testing the correlation of log2 expression values for RT–qPCR samples against microarray log2 expression ratios (Cy5/Cy3). (XLSX 75 kb)

Additional file 7: Table S4.

List of primer sequences for RT-qPCR. (XLSX 42 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Braden, L.M., Sutherland, B.J.G., Koop, B.F. et al. 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 18, 110 (2017). https://doi.org/10.1186/s12864-017-3520-1

Download citation

Keywords

  • Copepoda
  • Host susceptibility
  • Lepeophtheirus salmonis
  • Resistance
  • Salmon
  • Sea lice
  • Virulence factors

BMC Genomics

ISSN: 1471-2164

Contact us