RNA sequencing reveals distinct gene expression patterns during the development of parasitic larval stages of the salmon louse (Lepeophtheirus salmonis)

Salmon louse labial gland enzymes: implications for host settlement and immune modulation

Salmon louse (Lepeophtheirus salmonis) is a skin- and blood-feeding ectoparasite, infesting salmonids. While feeding, labial gland proteins from the salmon louse may be deposited on the Atlantic salmon (Salmo salar) skin. Previously characterized labial gland proteins are involved in anti-coagulation and may contribute to inhibiting Atlantic salmon from mounting a sufficient immune response against the ectoparasite. As labial gland proteins seem to be important in the host-parasite interaction, we have, therefore, identified and characterized ten enzymes localized to the labial gland. They are a large group of astacins named L. salmonis labial gland astacin 1-8 (LsLGA 1-8), one serine protease named L. salmonis labial gland serine protease 1 (LsLGSP1), and one apyrase named L. salmonis labial gland apyrase 1 (LsLGAp1). Protein domain predictions showed that LsLGA proteins all have N-terminal ShK domains, which may bind to potassium channels targeting the astacins to its substrate. LsLGA1 and -4 are, in addition, expressed in another gland type, whose secrete also meets the host-parasite interface. This suggests that LsLGA proteins may have an anti-microbial function and may prevent secondary infections in the wounds. LsLGAp1 is predicted to hydrolyze ATP or AMP and is, thereby, suggested to have an immune dampening function. In a knockdown study targeting LsLGSP1, a significant increase in IL-8 and MMP13 at the skin infestation site was seen under LsLGSP1 knockdown salmon louse compared to the control, suggesting that LsLGSP1 may have an anti-inflammatory effect. Moreover, most of the identified labial gland proteins are expressed in mature copepodids prior to host settlement, are not regulated by starvation, and are expressed at similar or higher levels in lice infesting the salmon louse-resistant pink salmon (Oncorhynchus gorbuscha). This study, thereby, emphasizes the importance of labial gland proteins for host settlement and their immune dampening function. This work can further contribute to anti-salmon louse treatment such as vaccine development, functional feed, or gene-edited salmon louse-resistant Atlantic salmon.

Identification and characterization of two salmon louse heme peroxidases and their potential as vaccine antigens

Salmon louse, Lepeophtheirus salmonis, represents major challenge for salmon farming. Current treatments impose welfare issues and are costly, whereas prophylactic measures are unavailable. Two salmon louse heme peroxidases (LsPxtl-1 and LsPxtl-2) were tested for their importance for parasite development and as potential vaccine candidates. LsPxtl-1 possesses two heme peroxidase domains and is expressed in ovaries and gut, whereas LsPxtl-2 encodes one domain and contains N-terminal signal peptide and an Eph receptor ligand-binding domain. LsPxtl-1, but not LsPxtl-2, knockdown in nauplius II stage caused poor swimming and death, indicating its importance for parasite development. Immunizations using single DNA plasmid injection encoding the peroxidases or heterologous prime (DNA) and boost (recombinant LsPxtl-2 protein) gave non-significant reduction in lice numbers. Single injection gave low specific antibody levels compared with the prime-boost. The findings suggest LsPxtl-1 is important for parasite development but formulations and vaccination modalities used did not significantly reduce lice infestation.

Repeated exposure affects susceptibility and responses of Atlantic salmon (Salmo salar) towards the ectoparasitic salmon lice (Lepeophtheirus salmonis)

Atlantic salmon (Salmo salar) is repeatedly exposed to and infected with ectoparasitic salmon lice (Lepeophtheirus salmonis) both in farms and in nature. However, this is not reflected in laboratory experiments where fish typically are infected only once. To investigate if a previous lice infection affects host response to subsequent infections, fish received 4 different experimental treatments; including 2 groups of fish that had previously been infected either with adult or infective salmon lice larvae (copepodids). Thereafter, fish in all treatment groups were infected with either a double or a single dose of copepodids originating from the same cohort. Fish were sampled when lice had developed into the chalimus, the pre-adult and the adult stage, respectively. Both the specific growth rate and cortisol levels (i.e. a proxy for stress) of the fish differed between treatments. Lice success (i.e. ability to infect and survive on the host) was higher in naïve than in previously infected fish (pre-adult stage). The expression of immune and wound healing transcripts in the skin also differed between treatments, and most noticeable was a higher upregulation early in the infection in the group previously infected with copepodids. However, later in the infection, the least upregulation was observed in this group, suggesting that previous exposure to salmon lice affects the response of Atlantic salmon towards subsequent lice infections.

