Bioinformatics analysis of potential key genes and mechanisms in type 2 diabetes mellitus

Type 2 diabetes mellitus (T2DM) is etiologically related to metabolic disorder. The aim of our study was to screen out candidate genes of T2DM and to elucidate the underlying molecular mechanisms by bioinformatics methods. Expression profiling by high throughput sequencing data of GSE154126 was downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) between T2DM and normal control were identified. And then, functional enrichment analyses of gene ontology (GO) and REACTOME pathway analysis was performed. Protein-protein interaction (PPI) network and module analyses were performed based on the DEGs. Additionally, potential miRNAs of hub genes were predicted by miRNet database. Transcription factors (TFs) of hub genes were detected by NetworkAnalyst database. Further, validations were performed by receiver operating characteristic curve (ROC) analysis and real-time polymerase chain reaction (RT-PCR). In total, 925 DEGs were identified in T2DM, including 447 up regulated genes and 478 down-regulated genes. Functional enrichment analysis results showed that up regulated DEGs were significantly enriched in defense response, neutrophil degranulation, cell adhesion and extracellular matrix organization. The top 10 hub genes, JUN, VCAM1, RELA, U2AF2, ADRB2, FN1, CDK1, TK1, A2M and ACTA2 were identified from the PPI network, modules, miRNA-hub gene regulatory network and TF-hub gene regulatory network. Furthermore, ROC analysis and RT-PCR revealed that JUN, VCAM1, RELA, U2AF2, ADRB2, FN1, CDK1, TK1, A2M and ACTA2 might serve as biomarkers in T2DM. Bioinformatics analysis is a useful tool to explore the molecular mechanism and pathogenesis of T2DM. The identified hub genes may participate in the onset and advancement of T2DM and serve as therapeutic targets.