Population and Transcriptomic Changes of the Tropical Fish Parasite Caligus confusus (Copepoda: Caligidae) with Seasonal Variations in Sea Temperature

Abstract

Fish–parasite systems could be subject to two scenarios under climate change: (i) increased water temperature might enhance parasite metabolism, allowing the parasite to spread rapidly; (ii) parasitism could decrease if the optimal temperature for growth and transmission is exceeded. Sea lice are parasitic copepods commonly found on marine fish in tropical regions, yet their biology remains poorly investigated. In this study, we analyzed the changes in infection levels and the transcriptomic response of the tropical sea louse Caligus confusus to two seasonal seawater temperatures (30 °C, “warm”, and 21 °C, “cold”). The prevalence of C. confusus was significantly higher in the colder water. A de novo transcriptomic analysis of C. confusus, the first for a tropical sea louse, revealed 426 over-expressed and 1402 down-expressed transcripts at the lower temperature. In particular, we observed over-expression of transcripts encoding vitellogenins (vit-1, vit-2, vit-4, and vit-6) and matrix metalloproteinases (mmp-2 and mmp-9), which are involved in reproduction and development. These results suggest that the cold tropical season physiologically favors C. confusus and that low temperature favors embryo development, which might ultimately lead to a higher prevalence. It is possible, therefore, that climate change could reduce some tropical sea lice populations during extreme warming events.

1. Introduction

Copepods are a diverse group of crustaceans composed of the free-living species that are the primary components of marine zooplankton, and symbiotic species, most of which have a direct life cycle and are commonly ectoparasitic on fish, although they have also been found in practically all marine animals [1]. Members of the family Caligidae, which includes 30 genera and more than 450 species, also called sea lice, are the most common parasitic copepods in marine fishes [2]. The most speciose genera in this family are Caligus and Lepeophtheirus, comprising 278 and 124 species, respectively [1]. Sea lice, mainly Caligus rogercresseyi and Lepeophtheirus salmonis, have become widely known due to their negative impact on fish aquaculture, causing skin irritation, ulcerations, anemia, lethargy, weight loss, secondary infections, and mortality in fish [3,4,5].

Variation in water temperature influences sea lice population dynamics. For instance, in temperate waters the prevalence of L. salmonis and C. elongatus can increase in spring (12–19 °C) and decline in winter (2–6 °C) [6]. This is because in warmer waters, the reproduction rate of sea lice is faster and the development time is shorter than in colder waters [7,8,9]. Conversely, some tropical species (e.g., L. simplex) thrive in warmer (30 °C) conditions [10], but others (e.g., C. omissus) do so in cooler (22 °C) environments [11]. Thus, global climate change can also be expected to perturb fish–parasite systems, initially by amplifying parasite metabolism, resulting in their more rapid spread; however, if the optimal temperature for growth and transmission is exceeded, parasite infections will decrease [12].

These general scenarios for parasitic copepods, especially at tropical latitudes, have been scarcely investigated, and the molecular mechanisms involved are not well understood. Whereas transcriptomic research has shed some light on the complex biological processes of C. rogercresseyi and L. salmonis [13,14,15,16], to date, transcriptomes are not available for tropical sea lice.

The aims of this study were to perform a de novo transcriptome assembly of C. confusus, and to document changes in infection levels and the transcriptomic responses of this parasite at contrasting seawater temperatures. This sea louse species is commonly found in fish of the family Carangidae (jacks) from the tropical waters of the Eastern Pacific and Indo-Pacific regions [17,18]. In particular, we focused on genes that are involved in reproduction and development, which might be associated with observable changes in the prevalence of C. confusus. Such data could be useful for better understanding the effect of climate change on sea lice, and parasitic copepods in general, in tropical regions.

2. Materials and Methods

2.1. Fish Sampling

Fish (Caranx caninus) were captured in the coastal waters of Mazatlán (23°12′48″ N, 106°26′20″ W), on the eastern side of the entrance to the Gulf of California (Mexico). The average sea surface temperature in this zone reaches a maximum of 31 °C between July and October and a minimum of 21 °C during January to March [19]. Thus, to determine the population and molecular changes occurring in C. confusus at contrasting temperatures, field samples were collected during October 2021 (warm season) and February 2022 (cold season).

2.2. Changes in Infection Levels

A total of 93 and 26 fish specimens were collected during October 2021 (average 30 °C) and February 2022 (average 21.8 °C), respectively. The fish were acquired directly from local fishermen while fresh. In the laboratory, the gills of the fish were removed and examined under a stereomicroscope (Motic, Richmond, BC, Canada) for the presence of sea lice, which were carefully transferred to labeled vials and preserved in 96% alcohol for later identification. Specimens of C. confusus were identified based on the morphological characteristics described by Ho and Lin [20]. To that end, copepods were cleared in lactic acid in order to observe their appendages in detail under a compound microscope (Leica DMLB, Leica Microsystems, Wetzlar, Germany). The prevalence and mean intensity of infection [21], with confidence intervals (CIs) at 95%, were calculated and compared between sampling months by using Fisher’s exact test and bootstrap t-testing, respectively, using the software Quantitative Parasitology (Qpweb 1.0.15) [22].

2.3. Transcriptomic Responses

2.3.1. Acquisition of Copepod Specimens

For molecular analysis, four fish (mean total length 56 cm) were captured on 25 October 2021 (seawater conditions: 30 °C, oxygen 90%, salinity 34.4 psu), and another four fish (mean total length 53) on 14 February 2022 (21 °C, oxygen 98%, salinity 34.6 psu). On these dates, the fish were caught together at the same site by local fishermen at a depth of 2 m, and immediately inspected for the purpose of recovering live sea lice. To that end, gills were removed, and sea lice were carefully collected under a stereomicroscope using sterile forceps. Specimens of C. confusus were transferred to a Petri dish containing filtered seawater (to remove excess mucus) for 1 min, preserved in RNAlater (Sigma–Aldrich) at room temperature for 24 h, and then stored at −20 °C until RNA extraction. Sixty female copepods were pooled for each sampling month, that is, October (warm) and February (cold).

2.3.2. RNA Extraction, Library Preparation, and Sequencing

