Molecular subtyping based on immune cell marker genes predicts prognosis and therapeutic response in patients with lung adenocarcinoma

Objective Lung adenocarcinoma (LA) is one of the most common malignancies and is responsible for the greatest number of tumor-related deaths. Our research aimed to explore the molecular subtype signatures of LA to clarify the correlation among the immune microenvironment, clinical outcomes, and therapeutic response. Methods The LA immune cell marker genes (LICMGs) identified by single-cell RNA sequencing (scRNA-seq) analysis were used to discriminate the molecular subtypes and homologous immune and metabolic traits of GSE72094 LA cases. In addition, the model-building genes were identified from 1441 LICMGs by Cox-regression analysis, and a LA immune difference score (LIDscore) was developed to quantify individual differences in each patient, thereby predicting prognosis and susceptibility to immunotherapy and chemotherapy of LA patients. Results Patients of the GSE72094 cohort were divided into two distinct molecular subtypes based on LICMGs: immune activating subtype (Cluster-C1) and metabolically activating subtype (cluster-C2). The two molecular subtypes have distinct characteristics regarding prognosis, clinicopathology, genomics, immune microenvironment, and response to immunotherapy. Among the LICMGs, LGR4, GOLM1, CYP24A1, SFTPB, COL1A1, HLA-DQA1, MS4A7, PPARG, and IL7R were enrolled to construct a LIDscore model. Low-LIDscore patients had a higher survival rate due to abundant immune cell infiltration, activated immunity, and lower genetic variation, but probably the higher levels of Treg cells in the immune microenvironment lead to immune cell dysfunction and promote tumor immune escape, thus decreasing the responsiveness to immunotherapy compared with that of the high-LIDscore patients. Overall, high-LIDscore patients had a higher responsiveness to immunotherapy and a higher sensitivity to chemotherapy than the low-LIDscore group. Conclusions Molecular subtypes based on LICMGs provided a promising strategy for predicting patient prognosis, biological characteristics, and immune microenvironment features. In addition, they helped identify the patients most likely to benefit from immunotherapy and chemotherapy. Keywords: Lung adenocarcinoma, Single-cell RNA sequencing, Prognosis, Molecular subtype, Immunotherapy response