Dataset Information

Construction of a five-gene prognostic model based on immune-related genes for the prediction of survival in pancreatic cancer.

ABSTRACT:

Purpose

To identify differentially expressed immune-related genes (DEIRGs) and construct a model with survival-related DEIRGs for evaluating the prognosis of patients with pancreatic cancer (PC).

Methods

Six microarray gene expression datasets of PC from the Gene Expression Omnibus (GEO) and Immunology Database and Analysis Portal (ImmPort) were used to identify DEIRGs. RNA sequencing and clinical data from The Cancer Genome Atlas Program-Pancreatic Adenocarcinoma (TCGA-PAAD) database were used to establish the prognostic model. Univariate, least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analyses were applied to determine the final variables of the prognostic model. The median risk score was used as the cut-off value to classify samples into low- and high-risk groups. The prognostic model was further validated using an internal validation set of TCGA and an external validation set of GSE62452.

Results

In total, 142 DEIRGs were identified from six GEO datasets, 47 were survival-related DEIRGs. A prognostic model comprising five genes (i.e., ERAP2, CXCL9, AREG, DKK1, and IL20RB) was established. High-risk patients had poor survival compared with low-risk patients. The 1-, 2-, 3-year area under the receiver operating characteristic (ROC) curve of the model reached 0.85, 0.87, and 0.93, respectively. Additionally, the prognostic model reflected the infiltration of neutrophils and dendritic cells. The expression of most characteristic immune checkpoints was significantly higher in the high-risk group versus the low-risk group.

Conclusions

The five-gene prognostic model showed reliably predictive accuracy. This model may provide useful information for immunotherapy and facilitate personalized monitoring for patients with PC.

SUBMITTER: Liu B

PROVIDER: S-EPMC8252190 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundPancreatic cancer patients with similar clinicopathological status exhibit substantially different therapeutic responses, which might be caused by the vast molecular heterogeneity of tumors. In this study, we attempted to identify specific molecular subgroups and construct a prognostic prediction model based on the expression level of immune-related genes in pancreatic cancer. The transcriptome profiling, single nucleotide variation, copy number variation, clinicopathological information, and follow-up data of pancreatic cancer patients were obtained from The Cancer Genome Atlas database. Thereafter, the immune-related genes with prognostic significance were identified for further consensus cluster analysis. The molecular characteristics and clinicopathological information were compared between the identified subgroups, and a weighted correlation network analysis was performed to identify the hub genes associated with the subgroups. Finally, the prognostic prediction model based on immune-related genes was established using the least absolute shrinkage and selection operator (LASSO) analysis.ResultsA total of 67 immune-relevant genes with prognostic significance were selected and used for the consensus cluster analysis. The total samples were divided into two groups, C1 and C2. The subgroup C1 had a significantly worse prognosis than C2, as well as lower levels of immune cell infiltration, which indicate an immunosuppressed state. The mutational rate of the cancer-related genes including KRAS, TP53, and RNF43 was higher in the C1 subgroup. The C1 subgroup was associated with more advanced tumor grade and T stage and with higher mortality. Using LASSO regression, we developed a prognostic prediction model based on the expression levels of 19 immune-related genes, which we validated in three external data sets. In addition, we identified four potential therapeutic and prognostic biomarkers (TNNT1, KCNN4, SH2D3A, and PHLDA2).ConclusionWe identified two novel molecular subgroups of pancreatic cancer and developed a prognostic prediction model based on the expression levels of immune-related genes, which could be used in a clinical setting and could aid in unraveling the molecular processes leading to the development of pancreatic cancer.

Project description:BackgroundThough immunotherapy has become one of the standard therapies for colon cancer, the overall effective rate of immunotherapy is very low. Constructing an immune-related genes prognostic index (IRGPI) model may help to predict the response to immunotherapy and clinical outcomes.MethodsDifferentially expressed immune-related genes (DEIRGs) between normal tissues and colon cancer tissues were identified and used to construct the co-expression network. Genes in the module with the most significant differences were further analyzed. Independent prognostic immune-related genes (IRGs) were identified by univariate and multivariate cox regression analysis. Independent prognostic IRGs were used to construct the IRGPI model using the multivariate cox proportional hazards regression model, and the IRGPI model was validated by independent dataset. ROC curves were plotted and AUCs were calculated to estimate the predictive power of the IRGPI model to prognosis. Gene set enrichment analysis (GSEA) was performed to screen the enriched KEGG pathways in the high-risk and low-risk phenotype. Correlations between IRGPI and clinical characteristic, immune checkpoint expression, TMB, immune cell infiltration, immune function, immune dysfunction, immune exclusion, immune subtype were analyzed.ResultsTotally 680 DEIRGs were identified. Three independent IRGs,NR5A2, PPARGC1A and LGALS4, were independently related to survival. NR5A2, PPARGC1A and LGALS4 were used to establish the IRGPI model. Survival analysis showed that patients with high-risk showed worse survival than patients in the low-risk group. The AUC of the IRGPI model for 1-year, 3-year and 5-year were 0.584, 0.608 and 0.697, respectively. Univariate analysis and multivariate cox regression analysis indicated that IRGPI were independent prognostic factors for survival. Stratified survival analysis showed that patients with IRGPI low-risk and low TMB had the best survival, which suggested that combination of TMB and IRGPI can better predict clinical outcome. Immune cell infiltration, immune function, immune checkpoint expression and immune exclusion were different between IRGPI high-risk and low-risk patients.ConclusionAn immune-related genes prognostic index (IRGPI) was constructed and validated in the current study and the IRGPI maybe a potential biomarker for evaluating response to immunotherapy and clinical outcome for colon cancer patients.

Project description:Aims. This study is aimed at identifying a prognostic signature for cervical cancer. Main Methods. The gene expression data and clinical information of cervical cancer and normal cervical tissues were acquired from The Cancer Genome Atlas and from three datasets of the Gene Expression Omnibus database. DESeq2 and Limma were employed to screen differentially expressed genes (DEGs). The overlapping DEGs among all datasets were considered the final DEGs. Then, the functional enrichment analysis was performed. Moreover, the Cox proportional hazards regression was performed to establish a prognostic signature of the DEGs. The Kaplan-Meier analysis was applied to test the model. Relationships between gene expression and clinicopathological parameters in cervical cancer, including age, HPV status, histology, stage, and lymph node metastasis, were analysed by the chi-square test. The somatic mutations of these prognostic genes were assessed through cBioPortal. The robustness of the model was verified in another two independent validation cohorts. Key Findings. In total, 169 overlapping upregulated genes and 29 overlapping downregulated genes were identified in cervical cancer compared with normal cervical tissues. Functional enrichment analysis indicated that the DEGs were mainly enriched in DNA replication, the cell cycle, and the p53 signalling pathway. Finally, a 5-gene- (ITM2A, DSG2, SPP1, EFNA1, and MMP1) based prognostic signature was built. According to this model, each patient was given a prognostic-related risk value. The Kaplan-Meier analysis showed that a higher risk was related to worse overall survival in cervical cancer, with an area under the receiver operating characteristic curve of 0.811 for 15 years. The validity of this model in the prediction of cervical cancer outcome was verified in another two independent datasets. In addition, our study also found that the low expression of ITM2A was associated with cervical adenocarcinoma. Interestingly, DSG2 was associated with the HPV status of cervical cancer. Significance. Our study constructed a prognostic model in cervical cancer and discovered two novel genes, ITM2A and DSG2, associated with cervical carcinogenesis and survival.

Project description:Long non-coding RNAs (lncRNAs) may play a role in oxidative stress by altering the tumor microenvironment, thereby affecting pancreatic cancer progression. There is currently limited information on oxidative stress-related lncRNAs as novel prognostic markers of pancreatic cancer. Gene expression and clinical data of patients with pancreatic cancer were downloaded from The Cancer Genome Atlas (TCGA-PAAD) and the International Cancer Genome Consortium (ICGC-PACA) database. A weighted gene co-expression network analysis (WGCNA) was constructed to identify genes that were differentially expressed between normal and tumor samples. Based on the TCGA-PAAD cohort, a prediction model was established using lasso regression and Cox regression. The TCGA-PAAD and ICGC-PACA cohorts were used for internal and external validation, respectively. Furthermore, a nomogram based on clinical characteristics was used to predict mortality of patients. Differences in mutational status and tumor-infiltrating immune cells between risk subgroups were also explored and model-based lncRNAs were analyzed for potential immune-related therapeutic drugs. A prediction model for 6-lncRNA was established using lasso regression and Cox regression. Kaplan-Meier survival curves and receiver operating characteristic (ROC) curves indicated that patients with lower risk scores had a better prognosis. Combined with Cox regression analysis of clinical features, risk score was an independent factor predicting overall survival of patients with pancreatic cancer in both the TCGA-PAAD and ICGC-PACA cohorts. Mutation status and immune-related analysis indicated that the high-risk group had a significantly higher gene mutation rate and a higher possibility of immune escape, respectively. Furthermore, the model genes showed a strong correlation with immune-related therapeutic drugs. A pancreatic cancer prediction model based on oxidative stress-related lncRNA was established, which may be used as a biomarker related to the prognosis of pancreatic cancer to evaluate the prognosis of pancreatic cancer patients.

Project description:BackgroundPancreatic cancer (PC) is a highly fatal and aggressive disease with its incidence and mortality quite discouraging. An effective prediction model is urgently needed for the accurate assessment of patients' prognosis to assist clinical decision-making.MethodsGene expression data and clinicopathological data of the samples were acquired from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Gene Expression Omnibus (GEO) databases. Differential expressed genes (DEGs) analysis, univariate Cox regression analysis, least absolute shrinkage and selection operator (LASSO) regression analysis, random forest screening and multivariate Cox regression analysis were applied to construct the risk signature. The effectiveness and independence of the model were validated by time-dependent receiver operating characteristic (ROC) curve, Kaplan-Meier (KM) survival analysis and survival point graph in training set, test set, TCGA entire set and GSE57495 set. The validity of the core gene was verified by immunohistochemistry and our own independent cohort. Meanwhile, functional enrichment analysis of DEGs between the high and low risk groups revealed the potential biological pathways. Finally, CMap database and drug sensitivity assay were utilized to identify potential small molecular drugs as the risk model-related treatments for PC patients.ResultsFour histone modification-related genes were identified to establish the risk signature, including CBX8, CENPT, DPY30 and PADI1. The predictive performance of risk signature was validated in training set, test set, TCGA entire set and GSE57495 set, with the areas under ROC curve (AUCs) for 3-year survival were 0.773, 0.729, 0.775 and 0.770 respectively. Furthermore, KM survival analysis, univariate and multivariate Cox regression analysis proved it as an independent prognostic factor. Mechanically, functional enrichment analysis showed that the poor prognosis of high-risk population was related to the metabolic disorders caused by inadequate insulin secretion, which was fueled by neuroendocrine aberration. Lastly, a cluster of small molecule drugs were identified with significant potentiality in treating PC patients.ConclusionsBased on a histone modification-related gene signature, our model can serve as a reliable prognosis assessment tool and help to optimize the treatment for PC patients. Meanwhile, a cluster of small molecule drugs were also identified with significant potentiality in treating PC patients.

Project description:BackgroundLung 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.MethodsAll 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.ResultsThrough 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.ConclusionThe immune lncRNA risk model established by us could better predict the prognosis of patients with LUAD.