Abstract: Immune changes and inflammatory responses have been identified as central events in the pathological process of spinal cord injury. They can greatly affect nerve regeneration and functional recovery. However, there is still limited understanding of the peripheral immune inflammatory response in spinal cord injury. In this study, we obtained microRNA expression profiles from the peripheral blood of patients with spinal cord injury using high-throughput sequencing. We also obtained the mRNA expression profile of spinal cord injury patients from the Gene Expression Omnibus (GEO) database (GSE151371). We identified 54 differentially expressed microRNAs and 1656 differentially expressed genes using bioinformatics approaches. Functional enrichment analysis revealed that various common immune and inflammation-related signaling pathways, such as neutrophil extracellular trap formation pathway, T cell receptor signaling pathway, and nuclear factor-κB signal pathway, were abnormally activated or inhibited in spinal cord injury patient samples. We applied an integrated strategy that combines weighted gene co-expression network analysis, LASSO logistic regression, and SVM-RFE algorithm and identified three biomarkers associated with spinal cord injury: ANO10, BST1, and ZFP36L2. We verified the expression levels and diagnostic performance of these three genes in the original training dataset and clinical samples through the receiver operating characteristic curve. Quantitative polymerase chain reaction results showed that ANO10 and BST1 mRNA levels were increased and ZFP36L2 mRNA was decreased in the peripheral blood of spinal cord injury patients. We also constructed a small RNA-mRNA interaction network using Cytoscape. Additionally, we evaluated the proportion of 22 types of immune cells in the peripheral blood of spinal cord injury patients using the CIBERSORT tool. The proportions of naïve B cells, plasma cells, monocytes, and neutrophils were increased while the proportions of memory B cells, CD8+ T cells, resting natural killer cells, resting dendritic cells, and eosinophils were markedly decreased in spinal cord injury patients increased compared with healthy subjects, and ANO10, BST1 and ZFP26L2 were closely related to the proportion of certain immune cell types. The findings from this study provide new directions for the development of treatment strategies related to immune inflammation in spinal cord injury and suggest that ANO10, BST1, and ZFP36L2 are potential biomarkers for spinal cord injury. The study was registered in the Chinese Clinical Trial Registry (registration No. ChiCTR2200066985, December 12, 2022).

Key words: bioinformatics analysis, biomarker, CIBERSORT, GEO dataset, LASSO, miRNA-mRNA network, RNA sequencing, spinal cord injury, SVM-RFE, weighted gene co-expression network analysis