Uncovering the gene regulatory network of type 2 diabetes through multi-omic data integration.

Background: Type 2 diabetes (T2D) onset is a complex, organized biological process with multilevel regulation, and its physiopathological mechanisms are yet to be elucidated. This study aims to find out the key drivers and pathways involved in the pathogenesis of T2D through multi-omics analysis.

Methods: The datasets used in the experiments comprise three groups: (1) genomic (2) transcriptomic, and (3) epigenomic categories. Then, a series of bioinformatics technologies including Marker set enrichment analysis (MSEA), weighted key driver analysis (wKDA) was performed to identify key drivers. The hub genes were further verified by the Receiver Operator Characteristic (ROC) Curve analysis, proteomic analysis, and Real-time quantitative polymerase chain reaction (RT-qPCR). The multi-omics network was applied to the Pharmomics pipeline in Mergeomics to identify drug candidates for T2D treatment. Then, we used the drug-gene interaction network to conduct network pharmacological analysis. Besides, molecular docking was performed using AutoDock/Vina, a computational docking program.

Results: Module-gene interaction network was constructed using MSEA, which revealed a significant enrichment of immune-related activities and glucose metabolism. Top 10 key drivers (PSMB9, COL1A1, COL4A1, HLA-DQB1, COL3A1, IRF7, COL5A1, CD74, HLA-DQA1, and HLA-DRB1) were selected by wKDA analysis. Among these, COL5A1, IRF7, CD74, and HLA-DRB1 were verified to have the capability to diagnose T2D, and expression levels of PSMB9 and CD74 had significantly higher in T2D patients. We further predict the co-expression network and transcription factor (TF) binding specificity of the key driver. Besides, based on module interaction networks and key driver networks, 17 compounds are considered to possess T2D-control potential, such as sunitinib.

Conclusions: We identified signature genes, biomolecular processes, and pathways using multi-omics networks. Moreover, our computational network analysis revealed potential novel strategies for pharmacologic interventions of T2D.