Mining Lepeophtheirus salmonis RNA-Seq data for qPCR reference genes and their application in Caligus elongatus

Lepeophtheirus salmonis and Caligus elongatus are two parasitic copepod species posing a significant threat to salmonid aquaculture. Consequently, several gene expression studies are executed each year to gain new knowledge and treatment strategies. Though, to enable accurate gene expression measurements by quantitative real time PCR, stable reference genes are needed. Previous studies have mainly focused on a few genes selected based on their function as housekeeping genes, as these are often stably expressed in various cells and tissues. In the present study, however, RNA-sequencing data from 127 L. salmonis samples from different life stages and diverse environmental conditions were used to identify new candidate reference genes displaying low variation. From this, six genes were selected, and the stability validated by qPCR on samples from different life stages. Since neither a genome nor comprehensive RNA sequencing data are available for C. elongatus, homologous genes to those identified for L. salmonis were identified within a C. elongatus transcriptome assembly and validated by qPCR in different life stages. Overall, the genes eukaryotic translation initiation factor 1A (EIF1A) and serine/threonine-protein phosphatase 1 (PP1) displayed the highest stability in L. salmonis, while the combination of PP1 and ribosomal protein S13 (RPS13) was found to have the highest stability in C. elongatus. These genes are well-suited reference genes for qPCR applications which allow for accurate normalization of target genes.

The effect of different intensities of the ectoparasitic salmon lice (Lepeophtheirus salmonis) on Atlantic salmon (Salmo salar)

The effect of different intensities of the ectoparasitic salmon lice (Lepeophtheirus salmonis) on stress, growth and the expression of immune and wound healing transcripts in the skin of Atlantic salmon (Salmo salar) was investigated. Lice infection success and survival were similar at the chalimus and preadult stage in the low and high dose group, but infection success and survival were significantly lower in the high than in the low dose group at the adult stage. The expression of investigated transcripts was not correlated to lice intensities, but several of them were significantly differently expressed locally in the skin at the site of lice attachment in infected fish compared to controls. This included an up-regulation of pro-inflammatory markers at the site of lice attachment (e.g., interleukin 1-beta, interleukin 8 and the acute phase protein serum amyloid A), a reduction of markers of adaptive immunity (cluster of differentiation 8-alpha and immunoglobulin M) and decreased expression of the anti-inflammatory cytokine interleukin 10.

Small, charged proteins in salmon louse (Lepeophtheirus salmonis) secretions modulate Atlantic salmon (Salmo salar) immune responses and coagulation

Little is known about glandular proteins secreted from the skin- and blood-feeding ectoparasite salmon louse (Lepeophtheirus salmonis). The labial gland has ducts extending into the oral cavity of the lice, and the present study aimed to identify novel genes expressed by this gland type and to investigate their role in modulation of host parameters at the lice feeding site. Five genes associated with labial gland function were identified and named Lepeophteirus salmonis labial gland protein (LsLGP) 1-4 and 1 like (LsLGP1L). All LsLGPs were predicted to be small charged secreted proteins not encoding any known protein domains. Functional studies revealed that LsLGP1 and/or LsLGP1L regulated the expression of other labial gland genes. Immune dampening functions were indicated for LsLGP2 and 3. Whereas LsLGP2 was expressed throughout the parasitic life cycle and found to dampen inflammatory cytokines, LsLGP3 displayed an increased expression in mobile stages and appeared to dampen adaptive immune responses. Expression of LsLGP4 coincided with moulting to the mobile pre-adult I stage where hematophagous feeding is initiated, and synthetic LsLGP4 decreased the clotting time of Atlantic salmon plasma. Results from the present study confirm that the salmon louse secretes immune modulating and anti-coagulative proteins with a potential application in new immune based anti-salmon louse treatments.

The salmon louse genome may be much larger than sequencing suggests

