Genomic resources for sea lice: analysis of ESTs and mitochondrial genomes

Caligus minimus Otto, 1821 (Copepoda: Caligidae): A commercially important but poorly described parasite of cultured European Sea Bass, Dicentrarchus labrax (Linnaeus, 1758)

Caligus minimus Otto, 1821 has been known for over two centuries and it is the second oldest of the approximately 275 species of Caligus O. F. Müller, 1985. Despite the numerous records of this species from European waters, it has never been fully described to modern standards. The lack of a comprehensive modern description has resulted in numerous misidentifications, even in recently published reports, and this is especially problematic for a species that is known to have a significant economic impact in aquaculture. This study presents a detailed description of both sexes and documents newly observed features of C. minimus collected from the buccal cavity of farmed European Sea Bass (ESB), Dicentrarchus labrax (Linnaeus, 1758). The morphology of C. minimus was examined using light microscope (LM), scanning electron microscope (SEM), and confocal laser scanning microscope (CLSM), and new details are revealed regarding the structure and ornamentation of the marginal membrane of the cephalothorax, maxilliped, antenna, sternal furca, abdomen, and legs 1, 3, 4, and 6. The ornamentation of the marginal membrane of the cephalothorax is unique and its impact on the functioning of the cephalothoracic sucker requires further investigation. Additionally, partial COI gene region sequences were obtained from four individuals of C. minimus and provided for future references. A phylogenetic analysis was conducted in conjunction with Caligus sequences available in the NCBI GenBank database.

Insights from Hi-C data regarding the Pacific salmon louse (Lepeophtheirus salmonis) sex chromosomes

Salmon lice, Lepeophtheirus salmonis (family Caligidae), are ectoparasites that have negatively impacted the salmon aquaculture industry and vulnerable wild salmon populations. Researchers have studied salmon lice to better understand their biology to develop effective control strategies. In this study, we updated the chromosome-level reference genome assembly of the Pacific subspecies of L. salmonis using Hi-C data. The previous version placed contigs/scaffolds using an Atlantic salmon louse genetic map. By utilizing Hi-C data from Pacific salmon lice, we were able to properly assign locations to contigs/scaffolds previously unplaced or misplaced. This resulted in a more accurate genome assembly and a more comprehensive characterization of the sex chromosome unique to females (W). We found evidence that the same ZW-ZZ mechanism is common in both Atlantic and Pacific subspecies of salmon lice using PCR assays. The W chromosome was approximately 800 kb in size, which is ∼30 times smaller than the Z chromosome (24 Mb). The W chromosome contained 61 annotated genes, including 32 protein-coding genes, 27 long noncoding RNA (lncRNA) genes, and 2 pseudogenes. Among these 61 genes, 39 genes shared homology to genes found on other chromosomes, while 20 were unique to the W chromosome. Two genes of interest on the W chromosome, prohibitin-2 and kinase suppressor of ras-2, were previously identified as potential sex-linked markers in the salmon louse. However, we prioritized the 20 unique genes on the W chromosome as sex-determining candidates. This information furthers our understanding of the biology of this ectoparasite and may help in the development of more effective management strategies.

Complete Mitogenome sequencing of the fish louse Argulus japonicus (Crustacea: Branchiura): Comparative analyses and phylogenetic implications

The fish louse Argulus japonicus, a branchiuran crustacean of the Argulidae family, is attracting increasing attention because of its parasitic tendencies and significant health threats to global fish farming. The mitogenomes can yield a foundation for studying epidemiology, genetic diversity, and molecular ecology and therefore may be used to assist in the surveillance and control of A. japonicus. In this study, we sequenced and assembled the complete mitogenome of A. japonicus to shed light on its genetic and evolutionary blueprint. Our investigation indicated that the 15,045-bp circular genome of A. japonicus encodes 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs) with significant AT and GC skews. Comparative genomics provided an evolutionary scenario for the genetic diversity of 13 PCGs: all were under purifying selection, with cox1 and nad6 having the lowest and highest evolutionary rates, respectively. Genome-wide phylogenetic trees established a close relationship between species of the families Argulidae (Arguloida) and Armilliferidae (Porocephalida) within Crustacea, and further, A. japonicus and Argulus americanus were determined to be more closely related to each other than to others within the family Argulidae. Single PCG-based phylogenies supported nad1 and nad6 as the best genetic markers for evolutionary and phylogenetic studies for branchiuran crustaceans due to their similar phylogenetic topologies with those of genome-based phylogenetic analyses. To sum up, these comprehensive mitogenomic data of A. japonicus and related species refine valuable marker resources and should contribute to molecular diagnostic methods, epidemiological investigations, and ecological studies of the fish ectoparasites in Crustacea.

Divergent RNA viruses infecting sea lice, major ectoparasites of fish

