Machine learning and bioinformatics identify and validate membrane protein biomarkers for type 2 diabetes mellitus

Tong Liu, Zhaoyuan Tang, Zulipikaer Aierken, Jie Zhang, Cheng Lei

Introduction Type 2 diabetes mellitus (T2DM) is among the most rapidly increasing metabolic disorders worldwide. Membrane proteins, integral components of biological membranes, are pivotal in insulin signal transduction and significantly contribute to the pathogenesis of T2DM. However, the systematic investigation of membrane proteins linked to T2DM remains insufficient. Methods The T2DM dataset and membrane protein-related genes were sourced from the GEO database and the Uniprot website, respectively. Bioinformatics methodologies, including Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and protein-protein interaction (PPI) network analysis, were employed to assess the differentially expressed membrane protein-related genes between the normal control group and the T2DM group. Subsequently, machine learning algorithms, including Gaussian Mixture Model (GMM), Random Forest (RF), and Support Vector Machine (SVM), were utilized to identify hub genes. Following this, a clinical diagnostic model was constructed, and the receiver operating characteristic (ROC) curve was plotted. Candidate genes were further examined in palmitic acid (PA)-induced NES2Y and HepG2 cell models and in high-fat diet/streptozotocin-induced T2DM mice. Finally, Functional effects were assessed by gene knockdown, reverse transcription quantitative PCR (RT-qPCR), western blotting (WB), immunofluorescence (IF), co-immunoprecipitation, histological examination, and lipid staining. Results A total of 1,856 DEGs were identified between T2DM and control samples, including 599 upregulated genes and 1,257 downregulated genes, among which 42 were membrane proteinrelated T2DM DEGs. Machine learning algorithms were then applied to pinpoint two target genes: ICAM1 and EZR, the latter of which encodes Ezrin, a membrane–cytoskeleton linker protein. The ROC curve analysis showed that the diagnostic model exhibited strong predictive capability, with an AUC value of 0.95. PA treatment increased lipid accumulation and EZR/ICAM1 mRNA and Ezrin/ICAM1 protein expression in the cell models. Ezrin co-immunoprecipitated with the insulin receptor. In PA-treated HepG2 cells, ICAM1 or EZR knockdown increased p85α and AKT phosphorylation and GLUT4 protein expression. Ezrin and ICAM1 were also elevated in the livers of T2DM mice. Discussion ICAM1 and EZR may serve as potential diagnostic biomarkers and candidate therapeutic targets associated with T2DM. Furthermore, ICAM1 and EZR may be associated with insulin resistance by reducing glucose transport efficiency, potentially through modulation of the PI3K-AKT signaling pathway.