Transcriptional profiling and multiple machine learning methods identified novel microglia biomarkers for Experimental Autoimmune Encephalomyelitis

Background: Microglia plays complex and crucial roles in multiple sclerosis (MS). This study aimed to explore the biological significance of microglia-associated genes in experimental autoimmune encephalomyelitis (EAE) . Methods: Differentially expressed genes (DEGs) were screened with six machine learning (ML) methods, which were also utilized to validate the microglia-associated DEGs in three public databases. ceRNA and Protein–protein interaction (PPI) network analyses were utilized to identify the interaction of the 6 novel biomarkers with other molecules. Then, CIBERSORT and single-sample gene set enrichment analysis (ssGSEA) were employed to quantify the relative abundance of each immune cell infiltration, respectively. qRT-PCR was performed to test the expression of key DEGs in murine models. Results: A total of 247 DEmRNA, 499 DElncRNAs and 269 DEcircRNAs were identified. With screening strategy of five ML algorithms, 6 DEmRNAs were obtained including NGP, HIST1H2BJ, PBLD1, MBLN3, CD180 and F10. Then the 6 DEmRNAs were used as a multigene signature to construct models to differentiate EAE from normal microglia, and AUC value for each model was greater than 0.8. The diagnostic value of these 6 DEmRNAs were identified and further verified by qRT-PCR. Then, differential expression for five out of these 6 DEmRNAs, namely NGP, HIST1H2BJ, PBLD1, MBLN3, and F10 were confirmed. Using PPI analysis, DEmRNAs frequently interacting with transcription factors (TFs), potential drugs and RBPs were identified. With immune cell infiltration analyses, we found EAE microglia presented high levels of immune infiltration, especially Nature Killer (NK) cells and CD8+ T cells. We also reported circRNA (circRNA_00638) was predicted to bind to 76 RBPs. Conclusions: We identified and validated 6 novel microglia related genes and developed a multigene signature with ML methods to confirm their ability to accurately diagnose and characterize biological alterations in EAE microglia. The six key DEmRNAs might also be latent targets for immunoregulatory therapy. To explore the transcriptional profiling as well as the interaction of key DEGs and other molecules, including circRNA, lncRNA, RBPs, immune cells of EAE microglia, we established EAE mice model and extracted brain microglia to compare with healthy cells. We then performed gene expression profiling analysis using data obtained from RNA-seq of 8 samples of EAE microglia and 8 samples of PBS control microglia