Sea lice, the major ectoparasites of fish, have significant economic impacts on wild and farmed finfish, and have been implicated in the decline of wild salmon populations. As blood-feeding arthropods, sea lice may also be reservoirs for viruses infecting fish. However, except for two groups of negative-strand RNA viruses within the order Mononegavirales, nothing is known about viruses of sea lice. Here, we used transcriptomic data from three key species of sea lice (Lepeophtheirus salmonis, Caligus clemensi, and Caligus rogercresseyi) to identify 32 previously unknown RNA viruses. The viruses encompassed all the existing phyla of RNA viruses, with many placed in deeply branching lineages that likely represent new families and genera. Importantly, the presence of canonical virus-derived small interfering RNAs (viRNAs) indicates that most of these viruses infect sea lice, even though in some cases their closest classified relatives are only known to infect plants or fungi. We also identified both viRNAs and PIWI-interacting RNAs (piRNAs) from sequences of a bunya-like and two qin-like viruses in C. rogercresseyi. Our analyses showed that most of the viruses found in C. rogercresseyi occurred in multiple life stages, spanning from planktonic to parasitic stages. Phylogenetic analysis revealed that many of the viruses infecting sea lice were closely related to those that infect a wide array of eukaryotes with which arthropods associate, including fungi and parasitic tapeworms, implying that over evolutionary time there has been cross-phylum and cross-kingdom switching of viruses between arthropods and other eukaryotes. Overall, this study greatly expands our view of virus diversity in crustaceans, identifies viruses that infect and replicate in sea lice, and provides evidence that over evolutionary time, viruses have switched between arthropods and eukaryotic hosts in other phyla and kingdoms.

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.

A genetic linkage map for the salmon louse (Lepeophtheirus salmonis): evidence for high male:female and inter-familial recombination rate differences

A salmon louse (Lepeophtheirus salmonis salmonis) genetic linkage map was constructed to serve as a genomic resource for future investigations into the biology of this important marine parasitic copepod species, and to provide insights into the inheritance patterns of genetic markers in this species. SNP genotyping of 8 families confirmed the presence of 15 linkage groups based upon the assignment of 93,773 markers. Progeny sample size weight adjusted map sizes in males (with the exception of SL12 and SL15) ranged in size from 96.50 cM (SL11) to 134.61 cM (SL06), and total combined map steps or bins ranged from 143 (SL09) to 203 (SL13). The SL12 male map was the smallest linkage group with a weight-averaged size of 3.05 cM with 6 recombination bins. Male:female specific recombination rate differences are 10.49:1 and represent one of the largest reported sex-specific differences for any animal species. Recombination ratio differences (M:F) ranged from 1.0 (SL12) to 29:1 (SL15). The number of markers exhibiting normal Mendelian segregation within the sex linkage group SL15 was extremely low (N = 80) in comparison to other linkage groups genotyped [range: 1459 (SL12)-10206 markers (SL05)]. Re-evaluation of Mendelian inheritance patterns of markers unassigned to any mapping parent according to hemizygous segregation patterns (models presented) identified matches for many of these markers to hemizygous patterns. The greatest proportion of these markers assigned to SL15 (N increased to 574). Inclusion of the hemizygous markers revised SL15 sex-specific recombination rate differences to 28:1. Recombination hot- and coldspots were identified across all linkage groups with all linkage groups possessing multiple peaks. Nine of 13 linkage groups evaluated possessed adjacent domains with hot-coldspot transitional zones. The most common pattern was for one end of the linkage to show elevated recombination in addition to internal regions. For SL01 and SL06, however, a terminal region with high recombination was not evident while a central domain possessing extremely high-recombination levels was present. High levels of recombination were weakly coupled to higher levels of SNP variation within domains, but this association was very strong for the central domains of SL01 and SL06. From the pooled paternal half-sib lots (several virgin females placed with 1 male), only 1 or two surviving family lots were obtained. Surviving families possessed parents where both the male and female possessed either inherently low or high recombination rates. This study provides insight into the organization of the sea louse genome, and describes large differences in recombination rate that exist among individuals of the same sex, and between the sexes. These differences in recombination rate may be coupled to the capabilities of this species to adapt to environmental and pharmaceutical treatments, given that family survivorship appears to be enhanced when parents have similar recombination levels.

High level efficacy of lufenuron against sea lice (Lepeophtheirus salmonis) linked to rapid impact on moulting processes

