Skip to main content

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



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.


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.


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.


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.

Peer Review reports


Lung cancer remains the leading cancer type that threatens human life expectancy and quality of life worldwide [1]. Non-small cell lung cancer (NSCLC) is the predominant pathological type of all lung cancers, accounting for approximately 85% of cases, while lung adenocarcinoma (LUAD) is the most common subtype of NSCLC [2]. LUAD is characterized by a high rate of recurrence and metastasis. Despite appropriate surgery, chemotherapy, radiotherapy, targeted therapy, and immunotherapy, the 5-year survival of lung cancer patients is still only 16.8%, the prognosis of most LUAD patients remains suboptimal [3]. In addition, the complex pathogenesis of LUAD leads to inaccurate prediction of patient prognosis according to the current TNM staging system [4]. Therefore, it is urgent to decipher the tumorigenesis of LUAD and discover effective biomarkers and potential therapeutic targets to predict the prognosis of LUAD patients.

Cellular metabolic pathways mainly include glycolysis, lipid metabolism, glutaminolysis and oxidative phosphorylation [5]. The metabolism of cancer cells differs from that of normal cells in that they have elevated levels of metabolism and maintain a high proliferation rate that is used to resist some cell death signals [6]. This phenomenon is known as the “Warburg effect” and is closely related to the discovery of new therapeutic targets and the development of new anti-cancer drugs [7]. Angiogenesis is an important feature of tumor growth, but as the tumor grows, it causes some of the tumor tissues to become more distant from the blood vessels, so the content of glucose, oxygen and other components in these tissues decreases [8, 9]. Both glycolytic and glutaminolytic pathways are enhanced in cancer cells. Glycolysis can meet the energy requirements of cancer cells, and glutaminolysis can provide cancer cells with precursors for the synthesis of various substances. It has been found that the level of glucose in tumor extracellular fluid is half that of the surrounding normal tissue [10], a result that is consistent with some previous metabolomics studies reporting reduced glucose levels in malignant tissues [11]. Therefore, in rapidly proliferating and metastasis- and recurrence-prone tumors like LUAD, cancer cells are even more in need of rapid adaptation to glucose deficiency [12]. To adapt to this change, and to compensate for pyruvate, cancer cells maintain the citric acid cycle by accelerating the rate of glutamine catabolism. This series of reactions is of interest to cancer cells, as it can produce glutathione, fatty acids for the citric acid cycle, as well as carbon for nucleotides and even nitrogen for many non-essential amino acids [13]. Glutamine is the most abundant circulating amino acid in blood and muscle, and previous studies have found high intake of glutamine in many cancers including pancreatic, ovarian and breast cancers [14,15,16]. Glutamine is a key amino acid that supports many essential cellular functions in cancer cells, so glutaminolysis is closely related to the development of cancer cells [17].

Inhibition of key enzymes in the glycolytic and glutaminolytic pathways is emerging as a popular area of cancer research, and inhibition of these pathways has been shown to be effective in suppressing cancer cell proliferation [17,18,19,20]. The aim of this study was to identify LUAD metabolic subtypes with different prognosis. Patients were classified into different subtypes based on the expression of based on the glycolysis- and glutaminolysis-related genes. We explored the differences in clinical characteristics including prognosis of LUAD patients with different metabolic subtypes, and screened for possible prognostic markers in LUAD thus establishing a clinically feasible prognostic model that is expected to guide and design targeted therapies for LUAD in the future. The detailed workflow of this study was visible in Fig. 1.

Fig. 1
figure 1

The detailed workflow of this study

Materials and methods

Data download and pre-processing

