Dataset Information

Identification and Validation of Three PDAC Subtypes and Individualized GSVA Immune Pathway-Related Prognostic Risk Score Formula in Pancreatic Ductal Adenocarcinoma Patients.

ABSTRACT:

Background

With the progress of precision medicine treatment in pancreatic ductal adenocarcinoma (PDAC), individualized cancer-related medical examination and prediction are of great importance in this high malignant tumor and tumor-immune microenvironment with changed pathways highly enrolled in the carcinogenesis of PDAC.

Methods

High-throughput data of pancreatic ductal adenocarcinoma were downloaded from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) database. After batch normalization, the enrichment pathway and relevant scores were identified by the enrichment of immune-related pathway signature using gene set variation analysis (GSVA). Then, cancerous subtype in TCGA and GEO samples was defined through the NMF methods by cancertypes packages in R software, respectively. Subsequently, the significance between the characteristics of each TCGA sample and cancer type and the significant prognosis-related pathway with risk score formula is calculated through t-test and univariate Cox analysis. Next, the prognostic value of gained risk score formula and each significant prognosis-related pathway were validated in TCGA and GEO samples by survival analysis. The pivotal hub genes in the enriched significant prognosis-related pathway are identified and validated, and the TIMER database was used to identify the potential role of hub genes in the PDAC immune environment. The potential role of hub genes is promoting the transdifferentiation of cancer-associated fibroblasts.

Results

The enrichment pathway and relevant scores were identified by GSVA, and 3 subtypes of pancreatic ductal adenocarcinoma were defined in TCGA and GEO samples. The clinical stage, tumor node metastasis classification, and tumor grade are strongly relative to the subtype above in TCGA samples. A risk formula about GSVA significant pathway "GSE45365_WT_VS_IFNAR_KO_CD11B_DC_MCMV_INFECTION_DN ∗ 0.80 + HALLMARK_GLYCOLYSIS ∗ 16.8 + GSE19888_CTRL_VS_T_CELL_MEMBRANES_ACT_MAST_CELL_DN ∗ 14.4" was identified and validated in TCGA and GEO samples through survival analysis with significance. DCN, VCAN, B4GALT7, SDC1, SDC2, B3GALT6, B3GAT3, SDC3, GPC1, and XYLT2 were identified as hub genes in these GSVA significant pathways and validated in silico.

Conclusions

Three pancreatic ductal adenocarcinoma subtypes are identified, and an individualized GSVA immune pathway score-related prognostic risk score formula with 10 hub genes is identified and validated. The predicted function of the 10 upregulated hub genes in tumor-immune microenvironment was promoting the infiltration of cancer-associated fibroblasts. These findings will contribute to the precision medicine of pancreatic ductal adenocarcinoma treatment and tumor immune-related basic research.

SUBMITTER: Zhang D

PROVIDER: S-EPMC8723862 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundThe epithelial-mesenchymal transition (EMT) is a pivotal process for fibrotic disease, embryonic development, and wound healing. Moreover, some evidence has proven that the disorder of EMT also plays an important role in carcinogenesis, especially invasion and metastasis of various tumors (Ritchie et al., 2015). Additionally, gastric adenocarcinoma (GAC) is a common gastrointestinal malignancy which is the fourth most commonly diagnosed tumor. Our study is aimed at identifying the prognostic value of EMT-related genes in gastric adenocarcinoma.MethodsFirstly, high-throughput and clinical data were downloaded from The Cancer Genome Atlas (TCGA) database. 99 differentially expressed EMT-related genes (ERGs) were obtained in these gastric adenocarcinoma data. Secondly, GO and KEGG enrichment analyses show that EMT may promote gastric carcinogenesis. Next, 10 ERGs associated with prognosis of gastric adenocarcinoma patients are screened out by univariate Cox regression, and 6 pivotal prognostic ERGs (MMP8, MMP11, TFDP3, MYB, F2, and CNTN1) are identified through multivariate Cox regression. These 6 genes are confirmed with significant prognostic value in gastric adenocarcinoma through overall survival (OS) analysis. Finally, a risk score formula is constructed and tested in another gastric adenocarcinoma cohort from GEO.Results99 differentially expressed EMT-related genes (ERGs) and their enriched pathways are identified. 10 ERGs are strongly related to the prognosis of GAC patients. A risk score formula of 6 prognosis-related ERGs used to predict the prognosis of gastric adenocarcinoma patients is identified and tested (risk score = 0.448115∗expression value of MMP8 + 0.378892∗expression value of MMP11 - 0.3226∗expression value of MYB + 1.322812∗expression value of TFDP3 + 0.325063∗expression value of F2 + 0.334197∗expression value of CNTN1).ConclusionThis study provides a potential prognostic signature for predicting prognosis of gastric adenocarcinoma patients and molecular insights of EMT in gastric adenocarcinoma, and the formula focusing on the prognosis of gastric adenocarcinoma can be effective.