Drug resistance in the salmon louse Lepeophtheirus salmonis is a global issue for Atlantic salmon aquaculture. Multiple resistance has been described across most available compound classes with the exception of the benzoylureas. To target this gap in effective management of L. salmonis and other species of sea lice (e.g. Caligus spp.), Elanco Animal Health is developing an in-feed treatment containing lufenuron (a benzoylurea) to be administered prior to seawater transfer of salmon smolts and to provide long-term protection of salmon against sea lice infestations. Benzoylureas disrupt chitin synthesis, formation, and deposition during all moulting events. However, the mechanism(s) of action are not yet fully understood and most research completed to date has focused on insects. We exposed the first parasitic stage of L. salmonis to 700 ppb lufenuron for three hours and observed over 90% reduction in survival to the chalimus II life stage on the host, as compared to vehicle controls. This agrees with a follow up in vivo administration study on the host, which showed >95% reduction by the chalimus I stage. Transcriptomic responses of salmon lice exposed to lufenuron included genes related to moulting, epithelial differentiation, solute transport, and general developmental processes. Global metabolite profiles also suggest that membrane stability and fluidity is impacted in treated lice. These molecular signals are likely the underpinnings of an abnormal moulting process and cuticle formation observed ultrastructurally using transmission electron microscopy. Treated nauplii-staged lice exhibited multiple abnormalities in the integument, suggesting that the coordinated assembly of the epi- and procuticle is impaired. In all cases, treatment with lufenuron had rapid impacts on L. salmonis development. We describe multiple experiments to characterize the efficacy of lufenuron on eggs, larvae, and parasitic stages of L. salmonis, and provide the most comprehensive assessment of the physiological responses of a marine arthropod to a benzoylurea chemical.

Characterisation of proteins in excretory/secretory products collected from salmon lice, Lepeophtheirus salmonis

BACKGROUND

The salmon louse, Lepeophtheirus salmonis, is an ectoparasitic copepod which feeds on the mucus, skin and blood of salmonid fish species. The parasite can persist on the surface of the fish without any effective control being exerted by the host immune system. Other ectoparasitic invertebrates produce compounds in their saliva, excretions and/or secretions which modulate the host immune responses allowing them to remain on or in the host during development. Similarly, compounds are produced in secretions of L. salmonis which are thought to be responsible for immunomodulation of the host responses as well as other aspects of crucial host-parasite interactions.

METHODS

In this study we have identified and characterised the proteins in the excretory/secretory (E/S) products of L. salmonis using LC-ESI-MS/MS.

RESULTS

In total 187 individual proteins were identified in the E/S collected from adult lice and pre-adult sea lice. Fifty-three proteins, including 13 serine-type endopeptidases, 1 peroxidase and 5 vitellogenin-like proteins were common to both adult and pre-adult E/S products. One hundred and seven proteins were identified in the adult E/S but not in the pre-adult E/S and these included serine and cysteine-type endopeptidases, vitellogenins, sphingomyelinase and calreticulin. A total of 27 proteins were identified in pre-adult E/S products but not in adult E/S.

CONCLUSIONS

The assigned functions of these E/S products and the potential roles they play in host-parasite interaction is discussed.

Mitochondrial genomes of the key zooplankton copepods Arctic Calanus glacialis and North Atlantic Calanus finmarchicus with the longest crustacean non-coding regions

We determined the nearly complete mitochondrial genomes of the Arctic Calanus glacialis and its North Atlantic sibling Calanus finmarchicus, which are key zooplankton components in marine ecosystems. The sequenced part of C. glacialis mitogenome is 27,342 bp long and consists of two contigs, while for C. finmarchicus it is 29,462 bp and six contigs, what makes them the longest reported copepod mitogenomes. The typical set of metazoan mitochondrial genes is present in these mitogenomes, although the non-coding regions (NCRs) are unusually long and complex. The mitogenomes of the closest species C. glacialis and C. finmarchicus, followed by the North Pacific C. sinicus, are structurally similar and differ from the much more typical of deep-water, Arctic C. hyperboreus. This evolutionary trend for the expansion of NCRs within the Calanus mitogenomes increases mitochondrial DNA density, what resulted in its similar density to the nuclear genome. Given large differences in the length and structure of C. glacialis and C. finmarchicus mitogenomes, we conclude that the species are genetically distinct and thus cannot hybridize. The molecular resources presented here: the mitogenomic and rDNA sequences, and the database of repetitive elements should facilitate the development of genetic markers suitable in pursuing evolutionary research in copepods.

Circadian signaling in the Northern krill Meganyctiphanes norvegica: In silico prediction of the protein components of a putative clock system using a publicly accessible transcriptome

The Northern krill Meganyctiphanes norvegica is a significant component of the zooplankton community in many regions of the North Atlantic Ocean. In the areas it inhabits, M. norvegica is of great importance ecologically, as it is both a major consumer of phytoplankton/small zooplankton and is a primary food source for higher-level consumers. One behavior of significance for both feeding and predator avoidance in Meganyctiphanes is diel vertical migration (DVM), i.e., a rising from depth at dusk and a return to depth at dawn. In this and other euphausiids, an endogenous circadian pacemaker is thought, at least in part, to control DVM. Currently, there is no information concerning the identity of the genes/proteins that comprise the M. norvegica circadian system. In fact, there is little information concerning the molecular underpinnings of circadian rhythmicity in crustaceans generally. Here, a publicly accessible transcriptome was used to identify the molecular components of a putative Meganyctiphanes circadian system. A complete set of core clock proteins was deduced from the M. norvegica transcriptome (clock, cryptochrome 2, cycle, period and timeless), as was a large suite of proteins that likely function as modulators of the core clock (e.g., doubletime), or serves as inputs to it (cryptochrome 1) or outputs from it (pigment dispersing hormone). This is the first description of a "complete" (core clock through putative output pathway signals) euphausiid clock system, and as such, provides a foundation for initiating molecular investigations of circadian signaling in M. norvegica and other krill species, including how clock systems may regulate DVM and other behaviors.