Five hundred ninety-five LUAD tissues and normal lung tissues, simple nucleotide variations (SNVs), and copy number variations (CNVs) were downloaded from the TCGA GDC portal ( Clinical data and survival data for LUAD were from UCSC Xena Portal. The gene of aerobic glycolytic pathways is derived from WP_AEROBIC_GLYCOLYSIS (n = 12) in the MSigDB database ( The gene (n = 7) of glutaminolytic pathways from a previous study [21]. A total of 18 genes were included in the analysis excluding genes whose expression level was 0 in 50% of the samples. In addition, we also obtained some LUAD microchip datasets with detailed survival information and transcription profile data from the GEO database ( These LUAD datasets are GSE29013 [22], GSE30219 [23], GSE31210 [24, 25], GSE37745 [26,27,28], GSE42127 [29, 30], GSE50081 [31], GSE68465 [32], GSE72094 [33] respectively. We merged the datasets originating from the same platform (GSE29013, GSE30219, GSE31210, GSE37745, GSE50081), and we used the ‘ComBat’ algorithm of ‘sva’ R package to reduce the batch effect due to the non-biotechnical bias [34].

Identification of different metabolic subtypes by consensus clustering

Based on the above 18 genes of the aerobic glycolytic pathway and glutaminolytic pathway, we performed consensus clustering of LUAD samples using the “ConsensusClusterPlus” R package. We refer to previous studies to perform specific practices [35,36,37]. Subsequently, we classified the LUAD samples into different metabolic subtypes based on the median of the two types of co-expressed genes. They were glycolytic type (glycolytic median > 0 and glutaminolytic median < = 0), glutaminolytic type (glycolytic median < = 0 and glutaminolytic median > 0), mixed type (glycolytic median > 0 and glutaminolytic median > 0) and quiescent type (glycolytic median < 0 and glutaminolytic median < 0).

Weighted gene co-expression network analysis (WGCNA) and functional annotations

We explored the genes differentially expressed between glycolytic type and glutaminolytic type tumor tissues using the limma R package. The threshold was set as:|log2FC| > 1 and FDR < 0.05. Differential expressions were presented in the form of volcano maps. Subsequently, a weighted gene co-expression network analysis (WGCNA) of the above differential genes was performed using the WGCNA R package. Accordingly, the correlation between each gene was obtained, and the correlation matrix and topological overlap matrices between the genes were constructed to measure the network connectivity of the genes and to determine the soft threshold size. Genes with similar expression levels were grouped into a gene module with linkage hierarchical clustering. The weight of each module in the dataset was also calculated, and the maximum weight dataset was filtered for subsequent analysis. Subsequently, we selected the modules most related to the glycolytic and glutaminolytic types, and functionally annotated this module genes using the Metascape website (

Screening of hub genes associated to glycolytic\glutaminolytic type

The genes within the significance module obtained above were imported into the STRING database ( to construct a protein-protein interaction (PPI) network, which was also visualized by using Cyctoscape software. The minimum interaction score > 0.5 between proteins was selected in the setup as the basis for the reliability of protein interactions, and the nodes outside the network linkage were hidden to finally obtain the grid map of protein interactions. The topological analysis of the PPI network was performed by CytoNCA app, using “Betweenness”, “Closeness”, “Degree”, “Eigenvector”, “LAC” and “Network” as standard reference, and the genes whose all standard scores were greater than their corresponding standard mean were selected as candidate targets, which were then subjected to the above topological analysis twice to obtain the candidate genes. Subsequently, we further explored the differential expression of these genes between tumor and normal tissues and the relationship with overall survival of patients. Next, we identified the genes that were differentially expressed between LUAD and normal lung tissues and significantly correlated with prognosis as Hub genes.

Construction and validation of a hub gene-based prognostic model

Here, a more prevalent machine learning algorithm - LASSO was used to construct a Hub gene-based prognostic model for LUAD (TCGA-LUAD, n = 500). The relevant parameters in this process were as follows: family = “cox”, maxit = 1200, and other parameters were set as the default values. The risk score for each patient was obtained by the following formula: Riskscore = CoefG1*ExpressionG1 + CoefG2*ExpressionG2 + ... + CoefGn*ExpressionGn. Using the median value of the risk score, patients were divided into two groups. and Kaplan-Meier survival analysis was performed to verify the effectiveness of this model in the training cohort (n = 500) and multiple independent validation cohorts, including GSE42127 cohort (n = 133), GSE68465 cohort (n = 442), GSE72094 cohort (n = 398), and merge-GEO cohort (GSE29013, GSE30219, GSE31210, GSE37745, GSE50081) (n = 574).

Cell culture and quantitative real-time polymerase chain reaction (qRT-PCR)

Lung cancer cell lines (H1299 and A549) and normal lung epithelial cell line (Beas-2B) were purchased from the American Type Culture Collection (ATCC, Manassas, VA). H1299 and A549 cells were cultured using RPMI 1640 medium, and Beas-2B cells were cultured in DMEM medium. These mediums contain 10% fetal bovine serum (FBS), 50 mg/mLstreptomycin and 50 IU/mL penicillin, and all cells were placed in an incubator conditioned at 37 °C, 5% CO2. All cell lines were tested and authenticated by short tandem repeat (STR) analysis. The Hub gene mRNA expression in cells was measured using qRT-PCR analysis. Applying the Trizol method, the Total RNA was extracted and subsequently used to synthesize cDNA and subjected to PCR reactions (all experimental procedures were performed strictly according to the instructions of the kit). The primer sequences of Hub genes, and GAPDH used in the RT–qPCR are listed in Table S1.GAPDH was used as the reference gene and relative gene expression was calculated by the 2 − ΔΔCT method.

Immunohistochemistry (IHC) and subcellular localization of the proteins encoded by the hub genes

In the above analysis, we identified the Hub genes. To further explore the expression levels of the Hub gene-encoded proteins in LUAD tumor tissues and normal lung tissues, we selected their immunohistochemical results for presentation. These IHC results can be obtained from an open-source database (The Human Protein Atlas (HPA): Subsequently, we were interested in the subcellular localization of their proteins, and for this purpose, we also obtained immunofluorescence and confocal images of the subcellular localization of the Hub genes in cancer cells (Hep G2) from this database.


Identification of the four metabolic subtypes of LUAD

Consensus clustering of 11 genes of the glycolytic pathway and 7 genes of the glutaminolytic pathway was performed using TCGA-LUAD data to screen for co-expression of glycolytic and glutaminolytic related genes. When K = 6, glycolytic and glutaminolytic genes were clustered together. As shown in Fig. 2a, the genes in C1 and C2 (defined as glycolytic co-expression genes) all belong to the glycolytic pathway, and these genes include ALDOA, ENO1, GAPDH, GPI, PKM, TPI1, LDHA, PGK1, SLC2A1. The genes in C3, C4 and C6 (defined as glutaminolytic co-expression genes) all belong to the glutaminolytic pathway, and these genes include GLS, GOT1, GPT, GLS2, GLUD1, and GLUD2. Subsequently, the median expression values of co-expressed glycolytic and glutaminolytic genes were z-scored, and four metabolic subgroups were identified based on median expression (Fig. 2b). The expression levels of these selected genes were visualized in the four subgroups, with glutaminolytic genes being generally highly expressed in the glutaminolytic and mixed types but relatively low in the Quiescent and glycolytic types. Glycolytic genes were highly expressed in the glycolytic and mixed types, but low in the Quiescent and glutaminolytic types (Fig. 2c). Next, we investigated the prognostic relationships among these four LUAD subtypes, and the results of Kaplan-Meier survival analysis indicated that OS was significantly different among the four subtypes. Even with multiple hypothesis tests with the Bonferroni method to correct the significance level, the overall survival of glycolytic patients remained significantly worse than the glutaminolytic patients (Fig. 2d). In addition, we also analyzed the distribution of clinical characteristics among the four different subtypes. Extensive similarities in clinical features exist among patients of the four metabolic subtypes, however, they differ in the M stage as well as in the anatomic neoplasm subdivision, as shown in Table 1.

Fig. 2
figure 2

Identification of the four metabolic subtypes of LUAD. A Consistent clustering of the aerobic glycolytic and glutaminolytic pathway related genes. B The four LUAD metabolic subtypes (Glycolytic, Glutaminolytic, Quiescent, and Mixed) were identified according to aerobic glycolytic and glutaminolytic pathway related gene expression levels. C The heatmap showing the expression levels of these selected genes in the four subgroups. D Overall survival time prognostic survival curves of the four molecular subtypes

Table 1 Distribution and comparison of clinical characteristics across the four LUAD metabolic subtypes

Identification of hub genes associated with glycolytic/glutaminolytic type

In the previous analysis, we found significantly different prognostic differences between patients with glycolytic and glutaminolytic types. For this reason, we further investigated the Hub genes significantly associated with the two subtypes. First, we performed differential expression analysis of the transcriptional profiles of the two subtypes (Fig. 3a) and obtained a total of 2772 differentially expressed genes, including 319 genes upregulated in the glycolytic type and 2453 genes upregulated in the glutaminolytic type. We performed WGCNA analysis of glycolytic and glutaminolytic types of samples based on these 2772 differentially expressed genes. To remove sample outliers, we clustered the LUAD glycolytic and glutaminolytic types of samples based on the gene expression matrix and constructed a clustering dendrogram, where the ordinate indicates each sample and the abscissa indicate the clustering distance. The results showed that none of these samples were significantly deviated, and no samples were rejected (Fig. 3b). We set the selection criterion of soft threshold as signed R2 > 0.82, selected a set of candidate thresholds and output the corresponding network parameters, as in Fig. 3c, when the soft threshold was 6, the gene network satisfies both high internal connectivity and high gene similarity. The gene co-expression network was constructed with a threshold value of 6, and the differentially expressed genes were clustered hierarchically according to the dissimilarity matrix, and a clustering dendrogram was constructed (Fig. 3d). The network modules were set to contain at least 50 genes, and the different gene modules were identified using the dynamic cut method, and the modules with high similarity were merged to finally obtain five different gene modules. These different gene modules were indicated by different colors, and the genes in the same color module had high similarity. To screen the modules with high correlation with LUAD glycolytic and glutaminolytic types, we first did principal component analysis (PCA) on the genes in each module separately, extracted the value of the first principal component as the module eigenvalue (ME), and then calculated the correlation coefficient between the ME and glycolytic\glutaminolytic type, and the correlation heat map was shown in Fig. 3e. We found that the turquoise color module was most correlated with glycolysis\glutaminolytic type. The genetic significance of the turquoise color module and the module membership relationships were shown in Fig. 3f. The values of these variables showed a strong positive correlation (cor = 0.83, p = 2.8e-155). In view of the above, we conductd the turquoise color module to Metascape analysis (Table 2). The results showed that these genes are mainly involved in biological processes such as naba matrisome-associated, response to hypoxia, regulation of hormone levels, and response to extracellular stimulus. Subsequently, in order to further screen out the genes related to glycolytic\glutaminolytic type, we input the above turquoise module genes into String database, selected human genes for PPI network analysis (related protein nodes were filtered by protein interaction score ≥ 0.5 criterion), the PPI network TSV format file was imported into Cytoscape 3.9.0 to obtain the PPI network map (Fig. 4a), and the PPI network was topologically analyzed twice by CytoNCA app (Fig. 4b-c), and finally nine candidate genes (CXCL8, CNR1, AGER, ALB, S100A7, SLC2A1, TH, SPP1, LEP) were obtained. Further exploring the expression levels of these nine genes between LUAD tumors and normal tissues and their relationship with prognosis (Fig. 5), we found that only three genes (SPP1, SLC2A1 and AGER) were both differentially expressed and associated with prognosis. Thus, they were identified as glycolytic\glutaminolytic-related Hub genes.

Fig. 3
figure 3

Differentially expressed genes and WGCNA between glycolytic and glutaminolytic samples. A The Volcano plots showing the distribution of differentially expressed genes between the two subtypes. B Cluster plot of the LUAD samples. C Screening of the optimal soft-threshold values. D Dendrogram of LUAD gene clustering. E Heatmap of the correlation of gene modules with the glycolytic / glutaminolytic type. F Correlation of the selected module membership with gene significance

Table 2 The metascape analysis of the module genes
Fig. 4
figure 4

The PPI network and the topology analysis based on CytoNCA App in Cytoscape 3.9.0. A The PPI for the differentially expressed genes between the glycolytic and glutaminolytic samples. B-C The screening of candidate genes by two topological analysis

Fig. 5
figure 5

Expression of the candidate genes and their relationship to the overall survival. A Comparison of the expression of the 11 candidate genes between LUAD tumor and normal tissues. B Effect of 11 candidate genes on LUAD overall survival

Estimation of a prognostic model

In the previous analysis, we identified 3 Hub genes that were closely associated with the glycolytic\glutaminolytic types. Considering their significant impact in terms of patient outcome (they were also considered as prognostic markers of LUAD), we constructed a risk model that can be used to assess the prognosis of LUAD patients based on the transcript expression profiles of these 3 genes using the LASSO method. The risk score for each patient was calculated by the following formula: Riskscore = 0.00236992044579479 * Expression (SPP1) + 0.205240153500597 * Expression (SLC2A1) + (− 0.0441022899873212) * Expression (AGER). Using the median of the patients’ risk scores as the cut-off point, patients were divided into high-risk and low-risk groups. In the TCGA cohort, we can clearly find that the prognosis of patients with high risk score was significantly worse than that of the low risk group (Fig. 6a). Patient survival status and Hub gene expression changed with risk score, more patients died in the high-risk group and the expression of SPP1 and SLC2A1 increased with increasing risk score, while the opposite was true for AGER (Fig. 6b). In addition, independent prognostic analysis also suggested that the risk score could also be used to assess patient prognosis independently of other factors (Table 3). To further validate the stability of the model, we also performed validation in several iindependent external validation cohorts, including the GSE42127 cohort (Fig. 6c), the GSE72094 cohort (Fig. 6d), the GSE68465 cohort (Fig. 6e), and the Meta-GEO cohort (Fig. 6f). Overall, these results suggested that the prognostic model constructed in this study can predict patient outcomes more consistently and accurately.

Fig. 6
figure 6

Validation of the LUAD risk assessment model. A Kaplan-Meier survival curves for patients at high and low risk in TCGA. B Patients of TCGA-LUAD were arranged in the same ascending order of the risk score. C Kaplan-Meier survival curves for patients at high and low risk in GSE42127. D Kaplan-Meier survival curves for patients at high and low risk in GSE72094. E Kaplan-Meier survival curves for patients at high and low risk in GSE68465. F Kaplan-Meier survival curves for patients at high and low risk in merge-GEO (GSE29013, GSE31210, GSE30219, GSE37745, GSE50081)

Table 3 Independent prognostic analysis based on clinical characteristics and risk score

SPP1, SLC2A1, AGER as key prognostic markers for LUAD

In the previous analysis, we identified 3 LUAD prognostic markers, and their expression levels in LUAD were initially revealed by bioinformatic methods. To further clarify this relationship, we performed qPCR analysis in human normal lung epithelial cell line (Beas-2B) and lung cancer cell lines (H1299 and A549) (Fig. 7a). The results showed that SPP1, SLC2A1 was significantly highly expressed in the lung cancer cell, while AGER was in the opposite. These results were consistent with the above findings. In addition, we also verified their expression in LUAD tumors and normal tissues from the protein level (Fig. 7b). Further, we used images obtained by immunofluorescence and confocal microscopy to point out the subcellular localization profiles of the proteins of these three prognostic markers in human cancer cells (Fig. 7c). These findings were beneficial to improve the understanding of these three prognostic markers.

Fig. 7
figure 7

The qPCR, immunofluorescence and confocal microscopy, and immunohistochemistry of the hub genes. A The qPCR revealed the expression of three genes in two LUAD cell lines and normal lung epithelial cells. B The immunohistochemistry of the hub genes in LUAD tissues and normal lung tissues. C Immunofluorescence and confocal microscopy of the hub genes in human cancer cells (Hep G2)


As one of the more heterogeneous malignant tumors, LUAD has its complex oncogenic mechanism [38]. Although the current study has largely increased awareness of LUAD, the treatments for LUAD are significantly diversified, and the prognosis is significantly improved, it is undeniable that there is still a large room for improvement in the prognosis of LUAD. Therefore, it is imperative to further refine the LUAD research and accelerate the individualized treatment of LUAD.

Glucose and glutamine are the main energy sources of tumor cells, and the metabolic abnormalities of both substances significantly affect the tumor cell fate. In view of this, we attempted to classify LUAD patients into differently characterized tumor subtypes based on the transcriptional profile data of the two metabolic genes involved in glycolysis and glutaminolysis. Our results showed that LUAD could be divided into four different metabolic subtypes: Glycolytic, Glutaminolytic, Mixed, and Quiescent. Overall survival remained significantly worse in glycolytic type patients than in glutaminolytic ones. In univariate analysis, we observed that M stage has a large impact on patient prognosis and is a poor prognostic factor. That could be expected. Also, the proportion of M1 patients varies among the different subtypes. Among them, M1 had the highest proportion of Mixed type at 9.45%. This also seems to correspond to a poor prognosis of the Mixed subtype. Interestingly, the glutaminolytic type had a significantly different prognosis from the glycolytic type, however, both did not differ significantly different in M stage. This further emphasizes the deficiency of M stage in the prognostic stratification of LUAD patients, and that the more rational molecular typing of LUAD can still be mined. This illustrates to some extent the importance of the four LUAD metabolic subtypes identified in this study. It is noteworthy that the proportion of patients with M1 was too small in the LUAD sample, which may also be an important factor influencing the prognostic stratification.

Furthermore, our study also identified nine hub genes closely related to the glycolytic/ glutaminolytic LUAD (CXCL8, CNR1, AGER, ALB, S100A7, SLC2A1, TH, SPP1, LEP), which were also further filtered for significant markers affecting LUAD prognosis, thus constructing a risk assessment model. Three key genes (SPP1, SLC2A1 and AGER) were selected for the model construction. The results showed that this risk prognostic model robustly divided LUAD into subgroups with significantly different prognosis, and that this accuracy was also validated in multiple LUAD cohorts. SPP1 is a secreted multifunctional phosphoprotein, also known as bone bridge protein-like protein or early T lymphocyte activation 1 protein, that specifically binds and activates matrix metalloproteinases (MMPs) in cancer [39]. Its main function is to participate in immune response and tissue remodeling, and it is also associated with cell growth, proliferation, migration, and apoptosis [40]. Previous findings have found that SPP1 shows high expression in many cancers and can be used to predict patient prognosis, including ovarian cancer, glioblastoma, hepatocellular carcinoma and gastric cancer [41,42,43], but no studies have been shown to explore the relevance of SPP1 to LUAD,therefore our study may prove to SPP1 in LUAD and its potential clinical value. Among the many genes associated with glucose metabolism, SLC2A1 is the gene encoding a glucose transporter protein that controls glucose uptake and plays a key role in the growth and proliferation of tumor cells [44, 45]. SLC2A1 has been reported to be aberrantly expressed in several cancer types and is closely associated with the development and progression of human cancer [46,47,48]. SLC2A1 was found to be significantly overexpressed in LUAD and closely correlated with overall survival (OS) of patients, which is consistent with the results we obtained. AGER is a highly polymorphic gene with polymorphisms or SNPs that may be responsible or co-responsible for disease development, is expressed primarily in the lung, and is involved in multiple pathways that initiate and maintain an unfavorable pro-inflammatory state [49]. AGER overexpression decreased proliferation, invasion and migration of LUAD cells H1299 and increased apoptosis, AGER may act as a potential molecular marker for LUAD [50], this is very much in line with our expections. Interestingly, another study found that blocking AGER could inhibit cervical squamous cell proliferation and migration, while overexpression of AGER could increase cell proliferation and migration, and inhibit cell apoptosis [51].

Although our study utilized powerful open-source data information to reveal four different features of LUAD metabolic subtypes and to construct a robust risk assessment model, some limitations of this study remain. First, the intrinsic molecular driving mechanisms of the four metabolic isoforms identified in this study were not explored by the underlying experiments. Secondly, the prognostic model we developed was based only on the transcriptional profile information of key genes, and other omics information, such as genome and proteomics information, was lacking. Our prognostic model can refer to more data features at follow-up to further improve the accuracy of the prognostic model.

Overall, our study is the first to identify four distinct LUAD metabolic isotypes based on gene transcription profiling data related to glycolysis and glutaminysis. Nine genes (CXCL8, CNR1, AGER, ALB, S100A7, SLC2A1, TH, SPP1, LEP) may play an important role in the subtype-intrinsic drive. We explored the differences in survival and other clinical characteristics of LUAD patients with different metabolic subtypes, and screened potential prognostic markers in LUAD, thus establishing a clinically feasible prognostic model that is expected to guide and design future targeted therapies for LUAD.

Availability of data and materials

The datasets generated and analysed during the current study are available in the TCGA GDC repository, (, GEO repository, (, Human Protein Atlas (HPA) database (


  1. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics. CA Cancer J Clin. 2022;72(2022):7–33.

    Article  PubMed  Google Scholar 

  2. Song C, Wu Z, Wang Q, Wang Y, Guo Z, Li S, et al. A combined two-mRNA signature associated with PD-L1 and tumor mutational burden for prognosis of lung adenocarcinoma. Front Cell Dev Biol. 2021;9:634697.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Rodriguez-Canales J, Parra-Cuentas E, Wistuba II. Diagnosis and molecular classification of lung Cancer. Cancer Treat Res. 2016;170:25–46.

    Article  PubMed  Google Scholar 

  4. Song C, Lu Z, Lai K, Li D, Hao B, Xu C, et al. Identification of an inflammatory response signature associated with prognostic stratification and drug sensitivity in lung adenocarcinoma. Sci Rep. 2022;12:10110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Feng Z, Ou Y, Hao L. The roles of glycolysis in osteosarcoma. Front Pharmacol. 2022;13:950886.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Tennant DA, Duran RV, Gottlieb E. Targeting metabolic transformation for cancer therapy. Nat Rev Cancer. 2010;10:267–77.

    Article  CAS  PubMed  Google Scholar 

  7. Li C, Zhang G, Zhao L, Ma Z, Chen H. Metabolic reprogramming in cancer cells: glycolysis, glutaminolysis, and Bcl-2 proteins as novel therapeutic targets for cancer. World J Surg Oncol. 2016;14:15.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Martin JD, Fukumura D, Duda DG, Boucher Y, Jain RK. Reengineering the tumor microenvironment to alleviate hypoxia and overcome Cancer heterogeneity. Cold Spring Harb Perspect Med. 2016;6:a027094.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Vaupel P. Tumor microenvironmental physiology and its implications for radiation oncology. Semin Radiat Oncol. 2004;14:198–206.

    Article  PubMed  Google Scholar 

  10. Sullivan MR, Danai LV, Lewis CA, Chan SH, Gui DY, Kunchok T, et al. Quantification of microenvironmental metabolites in murine cancers reveals determinants of tumor nutrient availability. Elife. 2019;8:e44235.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Garcia-Canaveras JC, Chen L, Rabinowitz JD. The tumor metabolic microenvironment: lessons from lactate. Cancer Res. 2019;79:3155–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Grasmann G, Smolle E, Olschewski H, Leithner K. Gluconeogenesis in cancer cells - repurposing of a starvation-induced metabolic pathway? Biochim Biophys Acta Rev Cancer. 1872;2019:24–36.

    Google Scholar 

  13. Xiang Y, Stine ZE, Xia J, Lu Y, O'Connor RS, Altman BJ, et al. Targeted inhibition of tumor-specific glutaminase diminishes cell-autonomous tumorigenesis. J Clin Invest. 2015;125:2293–306.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Fan J, Kamphorst JJ, Mathew R, Chung MK, White E, Shlomi T, et al. Glutamine-driven oxidative phosphorylation is a major ATP source in transformed mammalian cells in both normoxia and hypoxia. Mol Syst Biol. 2013;9:712.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Yang L, Moss T, Mangala LS, Marini J, Zhao H, Wahlig S, et al. Metabolic shifts toward glutamine regulate tumor growth, invasion and bioenergetics in ovarian cancer. Mol Syst Biol. 2014;10:728.

    Article  PubMed  PubMed Central  Google Scholar 

  16. van Geldermalsen M, Wang Q, Nagarajah R, Marshall AD, Thoeng A, Gao D, et al. ASCT2/SLC1A5 controls glutamine uptake and tumour growth in triple-negative basal-like breast cancer. Oncogene. 2016;35:3201–8.

    Article  PubMed  Google Scholar 

  17. Hensley CT, Wasti AT, DeBerardinis RJ. Glutamine and cancer: cell biology, physiology, and clinical opportunities. J Clin Invest. 2013;123:3678–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Mondesir J, Willekens C, Touat M, de Botton S. IDH1 and IDH2 mutations as novel therapeutic targets: current perspectives. J Blood Med. 2016;7:171–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Cho YS, Levell JR, Liu G, Caferro T, Sutton J, Shafer CM, et al. Discovery and evaluation of clinical candidate IDH305, a brain penetrant mutant IDH1 inhibitor. ACS Med Chem Lett. 2017;8:1116–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Shapiro RA, Clark VM, Curthoys NP. Inactivation of rat renal phosphate-dependent glutaminase with 6-diazo-5-oxo-L-norleucine. Evidence for interaction at the glutamine binding site. J Biol Chem. 1979;254:2835–8.

    Article  CAS  PubMed  Google Scholar 

  21. Yu TJ, Ma D, Liu YY, Xiao Y, Gong Y, Jiang YZ, et al. And Di GH, bulk and single-cell transcriptome profiling reveal the metabolic heterogeneity in human breast cancers. Mol Ther. 2021;29:2350–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Xie Y, Xiao G, Coombes KR, Behrens C, Solis LM, Raso G, et al. Robust gene expression signature from formalin-fixed paraffin-embedded samples predicts prognosis of non-small-cell lung cancer patients. Clin Cancer Res. 2011;17:5705–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Rousseaux S, Debernardi A, Jacquiau B, Vitte AL, Vesin A, Nagy-Mignotte H, et al. Ectopic activation of germline and placental genes identifies aggressive metastasis-prone lung cancers. Sci Transl Med. 2013;5:186ra66.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Okayama H, Kohno T, Ishii Y, Shimada Y, Shiraishi K, Iwakawa R, et al. Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. 2012;72:100–11.

    Article  CAS  PubMed  Google Scholar 

  25. Yamauchi M, Yamaguchi R, Nakata A, Kohno T, Nagasaki M, Shimamura T, et al. Epidermal growth factor receptor tyrosine kinase defines critical prognostic genes of stage I lung adenocarcinoma. PLoS One. 2012;7:e43923.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Botling J, Edlund K, Lohr M, Hellwig B, Holmberg L, Lambe M, et al. Biomarker discovery in non-small cell lung cancer: integrating gene expression profiling, meta-analysis, and tissue microarray validation. Clin Cancer Res. 2013;19:194–204.

    Article  CAS  PubMed  Google Scholar 

  27. Jabs V, Edlund K, Konig H, Grinberg M, Madjar K, Rahnenfuhrer J, et al. Integrative analysis of genome-wide gene copy number changes and gene expression in non-small cell lung cancer. PLoS One. 2017;12:e0187246.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Goldmann T, Marwitz S, Nitschkowski D, Krupar R, Backman M, Elfving H, et al. PD-L1 amplification is associated with an immune cell rich phenotype in squamous cell cancer of the lung. Cancer Immunol Immunother. 2021;70:2577–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Tang H, Xiao G, Behrens C, Schiller J, Allen J, Chow CW, et al. A 12-gene set predicts survival benefits from adjuvant chemotherapy in non-small cell lung cancer patients. Clin Cancer Res. 2013;19:1577–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Hight SK, Mootz A, Kollipara RK, McMillan E, Yenerall P, Otaki Y, et al. An in vivo functional genomics screen of nuclear receptors and their co-regulators identifies FOXA1 as an essential gene in lung tumorigenesis. Neoplasia. 2020;22:294–310.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Der SD, Sykes J, Pintilie M, Zhu CQ, Strumpf D, Liu N, et al. Validation of a histology-independent prognostic gene signature for early-stage, non-small-cell lung cancer including stage IA patients. J Thorac Oncol. 2014;9:59–64.

    Article  CAS  PubMed  Google Scholar 

  32. Shedden K, Taylor JM, Enkemann SA, Tsao MS, Yeatman TJ, Gerald WL, et al. Gene expression-based survival prediction in lung adenocarcinoma: a multi-site, blinded validation study. Nat Med. 2008;14:822–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Schabath MB, Welsh EA, Fulp WJ, Chen L, Teer JK, Thompson ZJ, et al. Differential association of STK11 and TP53 with KRAS mutation-associated gene expression, proliferation and immune surveillance in lung adenocarcinoma. Oncogene. 2016;35:3209–16.

    Article  CAS  PubMed  Google Scholar 

  34. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4:249–64.

    Article  PubMed  Google Scholar 

  35. Karasinska JM, Topham JT, Kalloger SE, Jang GH, Denroche RE, Culibrk L, et al. Altered gene expression along the glycolysis-cholesterol synthesis Axis is associated with outcome in pancreatic Cancer. Clin Cancer Res. 2020;26:135–46.

    Article  CAS  PubMed  Google Scholar 

  36. Zhu Z, Qin J, Dong C, Yang J, Yang M, Tian J, et al. Identification of four gastric cancer subtypes based on genetic analysis of cholesterogenic and glycolytic pathways. Bioengineered. 2021;12:4780–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Kim SS, Aprahamian ML, Lindert S. Improving inverse docking target identification with Z-score selection. Chem Biol Drug Des. 2019;93:1105–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Song C, Guo Z, Yu D, Wang Y, Wang Q, Dong Z, et al. A prognostic nomogram combining immune-related gene signature and clinical factors predicts survival in patients with lung adenocarcinoma. Front Oncol. 2020;10:1300.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Su X, Xu BH, Zhou DL, Ye ZL, He HC, Yang XH, et al. Polymorphisms in matricellular SPP1 and SPARC contribute to susceptibility to papillary thyroid cancer. Genomics. 2020;112:4959–67.

    Article  CAS  PubMed  Google Scholar 

  40. Zeng B, Zhou M, Wu H, Xiong Z. SPP1 promotes ovarian cancer progression via integrin beta1/FAK/AKT signaling pathway. Onco Targets Ther. 2018;11:1333–43.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Kijewska M, Kocyk M, Kloss M, Stepniak K, Korwek Z, Polakowska R, et al. The embryonic type of SPP1 transcriptional regulation is re-activated in glioblastoma. Oncotarget. 2017;8:16340–55.

    Article  PubMed  Google Scholar 

  42. Wang J, Hao F, Fei X, Chen Y. SPP1 functions as an enhancer of cell growth in hepatocellular carcinoma targeted by miR-181c. Am J Transl Res. 2019;11:6924–37.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. Song SZ, Lin S, Liu JN, Zhang MB, Du YT, Zhang DD, et al. Targeting of SPP1 by microRNA-340 inhibits gastric cancer cell epithelial-mesenchymal transition through inhibition of the PI3K/AKT signaling pathway. J Cell Physiol. 2019;234:18587–601.

    Article  CAS  PubMed  Google Scholar 

  44. Kunkel M, Reichert TE, Benz P, Lehr HA, Jeong JH, Wieand S, et al. Overexpression of Glut-1 and increased glucose metabolism in tumors are associated with a poor prognosis in patients with oral squamous cell carcinoma. Cancer. 2003;97:1015–24.

    Article  CAS  PubMed  Google Scholar 

  45. DeBerardinis RJ, Cheng T. Q's next: the diverse functions of glutamine in metabolism, cell biology and cancer. Oncogene. 2010;29:313–24.

    Article  CAS  PubMed  Google Scholar 

  46. Feng W, Cui G, Tang CW, Zhang XL, Dai C, Xu YQ, et al. Role of glucose metabolism related gene GLUT1 in the occurrence and prognosis of colorectal cancer. Oncotarget. 2017;8:56850–7.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Koh YW, Lee SJ, Park SY. Differential expression and prognostic significance of GLUT1 according to histologic type of non-small-cell lung cancer and its association with volume-dependent parameters. Lung Cancer. 2017;104:31–7.

    Article  PubMed  Google Scholar 

  48. Oh S, Kim H, Nam K, Shin I. Glut1 promotes cell proliferation, migration and invasion by regulating epidermal growth factor receptor and integrin signaling in triple-negative breast cancer cells. BMB Rep. 2017;50:132–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Serveaux-Dancer M, Jabaudon M, Creveaux I, Belville C, Blondonnet R, Gross C, et al. Pathological implications of receptor for advanced glycation end-product (AGER) gene polymorphism. Dis Markers. 2019;2019:2067353.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Wang Q, Zhu W, Xiao G, Ding M, Chang J, Liao H. Effect of AGER on the biological behavior of nonsmall cell lung cancer H1299 cells. Mol Med Rep. 2020;22:810–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Zhu X, Zhou L, Li R, Shen Q, Cheng H, Shen Z, et al. AGER promotes proliferation and migration in cervical cancer. Biosci Rep. 2018;38. PMID: 29298878.

Download references


This study largely benefited from the large amount of data provided by public databases including GEO, TCGA, HPA, GEPIA2 and so on. We are grateful for the efforts made by the resources and staff to expand and improve the databases.


This work was supported by grants from the Medical Research Project of Wuhan Health Commission (WX21D08).

Author information

Authors and Affiliations



Jinjin Zhang, Congkuan Song and Qi Li designed the study. Jinjin Zhang and Xiaopeng Wang consulted a large number of literature. Jinjin Zhang, and Congkuan Song performed the data analysis, designed and drawn the Figures and Tables in the manuscript. Jinjin Zhang drafted the manuscript. Congkuan Song and Qi Li reviewed and modified this manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Congkuan Song or Qi Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable (NA).

Consent for publication

Not applicable (NA).

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Primer information for the three Hub genes and the reference genes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, J., Wang, X., Song, C. et al. Identification of four metabolic subtypes and key prognostic markers in lung adenocarcinoma based on glycolytic and glutaminolytic pathways. BMC Cancer 23, 152 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: