Comprehensive analysis of differential expression profiles via transcriptome sequencing in SH-SY5Y cells infected with CVA16

Abstract

Coxsackievirus A16 (CV-A16) is one of the viruses that is most frequently associated with hand-foot-and-mouth disease (HFMD). Previous studies have shown that CV-A16 infections are mostly self-limiting, but in recent years, it has been gradually found that CV-A16 infections can also induce neurological complications and eventually cause death in children with HFMD. Moreover, no curative drugs or preventative vaccines have been developed for CV-A16 infection. Therefore, it is particularly important to investigate the mechanism of CVA16 infection-induced neuropathy. In the current study, transcriptome sequencing technology was used to identify changes in the transcriptome of SH-SY5Y cells infected with CVA16, which might hide the mechanism of CV-A16-induced neuropathology. The transcriptome profiling showed that 82,406,974, 108,652,260 and 97,753,565 clean reads were obtained in the Control, CV-A16-12 h and CV-A16-24 h groups, respectively. And it was further detected that a total of 136 and 161 differentially expressed genes in CV-A16-12 h and CV-A16-24 h groups, respectively, when compared with Control group. Then, to explore the mechanism of CV-A16 infection, we focused on the common differentially expressed genes at different time points of CV-A16 infection and found that there were 34 differentially expressed genes based on which clustering analysis and functional category enrichment analysis were performed. The results indicated that changes in oxidation levels were particularly evident in the GO term analysis, while only the "Gonadotropin-releasing hormone receptor pathway"was enriched in the KEGG pathway analysis, which might be closely related to the neurotoxicity caused by CV-A16 infection. Meanwhile, the ID2 closely related to nervous system has been demonstrated to be increased during CV-A16 infection. Additionally, the data on differentially expressed non-protein-coding genes of different types within the transcriptome sequencing results were analyzed, and it was speculated that these dysregulated non-protein-coding genes played a pivotal role in CV-A16 infection. Ultimately, qRT-PCR was utilized to validate the transcriptome sequencing findings, and the results of qRT-PCR were in agreement with the transcriptome sequencing data. In conclusion, transcriptome profiling was carried out to analyze response of SH-SY5Y cells to CVA16 infection. And our findings provide important information to elucidate the possible molecular mechanisms which were linked to the neuropathogenesis of CV-A16 infection.