Abstract
Tumor metastasis imposes metabolic requirements for escaping from primary tissues, producing vulnerability in treatment. This study aimed to explore the metabolic reprogramming relevant to lung adenocarcinoma (LUAD) metastasis and decode the underlying intercellular alterations. Using the gene expression profiles of 394 LUAD samples derived from The Cancer Genome Atlas (TCGA), we identified 11 metastasis-related metabolic genes involved in glycolysis and lipid metabolism, and defined three metabolic reprogramming phenotypes (MP-I, -II, and -III) using unsupervised clustering. MP-III with the highest glycolytic and lowest lipid metabolic levels exhibited the highest metastatic potency and poorest survival in TCGA and six independent cohorts totaling 1,235 samples. Genomic analyses showed that mutations in TP53 and KEAP1, and deletions in SETD2 and PBRM1 might drive metabolic reprogramming in MP-III. Single-cell RNA-sequencing data from LUAD validated a metabolic evolutionary trajectory from normal to MP-II and MP-III, through MP-I. The further intercellular communications revealed that MP-III interacted uniquely with endothelial cells and fibroblasts in the ANGPTL pathway, and had stronger interactions with endothelial cells in the VEGF pathway. Herein, glycolysis-lipid dysregulation patterns suggested metabolic reprogramming phenotypes relevant to metastasis. Further insights into the oncogenic drivers and microenvironmental interactions would facilitate the treatment of LUAD metastasis in the future.