Cryptocaryon irritans infection induces the acute phase response in Lates calcarifer: a transcriptomic perspective

Transcriptome Analysis Reveals the Immunosuppression in Tiger Pufferfish (Takifugu rubripes) under Cryptocaryon irritans Infection

The tiger pufferfish (Takifugu rubripes), also known as fugu, has recently suffered from severe C. irritans infections under aquaculture environment, yet the underlying immune mechanisms against the parasite remain poorly understood. In this study, we conducted a comprehensive transcriptome analysis of the gill tissue from infected and uninfected fish using PacBio long-read (one pooled sample each for seriously infected and healthy individuals, respectively) and Illumina short-read (three pools for mildly infected, seriously infected, and healthy individuals, respectively) RNA sequencing technologies. After aligning sequence data to fugu's reference genome, 47,307 and 34,413 known full-length transcripts were identified and profiled in healthy and infected fish, respectively. Similarly, we identified and profiled 1126 and 803 novel genes that were obtained from healthy and infected fish, respectively. Interestingly, we found a decrease in the number of alternative splicing (AS) events and long non-coding RNAs (lncRNAs) after infection with C. irritans, suggesting that they may be involved in the regulation of the immune response in fugu. There were 687 and 1535 differentially expressed genes (DEGs) in moderately and heavily infected fish, respectively, compared to uninfected fish. Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that immune-related DEGs in the two comparison groups were mainly enriched in cytokine-cytokine receptor interactions, ECM-receptor interactions, T-cell receptor signaling pathways, Th1 and Th2 cell differentiation, and Th17 cell differentiation pathways. Further analysis revealed that a large number of immune-related genes were downregulated in infected fish relative to uninfected ones, such as CCR7, IL7R, TNFRSF21, CD4, COL2A1, FOXP3B, and ITGA8. Our study suggests that C. irritans is potentially a highly efficient parasite that may disrupt the defense mechanisms of fugu against it. In addition, in combination of short-read RNA sequencing and previous genome-wide association analyses, we identified five key genes (NDUFB6, PRELID1, SMOX, SLC25A4, and DENND1B) that might be closely associated with C. irritans resistance. This study not only provides valuable resources of novel genic transcripts for further research, but also provides new insights into the immune mechanisms underlying C. irritans infection response in farmed fugu.

Transcriptomic analysis of large yellow croaker (Larimichthys crocea) reveals the suppression of the inflammatory response from Cryptocaryon irritans infection

Large yellow croaker (Larimichthys crocea) is the most productive marine fish in China. Cryptocaryon irritans is an extremely destructive parasite that causes great economic losses in large yellow croaker aquaculture industry. Therefore, it is very necessary to study the immune response of large yellow croaker in response to C. irritans infection. In this study, the transcriptomic profiles of large yellow croaker were sequenced and analyzed in the brain and head kidney at 72 h after C. irritans infection. Cytokines and chemokines related terms were significantly enriched based on the GO enrichment of down-regulated differentially expressed genes (DEGs) from the head kidney. Meanwhile, cytokine-cytokine receptor interaction was significantly enriched based on the KEGG enrichment of up-regulated DEGs from the brain and down-regulated DEGs from the head kidney, respectively. Moreover, the majority of inflammation-related DEGs were significantly up-regulated in the brain, but distinctly down-regulated in the head kidney. These results showed that the brain and head kidney might play different roles against C. irritans infection, and the inflammatory response of large yellow croaker may be restrained during C. irritans infection. Taken together, the transcriptomic analyses will be helpful to more comprehensively understand the immune mechanism of teleost against C. irritans infection, and provide a theoretical basis for the prevention and treatment of Cryptosporidiosis.

Growth, sex reversal pattern, and reproductive characteristics of Barramundi (Lates calcarifer) broodstock candidates reared in floating cages

Sex reversal of male to female is a characteristic of barramundi (Lates calcarifer), which is affected by several factors, thereby changing the broodstock population. A study was conducted in floating cages in Langkawi, Malaysia, to determine the weight point at the onset of the sex reversal phenomena. A total of 75 female and 55 male adult individuals (3-4 weeks of age) were sampled from the fish cultured in cages to ascertain their sex at different weights. The water temperature and salinity values were 29.82 °C and 33.12 ppt, respectively. The specimens were classified into twelve bodyweight classes (2.00-8.00 ± 0.5 kg intervals). Female specimen body weight distribution was highest in the 6.01-6.50 kg class (22.6%), followed by the 5.51-6.00 kg and 4.51-5.00 class (13.3%), while male specimen body weight distribution was highest in the 4.51-5.00 kg class (32.1%), followed by the 4.01-4.50 kg class (30.3%). Length-to-weight relationships for females and males of Asian Seabass indicated positive allometric growth. The correlation between body weight and GSI, using Pearson's correlation, for both sexes, for the male and female barramundi, there was a weak correlation between body weight and GSI, which was 37 and 30%, respectively. Based on the present study's findings, it can be concluded that sex reversal from male to female in Barramundi largely occurred at 4.57 kg body weight and 66.8 cm total length.

