Transcriptome analysis of the Larimichthys crocea liver in response to Cryptocaryon irritans

The large yellow croaker (Larimichthys crocea) is an economically important marine fish cultured in China and East Asian countries and is facing a serious threat from Cryptocaryon irritans, which is a protozoan ectoparasite that infects most reared marine fish species. To understand the molecular immune mechanisms underlying the response to C. irritans, we first performed a comparative gene transcription analysis using livers from C. irritans-immunized L. croceas and from a control group through RNA-Seq technology. After the removal of low-quality sequences and assembly, 51360 contigs were obtained, with an average length of 1066.93 bp. Further, a blast analysis indicates that 30747 contigs can be annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was used to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. Moreover, 14470 genes were found in 303 KEGG pathways. We used RSEM and EdgeR to determine that 3841 genes were significantly differentially expressed (FDR < 0.001), including 2129 up-regulated genes and 1712 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes revealed major immune-related pathways, including the toll-like receptor, complement and coagulation cascades, and chemokine signaling pathways. In addition, 28748 potential simple sequence repeats (SSRs) were detected from 12776 transcripts, and 62992 candidate single nucleotide polymorphisms (SNPs) were identified in the L. croceas liver transcriptome. This study characterized a gene expression pattern for normal and C. irritans-immunized L. croceas for the first time and not only sheds new light on the molecular mechanisms underlying the host-C. irritans interaction but also facilitates future studies on L. croceas gene expression and functional genomics.