Characterization of Critical Functions of Long Non-Coding RNAs and mRNAs in Rhabdomyosarcoma Cells and Mouse Skeletal Muscle Infected by Enterovirus 71 Using RNA-Seq

Enterovirus 71 (EV71) is the main pathogen of severe hand-foot-mouth disease (HFMD). Long non-coding RNAs (lncRNAs) are recognized as pivotal factors during the pathogenesis of viral infection. However, the critical functions of lncRNAs in EV71-host interactions have not been characterized. Here, for the first time, we performed global transcriptome analysis of lncRNA and mRNA expression profiles in EV71-infected human rhabdomyosarcoma (RD) cells and skeletal muscle of mice using second-generation sequencing. In our study, a total of 3801 novel lncRNAs were identified. In addition, 23 lncRNAs and 372 mRNAs exhibited remarkable differences in expression levels between infected and uninfected RD cells, while 104 lncRNAs and 2647 mRNAs were differentially expressed in infected skeletal muscle from neonatal mice. Comprehensive bioinformatics analysis included target gene prediction, lncRNA-mRNA co-expression network construction, as well as gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis mainly focused on differentially-expressed genes (DEGs). Our results suggest that lncRNAs may participate in EV71 infection-induced pathogenesis through regulating immune responses, protein binding, cellular component biogenesis and metabolism. The present study provides novel insights into the functions of lncRNAs and the possible pathogenic mechanism following EV71 infection.