The total RNA of each pool of copepods was extracted using Trizol reagent (Molecular Research Center, Inc., Cincinnati, OH, USA) according to the manufacturer’s instructions. The remaining genomic DNA was removed with a DNA-free kit (Life Technologies Corporation, South San Francisco, CA, USA). The quantity and quality of the total RNA were measured at 230, 260, and 280 nm with a Nano-Drop 2100 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). Total RNA was loaded onto a 1% agarose gel containing ethidium bromide. After electrophoresis, the gel was visualized using a Biorad Universal Hood II Gel Doc System (Bio-Rad Laboratories Segrate, Milan, Italy). The final quality of total RNA was verified by capillary electrophoresis using a 2100 Agilent Bioanalyzer (Agilent, Santa Clara, CA, USA) with a Nano Kit RNA chip. The RNA integrity number (RIN) was 7.5.

The RNA sequencing (RNA-seq) libraries were prepared using the MGI Easy RNA Library Prep Kit (MGI-Tech, MGI Tech Co., Ltd., Shenzhen, China). Sequencing of each library was performed on the MGSEQ-2000 platform (MGI Tech Co., Ltd., Shenzhen, China). Extraction of RNA, cDNA library construction, and sequencing were performed at the Genomic Services Laboratory (UGA-Langebio, CINVESTAV, Irapuato, Mexico). Raw sequence reads of C. confusus from cold and warm temperatures were deposited in the NCBI Short Read Archive (SRA) database under the accession numbers SRR25611633 and SRR25611632, respectively (BioProject PRJNA1004606).

2.3.3. De Novo Transcriptome Assembly and Functional Annotation

The quality of the raw RNA-seq library data was assessed with FastQC 0.11.9 and reports were gathered with MultiQC 1.14. [23]. The RNA-seq data showed neither adapter bias nor low quality at the 3′ end. No processing of the fragments was necessary. Phred scores per cycle for each library are reported in Supplementary File S1. In total, 893 million RNA-seq fragments from the two libraries (warm and cold) were used to assemble a de novo transcriptome in Trinity 2.6.6 [24] with default parameters for paired-end reads. Transcripts with 96% identity were clustered using the CD-HIT 4.6 software [25] to reduce redundant sequences. The open reading frames (ORFs) and putative proteins were predicted from the assembled transcriptomes using TransDecoder 5.5.0 [26], in combination with Swiss-Prot and Pfam searches (options: –retain_pfam_hits, –retain_blastp_hits). The functional annotation was performed with Trinotate [27] to gather the results of Blastn 2.6.0, Blastp 2.6.0 [28], Hmmer 3.1b2 [29], SignalP 4.1 [30], and TMHMM 2.0 [31] for each assembled transcript. This Transcriptome Shotgun Assembly project has been deposited at DDBJ/ENA/GenBank under the accession GKOC00000000. The version described in this paper is the first version, GKOC01000000.

2.3.4. Evaluation of Assembly

Although pooling many individuals is useful for obtaining sufficient genomic material from small organisms, a disadvantage is the increased likelihood of contamination, which compromises the quality of the assembly. Therefore, to confirm that the genomic material was representative of a single species, that is, C. confusus, the assembled transcriptome in this study was aligned against the 28s ribosomal gene of L. chilensis (NCBI ID: JX896325.2). Hit sequences were retrieved and aligned against the NCBI nucleotide collection (nt) database using Blastn, with default parameters. The results revealed the presence of sequences of more than one caligid species. It is possible that a few specimens other than C. confusus were inadvertently collected, as the observation of individuals was not conducted with absolute rigor to avoid subjecting them to stress related to handling. Because sequences cannot be filtered computationally between closely related species, we continued the posterior analysis using this transcriptome.

To eliminate sequences of fish and possible contamination by bacteria, virus, fungi, and another species, strict filtering was performed by aligning the putative proteins against the UniRef90 database (accessed on 19 June 2022) [32] using Diamond 2.0.15 [33] with the option “--sensitive” (e-value < 1 × 10⁻³). Similar to the findings of Caña-Bozada et al., the hit sequences that were the best matches with non-Protostomia species were considered contaminants and eliminated from the transcriptome and putative proteins [34]. Completeness of the transcriptome was assessed in BUSCO 5.3.2 [35] using the dataset of orthologous genes conserved across arthropods, which contains 1013 BUSCO groups.

2.3.5. Differentially Expressed Transcripts and Enrichment Analysis

For the differential expression analysis, first, the RNA-seq reads were mapped to the assembled transcriptome using the package Kallisto 0.42.4 [36], to obtain a matrix of transcript counts and abundances. The matrix of transcript abundances was used as an input in the Bioconductor package edgeR 3.40.2 [37] to identify differentially expressed transcripts (DETs) by comparing the transcripts obtained from the cold group against those from the warm. In this analysis, the transcripts with expression levels below 4 counts per million reads (cpm) were removed and transcript counts were normalized using the “calcNormFactors” function. Transcripts with the values log (fold change) ≥ 2.0 and p < 0.01 were considered DETs.

Functional enrichment analysis was performed using two approaches—Gene Ontology (GO) enrichment analysis and gene set enrichment analysis (GSEA). Enrichment analysis of the GO terms was performed with TopGop 2.52.0 [38], using Fisher’s exact testing and the weight algorithm. Term enrichment was considered at p < 0.01. To avoid redundancy, only those GO terms with fewer than 400 annotated proteins were considered.

GSEA was performed using Fast Gene Set Enrichment Analysis (FGSEA 1.24.0) [39]. This analysis was performed with the following parameters: nperm = 10,000; minSize = 14; maxSize = 400. Gene sets with p < 0.01 were considered as differentially enriched. GSEA results were collapsed with the collapsePathways function, implemented in the fgsea package, to reduce redundancy.

3. Results

3.1. Changes in Infection Levels

The prevalence of C. confusus in Ca. caninus fish was significantly higher in cold (42.3%, CI 30.4–63.7) than in warm (25.8%, CI 17.4–35.9; p < 0.05) tropical waters, whereas its mean intensity did not exhibit significant differences (p > 0.05), with values of 3.33 (CI 2.3–4.79) and 2.36 (CI 1.55–3.18) parasites per fish, respectively.

3.2. Transcriptomic Responses

De Novo Transcriptome Assembly of C. confusus

The de novo transcriptome of C. confusus combining the two MGI-seq libraries (warm and cold) generated 230,861 transcripts and 182,298 genes. The N10 and N50 results were 3841 and 1000, respectively, with an average length of 1012.3 bp (Table 1). We obtained an average of 35% mapping for all the libraries. A total of 101,640 ORF and protein sequences were predicted from the transcriptome. After contaminant sequences were filtered out, 55,530 ORF and protein sequences were retained. The ORF and protein sequences in fasta format are available in Supplementary File S2.