Effects of the vertically transmitted microsporidian Facilispora margolisi and the parasiticide emamectin benzoate on salmon lice (Lepeophtheirus salmonis)

Background

Microsporidia are highly specialized, parasitic fungi that infect a wide range of eukaryotic hosts from all major taxa. Infections cause a variety of damaging effects on host physiology from increased stress to death. The microsporidian Facilispora margolisi infects the Pacific salmon louse (Lepeophtheirus salmonis oncorhynchi), an economically and ecologically important ectoparasitic copepod that can impact wild and cultured salmonids.

Results

Vertical transmission of F. margolisi was demonstrated by using PCR and in situ hybridization to identify and localize microsporidia in female L. salmonis and their offspring. Spores and developmental structures of F. margolisi were identified in 77% of F₁ generation copepods derived from infected females while offspring from uninfected females all tested negative for the microsporidia. The transcriptomic response of the salmon louse to F. margolisi was profiled at both the copepodid larval stage and the pre-adult stage using microarray technology. Infected copepodids differentially expressed 577 transcripts related to stress, ATP generation and structural components of muscle and cuticle. The infection also impacted the response of the copepodid to the parasiticide emamectin benzoate (EMB) at a low dose of 1.0 ppb for 24 h. A set of 48 transcripts putatively involved in feeding and host immunomodulation were up to 8-fold underexpressed in the F. margolisi infected copepodids treated with EMB compared with controls or either stressor alone. Additionally, these infected lice treated with EMB also overexpressed 101 transcripts involved in stress resistance and signalling compared to the other groups. In contrast, infected pre-adult lice did not display a stress response, suggesting a decrease in microsporidian virulence associated with lice maturity. Furthermore, copepodid infectivity and moulting was not affected by the microsporidian infection.

Conclusions

This study demonstrated that F. margolisi is transmitted vertically between salmon louse generations and that biological impacts of infection differ depending on the stage of the copepod host. The infection caused significant perturbations of larval transcriptomes and therefore must be considered in future studies in which impacts to host development and environmental factors are assessed. Fitness impacts are probably minor, although the interaction between pesticide exposure and microsporidian infection merits further study.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4040-8) contains supplementary material, which is available to authorized users.

Host–parasite transcriptomics during immunostimulant-enhanced rejection of salmon lice (Lepeophtheirus salmonis) by Atlantic salmon (Salmo salar)

Salmon lice ( Lepeophtheirus salmonis) are important ectoparasites of wild and farmed salmonids and cause major losses to the salmon farming industry throughout the Northern Hemisphere. With the emergence of resistance to several commonly used parasiticides, novel control strategies and integration of multiple treatment options are needed, including host immunostimulation. Here, we investigate the effects of a functional feed containing a peptidoglycan and nucleotide formulation on L. salmonis infection of Atlantic salmon ( Salmo salar) by characterizing lice infection levels, the expression of several host immune genes, and the parasite transcriptomic response to the immunostimulated host. Although initial infection intensities were low, the low dose (LD) immunostimulant diet reduced the total lice burden by 50% relative to controls. Immunostimulant fed hosts up-regulated interleukin-1β in the skin and spleen. This gene has been implicated in successful responses of several salmonid species to salmon lice but is typically not observed in Atlantic salmon, suggesting a favorable influence on the immune response. Lice infecting LD immunostimulated salmon overexpressed genes putatively involved in parasite immunity, including carboxylesterases, and underexpressed genes putatively involved in feeding (e.g., proteases). These lice response genes further improve the characterization of the transcriptome of the non-model parasite by identifying genes potentially involved in evading host immunity.

Enhanced transcriptomic responses in the Pacific salmon louse Lepeophtheirus salmonis oncorhynchi to the non-native Atlantic Salmon Salmo salar suggests increased parasite fitness

Background

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

Results

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

Conclusions

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

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3520-1) contains supplementary material, which is available to authorized users.

Cypermethrin exposure induces metabolic and stress-related gene expression in copepodid salmon lice (Lepeophtheirus salmonis)

