Regulation of gene expression in rats with spinal cord injury based on microarray data

The present study aimed to investigate the molecular mechanisms of spinal cord injury (SCI) in rats. First, the differentially expressed genes (DGEs) were screened based on GSE45006 microarray data downloaded from Gene Expression Omnibus using the significant analysis of microarray (SAM) method. Screening was performed for DEGs which were negatively or possibly correlated with time and subsequently subjected to gene ontology (GO) functional annotation. Furthermore, pathway enrichment analysis using the Kyoto Encyclopedia of Genes and Genomes was also performed. In addition, a protein-protein interaction (PPI) network was constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins database. Finally, GeneCodis was used to seek transcription factors and microRNAs that are involved in the regulation of DEGs. A total of 806 DEGs were upregulated and 549 DEGs were downregulated in the rats with SCI. Cholesterol metabolism-associated genes (e.g. HMGCS1, FDFT1 and IDI1) were negatively correlated with time, while injury genes (e.g. SERPING1, C1S and RAB27A) were positively correlated with time after SCI. PCNA, MCM2, JUN and SNAP25 were the hub proteins of the PPI network. The transcription factors LEF1 and SP1 were observed to be associated with the regulation of two DEGs that were involved in the cholesterol-associated metabolism as well as injury responses. A number of microRNAs (e.g. miR210, miR-487b and miR-16) were observed to target cholesterol metabolism-associated DGEs. The hub genes PCNA, MCM2, JUN and SNAP25 presumably have critical roles in rats with SCI, and the transcription factors LEF1 and SP1 may be important for the regulation of cholesterol metabolism and injury responses following SCI.