Identification of key genes and pathways associated with spinal cord injury

Abstract

The present study was designed to identify key genes or significant signaling pathways associated with spinal cord injury (SCI), and to clarify the underlying molecular mechanisms of SCI. Data from the GSE45550 array were downloaded from the Gene Expression Omnibus database. A total of 6 control samples, 6 samples at 3 days post-SCI (SCI3d), 6 samples at 8 days post-SCI (SCI8d) and 6 samples at 14 days post-SCI (SCI14d) were included. The microarray data was preprocessed by the robust multi-array average algorithm. The differentially expressed genes (DEGs) were identified using the limma package. The overlapping DEGs among groups were analyzed using the Venny 2.0 online tool. Modules enriched by DEGs were selected by weighted gene co-expression network analysis. Gene Ontology annotation and the Kyoto Encyclopedia of Genes and Genomes pathways were identified for DEGs using the Database for Annotation, Visualization and Integrated Discovery. A total of 693 genes were obtained by combining the DEGs of the SCI3d, SCI8d and SCI14d groups. The pink module and green module with smaller P-values obtained from weighted gene co-expression network analysis module analyses of DEGs demonstrated a higher correlation with SCI. In addition, the peroxisome proliferator-activated receptor (PPAR) signaling pathway that the cluster of differentiation 36 (CD36) was significantly enriched in, was one of the significant pathways in the pink module. The p53 signaling pathway that Caspase-3 (Casp3) was significantly enriched in was one of the significant pathways in the green module. In conclusion, the results of the present study demonstrated that the PPAR and p53 signaling pathway may serve important roles in the progression of SCI. In addition, CD36 and Casp3 may be involved in the progression of SCI via the PPAR and p53 signaling pathways, respectively.