Identification of hub genes and key pathways in spinal cord injury via bioinformatics analysis

Abstract

This study aimed to explore the hub genes and related key pathways in Spinal Cord Injury (SCI) based on the bioinformatics analysis. Two microarray datasets (GSE45006, GSE45550) were obtained from the GEO database and were merged and batch-corrected. The differentially expressed genes (DEGs) in SCI were explored with the Limma, and the weighted gene co-expression network analysis (WGCNA) was conducted to explore the module genes. Functional enrichment analysis and Gene set variation analysis (GSVA) were used to investigate the biological functions and key pathways of the key genes related to SCI. Then the proteinprotein interaction (PPI) network was generated using the STING online tool, and the hub genes in SCI were identified. Receiver operating characteristic (ROC) curves were applied to assess the diagnostic value of the selected hub genes. We identified 554 DEGs in SCI, and 1236 key genes in SCI were selected via WGCNA. Totally 111 key genes related to SCI were discovered. Furthermore, the functional enrichment analysis showed that these key mRNAs were primarily enriched in the extracellular matrix (ECM)-related pathways and processes associated with wound healing and cell growth. The PPI network further filtered six hub genes (Cd44, Timp1, Loxl1, Col6a1, Col3a1, Col5a1) ranked by the degree, and the diagnostic value of the six hub genes was confirmed by the ROC curves. Six hub genes including Cd44, Timp1, Loxl1, Col6a1, Col3a1, and Col5a1 were identified in SCI, with differential expression and excellent diagnostic value, which might provide insight into the targeted therapy of SCI.