Molecular signatures of tumor progression in pancreatic adenocarcinoma identified by energy metabolism characteristics.
Ontology highlight
ABSTRACT: Background
In this study, we performed a molecular evaluation of primary pancreatic adenocarcinoma (PAAD) based on the comprehensive analysis of energy metabolism-related gene (EMRG) expression profiles.Methods
Molecular subtypes were identified by nonnegative matrix clustering of 565 EMRGs. An overall survival (OS) predictive gene signature was developed and internally and externally validated based on three online PAAD datasets. Hub genes were identified in molecular subtypes by weighted gene correlation network analysis (WGCNA) coexpression algorithm analysis and considered as prognostic genes. LASSO cox regression was conducted to establish a robust prognostic gene model, a four-gene signature, which performed better in survival prediction than four previously reported models. In addition, a novel nomogram constructed by combining clinical features and the 4-gene signature showed high-confidence clinical utility. According to gene set enrichment analysis (GSEA), gene sets related to the high-risk group participate in the neuroactive ligand receptor interaction pathway.Conclusions
In summary, EMRG-based molecular subtypes and prognostic gene models may provide a novel research direction for patient stratification and trials of targeted therapies.
SUBMITTER: Tan C
PROVIDER: S-EPMC9006543 | biostudies-literature | 2022 Apr
REPOSITORIES: biostudies-literature
Publications
Molecular signatures of tumor progression in pancreatic adenocarcinoma identified by energy metabolism characteristics.
Tan Cong C Wang Xin X Wang Xu X Weng Weiwei W Ni Shu-Juan SJ Zhang Meng M Jiang Hesheng H Wang Lei L Huang Dan D Sheng Weiqi W Xu Mi-Die MD
BMC cancer 20220413 1
<h4>Background</h4>In this study, we performed a molecular evaluation of primary pancreatic adenocarcinoma (PAAD) based on the comprehensive analysis of energy metabolism-related gene (EMRG) expression profiles.<h4>Methods</h4>Molecular subtypes were identified by nonnegative matrix clustering of 565 EMRGs. An overall survival (OS) predictive gene signature was developed and internally and externally validated based on three online PAAD datasets. Hub genes were identified in molecular subtypes b ...[more]
Similar Datasets
Project description:BackgroundLung adenocarcinoma (LUAD) is the major non-small-cell lung cancer pathological subtype with poor prognosis worldwide. Herein, we aimed to build an energy metabolism-associated prognostic gene signature to predict patient survival.MethodsThe gene expression profiles of patients with LUAD were downloaded from the TCGA and GEO databases, and energy metabolism (EM)-related genes were downloaded from the GeneCards database. Univariate Cox and LASSO analyses were performed to identify the prognostic EM-associated gene signatures. Kaplan-Meier and receiver operating characteristic (ROC) curves were plotted to validate the predictive effect of the prognostic signatures. A CIBERSORT analysis was used to evaluate the correlation between the risk model and immune cells. A nomogram was used to predict the survival probability of LUAD based on a risk model.ResultsWe constructed a prognostic signature comprising 13 EM-related genes (AGER, AHSG, ALDH2, CIDEC, CYP17A1, FBP1, GNB3, GZMB, IGFBP1, SORD, SOX2, TRH and TYMS). The Kaplan-Meier curves validated the good predictive ability of the prognostic signature in TCGA AND two GEO datasets (p<0.0001, p=0.00021, and p=0.0034, respectively). The area under the curve (AUC) of the ROC curves also validated the predictive accuracy of the risk model. We built a nomogram to predict the survival probability of LUAD, and the calibration curves showed good predictive ability. Finally, a functional analysis also unveiled the different immune statuses between the two different risk groups.ConclusionOur study constructed and verified a novel EM-related prognostic gene signature that could improve the individualized prediction of survival probability in LUAD.
| S-EPMC9047773 | biostudies-literature
Project description:BackgroundMultiple molecular subtypes with distinct clinical outcomes in pancreatic adenocarcinoma (PAAD) have been identified in recent years. Cuproptosis is a new form of cell death that likely involved in tumor progression. However, the cuproptosis-related molecular subtypes as well as its mediated tumor microenvironment (TME) cell infiltration characteristics largely remain unclear.MethodsExpression profiles of 10 cuproptosis-related genes (CRGs) and their association with patient survival, TME, cancer stemness and drug resistance were studied in 33 cancer types using the TCGA pan-cancer data. Using 437 PAAD samples from five cohorts (TCGA-PAAD cohort and four GEO cohorts), we explored the molecular subtypes mediated by CRGs, along with the associated TME cell infiltration. Unsupervised methods were utilized to perform cuproptosis subtype clustering. The cuproptosis score was constructed using the COX regression model with least absolute shrinkage and selection operator regression (LASSO) algorithm to quantify the cuproptosis characteristics of a single tumor.ResultsThe expression of 10 CRGs varies in different cancer types with striking inter- and intra- cancer heterogeneity. We integrated the genomic profiling of the CRGs and identified three distinct cuproptosis subtypes, and found that multi-layer CRG alterations were correlated with patient prognosis and TME cell infiltration characteristics. In addition, a cuproptosis score signature was constructed to predict prognosis, and its clinical impacts were characterized in the TCGA-PAAD cohort. The cuproptosis signature was significantly associated with prognosis, tumor subtypes, CD8 T-cell infiltration, response to immune checkpoint inhibitors (ICIs) and chemotherapeutic drug sensitivity. Furthermore, the expression patterns of CRGs in pancreatic cancer cells and normal controls were validated, which was almost consistent with the results from the public database. The expression level and prognostic predictive capability of DLAT were verified in 97 PAAD patients from our patient cohort.ConclusionsThese findings may help understand the roles of CRGs in PAAD and the molecular characterization of cuproptosis subtypes. In addition, the cuproptosis score could serve as a promising biomarker for predicting prognosis and response to immunotherapy in PAAD patients.
| S-EPMC9844034 | biostudies-literature
Project description:Molecular signatures have been suggested as biomarkers to classify pancreatic ductal adenocarcinoma (PDAC) into two, three, four or five subtypes. Since the robustness of existing signatures is controversial, we performed a systematic evaluation of four established signatures for PDAC stratification across nine publicly available datasets. Clustering revealed inconsistency of subtypes across independent datasets and in some cases a different number of PDAC subgroups than in the original study, casting doubt on the actual number of existing subtypes. Next, we built sixteen classification models to investigate the ability of the signatures for tumor subtype prediction. The overall classification performance ranged from ∼35% to ∼90% accuracy, suggesting instability of the signatures. Notably, permuted subtypes and random gene sets achieved very similar performance. Cellular decomposition and functional pathway enrichment analysis revealed strong tissue-specificity of the predicted classes. Our study highlights severe limitations and inconsistencies that can be attributed to technical biases in sample preparation and tumor purity, suggesting that PDAC molecular signatures do not generalize across datasets. How stromal heterogeneity and immune compartment interplay in the diverging development of PDAC is still unclear. Therefore, a more mechanistic or a cross-platform multi-omic approach seems necessary to extract more robust and clinically exploitable insights.
| S-EPMC9575186 | biostudies-literature
Project description:BackgroundPancreatic ductal adenocarcinoma (PDAC) is one of the deadliest solid tumors. The rapid progression of PDAC results in an advanced stage of patients when diagnosed. However, the dynamic molecular mechanism underlying PDAC progression remains far from clear.MethodsThe microarray GSE62165 containing PDAC staging samples was obtained from Gene Expression Omnibus and the differentially expressed genes (DEGs) between normal tissue and PDAC of different stages were profiled using R software, respectively. The software program Short Time-series Expression Miner was applied to cluster, compare, and visualize gene expression differences between PDAC stages. Then, function annotation and pathway enrichment of DEGs were conducted by Database for Annotation Visualization and Integrated Discovery. Further, the Cytoscape plugin DyNetViewer was applied to construct the dynamic protein-protein interaction networks and to analyze different topological variation of nodes and clusters over time. The phosphosite markers of stage-specific protein kinases were predicted by PhosphoSitePlus database. Moreover, survival analysis of candidate genes and pathways was performed by Kaplan-Meier plotter. Finally, candidate genes were validated by immunohistochemistry in PDAC tissues.ResultsCompared with normal tissues, the total DEGs number for each PDAC stage were 994 (stage I), 967 (stage IIa), 965 (stage IIb), 1027 (stage III), 925 (stage IV), respectively. The stage-course gene expression analysis showed that 30 distinct expressional models were clustered. Kyoto Encyclopedia of Genes and Genomes analysis indicated that the up-regulated DEGs were commonly enriched in five fundamental pathways throughout five stages, including pathways in cancer, small cell lung cancer, ECM-receptor interaction, amoebiasis, focal adhesion. Except for amoebiasis, these pathways were associated with poor PDAC overall survival. Meanwhile, LAMA3, LAMB3, LAMC2, COL4A1 and FN1 were commonly shared by these five pathways and were unfavorable factors for prognosis. Furthermore, by constructing the stage-course dynamic protein interaction network, 45 functional molecular modules and 19 nodes were identified as featured regulators for all PDAC stages, among which the collagen family and integrins were considered as two main regulators for facilitating aggressive progression. Additionally, the clinical relevance analysis suggested that the stage IV featured nodes MLF1IP and ITGB4 were significantly correlated with shorter overall survival. Moreover, 15 stage-specific protein kinases were identified from the dynamic network and CHEK1 was particularly activated at stage IV. Experimental validation showed that MLF1IP, LAMA3 and LAMB3 were progressively increased from tumor initiation to progression.ConclusionsOur study provided a view for a better understanding of the dynamic landscape of molecular interaction networks during PDAC progression and offered potential targets for therapeutic intervention.
| S-EPMC6303882 | biostudies-literature
Project description:Background: Pancreatic adenocarcinoma (PAAD) is a highly malignant gastrointestinal tumor and is associated with an unfavorable prognosis worldwide. Considering the effect of mitochondrial metabolism on the prognosis of pancreatic cancer has rarely been investigated, we aimed to establish prognostic gene markers associated with mitochondrial energy metabolism for the prediction of survival probability in patients with PAAD. Methods: Gene expression data were obtained from The Cancer Genome Atlas and Gene Expression Omnibus databases, and the mitochondrial energy metabolism-related genes were obtained from the GeneCards database. Based on mitochondrial energy metabolism score (MMs), differentially expressed MMRGs were established for MMs-high and MMs-low groups using ssGSEA. After the univariate Cox and least absolute and selection operator (LASSO) analyses, a prognostic MMRG signature was used in the multivariate Cox proportional regression model. Survival and immune cell infiltration analyses were performed. In addition, a nomogram based on the risk model was used to predict the survival probability of patients with PAAD. Finally, the expression of key genes was verified using quantitative polymerase chain reaction and immunohistochemical staining. Intro cell experiments were performed to evaluated the proliferation and invasion of pancreatic cancer cells. Results: A prognostic signature was constructed consisting of two mitochondrial energy metabolism-related genes (MMP11, COL10A1). Calibration and receiver operating characteristic (ROC) curves verified the good predictability performance of the risk model for the survival rate of patients with PAAD. Finally, immune-related analysis explained the differences in immune status between the two subgroups based on the risk model. The high-risk score group showed higher estimate, immune, and stromal scores, expression of eight checkpoint genes, and infiltration of M0 macrophages, which might indicate a beneficial response to immunotherapy. The qPCR results confirmed high expression of MMP11 in pancreatic cancer cell lines, and IHC also verified high expression of MMP11 in clinical pancreatic ductal adenocarcinoma tissues. In vitro cell experiments also demonstrated the role of MMP11 in cell proliferation and invasion. Conclusion: Our study provides a novel two-prognostic gene signature-based on MMRGs-that accurately predicted the survival of patients with PAAD and could be used for mitochondrial energy metabolism-related therapies in the future.
| S-EPMC10987750 | biostudies-literature
Project description:BackgroundExosomes play a role in intercellular communication and participate in the interaction between pancreatic ductal adenocarcinoma (PDAC) cells and immune cells. Macrophages can receive tumor cell-derived exosomes to polarize into M2-type macrophages, which can enhance the invasion and metastasis of pancreatic cancer, leading to poor prognosis. However, the mechanism by which pancreatic cancer cell-derived exosomes promote M2-type macrophages is still unclear.MethodsM2 macrophage-associated exosome-derived key module genes were identified by differentially expressed genes (DEGs) and weighted gene co-expression network analysis (WGCNA) analysis using exoRbase 2.0, The Cancer Genome Atlas (TCGA), and The International Cancer Genome Consortium (ICGC) databases. Multivariate Cox regression analysis was used to identify key prognostic genes and obtain regression coefficients to establish prognostic signature. Immune infiltration, tumor mutations, and GSEA among different risk groups were compared. exoRbase 2.0, Gene Expression Profiling Interactive Analysis 2 (GEPIA2), HPA, and TISCH2 databases were used to further analyze the expression pattern of S100A9 in pancreatic cancer. In vitro experiments, cell-derived exosome isolation, quantitative polymerase chain reaction (qPCR), western blot, flow cytometry analysis, cell transfection, transwell assay, and CCK-8 assay were applied to investigate the roles of S100A9 in macrophage M2 polarization and tumor progression.ResultsThe key genes of PDAC-derived exosomes promoting M2-type macrophage polarization were identified, and a risk score model was established. The risk score is related to the expression of common immune checkpoints, immune score, and stromal score, and the tumor mutational burden and biological function of high- and low-risk groups were also different. S100A9 was positively correlated with M2-type macrophage marker. In addition, scRNA-seq data from the TISCH2 database revealed that S100A9 is predominantly expressed in pancreatic cancer cells and mono/macrophage cells, suggesting that S100A9 in pancreatic cancer cells could be received by macrophages, thereby inducing macrophage polarization. In vitro, we used exosomes from BxPC-3 cell lines to coculture macrophages and found that macrophages were mainly polarized toward M2 type, which further promoted the proliferation and metastasis of PDAC.ConclusionsOur study established a reliable risk score model for PDAC-derived exosomes and M2 macrophages, identified the important role of S100A9 in macrophage M2 polarization, which provides a new strategy for the diagnosis and treatment of PDAC, and strengthened the understanding of the mechanism of tumor development and metastasis.
| S-EPMC10952989 | biostudies-literature
Project description:Molecular assessment using circulating tumor DNA (ctDNA) has not been well-defined. We recruited 61 pancreatic cancer (PC) patients who underwent initial computed tomography (CT) imaging study during first-line chemotherapy. Initial molecular assessment was performed using droplet digital PCR and defined as the change in KRAS-mutated ctDNA before and after treatments, which was classified into five categories: mNT, molecular negative; mCR, complete response; mPR, partial response; mSD, stable disease; mPD, progressive disease. Of 61 patients, 14 diagnosed with PD after initial CT imaging showed significantly worse therapeutic outcomes than 47 patients with disease control. In these 47 patients, initial molecular assessment exhibited significant differences in therapeutic outcomes between patients with and without ctDNA (mPD + mSD vs. mCR + mNT; 13.2 M vs. 21.7 M, P = 0.0029) but no difference between those with mPD and mSD + mCR + mNT, suggesting that the presence of ctDNA had more impact on the therapeutic outcomes than change in its number. Multivariate analysis revealed that it was the only independent prognostic factor (P = 0.0405). The presence of ctDNA in initial molecular assessment predicted early tumor progression and identified PC patients more likely to benefit from chemotherapy.
| S-EPMC10036464 | biostudies-literature
Project description:BackgroundMetabolic reprogramming is an emerging cancer feature that has recently drawn special attention since it promotes tumor cell growth and proliferation. However, the mechanism of the Warburg effect is still largely unknown. This research aimed to reveal the effects of BarH-like homeobox 2 (BARX2) in regulating tumor progression and glucose metabolism in lung adenocarcinoma (LUAD).MethodsExpression of BARX2 was measured by quantitative real-time polymerase chain reaction (qRT-PCR) in LUAD cell line and tissues, and the tumor-promoting function of BARX2 in LUAD cells was detected in vitro and in vivo xenograft models. The metabolic effects of BARX2 were examined by detecting glucose uptake, the production levels of lactate and pyruvate, and the extracellular acidification rate (ECAR). Chromatin immunoprecipitation (ChIP) assay and luciferase reporter gene assay were used to identify the underlying molecular mechanism of BARX2 regulation of HK2. Further studies showed that transcription factor FOXA1 directly interacts with BARX2 and promotes the transcriptional activity of BARX2.ResultsBARX2 was remarkably up-regulated in LUAD tissues and positively linked to advanced clinical stage and poor prognosis. In vitro and in vivo data indicated ectopic expression of BARX2 enhanced cell proliferation and tumorigenesis, whereas BARX2 knockdown suppressed these effects. Metabolic-related experiments showed BARX2 promoted the reprogramming of glucose metabolism. Mechanistically, the BARX2/FOXA1/HK2 axis promoted LUAD progression and energy metabolism reprogramming.ConclusionsIn summary, our research first defined BARX2 as a tumor-promoting factor in LUAD and that it may act as a novel prognostic biomarker and new therapeutic target for the disease.
| S-EPMC9359959 | biostudies-literature
Project description:BackgroundHerein, we investigate the relationship between pancreatic stem cell markers (PCSC markers), CD44, and epithelial-specific antigen (ESA), tumor stroma, and the impact on recurrence outcomes in pancreatic ductal adenocarcinoma (PDAC) patients.MethodsPDAC patients who underwent surgical resection between 01/2012-06/2014 were identified. CD44 and ESA expression was assessed by immunohistochemistry. Stroma was classified as loose, moderate, and dense based on fibroblast content. Overall survival (OS) and relapse-free survival (RFS) were estimated using the Kaplan-Meier method and compared between subgroups by log-rank test. The association between PCSC markers and stroma type was assessed by Fisher's exact test.ResultsN = 93 PDAC patients were identified. The number of PDAC patients with dense, moderate density, and loose stroma was 11 (12%), 51 (54%), and 31 (33%) respectively. PDAC with CD44+/ESA- had highest rate of loose stroma (63%) followed by PDAC CD44+/ESA+ (50%), PDAC CD44-/ESA+ (35%), CD44-/ESA- (9%) (p = 0.0033). Conversely, lack of CD44 and ESA expression was associated with the highest rate of moderate and dense stroma (91% p = 0.0033). No local recurrence was observed in patients with dense stroma and 9 had distant recurrence. The highest rate of cumulative local recurrence was observed in patients with loose stroma. No statistically significant difference in RFS and OS was observed among subgroups (P = 0.089).ConclusionsThese data indicate PCSCs may have an important role in stroma differentiation in PDAC. Our results further suggest that tumor stroma may influence the recurrence pattern in PDAC patients.
| S-EPMC8034174 | biostudies-literature
Project description:Patients with pancreatic ductal adenocarcinoma (PDAC) face a clinically intractable disease with poor survival rates, attributed to exceptionally high levels of metastasis. Epithelial-to-mesenchymal transition (EMT) is pronounced at inflammatory foci within the tumor; however, the immunological mechanisms promoting tumor dissemination remain unclear. It is well established that tumors exhibit the Warburg effect, a preferential use of glycolysis for energy production, even in the presence of oxygen, to support rapid growth. We hypothesized that the metabolic pathways utilized by tumor-infiltrating macrophages are altered in PDAC, conferring a pro-metastatic phenotype. We generated tumor-conditioned macrophages in vitro, in which human peripheral blood monocytes were cultured with conditioned media generated from normal pancreatic or PDAC cell lines to obtain steady-state and tumor-associated macrophages (TAMs), respectively. Compared with steady-state macrophages, TAMs promoted vascular network formation, augmented extravasation of tumor cells out of blood vessels, and induced higher levels of EMT. TAMs exhibited a pronounced glycolytic signature in a metabolic flux assay, corresponding with elevated glycolytic gene transcript levels. Inhibiting glycolysis in TAMs with a competitive inhibitor to Hexokinase II (HK2), 2-deoxyglucose (2DG), was sufficient to disrupt this pro-metastatic phenotype, reversing the observed increases in TAM-supported angiogenesis, extravasation, and EMT. Our results indicate a key role for metabolic reprogramming of tumor-infiltrating macrophages in PDAC metastasis, and highlight the therapeutic potential of using pharmacologics to modulate these metabolic pathways.
| S-EPMC5007961 | biostudies-literature