Dataset Information

Construction of a Prognostic Immune-Related LncRNA Risk Model for Lung Adenocarcinoma.

ABSTRACT:

Background

Lung adenocarcinoma (LUAD) originates mainly from the mucous epithelium and glandular epithelium of the bronchi. It is the most common pathologic subtype of non-small cell lung cancer (NSCLC). At present, there is still a lack of clear criteria to predict the efficacy of immunotherapy. The 5-year survival rate for LUAD patients remains low.

Methods

All data were downloaded from The Cancer Genome Atlas (TCGA) database. We used Gene Set Enrichment Analysis (GSEA) database to obtain immune-related mRNAs. Immune-related lncRNAs were acquired by using the correlation test of the immune-related genes with R version 3.6.3 (Pearson correlation coefficient cor = 0.5, P < 0.05). The TCGA-LUAD dataset was divided into the testing set and the training set randomly. Based on the training set to perform univariate and multivariate Cox regression analyses, we screened prognostic immune-related lncRNAs and given a risk score to each sample. Samples were divided into the high-risk group and the low-risk group according to the median risk score. By the combination of Kaplan-Meier (KM) survival curve, the receiver operating characteristic (ROC) (AUC) curve, the independent risk factor analysis, and the clinical data of the samples, we assessed the accuracy of the risk model. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed on the differentially expressed mRNAs between the high-risk group and the low-risk group. The differentially expressed genes related to immune response between two risk groups were analyzed to evaluate the role of the model in predicting the efficacy and effects of immunotherapy. In order to explain the internal mechanism of the risk model in predicting the efficacy of immunotherapy, we analyzed the differentially expressed genes related to epithelial-mesenchymal transition (EMT) between two risk groups. We extracted RNA from normal bronchial epithelial cell and LUAD cells and verified the expression level of lncRNAs in the risk model by a quantitative real-time polymerase chain reaction (qRT-PCR) test. We compared our risk model with other published prognostic signatures with data from an independent cohort. We transfected LUAD cell with siRNA-LINC0253. Western blot analysis was performed to observed change of EMT-related marker in protein level.

Results

Through univariate Cox regression analysis, 24 immune-related lncRNAs were found to be strongly associated with the survival of the TCGA-LUAD dataset. Utilizing multivariate Cox regression analysis, 10 lncRNAs were selected to establish the risk model. The K-M survival curves and the ROC (AUC) curves proved that the risk model has a fine predictive effect. The GO enrichment analysis indicated that the effect of the differentially expressed genes between high-risk and low-risk groups is mainly involved in immune response and intercellular interaction. The KEGG enrichment analysis indicated that the differentially expressed genes between high-risk and low-risk groups are mainly involved in endocytosis and the MAPK signaling pathway. The expression of genes related to the efficacy of immunotherapy was significantly different between the two groups. A qRT-PCR test verified the expression level of lncRNAs in LUAD cells in the risk model. The AUC of ROC of 5 years in the independent validation dataset showed that this model had superior accuracy. Western blot analysis verified the change of EMT-related marker in protein level.

Conclusion

The immune lncRNA risk model established by us could better predict the prognosis of patients with LUAD.

SUBMITTER: Li Y

PROVIDER: S-EPMC8044985 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Necroptosis is a new form of cell death. Since the discovery that long non-coding RNAs can affect the proliferation of lung adenocarcinoma, much has been learned about it, yet those of necroptosis-related long non-coding RNAs (NRlncRNAs) in lung adenocarcinoma (LUAD) remain enigmatic. This study aims to explore novel biomarkers and therapeutic targets for LUAD. The LUAD data was downloaded from The Cancer Genome Atlas, and necroptosis-related genes were retrieved from published literature. Co-expression analysis, univariate Cox analysis, least absolute shrinkage and selection operator regression analysis were used to identify necroptosis-related prognostic long non-coding RNAs. A comprehensive evaluation of tumor immunity for necrosis-related features was performed, and we identified a 9-NRlncRNA signature. Kaplan-Meier and Cox regression analyses confirmed that the signature was an independent predictor of LUAD outcome in the test and train sets (all P < 0.05). The areas of 1-, 2-, and 3-year overall survival under the time-dependent receiver operating characteristics (ROC) curve (AUC) were 0.754, 0.746, and 0.720, respectively. The GSEA results showed that 9 NRlncRNAs were associated with multiple malignancy-associated and immunoregulatory pathways. Based on this model, we found that the immune status and level of response to chemotherapy and targeted therapy were significantly different in the low-risk group compared with the high-risk group. qRT-PCR assay revealed that 9 NRlncRNAs were involved in the regulation of tumor cell proliferation and may affect the expression of programmed cell death 1 (PD1) and CD28 at human immune checkpoints. Our results indicated that the novel signature involving 9 NRlncRNAs (AL031600.2, LINC01281, AP001178.1, AL157823.2, LINC01290, MED4-AS1, AC026355.2, AL606489.1, FAM83A-AS1) can predict the prognosis of LUAD and are associated with the immune response. This will provide new insights into the pathogenesis and development of therapies for LUAD.