Of the 294,679 assembled transcripts, 8.3% were contaminant sequences. The final transcriptome with clean sequences resulted in 230,861 transcripts. A total of 1211 GO terms were assigned to cellular component (CC, 10.30%), 7903 to biological process (BP, 67.23%), and 2641 to molecular function (MF, 22.41%) (Figure 1). BUSCO analysis of the transcripts revealed a completeness of 94.7% (15% single-copy BUSCOs and 79.7% duplicated BUSCOs), 1.7% of fragmented genes, and 3.7% missing genes in the transcriptome of C. confusus (Figure 2).

3.3. DET and Enrichment Analysis

DET analysis was conducted on a subset of 20,640 genes that exhibited cpm values greater than 4. This analysis showed that, under contrasting temperatures, 1828 transcripts were differentially expressed (LFC ≥ 2.0 and p < 0.01); of these, 426 were over-expressed and 1402 were down-expressed in colder sea temperature (Supplementary File S3).

The GO enrichment analysis revealed that 66 processes were up-regulated in cold conditions, of which 35 are associated with BP, 17 with MF, and 14 with CC. Conversely, 111 processes were down-regulated, of which 62 relate to BP, 27 to MF, and 22 to CC. The most prominent processes are depicted in Figure 3.

FGSEA identified 59 enriched GO terms (p < 0.01), 41 of which correspond to BP, 15 to MF, and 3 to CC (Figure 4). Particularly in processes related to reproduction and development, we observed transcripts encoding vitellogenin (vit-1, vit-2, vit-4, and vit-6) in the GO terms’ nutrient reservoir activity, and metalloproteinases (mmp-2 and mmp-9) in the GO terms’ endodermal cell differentiation and positive regulation of vascular-associated smooth muscle cell proliferation. Other processes involving metalloproteinases were also found, including cellular response to UV-A, response to amyloid-beta, and negative regulation of cation channel activity (Figure 3 and Supplementary File S4).

4. Discussion

4.1. Changes in Infection Levels and Transcriptomic Response

The higher prevalence of C. confusus on Ca. caninus in cold seas is similar to that observed in another tropical sea louse—C. omissus on Scomberomorus sierra [11]. These findings contrast with the population dynamics of sea lice at higher latitudes. Lõhmus and Björklund reported that the occurrence of parasites in general in temperate regions has strong seasonality, with the occurrence and transmission being higher during warm periods [12]. This seasonality has been observed in the sea lice L. salmonis, C. epinepheli, and C. rogercresseyi from temperate latitudes [40,41,42]. This is because sea lice follow the temperature-dependent seasonal patterns common in ectotherms, with their development being faster in summer than in winter [7].

Conversely, in the tropics, it seems that the increase in temperature compromises the reproduction and infectivity of sea lice, as suggested by the low prevalence of C. confusus observed in the warm season. In general, animals in shallow tropical marine ecosystems subsist near their upper temperature limits and are highly vulnerable to increases in temperature [43]. This issue has been shown for the free-living copepod Pseudodiaptomus annandalei, which results in prolonged development, reduced size at maturity, smaller clutch sizes, lower hatching success, and reduced naupliar production under the extreme warming conditions (34 °C) that occur frequently during summer periods [44]. In sea lice larvae, body size and food reserves are reduced at higher temperatures, so that the infective copepodid stage remains viable for a shorter period of time, whereas at lower temperatures viability is higher because sea lice produce more and larger eggs and have more time to find a host [7,9,45].

It is possible that in the (tropical) cold season C. confusus produces larger eggs with greater energy reserves which leads to higher prevalence. This hypothesis is supported by our transcriptomic analysis, which revealed a positive regulation of vitellogenin and metalloproteinase transcripts in C. confusus sampled in the cold season relative to specimens from the warm season. Vitellogenins are a group of proteins that play a crucial role in the reproductive biology of many animals. In fish, warm temperatures can reduce vitellogenin levels, which might be reflected in poor survival of ova and infertility [46]. In particular, in egg-laying species such as copepods, vitellogenin serves as a reserve of nutrients for the embryos as they develop inside the egg and on which they survive after hatching until they feed autonomously [47,48]. Semmouri et al. observed that the expression of vitellogenins was positively correlated with the density of the free-living copepods [49]. Something similar could occur in C. confusus, in such a way that vitellogenins provide energy for embryonic development and food larval stages, leading to higher infection rates and, consequently, higher prevalence.

In C. confusus, the over-expression of the nutrient reservoir activity process was activated by vitellogenin genes vit-1, vit-2, vit-4, and vit-6. vit-1 and vit-2 have been previously identified and characterized in the sea lice L. salmonis and C. rogercresseyi [13,16,50], which is consistent with our results. Conversely, vit-4 and vit-6 have not been previously reported in copepods. In nematodes, increased expression of vit-4 and vit-6 is associated with aging [51,52]. It would be interesting to carry out more studies to validate and characterize other possible vitellogenin precursor genes in parasitic copepods, such as vit-4 and vit-6. The identification and understanding of vitellogenins might provide valuable insights into the reproductive biology and ecological dynamics of sea lice.

Other important transcripts over-expressed in C. confusus in colder tropical seas were matrix metalloproteinases (MMPs) mmp-2 and mmp-9, which, according to the GO annotation, are involved in different biological processes and molecular functions such as endodermal cell differentiation and regulation of vascular cells associated with smooth muscle cell proliferation. MMPs are involved in cell–cell adhesion and cell and basement membrane remodeling, and in the degradation of the extracellular matrix [53]. In invertebrates, MMPs are important for tissue remodeling [54]. To the best of our knowledge, to date there have been no reports of the role of metalloproteinases in copepods. In Drosophila melanogaster, only two MMPs genes have been identified—mmp-1 and mmp-2—which participate in various processes during embryogenesis by histolysis, tissue remodeling, nervous system development, and neuronal dendritic remodeling [54]. In other arthropods, such as the crab Eriocheir sinensis, molting processes and innate immune responses are regulated by MMPs [55], whereas in the silkworm, Bombyx mori, MMPs are functionally required for fat body cell dissociation and ovary development in female pupae [56]. In the parasitic nematode Haemonchus contortus, silencing of mmp-12 significantly reduced egg count, larval hatchability, and adult worm count and size [57]. Based on this evidence, and with the over-expression of MMPs observed herein, we suppose that the cold tropical season physiologically favors C. confusus, and that low sea temperatures favor embryo development, which might ultimately lead to higher prevalence.

