Gene expression profile changes in rat dorsal horn after sciatic nerve injury

Abstract

Objective: This study aims to investigate gene expression changes in rat dorsal horns after sciatic nerve injury (SNI).

Methods: The GSE18803 microarray data collected from young and adult rats were downloaded from GEO. After preprocessing, differentially expressed genes (DEGs) between SNI and sham-operated groups were selected using Limma package, in young and adult group, respectively, followed by Venn analysis. Then, enrichment analyses were performed for these DEGs using DAVID. The STRING database was used to identify protein–protein interactions (PPIs) among these DEGs, and the module network was further extracted using plugin ClusterONE. Finally, protein domain enrichment analysis of DEGs in each module was performed using InterPro database.

Results: Totally, 210 and 50 DEGs were identified in adult and young group, respectively. Among them, 160 were specific in adult group (e.g. CCL2, NF-κB1 and RAC2); 9 were specific in young group (e.g. ILF3 and LYVE1); and 41 were common in both two groups (e.g. FCER1G, C1QA, C1QB and C1QC). The up-regulated DEGs were mostly enriched in immune response-related biological processes, as well as 15 immune- and inflammation-related pathways. Then, two modules were identified in PPI network. CCL2 and NF-κB1 had high connectivity degrees in module 1, and RAC2, FCER1G and CD68 in module 2.

Conclusion: CCL2, NF-κB1, RAC2, FCER1G and C1Q may contribute to the generation of neuropathic pain after SNI via immune and defense pathways. Among the five genes, the first three are specific in adult population, while the latter two are age-independent. They all might function through involvement of these immune or inflammatory pathways.

Keywords:

• Sciatic nerve injury
• differentially expressed genes
• enrichment analysis
• protein–protein interaction
• module analysis
• adult
• age-independent

Acknowledgement

We thank Costigan very much, the author who deposited their microarray dataset, GSE18803, into the public GEO database.