An Immune-Related Signature Predicts Survival in Patients With Lung Adenocarcinoma

• Published in: Frontiers in oncology Vol. 9; p. 1314
• Main Authors: Zhang, Minghui, Zhu, Kaibin, Pu, Haihong, Wang, Zhuozhong, Zhao, Hongli, Zhang, Jinfeng, Wang, Yan
• Format: Journal Article
• Language: English
• Published: Switzerland Frontiers Media S.A 10.12.2019
• Subjects: hierarchal clustering; immunophenotypes; lung adenocarcinoma; Oncology; patient prognosis; riskscore; hierarchal clustering; patient prognosis; riskscore; lung adenocarcinoma; immunophenotypes
• ISSN: 2234-943X; 2234-943X
• DOI: 10.3389/fonc.2019.01314

We investigated the local immune status and its prognostic value in lung adenocarcinoma. In total, 513 lung adenocarcinoma samples from TCGA and ImmPort databases were collected and analyzed. The R package coxph was employed to mine immune-related genes that were significant prognostic indicators using both univariate and multivariate analyses. The R software package glmnet was then used for Lasso Cox regression analysis, and a prognosis prediction model was constructed for lung adenocarcinoma; clusterProfiler was selected for functional gene annotations and KEGG enrichment analysis. Finally, correlations between the RiskScore and clinical features or signaling pathways were established. Sixty-four immune-related genes remarkably correlated with patient prognosis and were further applied. Samples were hierarchically clustered into two subgroups. Accordingly, the LASSO regression algorithm was employed to screen the 14 most representative immune-related genes ( , and ) with respect to patient prognosis. Then, the prognosis prediction model for lung adenocarcinoma patients (namely, the RiskScore equation) was constructed, and the training set samples were incorporated to evaluate the efficiency of this model to predict and classify patient prognosis. Subsequently, based on functional annotations and KEGG pathway analysis, the 14 immune-related genes were mainly enriched in pathways closely associated with lung adenocarcinoma and its immune microenvironment, such as cytokine-cytokine receptor interaction and human T-cell leukemia virus 1 infection. Furthermore, correlations between the RiskScore and clinical features of the training set samples and signaling pathways (such as p53, cell cycle, and DNA repair) were also demonstrated. Finally, the test set sample data were employed for independent testing and verifying the model. We established a prognostic prediction RiskScore model based on the expression profiles of 14 immune-related genes, which shows high prediction accuracy and stability in identifying immune features. This could provide clinical guidance for the diagnosis and prognosis of different immunophenotypes, and suggest multiple targets for precise advanced lung adenocarcinoma therapy based on subtype-specific immune molecules.