# Machine learning-driven identification of shared and disease-specific mitochondria-related genes in COPD, NSCLC, and NSCLC with COPD

**Authors:** Siyu Wu, Zelin Chen, Tongxinwei Sun, Beibei Song, Xinxiu Liu, Liwen Zhang, Jing Li, Haoran Lu, Wenhui Song, Aihong Meng

PMC · DOI: 10.1016/j.isci.2026.114857 · 2026-01-29

## TL;DR

This study identifies blood-based mitochondrial genes that can distinguish COPD, NSCLC, and their coexistence, offering potential diagnostic and treatment options.

## Contribution

The study introduces novel mitochondrial gene signatures and drug repurposing candidates for COPD, NSCLC, and their comorbidity.

## Key findings

- Mitochondrial gene signatures (NDUFB6/BID/COX7A2) achieved high diagnostic accuracy (AUC >0.7) across COPD, NSCLC, and their overlap.
- Machine learning models revealed mitochondria-immune crosstalk in disease progression and immune cell composition changes.
- Metformin and ME-344 were identified as potential repurposed drugs targeting mitochondrial genes in these conditions.

## Abstract

Chronic obstructive pulmonary disease (COPD) and non-small-cell lung cancer (NSCLC) often coexist; here, the shared mitochondrial drivers were investigated. Serum from 30 subjects (seven controls, nine COPD, eight NSCLC, and six NSCLC with COPD) underwent RNA-seq, integrated with 1,136 MitoCarta 3.0-derived mitochondrial-related genes (MRGs). DESeq2 identified 25, 124, and 58 mitochondria-related differentially expressed genes (MR-DEGs) in COPD, NSCLC, and their comorbidity, respectively, with 15 and 58 overlapping genes in relevant pairs. SVM-RFE selected two biomarker sets (3-gene and 5-gene), showing excellent diagnostic performance via ROC (AUC 0.89–0.92) and accurate multivariate logistic regression models. GSEA highlighted immune-inflammatory and oxidative phosphorylation pathways; CIBERSORT revealed altered immune cell proportions (e.g., elevated monocytes in COPD) with biomarker-immune cell correlations. CTD linked BID/COX7A2 to NSCLC, and DGIdb identified metformin/ME-344 as potential drugs. These mitochondrial gene signatures, validated in blood as robust classifiers of COPD, NSCLC, and their overlap, simultaneously furnish diagnostic biomarkers and actionable therapeutic targets, underscoring the translational value of mitochondria-immune crosstalk.

•Blood-based mitochondrial gene signatures distinguish COPD, NSCLC and NSCLC with COPD•NDUFB6/BID/COX7A2 yield ROC-AUC >0.7 across all groups; validated by qPCR•Machine-learning models reveal mitochondria-immune crosstalk in COPD, NSCLC & NSCLC with COPD•Metformin & ME-344 nominated to target NDUFB6/NDUFA1 via drug-repurposing

Blood-based mitochondrial gene signatures distinguish COPD, NSCLC and NSCLC with COPD

NDUFB6/BID/COX7A2 yield ROC-AUC >0.7 across all groups; validated by qPCR

Machine-learning models reveal mitochondria-immune crosstalk in COPD, NSCLC & NSCLC with COPD

Metformin & ME-344 nominated to target NDUFB6/NDUFA1 via drug-repurposing

Health sciences; Medicine; Medical specialty; Internal medicine; Oncology; Respiratory medicine

## Linked entities

- **Genes:** NDUFB6 (NADH:ubiquinone oxidoreductase subunit B6) [NCBI Gene 4712], BID (BH3 interacting domain death agonist) [NCBI Gene 637], COX7A2 (cytochrome c oxidase subunit 7A2) [NCBI Gene 1347], NDUFA1 (NADH:ubiquinone oxidoreductase subunit A1) [NCBI Gene 4694]
- **Chemicals:** metformin (PubChem CID 4091), ME-344 (PubChem CID 68026984)
- **Diseases:** COPD (MONDO:0005002), non-small-cell lung cancer (MONDO:0005233), NSCLC (MONDO:0005233)

## Full-text entities

