Establishment and Evaluation of a 6-Gene Survival Risk Assessment Model Related to Lung Adenocarcinoma Microenvironment

Abstract

Objective. A survival risk assessment model associated with a lung adenocarcinoma (LUAD) microenvironment was established and evaluated to identify effective independent prognostic factors for LUAD. Methods. The public data were downloaded from the TCGA database, and ESTIMATE prediction software was used to score immune cells and stromal cells for tumor purity prediction. The samples were divided into the high-score group and the low-score group by the median value of the immune score (or stromal score). The Wilcoxon test was used for differential analysis. GO and KEGG enrichment analysis of differentially expressed genes (DEGs) was performed using "clusterProfiler" of R package. Meanwhile, univariate and multivariate regression analysis was performed on DEGs to construct a multivariate Cox risk regression model with variable gene expression levels as independent prognostic factors affecting a tumor microenvironment (TME) and tumor immunity. Results. This study found that LUAD patients with high immune cell (stromal cell) infiltration had better prognosis and were in earlier staging. Functional enrichment analysis revealed that most DEGs were related to the proliferation and activation of immune cells or stromal cells. A survival prediction model composed of 6 TME-related genes (CLEC17A, TAGAP, ABCC8, BCAN, FLT3, and CCR2) was established, and finally, the 6 feature genes closely related to the prognosis of LUAD were proved. The AUC value of the ROC curve in this model was 0.7, indicating that the model was reliable. Conclusion. Six genes related to the LUAD microenvironment have a predictive prognostic value in LUAD.