Abstract
Objectives: As a high-cost neurological disability, spinal cord injury (SCI) can result in permanent paralysis and loss of sensation. To identify the temporal genes involved in the pathogenesis of SCI, we analysed the expression profile of GSE45006.
Methods: GSE45006 was downloaded from Gene Expression Omnibus, including 20 SCI samples (samples at 1 day, 3 days, 1 week, 2 weeks and 8 weeks after injury, four repetitions for each time point) and 4 normal samples. The Bayesian Estimation of Temporal Regulation (BETR) and randomForest packages were used to screen the temporal genes and the top 100 temporal genes, respectively. Then, the gplots package and Pearson correlation analysis separately were used to perform hot map analysis and expression pattern clustering for the top 100 temporal genes. Using the clusterProfiler package and TargetMine tool, their potential functions were analysed by enrichment analyses. Moreover, interaction relationships between these temporal genes and pathways were investigated by pathway-gene crosslinking networks.
Results: In total, 1907 temporal genes were identified. The top 100 temporal genes were obtained and divided into six clusters. Most of the gene functions were enriched in biological process categories. ARG1 and NOS3 in cluster 4 were enriched in biological process of arginine catabolic process. TGF beta 2, TGF beta 3, ALDH2 and ALDH3A2 were correlated with numerous pathways in the pathway-gene crosslinking network. Pathways related to TGF beta 2 and TGF beta 3 were connected to pathways related to ARG1 and NOS3 via ARG1.
Conclusion: Several temporal genes, including TGF beta 2, TGF beta 3, ALDH2, ALDH3A2, ARG1 and NOS3, might be involved in SCI.