The genome size of organisms impacts their evolution and biology and is often assumed to be characteristic of a species. Here we present the first published estimates of genome size of the ecologically and economically important ectoparasite, Lepeophtheirus salmonis (Copepoda, Caligidae). Four independent L. salmonis genome assemblies of the North Atlantic subspecies Lepeophtheirus salmonis salmonis, including two chromosome level assemblies, yield assemblies ranging from 665 to 790 Mbps. These genome assemblies are congruent in their findings, and appear very complete with Benchmarking Universal Single-Copy Orthologs analyses finding > 92% of expected genes and transcriptome datasets routinely mapping > 90% of reads. However, two cytometric techniques, flow cytometry and Feulgen image analysis densitometry, yield measurements of 1.3-1.6 Gb in the haploid genome. Interestingly, earlier cytometric measurements reported genome sizes of 939 and 567 Mbps in L. salmonis salmonis samples from Bay of Fundy and Norway, respectively. Available data thus suggest that the genome sizes of salmon lice are variable. Current understanding of eukaryotic genome dynamics suggests that the most likely explanation for such variability involves repetitive DNA, which for L. salmonis makes up ≈ 60% of the genome assemblies.

Sex differences in the early life stages of the salmon louse Lepeophtheirus salmonis (Copepoda: Caligidae)

Salmon lice are ectoparasites on salmonids and feed on blood, mucus, and skin from their hosts. This causes high annual costs for treatment and control for the aquaculture industry. Salmon lice have a life cycle consisting of eight life stages. Sex determination by eye is only possible from the sixth stage onwards. A molecular sex determination has not been carried out so far, even though few individual sex-linked SNPs have been reported. In the present study, we used known sex-specific SNPs as a basis to sequence the complete sex-specific gene variants and used the sequence information to develop a sex determination assay. This assay could be used to determine the developmental speed of the two sexes already in the earliest life stages. Additionally, we sampled salmon lice in the nauplius II stage, determined the sex of each individual, pooled their RNA according to their sex, and used RNA sequencing to search for differences in gene expression and further sex-specific SNPs. We succeeded in developing a sex-determination assay that works on DNA or RNA from even the earliest larval stages of the salmon louse after hatching. At these early developmental stages, male salmon lice develop slightly quicker than females. We detected several previously unknown, sex-specific SNPs in our RNA-data seq, but only very few genes showed a differential expression between the sexes. Potential connections between SNPs, gene expression, and development are discussed.

A novel approach to co-expression network analysis identifies modules and genes relevant for moulting and development in the Atlantic salmon louse (Lepeophtheirus salmonis)

BACKGROUND

The salmon louse (Lepeophtheirus salmonis) is an obligate ectoparasitic copepod living on Atlantic salmon and other salmonids in the marine environment. Salmon lice cause a number of environmental problems and lead to large economical losses in aquaculture every year. In order to develop novel parasite control strategies, a better understanding of the mechanisms of moulting and development of the salmon louse at the transcriptional level is required.

METHODS

Three weighted gene co-expression networks were constructed based on the pairwise correlations of salmon louse gene expression profiles at different life stages. Network-based approaches and gene annotation information were applied to identify genes that might be important for the moulting and development of the salmon louse. RNA interference was performed for validation. Regulatory impact factors were calculated for all the transcription factor genes by examining the changes in co-expression patterns between transcription factor genes and deferentially expressed genes in middle stages and moulting stages.

RESULTS

Eight gene modules were predicted as important, and 10 genes from six of the eight modules have been found to show observable phenotypes in RNA interference experiments. We knocked down five hub genes from three modules and observed phenotypic consequences in all experiments. In the infection trial, no copepodids with a RAB1A-like gene knocked down were found on fish, while control samples developed to chalimus-1 larvae. Also, a FOXO-like transcription factor obtained highest scores in the regulatory impact factor calculation.

CONCLUSIONS

We propose a gene co-expression network-based approach to identify genes playing an important role in the moulting and development of salmon louse. The RNA interference experiments confirm the effectiveness of our approach and demonstrated the indispensable role of a RAB1A-like gene in the development of the salmon louse. We propose that our approach could be generalized to identify important genes associated with a phenotype of interest in other organisms.

The salmon louse genome: Copepod features and parasitic adaptations

