Abstract
Diabetic neuropathy (DN) is one of the chronic complications of diabetes which can cause severe harm to patients. In order to determine the key genes and pathways related to the pathogenesis of DN, we downloaded the microarray data set GSE27382 from Gene Expression Omnibus (GEO) and adopted bioinformatics methods for comprehensive analysis, including functional enrichment, construction of PPI networks, central genes screening, TFs-target interaction analysis, and evaluation of immune infiltration characteristics. Finally, we examined quantitative real- time PCR (qPCR) to validate the expression of hub genes. A total of 318 differentially expressed genes (DEGs) were identified, among which 125 upregulated DEGs were enriched in the mitotic nuclear division, extracellular region, immunoglobulin receptor binding, and p53 signaling pathway, while 193 downregulated DEGs were enriched in ion transport, membrane, synapse, sodium channel activity, and retrograde endocannabinoid signaling. GSEA plots showed that condensed nuclear chromosome kinetochore were the most significant enriched gene set positively correlated with the DN group. Importantly, we identified five central genes (Birc5, Bub1, Cdk1, Ccnb2, and Ccnb1), and KEGG pathway analysis showed that the five hub genes were focused on progesterone-mediated oocyte maturation, cell cycle, and p53 signaling pathway. The proportion of immune cells from DN tissue and normal group showed significant individual differences. In DN samples, T cells CD4 memory resting and dendritic cells resting accounted for a higher proportion, and macrophage M2 accounted for a lower proportion. In addition, all five central genes showed consistent correlation with immune cell infiltration levels. qPCR showed the same expression trend of five central genes as in our analysis. Our research identified key genes related to differential genes and immune infiltration related to the pathogenesis of DN and provided new diagnostic and potential therapeutic targets for DN.