Multiple sclerosis (MS) is the main reason for disability caused by autoimmunity. However, effective biomarkers for MS have not been identified. This study explored the molecular mechanisms and potential therapeutic targets for MS occurrence and progression using bioinformatics and next generation sequencing (NGS) data analysis. NGS data GSE138614, including patients with MS and 25 normal control samples, were obtained from Gene Expression Omnibus (GEO) database Differentially expressed genes (DEGs) were filtered and subjected to gene ontology (GO) and pathway enrichment analyses. Protein–protein interaction (PPI) network and module analyses were performed based on the DEGs. MiRNA-hub gene regulatory network and TF-hub gene regulatory network analysis were built by Cytoscape to predict the underlying microRNAs (miRNAs) and transcription factors (TFs) associated with hub genes. We also validated the identified hub genes via receiver operating characteristic (ROC) curve analysis. A total of 959 DEGs (479 up regulated genes and 480 down regulated genes) were detected. The GO terms and pathways of DEGs involve immune system process, developmental process, immune system and regulation of cholesterol biosynthesis by SREBP (SREBF). The network analysis revealed that hub genes include LCK, PYHIN1, SLAMF1, DOK2, TAB2, CFTR, RHOB, LMNA, EGLN3 and ERBB3 might be involved in the development of MS. The present investigation illustrates a characteristic gene profile in MS, which might contribute to the interpretation of the progression of MS and provide novel biomarkers and therapeutic targets for MS.