Skip to main content

Roles of three putative salmon louse (Lepeophtheirus salmonis) prostaglandin E2 synthases in physiology and host–parasite interactions



The salmon louse (Lepeophtheirus salmonis) is a parasite of salmonid fish. Atlantic salmon (Salmo salar) exhibit only a limited and ineffective immune response when infested with this parasite. Prostaglandins (PGs) have many biological functions in both invertebrates and vertebrates, one of which is the regulation of immune responses. This has led to the suggestion that prostaglandin E2 (PGE2) is important in the salmon louse host–parasite interaction, although studies of a salmon louse prostaglandin E2 synthase (PGES) 2 gene have not enabled conformation of this hypothesis. The aim of the present study was, therefore, to characterize two additional PGES-like genes.


Lepeophtheirus salmonis microsomal glutathione S-transferase 1 like (LsMGST1L) and LsPGES3L were investigated by sequencing, phylogenetics, transcript localization and expression studies. Moreover, the function of these putative PGES genes in addition to the previously identified LsPGES2 gene was analyzed in double stranded (ds) RNA-mediated knockdown (KD) salmon louse.


Analysis of the three putative LsPGES genes showed a rather constitutive transcript level throughout development from nauplius to the adult stages, and in a range of tissues, with the highest levels in the ovaries or gut. DsRNA-mediated KD of these transcripts did not produce any characteristic changes in phenotype, and KD animals displayed a normal reproductive output. The ability of the parasite to infect or modulate the immune response of the host fish was also not affected by KD.


Salmon louse prostaglandins may play endogenous roles in the management of reproduction and oxidative stress and may be a product of salmon louse blood digestions.

Graphic Abstract


The salmon louse (Lepeophtheirus salmonis) is a marine ectoparasite of salmonid fish, feeding on fish mucus, skin and blood [1, 2]. Salmon lice have a life-cycle consisting of eight stages, starting with three planktonic larval stages: nauplius 1, nauplius 2 and the infective copepodid stage. All further development takes place on the fish and consists of two chalimi, two preadult stages and the adult stage [3, 4]. Heavy infestations can cause severe problems to the host, including wounding, secondary infections and osmotic disturbances [5, 6]. Infestations of farmed fish have been difficult to manage due to the development of resistance against parasiticides [7], and the negative impact on wild salmonid fish through transmission of lice from farmed to wild fish stocks is a cause of environmental concern [8,9,10].

Prostaglandin E2 (PGE2) is an eicosanoid component with hormone-like functions, derived from arachidonic acid (AA) through three sequential enzymatic reactions [11]. Generally, phospholipase A2 releases AA from cell membranes, a cyclooxygenase (COX) converts AA to PGG2 and further to PGH2, which can be metabolized to PGE2 by PGE2 synthases (PGESs). The action of PGE2 is dependent on the binding of prostanoid EP receptor subtypes (EP1–4), and PGE2 has diverse roles in physiological processes, such as inflammation, reproduction and embryonic development, depending on cell type and EP receptor subtype [12,13,14,15,16].

PGESs are found in a large range of organisms where they play diverse physiological roles [17]. In mammals, three isotypes of PGES exist, two of which are membrane bound (PGES1 and -2) and one cytosolic (PGES3). These isotypes typically display different enzymatic properties, expression patterns, cellular localizations and functions. PGESs have also been described in many arthropods, including crustacean species, and similar to mammals, three genes encoding PGESs have been found in penaeid shrimp [18]. The functional property of a few invertebrate synthases has also been confirmed, where the conversion of PGH2 to PGE2 by recombinant Gammarus and Caprella sp. PGES2 has been demonstrated [19]. PGESs and their product PGE2 are involved in reproduction, including maturation of ovaries and vitellogenesis in prawn and crabs [18, 20,21,22]. PGES has also been detected in blood cells of crab, indicative of a role in immunity [23], and PGE2 has been detected in both tapeworms and crab nervous tissue where it is likely to be important for development and function [20, 24].

Parasites with a long-term interaction with their host need mechanisms to evade potential host immune reaction, and prostaglandins (PGs) such as PGE2 have been found to play important roles in many host–parasite protozoan and metazoan interactions [25, 26]. However, uncovering the role of PGs is complicated by their diverse roles and the fact that these compounds are produced by both the host and the parasite [27]. In blood-feeding arthropods such as ticks, PGs are suggested to be involved in vasodilation and immune modulation in the host, as high concentrations of PGE2 are found in tick saliva [28,29,30]. Both PGE2 and PGD2 have been detected in gland secretions from the fish tapeworm Diphyllobothrium dendriticum when incubated in serum from a host fish [24], and in vitro studies have suggested that PGE2 derived from the tapeworm Schistocephalus solidus have the potential to modulate host leukocyte viability and reactive oxygen species production [31].

Atlantic salmon (Salmo salar) display a limited and mainly local immune response to salmon louse, which does not eliminate the infection [32,33,34,35]. This has led to efforts to understand how the salmon louse modulates its host, and PGE2 has been put forward as a putative candidate for immune modulation. PGE2 has been detected in excretory/secretory products from dopamine-treated salmon louse, and it was demonstrated that the level of PGE2 quickly dropped after lice were removed from the host fish [36, 37]. Moreover, Atlantic salmon immune tissue express the EP4 prostaglandin E receptors [38], and PGE2 treatment has been shown to down-modulate the expression of MHC class I proteins and interleukin-1β in a salmon macrophage-like cell line [37]. However, in a more recent study in which excretory/secretory products were obtained by agitating lice in sea water, PGE2 was not identified [39]. Moreover, studies of a putative salmon louse gene, PGES2 (LsPGES2) were unable to reveal any involvement of LsPGES2 in the host–parasite interaction [40]. Instead, expression of LsPGES2 indicated a function in reproductive organs, although no visible alteration in reproductive phenotype was observed in knockdown (KD) animals.

In the present study, we hypothesized that this lack of phenotypic alterations is due to the action of additional salmon lice PGES genes beyond the already described LsPGES2 and that these may be involved in the host–parasite interaction in addition to endogenous processes. Two additional genes with putative PGES activity were identified and their functional role analyzed by double-stranded (ds) RNA-mediated KD and studies of host–parasite interactions, reproduction, oxidative stress and blood digestion.


Culture of salmon lice and source of tissue

Laboratory strains of the Atlantic salmon lice subspecies Lepeophtheirus salmonis salmonis [41] were maintained on farmed Atlantic salmon (Salmo salar) according to Hamre et al. [42]. The salmon were hand fed a commercial diet and reared in sea water with a salinity of 34.5 g/kg and a temperature of 10 °C. Eggs, nauplii and copepodids were kept in seawater from the same source. Nauplii were obtained from hatching eggs and kept in a flow-through single-well system [42].

