Dataset Information

Discovering a Four-Gene Prognostic Model Based on Single-Cell Data and Gene Expression Data of Pancreatic Adenocarcinoma.

ABSTRACT:

Background

Pancreatic cancer has a 5-year overall survival lower than 8%. Pancreatic adenocarcinoma (PAAD) is the most common type. This study attempted to explore novel molecular subtypes and a prognostic model through analyzing tumor microenvironment (TME).

Materials and methods

Single-cell RNA sequencing (scRNA-seq) data and expression profiles from public databases were downloaded. Three PAAD samples with single-cell data and 566 samples with gene expression data were included. Seurat was used to identify cell subsets. SVA merged and removed batch effects from multichip datasets. CIBERSORT was used to evaluate the components of different cells in transcriptome, ConsensusClusterPlus was used to identify molecular subtypes, and gene set enrichment analysis was used for functional enrichment analysis. LASSO Cox was performed to construct dimensionality reduction and prognosis model.

Results

Memory B cells (MBCs) were identified to be significantly with PAAD prognosis. Two immune subtypes (IS1 and IS2) with distinct overall survival were constructed. Forty-one DEGs were identified between IS1 and IS2. Four prognostic genes (ANLN, ARNTL2, SERPINB5, and DKK1) were screened to develop a prognostic model. The model was effective in classifying samples into high-risk and low-risk groups with distinct prognosis. Three subgroups of MBCs were identified, where MBC_0 and MBC_1 were differentially distributed between IS1 and IS2, high-risk and low-risk groups.

Conclusions

MBCs were closely involved in PAAD progression, especially MBC_0 and MBC_1 subgroups. The four-gene prognostic model was predictive of overall survival and could guide immunotherapy for patients with PAAD.

SUBMITTER: Huang W

PROVIDER: S-EPMC9253429 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background and aimsPancreatic adenocarcinoma (PAAD) is highly aggressive and characterized by a poor prognosis. Oxidative stress has great impacts on the occurrence and development of tumors. However, the predictive role of oxidative stress related genes on PAAD patients' prognosis remains unclear. In this study, we aimed to construct a prognostic model for PAAD based on oxidative stress genes and to evaluate its predictive value.MethodsThe Cancer Genome Atlas (TCGA) and three Gene Expression Omnibus (GEO) datasets were used to identify differentially expressed oxidative stress genes. Univariate Cox regression, Kaplan-Meier and multivariate Cox regression analysis were used to select genes and to construct a prognosis model. According to the median value of the model's risk score, patients were divided into high and low risk groups, and gene set enrichment analysis (GSEA), immune infiltration and immunotherapy effect, drug resistance and the expression of immune checkpoint related genes and synthetic driver genes of T cell proliferation were analyzed. Finally, the mRNA and protein levels of four genes in PAAD were verified by the clinical proteomic tumor analysis consortium (CPTAC) database and the immunostaining of patients' tissue.Results55 differentially expressed oxidative stress genes were identified, and four genes including MET, FYN, CTTN and CDK1 were selected to construct a prognosis model. GESA indicated that immune related pathways, metabolic pathways and DNA repair pathways were significantly enriched in the high risk group as compared to the low risk group. The frequency of genetic mutations was also significantly higher in high risk groups than that in low risk groups. Moreover, the infiltration level of 23 immune cells as well as the expression of immune checkpoint related and synthetic driver genes of T cell proliferation were significantly altered, with the better immunotherapy effect occurring in low risk group. In patient PAAD tissues, the mRNA and protein levels of these four genes were up-regulated.ConclusionWe have successfully constructed a four oxidative stress gene prognostic model that has important predictive value for PAAD patients, and this model might be a promising guidance for prognostic prediction and efficacy monitoring in clinical individualized therapy.

Project description:Objective: Pancreatic ductal adenocarcinoma (PDAC) is highly lethal. Although progress has been made in the treatment of PDAC, its prognosis remains unsatisfactory. This study aimed to develop novel prognostic genes related to glycolysis in PDAC and to apply these genes to new risk stratification. Methods: In this study, based on the Cancer Genome Atlas (TCGA) PAAD cohort, the expression level of glycolysis-related gene at mRNA level in PAAD and its relationship with prognosis were analyzed. Non-negative matrix decomposition (NMF) clustering was used to cluster PDAC patients according to glycolytic genes. Prognostic glycolytic genes, screened by univariate Cox analysis and LASSO regression analysis were established to calculate risk scores. The differentially expressed genes (DEGs) in the high-risk group and the low-risk group were analyzed, and the signal pathway was further enriched to analyze the correlation between glycolysis genes. In addition, based on RNA-seq data, CIBERSORT was used to evaluate the infiltration degree of immune cells in PDAC samples, and ESTIMATE was used to calculate the immune score of the samples. Results: A total of 319 glycolysis-related genes were retrieved, and all PDAC samples were divided into two clusters by NMF cluster analysis. Survival analysis showed that PDAC patients in cluster 1 had shorter survival time and worse prognosis compared with cluster 2 samples (P < 0.001). A risk prediction model based on 11 glycolysis genes was constructed, according to which patients were divided into two groups, with significantly poorer prognosis in high-risk group than in low-risk group (P < 0.001). Both internal validation and external dataset validation demonstrate good predictive ability of the model (AUC = 0.805, P < 0.001; AUC = 0.763, P < 0.001). Gene aggregation analysis showed that DEGs highly expressed in high-risk group were mainly concentrated in the glycolysis level, immune status, and tumor cell proliferation, etc. In addition, the samples in high-risk group showed immunosuppressed status and infiltrated by relatively more macrophages and less CD8+T cell. Conclusions: These findings suggested that the gene signature based on glycolysis-related genes had potential diagnostic, therapeutic, and prognostic value for PDAC.