Skip to main content

Transcriptome analysis of dormant tomonts of the marine fish ectoparasitic ciliate Cryptocaryon irritans under low temperature

Abstract

Background

Cryptocaryon irritans, a species of obligatory ciliate ectoparasite, can infect various species of marine teleost fish. Cryptocaryon irritans that fall to the seabed or aquarium bottom in winter can form “dormant tomonts” and wake up when the temperature rises the next year. Abundant studies and analyses on the dormant tomonts were carried out at the transcriptome level, in order to investigate the molecular mechanism of C. irritans tomonts entering the dormant state under low-temperature conditions.

Methods

The paired-end sequencing strategy was used to better assemble the entire transcriptome de novo. All clean sequencing reads from each of the three libraries (Group A: untreated blank control; Group B: treated for 24 h at 12 °C; and Group C: developed for 24 h at 25 °C) were respectively mapped back to the transcriptome assembly using the bioinformatics software.

Results

In this study, 25,695,034, 21,944,467, and 28,722,875 paired-end clean reads were obtained respectively from the three cDNA libraries of the C. irritans tomont by Illumina paired-end sequencing technology. A total of 25,925 unique transcript fragments (unigenes) were assembled, with an average length of 839 bp. Differentially expressed genes (DEGs) were scrutinized; in Group B/A pairwise comparison, 343 genes presented differential expression, including 265 up-regulated genes and 78 down-regulated genes in Group B; in Group C/A pairwise comparison, there were 567 DEGs, including 548 up-regulated genes and 19 down-regulated genes in Group C; and in Group B/C pairwise comparison, 185 genes showed differential expression, including 145 up-regulated genes and 40 down-regulated genes in Group B.

Conclusions

This is the first transcriptomic analytical study of the C. irritans tomonts under low temperature. It can be concluded that most of the genes required for its cell survival under low temperature, or for cell entry into a deeper dormancy state were discovered, and that they might be considered as candidate genes to develop the diagnostic and control measures for cryptocaryoniasis.

Background

Cryptocaryoniasis is categorized as a “Class II animal epidemic” in China, which is caused by the ciliate Cryptocaryon irritans inhabiting the body surface of marine teleosts. In recent years, this parasitic disease has often occurred in China’s coastal provinces and cities, bringing great losses to the fish farmers and relevant departments. In order to explore the pathogenesis of and effective control measures for C. irritans, researchers have carried out extensive studies from variable perspectives [13]. The regularity of parasite occurrence suggests that outbreaks of this disease are related to water temperature [4]. The incidence of this disease is relatively higher from April to November each year when the water temperature is between 20–30 °C than the rest of the year; after November, the incidence gradually reduces until complete disappearance; and after the 3–4-month low-temperature season, the disease would strike back again [5]. Researchers have revealed the 4-stage life-cycle of C. irritans: trophonts, protomonts, tomonts, and theronts. Trophonts mainly parasitize on the surfaces of the hosts, and leave their hosts after maturation, falling to seabed or aquarium bottom to form protomonts. Tomonts are formed after a brief phase of protomonts; and after continuous asymmetric division, 200–300 theronts are produced; theronts break the tomont wall and enter the water, swim very fast without food intake, and intrude into suitable hosts that they encounter, where they will develop into trophonts again [6]. Reports have proved that the development of the parasite is slow and can be even stopped at a low temperature. Tomonts can stay alive for 4–5 months at 12 °C. When water temperature rises, tomonts can wake up and start dividing again, producing infectious theronts. So it is believed that C. irritans that fall to seabed or aquarium bottom in winter can form “dormant tomonts” and wake up when the temperature rises next year [7].

Tomonts are a state of cells formed by ciliates and other protozoa after immobilizing from an active state, when they shrink gradually and lose some structures, followed by the formation of the tomont wall with secreted substances, forming a spherical or nearly spherical shape. Ciliate tomonts are classified into two types, namely dormant tomonts and proliferative tomonts; the former is a dormant state formed to withstand adverse environment, while the latter is a specific metamorphosis period of the life-cycle in which parasites in tomonts split into more daughter cells. It is therefore obvious that the formation and regulation of the two are different. Dormant tomonts are formed after a sudden change in temperature or food shortage, and they resume their normal activities after excystment once the environment is suitable. Current studies mainly focus on the processes of tomonts formation, dedifferentiation, and redifferentiation of free living ciliates [8]. The morphology of “proliferative tomonts” of C. irritans has also been described in great detail [1]. As for “dormant tomonts” described in this paper, however, they are formed when the cells stop dividing but keep alive after the water temperature decreases. This has been extensively investigated in other parasites and ciliates, e.g. Colpoda maupasi [9], C. steinii [10], dinoflagellate Scrippsiella hangoei [11], Alexandrium catenella and A. tamarense [12], etc. However, no study on the formation and regulation of dormant tomonts of C. irritans has been carried out yet.

A transcriptome represents all RNA transcripts in one cell or tissue, and reflects genes expressed in specific tissues in different life-cycle stages, physiological states, and environments [13]. Transcriptome studies can holistically exhibit functions and structures of genes and reveal the molecular mechanism of biological process and pathogenesis [14], thus transcriptomics has been widely applied in fundamental research, clinical diagnosis, drug development, and potential vaccine candidate proteins screening, etc. In recent years, RNA-sequencing has become a widely used approach in the studies on the development of ciliates, parasites, e.g. Leishmania donovani [15], salmon louse Caligus rogercresseyi [16], and Tetrahymena thermophila [17], etc. For the studies on cryptocaryoniasis, Lokanathan et al. [18] generated and analyzed ESTs of C. irritans tomonts to identify genes that encode surface proteins, excretory/secretory proteins and repeat-containing proteins; and this is the only report so far. In the present study, tomonts were induced to enter the state of dormancy at 12 °C and the changes in transcriptome of dormant tomonts were compared with RNA-seq technology to explore the molecular mechanism of C. irritans tomonts entering dormant state in the low-temperature season.

