Identifying gene expression profile of spinal cord injury in rat by bioinformatics strategy

Spinal cord injury (SCI) leads to the loss of sensory, motor, and autonomic function. We aimed to identify the therapeutic targets of-SCI by bioinformatics analysis. The gene expression profile of GSE20907 was downloaded from gene expression omnibus database. By comparing gene expression profiles with control samples, we screened out several differentially expressed genes (DEGs) in 3 days, 2 weeks and 1 month post-SCI. The pathway enrichment and protein-protein interaction (PPI) network analysis for the identified DEGs were performed. Then, transcription factors and microRNAs for DEGs were predicted. We found that up-regulated DEGs mainly participated in cell cycle, oxidative phosphorylation and immune-related pathways; while down-regulated DEGs were mainly involved in oxidative phosphorylation and central nervous system disease signaling pathways. In the constructed PPI network, Bub1, Vascular endothelial growth factor, Topoisomerase II alpha (TOP2a) and Cdc20 showed better correspondence with cell cycle, repair system and nerve system. Furthermore, the up-regulated genes (Arpc1b, CD74 and Brd2) significantly mapped to the target genes of transcription factors. The down-regulated genes of 3 days post-injury and the up-regulated genes of 2 weeks post-injury were significantly enriched as the target genes of microRNAs (miR-129 and miR-124). In conclusion, our results may provide guidelines to discuss the collaboration of PPI network in carcinogenesis of SCI.