Neutrophil extracellular traps (NETs) cause delayed wound closed up in type 2 diabetes mellitus (T2DM), but the specific regulatory mechanism of NETs-related genes (NETs-RGs) in T2DM is unclear.
We acquired GSE21321 and GSE15932 datasets from gene expression omnibus (GEO) database. First, differentially expressed genes (DEGs) between T2DM and control samples of GSE21321 dataset were sifted out by differential expression analysis. NETs scores were calculated for all samples in GSE21321 dataset, and key module genes associated with NETs scores were screened by constructing co-expression network. Then, DEGs and key module genes were intersected to yield intersection genes, and candidate genes were identified by constructing a protein protein interaction (PPI) network. Least absolute shrinkage and selection operator (LASSO) regression analysis was implemented on candidate genes to screen out diagnostic genes, and they were subjected to single sample gene set enrichment analysis (ssGSEA). Finally, immune characteristic analysis was carried out, and we constructed the gene-drug and transcription factor (TF)-miRNA-mRNA networks. Besides, we validated the expression of diagnostic genes by quantitative real-time polymerase chain reaction (qRT-PCR).
In total, 23 candidate genes were gained by PPI analysis. The 5 diagnostic genes, namely, inter-trypsin inhibitor heavy chain 3 (ITIH3), fibroblast growth factor 1 (FGF1), neuron cell adhesion molecule (NRCAM), advanced glycosylation end-product-specific receptor (AGER), and calcium voltage-gated channel subunit alpha1 C (CACNA1C), were identified via LASSO analysis, and they were involved in carboxylic acid transport, axonogenesis, etc. M2 Macrophage, Monocyte, Natural killer (NK) cell, and Myeloid dendritic cells (DC) were remarkably different between T2DM and control samples. Diagnostic genes had the strongest and the most significant positive correlation with B cells. The gene-drug network included CACNA1C-Isradipine, CACNA1C-Benidipine and other relationship pairs. Totally 76 nodes and 44 edges constituted the TF-miRNA-mRNA network, including signal transducer and activator of transcription 1(STAT1) -hsa-miR-3170-AGER, CCCTC-binding factor (CTCF)-hsa-miR-455-5p-CACNA1C, etc. Moreover, qRT-PCR suggested that the expression trends of FGF1 and AGER were in keeping with the results of bioinformatic analysis. FGF1 and AGER were markedly regulated downwards in the T2DM group.
Through bioinformatic analysis, we identified NETs-related diagnostic genes (ITIH3, FGF1, NRCAM, AGER, CACNA1C) in T2DM, and explored their mechanism of action from different aspects, providing new ideas for the studies related to diagnosis and treatment of T2DM.