Background: Insulin resistance (IR) is a major contributing factor to the pathogenesis of metabolic syndrome and type 2 diabetes mellitus (T2D). Adipocyte metabolism is known to play a crucial role in IR. Therefore, the aims of this study were to identify metabolism-related proteins that could be used as potential biomarkers of IR and to investigate the role of N6-methyladenosine (m6A) modification in the pathogenesis of this condition.
Methods: RNA-seq data on human adipose tissue were retrieved from the Gene Expression Omnibus database. The differentially expressed genes of metabolism-related proteins (MP-DEGs) were screened using protein annotation databases. Biological function and pathway annotations of the MP-DEGs were performed through Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses. Key MP-DEGs were screened, and a protein-protein interaction (PPI) network was constructed using STRING, Cytoscape, MCODE, and CytoHubba. LASSO regression analysis was used to select primary hub genes, and their clinical performance was assessed using receiver operating characteristic (ROC) curves. The expression of key MP-DEGs and their relationship with m6A modification were further verified in adipose tissue samples collected from healthy individuals and patients with IR.
Results: In total, 69 MP-DEGs were screened and annotated to be enriched in pathways related to hormone metabolism, low-density lipoprotein particle and carboxylic acid transmembrane transporter activity, insulin signaling, and AMPK signaling. The MP-DEG PPI network comprised 69 nodes and 72 edges, from which 10 hub genes (FASN, GCK, FGR, FBP1, GYS2, PNPLA3, MOGAT1, SLC27A2, PNPLA3, and ELOVL6) were identified. FASN was chosen as the key gene because it had the highest maximal clique centrality (MCC) score. GCK, FBP1, and FGR were selected as primary genes by LASSO analysis. According to the ROC curves, GCK, FBP1, FGR, and FASN could be used as potential biomarkers to detect IR with good sensitivity and accuracy (AUC = 0.80, 95% CI: 0.67-0.94; AUC = 0.86, 95% CI: 0.74-0.94; AUC = 0.83, 95% CI: 0.64-0.92; AUC = 0.78, 95% CI: 0.64-0.92). The expression of FASN, GCK, FBP1, and FGR was significantly correlated with that of IGF2BP3, FTO, EIF3A, WTAP, METTL16, and LRPPRC (p < 0.05). In validation clinical samples, the FASN was moderately effective for detecting IR (AUC = 0.78, 95% CI: 0.69-0.80), and its expression was positively correlated with the methylation levels of FASN (r = 0.359, p = 0.001).
Conclusion: Metabolism-related proteins play critical roles in IR. Moreover, FASN and GCK are potential biomarkers of IR and may be involved in the development of T2D via their m6A modification. These findings offer reliable biomarkers for the early detection of T2D and promising therapeutic targets.
Keywords: bioinformatics analysis; insulin resistance; m6A modification; type 2 diabetes.