- **Genes:** NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, CD4 (CD4 molecule) [NCBI Gene 920] {aka CD4mut, IMD79, Leu-3, OKT4D, T4}, NDUFA2 (NADH:ubiquinone oxidoreductase subunit A2) [NCBI Gene 4695] {aka B8, CIB8, MC1DN13}, PRDX5 (peroxiredoxin 5) [NCBI Gene 25824] {aka ACR1, AOEB166, B166, HEL-S-55, PLP, PMP20}, HAGH (hydroxyacylglutathione hydrolase) [NCBI Gene 3029] {aka GLO2, GLO2D, GLX2, GLXII, HAGH1}, COX5B (cytochrome c oxidase subunit 5B) [NCBI Gene 1329] {aka COXVB}, GAPDH (glyceraldehyde-3-phosphate dehydrogenase) [NCBI Gene 2597] {aka G3PD, GAPD, HEL-S-162eP}, DEK (DEK proto-oncogene) [NCBI Gene 7913] {aka D6S231E}, UQCR11 (ubiquinol-cytochrome c reductase, complex III subunit XI) [NCBI Gene 10975] {aka 0710008D09Rik, QCR10, UQCR}, MALAT1 (metastasis associated lung adenocarcinoma transcript 1) [NCBI Gene 378938] {aka HCN, LINC00047, NCRNA00047, NEAT2, PRO2853, miPEP-52}, ALAS2 (5'-aminolevulinate synthase 2) [NCBI Gene 212] {aka ALAS-E, ALASE, ANH1, ASB, SIDBA1, XLDPP}, CYCS (cytochrome c, somatic) [NCBI Gene 54205] {aka CYC, HCS, THC4}, STING1 (stimulator of interferon response cGAMP interactor 1) [NCBI Gene 340061] {aka ERIS, MITA, MPYS, NET23, SAVI, STING}, NDUFB6 (NADH:ubiquinone oxidoreductase subunit B6) [NCBI Gene 4712] {aka B17, CI}, NDUFA13 (NADH:ubiquinone oxidoreductase subunit A13) [NCBI Gene 51079] {aka B16.6, CDA016, CGI-39, GRIM-19, GRIM19, MC1DN28}, TP53 (tumor protein p53) [NCBI Gene 7157] {aka BCC7, BMFS5, LFS1, P53, TRP53}, COX7C (cytochrome c oxidase subunit 7C) [NCBI Gene 1350], Mfn2 (mitofusin 2) [NCBI Gene 170731] {aka D630023P19Rik, Fzo}, TENM1 (teneurin transmembrane protein 1) [NCBI Gene 10178] {aka ODZ1, ODZ3, TEN-M1, TEN1, TNM, TNM1}, NRF1 (nuclear respiratory factor 1) [NCBI Gene 4899] {aka ALPHA-PAL}, NDUFA1 (NADH:ubiquinone oxidoreductase subunit A1) [NCBI Gene 4694] {aka CI-MWFE, MC1DN12, MWFE, ZNF183}, STAT3 (signal transducer and activator of transcription 3) [NCBI Gene 6774] {aka ADMIO, ADMIO1, APRF, HIES}, COX7A2 (cytochrome c oxidase subunit 7A2) [NCBI Gene 1347] {aka COX7AL, COX7AL1, COXVIIAL, COXVIIa-L, VIIAL}, CASP8 (caspase 8) [NCBI Gene 841] {aka ALPS2B, CAP4, Casp-8, FLICE, MACH, MCH5}, STK11 (serine/threonine kinase 11) [NCBI Gene 6794] {aka LKB1, PJS, hLKB1}, IL6 (interleukin 6) [NCBI Gene 3569] {aka BSF-2, BSF2, CDF, HGF, HSF, IFN-beta-2}, Ern1 (endoplasmic reticulum to nucleus signalling 1) [NCBI Gene 78943] {aka 9030414B18Rik, Ire1a, Ire1alpha, Ire1p}, FIS1 (fission, mitochondrial 1) [NCBI Gene 51024] {aka CGI-135, TTC11}, IL17A (interleukin 17A) [NCBI Gene 3605] {aka CTLA-8, CTLA8, IL-17, IL-17A, IL17, ILA17}, GLRX5 (glutaredoxin 5) [NCBI Gene 51218] {aka C14orf87, FLB4739, GRX5, PR01238, PRO1238, PRSA}, GUK1 (guanylate kinase 1) [NCBI Gene 2987] {aka GMK, MTDPS21}, MTOR (mechanistic target of rapamycin kinase) [NCBI Gene 2475] {aka FRAP, FRAP1, FRAP2, RAFT1, RAPT1, SKS}, FKBP8 (FKBP prolyl isomerase 8) [NCBI Gene 23770] {aka FKBP38, FKBPr38}, Aim2 (absent in melanoma 2) [NCBI Gene 383619] {aka Gm1313, Ifi210}, SLC25A39 (solute carrier family 25 member 39) [NCBI Gene 51629] {aka CGI-69, CGI69}, BCL2 (BCL2 apoptosis regulator) [NCBI Gene 596] {aka Bcl-2, PPP1R50}, UQCRB (ubiquinol-cytochrome c reductase binding protein) [NCBI Gene 7381] {aka MC3DN3, QCR7, QP-C, QPC, UQBC, UQBP}, TNF (tumor necrosis factor) [NCBI Gene 7124] {aka DIF, IMD127, TNF-alpha, TNFA, TNFSF2, TNLG1F}, CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, Ndufb7 (NADH:ubiquinone oxidoreductase subunit B7) [NCBI Gene 66916] {aka 1110002H15Rik}, Higd1b (HIG1 domain family, member 1B) [NCBI Gene 75689] {aka 2310056K19Rik}, PRKAA1 (protein kinase AMP-activated catalytic subunit alpha 1) [NCBI Gene 5562] {aka AMPK, AMPK alpha 1, AMPKa1}, MIR34A (microRNA 34a) [NCBI Gene 407040] {aka MIRN34A, miRNA34A, mir-34, mir-34a}, TGFB1 (transforming growth factor beta 1) [NCBI Gene 7040] {aka CAEND1, CED, DPD1, IBDIMDE, LAP, TGF-beta1}, COX6B1 (cytochrome c oxidase subunit 6B1) [NCBI Gene 1340] {aka COX6B, COXG, COXVIb1, MC4DN7}, MAS1L (MAS1 proto-oncogene like, G protein-coupled receptor) [NCBI Gene 116511] {aka MAS-L, MRG, dJ994E9.2}, BID (BH3 interacting domain death agonist) [NCBI Gene 637] {aka FP497}, Fastk (Fas-activated serine/threonine kinase) [NCBI Gene 66587] {aka 0610011K02Rik, 610011K02Rik}, KRAS (KRAS proto-oncogene, GTPase) [NCBI Gene 3845] {aka 'C-K-RAS, C-K-RAS, CFC2, K-RAS2A, K-RAS2B, K-RAS4A}, Mtch2 (mitochondrial carrier 2) [NCBI Gene 56428] {aka 2310034D24Rik, 4930539J07Rik, HSPC032}, Ndufa7 (NADH:ubiquinone oxidoreductase subunit A7) [NCBI Gene 66416] {aka 14.5kDa, 2400007M02Rik, CI-B14.5a}
- **Diseases:** respiratory disorders (MESH:D012131), fatty liver disease (MESH:D005234), heart failure (MESH:D006333), autoimmune diseases (MESH:D001327), COPD (MESH:D029424), tuberculosis (MESH:D014376), pulmonary, cardiac, cerebral, hepatic or renal disease (MESH:D006331), pulmonary fibrosis (MESH:D011658), cor pulmonale (MESH:D011660), metabolic (MESH:D008659), immune dysregulation (OMIM:614878), airway obstruction (MESH:D000402), NSCLC (MESH:D002289), large cell carcinoma (MESH:D018287), adenosquamous carcinoma (MESH:D018196), squamous cell carcinoma (MESH:D002294), systemic diseases (MESH:D034721), airway inflammation (MESH:D007249), neurodegenerative diseases (MESH:D019636), injury (MESH:D014947), mitochondria (MESH:C564971), metastasis (MESH:D009362), immunodeficiency (MESH:D007153), SCLC (MESH:D055752), death (MESH:D003643), diabetic complications (MESH:D048909), Respiratory Diseases (MESH:D012140), mucosal immunity impairment (MESH:D020274), hypertension (MESH:D006973), mitochondrial damage (MESH:D028361), viral infections (MESH:D014777), Lung cancer (MESH:D008175), non-alcoholic fatty liver disease (MESH:D065626), cancer (MESH:D009369), adenocarcinoma (MESH:D000230), caspase (MESH:D056735), lung disease (MESH:D008171), diabetes mellitus (MESH:D003920), liver cirrhosis (MESH:D008103), toxicity (MESH:D064420), chronic cough (MESH:D003371), chronic kidney disease (MESH:D051436)
- **Chemicals:** ATP (MESH:D000255), ME-344 (MESH:C000597890), peroxynitrite (MESH:D030421), Metformin (MESH:D008687), acetyl-CoA (MESH:D000105), MitoQ (MESH:C429014), TRIzol (MESH:C411644), SDS (MESH:D012967), PVDF (MESH:C024865), cholesterol (MESH:D002784), magnesium (MESH:D008274), ROS (MESH:D017382), calcium (MESH:D002118), hydrogen peroxide (MESH:D006861), CHEMBL148831 (-), TERPESTACIN (MESH:C080890), dUTP (MESH:C027078), nitrogen (MESH:D009584), heme (MESH:D006418), EDTA (MESH:D004492)
- **Species:** Gammacoronavirus (genus) [taxon 694013], Mus musculus (house mouse, species) [taxon 10090], Homo sapiens (human, species) [taxon 9606]
- **Mutations:** 3G, S3E, 4 C

## Figures

10 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12925241/full.md

---
Source: https://tomesphere.com/paper/PMC12925241