A novel tumor mutation-related long non-coding RNA signature for predicting overall survival and immunotherapy response in lung adenocarcinoma.

ABSTRACT

Background: Immunotherapy has changed the treatment landscape for lung cancer. This study aims to construct a tumor mutation-related model that combines long non-coding RNA (lncRNA) expression levels and tumor mutation levels in tumor genomes to detect the possibilities of the lncRNA signature as an indicator for predicting the prognosis and response to immunotherapy in lung adenocarcinoma (LUAD). Methods: We downloaded the tumor mutation profiles and RNA-seq expression database of LUAD from The Cancer Genome Atlas (TCGA). Differentially expressed lncRNAs were extracted based on the cumulative number of mutations. Cox regression analyses were used to identify the prognostic lncRNA signature, and the prognostic value of the five selected lncRNAs was validated by using survival analysis and the receiver operating characteristic (ROC) curve. We used qPCR to validate the expression of five selected lncRNAs between human lung epithelial and human lung adenocarcinoma cell lines. The ImmuCellAI, immunophenoscore (IPS) scores and Tumor Immune Dysfunction and Exclusion (TIDE) analyses were used to predict the response to immunotherapy for this mutation related lncRNA signature. Results: A total of 162 lncRNAs were detected among the differentially expressed lncRNAs between the Tumor mutational burden (TMB)-high group and the TMB-low group. Then, five lncRNAs (PLAC4, LINC01116, LINC02163, MIR223HG, FAM83A-AS1) were identified as tumor mutation-related candidates for constructing the prognostic prediction model. Kaplan‒Meier curves showed that the overall survival of the low-risk group was significantly better than that of the high-risk group, and the results of the GSE50081 set were consistent. The expression levels of PD1, PD-L1 and CTLA4 in the low-risk group were higher than those in the high-risk group. The IPS scores and TIDE scores of patients in the low-risk group were significantly higher than those in the high-risk group. Conclusion: Our findings demonstrated that the five lncRNAs (PLAC4, LINC01116, LINC02163, MIR223HG, FAM83A-AS1) were identified as candidates for constructing the tumor mutation-related model which may serve as an indicator of tumor mutation levels and have important implications for predicting the response to immunotherapy in LUAD.