Identification of four metabolic subtypes and key prognostic markers in lung adenocarcinoma based on glycolytic and glutaminolytic pathways

Abstract

Abstract Background Glucose and glutamine are the main energy sources for tumor cells. Whether glycolysis and glutaminolysis play a critical role in driving the molecular subtypes of lung adenocarcinoma (LUAD) is unknown. This study attempts to identify LUAD metabolic subtypes with different characteristics and key genes based on gene transcription profiling data related to glycolysis and glutaminolysis, and to construct prognostic models to facilitate patient outcome prediction. Methods LUAD related data were obtained from the Cancer Genome Atlas and Gene Expression Omnibus, including TCGA-LUAD, GSE42127, GSE68465, GSE72094, GSE29013, GSE31210, GSE30219, GSE37745, GSE50081. Unsupervised consensus clustering was used for the identification of LUAD subtypes. Differential expression analysis, weighted gene co-expression network analysis (WGCNA) and CytoNCA App in Cytoscape 3.9.0 were used for the screening of key genes. The Cox proportional hazards model was used for the construction of the prognostic risk model. Finally, qPCR analysis, immunohistochemistry and immunofluorescence colocalization were used to validate the core genes of the model. Result This study identified four distinct characterized LUAD metabolic subtypes, glycolytic, glutaminolytic, mixed and quiescent types. The glycolytic type had a worse prognosis than the glutaminolytic type. Nine genes (CXCL8, CNR1, AGER, ALB, S100A7, SLC2A1, TH, SPP1, LEP) were identified as hub genes driving the glycolytic/glutaminolytic LUAD. In addition, the risk assessment model constructed based on three genes (SPP1, SLC2A1 and AGER) had good predictive performance and could be validated in multiple independent external LUAD cohorts. These three genes were differentially expressed in LUAD and lung normal tissues, and might be potential prognostic markers for LUAD. Conclusion LUAD can be classified into four different characteristic metabolic subtypes based on the glycolysis- and glutaminolysis-related genes. Nine genes (CXCL8, CNR1, AGER, ALB, S100A7, SLC2A1, TH, SPP1, LEP) may play an important role in the subtype-intrinsic drive. This metabolic subtype classification, provides new biological insights into the previously established LUAD subtypes.