Cypermethrin has been administered for decades to control salmon lice (Lepeophtheirus salmonis) infestations in Atlantic salmon farming regions globally. However, resistance to cypermethrin and other available therapeutants has threatened the sustainability of this growing industry. To better understand the effects of cypermethrin on L. salmonis, a 38K oligonucleotide microarray and RT-qPCR analyses were applied to pools of copepodid larvae exposed to 1.0ppb cypermethrin or seawater controls for 24h. Phenotypic assessments and global gene expression profiles showed a significant disruption of homeostasis in copepodid L. salmonis exposed to cypermethrin. Multiple degradative enzymes were overexpressed in cypermethrin-treated lice including five trypsin-like serine proteases and three cytochrome p450s CYP3a24 (p=0.03, fold change (FC)=3.8; GenBank accession no. JP326960.1), CYP6w1 (p=0.008, FC=5.3; GenBank accession no. JP317875.1), and CYP6d4 (p=0.01; FC=7.9; GenBank accession no. JP334550.1). These enzymes represent preliminary markers for understanding the physiological response of L. salmonis to cypermethrin exposure. A general stress response was also observed in cypermethrin-treated lice which included differential expression of cell signaling genes involved in the induction of cell growth, solute transport, and metabolism. Lastly, a consensus-based analysis was completed with two previously published L. salmonis transcriptome studies revealing genes that respond to cypermethrin, emamectin benzoate (another delousing agent) and hyposalinity. This included concordant differential expression of heat shock beta-1, ammonium transporter Rh types B, and 72kDa type IV collagenase across different L. salmonis studies. This is currently the most comprehensive transcriptome assessment of chemical exposure on the first infectious stage of L. salmonis, providing novel markers for studying drug resistance and general stress in this important parasite.

Sex-biased gene expression and sequence conservation in Atlantic and Pacific salmon lice (Lepeophtheirus salmonis)

BACKGROUND

Salmon lice, Lepeophtheirus salmonis (Copepoda: Caligidae), are highly important ectoparasites of farmed and wild salmonids, and cause multi-million dollar losses to the salmon aquaculture industry annually. Salmon lice display extensive sexual dimorphism in ontogeny, morphology, physiology, behavior, and more. Therefore, the identification of transcripts with differential expression between males and females (sex-biased transcripts) may help elucidate the relationship between sexual selection and sexually dimorphic characteristics.

RESULTS

Sex-biased transcripts were identified from transcriptome analyses of three L. salmonis populations, including both Atlantic and Pacific subspecies. A total of 35-43 % of all quality-filtered transcripts were sex-biased in L. salmonis, with male-biased transcripts exhibiting higher fold change than female-biased transcripts. For Gene Ontology and functional analyses, a consensus-based approach was used to identify concordantly differentially expressed sex-biased transcripts across the three populations. A total of 127 male-specific transcripts (i.e. those without detectable expression in any female) were identified, and were enriched with reproductive functions (e.g. seminal fluid and male accessory gland proteins). Other sex-biased transcripts involved in morphogenesis, feeding, energy generation, and sensory and immune system development and function were also identified. Interestingly, as observed in model systems, male-biased L. salmonis transcripts were more frequently without annotation compared to female-biased or unbiased transcripts, suggesting higher rates of sequence divergence in male-biased transcripts.

CONCLUSIONS

Transcriptome differences between male and female L. salmonis described here provide key insights into the molecular mechanisms controlling sexual dimorphism in L. salmonis. This analysis offers targets for parasite control and provides a foundation for further analyses exploring critical topics such as the interaction between sex and drug resistance, sex-specific factors in host-parasite relationships, and reproductive roles within L. salmonis.

Intraspecific Variation in Mitogenomes of Five Crassostrea Species Provides Insight into Oyster Diversification and Speciation

A large number of Crassostrea oysters are found in Asia-Pacific. While analyses of interspecific variation have helped to establish historical relationships among these species, studies on intraspecific variation are necessary to understand their recent evolutionary history and current forces driving population biology. We resequenced 18 and analyzed 31 mitogenomes of five Crassostrea species from China: Crassostrea gigas, Crassostrea angulata, Crassostrea sikamea, Crassostrea ariakensis, and Crassostrea hongkongensis. Our analysis finds abundant insertions, deletions, and single-nucleotide polymorphisms in all species. Intraspecific variation varies greatly among species with polymorphic sites ranging from 54 to 293 and nucleotide diversity ranging from 0.00106 to 0.00683. In all measurements, C. hongkongensis that has the narrowest geographic distribution exhibits the least sequence diversity; C. ariakensis that has the widest distribution shows the highest diversity, and species with intermediate distribution show intermediate levels of diversity. Low sequence diversity in C. hongkongensis may reflect recent bottlenecks that are probably exacerbated by human transplantation. High diversity in C. ariakensis is likely due to divergence of northern and southern China populations that have been separated without gene flow. The significant differences in mitogenome diversity suggest that the five sister species of Crassostrea have experienced different evolutionary forces since their divergence. The recent divergence of two C. ariakensis populations and the C. gigas/angulata species complex provides evidence for continued diversification and speciation of Crassostrea species along China's coast, which are shaped by unknown mechanisms in a north-south divide.

Next-Generation Transcriptome Profiling of the Salmon Louse Caligus rogercresseyi Exposed to Deltamethrin (AlphaMax™): Discovery of Relevant Genes and Sex-Related Differences