Climate change is causing an increase in ocean temperatures, and some studies propose that the effects of climate change might be more severe in coastal areas of tropical regions [19,43]. It is projected that by the year 2100, the sea temperatures will increase by 2 °C [58]. This increase can have negative consequences on sea lice populations, primarily because a 2 °C rise might exceed the tolerance limits of certain species [8,11,59]. Based on the present study, we could expect that climate change will reduce some tropical sea lice populations during extreme warming events. This could alter ecosystem processes in general because parasites are vital for regulating host abundance [60]; parasite extinctions might have unforeseen costs that impact the health and abundance of a large number of free-living species [61]. However, it would be interesting to investigate whether parasites can develop adaptive molecular mechanisms to counteract the detrimental effects of sea warming, as has been proposed for free-living species [62]. Thus, it is important to perform more detailed studies on the molecular biology of parasites.

4.2. De Novo Transcriptome Quality

This study presents a first draft of the de novo transcriptome for the tropical sea louse C. confusus. Only a few de novo transcriptomes are available for parasitic copepods (Table 2). The transcriptome of C. confusus presented herein has good quality metrics, one of which is the guanine–cytosine (GC) content of 36.96%. In free-living and parasitic copepods, GC content is typically around 30–40% [63,64]. The completeness (94.7%) indicated by BUSCO suggests that the majority of the expected genes were successfully recovered, indicating a high-quality assembly and accurate sequencing. Furthermore, the BUSCO score for C. confusus is higher than that of other assemblies of parasitic copepods such as C. rogercresseyi, L. salmonis, and Lernaeocera branchialis (Figure 2).

5. Conclusions

The present study provides the first transcriptome for a tropical sea louse species (C. confusus). The analysis of this transcriptome indicated over-expression of vitellogenin and metalloproteinase transcripts in specimens obtained from colder topical waters (21 °C) compared with those from warmer waters (30 °C). These proteins favor growth, development, and survival in the specimens obtained. This might help to explain the higher prevalence of C. confusus on its host Ca. caninus during February, one of the coldest months of the year in the coastal Mazatlán region.

Author Contributions

Conceptualization, F.N.M.-S.; methodology, F.N.M.-S., C.A.P.-A., V.H.C.-B. and J.M.O.-C.; software, C.O.-V. and V.H.C.-B.; formal analysis, C.A.P.-A., V.H.C.-B. and J.M.O.-C.; investigation, F.N.M.-S. and C.A.P.-A.; resources, C.O.-V. and F.N.M.-S.; data curation, C.O.-V., C.A.P.-A. and V.H.C.-B.; writing—original draft preparation, C.A.P.-A. and F.N.M.-S.; writing—review and editing, E.J.F.-Á.; C.O.-V. and V.H.C.-B.; visualization, V.H.C.-B.; supervision, F.N.M.-S. and E.J.F.-Á.; funding acquisition, E.J.F.-Á. and F.N.M.-S. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

This study did not consider experiments with live animals. All fish were obtained from commercial catches. The fish used to determine the prevalence and mean intensity of C. confusus were already deceased, and the fish from which C. confusus specimens were obtained for transcriptomics studies were handled according to the American Veterinary Medical Association (AVMA) Guidelines for the Euthanasia of Animals (AVMA 2020). None of the species involved are subject to conservation measures.

Data Availability Statement

