Dataset Information

Classification of Lung Adenocarcinoma Based on Immune Checkpoint and Screening of Related Genes.

ABSTRACT:

Aims

Lung adenocarcinoma (LUAD) cells could escape from the monitoring of immune cells and metastasize rapidly through immune escape. Therefore, we aimed to develop a method to predict the prognosis of LUAD patients based on immune checkpoints and their associated genes, thus providing guidance for LUAD treatment.

Methods

Gene sequencing data were downloaded from the Cancer Genome Atlas (TCGA) and analyzed by R software and R Bioconductor software package. Based on immune checkpoint genes, kmdist clustering in ConsensusClusterPlus R software package was utilized to classify LUAD. CIBERSORT was used to quantify the abundance of immune cells in LUAD samples. LM22 signature was performed to distinguish 22 phenotypes of human infiltrating immune cells. Gene set variation analysis (GSVA) was performed on immune checkpoint cluster and immune checkpoint score using GSVA R software package. The risk score was calculated by LASSO regression coefficient. Gene Ontology (GO), Hallmark, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed. PROC was performed to generate the ROC curve and calculate the area under the curve (AUC).

Results

According to the immune checkpoint, LUAD was classified into clusters 1 and 2. Survival rate, immune infiltration patterns, TMB, and immune score were significantly different between the two clusters. Functional prediction showed that the functions of cluster 1 focused on apoptosis, JAK/STAT signaling pathway, TNF-α/NFκB signaling pathway, and STAT5 signaling pathway. The risk score model was constructed based on nine genes associated with immune checkpoints. Survival analysis and ROC analysis showed that patients with high-risk score had poor prognosis. The risk score was significantly correlated with cancer status (with tumor), male proportion, status, tobacco intake, and cancer stage. With the increase of the risk score, the enrichment of 22 biological functions increased, such as p53 signaling pathway. The signature was verified in IMvigor immunotherapy dataset with excellent diagnostic accuracy.

Conclusion

We established a nine-gene signature based on immune checkpoints, which may contribute to the diagnosis, prognosis, and clinical treatment of LUAD.

SUBMITTER: Zhou T

PROVIDER: S-EPMC8337117 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background: In lung adenocarcinoma (LUAD), the predictive role of immune-related subgroup classification in immune checkpoint blockade (ICB) therapy remains largely incomplete. Methods: Transcriptomics analysis was performed to evaluate the association between immune landscape and ICB therapy in lung adenocarcinoma and the associated underlying mechanism. First, the least absolute shrinkage and selection operator (LASSO) algorithm and K-means algorithm were used to identify immune related subgroups for LUAD cohort from the Cancer Genome Atlas (TCGA) database (n = 572). Second, the immune associated signatures of the identified subgroups were characterized by evaluating the status of immune checkpoint associated genes and the immune cell infiltration. Then, potential responses to ICB therapy based on the aforementioned immune related subgroup classification were evaluated via tumor immune dysfunction and exclusion (TIDE) algorithm analysis, and survival analysis and further Cox proportional hazards regression analysis were also performed for LUAD. In the end, gene set enrichment analysis (GSEA) was performed to explore the metabolic mechanism potentially responsible for immune related subgroup clustering. Additionally, two LUAD cohorts from the Gene Expression Omnibus (GEO) database were used as validation cohort. Results: A total of three immune related subgroups with different immune-associated signatures were identified for LUAD. Among them, subgroup 1 with higher infiltration scores for effector immune cells and immune checkpoint associated genes exhibited a potential response to IBC therapy and a better survival, whereas subgroup 3 with lower scores for immune checkpoint associated genes but higher infiltration scores for suppressive immune cells tended to be insensitive to ICB therapy and have an unfavorable prognosis. GSEA revealed that the status of glucometabolic reprogramming in LUAD was potentially responsible for the immune-related subgroup classification. Conclusion: In summary, immune related subgroup clustering based on distinct immune associated signatures will enable us to screen potentially responsive LUAD patients for ICB therapy before treatment, and the discovery of metabolism associated mechanism is beneficial to comprehensive therapeutic strategies making involving ICB therapy in combination with metabolism intervention for LUAD.

Project description:BackgroundLung adenocarcinoma (LUAD), the major lung cancer histotype, represents 40% lung cancers. Currently, outcomes are remarkably different in LUAD patients with similar AJCC/UICC-TNM features. T cell proliferation-related regulator genes (TPRGs) relate to the proliferation, activity and function of T cells and tumor progression. The values of TPRGs in classifying LUAD patients and predicting outcomes remain unknown.MethodsGene expression profile and corresponding clinical data were downloaded from TCGA and the GEO databases. We systematically analyzed the expression profile characteristics of 35 TPRGs in LUAD patients and investigated the differences in overall survival (OS), biology pathway, immunity and somatic mutation between different TPRGs-related subtypes. Subsequently, we constructed a TPRGs-related risk model in TCGA cohort to quantify risk scores using LASSO cox regression analysis and then validated this risk model in two GEO cohorts. LUAD patients were divided into high- and low-risk subtypes according to the median risk score. We systematically compared the biology pathway, immunity, somatic mutation and drug susceptibility between the two risk subtypes. Finally, we validate biological functions of two TPRGs-encoded proteins (DCLRE1B and HOMER1) in LUAD cells A549.ResultsWe identified different TPRGs-related subtypes (including cluster 1/cluster A and its counterpart cluster 2/cluster B). Compared to the cluster 1/cluster A subtype, cluster 2/cluster B subtype tended to have a prominent survival advantage with an immunosuppressive microenvironment and a higher somatic mutation frequency. Then, we constructed a TPRGs-related 6-gene risk model. The high-risk subtype characterized by higher somatic mutation frequency and lower immunotherapy response had a worse prognosis. This risk model was an independent prognostic factor and showed to be reliable and accurate for LUAD classification. Furthermore, subtypes with different risk scores were significantly associated with drug sensitivity. DCLRE1B and HOMER1 suppressed cell proliferation, migration and invasion in LUAD cells A549, which was in line with their prognostic values.ConclusionWe construed a novel stratification model of LUAD based on TPRGs, which can accurately and reliably predict the prognosis and might be used as a predictive tool for LUAD patients.