Identification and expression analysis of lncRNAs reveal the immune mechanism of visceral white-nodules disease resistance in large yellow croaker

Long noncoding RNAs (lncRNAs) have several known functions in fish growth processes and signal transduction, but their possible roles in response to bacterial diseases remain largely unresolved. In this study, we report a comprehensive cold-water bacterial disease-responsive lncRNAs expression profile for defense against visceral white-nodules disease in large yellow croaker. A total of 2,534 high-confidence lncRNAs were identified by a rigorous filtering pipeline as basic sequence set for comparative transcriptional analysis. In addition, a total of 10,200 lncRNA-mRNA pairs with high correlation coefficients were identified by expressions level correlation analysis, including non-redundant 381 DE lncRNAs and 2,590 differential expressed genes, and MSTRG_11084_1 and MSTRG_20402_1 were considered as core lncRNA. We further revealed the conserved and idiosyncratic features of the disease response process between the resistant strain (RS) and the technical control strain (TCS). Immune-related pathways were enriched in GO terms and KEGG pathway, among which Cytokine-cytokine receptor interaction, MAPK signaling pathway, and NF-kappa B signaling pathway may play a key role in VWND resistance in large yellow croaker. Protein-protein interaction network (PPI) analysis revealed that immune-related target genes such as il-10, met, acta2, myc, cav1 and ntrk1, as well as growth and metabolism-related target genes such as pik3r2, igf1, sc5d, hmgcr and lss were considered the main hub genes. This study represents the first characterization of lncRNAs involved in VWND resistance in large yellow croaker and provides new clues for elucidating the disease response mechanism of large yellow croaker.