Project description:BackgroundAutophagy related genes (ARGs) regulate lysosomal degradation to induce autophagy, and are involved in the occurrence and development of a variety of cancers. The expression of ARGs in tumor tissues has a great prospect in predicting the survival of patients. The aim of this study was to construct a prognostic risk score model for lung adenocarcinoma (LUAD) based on ARGs.Methods5,786 ARGs were obtained from GeneCards database. Gene expression profiles and clinical data of 395 LUAD patients were collected from The Cancer Genome Atlas (TCGA) database. All ARGs expression data were extracted, and The ARGs differentially expressed were identified by R software. Survival analysis of differentially expressed ARGs was performed to screen for ARGs with prognostic value, and functional enrichment analysis was performed. The least absolute selection operator (LASSO) regression and Cox regression model were used to construct a prognostic risk scoring model for ARGs. The receiver operating characteristic (ROC) curve was drawn to obtain the optimal cut-off value of risk score. According to the cut-off value, the patients were divided into high-risk group and low-risk group. The area under curve (AUC) and the Kaplan-Meier survival curve was plotted to evaluate the model performance, which was verified in external data sets. Finally, univariate and multivariate Cox regression analysis was applied to evaluate the independent prognostic value of the model, and its clinical relevance was analyzed.ResultsSurvival analysis, Lasso regression and Cox regression analysis were used to construct a LUAD prognostic risk score model with five ARGs (ADAM12, CAMP, DKK1, STRIP2 and TFAP2A). The survival time of patients with low-risk score in this model was significantly better than that of patients with high-risk score (P<0.001). The model showed good prediction performance for LUAD in both the training set (AUCmax=0.78) and two external validation sets (AUCmax=0.88). Risk score was significantly associated with the prognosis of LUAD patients in univariate and multivariate Cox regression analyses, suggested that risk score could be a potential independent prognostic factor for LUAD. Correlation analysis of clinical characteristic showed that high risk score was closely associated with high T stage, high tumor stage and poor prognosis.ConclusionsWe constructed a LUAD risk score model consisting of five ARGs, which can provide a reference for predicting the prognosis of LUAD patients, and may be used in combination with tumor node metastasis (TNM) staging for prognosis prediction of LUAD patients in the future.

Project description:BackgroundSuppressive tumor microenvironment is closely related to the progression and poor prognosis of lung adenocarcinoma (LUAD). Novel individual and universal immune-related biomarkers to predict the prognosis and immune landscape of LUAD patients are urgently needed. Two-gene pairing patterns could integrate and utilize various gene expression data.MethodsThe RNA-seq and relevant clinicopathological data of the LUAD project from the TCGA and well-known immune-related genes list from the ImmPort database were obtained. Co-expression analysis followed by an analysis of variance was performed to identify differentially expressed immune-related lncRNA (irlncRNA) (DEirlncRNA) between tumor and normal tissues. Two arbitrary DEirlncRNAs (DEirlncRNAs pair) in a tumor sample underwent pairwise comparison to generate a score (0 or 1). Next, Univariate analysis, Lasso regression and Multivariate analysis were used to screen survival-related DEirlncRNAs pairs and construct a prognostic model. The Acak information standard (AIC) values of the receiver operating characteristic (ROC) curve for 3 years are calculated to determine the cut-off point for high- or low-risk score. Finally, we evaluated the relationship between the risk score and overall survival, clinicopathological features, immune landscape, and chemotherapy efficacy.ResultsData of 54 normal and 497 tumor samples of LUAD were enrolled. After a strict screening process, 15 survival-independent-related DEirlncRNA pairs were integrated to construct a prognostic model. The AUC value of the 3-year ROC curve was 0.828. Kaplan-Meier analysis showed that patients with low risk lived longer than patients with high risk (p <0.001). Univariate and Multivariate Cox analysis suggested that the risk score was an independent factor of survival. The risk score was negatively associated with most tumor-infiltrating immune cells, immune score, and microenvironment scores. The low-risk group was correlated with increased expression of ICOS. The high-risk group had a connection with lower half inhibitory centration (IC50) of most chemotherapy drugs (e.g., etoposide, paclitaxel, vinorelbine, gemcitabine, and docetaxel) and targeted medicine-erlotinib, but with higher IC50 of methotrexate.ConclusionThe established irlncRNA pairs-based model is a promising prognostic signature for LUAD patients. Furthermore, the prognostic signature has great potential in the evaluation of tumor immune landscape and guiding individualized treatment regimens.