Diabetes, Metabolic Syndrome and Obesity (Sep 2018)
Identification of key gene pathways and coexpression networks of islets in human type 2 diabetes
Abstract
Lu Li,1,* Zongfu Pan,2,* Si Yang,1 Wenya Shan,1 Yanyan Yang1 1Department of Pharmacy, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, Zhejiang, People’s Republic of China; 2Department of Pharmacy, Zhejiang Cancer Hospital, Hangzhou, Zhejiang, People’s Republic of China *These authors contributed equally to this work Purpose: The number of people with type 2 diabetes (T2D) is growing rapidly worldwide. Islet β-cell dysfunction and failure are the main causes of T2D pathological processes. The aim of this study was to elucidate the underlying pathways and coexpression networks in T2D islets. Materials and methods: We analyzed the differentially expressed genes (DEGs) in the data set GSE41762, which contained 57 nondiabetic and 20 diabetic samples, and developed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Protein–protein interaction (PPI) network, the modules from the PPI network, and the gene annotation enrichment of modules were analyzed as well. Moreover, a weighted correlation network analysis (WGCNA) was applied to screen critical gene modules and coexpression networks and explore the biological significance. Results: We filtered 957 DEGs in T2D islets. Then GO and KEGG analyses identified that key pathways like inflammatory response, type B pancreatic cell differentiation, and calcium ion-dependent exocytosis were involved in human T2D. Three significant modules were filtered from the PPI network. Ribosome biogenesis, extrinsic apoptotic signaling pathway, and membrane depolarization during action potential were associated with the modules, respectively. Furthermore, coexpression network analysis by WGCNA identified 13 distinct gene modules of T2D islets and revealed four modules, which were strongly correlated with T2D and T2D biomarker hemoglobin A1c (HbA1c). Functional annotation showed that these modules mainly enriched KEGG pathways such as NF-kappa B signaling pathway, tumor necrosis factor signaling pathway, cyclic adenosine monophosphate signaling pathway, and peroxisome proliferators-activated receptor signaling pathway. Conclusion: The results provide potential gene pathways and underlying molecular mechanisms for the prevention, diagnosis, and treatment of T2D. Keywords: type 2 diabetes, islet β cell, bioinformatics analysis, differentially expressed genes, WGCNA