Methods

Cryptocaryon irritans tomonts and collection

The C. irritans were derived from a naturally infected Larimichthys crocea, and L. crocea with an average body mass of 100 g were then used as animal models to establish the passage system [19]. The animal models were raised in a 1000 l aquarium (R × H: 60 × 60 cm), and were infected with a non-lethal concentration of theronts (≤ 10,000 theronts/fish) in 5 l of seawater per fish; 2 h after infection, fresh seawater was added. Four  days after infection, large numbers of tomonts were found to adhere to the bottom of aquarium. The fish were then transferred to another clean aquarium without tomonts and tomonts were collected by carefully discarding the debris and incubated in a 1 l beaker. Throughout the whole experiment, the water was oxygenated continuously and replaced to keep clean twice a day (09:00 and 15:00); the salinity, water temperature, light intensity, and photoperiod for aquaculture were 29–31 ‰, 26 ± 1 °C, 1000 lx, and 12 Light: 12 Dark, respectively.

Newly formed tomonts were collected within 10 h and divided into 3 groups: Group A, B and C. Group A was an untreated blank control group and placed in liquid nitrogen after drying; group B was treated for 24 h at 12 °C and placed in liquid nitrogen after drying; group C was set to develop normally for 24 h at room temperature (25 °C) without any extra-processing and placed in liquid nitrogen after drying. Previous studies have confirmed that at least a period of 36 h is required for newly formed tomonts to start dividing at about 25 °C, so the processing time in this study was set to 24 h to ensure that all cells should be about to divide but not yet. It was also confirmed microscopically that no cell division occurred in all three groups of tomonts described above before sampling.

Total RNA extraction and sample preparation for RNA-Seq

Total RNA was extracted from the tomonts using an RNeasy Plus Universal Mini Kit (QIAGEN, Maryland, USA), which is specific for genome DNA elimination, according to the manufacturer’s instructions. All RNA was processed by RNase free DNase I provided by the kit. The quantity, purity and integrity of RNA were measured on a 1.2 % (w/v) agarose gel and on a Nanodrop-1000 spectrophotometer (NanoDrop, USA). Samples with higher quality (absorbance ratios at 260 nm/280 nm > 1.8) were selected for high-throughput sequencing. The extracted total RNA was resuspended in RNase free water and stored at -80 °C until use [20]. After enrichment with oligo-dT-attached magnetic beads, the purified mRNA was fragmented by divalent cations under elevated temperature, and then considered as a template for first-strand cDNA synthesis by random primers and reverse transcriptase. The second-strand cDNA was synthesized with RNase H (Invitrogen, USA) and DNA polymerase I. A paired-end library was constructed from the cDNA synthesized by the Genomic Sample Prep Kit (Illumina). Multiplexed cDNA libraries were mixed in equal volumes with normalized 10 nM concentration (Agilent 2100) [21, 22]. Three normalized cDNA libraries were constructed with the RNA from the control group, group B and group C. The library was sequenced by Illumina Mi-Seq platform.

Assembly, comparative analysis and functional annotation of the transcriptome

