Transcriptome analysis reveals the temporal gene expression patterns in skin of large yellow croaker (Larimichthys crocea) in response to Cryptocaryon irritans infection.

Large yellow croaker (Larimichthys crocea) is one of the most important mariculture fish in China. In the past decades, cryptocaryonosis caused by Cryptocryon irritans has led to huge economic losses, posing great threat to the healthy and sustainable development of L. crocea mariculture industry. As the largest immunologically active mucosal organ in fish, skin provides the first defense line against external pathogens. To better understand the gene expression dynamics, the large yellow croakers were artificially infected with C. irritans and their skin tissues were collected at 0 h, 24 h, 48 h, 72 h and 96 h post infection. The total RNA in the skin tissues were extracted and the transcriptome were sequenced. After sequencing, a total of 1,131, 311, 140 million high quality RNA-seq reads were collected. A set of 215, 473, 968, 1055 differentially expressed genes were identified at 24 h, 48 h, 72 h and 96 h post infection respectively. Further analysis clustered these DEGs into six profiles and 75 hub genes for six profiles were identified. Among these hub genes, 18 immune related genes including TLR5, TOPK, NFKBIZ, MAPK14A were identified post C. irritans infection. Cytokine-cytokine receptor interaction was the only pathway that significantly enriched at four timepoints post infection. This study provides an in-depth understanding of skin transcriptome variance of large yellow croaker after C. irritans infection, which would be helpful for further understanding of the molecular mechanism of L. crocea in response to C. irritans infection.