Copepods encompass numerous ecological roles including parasites, detrivores and phytoplankton grazers. Nonetheless, copepod genome assemblies remain scarce. Lepeophtheirus salmonis is an economically and ecologically important ectoparasitic copepod found on salmonid fish. We present the 695.4 Mbp L. salmonis genome assembly containing ≈60% repetitive regions and 13,081 annotated protein-coding genes. The genome comprises 14 autosomes and a ZZ-ZW sex chromosome system. Assembly assessment identified 92.4% of the expected arthropod genes. Transcriptomics supported annotation and indicated a marked shift in gene expression after host attachment, including apparent downregulation of genes related to circadian rhythm coinciding with abandoning diurnal migration. The genome shows evolutionary signatures including loss of genes needed for peroxisome biogenesis, presence of numerous FNII domains, and an incomplete heme homeostasis pathway suggesting heme proteins to be obtained from the host. Despite repeated development of resistance against chemical treatments L. salmonis exhibits low numbers of many genes involved in detoxification.

The FTZ-F1 gene encodes two functionally distinct nuclear receptor isoforms in the ectoparasitic copepod salmon louse (Lepeophtheirus salmonis)

The salmon louse, Lepeophtheirus salmonis, is an ectoparasitic crustacean that annually inflicts substantial losses to the aquaculture industry in the northern hemisphere and poses a threat to the wild populations of salmonids. The salmon louse life cycle consists of eight developmental stages each separated by a molt. Fushi Tarazu Factor-1 (FTZ-F1) is an ecdysteroid-regulated gene that encodes a member of the NR5A family of nuclear receptors that is shown to play a crucial regulatory role in molting in insects and nematodes. Characterization of an FTZ-F1 orthologue in the salmon louse gave two isoforms named αFTZ-F1 and βFTZ-F1, which are identical except for the presence of a unique N-terminal domain (A/B domain). A comparison suggest conservation of the FTZ-F1 gene structure among ecdysozoans, with the exception of nematodes, to produce isoforms with unique N-terminal domains through alternative transcription start and splicing. The two isoforms of the salmon louse FTZ-F1 were expressed in different amounts in the same tissues and showed a distinct cyclical expression pattern through the molting cycle with βFTZ-F1 being the highest expressed isoform. While RNA interference knockdown of βFTZ-F1 in nauplius larvae and in pre-adult males lead to molting arrest, knockdown of βFTZ-F1 in pre-adult II female lice caused disruption of oocyte maturation at the vitellogenic stage. No apparent phenotype could be observed in αFTZ-F1 knockdown larvae, or in their development to adults, and no genes were found to be differentially expressed in the nauplii larvae following αFTZ-F1 knockdown. βFTZ-F1 knockdown in nauplii larvae caused both down and upregulation of genes associated with proteolysis and chitin binding and affected a large number of genes which are in normal salmon louse development expressed in a cyclical pattern. This is the first description of FTZ-F1 gene function in copepod crustaceans and provides a foundation to expand the understanding of the molecular mechanisms of molting in the salmon louse and other copepods.

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

BACKGROUND

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 E₂ (PGE₂) is important in the salmon louse host-parasite interaction, although studies of a salmon louse prostaglandin E₂ 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.

METHODS

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.

RESULTS

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.

CONCLUSIONS

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.

Host gill attachment causes blood-feeding by the salmon louse (Lepeophtheirus salmonis) chalimus larvae and alters parasite development and transcriptome

BACKGROUND

Blood-feeding is a common strategy among parasitizing arthropods, including the ectoparasitic salmon louse (Lepeophtheirus salmonis), feeding off its salmon host's skin and blood. Blood is rich in nutrients, among these iron and heme. These are essential molecules for the louse, yet their oxidative properties render them toxic to cells if not handled appropriately. Blood-feeding might therefore alter parasite gene expression.

METHODS

We infected Atlantic salmon with salmon louse copepodids and sampled the lice in two different experiments at day 10 and 18 post-infestation. Parasite development and presence of host blood in their intestines were determined. Lice of similar instar age sampled from body parts with differential access to blood, namely from gills versus lice from skin epidermis, were analysed for gene expression by RNA-sequencing in samples taken at day 10 for both experiments and at day 18 for one of the experiments.

RESULTS

