Supervised machine learning identifies impaired mitochondrial quality control in β cells with development of type 2 diabetes

In type 2 diabetes (T2D), molecular pathways driving β cell failure are difficult to resolve with standard single cell analysis. Here we developed an interpretable, supervised machine learning framework that couples sparse rule-based classification (SnakeClassifier), pathway constrained modelling (BlackSwanClassifier), and β cell mitochondrial fitness stratification (Kolmogorov-Arnold Neural Networks KANN), linking and integrating them into disease mechanisms in single cell RNA sequencing (scRNA-seq) from 52 human donors. SnakeClassifier trained on 50 genes accurately predicted T2D at single cell resolution, outperforming classical ensemble machine learning classifier models, and yielded donor level diabetes scores that correlated with chronic hyperglycemia. The clustering of β cell populations (β1-4) revealed a resilient non-diabetic (ND) β1 subtype characterized by preserved β cell identity genes and lower disease risk, whereas T2D β2-4 subtypes exhibited upregulation of genes involved in cellular and mitochondrial stress and suppression of genes promoting oxidative phosphorylation and insulin secretion. Mitophagy emerged as the dominant program linked to T2D and a mitophagy focused BlackSwanClassifier nominated PINK1, BNIP3 , and FUNDC1 as key regulators. PINK1 was enriched in ND β1, decreased with T2D disease score and connected sex stratified mitophagy. We generated a KANN derived mitochondrial fitness index (MFI) integrating mitophagy, mitochondrial proteostasis, biogenesis and oxidative phosphorylation into a single interpretable score (R ² = 0.934 vs module-based mitochondria quality index), which identified mitophagy PINK1, SQSTM1, PRKN and BNIP3 as top contributors to T2D progression. These transparent models unify prediction with T2D disease mechanism and identify the mitophagy receptor PINK1 as a central determinant of β cell metabolic fitness