Comprehensive transcriptomics and proteomics analysis of Carassius auratus gills in response to Aeromonas hydrophila

As one of the mucosal barriers, fish gills represent the first line of defense against pathogen infection. However, the exact mechanism of gill mucosal immune response to bacterial infection still needs further investigation in fish. Here, to investigate pathological changes and molecular mechanisms of the mucosal immune response in the gills of crucian carp (Carassius auratus) challenged by Aeromonas hydrophila, the transcriptomics and proteomics were performed by using multi-omics analyses of RNA-seq coupled with iTRAQ techniques. The results demonstrated gill immune response were mostly related to the activation of complement and coagulation cascades, antigen processing and presentation, phagosome, NOD-like receptor (NLR) and nuclear factor κB (NFκB) signaling pathway. Selected 21 immune-related DEGs (ie., Clam, nfyal, snrpf, acin1b, psme, sf3b5, rbm8a, rbm25, prpf18, g3bp2, snrpd3l, tecrem-2, cfl-A, C7, lysC, ddx5, hsp90, α-2M, C9, C3 and slc4a1a) were verified for their immune roles in the A. hydrophila infection via using qRT-PCR assay. Meanwhile, some complement (C3, C7, C9, CFD, DF and FH) and antigen presenting (HSP90, MHC Ⅱ, CALR, CANX and PSME) proteins were significantly participated in the process of defense against infections in gill tissues, and protein-protein interaction (PPI) network displayed the immune signaling pathways and interactions among these DEPs. The correlation analysis indicated that the iTRAQ and qRT-PCR results was significantly correlated (Pearson's correlation coefficient = 0.70, p < 0.01). To our knowledge, the transcriptomics and proteomics of gills firstly identified by multi-omics analyses contribute to understanding on the molecular mechanisms of local mucosal immunity in cyprinid species.

Investigating the etiologies behind emergent mass mortalities of farmed Liza carinata juveniles from coastal farms at Damietta, Egypt

This study aimed to identify the mortality present in private fish farm Amyloodinium ocellatum and Cryptocaryon irritans were isolated from this outbreak affecting Liza carinata fingerlings at an earthen-based aquaculture facility in Damietta, Egypt. A total of 140 moribunds, L. carinata, were collected from the fish ponds during the mortality events. Physico-chemical analysis of water was analyzed. The skin, fins, gills, and eyes of each fish specimen were scraped gently onto slides in areas over 2 cm area. All smears were examined separately under the light microscope. Molecular identification of the parasites using analysis of ITS rDNA regions flanking both 18S and 28S rDNA genes of Amyloodinium protozoa and C. irritans. Identities of the detected parasites were confirmed by gene sequence and phylogenetic analysis. The majority of the examined fish (90%) were infected, 66.42% had a mixed infection, and 23.57% had a single infection either with A. ocellatum (10.71%) or C. irritans (12.85%).The mean intensity of A. ocellatum was 16.5 ± 2.03 in the skin and 13.18 ± 1.90 in the gills of infected fish, while that of C. irritans was 4.75 ± 1.05 in gills and 7.43 ± 1.45 in the skin, respectively. To control the emergent mortalities, affected ponds were treated using copper sulfate pentahydrate, hydrogen peroxides solutions, and amprolium hydrochloride powder in feed. Fish across the treated ponds were gradually improved with low morbidity and mortalityrates during the treatment period. The clinical disease was almost diminished at the end of the second week of treatment. Coinciding with the clinical improvement of the treated juveniles, microscopical examination of skin/gill scraps exhibited a marked decline in the number of protozoan parasites at the end of the second week of treatment.

Outbreak of Cryptocaryon irritans infection in silver pomfret Pampus argenteus cultured in China

Silver pomfret Pampus argenteus is a major cultivated marine fish species with a high market value. In summer 2021, Cryptocaryon irritans, a ciliate parasite, infected the cultured silver pomfret in aquaculture ponds in Ningbo, Zhejiang Province, China. The symptoms of infected fish include white spots on the skin and fins, increased body surface mucus, loss of appetite, irritability, and shedding of scales. After collecting white spots from moribund fish, the 18S ribosomal RNA sequence of the pathogen on the fish skin was amplified by PCR; phylogenetic analysis showed that it was closely related to C. irritans strains from Ningde, Fujian, China. Four groups of silver pomfret were tested in an artificial infection experiment over the course of 72 h, consisting of 3 infected groups (1600, 4000, and 8000 theronts fish-1) and 1 healthy group. White spots were observed on the skin and fins of the infected fish, but not on their gills. Samples were taken from the gills, liver, kidney, and spleen of both infected and healthy fish and were compared to evaluate any significant histopathological differences. As the dose of infection increased, symptoms became more pronounced. At 72 h, mortality rates were 8.3, 50, and 66.7% for the 3 different concentrations, respectively. The median lethal concentration was calculated to be 366 theronts g-1 at 72 h, 298 theronts g-1 at 84 h, and 219 theronts g-1 at 96 h. This study emphasizes the importance of developing early diagnosis methods and appropriate prevention strategies to decrease the impact of C. irritans infection in the silver pomfret aquaculture industry.

Severe Natural Outbreak of Cryptocaryon irritans in Gilthead Seabream Produces Leukocyte Mobilization and Innate Immunity at the Gill Tissue

The protozoan parasite Cryptocaryon irritans causes marine white spot disease in a wide range of fish hosts, including gilthead seabream, a very sensitive species with great economic importance in the Mediterranean area. Thus, we aimed to evaluate the immunity of gilthead seabream after a severe natural outbreak of C. irritans. Morphological alterations and immune cell appearance in the gills were studied by light microscopy and immunohistochemical staining. The expression of several immune-related genes in the gills and head kidney were studied by qPCR, including inflammatory and immune cell markers, antimicrobial peptides (AMP), and cell-mediated cytotoxicity (CMC) molecules. Serum humoral innate immune activities were also assayed. Fish mortality reached 100% 8 days after the appearance of the C. irritans episode. Gill filaments were engrossed and packed without any space between filaments and included parasites and large numbers of undifferentiated and immune cells, namely acidophilic granulocytes. Our data suggest leukocyte mobilization from the head kidney, while the gills show the up-regulated transcription of inflammatory, AMPs, and CMC-related molecules. Meanwhile, only serum bactericidal activity was increased upon infection. A potent local innate immune response in the gills, probably orchestrated by AMPs and CMC, is triggered by a severe natural outbreak of C. irritans.

Characteristics delineation of piscidin 5 like from Larimichthys crocea with evidence for the potent antiparasitic activity

Several researches reported that piscidin members of teleosts owned strong antiparasitic activity. Cryptocaryon irritans, a type of ectoparasite, could infect most of the marine teleosts. Larimichthys crocea could severely suffer from marine white spot disease caused by C. irritans, and their mortality rate was significantly high. Concentrating on this problem, we have done many related works. Piscidin 5 like (termed Lc-P5L) was another piscidin member isolated from a comparative transcriptome of C. irritans-immuned L. crocea. In the paper, quantitative Real-time PCR (qRT-PCR) showed Lc-P5L was upregulated in examined tissues, including gill, head kidney, muscle, liver, spleen and intestine after challenged by C. irritans, the significant upregulation time was in accordance to key developmental stages of C. irritans, which implied different infection stages could result in host immune response. Furthermore, using microscope techniques, we observed theronts or trophonts became weakly motile, cilia became detached, cells were out of shape, membranes eventually lysed in different cell positions and cytoplasmic contents leaked. Laser confocal scanning microscope (LCSM) observed theronts macronucleus grew swell and depolymerized after treated by recombinant Lc-P5L (rLc-P5L). Data suggested rLc-P5L was significantly lethal to C. irritans, and the death state of the parasite incubated with rLc-P5L was remarkably similar to other piscidin members or other antiparasitic peptides (APPs). Thus, these data provided new insights into L. crocea immunity against C. irritans and potential of rLc-P5L as a therapeutic agent against pathogen invasion.

Molecular characterization and antiparasitic activity analysis of a novel piscidin 5-like type 4 from Larimichthys crocea

Cryptocaryon irritans is an obligate parasitic ciliate protozoan that can infect various commercially important mariculture teleosts and cause high lethality and economic loss, especially Larimichthys crocea. Current methods of controlling or preventing this parasite with chemicals or antibiotics are widely considered to be environmentally harmful. The antiparasitic activity of some antimicrobial peptides (AMPs) attracted extensive attention of scholars. In the study, a novel piscidin 5-like type 4 (termed Lc-P5L4) excavated from comparative transcriptome of C. irritans - immuned L. crocea was identified and characterized. Sequence analysis shows the full-length cDNA of Lc-P5L4 is 539 bp containing an open reading frame (ORF) of 198 bp which encodes a peptide of 65 amino acid residues. The genome consists of three exons and two introns which exist in its ORF, and all the exon-intron boundaries are in accordance with classical GT-AG rule (GT/intron/AG). Multiple alignments indicate the signal peptides share highly conserved identity, while mature peptides are more diverse. Phylogenetic analysis displays Lc-P5L4 clusters together with other members of piscidin 5-like family. Next, quantitative Real-time PCR (qRT-PCR) detection found C. irritans infection could upregulate Lc-P5L4 expression level in all tested tissues significantly, it appeared earliest upregulation in the theronts infection stage in the head kidney; the expression contents reached to maximum level in the intestine, gill and muscle during trophonts falling off stage; while it was just upregulated during secondary bacterial infection stage in the liver and spleen. The data showed Lc-P5L4 upregulation time points were in accordance with different infection stages. With recombinant Lc-P5L4 (rLc-P5L4) obtained through Escherichia coli system, in vitro assay showed rLc-P5L4 could cause cilia deactivation, cell bodiesclumping and sticking to each other, then cell membrane rupture and contents leakage. The data illustrated Lc-P5L4 played critical roles in the immune defense against C. irritans infection, and provided another proof that piscidins exhibit multiple anti- C. irritans features.

Infectivity and genes differentially expressed between young and aging theront cells of the marine fish parasite Cryptocaryon irritans

The ciliated protozoan Cryptocaryon irritans infects a wide range of marine fish and causes the highly lethal white spot disease. This parasite possesses three morphologically and physiologically distinct life stages: an infectious theront, a parasitic trophont, and an asexually reproductive tomont. In the past few years, several attempts have been made to help elucidate how C. irritans transforms from one stage to another using transcriptomic or proteomic approaches. However, there has been no research studying changes in transcription profiles between different time points of a single C. irritans life stage—the development of this parasite. Here we use RNA-seq and compare gene expression profiles of theront cells collected by 1 and 10 hrs after they emerged from tomonts. It has been shown that infectivity of theront cells declines 6–8 hours post-emergence, and we used this characteristic as a physiological marker to confirm the aging of theront cells. We identified a total of 41 upregulated and 90 downregulated genes that were differentially expressed between young and aging theront cells. Using Blast2Go to further analyze functions of these genes, we show that genes related to energy production are downregulated, but quite surprisingly many genes involved in transcription/translation processes are upregulated. We also show that expression of all nine detectable agglutination/immobilization antigen genes, with great sequence divergence, is invariably downregulated. Functions of other differentially expressed genes and indications are also discussed in our study.

Functional Analysis of IRF1 Reveals its Role in the Activation of the Type I IFN Pathway in Golden Pompano, Trachinotus ovatus (Linnaeus 1758)

Interferon (IFN) regulatory factor 1 (IRF1), a transcription factor with a novel helix–turn–helix DNA-binding domain, plays a crucial role in innate immunity by regulating the type I IFN signaling pathway. However, the regulatory mechanism through which IRF1 regulates type I IFN in fish is not yet elucidated. In the present study, IRF1 was characterized from golden pompano, Trachinotus ovatus (designated ToIRF1), and its immune function was identified to elucidate the transcriptional regulatory mechanism of ToIFNa3. The full-length complementary DNA (cDNA) of IRF1 is 1763 bp, including a 900-bp open reading frame (ORF) encoding a 299-amino-acid polypeptide. The putative protein sequence has 42.7–71.7% identity to fish IRF1 and possesses a representative conserved domain (a DNA-binding domain (DBD) at the N-terminus). The genomic DNA sequence of ToIRF1 consists of eight exons and seven introns. Moreover, ToIRF1 is constitutively expressed in all examined tissues, with higher levels being observed in immune-relevant tissues (whole blood, gill, and skin). Additionally, Cryptocaryon irritans challenge in vivo increases ToIRF1 expression in the skin as determined by Western blotting (WB); however, protein levels of ToIRF1 in the gill did not change significantly. The subcellular localization indicates that ToIRF1 is localized in the nucleus and cytoplasm with or without polyinosinic/polycytidylic acid (poly (I:C)) induction. Furthermore, overexpression of ToIRF1 or ToIFNa3 shows that ToIRF1 can notably activate ToIFNa3 and interferon signaling molecule expression. Promoter sequence analysis finds that several interferon stimulating response element (ISRE) binding sites are present in the promoter of ToIFNa3. Additionally, truncation, point mutation, and electrophoretic mobile shift (EMSA) assays confirmed that ToIRF1 M5 ISRE binding sites are functionally important for ToIFNa3 transcription. These results may help to illuminate the roles of teleost IRF1 in the transcriptional mechanisms of type I IFN in the immune process.

Immune response of mollusk Onchidium struma to extremely low-frequency electromagnetic fields (ELF-EMF, 50 Hz) exposure based on immune-related enzyme activity and De novo transcriptome analysis

Along with rapid offshore and onshore wind power development in modern society, extremely low frequency electromagnetic fields (ELF-EMF) is produced extensively in the habits of aquatic organisms. However, the biological effects of ELF-EMF on aquatic organisms are almost sparse. In this study, Onchidium struma without shell was chosen to aim whether ELF-EMF can elicit immune response of mollusk based on immune-related enzyme activities and gene expression through high-throughput transcriptome sequencing. Three experimental groups, i.e. ELF-EMF unexposed control group (C), ELF-EMF (50 Hz, 100 μT) exposed E1 group, and ELF-EMF (50 Hz, 500 μT) exposed E2 group, were set, and coelomocytes were collected to analyze. The results showed that total coelomocyte and spherulocyte density in E1 group increased significantly compared to groups C and E2 (P < 0.05). There were no significant differences on amoebocyte and chromatocyte density among groups C, E1 and E2. ELF-EMF exposure could significantly increase immune-related enzyme activities in coelomic fluid of O. struma, including acidic phosphatase, alkaline phosphatase, antioxidative capacity, catalase, superoxide dismutase, and polyphenol oxidase (P < 0.05). A total of 54.32 Mb and 55.27 Mb raw reads with average length of 1520 bp were obtained from coelomocytes of O. struma in unexposed and exposed groups, respectively. There were 341 differentially expressed genes (DGEs) between unexposed and exposed groups, including 209 up-regulated and 132 down-regulated unigenes. All the DGEs were allocated to 14 Kyoto Encyclopedia of Genes and Genomes pathways, and five pathways were associated with immune response, including TLR/TNF/NOD-like receptor/MAPK/Fc epsilon RI signaling pathways. Altogether, short-term (to one week) exposure of O. struma to lower luxy density ELF-EMF (<500 μT) could elicit the immune response, and antioxidant system is recommended as indicators of immunological effects. Hopefully, this study will further provide insights into exploring biomarker for evaluation of the effect of ELF-EMF exposure on aquatic organisms regarding to field density, frequency and exposure duration, and provide good guidance for exploitation and utilization of renewable energy.

Skin metabolome reveals immune responses in yellow drum Nibea albiflora to Cryptocaryon irritans infection

The yellow drum Nibea albiflora is less susceptible to Cryptocaryon irritans infection than is the case with other marine fishes such as Larimichthys crocea, Lateolabrax japonicus, and Pagrus major. To investigate further their resistance mechanism, we infected the N. albiflora with the C. irritans at a median lethal concentration of 2050 theronts/g fish. The skins of the infected and the uninfected fishes were sampled at 24 h and 72 h followed by an extensive analysis of metabolism. The study results revealed that there were 2694 potential metabolites. At 24 h post-infection, 12 metabolites were up-regulated and 17 were down-regulated whereas at 72 h post-infection, 22 metabolites were up-regulated and 26 were down-regulated. Pathway enrichment analysis shows that the differential enriched pathways were higher at 24 h with 22 categories and 58 subcategories (49 up, 9 down) than at 72 h whereby the differential enriched pathways were 6 categories and 8 subcategories (4 up, 4 down). In addition, the principal component analysis (PCA) plot shows that at 24 h the metabolites composition of infected group were separately clustered to uninfected group while at 72 h the metabolites composition in infected group were much closer to uninfected group. This indicated that C. irritans caused strong metabolic stress on the N. albiflora at 24 h and restoration of the dysregulated metabolic state took place at 72 h of infection. Also, at 72 h post infection a total of 17 compounds were identified as potential biomarkers. Furthermore, out of 2694 primary metabolites detected, 23 metabolites could be clearly identified and semi quantified with a known identification number and assigned into 66 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Most of the enriched KEGG pathways were mainly from metabolic pathway classes, including the metabolic pathway, biosynthesis of secondary metabolites, taurine and hypotaurine metabolism, purine metabolism, linoleic acid metabolism, phenylalanine, tyrosine and tryptophan biosynthesis. Others were glyoxylate and dicarboxylate metabolism, glutathione metabolism, and alanine, aspartate, and glutamate metabolism. Moreover, out of the identified metabolites, only 6 metabolites were statistically differentially expressed, namely, L -glutamate (up-regulated) at 24 h was important for energy and precursor for other glutathiones and instruments of preventing oxidative injury; 15-hydroxy- eicosatetraenoic acid (15-HETE), (S)-(-)-2-Hydroxyisocaproic acid, and adenine (up-regulated) at 72 h were important for anti-inflammatory and immune responses during infection; others were delta-valerolactam and betaine which were down-regulated compared to uninfected group at 72 h, might be related to immure responses including stimulation of immune system such as production of antibodies. Our results therefore further advance our understanding on the immunological regulation of N. albiflora during immune response against infections as they indicated a strong relationship between skin metabolome and C. irritans infection.

Integrating Genomic and Morphological Approaches in Fish Pathology Research: The Case of Turbot (Scophthalmus maximus) Enteromyxosis

Enteromyxosis, caused by Enteromyxum scophthalmi, is one of the most devastating diseases stemming from myxozoan parasites in turbot (Scophthalmus maximus L.), being a limiting factor for its production. The disease develops as a cachectic syndrome, associated to catarrhal enteritis and leukocytic depletion, with morbidity and mortality rates usually reaching 100%. To date, no effective treatment exists and there are different unknown issues concerning its pathogenesis. The gross and microscopic lesions associated to enteromyxosis have been thoroughly described, and several morphopathological studies have been carried out to elucidate the mechanisms of this host-parasite interaction. More recently, efforts have been focused on a multidisciplinary approach, combining histopathology and transcriptome analysis, which has provided significant advances in the understanding of the pathogenesis of this parasitosis. RNA-Seq technology was applied at early and advanced stages of the disease on fishes histologically evaluated and classified based on their lesional degree. In the same way, the transcriptomic data were analyzed in relation to the morphopathological picture and the course of the disease. In this paper, a comprehensive review of turbot enteromyxosis is presented, starting from the disease description up to the most novel information extracted by an integrated approach on the infection mechanisms and host response. Further, we discuss ongoing strategies toward a full understanding of host-pathogen interaction and the identification of suitable biomarkers for early diagnosis and disease management strategies.

Transcriptome analysis and discovery of genes involved in immune pathways from coelomocytes of Onchidium struma after bacterial challenge

Onchidium struma widely distributes in subtidal and low-tidal zones, which is considered to be an economical species with rich nutrition, a valuable biomonitor for heavy metal pollution and a representative species for evolution from ocean to land. However, there is limited genetic information available for O. struma development. This study compared transcriptomic profiles of coelomocytes from normal and bacteria infected O. struma by Illumina-based paired-end sequencing to explore the molecular immune mechanism of O. struma against bacterial infection. After assembly, a total of 92,450 unigenes with an average length of 1019 bp were obtained. Approximately 34,964 (37.82%) unigenes were annotated in the Nr NCBI database and 40.1% of unigenes were similar with that of Aplysia californica. Among them, 7609 unigenes were classified into three Gene Ontology (GO) categories: biological process (3250 unigenes, 42.7%), cellular component (2,281, 30.0%) and molecular function (2078 unigenes, 27.3%). A total of 22,776 unigenes were aligned to the Clusters of Orthologous Groups (COG) of proteins and classified into 25 functional categories. Following bacterial infection, 10,623 differently expressed unigenes (DEGs) were identified, including 7644 up-regulated and 2979 down-regulated unigenes. Further KEGG analysis annotated 11,681 DEGs to 42 pathways, and 11 pathways were identified to be related with diseases and immune system. To our knowledge, it was first time to analyze transcriptome profiles of O. struma. Results of the present study will provide valuable theoretical resources for future genetic and genomic research on O. struma. The research results will be helpful for improving the efficiency and quality of artificial breeding, establishing genetic linkage map, and enhancing health management for this species.

Transcriptomic analysis reveals the key immune-related signalling pathways of Sebastiscus marmoratus in response to infection with the parasitic ciliate Cryptocaryon irritans

Background

False kelpfish (Sebastiscus marmoratus) is one of the target species in artificial breeding in China, and is susceptible to infection by Cryptocaryon irritans, which is an obligate parasitic ciliate that lives in the epithelium of the fish gills, skin and fins. Here, we sought to understand the mechanisms of molecular immunity of S. marmoratus against C. irritans infection.

Methods

We carried out an extensive analysis of the transcriptome of S. marmoratus immune-related tissues. A paired-end library was constructed from the cDNA synthesized using a Genomic Sample Prep Kit. Five normalized cDNA libraries were constructed using RNA from the control group and the four groups of C. irritans-infected fish. The libraries were sequenced on an Illumina Mi-Seq platform, and functional annotation of the transcriptome was performed using bioinformatics software.

Results

The data produced a total of 149,983,397 clean reads from five cDNA libraries constructed from S. marmoratus immune-related tissues. A total of 33,291 unigenes were assembled with an average length of 1768 bp. In eggNOG (Evolutionary Genealogy of Genes: non-supervised orthologous groups) categories, 333 unigenes (0.94%) were assigned to defense mechanisms. In the immune system process sub-categories of gene ontology (GO) enrichment analysis, with the passage of time post-infection, the number of differentially expressed genes (DEGs) was reduced from 24 h to 48 h but then increased from 72 h to 96 h. Specifically, the immune-related differentially expressed genes (IRDEGs), which belong to the KEGG (Kyoto encyclopedia of genes and genomes) pathways, such as the complement and coagulation cascades, chemokine signalling pathways and toll-like receptor signalling pathways were mainly observed at 24 h post-infection.

Conclusions

Infection with C. irritans resulted in a large number of DEGs in the immune-related tissues of S. marmoratus. The rapid and significant response of the S. marmoratus immune signalling pathways following C. irritans infection may be associated with their involvement in the immune process.

Electronic supplementary material

The online version of this article (10.1186/s13071-017-2508-7) contains supplementary material, which is available to authorized users.

Effect of high temperature on immune response of grass carp (Ctenopharyngodon idellus) by transcriptome analysis

Exposure to high temperature is an inherent feature of grass carp culture in southern China and juvenile grass carps are predisposed to infectious disease in this condition. To understand how high temperature impacts the immune response to pathogens in grass carp, the transcriptomic profiles of the spleens from immune injected grass carp groups undergoing heat stress and normal temperature were investigated. An average of 72 million clean reads per library was obtained, and approximate 80% of these genes were successfully mapped to the reference genome. A total of 2287 up-regulated and 1068 down-regulated genes were identified. 10 immune-related categories involving 90 differently expressed genes were scrutinized. Expression patterns of 18 differentially expressed genes involving in immune response were validated by quantitative real-time RT-PCR. These results provide further significant insights into the influence mechanism of high temperature to immune response in grass carp.

RNA-seq analysis of early enteromyxosis in turbot (Scophthalmus maximus): new insights into parasite invasion and immune evasion strategies

Enteromyxum scophthalmi, an intestinal myxozoan parasite, is the causative agent of a threatening disease for turbot (Scophthalmus maximus, L.) aquaculture. The colonisation of the digestive tract by this parasite leads to a cachectic syndrome associated with high morbidity and mortality rates. This myxosporidiosis has a long pre-patent period and the first detectable clinical and histopathological changes are subtle. The pathogenic mechanisms acting in the early stages of infection are still far from being fully understood. Further information on the host-parasite interaction is needed to assist in finding efficient preventive and therapeutic measures. Here, a RNA-seq-based transcriptome analysis of head kidney, spleen and pyloric caeca from experimentally-infected and control turbot was performed. Only infected fish with early signs of infection, determined by histopathology and immunohistochemical detection of E. scophthalmi, were selected. The RNA-seq analysis revealed, as expected, less intense transcriptomic changes than those previously found during later stages of the disease. Several genes involved in IFN-related pathways were up-regulated in the three organs, suggesting that the IFN-mediated immune response plays a main role in this phase of the disease. Interestingly, an opposite expression pattern had been found in a previous study on severely infected turbot. In addition, possible strategies for immune system evasion were suggested by the down-regulation of different genes encoding complement components and acute phase proteins. At the site of infection (pyloric caeca), modulation of genes related to different structural proteins was detected and the expression profile indicated the inhibition of cell proliferation and differentiation. These transcriptomic changes provide indications regarding the mechanisms of parasite attachment to and invasion of the host. The current results contribute to a better knowledge of the events that characterise the early stages of turbot enteromyxosis and provide valuable information to identify molecular markers for early detection and control of this important parasitosis.

Transcriptome profiling of grass carp (Ctenopharyngodon idellus) infected with Aeromonas hydrophila

Aeromonas hydrophila is the causative pathogen of intestinal hemorrhage which has caused great economic loss in grass carp aquaculture. In order to understand the immunological response of grass carp to infection by A. hydrophila, the transcriptomic profiles of the spleens from infected and non-infected grass carp groups were obtained using HiSeq™ 2500 (Illumina). An average of 63 million clean reads per library was obtained, and approximately 80% of these genes were successfully mapped to the reference genome. A total of 1591 up-regulated and 530 down-regulated genes were identified. Eight immune-related categories involving 105 differently expressed genes were scrutinized. 16 of the differently expressed genes involving immune response were further validated by qRT-PCR. Our results provide valuable information for further analysis of the mechanisms of grass carp defense against A. hydrophila invasion.

Transcriptome and analysis on the complement and coagulation cascades pathway of large yellow croaker (Larimichthys crocea) to ciliate ectoparasite Cryptocaryon irritans infection

Large yellow croaker (Larimichthys crocea) is one of the most valuable marine fish in southern China. Given to the rapid development of aquaculture industry, the L. crocea was subjected to ciliate ectoparasite Cryptocaryon irritans. It therefore is indispensable and urgent to understand the mechanism of L. crocea host defense against C. irritans infection. In the present study, the extensively analysis at the transcriptome level for Cryptocaryoniasis in L. crocea was carried out. These results showed that 15,826,911, 16,462,921, and 15,625,433 paired-end clean reads were obtained from three cDNA libraries (A: 0 theronts/fish, B: 12,000 theronts/fish, and C: 24,000 theronts/fish) of the L. crocea immune-related tissues by Illumina paired-end sequencing technology. Totally, 30,509 unique transcript fragments (unigenes) were assembled, with an average length of 1715 bp. In B/A, C/A, and C/B pairwise comparison, 972, 900, and 1126 genes showed differential expression respectively. Differently expressed immune-related genes (DEIGs) were scrutinized, in B/A pairwise comparison, 48 genes showed differential expression, including 26 up-regulated genes and 22 down-regulated genes in B; in C/A pairwise comparison, there were 39 DEIGs, including 7 up-regulated genes and 32 down-regulated genes in C; in C/B pairwise comparison, 40 genes showed differential expression, including 11 up-regulated genes and 29 down-regulated genes in C. There were 16 DEIGs enriched KEGG pathways, in which the complement and coagulation cascades pathway was the top most DEIGs enriched pathway (B:A = 42; C:A = 28; C:B = 42). The coagulation and fibrinolytic system was in a highly active state after infected by C. irritans with non-lethal concentration; the alternative complement pathway may play an important role in the early stages of C. irritans infection. These results demonstrated that low-concentration infection can significantly induce the immunological response in fishes, however, when fishes were in fatal conditions, the immunity was suppressed.

Differential modulation of host genes in the kidney of brown trout Salmo trutta during sporogenesis of Tetracapsuloides bryosalmonae (Myxozoa)

Tetracapsuloides bryosalmonae (Myxozoa) is the causative agent of proliferative kidney disease in various species of salmonids in Europe and North America. In Europe, spores of T. bryosalmonae develop in the kidney of infected brown trout Salmo trutta and are released via urine to infect the freshwater bryozoan Fredericella sultana. The transcriptomes of kidneys of infected and non-infected brown trout were compared by suppressive subtractive hybridization. Differential screening and a subsequent NCBI BLAST analysis of expressed sequence tags revealed 21 transcripts with functions that included cell stress and cell growth, ribonucleoprotein, signal transduction, ion transporter, immune response, hemoglobin and calcium metabolisms. Quantitative real time PCR was used to verify the presence of these selected transcripts in brown trout kidney at sporogonic stages of T. bryosalmonae development. Expression of cold-inducible RNA-binding protein, cyclin-dependent kinase inhibitor 2A, prothymosin alpha, transforming protein RhoA, immunoglobulin light chain and major histocompatibility complex class I were up-regulated significantly in infected brown trout. Expression of both the hemoglobin subunit beta and stanniocalcin precursor were down-regulated significantly in infected brown trout. This study suggests that cell stress and cell growth processes, signal transduction activities, erythropoiesis and calcium homeostasis of the host are modulated during sporogonic stages of parasite development, which may support the sporogenesis of T. bryosalmonae in the kidney of brown trout.

Interleukin gene expression is strongly modulated at the local level in a fish-parasite model

The goal of this work was to identify interleukin (IL)-related genes in the gilthead sea bream (GSB) (Sparus aurata L.) and how they are modulated by the parasite Enteromyxum leei, a myxozoan that causes severe enteritis with a strong inflammatory response. A Blast-X search of our transcriptomic GSB database (www.nutrigroup-iats.org/seabreamdb) identified 16 new sequences encompassing seven ILs (IL-7, IL-8, IL-10, IL-12β, IL-15, IL-18, and IL-34), the interleukin enhancer-binding factor 2 (ILF2), and eight IL receptors (IL-R); IL-R1, IL-6RA, IL-6RB, IL-8RA, IL-10RA, IL-10RB, IL-18R1, and IL-22R. Except for ILF2, their expression, plus that of IL-1β, IL-1R2, IL-6, and TNF-α (from public repositories), were analysed by 96-well PCR array of samples of blood, spleen, head kidney, and intestine of GSB that were anally intubated with E. leei (recipient group, RCPT). Only the expression profile of the intestine of RCPT fish showed significant difference as compared to samples from PBS-inoculated fish. At 17 days post intubation (dpi), the expression of key pro-inflammatory ILs, such as IL-8, IL-8R, IL-12β, and TNFα was significantly up-regulated, whereas at 64 dpi, anti-inflammatory IL expression (IL-6, IL-6RB, IL-7, IL-10, IL-10RA, and IL-15) was predominant. These results indicate a modification of the IL expression at late times post infection, probably to protect the fish intestine from the parasite and damage inflicted by an excessive inflammatory response. Furthermore, the response is mainly mediated at the local level as no significant changes were detected in blood, spleen and head kidney.

Expression profile of immune-related genes in Lates calcarifer infected by Cryptocaryon irritans

Cryptocaryon irritans causes Cyptocaryonosis or white spot disease in a wide range of marine fish including Lates calcarifer (Asian seabass). However, the immune response of this fish to the parasite is still poorly understood. In this study, quantitative polymerase chain reaction (qPCR) was performed to assess the expression profile of immune-related genes in L. calcarifer infected by C. irritans. A total of 21 immune-related genes encoding various functions in the fish immune system were utilized for the qPCR analysis. The experiment was initiated with the infection of juvenile fish by exposure to theronts from 200 C. irritans cysts, and non-infected juvenile fish were used as controls. Spleen, liver, gills and kidney tissues were harvested at three days post-infection from control and infected fish. In addition, organs were also harvested on day-10 post-infection from fish that had been allowed to recover from day-4 up to day-10 post-infection. L. calcarifer exhibited pathological changes on day-3 post-infection with the characteristic presence of white spots on the entire fish body, excessive mucus production and formation of a flap over the fish eye. High quality total RNA was extracted from all tissues and qPCR was performed. The qPCR analysis on the cohort of 21 immune-related genes of the various organs harvested on day-3 post-infection demonstrated that most genes were induced significantly (p < 0.05) in all tissues, particularly liver (11/21 genes) and kidney (11/21). The expression profile demonstrated that induction of the MHC Class IIα gene was the highest compared to the other genes followed by serum amyloid A, CC chemokine and hepcidin-2 precursor genes. In fish that were allowed to recover from the C. irritans infection (10 days post-infection), expression of the immune-related genes was down-regulated to levels similar to the control fish. These results provide insights into the interaction between C. irritans and L. calcarifer and suggest that the innate immune system plays an important role in early defence against parasite infection allowing the fish to eventually recover from the infection.