Dataset Information

A novel immune-related ceRNA network that predicts prognosis and immunotherapy response in lung adenocarcinoma

ABSTRACT:

Background

The tumor microenvironment plays an important role in the progression and malignancy of lung adenocarcinoma and affects the immunotherapy response. There is increasing evidence that long non-coding RNAs (lncRNAs) as competing endogenous RNAs (ceRNAs) have significant functions in the development and treatment response of various kinds of cancer. This study aimed to explore the association between immune-related lncRNA-microRNA (miRNA)-messenger RNA (mRNA)-ceRNA networks, and the prognosis of and immunotherapy response in lung adenocarcinoma.

Methods

RNA-sequencing (RNA-seq) and miRNA-seq data from The Cancer Genome Atlas (TCGA) were used to evaluate the infiltration of immune cells in lung adenocarcinoma samples by undertaking a single-sample gene set enrichment analysis (ssGSEA) to divide the cells into high and low immune cell infiltration groups. The differentially expressed mRNA (DEmRNA) was further analyzed by a weighted gene co-expression network analysis (WGCNA), search tool for recurring instances of neighboring genes (STRING), and Cytoscape to select hub genes. The ceRNA network was constructed using Cytoscape. Additionally, survival analyses were conducted to screen out prognostic candidate genes.

Results

Seven thousand five hundred and thirty-eight mRNAs, 540 lncRNAs, and 138 miRNAs were found to be differentially expressed between the high and low immune cell infiltration groups. The two DEmRNA modules most significantly associated with immune cell infiltration were further analyzed, and four clusters, including 179 DEmRNAs, were selected based on Molecular Complex Detection (MCODE) scores. The selected DEmRNAs in the four clusters were mainly enriched in pathways involved in regulating the immune response. Ultimately, a ceRNA network of SNHG6-hsa-miR-30e-5p-CYSLTR1 was identified as being associated with the prognosis of and immunotherapy response in lung adenocarcinoma.

Conclusions

The present study extends understandings of immune-related lncRNA-miRNA-mRNA-ceRNA networks and identifies novel targets and a regulatory pathway for anti-tumor immunotherapy.

SUBMITTER: Gong W

PROVIDER: S-EPMC8506752 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundLung adenocarcinoma (LUAD) is the most common form of lung cancer, accounting for 30% of all cases and 40% of all non-small-cell lung cancer cases. Immune-related genes play a significant role in predicting the overall survival and monitoring the status of the cancer immune microenvironment. The present study was aimed at finding an immune-related gene signature for predicting LUAD patient outcomes.MethodsFirst, we chose the TCGA-LUAD project in the TCGA database as the training cohort for model training. For model validating, we found the datasets of GSE72094 and GSE68465 in the GEO database and took them as the candidate cohorts. We obtained 1793 immune-related genes from the ImmPort database and put them into a univariate Cox proportional hazard model to initially look for the genes with potential prognostic ability using the data of the training cohort. These identified genes then entered into a random survival forests-variable hunting algorithm for the best combination of genes for prognosis. In addition, the LASSO Cox regression model tested whether the gene combination can be further shrinkage, thereby constructing a gene signature. The Kaplan-Meier, Cox model, and ROC curve were deployed to examine the gene signature's prognosis in both cohorts. We conducted GSEA analysis to study further the mechanisms and pathways that involved the gene signature. Finally, we performed integrating analyses about the 22 TICs, fully interpreted the relationship between our signature and each TIC, and highlighted some TICs playing vital roles in the signature's prognostic ability.ResultsA nine-gene signature was produced from the data of the training cohort. The Kaplan-Meier estimator, Cox proportional hazard model, and ROC curve confirmed the independence and predictive ability of the signature, using the data from the validation cohort. The GSEA analysis results illustrated the gene signature's mechanism and emphasized the importance of immune-related pathways for the gene signature. 22 TICs immune infiltration analysis revealed resting mast cells' key roles in contributing to gene signature's prognostic ability.ConclusionsThis study discovered a novel immune-related nine-gene signature (BTK, CCR6, S100A10, SEMA3C, GPI, SCG2, TNFRSF11A, CCL20, and DKK1) that predicts LUAD prognosis precisely and associates with resting mast cells strongly.

Project description:Lung adenocarcinoma (LUAD), a malignancy with high incidence and mortality rates worldwide, contains multiple genomic and epigenomic abnormalities. And the useful tumor markers associated with these abnormalities need further investigation. Whereas apoptosis is a form of programmed cell death, the expression of apoptosis-related genes in LUAD and its relationship with prognosis is unclear. In the present study, we identified 64 differentially expressed apoptosis-related genes (DEARGs) that were differentially expressed between LUAD tissue and normal lung tissue. Based on these DEARGs, all LUAD cases were classified into two subtypes using The Cancer Genome Atlas (TCGA) cohort to assess the prognostic value of apoptosis-related genes for survival. An 11-gene signature was established by applying the Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression method to construct a multigene prediction model and classify all LUAD patients in the TCGA cohort into high or low AS-score groups. Patients in the low AS-score group had significantly higher survival and prognosis than those in the high AS-score group. Taking the median risk score of the AS-score, LUAD patients in the GSE68465 cohort were divided into two risk groups, low and high. The overall survival (OS) time was longer in the low AS-score group. Combined with clinical characteristics, the AS-score was an independent predictor of LUAD patients. Gene ontology (GO) and Kyoto Encylopedia of Genes and Genomes (KEGG) analyses showed that the differential genes between the two groups were mainly enriched in cellular immunity. Further analysis revealed higher immune checkpoint protein expression and higher tumor mutational burden (TMB) in the high AS-score group, suggesting better efficacy of immunotherapy in the high AS-score group than the low AS-score group. And the high AS-score group was better in chemotherapy and targeted therapy efficiency. In conclusion, the AS-score constructed based on apoptosis-related genes can predict the prognosis of LUAD patients and provide some guidance for the antitumor treatment of LUAD patients.

Project description:BackgroundLung adenocarcinoma (LUAD) is the most common type of non-small cell lung cancer (NSCLC) with poor survival in advanced stage. Nowadays the rate of nonsmoking patients has dramatically increased and may be associated with the presence of driver mutations. Better understanding of the mutation profile data of nonsmoking LUAD patients are critical to predict survival and provide greater benefits to more patients. The apolipoprotein B mRNA editing enzyme catalytic polypeptide-like (APOBEC) has been shown to play an important role in molecular tumorigenesis of NSCLC. However, the clinical relevance of APOBEC in nonsmoking LUAD remains to be understood.MethodsLUAD patients with somatic mutation and RNA sequencing data obtained from The Cancer Genome Atlas (TCGA) were assessed and screened in the Gene Expression Omnibus. Transcriptome data and mutational signatures were analyzed using R package. Then, we used the least absolute shrinkage and selection operator (LASSO) regression model to construct the APOBEC3 score (APOBEC3 score) model. The prognostic value was evaluated using Kaplan-Meier analysis. Finally, the functional enrichment analysis of differential expressed genes (DEGs) and the immune-related features were also estimated using R package.ResultsBy analyzing the mutational profile data of NSCLC in the TCGA database, we found that different mutation patterns existed between smoking and nonsmoking patients, and the APOBEC3 family played an important role in the mutation pattern of nonsmoking patients with LUAD. We established an APOBEC3 score and found that TCW (W = A or T) mutation counts were significantly greater in the high APOBEC3 score group than in the low APOBEC3 score group. Furthermore, there were different immune feathers and prognostic values between the high and low APOBEC3 score patients, suggesting an independent prognostic factor of APOBEC3 in nonsmoking LUAD patients.ConclusionsWe established a comprehensive view of APOBEC3 mutations in nonsmoking LUAD patients. Our review provides new insights into using the APOBEC3 mutation to predict prognosis and improve the immunotherapy response for future applications.

Project description:BackgroundThere is increasing evidence of the effectiveness of immune checkpoint blockade (ICB) therapy for the treatment of lung adenocarcinoma (LUAD). However, the benefits of ICB therapy vary among LUAD patients. Due to the research dimension, existing biomarkers, such as programmed death-ligand 1 (PD-L1) expression and tumor mutation burden (TMB), could not reflect the complex tumor environment, and had low prediction accuracy of ICB. Therefore, we aimed to uncover a prognostic biomarker that could also predict whether a patient would benefit from ICB therapy and other common treatments from multiple dimensions, so as to improve the prediction accuracy of pre-treatment patients.MethodsBased on the LUAD dataset retrieved from The Cancer Genome Atlas (TCGA) database, 50 immune-related hub genes were identified using weighted gene co-expression network analysis and univariate Cox regression analyses. An immune-related gene prognostic index (IRGPI) was constructed using a Cox proportional-hazards model based on 15 genes and validated using GSE72094 dataset. We tested its prognostic accuracy by Kaplan-Meier (K-M) survival curves of the two datasets and assessed its predictive power by comparing area under curve (AUC) of IRGPI with existing biomarkers. Subsequently, we analyzed the molecular and immune characteristics, and evaluated the benefits of ICB by PD-L1 expression and Tumor Immune Dysfunction and Exclusion (TIDE) analysis, predicted the inhibitory concentration 50 of common treatments drugs for two IRGPI score-related subgroups.ResultsPatients in the IRGPI-high subgroup had lower overall survival (OS) than patients in the IRGPI-low subgroup in K-M survival curve in two cohorts. And IRGPI has AUC values of 0.715, 0.724, and 0.743 in 1, 2, and 3 years, respectively. A higher tumor mutation burden and PD-L1 expression and the tumor microenvironment (TME) landscape demonstrated that IRGPI-high subgroup patients may respond better to ICB therapy. Genomics of Drug Sensitivity in Cancer (GDSC) analysis indicated that the IRGPI-high subgroup showed greater sensitivity to chemotherapy.ConclusionsIRGPI is a prospective biomarker for evaluating whether a patient will benefit from ICB therapy and other treatments, and distinguishing patients with different molecular and immune characteristics.