Cross-talk between necroptosis-related lncRNAs to construct a novel signature and predict the immune landscape of lung adenocarcinoma patients

Abstract

Background: As a new style of cell death, necroptosis plays a crucial role in tumor immune microenvironment. LncRNAs have been identified to act as competitive RNAs to influence genes involved in necroptosis. Therefore, we aim to create a signature based on necroptosis-related lncRNAs to predict the prognosis and immune landscape of lung adenocarcinoma (LUAD) patients in this study. Methods: TCGA database was used to acquire RNA sequencing (RNA-Seq) data and clinical information for 59 lung normal samples and 535 lung adenocarcinoma samples. The Pearson correlation analysis, univariate cox regression analysis and least absolute shrinkage and selection operator (LASSO) cox regression were performed to construct the prognostic NRlncRNAs signature. Then we used Kaplan-Meier (K-M) analysis, time-dependent ROC curves, univariate and multivariate cox regression analysis, and nomogram to validate this signature. In addition, GO, KEGG, and GSVA were analyzed to investigate the potential molecular mechanism. Moreover, we analyzed the relationship between our identified signature and immune microenvironment, TMB, and some clinical characteristics. Finally, we detected the expression of the six necroptosis-related lncRNAs in cells and tissues. Results: We constructed a NRlncRNAs signature consisting of six lncRNAs (FRMD6-AS1, LINC01480, FAM83A-AS1, FRMD6-AS1, MED4-AS1, and LINC01415) in LUAD. LUAD patients with high risk scores had lower chance of survival with an AUC of 0.739, 0.709, and 0.733 for 1-year, 3-year, and 5-year respectively. The results based on GO, KEGG, and GSVA enrichment analysis demonstrated that NRlncRNAs signature-related genes were mainly correlated with immune pathways, metabolic-and cell growth-related pathways, cell cycle, and apoptosis. Moreover, the risk score was correlated with the immune status of LUAD patients. Patients with higher risk scores had lower ESTIMATE scores and higher TIDE scores. The risk score was positively correlated with TMB. LINC01415, FRMD6-AS1 and FAM83A-AS1 were significantly overexpressed in lung adenocarcinoma, while the expression levels of MED4-AS1 and LINC01480 were lower in lung adenocarcinoma. Conclusion: Overall, an innovative prognostic signature based on NRlncRNAs was developed for LUAD through comprehensive bioinformatics analysis, which can act as a predictor of immunotherapy and may provide guidance for clinicians.