For ontogenetic expression analysis of the two PGES-like genes, LsMGST1L and LsPGES3L, we collected a time-series of all lice stages from egg to adult lice, in pentaplicate samples, as follows: (i) eggs: 1 egg sac (string) (containing approximately 200 eggs; (ii) nauplius 1, nauplius 2 and copepodids (free-living) (approximately 100 larvae); copepodids, 2 days post-infestation (dpi) and 4 dpi (60 larvae); (iv) chalimus: chalimus 1 (30 animals) and chalimus 2 (20 animals); (v) preadult and adult stages: single animals (adult females were defined as young before the first egg extrusion had taken place and as mature when producing eggs).

To analyze transcripts during the reproductive cycle, adult females were collected in RNAlater® (Qiagen, Hilden, Germany) while their egg strings were removed and kept in a flow-through single-well system at 10 °C [42]. Adult female lice continuously produce eggs that are deposited in batches approximately every 10 days at 10 °C. After deposition of the eggs carried externally by the female until hatching, a new batch of eggs (300–700) is matured and undergoes vitellogenesis in the genital segment [43]. To determine where in the reproductive cycle the sampled females were at sampling, the hatching day of the egg strings after sampling was recorded. The sampled females were further grouped according to days prior to hatching of their egg strings.

To analyze PGE2 levels in adult females in relation to blood-feeding, adult female lice were removed from a host fish and either frozen directly on dry ice or kept in saltwater in flow-through single-well incubators [42] for 3 h after removal from the host at 10 °C before sampling on dry ice. The animals were kept at − 80 °C until further analysis.

RNA purification and cDNA synthesis

All samples for RNA isolation were collected in RNAlater® (LifeTechnologies, Carlsbad, CA, USA), kept at 4 °C overnight and stored at − 20 °C. Total RNA was isolated with a combined TRI Reagent® (Sigma-Aldrich, St. Louis, MO, USA) and RNeasy (Qiagen) method, as previously described [44]. Extracted total RNA was either frozen at − 80 °C or used directly for cDNA synthesis. For real-time reverse transcriptase (RT)-PCR, cDNA synthesis was carried out using the AffinityScript qPCR cDNA Synthesis Kit (Stratagene, San Diego, CA, USA) with oligo(dT) as primers according to the supplier recommendations, adding 200 ng lice total RNA or 1000 ng salmon total RNA. The cDNA samples were diluted five times and stored at − 20 °C until use.

To produce a template for PCR, the qScript cDNA SuperMix (Quanta BioSciences, Gaithersburg, MD, USA) was used, applying 1 µg total RNA from adult female salmon louse.

PCR, cloning and sequencing

Candidate PGES genes were identified within the salmon louse genome by BLAST search in LiceBase ( Gene-specific primers were designed (Table 1), and rapid amplification of cDNA ends (RACE; (5′ and 3′) was performed using the SMARTer™ RACE cDNA amplification kit (Clontech Laboratories, Inc. Mountain View, CA, USA) using 1 μg total RNA from adult female lice.

Table 1 Primers used for rapid amplification of cDNA ends, in situ hybridization, double-stranded RNA synthesis for knockdown studies and quantitative PCR of the target gene and reference gene for Lepeophtheirus salmonis and Salmo salar

The RACE products were further cloned using the TOPO TA Cloning® Kit for sequencing (Invitrogen, Carlsbad, CA, USA). Colonies were used as templates for the PCR reaction with M13 forward and reverse primers, and the PCR products were further purified by ExoSAP-it (Affymetrix; Thermo Fisher Scientific, Waltham, MA, USA) prior to sequencing. The purified PCR products were sequenced using the BigDye Terminator v3.1 Cycle Sequencing Kit from Applied Biosystems and run by an ABI PRISM 7700 automated sequencing detecting system (Applied Biosystems, Foster City, CA, USA) at the University of Bergen sequencing facility.

Bioinformatic analysis

Sequences were assembled and translated using Vector NTI Advance 10 software (Invitrogen). The open reading frame (ORF) was blasted using ExPASy BLAST form ( and NCBI blast form (, and aligned in Clustal Omega ( Location of domains was predicted by InterPro ( Searches for orthologous sequences to PGE synthases were conducted using BlastP in the predicted L. salmonis proteome in Ensembl Metazoa (release 48) and TBlastN against the representative LSalAtl2s genome assembly ( Blast searches were executed both on the command line under Linux (NCBI Blast suite version 2.9.0+), Ensembl Metazoa, and through using LsPGES2 (EMLSAP00000000441) and human and Drosophila melanogaster mGST and PGES-like protein sequences as queries (UniProtKB accessions: Q9H7Z7, O14684, Q15185, B3KPZ2, Q8SY19, Q9V420) with a Blast E-value < 1e−6. In addition, we searched Ensembl (release 101) and Ensembl Metazoa for precomputed orthologues of D. melanogaster and L. salmonis having GO:0050220 (PGES activity) as Gene Ontology (GO) term annotation or InterPro annotation IPR023352 (Membrane associated eicosanoid/glutathione metabolism-like domain superfamily).

The expression of the three putative PGES was also compared to high-throughput data from previous investigations. For tissue-specific analysis, the expression measured by microarray of five different tissues (intestine female and male, ovary, testis, subcuticular tissue and frontal neuronal and gland enriched tissue) [45] was used. For each of the three genes, one oligo was shown (EMLSAT00000000441: CUST_3823_PI425513912; EMLSAT00000006733: CUST_21227_PI425553006; EMLSAT00000012943: CUST_11144_PI425553006), as identified in LiceBase GBrowse. For investigation of expression changes during development with respect to molting, RNA sequencing data from a series of chalimus and preadult 1 larvae of different instar age with respect to molting were used [46]. In this study, chalimus 1, 2 and preadult 1 salmon louse larvae were divided by size measurements into different sex and age groups (young, middle, old and molt, according to their instar age). The young group included lice that had recently molted to the respective stage; the middle group included lice approximately in the middle of the stage; the old group contained lice that were about to molt to the next stage; and the molt group contained the lice showing morphological signs of molting.

Phylogenetic analysis

Orthologous amino acid sequences for EMLSAP00000006733 and EMLSAP00000012943 were obtained by Blast searches versus UniprotKB and GeneBank [47, 48]. Sequence alignments and accessions used in the analysis are provided in Additional files 1 and 2. Multiple sequence alignments (MSA) were computed using MUSCLE via the MUSCLE web-service with default parameter [49]. MSAs were inspected and manually end-clipped to the first and last conserved column using Jalview [50]; no further editing of MSAs was applied. Phylogenetic trees were constructed using the Markov Chain Monte Carlo method (MCMC) in MrBayes (v3.2.6) with 5 million generations, eight chains and four rate categories for the Gamma distribution [51]. RAxML (v8.2.9) was used for maximum likelihood (ML) analysis with 1000 bootstrap iterations [52]. Substitution models for ML analysis were determined by prottest3 for each set of orthologous sequences based on the two best fit model under the Akaike information criterion (AIC) [53]; these were LG + I + Γ + F and LG + Γ + F for both groups of orthologues. MrBayes was first applied in mixed mode, and the resulting best fit model was selected for a second MCMC run, always allowing for invariant sites, resulting in the WAG + I + Γ model for both groups. Phylograms were inspected and rendered in FigTree v1.4.4 ( and figures were edited in Inkscape.

In situ hybridization

Salmon lice were fixed in 4% paraformaldehyde in 4% phosphate buffer (pH 7.4) for 24 h at 4 °C, and the specimens were then processed in the Reichert-Jung Histokinette 2000 rotary tissue processor (Leica Microsystems GmbH, Wetzlar, Germany) where they were washed in phosphate-buffered saline, dehydrated through a graded ethanol series and embedded in paraffin wax. Sections, 4.0 μm thick, were cut with a Leica RM 225 microtome (Leica Microsystems). Digoxigenin (DIG)-labeled antisense and sense RNA probes of the two putative PGES were prepared by in vitro transcription using the DIG RNA Labeling Kit (Roche AG, Basel, Switzerland) on purified PCR products (Table 1) that included T7 promoters in the templates. In situ hybridization was performed according to [54], with some modifications as earlier described [55]. Hybridizations with sense probes were carried out as negative control.

Real-time PCR

Real-time PCR was performed and the products analyzed as described below using the established salmon louse elongation factor 1 alfa (eEF1α) and adenine nucleotide translocator 3 (ADT3) gene as reference genes when analyzing developmental stages and animals on which RNA interference (RNAi) had been performed [56, 57]. In animals where RNAi was performed targeting LsHPX1 (see next section), 18S was used as reference gene as eEF1α and ADT3 previously have been found to be affected by LsHPX1 KD [57]. Primers used for real-time RT-PCR are listed in Table 1. Real-time PCR was performed with the 1× PowerUp™ SYBR Green Master Mix (Thermo Fisher Scientific), 500 nM forward and reverse primers and 2 µl diluted cDNA in 10-µl reaction volumes. Samples were run in duplicate on the Applied Biosystems 7500 Real-Time PCR System under standard conditions (50 °C, 2 min; then 95 °C/2 min; 95 °C/15 s and 60 °C/1 min for 40 cycles; followed by a melt curve analysis at 60–95 °C. A five-point standard curve of fourfold dilutions was made for each assay to calculate PCR efficiencies, given by the equation E% = (101/slope − 1) × 100 [58]. The relative differences in quantification cycle (Cq) between the target gene and the reference genes (ΔCq) and expression relative to a calibrator (ΔΔCq) were calculated, transformed by the equation 2−ΔΔCq [59]. T-tests were used to determine differently expressed genes with a threshold P value of 0.05.

RNA interference

Genes selected for KD were analyzed to identify sequences suitable for specific KD by treatment with dsRNA. Primers were chosen within regions that showed little similarity to other salmon louse genes. Moreover, the chosen fragment was aligned against the salmon louse genome by BlastN to exclude similarity between DICER-produced 19-mers and any salmon louse gene. RNAi was performed as previously described in nauplius [60]. In short, long dsRNA for the CPY control gene (cod trypsin) and PGES genes were produced using the Ambion MEGAscript® RNAi Kit (Thermo Fisher Scientific) according to supplier’s instructions using the primers listed in Table 1. Additionally, dsRNA targeting LsHPX1 transcripts was performed to create animals with increased levels of oxidative stress [57]. For soaking of nauplius, pools of 20–60 nauplius 1 larvae from the same egg string were incubated in 10 ng/µl dsRNA for each gene and incubated overnight (17 h). Thereafter, all animals from the same treatment group were pooled and transferred to a flow-through well, where they were kept until sampling.

For RNAi in preadult 2 females, a modified protocol based on [61] was used. Before injection of dsRNA (600 ng/ul), 1 μl of saturated and filtered bromophenol blue was added to 50 μl dsRNA solution; this concentration of bromophenol blue has been found to be well tolerated by the salmon louse [61]. The needle was placed dorsally into the haemocoel of the cephalothorax close to the ovaries. The custom-made needles were pulled by utilizing a 1-mm Borosilicate glass tube with an inner diameter of 0.5 mm (Sutter Instrument, Novato, CA, USA) on the P-2000 laser-based micropipette puller system (Sutter Instrument). The dsRNA solution was added to the needle using a microloader tip and then coupled to a HI-7 injection holder (Narishige Ltd., London, UK). The tip of the needle was removed with a scalpel blade. By blowing air into the needle, less than 1 µl solution was injected and quickly spread in the louse, as observed by the dispersal of blue color. After injection and a brief 2-h period in seawater, the lice were placed back onto fish to observe their development to adult individuals and measure their reproductive output. The gene expression was analyzed in control and KD animals.

Photographs were taken of lice from the RNAi experiments using a Canon EOS 600D camera (Canon Inc., Tokyo, Japan) mounted with an adaptor (LMscope) to an Olympus SZX9 dissecting microscope (Olympus Corp., Tokyo, Japan).

Infection trials

Copepodites treated with dsRNA targeting control fragments or a combination of the putative PGES genes were used to infect fish that were kept in single tanks (six fish/treatment group). Lice sampled to analyze KD by real-time RT-PCR were harvested at 2 and 14 dpi of fish, when the lice were still copepodites or had developed into the preadult 1 stage, respectively. From the fish, two different sample types of salmon skin, including the dermis, epidermis and scales, were collected on the same days: one sample was taken directly under a louse and one from nearby non-parasitized skin. Skin was obtained only from areas of the fish with scales since distinct expression of Atlantic salmon immune genes after salmon lice infestation has been observed in scaled skin compared to scaleless skin on the head region [62]. Moreover, skin at or above the lateral line was avoided to ensure more homogenous samples among individuals. The immune response in skin samples was analyzed with the primers listed in Table 1.

Prostaglandin E2 and F analysis

The Prostaglandin E Metabolite EIA Kit and the 13,14-dihydro-15-keto Prostaglandin F2α ELISA Kit (Cayman Chemical Company, Ann Arbor, MI, USA) were used to measure the levels of PGE2 and PGF metabolites in LsHPX1 KD copepodids and adult lice homogenates, according to supplier’s instructions with some modification as follows. For the PGE2 analysis, around 200 copepodids were pooled or one female was homogenized in 75 µl ddH2O using 1.4-mm zirconium oxide beads (Precellys 24; Thermo Fisher Scientific) and a TissueLyser LT (Qiagen) for 4 min at 50 Hz. Samples were centrifuged for 5 min at 1.5 g, and 55 µl of the supernatant was added to 16.5 µl carbonate buffer. After an overnight incubation on 37 °C, 22 µl of phosphate buffer and 16.5 µl EIA buffer were added, and 50 µl of the sample was added to two wells (two technical replicates). For the PGF analysis, around 150 copepodids were pooled in 100 µl ddH2O (N = 3) and homogenized as described above. Samples were centrifuged for 5 min at 1.5 g, and 50 µl of the supernatant was added to the ELISA wells. Two technical replicates were run for each of the three biological replicates.


Sequence analysis

Searches in the salmon louse genome ( revealed two putative PGE synthases in addition to the already published LsPGES2 [40]. Based on similarity to other genes, the putative PGES were named Lepeophtheirus salmonis microsomal glutathione S-transferase 1 like (LsMGST1L) and prostaglandin E synthase 3 like (LsPGES3L).

The sequenced LsMGST1L cDNA sequence (Accession no. MW495051) consists of 607 bp, containing an ORF of 432 bp followed by a 135-bp 3′-untranslated region (UTR). The ORF translates into a 146-amino acid (aa) protein identified as a microsomal glutathione S-transferase 1 like (MGST1) protein holding three transmembrane helixes. No other proteins within the membrane-associated proteins in eicosanoid and glutathione metabolism (MAPEG) superfamily were found within the salmon louse genome. Blast searches (BlastP, NCBI) with the LsMGST1L aa sequence showed the highest resemblance to genes annotated as MGST1 and PGES1 in both arthropod and vertebrate species. The PGE2 synthase activity of human PGES1 is glutathione dependent [63], and amino acids important for glutathione binding and enzymatic activity [64, 65] were found to be conserved in LsMGST1L.

The LsPGES3L cDNA sequence (Accession no. MW495052) consists of 952 bp, containing an ORF of 537 bp and a 305-bp 3′-UTR. The ORF translates into a 179-aa protein containing an alpha crystalline-Hsps-p23 like superfamily domain, predicted to have a cytosolic localization. Blast searches (BlastP, NCBI) with LsPGES3L showed the highest resemblance to genes annotated as prostaglandin E synthase 3 from arthropods, especially ants and termites, but also vertebrates and even mammalian PGES3.

Phylogenetic analysis

Using Bayesian (Fig. 1) and ML phylogenies (Additional file 3: Figure S1), we assessed whether homologous protein sequences of genes similar to PGES were clustered correctly with homologous arthropod sequences of similar annotation. Both phylogenies reconstructed large arthropod clades and vertebrate outgroups correctly, even though for some branches MrBayes was unable to determine bifurcation. ML analysis by RAxML produced bifurcations, but with low branch support. Overall, the ML phylogenies were less consistent with species relationships. Both, mGST-like and PGES3-like sequences largely form clades with sequences of the same annotation, and salmon louse sequences are well embedded within clades of closely related taxa. The only outlier with respect to annotation is Caligus rogercresseyi PGES, which is the sister taxon to LsMGST1L and is embedded in a cluster of mGST-like sequences. Notably, all Bayesian phylograms reproduce the major taxonomic groups that contain multiple sequences, such as Crustacea, into which the salmon louse is embedded. As expected, deeper branches indicate higher sequence diversity in invertebrates than in vertebrates.

Fig. 1

Phylogenetic trees for orthologues to EMLSAP00000006733 (MGST1L) (a) and EMLSAP00000012943 (PGES3L) (b). Branch labels indicate branch support in percentage, scale bars correspond to 0.3 (a) and 0.2 (b) substitutions per site. Sequence annotations are given as found. Both phylograms were rooted using mammalian (a) and vertebrate (b) sequences as outgroups. Both salmon lice sequences are well embedded within taxa with the same annotation, except for Caligus rogercresseyi PGE2 synthase (PGES). Major clades (Hymenoptera, Diptera, Nematoda and Crustacea) are indicated by vertical bars

Developmental expression of LsPGESs

LsMGST1L and LsPGES3L transcript levels were determined in all salmon louse developmental stages (Fig. 2a, b). Constitutive expression of both genes was detected throughout lice development, with small variations observed between stages. LsMGST1L did not show a transcriptional response to the shift from free-living stages to parasitic stages or major differences between male and females. Expression of LsPGES3L was similar, except for low levels of transcription in planktonic copepodites and adult males. These results were confirmed by inspection of transcriptional patterns based on RNA sequencing (Additional file: S1; Source:, EMLSAG00000006733 and EMLSAG00000012943).

Fig. 2

Transcriptional profile of LsMGST1L, LsPGES2 and LsPGES3L during development measured by qPCR (ac) and RNA sequencing (df). a, b Transcription level of LsMGST1L and LsPGES3L in all life stages: egg (fertilized), nauplius (nau), copepodites (cop), planktonic copepodites (cop free), cop days post-infestation of fish (cop pdi), chalimus (cha), readult (pad), adult (ad) and mature adult (mat ad). Expression is given as average 2−∆∆Cq ± standard deviation (SD) (N = 5). c Transcription of LsMGST1L, LsPGES2 and LsPGES3L during the reproductive cycle of adult females given as average 2−∆∆Cq ± SD (N = 5). The x-axis indicates the number of days before eggs carried by the female will hatch and new egg extrusion will take place. df Detailed ontogenetic expression profile during the chalimus (Cha1, Chal2) and preadult 1 (Pad1) stages: LsMGST1L (d), LsPGES3L (e) and LsPGES2 (f). Data are presented as average counts per million (CPM) ± SD (N = 3). y Young, m middle, o old, mo molting

The expression patterns of all three LsPGES genes were also investigated during the reproductive cycle of adult females (Fig. 2c) and in a more detailed ontogenetic analysis of animals with different stage-age (Fig. 2d–f). Analysis of LsPGES transcripts during the three molt cycles of chalimus 1, 2 and preadult 1 revealed different patterns of expression. Whereas LsPGES2 was expressed at similar levels throughout the molt cycle of chalimi and the preadult 1 stage, LsMGST1L exhibited a cyclic expression pattern. Higher levels of LsMGST1L were detected in newly molted lice and lice in the middle of the stage, while lower levels were detected in lice directly before the next molt. This cyclic pattern was also evident for LSPGES3L, especially in the expression in chalimus 2 larvae.

Adult female lice produce eggs continuously that are deposited in batches extruded into egg strings shortly after hatching of the previous batch. At 10 °C, this reproductive cycle lasts approximately 10 days [43]. To further determine the involvement of the putative LsPGES genes in reproduction, the transcriptional levels were analyzed during an entire reproductive cycle. The LsPGES transcript levels were, however, not found to change significantly during this cycle (Fig. 2c), although a large variation in the expression of LsPGES3L was evident in females directly after egg string extrusion.

Tissue specific expression of LsPGESs

To investigate the localization of LsMGST1L and LsPGES3L in salmon louse tissues, we performed in situ hybridization in sections of adult females. LsPGES3 transcripts were detected in the ovaries and developing oocytes within the genital segment and in some cells of the subepidermal tissue (Fig. 3). Parallel sections were incubated with the sense probe as a negative control, which displayed only unspecific staining of the cuticula (results not shown).

Fig. 3

Localization of LsPGES3 by in situ analysis of an adult female salmon louse (Lepeophtheirus salmonis). Dark coloring indicates the presence of transcripts. a Overview of the louse with positive staining seen in ovaries, oocytes and subepidermal tissue. Boxes denoted B, C and D indicate the position of the high magnification photos. Scale bar (a): 500 μm. b Unfertilized eggs in genital segment; scale bar: 50 μm. c Subepidermal tissue; scale bar: 50 μm. d Ovary; scale bar: 20 μm

Probably due to lower levels of LsMGST1L in individual cells, LsMGST1L was not successfully localized by in situ hybridization. Therefore, we analyzed the transcript level in the RNA sequencing data from a range of tissues and body segments (Source:, EMLSAG00000012943, EMLSAG00000000441 and EMLSAG00000006733) and in previously published microarray data [45]. LsPGES2 and LsPGES3L exhibited a similar pattern of expression, with the highest transcript level detected in the ovaries, especially in the RNA sequencing data (Fig. 4). LsMGST1L exhibited a slightly different pattern of expression, with the highest transcript level detected in louse intestine. Moreover, none of the genes showed high expression in the sample with a high content of salivary gland or in the leg sample where tegmental type 1–3 glands are present [66].

Fig. 4

mRNA level of the three putative PGES genes in different tissues/body segments analyzed by RNA sequencing (N = 1) (ac) and microarray (N = 4) (df). Error bars indicate standard deviation

PGES KD in larvae

To investigate the functional role of LsPGESs in salmon lice, RNAi was induced in nauplius 1 larvae. KD was performed using dsRNA fragments targeting LsMGST1l, LsPGES2 and LsPGES3L, and the effects on the animals were evaluated as the nauplii developed to the infective copepodid stage. Visual inspection of the animals revealed no abnormal morphology and pattern of movement or molting abnormalities compared to control animals (results not shown).

Due to the suggested role of lice-derived PGE2 in the immunomodulation of susceptible salmonid hosts, LsMGST1L, LsPGES2 and LsPGES3L KD copepodids were used to infest fish with subsequent analysis of the local immune response. At 2 dpi, the levels of all three transcripts were significantly reduced by 95% (LsMGST1L), 82% (LsPGES2) and 75% (LsPGES3L) compared to control animals, but at 14 dpi when the lice were in the chalimus 2 stage, they had returned to the level of the control animals (Fig. 5a). Sampling of fish skin was performed in parallel with sampling of the lice, and the expression of immune-related transcripts was analyzed. The transcript levels of the cytokines interleukin (IL), IL4/13a and IL8 and of the non-specific cytotoxic cell receptor P1 (NCCRP-1) in skin samples was significantly increased, especially at 14 dpi at the site of infestation compared to skin sections taken away from the louse (unaffected) and untreated control fish (Fig. 5b–o). Transcripts of T-cell markers cluster of differentiation 4 (CD4), CD8α, CD8β, major histocompatibility complex 2 (MHC2) and immunoglobulin T (IgT) were all expressed at a lower level when compared to untreated fish. Similar to the trend observed in upregulated transcripts, the response was stronger and—for some genes—only significant at 14 dpi. There were, however, no significant differences between control and KD lice for any of the analyzed genes, and infestation success was also similar between the two groups (Fig. 5p).

Fig. 5

a Knockdown (KD) of MGST1L, PGES2 and PGES3L in copepodids at 2 and chalimus 2 at 14 dpi. Expression level in KD lice were related to that of control lice (2−ΔΔCq ± SD, N = 4). Asterisk denotes significant KD compared to control lice. bo mRNA level of selected immune genes in salmon (S. salar) skin infested with control and KD lice (L. salmonis). Skin was sampled 2 and 14 dpi at the site of lice attachment (infested) and in unaffected sites. The expression level was calculated as 2−∆∆Cq ± SD (N = 5) related to untreated control fish. Lowercase letters above bars indicate statistical significance from the untreated control fish (a) and from unaffected skin (b). p Salmon louse infestation success after KD, given as average number of lice/fish ± SD. IL Interleukin, TNF tumor necrosis factor, NCCRP non-specific cytotoxic cell receptor P, IFN interferon, CD cluster of differentiation, MHC major histocompatibility complex, Ig immunoglobulin, C complement

PGES KD in preadult and adult stages

To investigate the functional role of the putative PGES genes in mobile lice, RNAi was induced in preadult 2 females and followed during their development to adults carrying their second egg string. KD was performed using two mixtures of dsRNA fragments targeting only LsPGES2 and LsPGES3L (RNAi1) or all three putative PGESs (RNAi2). Even though LsMGST1l, LsPGES2 and LsPGES3 were significantly reduced in KD animals (t-test, P < 0.001), all injected lice were successful at re-establishing on the host fish and underwent further development at a normal rate (Table 2). Visual inspection and morphological measurements of lice revealed no significant differences in survival on host fish, egg production, size or morphology compared to control lice (t-test, P > 0.05). Due to the mobile nature of preadult and adult lice, local immune responses were not assessed.

Table 2 Overview of interference RNA experiment injecting preadult 2 female salmon louse (Lepeophtheirus salmonis)

Prostaglandin E2 and blood feeding

Prostaglandin E2 in parasitic lice can either originate from an endogenous source or from the fish blood or skin. To determine whether the content of PGE2 in the louse was affected by blood-meal digestion, the amount of PGE2 metabolites was measured in adult females. Measurements were performed in four groups: females with blood-filled gut and females with a transparent gut with no visible blood; these were subdivided into one group measured just after collection from the fish (0 h) and one group measured 3 h later (Fig. 6). At 0 h, the amount of PGE2 metabolites was equal in both groups (with or without blood). After 3 h, significantly lower levels (fivefold decrease) of PGE2 metabolites were detected in females with transparent guts. A decrease in the amount of PGE2 metabolites from 0 to 3 h was also observed in females with a blood-filled gut, but this decrease was smaller and not significant.

Fig. 6

Indirect measurement of the prostaglandin E2 (PGE2) level (mean ± SD, N = 5) in adult lice with transparent or blood-filled gut directly after sampling (0 h) and at 3 h post-sampling. Asterisk indicates a significant difference in the amount of PGE2 metabolites, determined by t-tests with a threshold P value of 0.05

Oxidative stress and expression of LsMGST1L, LsPGES2 and LsPGES3L

The involvement of prostaglandins and the LsMGST1L, LsPGES2 and LsPGES3L genes during oxidative stress in salmon louse was further analyzed. The transcript levels of the three putative PGES genes in addition to the level of PGE2 and PGF metabolites were measured in LsHPX1 KD copepodids previously shown to have increased levels of several oxidative stress-related genes, such as catalase, glutathione reductase, glutathione S-transferase (GST) mu and zeta and glutathione peroxidase 1 and 4 [57]. Transcription of LsMGST1L, LsPGES2 and LsPGES3L was measured daily in KD animals as the larvae developed from nauplius 2 (3–4 days post hatching [dph]) to copepodids (5–8 dph). The expression of all three transcripts exhibited relatively stable expression in the nauplius stage followed by a peak just after molt to the copepodid stage and a subsequent sharp drop in expression levels at 8 dph (Fig. 7a–c). This drop was not evident in KD animals, and transcription of LsMGST1L, LsPGES2 and LsPGES3L was significantly higher in KD animals compared to controls in 8-day old copepodids. In these larvae, the amount of PGE2 metabolites and PGF was also increased significantly at 8 dph (Fig. 7d, e).

Fig. 7

ac Relative mRNA level of the three putative PGES genes in control and LsHPX1 KD animals. Expression was calculated relative to 18S expression, and the average mRNA level ± SD is shown (N = 3). Arrows indicate the approximate time of molting from nauplius 2 to copepodids. d Analysis of PGE2 metabolites. e PGF metabolites in control and KD animals at 8 days post-hatching (dph). Average metabolite amount ± SD is given (N = 3). Asterisks indicate statistical difference between control and KD animals, determined by t-tests with a threshold P value of 0.05


Salmon lice infestations of farmed Atlantic salmon are a major challenge to the industry and hence to future sustainable production. Scientific approaches to understand the interaction between the fish and this resilient ectoparasite has included searches for parasite substances that manipulate the host. Following the identification of PGE2 in lice secretory/excretory products, a pivotal role of PGE2 as a salmon louse immune-dampening substance was suggested [36] based on its role in immunomodulation by other endo- and ectoparasites [24, 26, 28,29,30, 67]. However, the exact role of PGE2 in the salmon louse host–parasite interaction remains elusive due to technical difficulties obtaining pure secretions from salmon lice glands. Moreover, the presence of additional salmon louse PGES genes might have further complicated functional studies. Hence, the aim of the present study was to identify additional salmon louse PGES genes, and the bioinformatic analysis of the complete salmon louse genome revealed two additional PGES candidates. The results from sequence and phylogenetic analyses support the annotation of LsMGST1L and LsPGES3L, with both sequences well embedded in their respective arthropod clades. However, the enzymatic activity of arthropod PGES has been poorly evaluated; for example, D. melanogaster mGST has been found to have GST but not PGES activity [68]. The present study has also not confirmed the enzymatic activity of the putative LsPGESs. Nevertheless, it has, based on present knowledge, identified those genes within the salmon louse genome that have the highest probability of such activity as they are the genes within the genome that show the highest resemblance to characterized PGES1 and PGES3 genes of both invertebrates and vertebrates.

Recently, functional studies of the salmon louse LsPGES2 gene suggested an endogenous role of PGE2 because an involvement in the host–parasite interaction could not be demonstrated [40]. The current work was initiated to challenge this role, but instead the results presented here further support these findings. If salmon louse-derived PGE2 were to play a major role in host immunomodulation, an upregulation of LsPGES transcripts would be expected at the transition from free-living to parasitic stages or when the louse shifts from the attached chalimus 2 stage to the mobile preadult stage. This latter shift allows the lice to move freely over the fish skin, a phase that can be associated with host mortality [6]. However, our expression analysis of LsMGST1L and LsPGES3 during all developmental stages revealed that both transcripts have a stable expression in all developmental stages, with little difference between free-living and parasitic stages, as previously shown for LsPGES2 [40]. Tissue expression also did not point in the direction of an immune modulative role of salmon louse-derived PGE2, as lower expression was detected in samples enriched with exocrine glands and in situ hybridization localized PGES3L to reproductive organs only. Moreover, KD studies in both copepodids and preadults further supported this conclusion. In mammals, four known PGE2 receptors are important for the downstream effects of PGE2, eliciting both pro-inflammatory and immunosuppressive responses, particularly in the innate immune cells (reviewed in [69]). This complexity of PGE2 signaling is possible not only due to the diversity of receptors that activate different signaling pathways, but also due to different affinities for PGE2 and that they are expressed on a variety of cell types in different combinations. This enables PGE2 to act as an adaptable signaling molecule in a wide range of cell types. While tick-derived PGE2 is believed to have immunosuppressive activity in the host skin [70, 71], PGE2 secreted by the enteric amoeba Entamoeba histolytica is believed to have a pro-inflammatory effect on colonic epithelial cells [72]. In the present study, a clear but dampened immune response was seen at the site of infestation, especially at the later time point when lice had developed to the chalimus 2 stage. This confirms earlier reports that the immune response to lice in fish skin is mainly limited to a small area surrounding the site of infection [35, 73,74,75]. Nevertheless, the local level of key pro-inflammatory transcripts expected to be regulated by salmon louse-derived PGE2 was induced at the same level in control and dsRNA-treated lice, even though a decrease in LsMGST1L, LsPGES2 and LsPGES3L transcripts of between 75–95% was detected. KD copepodids also successfully infested and developed on Atlantic salmon for 14 days at the same rate as the control lice. Likewise, KD in mobile lice stages targeting all three transcripts did not compromise the ability of the louse to parasitize or reproduce on the fish. Taken together, the presented results indicate that LsMGST1L, LsPGES2 and LsPGES3L play insignificant roles in immunomodulation; alternatively, the concentration of residual proteins was sufficient to maintain normal function through the experimental period.

Earlier studies have demonstrated that the amount of the dopamine-induced release of PGE2 in salmon lice secretory/excretory products decreased in the hours after removal from the fish [36], and this has been one of the arguments for an exogenous immunomodulatory role of PGE2. The same decrease of PGE2 was observed in the present study in whole lice homogenates, especially in lice with a transparent empty-looking gut. Interestingly, lice with a blood-filled intestine at the time of removal did not show the same level of PGE2 decrease after 3 h off the host, suggesting that PGE2 may be produced during blood digestion. Both salmon skin and blood cells consumed by the louse can represent a source of PGE2 as these cells are likely to express all enzymes necessary for PGE2 synthesis [11]. Moreover, lysis of blood cells will release toxic components, such as heme and iron, and allow for lipid peroxidation of blood-cell membranes. Thus, digestion of blood leads to an increase in reactive oxidative species, and blood-feeding parasites, including the salmon louse, have evolved mechanisms to handle this stress [76, 77]. Several salmon louse antioxidant genes, such as catalase (EMLSAT00000007315), glutathione reductase (EMLSAT00000000851), GST mu (EMLSAT00000005474) and glutathione peroxidase 1 (EMLSAT00000003186) display relatively high transcript levels within the lice gut ( LsMGST1l also has the highest expression in the lice gut and may act as an antioxidant here.

While non-enzymatic lipid peroxidation typically produces isoprostanes like 8-iso-PGF, enzymatic lipid peroxidation produces PGE2 and especially PGF [78, 79]. We therefore wanted to further explore the involvement of the three putative LsPGES in the production of both PGE2 and PGF metabolites during oxidative stress in the salmon louse. The transcription levels of LsMGST1L, LsPGES2 and LsPGES3L increased in copepodids with increased oxidative stress. Accordingly, the amounts of PGE2 and PGF metabolites were also elevated, indicating a role of prostaglandins in the salmon louse oxidative stress response. This observation is supported by studies of mud crab (Scylla serrata) where both lipid peroxidation and antioxidants werre found to vary during the molt cycle [80, 81]. In salmon louse, transcript levels of PGES genes, particularly LsMGSTL, also varied in different instar ages, as seen for some antioxidant genes, including the above-mentioned catalase, glutathione reductase and glutathione peroxidase 1 [46]. Further studies linking prostaglandins, molt and oxidative stress in salmon louse are thus warranted.

The results of the present study combined with previous findings by [40] also indicate a role of PGE2 in salmon louse reproduction. Tissue-specific analysis revealed the highest level of both LsPGES2 and LsPGES3L transcripts in the ovaries and both transcripts are localized to the ovaries and maturing oocytes within the salmon louse genital segment (Fig. 3 and [40]). PGs are well-known to be involved in ovarian development and egg laying behavior in invertebrates (reviewed by [22, 82, 83]), including crustaceans [18, 20, 21]. PGE2 is also found to be involved in vitellogenesis in the giant freshwater prawn (Macrobrachium rosenbergii), with increased vitellogenin levels in the hemolymph after PGE2 administration [18]. Interestingly, LsPGES3L transcripts were localized to selected cells of the subepidermal tissue, the site of vitellogenin production in salmon louse [43]. However, while KD of a PGES2 ortholog in the lepidopteran insect Spodoptera exigua reduced the egg laying behavior [84], simultaneous KD of LsPGES2 and LsPGES3L did not lead to any measurable decrease in salmon louse reproduction. Furthermore, in animals in which all three putative LsPGES genes were targeted, phenotypic alterations were not detected in either reproduction or development, and none of the transcripts were shown to be regulated during the salmon louse reproductive cycle. In the silkworm Bombyx mori, egg-laying behavior was not affected when virgin individuals were injected with PGE2, showing that PGE2 is not an oviposition-stimulating substance of silkworm [85]. The present study suggests a similar situation in salmon louse, with LsPGES2 and LsPGES3L mRNAs possibly only deposited in the salmon louse oocytes for further downstream application after hatching, although PGE2 biosynthesis has been confirmed in the ovaries of both the black tiger shrimp (Penaeus monodon) and the giant freshwater prawn [18, 21]. It should, however, be mentioned that in B. mori and the parasitic nematode Ascaridia gallia, a sigma-class GST has been shown to have PGES activity [86, 87], and the presence of additional salmon louse PGES genes could also secure reproduction in KD animals if PGE2 has a vital function in salmon louse reproduction. Searches within the salmon louse genome reveal at least 14 genes encoding GSTs (, although none of them annotates as a sigma-class GST. Nevertheless, the search for additional PGES genes in salmon louse should continue to resolve the involvement of PGE2 in both endogenous and exogenous processes.


Putative candidates for immunomodulation of the fish host by salmon louse have attracted interest and speculations, both to understand the molecular interaction but also in the hope of finding candidates for preventive measures. This ectoparasite remains on the fish host for prolonged periods of time with little resulting immunological response. Consequently, it appears prudent to argue that the louse somehow needs to manipulate the host immune response. To date, PGE2 has been the main candidate suggested for this role. However, if the present study has in fact succeeded to address all salmon louse PGESs, a role for PGE2 as immune modulator seems unlikely. This conclusion is based on three major points: (i) expressional analysis did not indicate a role in immunomodulation; for example, upregulation of transcripts in response to the shift to parasitic stages or high expression in samples with an increased content of exocrine glands were not detected; (ii) functional downregulation of all three transcripts by RNAi did not affect the fish immune response nor the ability of the lice to successfully parasitize the fish; (iii) instead, further experiments supported multiple endogenous roles of PGES in reproduction and in handling of oxidative stress associated with molting and hematophagy. However, the multiple endogenous roles of PGE2 together with the presence of multiple PGESs complicates functional studies, and this work could therefore not be regarded as conclusive. Nevertheless, future studies addressing compounds produced by salmon louse exocrine glands are warranted as they are more likely to provide new candidates for immunomodulation of the host.

Availability of data and materials

Data supporting the conclusions of this article are fully available without restriction upon reasonable request.


  1. 1.

    Brandal PO, Egidius E, Romslo I. Host blood: major food component for parasitic copepod Lepeophtheirus salmonis Krøyer, 1838 (Crustacea-Caligidae). Norw J Zool. 1976;24(4):341–3.

    Google Scholar 

  2. 2.

    Kabata Z. Mouth and mode of feeding of caligidae (Copepoda), parasites of fishes, as determined by light and scanning electron microscopy. J Fish Res Board Can. 1974;31(10):1583–8.

    Article  Google Scholar 

  3. 3.

    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.

    Article  CAS  Google Scholar 

  4. 4.

    Johnson SC, Albright LJ. The developmental stages of Lepeophtheirus salmonis (Krøyer, 1837) (Copepoda, Caligidae). Can J Zool. 1991;69(4):929–50.

    Article  Google Scholar 

  5. 5.

    Finstad B, Bjorn PA, Grimnes A, Hvidsten NA. Laboratory and field investigations of salmon lice Lepeophtheirus salmonis (Krøyer) infestation on Atlantic salmon (Salmo salar L.) post-smolts. Aquac Res. 2000;31(11):795–803.

    Article  Google Scholar 

  6. 6.

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

    Article  Google Scholar 

  7. 7.

    Jensen EM, Horsberg TE, Sevatdal S, Helgesen KO. Trends in de-lousing of Norwegian farmed salmon from 2000–2019—consumption of medicines, salmon louse resistance and non-medicinal control methods. PLoS ONE. 2020;15:10e0240894.

    Google Scholar 

  8. 8.

    Bøhn T, Gjelland KO, Serra-Llinares RM, Finstad B, Primicerio R, Nilsen R, et al. Timing is everything: survival of Atlantic salmon Salmo salar postsmolts during events of high salmon lice densities. J Appl Ecol. 2020;57:1149–60.

    Article  Google Scholar 

  9. 9.

    Fjørtoft HB, Nilsen F, Besnier F, Stene A, Bjørn PA, Tveten AK, et al. Salmon lice sampled from wild Atlantic salmon and sea trout throughout Norway display high frequencies of the genotype associated with pyrethroid resistance. Aquacult Env Interac. 2019;11:459–68.

    Article  Google Scholar 

  10. 10.

    Forseth T, Barlaup BT, Finstad B, Fiske P, Gjoaester H, Falkegard M, et al. The major threats to Atlantic salmon in Norway. ICES J Mar Sci. 2017;74(6):1496–513.

    Article  Google Scholar 

  11. 11.

    Kudo I, Murakami M. Prostaglandin E synthase, a terminal enzyme for prostaglandin E-2 biosynthesis. J Biochem Mol Biol. 2005;38(6):633–8.

    CAS  PubMed  Google Scholar 

  12. 12.

    Frolov A, Yang LH, Dong H, Hammock BD, Crofford LJ. Anti-inflammatory properties of prostaglandin E2: Deletion of microsomal prostaglandin E synthase-1 exacerbates non-immune inflammatory arthritis in mice. Prostag Leukotr Ess. 2013;89(5):351–8.

    CAS  Article  Google Scholar 

  13. 13.

    Kawahara K, Hohjoh H, Inazumi T, Tsuchiya S, Sugimoto Y. Prostaglandin E2 induced inflammation: relevance of prostaglandin E receptors. Biochem Biophys Acta. 2015;1851(4):414–21.

    CAS  PubMed  Google Scholar 

  14. 14.

    Nakanishi M, Rosenberg DW. Multifaceted roles of PGE2 in inflammation and cancer. Semin Immunopathol. 2013;35(2):123–37.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Niringiyumukiza JD, Cai HC, Xiang WP. Prostaglandin E2 involvement in mammalian female fertility: ovulation, fertilization, embryo development and early implantation. Reprod Biol Endocrin. 2018;16:43.

    Article  CAS  Google Scholar 

  16. 16.

    Yao CC, Sakata D, Esaki Y, Li YX, Matsuoka T, Kuroiwa K, et al. Prostaglandin E2-EP4 signaling promotes immune inflammation through Th1 cell differentiation and Th17 cell expansion. Nat Med. 2009;15(6):633.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Seo MJ, Oh DK. Prostaglandin synthases: molecular characterization and involvement in prostaglandin biosynthesis. Prog Lipid Res. 2017;66:50–68.

    CAS  PubMed  Article  Google Scholar 

  18. 18.

    Wimuttisuk W, Tobwor P, Deenarn P, Danwisetkanjana K, Pinkaew D, Kirtikara K, et al. Insights into the prostanoid pathway in the ovary development of the penaeid shrimp Penaeus monodon. PLoS ONE. 2013;8(10):e76934.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Hansen K, Varvas K, Jarving I, Samel N. Novel membrane-associated prostaglandin E synthase-2 from crustacean arthropods. Comp Biochem Phys B. 2014;174:45–52.

    CAS  Article  Google Scholar 

  20. 20.

    Duangprom S, Ampansri W, Suwansa-Ard S, Chotwiwatthanakun C, Sobhon P, Kornthong N. Identification and expression of prostaglandin E synthase (PGES) gene in the central nervous system and ovary during ovarian maturation of the female mud crab Scylla olivacea. Anim Reprod Sci. 2018;198:220–32.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Sumpownon C, Engsusophon A, Siangcham T, Sugiyama E, Soonklang N, Meeratana P, et al. Variation of prostaglandin E2 concentrations in ovaries and its effects on ovarian maturation and oocyte proliferation in the giant fresh water prawn Macrobrachium rosenbergii. Gen Comp Endocr. 2015;223:129–38.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Tahara D, Yano I. Maturation-related variations in prostaglandin and fatty acid content of ovary in the kuruma prawn (Marsupenaeus japonicus). Comp Biochem Phys A. 2004;137(4):631–7.

    Article  CAS  Google Scholar 

  23. 23.

    Hampson AJ, Rowley AF, Barrow SE, Steadman R. Biosynthesis of eicosanoids by blood-cells of the crab Carcinus maenas. Biochim Biophys Acta. 1992;1124(2):143–50.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Kutyrev IA, Biserova NM, Olennikov DN, Korneva JV, Mazur OE. Prostaglandins E2 and D2-regulators of host immunity in the model parasite Diphyllobothrium dendriticum: an immunocytochemical and biochemical study. Mol Biochem Parasit. 2017;212:33–45.

    CAS  Article  Google Scholar 

  25. 25.

    Daugschies A, Joachim A. Eicosanoids in parasites and parasitic infections. Adv Parasit. 2000;46:181–240.

    CAS  Article  Google Scholar 

  26. 26.

    de Almeida PE, Toledo DAM, Rodrigues GSC, D’Avila H. Lipid bodies as sites of prostaglandin E2 synthesis during chagas disease: Impact in the parasite escape mechanism. Front Microbiol. 2018;9:499.

    PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Kubata BK, Duszenko M, Martin S, Urade Y. Molecular basis for prostaglandin production in hosts and parasites. Trends Parasitol. 2007;23(7):325–31.

    CAS  PubMed  Article  Google Scholar 

  28. 28.

    Bakshi M, Kim TK, Porter L, Mwangi W, Mulenga A. Amblyomma americanum ticks utilizes countervailing pro and anti-inflammatory proteins to evade host defense. Plos Pathog. 2019;15(11):e1008128.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  29. 29.

    Bowman AS, Dillwith JW, Sauer JR. Tick salivary prostaglandins: presence, origin and significance. Parasitol Today. 1996;12(10):388–96.

    CAS  PubMed  Article  Google Scholar 

  30. 30.

    Esteves E, Bizzarro B, Costa FB, Ramirez-Hernandez A, Peti APF, Cataneo AHD, et al. Amblyomma sculptum salivary PGE2 modulates the dendritic cell-Rickettsia rickettsii interactions in vitro and in vivo. Front Immunol. 2019;10:118.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Kutyrev IA, Franke F, Buscher J, Kurtz J, Scharsack JP. In vitro effects of prostaglandin E2 on leucocytes from sticklebacks (Gasterosteus aculeatus) infected and not infected with the cestode Schistocephalus solidus. Fish Shellfish Immun. 2014;41(2):473–81.

    CAS  Article  Google Scholar 

  32. 32.

    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.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  33. 33.

    Jones MW, Sommerville C, Bron J. The histopathology associated with the juvenile stages of Lepeophtheirus salmonis on the Atlantic Salmon Salmo salar. L J Fish Dis. 1990;13(4):303–10.

    Article  Google Scholar 

  34. 34.

    Jones SRM. Lepeophtheirus salmonis. J Aquac Res Dev. 2011.

    Article  Google Scholar 

  35. 35.

    Øvergård AC, Hamre LA, Grotmol S, Nilsen F. Salmon louse rhabdoviruses: impact on louse development and transcription of selected Atlantic salmon immune genes. Dev Comp Immunol. 2018;86:86–95.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  36. 36.

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

    CAS  PubMed  Article  Google Scholar 

  37. 37.

    Fast MD, Ross NW, Johnson SC. Prostaglandin E2 modulation of gene expression in an Atlantic salmon (Salmo salar) macrophage-like cell line (SHK-1). Dev Comp Immunol. 2005;29(11):951–63.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Gamil AAA, Guo TC, Konig M, Evensen O. Distribution of EP4 receptor in different Atlantic salmon (Salmo salar L.) tissues. Dev Comp Immunol. 2015;48(1):143–50.

    CAS  PubMed  Article  Google Scholar 

  39. 39.

    Hamilton S, McLean K, Monaghan SJ, McNair C, Inglis NF, McDonald H, et al. Characterisation of proteins in excretory/secretory products collected from salmon lice Lepeophtheirus salmonis. Parasites Vectors. 2018;11:294.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  40. 40.

    Eichner C, Overgard AC, Nilsen F, Dalvin S. Molecular characterization and knock-down of salmon louse (Lepeophtheirus salmonis) prostaglandin E synthase. Exp Parasitol. 2015;159:79–93.

    CAS  PubMed  Article  Google Scholar 

  41. 41.

    Skern-Mauritzen R, Torrissen O, Glover KA. Pacific and Atlantic Lepeophtheirus salmonis (Kroyer, 1838) are allopatric subspecies: Lepeophtheirus salmonis salmonis and L. salmonis oncorhynchi subspecies novo. BMc Genet. 2014;15:32.

    PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Hamre LA, Glover KA, Nilsen F. Establishment and characterisation of salmon louse (Lepeophtheirus salmonis (Krøyer 1837)) laboratory strains. Parasitol Int. 2009;58(4):451–60.

    PubMed  Article  Google Scholar 

  43. 43.

    Dalvin S, Frost P, Loeffen P, Skern-Mauritzen R, Baban J, Ronnestad I, et al. Characterisation of two vitellogenins in the salmon louse Lepeophtheirus salmonis: molecular, functional and evolutional analysis. Dis Aquat Organ. 2011;94(3):211–24.

    CAS  PubMed  Article  Google Scholar 

  44. 44.

    Øvergård AC, Nerland AH, Patel S. Evaluation of potential reference genes for real time RT-PCR studies in Atlantic halibut (Hippoglossus hippoglossus L.); during development, in tissues of healthy and NNV-injected fish, and in anterior kidney leucocytes. BMC Mol Biol. 2010;11:36.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  45. 45.

    Edvardsen RB, Dalvin S, Furmanek T, Malde K, Maehle S, Kvamme BO, et al. Gene expression in five salmon louse (Lepeophtheirus salmonis Krøyer 1837) tissues. Mar Genom. 2014;18:39–44.

    Article  Google Scholar 

  46. 46.

    Eichner C, Dondrup M, Nilsen F. RNA sequencing reveals distinct gene expression patterns during the development of parasitic larval stages of the salmon louse (Lepeophtheirus salmonis). J Fish Dis. 2018;41(6):1005–29.

    CAS  PubMed  Article  Google Scholar 

  47. 47.

    Agarwala R, Barrett T, Beck J, Benson DA, Bollin C, Bolton E, et al. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2016;44(D1):D7–19.

    CAS  Article  Google Scholar 

  48. 48.

    Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.

    CAS  Article  Google Scholar 

  49. 49.

    Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  50. 50.

    Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ. Jalview Version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25(9):1189–91.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, et al. MrBayes 32: efficient bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61(3):539–42.

    PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  53. 53.

    Darriba D, Taboada GL, Doallo R, Posada D. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011;27(8):1164–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  54. 54.

    Dalvin S, Nilsen F, Skern-Mauritzen R. Localization and transcription patterns of LsVasa, a molecular marker of germ cells in Lepeophtheirus salmonis (Krøyer). J Nat Hist. 2013;47(5–12):889–900.

    Article  Google Scholar 

  55. 55.

    Trösse C, Nilsen F, Dalvin S. RNA interference mediated knockdown of the KDEL receptor and COPB2 inhibits digestion and reproduction in the parasitic copepod Lepeophtheirus salmonis. Comp Biochem Phys B. 2014;170:1–9.

    Article  CAS  Google Scholar 

  56. 56.

    Frost P, Nilsen F. Validation of reference genes for transcription profiling in the salmon louse, Lepeophtheirus salmonis, by quantitative real-time PCR. Vet Parasitol. 2003;118(1–2):169–74.

    CAS  PubMed  Article  Google Scholar 

  57. 57.

    Øvergård AC, Eichner C, Nilsen F, Dalvin S. Molecular characterization and functional analysis of a salmon louse (Lepeophtheirus salmonis, Krøyer 1838) heme peroxidase with a potential role in extracellular matrixes. Comp Biochem Phys A. 2017;206:1–10.

    Article  CAS  Google Scholar 

  58. 58.

    Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A. Guideline to reference gene selection for quantitative real-time PCR. Biochem Bioph Res Co. 2004;313(4):856–62.

    CAS  Article  Google Scholar 

  59. 59.

    Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9):e45.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  60. 60.

    Eichner C, Nilsen F, Grotmol S, Dalvin S. A method for stable gene knock-down by RNA interference in larvae of the salmon louse (Lepeophtheirus salmonis). Exp Parasitol. 2014;140:44–51.

    CAS  PubMed  Article  Google Scholar 

  61. 61.

    Dalvin S, Frost P, Biering E, Hamre LA, Eichner C, Krossoy B, et al. Functional characterisation of the maternal yolk-associated protein (LsYAP) utilising systemic RNA interference in the salmon louse (Lepeophtheirus salmonis) (Crustacea: Copepoda). Int J Parasitol. 2009;39(13):1407–15.

    CAS  PubMed  Article  Google Scholar 

  62. 62.

    Holm HJ, Skugor S, Bjelland AK, Radunovic S, Wadsworth S, Koppang EO, et al. Contrasting expression of immune genes in scaled and scaleless skin of Atlantic salmon infected with young stages of Lepeophtheirus salmonis. Dev Comp Immunol. 2017;67:153–65.

    CAS  PubMed  Article  Google Scholar 

  63. 63.

    Jakobsson PJ, Thoren S, Morgenstern R, Samuelsson B. Identification of human prostaglandin E synthase: a microsomal, glutathione-dependent, inducible enzyme, constituting a potential novel drug target. Proc Natl Acad Sci USA. 1999;96(13):7220–5.

    CAS  PubMed  Article  Google Scholar 

  64. 64.

    Jegerschöld C, Pawelzik SC, Purhonen P, Bhakat P, Gheorghe KR, Gyobu N, et al. Structural basis for induced formation of the inflammatory mediator prostaglandin E2. Proc Natl Acad Sci USA. 2008;105(32):11110–5.

    PubMed  Article  Google Scholar 

  65. 65.

    Sjögren T, Nord J, Ek M, Johansson P, Liu G, Geschwindner S. Crystal structure of microsomal prostaglandin E2 synthase provides insight into diversity in the MAPEG superfamily. Proc Natl Acad Sci USA. 2013;110(10):3806–11.

    PubMed  Article  Google Scholar 

  66. 66.

    Øvergård AC, Hamre LA, Harasimczuk E, Dalvin S, Nilsen F, Grotmol S. Exocrine glands of Lepeophtheirus salmonis (Copepoda: Caligidae): distribution, developmental appearance, and site of secretion. J Morphol. 2016;277(12):1616–30.

    PubMed  Article  CAS  Google Scholar 

  67. 67.

    Daugschies A, Joachim A. Eicosanoids in parasites and parasitic infections. Adv Parasitol. 2000;46:181–240.

    CAS  PubMed  Article  Google Scholar 

  68. 68.

    Bresell A, Weinander R, Lundqvist G, Raza H, Shimoji M, Sun TH, et al. Bioinformatic and enzymatic characterization of the MAPEG superfamily. FEBS J. 2005;272(7):1688–703.

    CAS  PubMed  Article  Google Scholar 

  69. 69.

    Agard M, Asakrah S, Morici LA. PGE2 suppression of innate immunity during mucosal bacterial infection. Front Cell Infect Microbiol. 2013;3:45.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  70. 70.

    Oliveira CJF, Sa-Nunes A, Francischetti IMB, Carregaro V, Anatriello E, Silva JS, et al. Deconstructing tick saliva non-protein molecules with potent iImmunomodulatory properties. J Biol Chem. 2011;286(13):10960–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  71. 71.

    Poole NM, Mamidanna G, Smith RA, Coons LB, Cole JA. Prostaglandin E2 in tick saliva regulates macrophage cell migration and cytokine profile. Parasites Vectors. 2013;6(1):261.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  72. 72.

    Dey I, Chadee K. Prostaglandin E2 produced by Entamoeba histolytica binds to EP4 receptors and stimulates interleukin-8 production in human colonic cells. Infect Immun. 2008;76(11):5158–63.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  73. 73.

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

    CAS  Google Scholar 

  74. 74.

    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(1):178–91.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  75. 75.

    Dalvin S, Jørgensen LVG, Kania PW, Buchmann K, Overgard AC. Rainbow trout Oncorhynchus mykiss skin responses to salmon louse Lepeophtheirus salmonis: From copepodid to adult stage. Fish Shellfish Immun. 2020;103:200–10.

    CAS  Article  Google Scholar 

  76. 76.

    Heggland EI, Dondrup M, Nilsen F, Eichner C. Host gill attachment causes blood-feeding by the salmon louse (Lepeophtheirus salmonis) chalimus larvae and alters parasite development and transcriptome. Parasite Vector. 2020;13:1.

    Article  CAS  Google Scholar 

  77. 77.

    Hernandez EP, Kusakisako K, Talactac MR, Galay RL, Hatta T, Matsuo T, et al. Characterization and expression analysis of a newly identified glutathione S-transferase of the hard tick Haemaphysalis longicornis during blood-feeding. Parasites Vectors. 2018;11:91.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  78. 78.

    Basu S. Carbon tetrachloride-induced lipid peroxidation: eicosanoid formation and their regulation by antioxidant nutrients. Toxicology. 2003;189(1–2):113–27.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  79. 79.

    Korbecki J, Baranowska-Bosiacka I, Gutowska I, Chlubek D. The effect of reactive oxygen species on the synthesis of prostanoids from arachidonic acid. J Physiol Pharmacol. 2013;64(4):409–21.

    CAS  PubMed  Google Scholar 

  80. 80.

    Salaenoi J, Setthamongkol P, Srimeetian P, Satapornvit K. Lipid peroxidation and antioxidant activity in molting stages of mud crab (Scylla serrata). Kasetsart J Nat Sci. 2015;49:422–32.

    CAS  Google Scholar 

  81. 81.

    Salaenoi J, Srimeetian P, Mingmuang M. Variations of catalase and glutathione activities in molting cycle of mud crab (Scylla serrata). Kasetsart J Nat Sci. 2014;48:64–71.

    CAS  Google Scholar 

  82. 82.

    Reddy PS, Reddy PR, Nagaraju GPC. The synthesis and effects of prostaglandins on the ovary of the crab Oziotelphusa senex senex. Gen Comp Endocr. 2004;135(1):35–41.

    CAS  PubMed  Article  Google Scholar 

  83. 83.

    Stanley-Samuelson DW, Pedibhotla VK. What can we learn from prostaglandins and related eicosanoids in insects? Insect Biochem Mol. 1996;26(3):223–34.

    CAS  Article  Google Scholar 

  84. 84.

    Ahmed S, Stanley D, Kim Y. An insect prostaglandin E2 synthase acts in immunity and reproduction. Front Physiol. 2018;9:1231.

    PubMed  PubMed Central  Article  Google Scholar 

  85. 85.

    Yamauchi M, Fugo H, Dedos SG. Prostaglandins do not release egg-laying behaviour in the silkmoth Bombyx mori. Zool Sci. 1997;14(1):135–40.

    CAS  Article  Google Scholar 

  86. 86.

    Meyer DJ, Muimo R, Thomas M, Coates D, Isaac RE. Purification and characterization of prostaglandin-H E-isomerase, a sigma-class glutathione S-transferase, from Ascaridia galli. Biochem J. 1996;313:223–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  87. 87.

    Yamamoto K, Higashiura A, Suzuki M, Aritake K, Urade Y, Uodome N, et al. Crystal structure of a Bombyx mori sigma-class glutathione transferase exhibiting prostaglandin E synthase activity. Biochem Biophys Acta. 2013;1830(6):3711–8.

    CAS  PubMed  Article  Google Scholar 

Download references


We thank Heidi Kongshaug, Lars Are Hamre and Per Gunnar Espedal for help in the laboratory and during fish trials, and Sindre Grotmol for discussions.


This research was funded by the Research Council Norway, SFI-Sea Lice Research Centre, Grant Number 203513/O30.

Author information




SD took part in designing the work, acquisition and interpretation of data and drafted the manuscript. CE took part in designing the work, acquisition, analysis and interpretation of data and substantially revised the manuscript. MD took part in designing the work, analysis and interpretation of data and substantially revised the manuscript. ACØ took part in designing the work, acquisition, analysis and interpretation of data and took part in drafting the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Aina-Cathrine Øvergård.

Ethics declarations

Ethics approval and consent to participate

All experimental procedures were performed in accordance with national legislation for animal welfare, approved by the Animal Ethics Committee by the Norwegian Food Safety Authority (approval number 4538) at University of Bergen.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

 Sequence alignments and accession numbers used in the phylogenetic analysis of EMLSAT00000006733 in FASTA format, gzipped.

Additional file 2.

 Sequence alignments and accession numbers used in the phylogenetic analysis of EMLSAT00000012943 in FASTA format, gzipped.

Additional file 3: Figure S1.

 Maximum Likelihood phylogenies of orthologues to EMLSAT00000006733 (LsMGST1L) (a) and EMLSAT00000012943 (LsPGES3L) (b). Branch labels indicate branch support in percentage, scale bars correspond to 0.3 substitutions per site. Sequence annotations are given as found, LsMGST1L = L_salmonis_hPGES in a. Both phylograms were rooted using mammalian (a) and vertebrate (b) sequences as outgroups.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Dalvin, S., Eichner, C., Dondrup, M. et al. Roles of three putative salmon louse (Lepeophtheirus salmonis) prostaglandin E2 synthases in physiology and host–parasite interactions. Parasites Vectors 14, 206 (2021).

Download citation


  • Invertebrate
  • Arthropod
  • Copepod
  • Blood-feeding
  • RNA interference
  • ROS
  • Immune response