- Open Access
Comparative transcriptional profile of the fish parasite Cryptocaryon irritans
Parasites & Vectorsvolume 9, Article number: 630 (2016)
Cryptocaryon irritans is an obligate ectoparasitic ciliate pathogen of marine fishes. It can infect most marine teleosts and cause heavy economic losses in aquaculture. There is currently no effective method of controlling this disease, and little information is available regarding the genes involved in its development and virulence. We aimed to investigate the distinct features of the three major life-cycle stages of C. irritans in terms of gene transcription level, and identify candidate vaccines/drug targets. We established a reference transcriptome of C. irritans by RNA-seq.
Three cDNA libraries using total poly(A)+ mRNA isolated from trophonts, tomonts, and theronts was constructed and sequenced, respectively. Clean reads from the three stages were de novo assembled to generated unigene. Annotation of unigenes and transcriptomic comparison of three stages was performed.
Totals of 73.15, 62.23, and 109.57 million clean reads were generated from trophont, tomont, and theront libraries, respectively. After de novo assembly, 49,104 unigenes were obtained, including 9,253 unigenes with significant similarities to proteins from other ciliates. Transcriptomic comparisons revealed that 2,470 genes were differentially expressed among the three stages, including 2,011, 1,404, and 1,797 genes that were significantly differentially expressed in tomont/theront, tomont/trophont, and theront/trophont pairwise comparisons, respectively. Based on the results of hierarchical clustering, all differentially expressed genes (DEGs) were located in five major clusters. DEGs in clusters 1 and 2 were more highly expressed in tomonts than in other stages, DEGs in cluster 3 were dominant in the tomont and trophont stages, whereas clusters 4 and 5 included genes upregulated in the theront stage. In addition, Immobilization antigens (I-antigens) and proteases have long been considered major targets for vaccine development and potential drug targets in parasites, respectively. In the present study, nine putative I-antigens transcripts and 161 protease transcripts were found in the transcriptome of C. irritans.
It was concluded that DEGs enriched in tomonts were involved in cell division, to increase the number of theronts and ensure parasite continuity. DEGs enriched in theronts were associated with response to stimuli, whereas genes enriched in trophonts were related to nutrient accumulation and cell growth. In addition, the I-antigen and protease transcripts in our transcriptome could contribute to the development of vaccines or targeted drugs. Together, the results of the present study provide novel insights into the physiological processes of a marine parasitic ciliate.
The protozoan Cryptocaryon irritans is one of the most common parasitic ciliates. It can infect many species of marine fish in tropical and subtropical regions, causing cryptocaryonosis [1–4]. Although the parasite does not usually cause severe infections in wild fish, it has become a major problem in aquaria and aquaculture industries. In South China, direct economic losses due to cryptocaryoniasis have amounted to over 16 million US dollars just in the Guangdong Province each year (data not published). Cryptocaryon irritans has four main life-cycle stages (Fig. 1): parasitic trophont, off-host protomont, reproductive tomont and infective theront . The trophont lives within the host epithelium, and the typical clinical signs of cryptocaryonosis, including pinhead-sized white nodules covering the surface of the skin, gills and eyes, can be observed by the naked eye during this stage. The mature trophont leaves the host and becomes a protomont before encysting and transforming into a tomont, which then attaches to the substrate, develops, and divides into numerous daughter tomites by asymmetric binary fission. These tomites subsequently leave the cyst as theronts, which actively seek new fish hosts and burrow into the epithelium to start the next life-cycle. An entire life-cycle takes one week in Trachinotus ovatus hosts used to maintain C. irritans at 27 °C in artificial conditions .
Despite major losses caused by cryptocaryonosis, there is currently no effective means of controlling this disease . However, previous studies indicated that C. irritans can regulate the expression of host immune-related genes [7–12], and immunization with C. irritans can confer protection [13–15], implicating immune prophylaxis as a potential alternative method of control. However, it is impractical to culture C. irritans continuously in vitro in large-scale industrial operations, and using fish as hosts to passage the parasite is both time-consuming and costly. Surface proteins called immobilization antigens (I-antigens) were recently identified in this parasite, and a modified DNA vaccine based on these antigens was demonstrated to protect fish against infection by C. irritans . However, no other potential antigen proteins have yet been identified.
Information on the transcriptome of a parasite can help to discover genes and understand the molecular processes involved in parasite development, reproduction and host interactions, as well as facilitating the search for potential vaccine candidates and drug targets [17–22]. Comparative analysis of Ichthyophthirius multifiliis trophonts and theronts showed that the most abundant transcripts were highly stage-specific and coincided with metabolic activities . In 2010, Lokanathan et al.  identified 2,659 expressed sequence tags (ESTs) from a cDNA library of C. irritans tomonts, of which just 1,611 matched proteins with known functions. Although these results offered insight into the genomics of C. irritans, they only identified gene sequences, with little information about the transcript levels of each gene, or the genes involved in different stages of the life-cycle. The present study therefore aimed to obtain a more complete transcript picture of C. irritans by sequencing mRNA from all three life-cycle stages (theront, trophont, and tomont) using the Illumina RNA-seq method. A total of 9,253 high-quality unigenes were identified that significantly matched other ciliate proteins. We analyzed the expression profiles of these unigenes and described their differential expression patterns in the three stages. The results of this study improve our understanding of C. irritans biology and will help to further the development of effective methods for controlling cryptocaryonosis.
Cryptocaryon irritans for sequencing were originally isolated from an infected T. ovatus obtained from a local farm in Daya Bay, Guangdong Province, China, and maintained by serial passage using T. ovatus as hosts, as described previously . To collect parasites at different life stages, mature trophonts were gently scraped from infected T. ovatus gills (3 days post-infection) and carefully washed to remove contaminating tissue debris. The trophonts were then incubated in sterilized seawater at 28 °C for 12 h to develop into tomonts, or for 3 days to develop into theronts. After centrifugation at 800× g for 5 min, trophonts and tomonts were collected. Theronts were cooled in an ice-bath for 30 min, and then harvested after centrifugation at 4,000× g for 5 min. All C. irritans were immediately cryopreserved in liquid nitrogen until RNA isolation.
RNA was extracted from trophonts, tomonts, and theronts using an RNeasy Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. The quantity and integrity of the RNA were determined using a NanoDropND-1000 spectrophotometer (Thermo, Waltham, USA) and Agilent 2100 Bioanalyzer (Agilent Technologies, Richardson, USA), respectively, and RNAs with a RNA Integrity Number (RIN) > 8, 28S/18S > 0.7, and A260/280 values of about 2.0 were used to prepare the libraries.
Library construction and sequencing
Sequencing libraries were constructed according to the TruSeq® RNA Sample Preparation v2 manuals (Illumina, San Diego, USA). Briefly, poly(A)+ mRNA was purified from the total RNA with oligo(dT) attached magnetic beads, and cleaved into small fragments. First-strand cDNA was synthesized using the fragmented RNA as template and random hexamers as primers, and second-strand cDNA was obtained using Second Strand Master Mix (Thermo, Waltham, USA). After repairing the overhang ends into blunt ends with End Repair Mix, a single “A” nucleotide was added to the 3′ end of the double-stranded cDNA, followed by ligation to an adapter containing a single “T”. Polymerase chain reaction (PCR) was used to enrich the DNA fragments with adapter molecules on both ends and amplify the amount of DNA in the library. The concentrations and sizes of the libraries were measured with a Qubit®2.0 Fluorometer (Life Technologies, Camarillo, USA) and Agilent 2100 Bioanalyzer, respectively. The libraries were then sequenced and analyzed using an Illumina HiSeq™ 2500 instrument (Illumina, San Diego, USA) with paired end 2 × 100 nucleotide multiplex, according to the manufacturer’s instructions  (Fig. 1).
De novo assembly and gene annotation
Adaptors, low-quality reads (> 50% bases with quality (Q) value ≤ 10), ambiguous nucleotide reads (> 10% ‘N’ rate, where ‘N’ represents ambiguous bases in reads), and sequence reads < 20 bp were removed. All clean reads from the three stages served as a pool reads and were processed with CLC Genomics Workbench software [24–26] using the scaffolding contig program (word-size = 45, minimum contig length ≥ 300) and CAP3 EST software to generate de novo assembled unigenes. All unigenes were first compared with the T. ovatus transcriptome database (http://www.ncbi.nlm.nih.gov/sra/SRX534464) to remove host-cell contamination. A bacterial database (total 222,262,566 proteins) was then downloaded from a non-redundant (NR) protein database to remove prokaryote contamination, with a cut-off E-value for contaminating contigs of < 1e-5 and identity ≥ 70% in the above analysis. Finally, the clean unigenes were annotated using a high-throughput BLASTx program against the NR protein database, and protein databases of Tetrahymena thermophila (http://www.ciliate.org/system/downloads/T_thermophila_June2014_proteins.fasta), Tetrahymena borealis (http://www.ciliate.org/system/downloads/T_borealis_oct2012_proteins.fasta), Tetrahymena elliotti (http://www.ciliate.org/system/downloads/T_elliotti_oct2012_proteins.fasta), Tetrahymena malaccensis (http://www.ciliate.org/system/downloads/T_malaccensis_oct2012_proteins.fasta), Paramecium tetraurelia (http://paramecium.cgm.cnrs-gif.fr/download/fasta/Ptetraurelia_peptides_cur.fasta), and I. multifiliis (http://ich.ciliate.org/system/downloads/img1_0407.aa.fsa) (E-value < 1e-5, identity ≥ 30% and coverage ≥ 50%). Potential C-terminal glycosylphosphatidylinositol (GPI)-modification sites were predicted using big-II predictor (http://mendel.imp.ac.at/sat/gpi/gpi_server.html) . Gene ontology (GO) annotation was performed using Blast2GO and GO enrichment analysis using top GO package. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were obtained from KEGG databases.
DEGs and cluster analysis
The relative expression levels of unigenes were calculated as reads per kilobases per million reads (RPKM) . To compare the pronounced active transcripts among the three stages of C. irritans, the edegR package was used to define significantly upregulated or downregulated genes with a threshold value of fold change ≥ 4 and false discovery rate (FDR) < 0.001 . Based on the expression pattern, cluster analysis of differentially expressed genes (DEGs) in the three stages was performed using Cluster package. A heat map was generated using R, and a GO annotation plot was generated using WEGO (http://wego.genomics.org.cn/cgi-bin/wego/index.pl) .
Experimental validation of transcription levels
To confirm the RNA-seq results, gene expression levels were determined using the LightCycler480 real-time PCR system (Roche, Mannheim, Germany) using SYBR Green Realtime PCR Master Mix (Thermo, Waltham, USA), according to the manufacturer’s instructions, as described previously . Total RNA was extracted from each parasite stage using an RNeasy Mini Kit, as described above, and stored at -80 °C. cDNA was synthesized from total RNA for each parasite stage using a ReverTra Ace-α-Kit (Toyobo, Osaka, Japan), according to the manufacturer’s protocol. Eight gene-specific primers were designed based on the transcriptome library (Additional file 1: Table S1). Elongation factor-1β (EF-1β) primers (EF-1β RTF/EF-1β RTR) were introduced as the reference gene. The cycling protocol was as follows: 94 °C for 2 min, and 94 °C for 15 s, 58 °C for 15 s, and 72 °C for 20 s for 40 cycles. Melting-curve analysis was used to detect the specificity of the PCR products, which was verified by sequencing. All samples were analyzed in triplicate. The expression of each target gene was normalized to the housekeeping gene EF-1β, calculated according to the 2-ΔΔCt method . Tomont samples were used as the calibration control. Thus, ΔΔCt = [(Cttarget – CtEF-1β)]trophont or theront stage – [(Cttarget – CtEF-1β)] tomont stage. All data were analyzed using SPSS (version 16.0) software. Correlations between quantitative PCR (qPCR) and RNA-seq results were determined using Pearson’s r.
Sequencing and de novo assembly
To examine the transcriptional profile of C. irritans, we constructed three cDNA libraries using total poly(A)+ mRNA isolated from trophonts, tomonts, and theronts, respectively. Totals of 79.35, 66.42, and 123.62 million raw reads, respectively, were produced by Illumina Hiseq 2500 sequencing (Table 1). All raw reads were submitted to the Sequence Read Archive database at NCBI (SUB1416064, SUB1416075, and SUB1416142). After removal of adaptors, low-quality reads, ambiguous reads, and reads < 20 bp, 73.15, 62.23, and 109.57 million clean reads were generated from the trophont, tomont, and theront libraries, respectively, with a clean-read ratio of 90.92%. A total of 49,104 unigenes were then obtained after de novo assembly, with a total length of 50,174,026 bp, an average length of 1,022 bp, and N50 of 1,338 bp.
To obtain more C. irritans-derived genes, the unigenes were first compared with a T. ovatus transcriptome database to obtain 31,893 non-host contigs, followed by a bacterial database to remove the prokaryote contamination. This left 9,253 contigs (Coding sequence ≥ 100 amino acids) (Table 1), of which 8,569 contigs matched entries in the NR protein database, and all 9,253 contigs (Additional file 1: Tables S2 and S3) matched to the six ciliate databases, corresponding to the number of I. multifiliis annotations [33, 34]. A Venn diagram summarizing C. irritans unigenes compared with peptide sequences from Tetrahymena (four species in total), I. multifiliis, and P. tetraurelia is shown in Fig. 2. A total of 7,041 (76.1%) unigenes were shared by all ciliate databases.
GO function annotation
For functional predictions and categories, all 9,253 unigenes were annotated with 43,436 GO terms and assigned into three functional GO terms, including cellular component (14 sub-categories), molecular function (13 sub-categories), and biological process categories (22 sub-categories) (Fig. 3 and Additional file 1: Table S4). Cell part (66.3%) was the largest sub-category in the cellular component category, while the main sub-categories in molecular function were catalytic activity (47.6%) and binding (43.3%), and the major biological processes were cellular process (60.7%) and metabolic process (47.3%).
KEGG enrichment analysis
Mapping the reference pathways in the KEGG database identified a total of 230 KEGG pathways (Additional file 1: Table S5). The pathways with most genes were the ribosome, spliceosome, carbon metabolism, and purine metabolism pathways, which was in accord with the previous transcriptome analysis of C. irritans tomonts . Comparisons with the T. ovatus transcriptome  and other ciliate genomes indicated that the overall metabolism of ciliates was largely identical, with only minor differences.
DEGs and cluster analysis
In the present study, we used a threshold of FDR < 0.001 and absolute fold change ≥ 4 to define significantly upregulated or downregulated genes in pairwise comparisons. A total of 2,470 DEGs were identified across the three stages, including 2,011 genes that were significantly differentially expressed between tomonts and theronts, comprising 1,103 upregulated and 908 downregulated genes. A total of 1,404 DEGs, including 631 upregulated and 773 downregulated genes in the tomont/trophont pairwise comparison, and 1,797 DEGs, including 805 upregulated and 992 downregulated genes in the theront/trophont pairwise comparison. Based on hierarchical clustering results, we generated a heat map to illustrate the differential gene expression patterns in the three stages (Fig. 4). DEGs formed five major clusters (Fig. 5 and Additional file 1: Table S6). DEGs in clusters 1 and 2 were more highly expressed in tomonts than in the other stages, while DEGs in cluster 3 had lower expression levels in theronts but relatively higher expression levels in tomonts and trophonts, while clusters 4 and 5 showed increased expression in the theront stage.
GO enrichment analysis (Fig. 6) revealed that DEGs in clusters 1 and 2, including the genes encoding histone H2A variant 3, protein kinase 2, Cyclic Adenosine monophosphate (cAMP)-dependent protein kinase regulatory subunit, cat-eye syndrome critical region protein 2, translationally-controlled tumor protein homolog, calpain-type cysteine protease DEK1, and calpain-type cysteine protease ADL1 were enriched in biological processes of reproduction and regulation of cell growth. DEGs such as those encoding phosphatidylinositol-4-phosphate 5-kinase 1, phosphatidylinositol-4-phosphate 5-kinase 4, phosphatidylinositol-4-phosphate 5-kinase 6, phosphatidylinositol-4-phosphate 5-kinase 9, and SNF1-related protein kinase catalytic subunit alpha KIN10, which were involved in lipid metabolism, were also enriched in those clusters.
DEGs in cluster 3, such as those involved in glycometabolism (glycogen phosphorylase 1, glycogen debranching enzyme, phosphoglucomutase-2, glucose-6-phosphate isomerase, cytosolic 1), lipid metabolism (acetyl-CoA acetyltransferase, hormone-sensitive lipase, phosphoglycerate kinase), amino acid metabolism (aspartate aminotransferase, alanine aminotransferase, glutamine-fructose-6-phosphate aminotransferase) were largely enriched in the metabolism of nutrient substances. In addition, genes such those encoding 40S ribosomal protein S3-3, aconitate hydratase 2, aldehyde dehydrogenase family 7 member B4, endoplasmin homolog, proteasome subunit beta type-4, shaggy-related protein kinase epsilon, and sodium/hydrogen exchanger 7, which were involved in the biological process of response to abiotic stimulus, salt stress and osmotic stress, were also increased in cluster 3.
DEGs in clusters 4 and 5, including those for calpain-D, E3 ubiquitin-protein ligase CHIP, SNW/SKI-interacting protein, cation/calcium exchanger 3, mitogen-activated protein kinase 2 (MAPK2), MAPK 3, MAPK 4 and serine/threonine-protein kinase AtPK1/AtPK6 were enriched in response to osmotic stress, abiotic stimuli, or salt stress. Genes for aspartic proteinase-like protein 1, aspartic proteinase-like protein 2, serine carboxypeptidase-like 47, serine carboxypeptidase-like 49, serine/threonine-protein kinase AFC1, ubiquitin carboxyl-terminal hydrolase 4 and methionine aminopeptidase 2B, which were involved in proteolysis and protein modification process, were also increased in cluster 4.
Identification of I-antigens and proteases
I-antigens and proteases are recognized as major targets for vaccine development and potential drug targets in parasites, respectively. A total of nine putative I-antigens were found in our transcriptome (Additional file 1: Table S3). Two I-antigen transcripts were highly expressed in all three stages of C. irritans (RPKM > 200), and shared only 57% similarity with each other. In addition, 122 proteases were found in the transcriptome, most of which were calpain family cysteine proteases (82). Six cysteine proteases were highly expressed (RPKM > 100) in trophonts.
Eight genes were randomly selected from the transcriptome library to verify the transcription levels in each stage by qPCR. After normalizing to the EF-1β gene, the relative expression levels were presented as the fold-change relative to the tomont stage. The correlation coefficient (r) between RNA-Seq and qPCR ranged from 0.8188 to 0.9906, depending on the specific gene and reference gene (Fig. 7).
Understanding the comparative transcriptional profile of C. irritans is an essential goal for parasite physiology and the prevention of cryptocaryonosis. In this study, we performed comparative gene transcription analysis among the three life-cycle stages of C. irritans. To obtain a more complete reference transcriptome, all raw reads from the three stages were first pooled and de novo assembled. A more specific C. irritans-derived transcriptome was obtained by removing contamination from host fish and prokaryotes by comparison with a T. ovatus transcriptome database and a bacterial database. A total of 9,253 unigenes were obtained, 7,041 of which were shared by other ciliate databases. A comprehensive repertoire of gene annotation, GO function annotation, KEGG analysis, DEG, and cluster analysis was conducted.
The macronuclear genome of I. multifiliis was published in 2011  and is available as a reference database for I. multifiliis mapping. However, there is currently no available C. irritans reference database for mapping, and we therefore used the de novo method to assemble the C. irritans transcriptome. Moreover, C. irritans must be cultured in host fish, and we recently demonstrated that nonspecific cytotoxic cells receptor protein+ (NCCRP+) cells can be found inside food vacuoles in trophonts . In addition, C. irritans cannot be cultured axenically, and contamination by host RNA and/or RNA from endosymbiotic bacterium cannot be completely eliminated when preparing samples. Although oligo (dT) magnetic beads were used to enrich the eukaryotic mRNA, almost half the assembled unigenes still matched to prokaryotic genes (22,640). A similar study in I. multifiliis showed that 17.2% of ESTs in the I. multifiliis transcriptome strongly matched the genome sequence of an endosymbiotic bacterium [34, 37]. A previous report of C. irritans transcriptome data also showed that only 57% of unique transcripts matched to other ciliate species . We cannot fully explain the reason for the high proportion of prokaryote transcripts in our transcriptome, but focused on those unigenes that significantly matched to other ciliate species.
Tomonts represent a crucial stage in the life-cycle of C. irritans, ensuring continuity of the parasite. In this stage, asynchronous theronts are released from tomonts post-encystment, and the maximum yield of theronts was 122-fold in in vitro propagation under experimental conditions . DEGs in clusters 1 and 2 were more highly expressed in tomonts. Consistent with a previous report in C. irritans , GO enrichment analysis showed that DEGs in clusters 1 and 2 were enriched in the biological processes of reproduction and regulation of cell growth, which are the major biological activities in cells during the tomont stage. Some genes involved in lipid metabolism were also enriched in clusters 1 and 2. Given that tomonts cannot obtain energy substrates from the external environment, the role of lipids as the major energy substrate in tomonts needs further investigation.
Trophonts are the only parasitic stage during the life-cycle of C. irritans, and need to acquire enough energy for the forthcoming reproductive tomont stage. DEGs in cluster 3 were highly expressed in the trophont and tomont stages, and were largely enriched in the utilization of nutrient substances. For protozoan parasites, the parasitic stage is an important state for acquiring nutrition from their host, and genes involved in metabolic processes thus dominate cellular activity. A previous transcriptional analysis of the model protozoan T. thermophila demonstrated that 60 of 148 selected representative genes were upregulated in the growth stage and were involved in metabolic processes . Approximately 30% of sequences in the parasite Trypanosoma vivax showed significantly high expression levels in the bloodstream stage, and were confirmed to be genes related to metabolic processes . Transcriptomic analysis of I. multifiliis revealed that transcripts related to metabolic enzymes, including succinyl-CoA ligase, aldo/keto reductase, and glutamine synthetase, were significantly upregulated in the trophont stage . In the present study, those metabolic enzymes were also highly expressed in the trophont stage. Highly expressed unigenes (RPKM > 100, 352 unigenes) encoding metabolic enzymes may relate to the importance of energy metabolism and material exchange in the cellular activity of parasitic trophonts. The oral route is a convenient way to apply medicines, and given that trophonts are the only parasitic stage during the life-cycle of C. irritans, they represent the ideal stage at which to administer oral drugs [40, 41]. Inhibitors of the proteins involved in trophont metabolic pathways may thus provide a means of controlling this parasite infection. In addition, DEGs involved in the biological processes of response to abiotic stimuli, salt stress, and osmotic stress were also increased in cluster 3, indicating that those genes may relate to the adjustment of permeation pressure during the transition from host body to saltwater, when the trophont leaves the host fish and develops into a tomont.
Among the three life-cycle stages of C. irritans, theronts are responsible for host invasion. DEGs in clusters 4 and 5 were more highly expressed in the theront stage. DEGs in cluster 4 were enriched in response to osmotic stress, abiotic stimuli, and salt stress. The active transcription of these genes in cluster 4 may reflect a series of stress reactions in response to exposure to saltwater, and may be involved in acclimating theronts to the change in permeation pressure during the transition from tomont to saltwater. Notably, we cooled the theronts in an ice-bath for 30 min for sample preparation, which may have affected transcription. However, genes related to response to temperature stimuli did not appear to be enriched in theronts. Protein metabolism was also active in cluster 4, and DEGs involved in proteolysis and protein modification were more highly transcribed in this cluster, suggesting that efficient protein metabolism is crucial for theronts exposed to saltwater. In addition, transcripts such as those for patatin family phospholipases and mitochondrial carrier proteins, which are reportedly involved in the balance of energy usage/storage [42, 43], were highly transcribed in clusters 4 and 5. Given that theronts cannot survive for long without finding a host, it seems that theronts require a careful energy balance during host finding.
Along with the DEGs described above, some transcripts highly expressed in all three stages of C. irritans merit additional consideration. I-antigens are well-studied surface proteins in some ciliates. Although their functions remain unknown, I-antigens in I. multifiliis have been shown to elicit protective immunity in fish and have long been considered major targets for vaccine development [44, 45]. Cryptocaryon irritans has previously been shown to elicit effective host systematic and skin humoral immune responses [46–48], and some putative I-antigens have been identified in C. irritans [49–51]. Vaccination with DNA or recombinant I-antigen vaccine improved post-challenge survival by 46% . In this study, we found nine putative I-antigen transcripts in our C. irritans transcriptome, which were expressed at various levels in the three stages (Additional file 1: Table S3). Two I-antigen transcripts were highly expressed in all three stages of C. irritans (RPKM > 200), and shared only 57% similarity with each other, though a potential C-terminal GPI-modification site was predicted in both transcripts. Further studies are needed to determine if these proteins share similar functions. These newly identified I-antigens provide more potential targets for vaccine development against C. irritans.
Proteases play crucial roles in parasite infection and development. Because of the feasibility of designing specific inhibitors, proteases have long been recognized as potential drug targets in parasites [33, 52–54]. For example, the synthetic peptide GlcA-Val-Leu-Gly-Lys-NHC2H5 effectively inhibited Plasmodium falciparum schizont cysteine protease Pf 68 in vitro . α2-Macroglobulin (α2M), a non-specific protease inhibitor of endogenous and exogenous proteases, was found to be involved in the immune response to I. multifiliis in infected carp . We previously established a transcriptome database of skin, gill, spleen and head kidney for the grouper Epinephelus coioides at different time points after C. irritans infection [54, 31]. We found three α2M isoforms in groupers, all of which were significantly upregulated in grouper gill and spleen post-infection (data unpublished). However, the role of α2M in the defense against C. irritans infection remains to be demonstrated. Nevertheless, these results suggest that anti-protease drugs could be exploited to control this parasite. We identified 122 protease transcripts in our transcriptome, most of which (82 transcripts) were calpain family cysteine proteases, thus significantly expanding the range of potential therapeutic targets for protease inhibitors. Calpain family cysteine proteases appear to play important roles in C. irritans, suggesting that they should be considered as the primary protease inhibitor targets. In addition, six proteases were highly expressed (RPKM > 100) in trophonts, and only shared about 30% amino acid similarity with teleosts. This suggests that specific inhibitors against C. irritans proteases could be designed based on those unigenes.
We established a reference transcriptome of C. irritans and characterized the genes and their functional categories at each stage of the parasite’s life-cycle. A total of 9,253 unigenes were generated by de novo assembly followed by removal of contamination and multiple annotations. A total of 2,470 DEGs were identified across the three life-cycle stages. DEGs enriched in tomonts were associated with cell division, while DEGs enriched in theronts were involved in response to stimuli, and DEGs enriched in trophonts were related to the accumulation of nutrients and cell growth. In addition, the identification of proteases and putative I-antigen transcripts in the transcriptome could contribute to vaccine development or drug targeting. Further studies are planned focusing on functional verification of the important DEGs in each life-cycle stage of C. irritans, with a view to developing effective vaccines or drugs to control this disease.
Differentially expressed gene
Elongation factor-1 β
False discovery rate
Mitogen-activated protein kinase
Reads Per Kilo bases per Million reads
Luo XC, Xie MQ, Zhu XQ, Li AX. Some characteristics of host-parasite relationship for Cryptocaryon irritans isolated from South China. Parasitol Res. 2008;102(6):1269–75.
Burgess PJ, Matthews RA. Fish host range of seven isolates of Cryptocaryon irritans (Ciliophora). J Fish Biol. 1995;46:727–9.
Colorni A, Burgess P. Cryptocaryon irritans Brown 1951, the cause of ‘white spot disease’ in marine fish: an update. Aquar Sci Conserv. 1997;1:217–38.
Wang FH, Xie MQ, Li AX. A novel protein isolated from the serum of rabbitfish (Siganus oramin) is lethal to Cryptocaryon irritans. Fish Shellfish Immunol. 2010;29:32–41.
Dan XM, Li AX, Lin XT, Teng N, Zhu XQ. A standardized method to propagate Cryptocaryon irritans on a susceptible host pompano Trachinotus ovatus. Aquaculture. 2006;258:127–33.
Dan XM, Zhang TW, Li YW, Li AX. Immune responses and immune-related gene expression profile in orange-spotted grouper after immunization with Cryptocaryon irritans vaccine. Fish Shellfish Immunol. 2013;34(3):885–91.
Mo ZQ, Li YW, Zhou L, Li AX, Luo XC, Dan XM. Grouper (Epinephelus coioides) IL-34/MCSF2 and MCSFR1/MCSFR2 were involved in mononuclear phagocytes activation against Cryptocaryon irritans infection. Fish Shellfish Immunol. 2015;43(1):142–9.
Li YW, Mo XB, Zhou L, Li X, Dan XM, Luo XC, et al. Identification of IRAK-4 in grouper (Epinephelus coioides) that impairs MyD88-dependent NF-κB activation. Dev Comp Immunol. 2014;45(1):190–7.
Li YW, Li X, Xiao XX, Zhao F, Luo XC, Dan XM, et al. Molecular characterization and functional analysis of TRAF6 in orange-spotted grouper (Epinephelus coioides). Dev Comp Immunol. 2014;44(1):217–25.
Li YW, Luo XC, Dan XM, Qiao W, Huang XZ, Li AX. Molecular cloning of orange-spotted grouper (Epinephelus coioides) TLR21 and expression analysis post Cryptocaryon irritans infection. Fish Shellfish Immunol. 2012;32:476–81.
Li YW, Luo XC, Dan XM, Huang XZ, Qiao W, Zhong ZP, et al. Orange-spotted grouper (Epinephelus coioides) TLR2, MyD88 and IL-1β involved in anti-Cryptocaryon irritans response. Fish Shellfish Immunol. 2011;30(6):1230–40.
Li YW, Dan XM, Zhang TW, Luo XC, Li AX. Immune-related genes expression profile in orange-spotted grouper during exposure to Cryptocaryon irritans. Parasite Immunol. 2011;33(12):679–987.
Bai JS, Xie MQ, Zhu XQ, Dan XM, Li AX. Comparative studies on the immunogenicity of theronts, tomonts and trophonts of Cryptocaryon irritans in grouper. Parasitol Res. 2008;102(2):307–13.
Luo XC, Xie MQ, Zhu XQ, Li AX. Protective immunity in grouper (Epinephelus coioides) following exposure to or injection with Cryptocaryon irritans. Fish Shellfish Immunol. 2007;22(4):427–32.
Misumi I, Lewis TD, Takemura A, Leong JA. Elicited cross-protection and specific antibodies in Mozambique tilapia (Oreochromis mossambicus) against two different immobilization serotypes of Cryptocaryon irritans isolated in Hawaii. Fish Shellfish Immunol. 2011;30(4-5):1152–8.
Jose Priya TA, Lin YH, Wang YC, Yang CS, Chang PS, Song YL. Codon changed immobilization antigen (iAg), a potent DNA vaccine in fish against Cryptocaryon irritans infection. Vaccine. 2012;30(5):893–903.
Fritz HM, Buchholz KR, Chen X, Durbin-Johnson B, Rocke DM, Conrad PA, Boothroyd JC. Transcriptomic analysis of Toxoplasma development reveals many novel functions and structures specific to sporozoites and oocysts. PLoS One. 2012;7(2):e29998.
Heizer E, Zarlenga D, Rosa B, Gao X, Gasser R, De Graef J, Geldhof P, Mitreva M. Transcriptome analyses reveal protein and domain families that delineate stage-related development in the economically important parasitic nematodes, Ostertagia ostertagi and Cooperia oncophora. BMC Genomics. 2013;14(1):118.
Morin-Adeline V, Lomas R, O'Meally D, Stack C, Conesa A, Slapeta J. Comparative transcriptomics reveals striking similarities between the bovine and feline isolates of Tritrichomonas foetus: consequences for in silico drug-target identification. BMC Genomics. 2014;15(1):955.
Chaudhary K, Roos DS. Protozoan genomics for drug discovery. Nat Biotechnol. 2005;23(9):1089–91.
Cassidy-Hanley DM, Cordonnier-Pratt MM, Pratt LH, Devine C, Mozammal Hossain M, Dickerson HW, Clark TG. Transcriptional profiling of stage specific gene expression in the parasitic ciliate Ichthyophthirius multifiliis. Mol Biochem Parasitol. 2011;178(1-2):29–39.
Lokanathan Y, Mohd-Adnan A, Wan KL, Nathan S. Transcriptome analysis of the Cryptocaryon irritans tomont stage identifies potential genes for the detection and control of cryptocaryonosis. BMC Genomics. 2010;11:76.
Xiong J, Lu X, Zhou Z, Chang Y, Yuan D, Tian M, et al. Transcriptome Analysis of the model protozoan, Tetrahymena thermophila, using deep RNA sequencing. PLoS One. 2012;7(2):e30630.
Brautigam A, Mullick T, Schliesky S, Weber AP. Critical assessment of assembly strategies for non-model species mRNA-Seq data and application of next-generation sequencing to the comparison of C(3) and C(4) species. J Exp Bot. 2011;62(9):3093–102.
Garg R, Patel RK, Tyagi AK, Jain M. De novo assembly of chickpea transcriptome using short reads for gene discovery and marker identification. DNA Res. 2011;18(1):53–63.
Su CL, Chao YT, Alex Chang YC, Chen WC, Chen CY, et al. De novo assembly of expressed transcripts and global analysis of the Phalaenopsis aphrodite transcriptome. Plant Cell Physiol. 2011;52(9):1501–14.
Eisenhaber B, Bork P, Yuan Y, Loffler G, Eisenhaber F. Automated annotation of GPI anchor sites: case study C. elegans. Trends Biochem Sci. 2000;25(7):340–1.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Meth. 2008;5(7):621–8.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Chetvernin V, et al. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2007;35(Database issue):D5–12.
Mo ZQ, Yang M, Wang HQ, Xu Y, Huang MZ, Lao GF, et al. Grouper (Epinephelus coioides) BCR signaling pathway was involved in response against Cryptocaryon irritans infection. Fish Shellfish Immunol. 2016;57:198–205.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.
Coyne RS, Hannick L, Shanmugam D, Hostetler JB, Brami D, Joardar VS, et al. Comparative genomics of the pathogenic ciliate Ichthyophthirius multifiliis, its free-living relatives and a host species provide insights into adoption of a parasitic lifestyle and prospects for disease control. Genome Biol. 2011;12(10):R100.
Abernathy J, Xu DH, Peatman E, Kucuktas H, Klesius P, Liu Z. Gene expression profiling of a fish parasite Ichthyophthirius multifiliis: Insights into development and senescence-associated avirulence. Comp Biochem Physiol Part D Genomics Proteomics. 2011;6(4):382–92.
Sun HY, Noe J, Barber J, Coyne RS, Cassidy-Hanley D, Clark TG, et al. Endosymbiotic bacteria in the parasitic ciliate Ichthyophthirius multifiliis. Appl Environ Microbiol. 2009;75(23):7445–52.
Zhenzhen X, Ling X, Dengdong W, Chao F, Qiongyu L, Zihao L, et al. Transcriptome analysis of the Trachinotus ovatus: identification of reproduction, growth and immune-related genes and microsatellite markers. PLoS One. 2014;9(10):e109419.
Greif G, Ponce de Leon M, Lamolle G, Rodriguez M, Pineyro D, Tavares-Marques LM, et al. Transcriptome analysis of the bloodstream stage from the parasite Trypanosoma vivax. BMC Genomics. 2013;14:149.
Kawano F, Hirazawa N. Antiparasitic effect of in-feed inhibitors of folic acid synthesis and dihydrofolate reductase against ciliate Cryptocaryon irritans infection in the red sea bream Pagrus major and against ciliate Ichthyophthirius multifiliis infection in black pop-eyed goldfish Carassius auratus. Aquaculture. 2012;330:1–7.
Hirazawa N, Oshima S-i, Hara T, Mitsuboshi T, Hata K. Antiparasitic effect of medium-chain fatty acids against the ciliate Cryptocaryon irritans infestation in the red sea bream Pagrus major. Aquaculture. 2001;198(3-4):219–28.
Kienesberger PC, Oberer M, Lass A, Zechner R. Mammalian patatin domain containing proteins: a family with diverse lipolytic activities involved in multiple biological functions. J Lipid Res. 2009;50(Suppl):S63–68.
Kuan J, Saier Jr MH. The mitochondrial carrier family of transport proteins: structural, functional, and evolutionary relationships. Crit Rev Biochem Mol Biol. 1993;28(3):209–33.
Wang X, Dickerson HW. Surface immobilization antigen of the parasitic ciliate Ichthyophthirius multifiliis elicits protective immunity in channel catfish (Ictalurus punctatus). Clin Diagn Lab Immunol. 2002;9(1):176–81.
Lin TL, Dickerson HW. Purification and partial characterization of immobilization antigens from Ichthyophthirius multifiliis. J Protozool. 1992;39(4):457–63.
Burgess PJ, Matthews RA. Cryptocaryon irritans (Ciliophora): acquired protective immunity in the thick-lipped mullet, Chelon labrosus. Fish Shellfish Immunol. 1995;5(6):459–68.
Yoshinaga T, Nakazoe J-I. Acquired protection and production of immunobilization antibody against Cryptocaryon irritans (Ciliophora, Hymenostomatida) in mummichog (Fundulus heteroclitus). Fish Pathology. 1997;32(4):229–30.
Hatanaka A, Umeda N, Yamashita S, Hirazawa N. Identification and characterization of a putative agglutination/immobilization antigen on the surface of Cryptocaryon irritans. Parasitology. 2007;134(Pt 9):1163–74.
Hatanaka A, Umeda N, Hirazawa N. Molecular cloning of a putative agglutination/immobilization antigen located on the surface of a novel agglutination/immobilization serotype of Cryptocaryon irritans. Parasitology. 2008;135(9):1043–52.
Huang X, Sun Z, Guo G, Zheng C, Xu Y, Yuan L, Liu C. Cloning and characterization of a surface antigen CiSA-32.6 from Cryptocaryon irritans. Exp Parasitol. 2012;130(3):189–94.
Atkinson HJ, Babbitt PC, Sajid M. The global cysteine peptidase landscape in parasites. Trends Parasitol. 2009;25(12):573–81.
Wu Y, Wang X, Liu X, Wang Y. Data-mining approaches reveal hidden families of proteases in the genome of malaria parasite. Genome Res. 2003;13(4):601–16.
Rosenthal PJ. Proteases of malaria parasites: new targets for chemotherapy. Emerg Infect Dis. 1998;4(1):49–57.
Mayer R, Picard I, Lawton P, Grellier P, Barrault C, Monsigny M, Schrevel J. Peptide derivatives specific for a Plasmodium falciparum proteinase inhibit the human erythrocyte invasion by merozoites. J Med Chem. 1991;34(10):3029–35.
Onara DF, Forlenza M, Gonzalez SF, Rakus KL, Pilarczyk A, Irnazarow I, Wiegertjes GF. Differential transcription of multiple forms of alpha-2-macroglobulin in carp (Cyprinus carpio) infected with parasites. Dev Comp Immunol. 2008;32(4):339–47.
Mo ZQ, Chen RA, Li YW, Huang XZ, Li AX, Luo XC, et al. Characterization and expression analysis of two novel CCR6 chemokine receptors and their three potential ligands CCL20Ls of grouper (Epinephelus coioides) post Cryptocaryon irritans infection. Fish Shellfish Immunol. 2015;47(1):280–8.
We acknowledge support of the Special Funds for Marine Fishery Science and Technology Promotion Project of Guangdong Province (grant No. A201501B11) to Dr. Xue-Ming Dan, Finance Supporting Specific Projects for Fish Diseases Control from the Government of Guangdong Province (grant No. YCN11H) to Dr. Yan-Wei Li and National Natural Science Foundation of China (grant No. 31272681) to Dr. An-Xing Li.
Special Funds for Marine Fishery Science and Technology Promotion Project of Guangdong Province (grant No. A201501B11) to Dr. Xue-Ming Dan. Finance Supporting Specific Projects for Fish Diseases Control from the Government of Guangdong Province (grant No. YCN11H) to Dr. Yan-Wei Li. National Natural Science Foundation of China (grant No. 31272681) to Dr. An-Xing Li.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files. All raw reads were submitted to the Sequence Read Archive database at NCBI (SUB1416064, SUB1416075, and SUB1416142).
Conceived and designed the experiments: XMD, XCL, AXL, ZQM. Conducted the experiment: ZQM, JLW, LYN, MY, GFL. Conducted the analysis: ZQM, YWL, HQW. Wrote the paper: ZQM, YWL, YZH, YZM. All authors read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Primers used in this study. Table S2. Gene annotation using ciliate databases. Table S3. Gene annotation using NR database. Table S4. GO function annotation. Table S5. KEGG annotation. Table S6. Number of DEGs in clusters. (XLSX 2815 kb)