Sea lice are one of the main parasites affecting the salmon aquaculture industry, causing significant economic losses worldwide. Increased resistance to traditional chemical treatments has created the need to find alternative control methods. Therefore, the objective of this study was to identify the transcriptome response of the salmon louse Caligus rogercresseyi to the delousing drug deltamethrin (AlphaMax™). Through bioassays with different concentrations of deltamethrin, adult salmon lice transcriptomes were sequenced from cDNA libraries in the MiSeq Illumina platform. A total of 78 million reads for females and males were assembled in 30,212 and 38,536 contigs, respectively. De novo assembly yielded 86,878 high-quality contigs and, based on published data, it was possible to annotate and identify relevant genes involved in several biological processes. RNA-seq analysis in conjunction with heatmap hierarchical clustering evidenced that pyrethroids modify the ectoparasitic transcriptome in adults, affecting molecular processes associated with the nervous system, cuticle formation, oxidative stress, reproduction, and metabolism, among others. Furthermore, sex-related transcriptome differences were evidenced. Specifically, 534 and 1033 exclusive transcripts were identified for males and females, respectively, and 154 were shared between sexes. For males, estradiol 17-beta-dehydrogenase, sphingolipid delta4-desaturase DES1, ketosamine-3-kinase, and arylsulfatase A, among others, were discovered, while for females, vitellogenin 1, glycoprotein G, transaldolase, and nitric oxide synthase were among those identified. The shared transcripts included annotations for tropomyosin, γ-crystallin A, glutamate receptor-metabotropic, glutathione S-transferase, and carboxipeptidase B. The present study reveals that deltamethrin generates a complex transcriptome response in C. rogercresseyi, thus providing valuable genomic information for developing new delousing drugs.

Characterization of a novel RXR receptor in the salmon louse (Lepeophtheirus salmonis, Copepoda) regulating growth and female reproduction

BACKGROUND

Nuclear receptors have crucial roles in all metazoan animals as regulators of gene transcription. A wide range of studies have elucidated molecular and biological significance of nuclear receptors but there are still a large number of animals where the knowledge is very limited. In the present study we have identified an RXR type of nuclear receptor in the salmon louse (Lepeophtheirus salmonis) (i.e. LsRXR). RXR is one of the two partners of the Ecdysteroid receptor in arthropods, the receptor for the main molting hormone 20-hydroxyecdysone (E20) with a wide array of effects in arthropods.

RESULTS

Five different LsRXR transcripts were identified by RACE showing large differences in domain structure. The largest isoforms contained complete DNA binding domain (DBD) and ligand binding domain (LBD), whereas some variants had incomplete or no DBD. LsRXR is transcribed in several tissues in the salmon louse including ovary, subcuticular tissue, intestine and glands. By using Q-PCR it is evident that the LsRXR mRNA levels vary throughout the L. salmonis life cycle. We also show that the truncated LsRXR transcript comprise about 50% in all examined samples. We used RNAi to knock-down the transcription in adult reproducing female lice. This resulted in close to zero viable offspring. We also assessed the LsRXR RNAi effects using a L. salmonis microarray and saw significant effects on transcription in the female lice. Transcription of the major yolk proteins was strongly reduced by knock-down of LsRXR. Genes involved in lipid metabolism and transport were also down regulated. Furthermore, different types of growth processes were up regulated and many cuticle proteins were present in this group.

CONCLUSIONS

The present study demonstrates the significance of LsRXR in adult female L. salmonis and discusses the functional aspects in relation to other arthropods. LsRXR has a unique structure that should be elucidated in the future.

Transcriptomic responses to emamectin benzoate in Pacific and Atlantic Canada salmon lice Lepeophtheirus salmonis with differing levels of drug resistance

Salmon lice Lepeophtheirus salmonis are an ecologically and economically important parasite of wild and farmed salmon. In Scotland, Norway, and Eastern Canada, L. salmonis have developed resistance to emamectin benzoate (EMB), one of the few parasiticides available for salmon lice. Drug resistance mechanisms can be complex, potentially differing among populations and involving multiple genes with additive effects (i.e., polygenic resistance). Indicators of resistance development may enable early detection and countermeasures to avoid the spread of resistance. Here, we collect sensitive Pacific L. salmonis and sensitive and resistant Atlantic L. salmonis from salmon farms, propagate in laboratory (F1), expose to EMB in bioassays, and evaluate either baseline (Atlantic only) or induced transcriptomic differences between populations. In all populations, induced responses were minor and a cellular stress response was not identified. Pacific lice did not upregulate any genes in response to EMB, but downregulated degradative enzymes and transport proteins at 50 ppb EMB. Baseline differences between sensitive and now resistant Atlantic lice were much greater than responses to exposures. All resistant lice overexpressed degradative enzymes, and resistant males, the most resistant group, overexpressed collagenases to the greatest extent. These results indicate an accumulation of baseline expression differences related to resistance.

Comprehensive transcriptome study to develop molecular resources of the copepod Calanus sinicus for their potential ecological applications

