Abstract
BACKGROUND: Non-small-cell lung cancer (NSCLC) exhibits pronounced molecular heterogeneity, and current predictive models rarely incorporate mitochondrial quality-control programs such as mitophagy. We hypothesize that an artificial intelligence model based on mitophagy-related genes (MRGs) at single-cell resolution could improve prediction of survival and immunotherapy benefit. METHODS: We analyzed single-cell RNA sequencing data from treatment-naïve NSCLC tumors to evaluate the activity of MRGs and identify genes exhibiting differential expression between cells with high versus low mitophagy levels. These differentially expressed genes were then cross-referenced with mitophagy gene sets to pinpoint candidate prognostic markers. Using LASSO regression combined with multiple machine learning classifiers, we constructed a risk model, which was validated in both internal and external cohorts, including a clinical immunotherapy trial. We further examined the relationship between the risk model, immune cell infiltration, and drug sensitivity in silico. The key MRGs were then experimentally validated in A549 cells using qRT-PCR, Western blotting, immunofluorescence, and functional assays for cell migration and wound healing. RESULTS: We quantified the mitochondrial autophagy activity of 18,167 single cells. Differential expression yielded 1,668 genes; intersection with the MRG list produced 39 candidates. A six-gene panel (FOS, CANX, EIF4G1, CALCOCO2, HSP90AB1, and PRKAR1A) emerged from LASSO. Gradient boosting machine (GBM) achieved the optimal cross-validated performance (testing set: AUC = 0.80, validation set: AUC = 0.72). SHAP analysis ranked PRKAR1A and CALCOCO2 as the top risk contributors. Patients classified into the high-MRG-score group exhibited consistently shorter overall survival (OS) across all datasets (HR = 3.66, 95% CI 1.72 − 7.81, P < 0.001). Low-MRG tumors displayed elevated immune and ESTIMATE scores with reduced tumor purity, and achieved a significantly higher objective response rate (35% vs 22%, P < 0.05) and prolonged OS (HR = 1.48, 95% CI 1.07 − 2.05, P = 0.017) with immune checkpoint blockade in the clinical-trial setting. Experimental results showed that knockdown of CALCOCO2 or overexpression of PRKAR1A significantly inhibited A549 cell proliferation and reduced mitochondrial membrane potential (P < 0.05), thereby affecting mitophagy. CONCLUSION: We developed an MRG-based model that reliably stratifies NSCLC patients according to prognosis and identifies those most likely to respond to immune checkpoint inhibitors, providing a framework for integrating tumor metabolic characteristics into personalized therapeutic decisions. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12931-026-03601-w.