Raw sequence data is available on the NCBI website under project number PRJNA1004606. Accession numbers of Short Read Archive (SRA) accession numbers SRR25611633 and SRR25611632 (https://www.ncbi.nlm.nih.gov/sra/?term=PRJNA1004606). Transcriptome Shotgun Assembly project has been deposited at DDBJ/ENA/GenBank under the accession GKOC00000000. The version described in this paper is the first version, GKOC01000000 (https://www.ncbi.nlm.nih.gov/nuccore/GKOC00000000).

Acknowledgments

The temperature and dissolved oxygen data were provided by the Global and Climate Change Laboratory of the Mazatlán Academic Unit of the Institute of Marine Sciences and Limnology at the National Autonomous University of Mexico. In particular, we wish to acknowledge the work carried out at the Coastal Observatory for Global Change of Mazatlan by Joan Albert Sánchez Cabeza (lead), Ana Carolina Ruiz Fernández, Marcela Guillermina Fregoso López, Benjamín Yáñez Chávez, Beatriz Yáñez Rivera, Sergio Rendón Rodríguez, Roberto Cruz García, Libia Hascibe Pérez Bernal, and León Felipe Álvarez Sánchez (UNINMAR). We also thank Jose Sálgado Barragán (ICML, UNAM) for his help in the laboratory work.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 1. Gene Ontology processes of annotated transcripts of C. confusus de novo assembly. BP, biological process; CC, cellular component; MF, molecular function. # on the X axis; number.

Figure 2. BUSCO completeness analysis of C. confusus transcriptome and other copepods of the order Siphonostomatoida. T, transcriptome; EST, expressed sequence tag. Completeness was assessed using the dataset of orthologous genes conserved across arthropods, which contains 1013 BUSCO groups.

Figure 3. Results of the GO terms enrichment analysis. The graphs report some GO terms assigned to biological process (BP), molecular function (MF), or cellular component (CC) and up-regulated (A) or down-regulated (B) in C. confusus taken from cold waters.

Figure 4. FGSEA results. The graph reports the 24 most enriched GO terms in C. confusus taken from cold waters.

Table 1

Summary statistics for de novo assembly and annotation of the Caligus confusus transcriptome.

Descriptors	De Novo Assembly
Reads (millions)	893
Trinity transcripts (n)	230,861
Trinity predicted genes (n)	182,298
Assembled bases	233,706,455
Average contig length (bp)	1012.3
Median contig length (bp)	732
N10 (bp)	3841
N50 (bp)	1000
GC content (%)	36.96
Overall mapping (%)	35
Annotation (total)
Predicted ORF	55,530
Swiss-Prot top BLASTX hit	44,902
Gene Ontology in Blastx	44,280
Kegg	39,021
Swiss-Prot top BLASTP hit	31,478
Pfam	30,596
Gene Ontology in blastp	31,050
Gene Ontology in Pfam	18,751
SignalP	5504
Eggnog	479

bp, base pairs.

Table 2

Metrics of the de novo assembly transcriptomes of selected parasitic copepods.

Copepod Species	Caligus confusus	Tracheliastes polycolpus	Lepeophtheirus salmonis	Lernaea cyprinacea	Caligus rogercresseyi
Reference	This study	[65]	[66]	[67]	[68]
Sequencing platform	MGIseq 2000	Illumina HiSeq 2000	Illumina HiSeq 2000	Illumina HiSeq 2000	Illumina MiSeq
Assembly software	Trinity 2.6.6	Trinity version 2014-07-17	Trinityversion 2012-06-08	Velvet_1.2.10 and Oases_0.2.08	CLC GenomicsWorkbench
Raw reads (millions)	433.7	420	389.92	31.67	151.78
Assembled bases	233,706,455	20,355,465	*	50,792	83,444
Transcripts filtered	230,861	17,157	33,933	57,349	83,444
Average contig length (bp)	1012.3	1186	1469	1071.2	819
Median contig length	732	*	979	23,727	*
N50 (bp)	1000	1599	2100	1750	*
Overall mapping (%)	35	65.8	*	*	*
GC content (%)	36.96	38.55	36.9	*	*
Transcripts with GO terms	32,152	12,098	28,547	2439	*
Transcripts with KEGG terms	39,021	10,924	*	*	*
Pfam	30,596	11,751	*	*	*

* Data not provided.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/fishes8100475/s1, Supplementary File S1: Phred scores per cycle for each library; Supplementary File S2: The ORF and protein sequences in fasta format; Supplementary File S3: Results obtained from differential expression analysis; Supplementary File S4: Results of FGSEA analysis.

References

1. Walter, T.C.; Boxshall, G. World of Copepods Database. Available online: https://www.marinespecies.org/copepoda (accessed on 11 May 2023).

2. Boxshall, G.A.; Halsey, S.H. An Introduction to Copepod Diversity; Natural History Museum: London, UK, 2004; Volume 15, 996.

3. Abolofia, J.; Asche, F.; Wilen, J.E. The Cost of Lice: Quantifying the Impacts of Parasitic Sea Lice on Farmed Salmon. Mar. Resour. Econ.; 2017; 32, pp. 329-349. [DOI: https://dx.doi.org/10.1086/691981]

4. Costello, M.J. The Global Economic Cost of Sea Lice to the Salmonid Farming Industry. J. Fish Dis.; 2009; 32, pp. 115-118. [DOI: https://dx.doi.org/10.1111/j.1365-2761.2008.01011.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19245636]

5. Hemmingsen, W.; MacKenzie, K.; Sagerup, K.; Remen, M.; Bloch-Hansen, K.; Dagbjartarson Imsland, A.K. Caligus elongatus and other Sea Lice of the Genus Caligus as Parasites of Farmed Salmonids: A Review. Aquaculture; 2020; 522, 735160. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2020.735160]

6. Schram, T. Seasonal Occurrence of Lepeophtheirus salmonis and Caligus elongatus (Copepoda: Caligidae) on Sea Trout (Salmo Trutta), off Southern Norway. ICES J. Mar. Sci.; 1998; 55, pp. 163-175. [DOI: https://dx.doi.org/10.1006/jmsc.1997.0357]

7. Samsing, F.; Oppedal, F.; Dalvin, S.; Johnsen, I.; Vågseth, T.; Dempster, T. Salmon Lice (Lepeophtheirus salmonis) Development Times, Body Size, and Reproductive Outputs Follow Universal Models of Temperature Dependence. Can. J. Fish. Aquat. Sci.; 2016; 73, pp. 1841-1851. [DOI: https://dx.doi.org/10.1139/cjfas-2016-0050]

8. Hamre, L.; Bui, S.; Oppedal, F.; Skern-Mauritzen, R.; Dalvin, S. Development of the Salmon Louse Lepeophtheirus salmonis Parasitic Stages in Temperatures Ranging from 3 to 24 °C. Aquac. Environ. Interact.; 2019; 11, pp. 429-443. [DOI: https://dx.doi.org/10.3354/aei00320]

9. Montory, J.A.; Cumillaf, J.P.; Gebauer, P.; Urbina, M.; Cubillos, V.M.; Navarro, J.M.; Marín, S.L.; Cruces, E. Early Development and Metabolic Rate of the Sea Louse Caligus rogercresseyi under Different Scenarios of Temperature and PCO2. Mar. Environ. Res.; 2020; 162, 105154. [DOI: https://dx.doi.org/10.1016/j.marenvres.2020.105154]

10. Morales-Serna, F.N.; Rubio-Godoy, M.; Gómez, S. Seasonality of Parasitic Copepods on Bullseye Puffer, Sphoeroides annulatus (Pisces: Tetraodontidae), From the Northwestern Coast of Mexico. J. Parasitol.; 2011; 97, pp. 565-573. [DOI: https://dx.doi.org/10.1645/GE-2638.1]

11. Morales-Serna, F.N.; Olivas-Padilla, L.; Marín-Enriquez, E.; Osuna-Cabanillas, J.M.; Aguirre-Villaseñor, H.; Hernández-Covarrubias, V. Copepod and Monogenean Infection of Fish under Sea Surface Temperature Anomalies. J. Sea Res.; 2021; 175, 102098. [DOI: https://dx.doi.org/10.1016/j.seares.2021.102098]

12. Lõhmus, M.; Björklund, M. Climate Change: What Will It Do to Fish-Parasite Interactions?. Biol. J. Linn. Soc.; 2015; 116, pp. 397-411. [DOI: https://dx.doi.org/10.1111/bij.12584]

13. Eichner, C.; Frost, P.; Dysvik, B.; Jonassen, I.; Kristiansen, B.; Nilsen, F. Salmon Louse (Lepeophtheirus salmonis) Transcriptomes during Post Molting Maturation and Egg Production, Revealed Using EST-Sequencing and Microarray Analysis. BMC Genom.; 2008; 9, 126. [DOI: https://dx.doi.org/10.1186/1471-2164-9-126]

14. Chávez-Mardones, J.; Gallardo-Escárate, C. Next-Generation Transcriptome Profiling of the Salmon Louse Caligus Exposed to Deltamethrin (AlphaMax^TM): Discovery of Relevant Genes and Sex-Related Differences. Mar. Biotechnol.; 2015; 17, pp. 793-810. [DOI: https://dx.doi.org/10.1007/s10126-015-9661-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26307019]

15. Poley, J.D.; Igboeli, O.O.; Fast, M.D. Towards a Consensus: Multiple Experiments Provide Evidence for Constitutive Expression Differences among Sexes and Populations of Sea Lice (Lepeophtheirus salmonis) Related to Emamectin Benzoate Resistance. Aquaculture; 2015; 448, pp. 445-450. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2015.06.026]

16. Gallardo-Escárate, C.; Valenzuela-Muñoz, V.; Boltaña, S.; Nuñez-Acuña, G.; Valenzuela-Miranda, D.; Gonçalves, A.T.; Détrée, C.; Tarifeño-Saldivia, E.; Farlora, R.; Roberts, S. The Caligus Rogercresseyi MiRNome: Discovery and Transcriptome Profiling during the Sea Lice Ontogeny. Agri Gene; 2017; 4, pp. 8-22. [DOI: https://dx.doi.org/10.1016/j.aggene.2017.03.002]

17. Morales-Serna, F.N.; Medina-Guerrero, R.M.; Fajer-Avila, E.J. Sea Lice (Copepoda: Caligidae) Parasitic on Fishes Reported from the Neotropical Region. Neotrop. Biodivers.; 2016; 2, pp. 141-150. [DOI: https://dx.doi.org/10.1080/23766808.2016.1236313]

18. Soler-Jiménez, L.C.; Morales-Serna, F.N.; Aguirre-Macedo, M.L.; Mclaughlin, J.P.; Jaramillo, A.G.; Shaw, J.C.; James, A.K.; Hechinger, R.F.; Kuris, A.M.; Lafferty, K.D. et al. Parasitic Copepods (Crustacea, Hexanauplia) on Fishes from the Lagoon Flats of Palmyra Atoll, Central Pacific. ZooKeys; 2019; 833, pp. 85-106. [DOI: https://dx.doi.org/10.3897/zookeys.833.30835]

19. Sanchez-Cabeza, J.-A.; Herrera-Becerril, C.A.; Carballo, J.L.; Yáñez, B.; Álvarez-Sánchez, L.F.; Cardoso-Mohedano, J.-G.; Ruiz-Fernández, A.C. Rapid Surface Water Warming and Impact of the Recent (2013–2016) Temperature Anomaly in Shallow Coastal Waters at the Eastern Entrance of the Gulf of California. Prog. Oceanogr.; 2022; 202, 102746. [DOI: https://dx.doi.org/10.1016/j.pocean.2022.102746]

20. Ho, J.S.; Lin, C.L. Sea Lice of Taiwan; The Sueichan Press: Keelung, Taiwan, 2004.

21. Bush, A.O.; Lafferty, K.D.; Lotz, J.M.; Shostak, A.W. Parasitology Meets Ecology on its Own Terms: Margolis et al. Revisited. J. Parasitol.; 1997; 83, pp. 575-583. [DOI: https://dx.doi.org/10.2307/3284227]

22. Reiczigel, J.; Marozzi, M.; Fábián, I.; Rózsa, L. Biostatistics for Parasitologists–A Primer to Quantitative Parasitology. Trends Parasitol.; 2019; 35, pp. 277-281. [DOI: https://dx.doi.org/10.1016/j.pt.2019.01.003]

23. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. 2010; Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 15 June 2023).

24. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q. et al. Full-Length Transcriptome Assembly from RNA-Seq Data without a Reference Genome. Nat. Biotechnol.; 2011; 29, pp. 644-652. [DOI: https://dx.doi.org/10.1038/nbt.1883]

25. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for Clustering the next-Generation Sequencing Data. Bioinformatics; 2012; 28, pp. 3150-3152. [DOI: https://dx.doi.org/10.1093/bioinformatics/bts565]

26. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M. et al. De Novo Transcript Sequence Reconstruction from RNA-Seq Using the Trinity Platform for Reference Generation and Analysis. Nat. Protoc.; 2013; 8, pp. 1494-1512. [DOI: https://dx.doi.org/10.1038/nprot.2013.084]

27. Bryant, D.M.; Johnson, K.; DiTommaso, T.; Tickle, T.; Couger, M.B.; Payzin-Dogru, D.; Lee, T.J.; Leigh, N.D.; Kuo, T.-H.; Davis, F.G. et al. A Tissue-Mapped Axolotl De Novo Transcriptome Enables Identification of Limb Regeneration Factors. Cell Rep.; 2017; 18, pp. 762-776. [DOI: https://dx.doi.org/10.1016/j.celrep.2016.12.063]

28. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic Local Alignment Search Tool. J. Mol. Biol.; 1990; 215, pp. 403-410. [DOI: https://dx.doi.org/10.1016/S0022-2836(05)80360-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/2231712]

29. Wheeler, T.J.; Eddy, S.R. Nhmmer: DNA Homology Search with Profile HMMs. Bioinformatics; 2013; 29, pp. 2487-2489. [DOI: https://dx.doi.org/10.1093/bioinformatics/btt403]

30. Petersen, T.N.; Brunak, S.; Heijne, G.; Nielsen, H. SignalP 4.0: Discriminating Signal Peptides from Transmembrane Regions. Nat. Methods; 2011; 8, pp. 785-786. [DOI: https://dx.doi.org/10.1038/nmeth.1701] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21959131]

31. Krogh, A.; Larsson, B.; Von Heijne, G.; Sonnhammer, E.L.L. Predicting Transmembrane Protein Topology with a Hidden Markov Model: Application to Complete Genomes. J. Mol. Biol.; 2001; 305, pp. 567-580. [DOI: https://dx.doi.org/10.1006/jmbi.2000.4315] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11152613]

32. The UniProt Consortium. UniProt: The universal protein knowledgebase in 2021. Nucleic Acids Res.; 2021; 49, pp. D480-D489. [DOI: https://dx.doi.org/10.1093/nar/gkaa1100]

33. Buchfink, B.; Xie, C.; Huson, D.H. Fast and Sensitive Protein Alignment Using DIAMOND. Nat. Methods; 2015; 12, pp. 59-60. [DOI: https://dx.doi.org/10.1038/nmeth.3176]

34. Caña-Bozada, V.; Morales-Serna, F.N.; Fajer-Ávila, E.J.; Llera-Herrera, R. De Novo Transcriptome Assembly and Identification of G-Protein-Coupled-Receptors (GPCRs) in Two Species of Monogenean Parasites of Fish. Parasite; 2022; 29, 51. [DOI: https://dx.doi.org/10.1051/parasite/2022052]

35. Manni, M.; Berkeley, M.R.; Seppey, M.; Zdobnov, E.M. BUSCO: Assessing Genomic Data Quality and Beyond. Curr. Protoc.; 2021; 1, e323. [DOI: https://dx.doi.org/10.1002/cpz1.323] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34936221]

36. Bray, N.L.; Pimentel, H.; Melsted, P.; Pachter, L. Near-Optimal Probabilistic RNA-Seq Quantification. Nat. Biotechnol.; 2016; 34, pp. 525-527. [DOI: https://dx.doi.org/10.1038/nbt.3519] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27043002]

37. Baldoni, P.L.; Chen, Y.; Hediyeh-zadeh, S.; Liao, Y.; Dong, X.; Ritchie, M.E.; Shi, W.; Smyth, G.K. Dividing out Quantification Uncertainty Allows Efficient Assessment of Differential Transcript Expression. bioRxiv; 2023; [DOI: https://dx.doi.org/10.1101/2023.04.02.535231]

38. Alexa, A.; Rahnenführer, J. Gene Set Enrichment Analysis with TopGO. Bioconductor Improv.; 2009; 27, pp. 1-26.

39. Korotkevich, G.; Sukhov, V.; Budin, N.; Shpak, B.; Artyomov, M.N.; Sergushichev, A. Fast Gene Set Enrichment Analysis. bioXiv; 2016; [DOI: https://dx.doi.org/10.1101/060012]

40. Rikardsen, A.H. Seasonal Occurrence of Sea Lice Lepeophtheirus salmonis on Sea Trout in Two North Norwegian Fjords. J. Fish Biol.; 2004; 65, pp. 711-722. [DOI: https://dx.doi.org/10.1111/j.0022-1112.2004.00478.x]

41. Khamees, N.R.; Adday, T.K. Occurrence of Sea Lice Caligus Epinepheli Yamaguti, 1936 (Copepoda: Siphonostomatoida) on Gills of Nemipterus japonicus (Bloch, 1775) from Northwest of the Arabian Gulf. Basrah J. Agric. Sci.; 2013; 26, pp. 1-14. [DOI: https://dx.doi.org/10.33762/bagrs.2013.111443]

42. Montes, R.M.; Quiñones, R.A.; Gallardo-Escárate, C. Disentangling the Effect of Sea Temperature and Stocking Density on Sea Louse Abundance Using Wavelets in a Highly Infested Salmon Farming Area. Aquaculture; 2022; 546, 737246. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2021.737246]

43. Nguyen, K.D.T.; Morley, S.A.; Lai, C.-H.; Clark, M.S.; Tan, K.S.; Bates, A.E.; Peck, L.S. Upper Temperature Limits of Tropical Marine Ectotherms: Global Warming Implications. PLoS ONE; 2011; 6, e29340. [DOI: https://dx.doi.org/10.1371/journal.pone.0029340]

44. Doan, N.X.; Vu, M.T.T.; Pham, H.Q.; Wisz, M.S.; Nielsen, T.G.; Dinh, K.V. Extreme Temperature Impairs Growth and Productivity in a Common Tropical Marine Copepod. Sci. Rep.; 2019; 9, 4550. [DOI: https://dx.doi.org/10.1038/s41598-019-40996-7]

45. Bravo, S.; Erranz, F.; Lagos, C. A Comparison of Sea Lice, Caligus rogercresseyi, Fecundity in Four Areas in Southern Chile. J. Fish Dis.; 2009; 32, pp. 107-113. [DOI: https://dx.doi.org/10.1111/j.1365-2761.2008.01012.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19245635]

46. King, H.R.; Pankhurst, N.W.; Watts, M.; Pankhurst, P.M. Effect of Elevated Summer Temperatures on Gonadal Steroid Production, Vitellogenesis and Egg Quality in Female Atlantic Salmon: Vitellogenesis in atlantic salmon. J. Fish Biol.; 2003; 63, pp. 153-167. [DOI: https://dx.doi.org/10.1046/j.1095-8649.2003.00137.x]

47. Lee, C.S.; O’Bryen, P.J.; Marcus, N.H. Copepods in Aquaculture; John Wiley & Sons: Hoboken, NJ, USA, 2008.

48. Tarrant, A.M.; Nilsson, B.; Hansen, B.W. Molecular Physiology of Copepods-from Biomarkers to Transcriptomes and Back Again. Comp. Biochem. Physiol. Part D Genom. Proteom.; 2019; 30, pp. 230-247. [DOI: https://dx.doi.org/10.1016/j.cbd.2019.03.005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30921652]

49. Semmouri, I.; Asselman, J.; Van Nieuwerburgh, F.; Deforce, D.; Janssen, C.R.; De Schamphelaere, K.A.C. The Transcriptome of the Marine Calanoid Copepod Temora longicornis under Heat Stress and Recovery. Mar. Environ. Res.; 2019; 143, pp. 10-23. [DOI: https://dx.doi.org/10.1016/j.marenvres.2018.10.017]

50. Dalvin, S.; Frost, P.; Loeffen, P.; Skern-Mauritzen, R.; Baban, J.; Rønnestad, I.; Nilsen, F. Characterisation of Two Vitellogenins in the Salmon Louse Lepeophtheirus salmonis: Molecular, Functional and Evolutional Analysis. Dis. Aquat. Organ.; 2011; 94, pp. 211-224. [DOI: https://dx.doi.org/10.3354/dao02331]

51. Zhai, C.; Zhang, N.; Li, X.; Chen, X.; Sun, F.; Dong, M. Fusion and Expansion of Vitellogenin Vesicles during Caenorhabditis elegans Intestinal Senescence. Aging Cell; 2022; 21, e13719. [DOI: https://dx.doi.org/10.1111/acel.13719]

52. Tang, J.; Ma, Y.-C.; Chen, Y.-L.; Yang, R.-Q.; Liu, H.-C.; Wang, X.; Ni, B.; Zou, C.-G.; Zhang, K.-Q. Vitellogenin Accumulation Leads to Reproductive Senescence by Impairing Lysosomal Function. Sci. China Life Sci.; 2023; 66, pp. 439-452. [DOI: https://dx.doi.org/10.1007/s11427-022-2242-8]

53. Luchian, I.; Goriuc, A.; Sandu, D.; Covasa, M. The Role of Matrix Metalloproteinases (MMP-8, MMP-9, MMP-13) in Periodontal and Peri-Implant Pathological Processes. Int. J. Mol. Sci.; 2022; 23, 1806. [DOI: https://dx.doi.org/10.3390/ijms23031806]

54. Pagemccaw, A. Remodeling the Model Organism: Matrix Metalloproteinase Functions in Invertebrates. Semin. Cell Dev. Biol.; 2008; 19, pp. 14-23. [DOI: https://dx.doi.org/10.1016/j.semcdb.2007.06.004]

55. Li, R.; Meng, Q.; Huang, J.; Wang, S.; Sun, J. MMP-14 Regulates Innate Immune Responses to Eriocheir Sinensis via Tissue Degradation. Fish Shellfish Immunol.; 2020; 99, pp. 301-309. [DOI: https://dx.doi.org/10.1016/j.fsi.2020.02.021]

56. Jia, Q.; Chen, X.; Wu, L.; Ruan, Z.; Li, K.; Li, S. Matrix Metalloproteinases Promote Fat Body Cell Dissociation and Ovary Development in Bombyx Mori. J. Insect Physiol.; 2018; 111, pp. 8-15. [DOI: https://dx.doi.org/10.1016/j.jinsphys.2018.10.002]

57. Naqvi, M.A.-H.; Li, H.; Gao, W.; Naqvi, S.Z.; Jamil, T.; Aimulajiang, K.; Xu, L.; Song, X.; Li, X.; Yan, R. Haemonchus contortus: SiRNA Mediated Knockdown of Matrix Metalloproteinase 12A (MMP-12) Results in Reduction of Infectivity. Parasit. Vectors; 2020; 13, 151. [DOI: https://dx.doi.org/10.1186/s13071-020-04025-1]

58. Widlansky, M.J.; Long, X.; Schloesser, F. Increase in Sea Level Variability with Ocean Warming Associated with the Nonlinear Thermal Expansion of Seawater. Commun. Earth Environ.; 2020; 1, 9. [DOI: https://dx.doi.org/10.1038/s43247-020-0008-8]

59. Dalvin, S.; Are Hamre, L.; Skern-Mauritzen, R.; Vågseth, T.; Stien, L.; Oppedal, F.; Bui, S. The Effect of Temperature on Ability of Lepeophtheirus Salmonis to Infect and Persist on Atlantic Salmon. J. Fish Dis.; 2020; 43, pp. 1519-1529. [DOI: https://dx.doi.org/10.1111/jfd.13253] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32882750]

60. Sato, T.; Egusa, T.; Fukushima, K.; Oda, T.; Ohte, N.; Tokuchi, N.; Watanabe, K.; Kanaiwa, M.; Murakami, I.; Lafferty, K.D. Nematomorph Parasites Indirectly Alter the Food Web and Ecosystem Function of Streams through Behavioural Manipulation of Their Cricket Hosts. Ecol. Lett.; 2012; 15, pp. 786-793. [DOI: https://dx.doi.org/10.1111/j.1461-0248.2012.01798.x]

61. Dobson, A.; Lafferty, K.D.; Kuris, A.M.; Hechinger, R.F.; Jetz, W. Homage to Linnaeus: How Many Parasites? How Many Hosts?. Proc. Natl. Acad. Sci. USA; 2008; 105, (Suppl. S1), pp. 11482-11489. [DOI: https://dx.doi.org/10.1073/pnas.0803232105]

62. Bowman, L.L.; Kondrateva, E.S.; Timofeyev, M.A.; Yampolsky, L.Y. Temperature Gradient Affects Differentiation of Gene Expression and SNP Allele Frequencies in the Dominant Lake Baikal Zooplankton Species. Mol. Ecol.; 2018; 27, pp. 2544-2559. [DOI: https://dx.doi.org/10.1111/mec.14704] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29691934]

63. Guo, Y.; Ye, F.; Sheng, Q.; Clark, T.; Samuels, D.C. Three-Stage Quality Control Strategies for DNA Re-Sequencing Data. Brief. Bioinform.; 2014; 15, pp. 879-889. [DOI: https://dx.doi.org/10.1093/bib/bbt069] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24067931]

64. Jørgensen, T.S.; Petersen, B.; Petersen, H.C.B.; Browne, P.D.; Prost, S.; Stillman, J.H.; Hansen, L.H.; Hansen, B.W. The Genome and MRNA Transcriptome of the Cosmopolitan Calanoid Copepod Acartia tonsa Dana Improve the Understanding of Copepod Genome Size Evolution. Genome Biol. Evol.; 2019; 11, pp. 1440-1450. [DOI: https://dx.doi.org/10.1093/gbe/evz067]

65. Mathieu-Bégné, E.; Loot, G.; Blanchet, S.; Toulza, E.; Genthon, C.; Rey, O. De Novo Transcriptome Assembly for Tracheliastes polycolpus, an Invasive Ectoparasite of Freshwater Fish in Western Europe. Mar. Genomics; 2019; 46, pp. 58-61. [DOI: https://dx.doi.org/10.1016/j.margen.2018.12.001]

66. Carmona-Antoñanzas, G.; Carmichael, S.N.; Heumann, J.; Taggart, J.B.; Gharbi, K.; Bron, J.E.; Bekaert, M.; Sturm, A. A Survey of the ATP-Binding Cassette (ABC) Gene Superfamily in the Salmon Louse (Lepeophtheirus salmonis). PLoS ONE; 2015; 10, e0137394. [DOI: https://dx.doi.org/10.1371/journal.pone.0137394] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26418738]

67. Pallavi, B.; Shankar, K.M.; Abhiman, P.B.; Iqlas, A. Identification of Putative Genes Involved in Parasitism in the Anchor Worm, Lernaea cyprinacea by de Novo Transcriptome Analysis. Exp. Parasitol.; 2015; 153, pp. 191-197. [DOI: https://dx.doi.org/10.1016/j.exppara.2015.03.014] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25816972]

68. Gallardo-Escárate, C.; Valenzuela-Muñoz, V.; Nuñez-Acuña, G. RNA-Seq Analysis Using De Novo Transcriptome Assembly as a Reference for the Salmon Louse Caligus rogercresseyi. PLoS ONE; 2014; 9, e92239. [DOI: https://dx.doi.org/10.1371/journal.pone.0092239] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24691066]