We found that lice started feeding on blood when becoming mobile preadults if sitting on the fish body; however, they may initiate blood-feeding at the chalimus I stage if attached to gills. Lice attached to gills develop at a slower rate. By differential expression analysis, we found 355 transcripts elevated in lice sampled from gills and 202 transcripts elevated in lice sampled from skin consistent in all samplings. Genes annotated with "peptidase activity" were among the ones elevated in lice sampled from gills, while in the other group genes annotated with "phosphorylation" and "phosphatase" were pervasive. Transcripts elevated in lice sampled from gills were often genes relatively highly expressed in the louse intestine compared with other tissues, while this was not the case for transcripts elevated in lice sampled from skin. In both groups, more than half of the transcripts were from genes more highly expressed after attachment.

CONCLUSIONS

Gill settlement results in an alteration in gene expression and a premature onset of blood-feeding likely causes the parasite to develop at a slower pace.

Heavy and light chain homologs of ferritin are essential for blood-feeding and egg production of the ectoparasitic copepod Lepeophtheirus salmonis

The salmon louse, Lepeophtheirus salmonis, is a hematophagous ectoparasite of salmonid fish. Due to its blood-feeding activity, the louse is exposed to great amounts of iron, which is an essential, yet potentially toxic mineral. The major known iron storage protein is ferritin, which the salmon louse encodes four genes of (LsFer1-4). Two of the ferritins are predicted to be secreted. These are one of the heavy chain homologs (LsFer1) and the light chain homolog (LsFer2). Here, we perform functional studies and characterize the two secreted ferritins. Our results show that knocking down LsFer1 and LsFer2 both negatively affect the parasite's physiology, as it is not able to properly feed and reproduce. In a starvation experiment, the transcript levels of both LsFer1 and LsFer2 decrease during the starvation period. Combined, these results demonstrate the importance of these genes for the normal parasite biology, and they could thus potentially be targets for pest management.

Sex-biased gene expression and evolution of candidate reproductive transcripts in adult stages of salmon lice (Lepeophtheirus salmonis)

The salmon louse Lepeophtheirus salmonis (Krøyer 1837) displays numerous sexually dimorphic characteristics. Insights into their underlying molecular components have only recently been explored, which serve to better understand both the basic biology of the louse, associated impacts on drug sensitivity, and evolution of resistance. Expression of 16 L. salmonis genes putatively involved in sexual dimorphism and reproduction were used to determine differences between sexes and better understand responses to mating using RT-qPCR of pre-adult and adult lice. Analysis of these genes revealed the dynamic nature of sex-biased expression across stages. However, female reception of a spermatophore did not appear to impact the expression of these particular genes. Furthermore six of these transcripts and 84 others identified previously in a large-scale louse transcriptomics experiment were used to estimate differences in evolutionary rates and codon-usage bias of sex-related genes using phylogenetic analysis by maximum likelihood (PAML) and codonW. Results suggest male-biased genes are evolving at significantly greater rates than female-biased and unbiased genes as evidenced by higher rates of non-synonymous substitutions and lower codon-usage bias in these genes. These analyses expand our understanding of interactions of sex-biased expression across the pre-adult and adult life stages and provide foundations for better understanding evolutionary differences in sex-biased genes of L. salmonis.

Norwegian e-Infrastructure for Life Sciences (NeLS)

The Norwegian e-Infrastructure for Life Sciences (NeLS) has been developed by ELIXIR Norway to provide its users with a system enabling data storage, sharing, and analysis in a project-oriented fashion. The system is available through easy-to-use web interfaces, including the Galaxy workbench for data analysis and workflow execution. Users confident with a command-line interface and programming may also access it through Secure Shell (SSH) and application programming interfaces (APIs).  NeLS has been in production since 2015, with training and support provided by the help desk of ELIXIR Norway. Through collaboration with NorSeq, the national consortium for high-throughput sequencing, an integrated service is offered so that sequencing data generated in a research project is provided to the involved researchers through NeLS. Sensitive data, such as individual genomic sequencing data, are handled using the TSD (Services for Sensitive Data) platform provided by Sigma2 and the University of Oslo. NeLS integrates national e-infrastructure storage and computing resources, and is also integrated with the SEEK platform in order to store large data files produced by experiments described in SEEK.   In this article, we outline the architecture of NeLS and discuss possible directions for further development.