To better assemble the entire transcriptome de novo, a paired-end (PE) sequencing strategy was used in this study. Raw PE reads with an average length of 250 bp were generated, and all sequences were examined for possible sequencing errors. Adaptor sequences and low quality sequences (i.e. the percentage of bases of quality value ≤ 5 exceeds 50 % in the read) were trimmed with the FastQC program (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Short sequences (< 50 bp) were removed by a custom Perl program. The resulting high-quality sequences were de novo assembled into contigs and transcripts with Trinity software (https://github.com/trinityrnaseq/trinityrnaseq/wiki) [23]. To reduce data redundancy, transcripts with a minimum length of 200 bp were assembled and clustered by TGICL with default parameters. The longest sequences in each cluster were reserved and designated as unigenes [24]. All unigenes were used as queries in searching Nr, SwissProt and functionally annotated by GO analysis with Blast2GO software (E-value < 10-5) (https://www.blast2go.com/) [25]. KO analyses of differentially expressed unigenes (DEGs) in pathways were conducted on website (http://www.genome.jp/kegg/tool/map_pathway2.html). KEGG pathway analysis was performed using KEGG Automatic Annotation Server (KASS) with default parameters [26].

Identification of DEGs

All clean sequencing reads from each of the three libraries (A, B and C) were, respectively, mapped back to the transcriptome assembly by the software Bowtie2 with default setting. The number of reads aligned to each unigene in the alignment file was counted for each sample. These read counts were normalized as RPKM (Reads per kilobase of transcripts per million fragments mapped) values [27] and further analysis of identify DEGs among different groups was conducted by a web tool DESeq (http://www-huber.embl.de/users/anders/DESeq). The P-value was applied to determine DEGs, and the FDR (false discovery rate) method was applied to determine the threshold P-value in multiple tests to judge the significance of the differences of gene expression. In our analysis, if the FDR was < 0.05 and RPKM values showing at least a 2-fold change among samples, these unigenes were considered as significant DEGs. Finally, DEGs were tested by GO functional analysis and KEGG pathway analysis [28].

Results and discussion

Sequencing and de novo transcriptome assembly

To obtain the C. irritans tomont transcriptome expression profile during the dormant phase, the three cDNA libraries were constructed using tomonts from groups A, B and C. A total of 80,847,073 paired-end raw reads with an average length of 126 bp and a Q20 percentage higher than 96.02 % were generated (Additional file 1: Table S1). A total of 76,362,376 clean reads were obtained for subsequent analysis after eliminating low-quality sequences and adaptor sequences from the original data sequence by quality analysis (Additional file 1: Table S1). Transcript de novo assembly was performed for the clean reads by Trinity. A summary of all contigs, transcripts and unigenes assembly is presented in Table 1. The total length and number of contigs were 92,485,665 bp and 273,199, respectively. The maximum contig  length was 17,250 bp with an average length of 338.53 bp (N50:409), with the GC% of 46.05 %. The total length and number of transcripts were 84,205,775 bp and 162,496, respectively. The maximum length of transcript was 17,250 bp with an average length of 518 bp (N50:634), with GC% of 43.60 %. The total length and number of unigenes were 21,739,745 bp and 25,925, respectively. The maximal length of unigenes was 17,250 bp with an average length of 839 bp (N50:1250), and the GC% was 39.51 % (Table 1). Figure 1 demonstrated the distribution of the unigene lengths. The most abundant unigenes were clustered in a group with 200–299 bp in length. There were 742 unigenes with length of more than 3000 bp.

Table 1 Summary of the de novo assembly of transcriptomic profiles of C. ittitans tomonts
Fig. 1
figure 1

Transcriptome sequence length distributions of C. irritans tomont unigenes. The x-axis indicates unigene size and the y-axis indicates the number of unigenes with different lengths

EggNOG

For further functional prediction and classifications, unigenes were aligned to the eggNOG database. The eggNOG IDs of all genes were listed and these IDs were assigned into appropriated eggNOG categories, by which classified statistic was conducted to determine the functions of all genes, so as to learn the distribution patterns of gene functions of this species. A total of 23,670 (91.30 %) hits were annotated into 25,925 NR top hit unigenes (Table 2), which formed 25 classifications (Fig. 2). Among the functional classes, with the exception of the general function prediction only (3423, 11.84 %) and Function unknown (10.85 %), the largest three groups were post-translational modification, protein turnover and chaperones (10.78 %), signal transduction mechanisms (10.39 %), and translation, ribosomal structure and biogenesis (9.46 %). The smallest three groups were cell motility (0.12 %), nuclear structure (0.21 %), and extracellular structures (0.21 %) (Fig. 2).

Table 2 Annotation of unigenes of transcriptomic profiles of C. irritans tomonts
Fig. 2
figure 2

eggNOG function classification of C. irritans tomonts unigenes. A total of 23,670 hits were classified into 25 categories

Figure 2 presented the eggNOG analysis of DEGs from the three groups. In group B/A pairwise comparison, the up-regulated DEGs from eggNOG were functionally classified into 22 categories in which the three most enriched terms were J (99, Translation, ribosomal structure and biogenesis), O (26, Post-translational modification, protein turnover, chaperones) and R (26, General function prediction only); the down-regulated DEGs from eggNOG were functionally classified into 17 categories in which the three most enriched terms were B (14, Chromatin structure and dynamics), O (11) and S (11, Function unknown). In group B/C pairwise comparison, the up-regulated DEGs from eggNOG were functionally classified into 16 categories in which the 3 most enriched terms were the J (9), F (8, Nucleotide transport and metabolism), and R (4); the down-regulated DEGs from eggNOG were functionally classified into 22 categories in which the three most enriched terms were J (47), O (22), and S (14). In group C/A pairwise comparison, the regulated DEGs from eggNOG were functionally classified into 24 categories in which the three most enriched terms were J (197), O (62) and Z (57, Cytoskeleton); the down-regulated DEGs from eggNOG were functionally classified into 13 categories in which the three most enriched terms were O (8), S (4), and U (3, Intracellular trafficking, secretion and vesicular transport). It is found that some genes related to normal life activity such as protein synthesis regulations and cytoskeleton assembly (Fig. 2), were significantly inhibited by low-temperatures, which is different from ciliates originated from Antarctic region: the latter has developed a strong adaptability that allow tubulin post-translational modifications at low temperature [29].

Identification of GO enrichment analysis

Gene Ontology (GO) is a standardized gene functional classification system. In this study, a total of 21,510 unigenes (82.97 %) (Table 2) were assigned into three major functional GO terms (biological process category, cellular component category, and molecular function category), which were further summarized into 103 sub-categories using Blast2GO (Fig. 3).

Fig. 3
figure 3

Gene ontology (GO) enrichment analysis of the differently expressed genes in pairwise comparison groups B vs A, B vs C and C vs A. A: untreated blank control, B: treated for 24 h at 12 °C, and C: developed for 24 h at 25 °C

Figure 3 represents pie charts of GO enrichment analysis of DEGs from above three groups. In C/A pairwise comparison, the up-regulated DEGs from GO were functionally classified into 61 sub-categories in which the three most enriched sub-categories were the cell (289, GO:0005623), intracellular (289, GO:0005622), and organelle (261, GO:0043226). The previous study on C. irritans tomont transcriptome indicated that proteins were annotated as the external encapsulating structure, cell projection proteins and proteinaceous extracellular matrix under the GO category cellular component were those with potential as serodiagnostic markers of the tomont stage parasites [18]. In groups B/A pairwise comparison, the up-regulated DEGs from GO were functionally classified into 56 sub-categories in which the three most enriched sub-categories were the molecular function (129, GO:0003674), biosynthetic process (127, GO:0009058), and cell (120); the down-regulated DEGs from GO were functionally classified into 48 sub-categories in which the three most enriched sub-categories were the molecular function (40), biological process (33, GO:0008150), and cellular nitrogen compound metabolic process (32, GO:0034641). In groups B/C pairwise comparison, the up-regulated DEGs from GO were functionally classified into 38 sub-categories in which the three most enriched sub-categories were the molecular function (29), ion binding (26, GO:0043167), and biological process (22); the down-regulated DEGs from GO were functionally classified into 48 sub-categories in which the three most enriched sub-categories were the molecular function (72), cell (63), intracellular (63). In a study on Aurantiochytrium sp. under cold stress, “cellular process”, “binding” and “metabolic process” were also the largest proportion sub-categories in the three categories, and fatty acid biosynthetic processes were also influenced [30]. The results above revealed the importance of biological processes and molecular function of the protozoans under low temperature.

Functional classification by KEGG

To identify the biological pathways in C. irritans dormant tomonts, the unigenes were mapped to the reference pathways recorded in the KEGG database. A total of 6607 unigenes (25.49 %) (Table 2) were further annotated by KEGG and classified into six categories with 42 sub-classes (Fig. 4) (419 known KEGG pathways) (Additional files 2 and 3: Tables S2, S3).

Fig. 4
figure 4

KEGG pathway enrichment analysis of the differently expressed genes in pairwise comparison groups B vs A, B vs C and C vs A. The x-axis represents KEGG pathway classification and y-axis represents the significance P-value of enrichment in hypergeometric distribution. The red line represents P = 0.05. A: Metabolism; B: Genetic information processing; C: Environmental information processing; D: Cellular processes; E: Organismal systems; F: Human diseases

Figure 4 shows the histogram of KEGG enrichment analysis of DEGs from above three groups. Compared with the whole genome expression, with -log10 (P-value) used to represent the enrichment degree, the most differential expressed enriched subclasses were as follows: Translation (B:A = 11.19*; C:A = 13.25*; C:B = 1.64*), Nucleotide metabolism (7.63*; 0.02; 0.72), Substance dependence (5.40*; 1.17; 1.05), Immune diseases (4.62*; 0.38; 0), Neurodegenerative diseases (1.11; 2.46*; 1.30*), Sensory system (0; 1.99*; 3.77*), Digestive system (0.24; 1.89*; 1.97*), Endocrine system (0.30; 1.37*; 2.41*), Xenobiotics biodegradation and metabolism (2.04*; 0; 1.14), Amino acid metabolism (1.49*; 0; 0.09), Signal transduction (0.25; 0.53; 1.34*), Cell communication (0; 1.33*; 1.14). “*” indicated significant difference. The ribosome pathway was the first of the most DEGs-enriched pathway (Additional file 3: Table S3). This result was similar to the results of the previous study on C. irritans tomont transcriptome [18].

It can also be seen that although many genes encoding ribosomal proteins were up-regulated in both low temperature treated group (B, 63) and normally developed group (C, 139) compared with group A, the type and number of genes show great difference between groups B/A and groups C/A (Additional file 3: Table S3, Fig. 5). It is assumed that tomont cells under a different physiological condition would require different ribosomal proteins for the synthesis of a different set of proteins, and that these ribosomal proteins can also play various roles such as transcription, signal recognition, apoptosis and nuclear transport protein synthesis [18, 31].

Fig. 5
figure 5

Gene list involved in ribosome pathway generated by KEGG of the DEGs in pairwise comparison groups B vs A, B vs C and C vs A. Green indicates significantly increased expression; red, significantly decreased expression; and pink, unchanged expression. Blue denotes genes that were not identified in the expression profile analysis

DEGs

Water temperature influences the outbreak of cryptocaryoniasis, via promoting the growth or dormancy of C. irritans cells [4, 7]. To identify the DEGs involved in C. irritans response to low temperature treatment, pairwise comparisons for differential expression analysis were carried out among the three groups (A, B and C) (Fig. 6). Figure 7 demonstrates a heat map obtained from the hierarchical clustering of these genes.

Fig. 6
figure 6

Pairwise comparisons of differentially expressed genes in pairwise comparison groups B vs A, B vs C and C vs A. The red, blue, and yellow circles indicate DEGs in groups B vs A, C vs A and B vs C, respectively

Fig. 7
figure 7

Heat map obtained from the hierarchical clustering of these genes. Shades of red to green represent the up- and down- regulation of genes within each row, where red denotes up-regulation and green denotes down-regulation

In the group C/A pairwise comparison, there were 567 DEGs, including 548 up-regulated genes and 19 down-regulated genes in group C (Fig. 6, Additional file 4: Table S4 C/A). 96.65 % of the DEGs were up-regulated, indicating that despite absence of visible cell division in tomonts at this moment, numerous genes inside have been mobilized, preparing for cell division [32]. The top ten of DEGs identified as preferentially up-regulated at group C included Tubulin/FtsZ family, XP_003884402.1, chaperone protein, collagen triple helix protein, AF466826_1histone H3p, Ribosomal protein L18, myo-inositol 1-phosphate synthase, Pyridoxal-dependent decarboxylase conserved domain containing protein, 60S ribosomal protein L10, and Ribosomal protein L37a. The 19 down-regulated genes are presented in Additional file 4: Table S4 C/A. It is indicated that these 19 genes represent the remaining molecules [33], and may gradually become unnecessary as the cells approach division. Other proteins involved in the normal development of C. irritans tomonts will be described in detail in separate papers; in the present study these were mainly used to help identify proteins controlling tomont dormancy at low temperature.

In the group B/A pairwise comparison, 343 genes were differentially expressed in group B (Fig. 6, Additional file 4: Table S4 B/A columns A, L), i.e. 60.49 % (343/567) of these in the pairwise comparison group C/A, suggesting that tomonts of group B were still in a relatively active physiological condition after 24 h development at low temperature. This is because the tomonts were in a dynamic phase with rich internal life activities still going on instead of entering a fixed state immediately after the formation of tomonts. This was consistent with the study of Verni & Rosati [34]. However, most of the above DEGs (295/343) did not show significant differences in group C/A (Additional file 4: Table S4 B/A columns J, U), indicating a role of these genes in the low-temperature response, i.e. entering a dormant state.

Of all DEGs in the pairwise comparison group B/A, up to 77.26 % were up-regulated; of these, 40 were identical to the up-regulated DEGs in the pairwise comparison group C/A (Fig. 6, Additional file 4: Table S4 B/A), possibly involved in preparing for the restart of cell division and waking up when the water temperature rises again. Previous studies have confirmed that a mRNA pool appears in the dormant state of ciliates [35]. The stored mRNA such as actin, α-tubulin, HS70, glyceraldehydes-3-phosphate dehydrogenase, and metallothioneins in microbial cryptobiotic states might be required during the excystment process [33]. In the present study, the top ten out of 40 up-regulated DEGs included Ribosomal protein genes, Histone H3, and hypothetical and unnamed DEGs: XP_003884402.1, XP_004036516.1, XP_001426579.1, CCF73045.1. Among these DEGs, the highly expressed Ribosomal protein L18 (related) participates in stress adaptation by potentiating the cellular translation machinery to achieve a robust cytosolic stress response [36]. 60S ribosomal protein L10 (alpha/beta hammerhead) genes are involved in development and translation under stress, such as ultraviolet B (UV-B) [37]. Histone H3 was related to the aggregation of chromosomes [38]. 40S ribosomal protein may be a subject to transcriptional regulation in different parasite life-cycle stages [39]; and acidic ribosomal P0 protein is involved in cold-adaptation in ciliates [40]. Other up-regulated genes (Fig. 6, Additional file 4: Table S4 B/A) may be necessary for cells to survive at low temperature, or to enter a deeper level of dormancy. Of these, the highly expressed ribosomal protein L7 could be involved in the suppression of global protein synthesis and the cell cycle under cold conditions [41]; Bifunctional purine biosynthesis protein (putative) catalyzes the second and the fifth steps of de novo purine biosynthesis pathway [42]; Vacuolar ATP synthase in C. irritans theronts with higher expression might be involved in the process of host infection [2]; Heat shock protein 70 gene can improve protein folding together with some other chaperones [43]; 40S ribosomal protein S18 (putative) is involved in the initiation of translation [44]; Glutamine synthetaseis is involved in nitrogen metabolism, recycling of the neurotransmitter glutamate and glutamine biosynthesis for the production of amino acids, sugars, and glucosamine-6-phosphate [45]; Eukaryotic initiation factor 4A is a protein complex that mediates recruitment of ribosomes to mRNA [46]; Ubiquitin-conjugating enzyme (putative) is involved in selective protein degradation, DNA repair, cell cycle control, and possibly the regulation of chromatin structure [47]. The mechanism of these DEGs still requires further investigation, e.g. to determine the dynamic change in expression of the above-mentioned genes at different time points during low-temperature dormancy and waking up upon temperature rising.

We found that 78.38 % (145/185) of the DEGs were down-regulated in group B compared with group C (Fig. 6, Additional file 4: Table S4 B/C column L); of these, 135 (93.10 %) genes were up-regulated in pairwise comparison group C/A (Additional file 4: Table S4 B/C column T). In other words, all 548 up-regulated genes in group C/A (Additional file 4: Table S4 C/A column A) were down-regulated (135) or not significantly changed in group B/C (413) (Additional file 4: Table S4 C/A column I) suggesting that the genes associated with normal development and metabolism were significantly inhibited in group B after low-temperature treatment. This vital phenomenon, dormancy, is an important survival strategy in protozoans [34]. Every microbial cryptobiotic state is a direct consequence of the opening and closing of specific genes [33]. In the present study, the top ten down-regulated DEGs at group B/C included Tubulin beta chain, Tubulin/FtsZ family, chaperone protein, cathepsin L, calmodulin, heat shock protein 70, EWS75591.1, cathepsin b, XP_003882889.1, and cut up CG6998-PA (Additional file 4: Table S4). Tubulins are the major components of microtubule cytoskeletons in eukaryotic cells. Cytoskeleton is related to the morphological and structural changes of the cells, intracellular material transport, and dynamical system. In this study, Alpha/beta-tubulin and Tubulin/FtsZ family were highly up-regulated in the normally developed cells in group C because a large number of matrix proteins are essential for cells to form new cytoskeletons during cell division [2] but down-regulated by 999.28 and 590.19 times respectively in group B compared with group C. This confirmed that the activity associated with cell division was substantially inhibited. In addition, proteins associated with cytoskeleton assembly and regulation were also significantly down-regulated, such as calmodulin, cut up CG6998-PA, dynein gamma chain (flagellar outer arm), and ciliary dynein heavy chain, etc. And cathepsin is a type of protease related to the development, pathogenicity, as well as immune escape of the parasites [48].

In many ciliates, when environmental stress appears, this is identified by membrane receptors and opens a specific genetic programme that leads to dormant cyst formation [33]. In the group B/C pairwise comparison, eight of the 40 up-regulated genes were not significantly changed in group B compared with Group A (Additional file 4: Table S4 B/C column J), indicating a role in the regulation of tomonts entering a dormant state, including Ion (cation) channel family protein (ICFP), Succinate dehydrogenase iron-sulfur subunit (SDIS), Triosephosphate isomerase (TI), Malonyl coa-acyl carrier protein transacylase (MCCPT), Ribosomal protein L12 (RPL12); 50S ribosomal protein L13 (RPL13), XP_001443344.1 and XP_672154.1. SDIS is a primary dehydrogenase and donates electrons to the aerobic and energy-generating respiratory chain in eukaryotic mitochondria and numerous prokaryotes [49]; TI is a crucial enzyme in the glycolytic pathway, which catalyzes the interconversion between dihydroxyacetone phosphate (DHAP) and D-glyceraldehyde 3-phosphate (GAP) by an isomerization reaction without any cofactor or metal ion [50]; MCCPT is one of the essential enzymes in the FAS II system [51]; RPL12 is required for the ribosome to function with Elongation factor Tu, Initiation factor 2, and Reversing factor [52]; RPL13 is essential for the formation of the RI*50 intermediate particle [53]. However, the dynamic changes and functions of these genes during low-temperature dormancy and recovery of tomonts remain to be further studied.

The previous study on C. irritans transcriptomic analysis identified a total of 690 unique transcripts encoding excretory/secretory proteins and membrane and membrane-associated proteins [18]; the proteomic analysis considered α-tubulin, actin, enolase, vacuolar ATP synthase catalytic subunit, and hsp70 potential vaccine antigen candidates [2]. In the present results, abundant known or unknown potential genes for vaccine development were provided (Additional file 4: Table S4), and further studies should be conducted to certify their functions.

Conclusions

To our knowledge, this is the first transcriptomic study of C. irritans tomonts under low temperature conditions. It was concluded that most of the genes associated with normal cell development were inhibited at low temperature, but genes required for cell survival under low temperature, or for cell entry into a deeper dormancy state were significantly up-regulated. These genes might be considered as candidate genes to develop diagnostic and control measures for cryptocaryonosis. The results can serve as a reference for further studies on the mechanisms of low-temperature dormancy and recovery of parasites.

Ethics statement

Experiments on fish and parasites were performed according to the regulations of local and central government. All experiments were approved by the Institutional Animal Care and Use Committee of the East China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences.

References

  1. Colorni A, Diamant A. Ultrastructural features of Cryptocaryon irritans, a ciliate parasite of marine fish. Eur J Protistol. 1993;29(4):425–34.

    Article  CAS  PubMed  Google Scholar 

  2. Mai YZ, Li YW, Li RJ, Li W, Huang XZ, Mo ZQ, Li AX. Proteomic analysis of differentially expressed proteins in the marine fish parasitic ciliate Cryptocaryon irritans. Vet Parasitol. 2015;211(1–2):1–11.

    Article  CAS  PubMed  Google Scholar 

  3. How KH, Zenke K, Yoshinaga T. Dynamics and distribution properties of theronts of the parasitic ciliate Cryptocaryon irritans. Aquaculture. 2015;438:170–5.

    Article  Google Scholar 

  4. Cheung PJ, Nigrelli RF, Ruggieri GD. Studies on cryptocaryoniasis in marine fish: effect of temperature and salinity on the reproductive cycle of Cryptocaryon irritans Brown, 1951. J Fish Dis. 1979;2(2):93–7.

    Article  Google Scholar 

  5. Chen AP, Jiang YL, Qian D, Chen CF, Li AX, Huang J, Yang B. [Cryptocaryoniasis. China Fish.] 2011;8:39–40 (In Chinese).

    Google Scholar 

  6. Colorni A, Burgess P. Cryptocaryon irritans Brown 1951, the cause of ‘white spot disease’ in marine fish: an update. Aquar Sci Conserv. 1997;1(4):217–38.

    Article  Google Scholar 

  7. Dan XM, Lin XT, Yan YX, Teng N, Tan ZL, Li AX. A technique for the preservation of Cryptocaryon irritans at low temperatures. Aquaculture. 2009;297(1–4):112–5.

    Article  Google Scholar 

  8. Gu f, Ni B. [An ultrastructural study on resting cyst of Euplotes encysticus]. Acta Biologica Experimentalis Sinica. 1995;28(2):163–71 (In Chinese with English abstract).

  9. Uspenska Z. Drying and freezing rate related to effect of UV radiation of Colpoda-maupasi resting cysts under low-temperatures. Tsitologiya. 1972;14(4):487–93.

    Google Scholar 

  10. Shatilovich A, Stoupin D, Rivkina E. Ciliates from ancient permafrost: assessment of cold resistance of the resting cysts. Eur J Protistol. 2015;51(3):230–40.

    Article  PubMed  Google Scholar 

  11. Kremp A, Anderson DM. Factors regulating germination of resting cysts of the spring bloom dinoflagellate Scrippsiella hangoei from the northern Baltic Sea. J Plankton Res. 2000;22(7):1311–27.

    Article  Google Scholar 

  12. Genovesi B, Laabir M, Masseret E, Collos Y, Vaquer A, Grzebyk D. Dormancy and germination features in resting cysts of Alexandrium tamarense species complex (Dinophyceae) can facilitate bloom formation in a shallow lagoon (Thau, southern France). J Plankton Res. 2009;31(10):1209–24.

    Article  Google Scholar 

  13. Weisheit S, Villar M, Tykalová H, Popara M, Loecherbach J, Watson M, Růžek D, Grubhoffer L, de la Fuente J, Fazakerley JK et al. Ixodes scapularis and Ixodes ricinus tick cell lines respond to infection with tick-borne encephalitis virus: transcriptomic and proteomic analysis. Parasit Vectors. 2015;8(1):1–26.

    Article  Google Scholar 

  14. Valenzuela-Muñoz V, Sturm A, Gallardo-Escárate C. Transcriptomic insights on the ABC transporter gene family in the salmon louse Caligus rogercresseyi. Parasit Vectors. 2015;8(1):1–14.

    Article  Google Scholar 

  15. Kraeva N, Butenko A, Hlavacova J, Kostygov A, Myskova J, Grybchuk D, Lestinova T, Votypka J, Volf P, Opperdoes F et al. Leptomonas seymouri: adaptations to the dixenous life cycle analyzed by genome sequencing, transcriptome profiling and co-infection with leishmania donovani. PLoS Pathog. 2015;11(8):e1005127.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Maldonado-Aguayo W, Gallardo-Escarate C. Increasing transcriptome response of serpins during the ontogenetic stages in the salmon louse Caligus rogercresseyi (Copepoda: Caligidae). Marine Genomics. 2014;15:55–64.

    Article  CAS  PubMed  Google Scholar 

  17. Xiong J, Lu XY, Zhou ZM, Chang Y, Yuan DX, Tian M, Zhou ZG, Wang L, Fu CJ, Orias E et al. Transcriptome analysis of the model protozoan, Tetrahymena thermophila, using deep RNA sequencing. PLoS One. 2012; 7(2):e30630.

  18. Lokanathan Y, Mohd-Adnan A, Wan K-L, Nathan S. Transcriptome analysis of the Cryptocaryon irritans tomont stage identifies potential genes for the detection and control of cryptocaryonosis. BMC Genomics. 2010;11(1):1–15.

    Article  Google Scholar 

  19. Yin F, Gong H, Ke QZ, Li AX. Stress, antioxidant defence and mucosal immune responses of the large yellow croaker Pseudosciaena crocea challenged with Cryptocaryon irritans. Fish Shellfish Immunol. 2015;47:344-51.

  20. Li M, Pang Z, Xiao W, Liu X, Zhang Y, Yu D, Yang M, Yang Y, Hu J, Luo K. A transcriptome analysis suggests apoptosis-related signaling pathways in hemocytes of Spodoptera litura after parasitization by Microplitis bicoloratus. PLoS One. 2014;9(10):e110967.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Hu JW, You F, Wang Q, Weng SD, Liu H, Wang LJ, Zhang PJ, Tan XG. Transcriptional responses of olive flounder (Paralichthys olivaceus) to low temperature. PLoS One. 2014; 9(10):e108582.

  22. Zhang X, Wang S, Chen S, Chen Y, Liu Y, Shao C, Wang Q, Lu Y, Gong G, Ding S et al. Transcriptome analysis revealed changes of multiple genes involved in immunity in Cynoglossus semilaevis during Vibrio anguillarum infection. Fish Shellfish Immunol. 2015;43(1):209–18.

    Article  CAS  PubMed  Google Scholar 

  23. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, 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 Biotech. 2011;29(7):644–52.

    Article  CAS  Google Scholar 

  24. Huang L, Li GY, Mo ZL, Xiao P, Li J, Huang J. De novo assembly of the Japanese Flounder (Paralichthys olivaceus) spleen transcriptome to identify putative genes involved in immunity. PLoS One. 2015; 10(2):e0131146.

  25. Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.

    Article  CAS  PubMed  Google Scholar 

  26. Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M. KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010;38(Database issue):D355–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. 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.

    Article  CAS  Google Scholar 

  28. Zhang F, Guo ZN, Zhong H, Wang S, Yang WQ, Liu YF, Wang SH. RNA-Seq-based transcriptome analysis of aflatoxigenic aspergillus flavus in response to water activity. Toxins. 2014;6(11):3187–207.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Pucciarelli S, Ballarini P, Miceli C. Cold-adapted microtubules: characterization of tubulin posttranslational modifications in the Antarctic ciliate Euplotes focardii. Cell Motil Cytoskeleton. 1997;38(4):329–40.

    Article  CAS  PubMed  Google Scholar 

  30. Ma Z, Tan Y, Cui G, Feng Y, Cui Q, Song X. Transcriptome and gene expression analysis of DHA producer Aurantiochytrium under low temperature conditions. Sci Rep. 2015;5:14446.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Ojha S, Ahamad J, Bhattacharya A, Bhattacharya S. Ribosomal RNA and protein transcripts persist in the cysts of Entamoeba invadens. Mol Biochem Parasitol. 2014;195(1):6–9.

    Article  CAS  PubMed  Google Scholar 

  32. Breyne P, Dreesen R, Vandepoele K, De Veylder L, Van Breusegem F, Callewaert L, Rombauts S, Raes J, Cannoot B, Engler G et al. Transcriptome analysis during cell division in plants. Proc Natl Acad Sci U S A. 2002;99(23):14825–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Gutiérrez J, Callejas S, Borniquel S, Benítez L, Martín-González A. Ciliate cryptobiosis: a microbial strategy against environmental starvation. Int Microbiol Offic J Spanish Soc Microbiol. 2001;4(3):151–7.

    Google Scholar 

  34. Verni F, Rosati G. Resting cysts: a survival strategy in Protozoa Ciliophora. Italian J Zool. 2011;78(2):134–45.

    Article  Google Scholar 

  35. Tourancheau AB, Morin L, Yang T, Perasso R. Messenger RNA in dormant cells of Sterkiella histriomuscorum (Oxytrichidae): identification of putative regulatory gene transcripts. Protist. 1999;150(2):137–47.

    Article  CAS  PubMed  Google Scholar 

  36. Zhang X, Gao X, Coots RA, Conn CS, Liu B, Qian S-B. Translational control of the cytosolic stress response by mitochondrial ribosomal protein L18. Nat Struct Mol Biol. 2015;22(5):404–10.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Ferreyra MLF, Pezza A, Biarc J, Burlingame AL, Casati P. Plant L10 ribosomal proteins have different roles during development and translation under ultraviolet-B stress. Plant Physiol. 2010;153(4):1878–94.

    Article  CAS  PubMed Central  Google Scholar 

  38. Wei Y, Mizzen CA, Cook RG, Gorovsky MA, Allis CD. Phosphorylation of histone H3 at serine 10 is correlated with chromosome condensation during mitosis and meiosis in Tetrahymena. Proc Natl Acad Sci U S A. 1998;95(13):7480–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Oura CAL, Kinnaird J, Tait A, Shiels BR. Identification of a 40S ribosomal protein (S17) that is differentially expressed between the macroschizont and piroplasm stages of Theileria annulata. Int J Parasitol. 2002;32(1):73–80.

    Article  CAS  PubMed  Google Scholar 

  40. Pucciarelli S, Marziale F, Di Giuseppe G, Barchetta S, Miceli C. Ribosomal cold-adaptation: Characterization of the genes encoding the acidic ribosomal P0 and P2 proteins from the Antarctic ciliate Euplotes focardii. Gene. 2005;360(2):103–10.

    Article  CAS  PubMed  Google Scholar 

  41. Wu S, De Croos JNA, Storey KB. Cold acclimation-induced up-regulation of the ribosomal protein L7 gene in the freeze tolerant wood frog, Rana sylvatica. Gene. 2008;424(1–2):48–55.

    Article  CAS  PubMed  Google Scholar 

  42. Stoyanov A, Petrova P, Lyutskanova D, Lahtchev K. Structural and functional analysis of PUR2,5 gene encoding bifunctional enzyme of de novo purine biosynthesis in Ogataea (Hansenula) polymorpha CBS 4732 T. Microbiol Res. 2014;169(5–6):378–87.

    Article  CAS  PubMed  Google Scholar 

  43. Zininga T, Makumire S, Gitau GW, Njunge JM, Pooe OJ, Klimek H, Scheurr R, Raifer H, Prinsloo E, Przyborski JM et al. Plasmodium falciparum Hop (PfHop) interacts with the Hsp70 chaperone in a nucleotide-dependent fashion and exhibits ligand selectivity. PLoS One. 2015;10(8):e0135326.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Tanaka T, Tanaka M, Mitsui Y. The cloning and sequencing of ribosomal protein S18 of parasitic protozoa, Entamoeba histolytica. DNA Seq. 1998;8(4):267–70.

    Article  CAS  PubMed  Google Scholar 

  45. Smartt CT, Kiley LM, Hillyer JF, Dasgupta R, Christensen BM. Aedes aegypti glutamine synthetase: expression and gene structure. Gene. 2001;274(1–2):35–45.

    Article  CAS  PubMed  Google Scholar 

  46. Gingras AC, Raught B, Sonenberg N. eIF4 initiation factors: effectors of mRNA recruitment to ribosomes and regulators of translation. Ann Rev Biochem. 1999;68:913–63.

    Article  CAS  PubMed  Google Scholar 

  47. Okano S, Tokushima H, Nakaoka Y, Shimizu K. Cloning of a novel ubiquitin-conjugating enzyme (E2) gene from the ciliate Paramecium tetraurelia. Febs Lett. 1996;391(1–2):1–4.

    Article  CAS  PubMed  Google Scholar 

  48. Shin SP, Han SY, Han JE, Jun JW, Kim JH, Park SC. Expression and characterization of cathepsin L-like cysteine protease from Philasterides dicentrarchi. Parasitol Int. 2014;63(2):359–65.

    Article  CAS  PubMed  Google Scholar 

  49. Yu J, Wang S, Li W, Dai G, Xu S, He Z, Peng X, Zhou S, Liu X. Cloning, expression and protective immunity evaluation of the full-length cDNA encoding succinate dehydrogenase iron-sulfur protein of Schistosoma japonicum. Sci China Ser C. 2007;50(2):221–7.

    Article  CAS  Google Scholar 

  50. Kurkcuoglu Z, Findik D, Akten Ebru D, Doruker P. How an inhibitor bound to subunit interface alters triosephosphate isomerase dynamics. Biophys J. 2015;109(6):1169–78.

    Article  CAS  PubMed  Google Scholar 

  51. Sun M, Zhu G, Qin Z, Wu C, Lv M, Liao S, Qi N, Xie M, Cai J. Functional characterizations of malonyl-CoA:acyl carrier protein transacylase (MCAT) in Eimeria tenella. Mol Biochem Parasitol. 2012;184(1):20–8.

    Article  CAS  PubMed  Google Scholar 

  52. Brot N, Weissbach H. Chemistry and biology ofE. coli ribosomal protein L12. Mol Cell Biochem. 1981;36(1):47–63.

    Article  CAS  PubMed  Google Scholar 

  53. Metaxas A, Tzartos S, Liakopoulou-Kyriakides M. The production of anti-hexapeptide antibodies which recognize the S7, L6 and L13 ribosomal proteins of Escherichia coli. J Pept Sci. 2002;8(3):118–24.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was funded by the National Natural Science Foundation of China (grant no. 31101932) and the Special Scientific Research Funds Central Non-profit Institutes of China (grant no. 2014 M01, East China Sea Fisheries Research Institute; grant no. 2015B05XK01, Chinese Academy of Fishery Sciences).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Fei Yin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

FY conceived and designed the study; FY, PS, and JTW prepared the first draft of the manuscript; FY, JTW, and QXG provided strategic advice and assisted with editing of the manuscript. All authors read and approved the final version of the manuscript.

Additional files

Additional file 1: Table S1.

C. irritans tomont transcriptome expression profile under low temperature. (DOCX 15 kb)

Additional file 2: Table S2.

KEGG pathway enrichment analysis of the differently expressed genes. (DOCX 40 kb)

Additional file 3: Table S3.

KEGG pathway enrichment analysis of the differently expressed genes among low temperature treated (B) and untreated (A, C) C. irritans tomonts. (XLSX 16 kb)

Additional file 4: Table S4.

Differently expressed genes among low temperature treated (B) and untreated (A, C) C. irritans tomonts. (XLSX 188 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yin, F., Sun, P., Wang, J. et al. Transcriptome analysis of dormant tomonts of the marine fish ectoparasitic ciliate Cryptocaryon irritans under low temperature. Parasites Vectors 9, 280 (2016). https://0-doi-org.brum.beds.ac.uk/10.1186/s13071-016-1550-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13071-016-1550-1

Keywords