Calanus sinicus Brodsky (Copepoda, Crustacea) is a dominant zooplanktonic species widely distributed in the margin seas of the Northwest Pacific Ocean. In this study, we utilized an RNA-Seq-based approach to develop molecular resources for C. sinicus. Adult samples were sequenced using the Illumina HiSeq 2000 platform. The sequencing data generated 69,751 contigs from 58.9 million filtered reads. The assembled contigs had an average length of 928.8 bp. Gene annotation allowed the identification of 43,417 unigene hits against the NCBI database. Gene ontology (GO) and KEGG pathway mapping analysis revealed various functional genes related to diverse biological functions and processes. Transcripts potentially involved in stress response and lipid metabolism were identified among these genes. Furthermore, 4,871 microsatellites and 110,137 single nucleotide polymorphisms (SNPs) were identified in the C. sinicus transcriptome sequences. SNP validation by the melting temperature (T_m)-shift method suggested that 16 primer pairs amplified target products and showed biallelic polymorphism among 30 individuals. The present work demonstrates the power of Illumina-based RNA-Seq for the rapid development of molecular resources in nonmodel species. The validated SNP set from our study is currently being utilized in an ongoing ecological analysis to support a future study of C. sinicus population genetics.

RNA-Seq analysis using de novo transcriptome assembly as a reference for the salmon louse Caligus rogercresseyi

Despite the economic and environmental impacts that sea lice infestations have on salmon farming worldwide, genomic data generated by high-throughput transcriptome sequencing for different developmental stages, sexes, and strains of sea lice is still limited or unknown. In this study, RNA-seq analysis was performed using de novo transcriptome assembly as a reference for evidenced transcriptional changes from six developmental stages of the salmon louse Caligus rogercresseyi. EST-datasets were generated from the nauplius I, nauplius II, copepodid and chalimus stages and from female and male adults using MiSeq Illumina sequencing. A total of 151,788,682 transcripts were yielded, which were assembled into 83,444 high quality contigs and subsequently annotated into roughly 24,000 genes based on known proteins. To identify differential transcription patterns among salmon louse stages, cluster analyses were performed using normalized gene expression values. Herein, four clusters were differentially expressed between nauplius I–II and copepodid stages (604 transcripts), five clusters between copepodid and chalimus stages (2,426 transcripts), and six clusters between female and male adults (2,478 transcripts). Gene ontology analysis revealed that the nauplius I–II, copepodid and chalimus stages are mainly annotated to aminoacid transfer/repair/breakdown, metabolism, molting cycle, and nervous system development. Additionally, genes showing differential transcription in female and male adults were highly related to cytoskeletal and contractile elements, reproduction, cell development, morphogenesis, and transcription-translation processes. The data presented in this study provides the most comprehensive transcriptome resource available for C. rogercresseyi, which should be used for future genomic studies linked to host-parasite interactions.

The phylogenetic distribution and evolution of enzymes within the thymidine kinase 2-like gene family in metazoa

Deoxyribonucleoside kinases (dNKs) carry out the rate-determining step in the nucleoside salvage pathway within all domains of life where the pathway is present, and, hence, are an indication on whether or not a species/genus retains the ability to salvage deoxyribonucleosides. Here, a phylogenetic tree is constructed for the thymidine kinase 2-like dNK gene family in metazoa. Each enzyme class (deoxycytidine, deoxyguanosine, and deoxythymidine kinases, as well as the multisubstrate dNKs) falls into a monophyletic clade. However, in vertebrates, dCK contains an apparent duplication with one paralog lost in mammals, and a number of crustacean genomes (like Caligus rogercresseyi and Lepeophtheirus salmonis) unexpectedly contain not only the multisubstrate dNKs, related to Drosophila multisubstrate dNK, but also a TK2-like kinase. Additionally, crustaceans (Daphnia, Caligus, and Lepeophtheirus) and some insects (Tribolium, Danaus, Pediculus, and Acyrthosiphon) contain several multisubstrate dNK-like enzymes which group paraphyletically within the arthropod clade. This might suggest that the multisubstrate dNKs underwent multiple rounds of duplications with differential retention of duplicate copies between insect families and more complete retention within some crustaceans and insects. Genomes of several basal animalia contain more than one dNK-like sequence, some of which group outside the remaining eukaryotes (both plants and animals) and/or with bacterial dNKs. Within the vertebrates, the mammalian genomes do not contain the second dCK, while birds, fish, and amphibians do retain it. Phasianidae (chicken and turkey) have lost dGK, while it has been retained in other bird lineages, like zebra finch. Reconstruction of the ancestral sequence between the multisubstrate arthropod dNKs and the TK2 clade of vertebrates followed by homology modeling and discrete molecular dynamics calculations on this sequence were performed to examine the evolutionary path which led to the two different enzyme classes. The structural models showed that the carboxyl terminus of the ancestral sequence is more helical than dNK, in common with TK2, although any implications of this for enzyme specificity will require biochemical validation. Finally, rate-shift and conservation-shift analysis between clades with different specificities uncovered candidate residues outside the active site pocket which may have contributed to differentiation in substrate specificity between enzyme clades.

Salmon lice (Lepeophtheirus salmonis) showing varying emamectin benzoate susceptibilities differ in neuronal acetylcholine receptor and GABA-gated chloride channel mRNA expression

Background

