Analysis of gene expression profiles in two spinal cord injury models

Abstract Objectives To analyze the changes of gene expression at different timepoints after spinal cord injury (SCI) with tenth segment thoracic injury. Methods Two SCI models, the complete paraplegia (H) and Allen’s strike (D) methods were applied to induce SCI in rats, and transcriptome sequencing was performed 1, 3, 7, 14, 56, and 70 days after SCI, respectively. Principal component analysis, differentially expressed gene analysis, and hierarchical clustering analysis were applied to analyze the differentially expressed genes (DEGs). Gene Ontology GO enrichment analysis, Kyoto Encyclopedia of Genes and Genomes enrichment analysis, and Gene Set Enrichment Analysis revealed the pathway of gene enrichment. Results There were 1,907, 3,120, 3,728, 978, 2,319, and 3,798 DEGs in the complete paraplegia group and 2,380, 878, 1,543, 6,040, 1,945, and 3,850 DEGs in the Allen’s strike method group and after SCI at 1, 3, 7, 14, 56, and 70 days, respectively. The transcriptome contours of D1, H1, D3, and H14 were clustered with C; the H56, D56, H70, and D70 transcriptome contours were similar and clustered together. H3, D7, and H7 were clustered together, and D14 was clustered separately. The transcriptome differences of the two SCI models were mainly concentrated during the first 2 weeks after SCI. The DEGs after SCI in the complete paraplegia group were more concentrated. Most of the early transcriptional regulation stabilized within 2 weeks after injury. Conclusions There were DEGs between the two SCI models. Through the gene changes and pathway enrichment of the entire time period after SCI, the molecular mechanism of SCI repair was revealed in depth, which provided a reference for SCI treatment in the future.