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.

Single-cell RNA sequencing uncovers the cell type-dependent transcriptomic changes in the retrosplenial cortex after peripheral nerve injury

Single-cell RNA sequencing uncovers the excitatory/inhibitory synaptic unbalance in the retrosplenial cortex after peripheral nerve injury

Single-cell RNA Sequencing in the Context of Neuropathic Pain: Progress, Challenges, and Prospects.

Single-cell RNA-seq reveals distinct injury responses in different types of DRG sensory neurons

Dorsal Root Ganglion Transcriptome Analysis Following Peripheral Nerve Injury In Mice

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

Gene Transcript Alterations in the Spinal Cord, Anterior Cingulate Cortex, and Amygdala in Mice Following Peripheral Nerve Injury

Single-Nuclei RNA Sequencing Identifies Type C Low-Threshold Mechanoreceptors as Key Players in Paclitaxel-Induced Peripheral Neuropathy

Single-cell profiling of cellular changes in the somatic peripheral nerves following nerve injury

Major differences in transcriptional alterations in dorsal root ganglia between spinal cord injury and peripheral neuropathic pain models

RNA Sequencing of Peripheral Blood Revealed that the Neurotropic TRK Receptor Signaling Pathway Shows Apparent Correlation in Recovery Following Spinal Cord Injury at Small Cohort

Comparative transcriptional analysis of satellite glial cell injury response

Convergent and divergent transcriptional reprogramming of motor and sensory neurons underlying response to peripheral nerve injury

Single-cell analysis of spinal cord injury reveals functional heterogeneity of oligodendrocyte lineage cells

Schwann cells regulate sensory neuron gene expression before and after peripheral nerve injury

Temporal Changes in the Spinal Cord Transcriptome after Peripheral Nerve Injury

Single-cell analysis reveals region-heterogeneous responses in rhesus monkey spinal cord with complete injury

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

Global Analysis of Transcriptome in Dorsal Root Ganglia Following Peripheral Nerve Injury in Rats

The Temporal and Spatial Signature of Microglial Transcriptome in Neuropathic Pain

Temporal and spatial cellular and molecular pathological alterations with single-cell resolution in the adult spinal cord after injury