Caligid copepods, also called sea lice, are fish ectoparasites, some species of which cause significant problems in the mariculture of salmon, where the annual cost of infection is in excess of €300 million globally. At present, caligid control on farms is mainly achieved using medicinal treatments. However, the continued use of a restricted number of medicine actives potentially favours the development of drug resistance. Here, we report transcriptional changes in a laboratory strain of the caligid Lepeophtheirus salmonis (Krøyer, 1837) that is moderately (~7-fold) resistant to the avermectin compound emamectin benzoate (EMB), a component of the anti-salmon louse agent SLICE® (Merck Animal Health).

Results

Suppression subtractive hybridisation (SSH) was used to enrich transcripts differentially expressed between EMB-resistant (PT) and drug-susceptible (S) laboratory strains of L. salmonis. SSH libraries were subjected to 454 sequencing. Further L. salmonis transcript sequences were available as expressed sequence tags (EST) from GenBank. Contiguous sequences were generated from both SSH and EST sequences and annotated. Transcriptional responses in PT and S salmon lice were investigated using custom 15 K oligonucleotide microarrays designed using the above sequence resources. In the absence of EMB exposure, 359 targets differed in transcript abundance between the two strains, these genes being enriched for functions such as calcium ion binding, chitin metabolism and muscle structure. γ-aminobutyric acid (GABA)-gated chloride channel (GABA-Cl) and neuronal acetylcholine receptor (nAChR) subunits showed significantly lower transcript levels in PT lice compared to S lice. Using RT-qPCR, the decrease in mRNA levels was estimated at ~1.4-fold for GABA-Cl and ~2.8-fold for nAChR. Salmon lice from the PT strain showed few transcriptional responses following acute exposure (1 or 3 h) to 200 μg L^-1 of EMB, a drug concentration tolerated by PT lice, but toxic for S lice.

Conclusions

Avermectins are believed to exert their toxicity to invertebrates through interaction with glutamate-gated and GABA-gated chloride channels. Further potential drug targets include other Cys-loop ion channels such as nAChR. The present study demonstrates decreased transcript abundances of GABA-Cl and nAChR subunits in EMB-resistant salmon lice, suggesting their involvement in avermectin toxicity in caligids.

Complete mitochondrial genome of Concholepas concholepas inferred by 454 pyrosequencing and mtDNA expression in two mollusc populations

Despite the great relevance of mitochondrial genome analysis in evolutionary studies, there is scarce information on how the transcripts associated with the mitogenome are expressed and their role in the genetic structuring of populations. This work reports the complete mitochondrial genome of the marine gastropod Concholepas concholepas, obtained by 454 pryosequencing, and an analysis of mitochondrial transcripts of two populations 1000 km apart along the Chilean coast. The mitochondrion of C. concholepas is 15,495 base pairs (bp) in size and contains the 37 subunits characteristic of metazoans, as well as a non-coding region of 330 bp. In silico analysis of mitochondrial gene variability showed significant differences among populations. In terms of levels of relative abundance of transcripts associated with mitochondrion in the two populations (assessed by qPCR), the genes associated with complexes III and IV of the mitochondrial genome had the highest levels of expression in the northern population while transcripts associated with the ATP synthase complex had the highest levels of expression in the southern population. Moreover, fifteen polymorphic SNPs were identified in silico between the mitogenomes of the two populations. Four of these markers implied different amino acid substitutions (non-synonymous SNPs). This work contributes novel information regarding the mitochondrial genome structure and mRNA expression levels of C. concholepas.

Transcriptomics of coping strategies in free-swimming Lepeophtheirus salmonis (Copepoda) larvae responding to abiotic stress

The salmon louse Lepeophtheirus salmonis is a marine ectoparasite of wild and farmed salmon in the Northern Hemisphere. Infections of farmed salmon are of economic and ecological concern. Nauplius and copepodid salmon lice larvae are free-swimming and disperse in the water column until they encounter a host. In this study, we characterized the sublethal stress responses of L. salmonis copepodid larvae by applying a 38K oligonucleotide microarray to profile transcriptomes following 24 h exposures to suboptimal salinity (30–10 parts per thousand (‰)) or temperature (16–4 °C) environments. Hyposalinity exposure resulted in large-scale gene expression changes relative to those elicited by a thermal gradient. Subsequently, transcriptome responses to a more finely resolved salinity gradient between 30 ‰ and 25 ‰ were profiled. Minimal changes occurred at 29 ‰ or 28 ‰, a threshold of response was identified at 27 ‰, and the largest response was at 25 ‰. Differentially expressed genes were clustered by pattern of expression, and clusters were characterized by functional enrichment analysis. Results indicate larval copepods adopt two distinct coping strategies in response to short-term hyposaline stress: a primary response using molecular chaperones and catabolic processes at 27 ‰; and a secondary response up-regulating ion pumps, transporters, a different suite of chaperones and apoptosis-related transcripts at 26 ‰ and 25 ‰. The results further our understanding of the tolerances of L. salmonis copepodids to salinity and temperature gradients and may assist in the development of salmon louse management strategies.