Identification of crucial genes associated with rat traumatic spinal cord injury

The aim of the present study was to investigate the key genes associated with traumatic spinal cord injuries (TSCI). The dataset GSE52763 was downloaded from the Gene Expression Omnibus, for which lumbar spinal cord samples were obtained from rats at 1 and 3 weeks following contusive spinal cord injury and 1 week subsequent to a sham laminectomy, and used to identify differentially expressed genes (DEGs). Functional enrichment analysis, co‑expression analysis and transcription factor (TF) identification were performed for DEGs common to the 1 and 3 week injury samples. In total, 234 upregulated and 51 downregulated DEGs were common to the 1 and 3 week injury samples. The upregulated DEGs were significantly enriched in Gene Ontology terms concerning immunity (e.g. Itgal and Ccl2) and certain pathways, including natural killer cell mediated cytotoxicity [e.g. Ras‑related C3 botulinum toxin substrate 2 (Rac2) and TYRO protein tyrosine kinase binding protein (Tyrobp)]. The downregulated DEGs were highly enriched in female gonad development [e.g. progesterone receptor (Pgr)], and the steroid biosynthesis pathway. A total of 139 genes had co‑expression associations and the majority of them were upregulated genes. The upregulated co‑expressed genes were predominantly enriched in biological regulation, including TGFB induced factor homeobox 1 (Tgif1) and Rac2. The downregulated co‑expressed genes were enriched in anatomical structure development (e.g. Dnm3). A total of 92 co‑expressed genes composed the protein‑protein interaction network. Additionally, 9 TFs (e.g. Pgr and Tgif1) were identified from the DEGs. It was hypothesized that the genes including Tgif1, Rac2, Tyrobp, and Pgr may be closely associated with TSCI.