Metabolic heterogeneity in early-stage lung adenocarcinoma revealed by RNA-seq and scRNA-seq

Purpose Cancer cells maintain cell growth, division, and survival through altered energy metabolism. However, research on metabolic reprogramming in lung adenocarcinoma (LUAD) is limited Methods We downloaded TCGA and GEO sequencing data. Consistent clustering with the ConsensusClusterPlus package was employed to detect the scores for four metabolism-related pathways. The LUAD samples in the TCGA dataset were clustered with ConsensusClusterPlus, and the optimal number of clusters was determined according to the cumulative distribution function (CDF). The cell score for each sample in the TCGA dataset was calculated using the MCPcounter estimate function of the MCPcounter package. Results We identified two subtypes by scoring the samples based on the 4 metabolism-related pathways and cluster dimensionality reduction. The prognosis of cluster B was obviously poorer than that of cluster A in patients with LUAD. The analysis of single-nucleotide variation (SNV) data showed that the top 15 genes in the four metabolic pathways with the most mutations were TKTL2, PGK2, HK3, EHHADH, GLUD2, PKLR, TKTL1, HADHB, CPT1C, HK1, HK2, PFKL, SLC2A3, PFKFB1, and CPT1A. The IFNγ score of cluster B was significantly higher than that of cluster A. The immune T-cell lytic activity score of cluster B was significantly higher than that of cluster A. We further identified 5 immune cell subsets from single-cell sequencing data. The top 5 marker genes of B cells were IGHM, JCHAIN, IGLC3, IGHA1, and IGKC. The C0 subgroup of monocytes had a higher pentose phosphate pathway (PPP) score than the C6 subgroup. Conclusions Metabolism-related subtypes could be potential biomarkers in the prognosis prediction and treatment of LUAD.