Abstract:Background: Spinal cord injury (SCI) is one of the most devastating diseases with a high incidence rate around the world. SCI-related neuropathic pain (NeP) is a common complication, whereas its pathomechanism is still unclear. The purpose of this study is to identify key genes and cellular components for SCI-related NeP by an integrated transcriptome bioinformatics analysis. Methods: The gene expression profile of 25 peripheral blood samples from chronic phase SCI patients (E-GEOD-69901) and 337 normal peripheral blood samples were downloaded from ArrayExpress and Genotype-Tissue Expression Portal (GTEx), respectively. A total of 3,368 normal peripheral blood mononuclear cells (PBMC) were download from Sequence Read Archive (SRA713577). Non-parametric tests were used to evaluate the association between all of differential expression genes (DEGs) and SCI-related NeP. CellPhoneDB algorithm was performed to identify the ligand-receptor interactions and their cellular localization among single PBMCs. Transcription factor (TF) enrichment analysis and Gene Set Variation Analysis (GSVA) were used to identify the potential upstream regulatory TFs and downstream signaling pathways, respectively. Co-expression analysis among significantly enriched TFs, key cellular communication genes and differentially expressed signaling pathways were performed to identify key genes and cellular components for SCI-related NeP. Results: A total of 2,314 genes were identified as DEGs between the experimental and the control group. Five proteins (ADRB2, LGALS9, PECAM1, HAVCR2, LRP1) were identified in the overlap of proteins in the significant ligand-receptor interactions of PBMCs and protein-protein interaction (PPI) network based on the DEGs. Only HAVCR2 was significantly associated with NeP (P = 0.005). Besides, the co-expression analysis revealed that TF YY1 had significantly co-expression pattern with cellular communication receptor HAVCR2 (R = -0.54, P < 0.001) in NK cells while HAVCR2 was also co-expressed with mTOR signaling pathway (R = 0.57, P < 0.001). The results of RT-qPCR and external dataset validation supported the signaling axis with the most significant co-expression patterns. Conclusion: In peripheral blood of chronic SCI, HAVCR2 might act as a key receptor on the surface of NK cells and interact with ligand LGALS9 secreted by CD14+ monocytes, inhibiting NK cells through mTOR signaling pathway and ultimately predicting the occurrence of SCI-related NeP. This hypothetical signaling axis may provide prognostic biomarkers and therapeutic targets for SCI-related NeP.

Functional investigation and two-sample Mendelian randomization study of neuropathic pain hub genes obtained by WGCNA analysis

Comprehensive Analysis of Differentially Expressed Genes and Key Pathways in Neuropathic Pain by Spinal Nerve Ligation

Identification of the Key Genes Associated with Neuropathic Pain.

Bioinformatic Analysis of Neuroimmune Mechanism of Neuropathic Pain

Identification of the Hub Genes Related to Nerve Injury-Induced Neuropathic Pain.

Jun, Gal, Cd74, and C1qb as potential indicator for neuropathic pain

Bidirectional two-sample Mendelian randomization analysis reveals a causal effect of interleukin-18 levels on postherpetic neuralgia risk

Identification of Neuropathic Pain-Associated Genes and Pathways Via Random Walk with Restart Algorithm

The Integrated Transcriptome Bioinformatics Analysis Identifies Key Genes and Cellular Components for Spinal Cord Injury-Related Neuropathic Pain

Temporal Changes of Spinal Transcriptomic Profiles in Mice With Spinal Nerve Ligation.

Bioinformatics Analysis Identifies P53 As a Candidate Prognostic Biomarker for Neuropathic Pain

Gene Expression Profiles Reveal Key Pathways and Genes Associated with Neuropathic Pain in Patients with Spinal Cord Injury

Identification and Validation of Core Genes Associated with Intracranial Aneurysms Through Bioinformatics Analysis and Mendelian Randomization

Gene expression profiles reveal key pathways and genes associated with neuropathic pain in patients with spinal cord injury.

Systemic Inflammatory Regulators Are Associated with Two Common Types of Neuropathic Pain: A Bidirectional Mendelian Randomization Study

Shared Genetic Regulatory Networks Contribute to Neuropathic and Inflammatory Pain: Multi-Omics Systems Analysis

Co‐expression gene modules involved in cisplatin‐induced peripheral neuropathy according to sensitivity, status, and severity

Effect of Surgical Damage to Spinal Nerve on Dorsal-Root Ganglion Gene Expressions: Comprehensive Analysis of Differentially Expressed Genes.

Bioinformatic Analysis of Potential Biomarkers for Spinal Cord-injured Patients with Intractable Neuropathic Pain.

Predicting pain genes: multi-modal data integration using probabilistic classifiers and interaction networks

Genetics of causal relationships between circulating inflammatory proteins and postherpetic neuralgia: a bidirectional Mendelian randomization study