A novel messenger RNA signature as a prognostic biomarker for predicting relapse in pancreatic ductal adenocarcinoma

An Immunity-Related Gene Model Predicts Prognosis in Cholangiocarcinoma

The prognosis of patients with cholangiocarcinoma (CCA) is closely related to both immune cell infiltration and mRNA expression. Therefore, we aimed at conducting multi-immune-related gene analyses to improve the prediction of CCA recurrence. Immune-related genes were selected from the Gene Expression Omnibus (GEO), The Cancer Genome Atlas (TCGA), and the Immunology Database and Analysis Portal (ImmPort). The least absolute shrinkage and selection operator (LASSO) regression model was used to establish the multi-gene model that was significantly correlated with the recurrence-free survival (RFS) in two test series. Furthermore, compared with single genes, clinical characteristics, tumor immune dysfunction and exclusion (TIDE), and tumor inflammation signature (TIS), the 8-immune-related differentially expressed genes (8-IRDEGs) signature had a better prediction value. Moreover, the high-risk subgroup had a lower density of B-cell, plasma, B-cell naïve, CD8+ T-cell, CD8+ T-cell naïve, and CD8+ T-cell memory infiltration, as well as more severe immunosuppression and higher mutation counts. In conclusion, the 8-IRDEGs signature was a promising biomarker for distinguishing the prognosis and the molecular and immune features of CCA, and could be beneficial to the individualized immunotherapy for CCA patients.

Overexpression of abnormal spindle-like microcephaly-associated (ASPM) increases tumor aggressiveness and predicts poor outcome in patients with lung adenocarcinoma

Background

Cumulative evidence points to abnormal spindle-like microcephaly-associated (ASPM) protein being overexpressed in various cancers, and the aberrant expression of ASPM has been shown to promote cancer tumorigenicity and progression. However, its role and clinical significance in lung adenocarcinoma (LUAD) remains unclear. This study aimed to determine the expression patterns of ASPM and its clinical significance in LUAD.

Methods

In total, 4 original worldwide LUAD microarray mRNA expression datasets (N=1,116) with clinical and follow-up annotations were downloaded from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. The expression of ASPM protein in LUAD patients was detected by immunohistochemistry. Survival analysis and Cox regression analysis were used to examine the prognostic value of ASPM expression. Gene set enrichment analysis (GSEA) was performed to investigate the relationship between ASPM and LUAD.

Results

Dataset analyses and immunohistochemistry revealed that ASPM expression was significantly higher in the LUAD tissues compared with normal lung tissues, especially in the advanced tumor stage. Additionally, overexpression of ASPM was significantly correlated with shorter overall survival (OS) and relapse-free survival (RFS) in LUAD. Univariate and multivariate Cox regression analyses revealed that the overexpression of ASPM was a potential independent predictor of poor OS and RFS. However, ASPM overexpression was not significantly associated with predicting OS in lung squamous cell carcinoma. GSEA analysis demonstrated that ASPM was significantly enriched in the cell cycle, DNA replication, homologous recombination, RNA degradation, mismatch repair, and p53 signaling pathways.

Conclusions

These findings demonstrate the important role of ASPM in the tumorigenesis and progression of LUAD.

Prediction of early recurrence and response to adjuvant Sorafenib for hepatocellular carcinoma after resection

Background

Early recurrence of hepatocellular carcinoma (HCC) is a major obstacle to improving the prognosis, and no widely accepted adjuvant therapy guideline for patients post-liver resection is available. Currently, all available methods and biomarkers are insufficient to accurately predict post-operation HCC patients’ risk of early recurrence and their response to adjuvant therapy.

Methods

In this study, we downloaded four gene expression datasets (GSE14520, GSE54236, GSE87630, and GSE109211) from the Gene Expression Omnibus database and identified 34 common differentially expressed genes associated with HCC dysregulation and response to adjuvant sorafenib. Then, we constructed a novel 11-messenger RNA predictive model by using ROC curves analysis, univariate Cox regression analysis, and LASSO Cox regression analysis. Furthermore, we validated the predictive values of the risk model in GSE14520 and TCGA-LIHC cohorts by using Kaplan–Meier survival analysis, multivariable Cox regression analysis, and decision curve analysis, respectively.

Results

The risk score model could identify patients with a high risk of HCC recurrence at the early stage and could predict the response of patients to adjuvant sorafenib. Patients with a high risk score had a worse recurrence rate in training cohorts (2-year: p < 0.0001, hazard ratio (HR): 4.658, confidence interval 95% CI [2.895–7.495]; 5-year: p < 0.0001, HR: 3.251, 95% CI [2.155–4.904]) and external validation cohorts (2-year: p < 0.001, HR: 3.65, 95% CI [2.001–6.658]; 5-year: p < 0.001, HR: 3.156, 95% CI [1.78–5.596]). The AUC values of the risk score model for predicting tumor early recurrence were 0.746 and 0.618, and that of the risk score model for predicting the response to adjuvant sorafenib were 0.722 and 0.708 in the different cohort, respectively. Multivariable Cox regression analysis and decision curve analysis also showed that the risk score model was superior to and independent of other clinicopathologic characteristics. Moreover, the risk score model had excellent abilities to predict the overall survival and HCC recurrence of patients with the same tumor stage category.

Conclusions

Our risk model is a reliable and superior predictive tool. With this model, we could optimize the risk stratification based on early tumor recurrence and could evaluate the response of patients to adjuvant sorafenib after liver resection.

A Novel Ferroptosis-Related Gene Signature Predicts Recurrence in Patients With Pancreatic Ductal Adenocarcinoma

Background: Recurrence after surgery is largely responsible for the extremely poor outcomes for patients with pancreatic ductal adenocarcinoma (PDAC). Ferroptosis is implicated in chemotherapy sensitivity and tumor recurrence, we aimed to find out survival-associated ferroptosis-related genes and use them to build a practical risk model with the purpose to predict PDAC recurrence. Methods: Univariate Cox regression analysis was conducted to obtain prognostic ferroptosis-related genes in The Cancer Genome Atlas (TCGA, N = 140) cohort. Multivariate Cox regression analysis was employed to construct a reliable and credible gene signature. The prognostic performance was verified in a MTAB-6134 (N = 286) validation cohort and a PACA-CA (N = 181) validation cohort. The stability of the signature was tested in TCGA and MTAB-6134 cohorts by ROC analyses. Pathway enrichment analysis was adopted to preliminary illuminate the biological relevance of the gene signature. Results: Univariate and multivariate Cox regression analyses identified a 5-gene signature that contained CAV1, DDIT4, SLC40A1, SRXN1 and TFAP2C. The signature could efficaciously stratify PDAC patients with different recurrence-free survival (RFS), both in the training and validation cohorts. Results of subgroup receiver operating characteristic curve (ROC) analyses confirmed the stability and the independence of this signature. Our signature outperformed clinical indicators and previous reported models in predicting RFS. Moreover, the signature was found to be closely associated with several cancer-related and drug response pathways. Conclusion: This study developed a precise and concise prognostic model with the clinical implication in predicting PDAC recurrence. These findings may facilitate individual management of postoperative recurrence in patients with PDAC.

A novel 20-gene prognostic score in pancreatic adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is among the most lethal cancers. Known risk factors for this disease are currently insufficient in predicting mortality. In order to better prognosticate patients with PDAC, we identified 20 genes by utilizing publically available high-throughput transcriptomic data from GEO, TCGA and ICGC which are associated with overall survival and event-free survival. A score generated based on the expression matrix of these genes was validated in two independent cohorts. We find that this “Pancreatic cancer prognostic score 20 –PPS20” is independent of the confounding factors in multivariate analyses, is dramatically elevated in metastatic tissue compared to primary tumor, and is higher in primary tumors compared to normal pancreatic tissue. Transcriptomic analyses show that tumors with low PPS20 have overall more immune cell infiltration and a higher CD8 T cell/Treg ratio when compared to those with high PPS20. Analyses of proteomic data from TCGA PAAD indicated higher levels of Cyclin B1, RAD51, EGFR and a lower E-cadherin/Fibronectin ratio in tumors with high PPS20. The PPS20 score defines not only prognostic and biological sub-groups but can predict response to targeted therapy as well. Overall, PPS20 is a stronger and more robust transcriptomic signature when compared to similar, previously published gene lists.

Identification of a robust functional subpathway signature for pancreatic ductal adenocarcinoma by comprehensive and integrated analyses

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal malignancy and its mortality continues to rise globally. Because of its high heterogeneity and complex molecular landscapes, published gene signatures have demonstrated low specificity and robustness. Functional signatures containing a group of genes involved in similar biological functions may display a more robust performance.

METHODS

The present study was designed to excavate potential functional signatures for PDAC by analyzing maximal number of datasets extracted from available databases with a recently developed method of FAIME (Functional Analysis of Individual Microarray Expression) in a comprehensive and integrated way.

RESULTS

Eleven PDAC datasets were extracted from GEO, ICGC and TCGA databases. By systemically analyzing these datasets, we identified a robust functional signature of subpathway (path:00982_1), which belongs to the drug metabolism-cytochrome P450 pathway. The signature has displayed a more powerful and robust capacity in predicting prognosis, drug response and chemotherapeutic efficacy for PDAC, particularly for the classical subtype, in comparison with published gene signatures and clinically used TNM staging system. This signature was verified by meta-analyses and validated in available cell line and clinical datasets with chemotherapeutic efficacy.

CONCLUSION

The present study has identified a novel functional PDAC signature, which has the potential to improve the current systems for predicting the prognosis and monitoring drug response, and to serve a linkage to therapeutic options for combating PDAC. However, the involvement of path:00982_1 subpathway in the metabolism of anti-PDAC chemotherapeutic drugs, particularly its biological interpretation, requires a further investigation. Video Abstract.

Prognostic values of a novel multi-mRNA signature for predicting relapse of cholangiocarcinoma

Cholangiocarcinoma (CCA) is an epithelial cancer and has high death and recurrence rates, current methods cannot satisfy the need for predicting cancer relapse effectively. So, we aimed at conducting a multi-mRNA signature to improve the relapse prediction of CCA. We analyzed mRNA expression profiling in large CCA cohorts from the Gene Expression Omnibus (GEO) database (GSE76297, GSE32879, GSE26566, GSE31370, and GSE45001) and The Cancer Genome Atlas (TCGA) database. The Least absolute shrinkage and selection operator (LASSO) regression model was used to establish a 7-mRNA-based signature that was significantly related to the recurrence-free survival (RFS) in two test series. Based on the 7-mRNA signature, the cohort TCGA patients could be divided into high-risk or low-risk subgroups with significantly different RFS [p < 0.001, hazard ratio (HR): 48.886, 95% confidence interval (CI): 6.226-3.837E+02]. Simultaneously, the prognostic value of the 7-mRNA signature was confirmed in clinical samples of Ren Ji hospital (p < 0.001, HR: 4.558, 95% CI: 1.829-11.357). Further analysis including multivariable and sub-group analyses revealed that the 7-mRNA signature was an independent prognostic value for recurrence of patients with CCA. In conclusion, our results might provide an efficient tool for relapse prediction and were beneficial to individualized management for CCA patients.

The differential distributions of ASPM isoforms and their roles in Wnt signaling, cell cycle progression, and pancreatic cancer prognosis

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive and treatment‐resistant malignancy. The lack of pathway‐informed biomarkers hampers the development of rational diagnostics or therapies. Recently, the protein abnormal spindle‐like microcephaly‐associated (ASPM) was identified as a novel Wnt and stemness regulator in PDAC, while the pathogenic roles of its protein isoforms remain unclarified. We developed novel isoform‐specific antibodies and genetic knockdown (KD) of putative ASPM isoforms, whereby we uncovered that the levels of ASPM isoform 1 (iI) and ASPM‐iII are variably upregulated in PDAC cells. ASPM isoforms show remarkably different subcellular locations; specifically, ASPM‐iI is exclusively localized to the cortical cytoplasm of PDAC cells, while ASPM‐iII is predominantly expressed in cell nuclei. Mechanistically, ASPM‐iI co‐localizes with disheveled‐2 and active β‐catenin as well as the stemness marker aldehyde dehydrogenase‐1 (ALDH‐1), and its expression is indispensable for the Wnt activity, stemness, and the tumorigenicity of PDAC cells. By contrast, ASPM‐iII selectively regulates the expression level of cyclin E and cell cycle progression in PDAC cells. The expression of ASPM‐iI and ASPM‐iII displays considerable intratumoral heterogeneity in PDAC tissues and only that of ASPM‐iI was prognostically significant; it outperformed ALDH‐1 staining and clinico‐pathological variables in a multivariant analysis. Collectively, the distinct expression patterns and biological functions of ASPM isoforms may illuminate novel molecular mechanisms and prognosticators in PDAC and may pave the way for the development of therapies targeting this novel oncoprotein. © 2019 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Identification of novel genes associated with a poor prognosis in pancreatic ductal adenocarcinoma via a bioinformatics analysis

Pancreatic ductal adenocarcinoma (PDAC) is a class of the commonest malignant carcinomas. The present study aimed to elucidate the potential biomarker and prognostic targets in PDAC. The array data of GSE41368, GSE43795, GSE55643, and GSE41369 were downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) and differentially expressed microRNAs (DEmiRNAs) in PDAC were obtained by using GEO2R, and overlapped DEGs were acquired with Venn Diagrams. Functional enrichment analysis of overlapped DEGs and DEmiRNAs was conducted with Metascape and FunRich, respectively. The protein-protein interaction (PPI) network of overlapped DEGs was constructed by STRING and visualized with Cytoscape. Overall survival (OS) of DEmiRNAs and hub genes were investigated by Kaplan-Meier (KM) plotter (KM plotter). Transcriptional data and correlation analyses among hub genes were verified through GEPIA and Human Protein Atlas (HPA). Additionally, miRNA targets were searched using miRTarBase, then miRNA-DEG regulatory network was visualized with Cytoscape. A total of 32 DEmiRNAs and 150 overlapped DEGs were identified, and Metascape showed that DEGs were significantly enriched in cellular chemical homeostasis and pathways in cancer, while DEmiRNAs were mainly enriched in signal transduction and Glypican pathway. Moreover, seven hub genes with a high degree, namely, V-myc avian myelocytomatosis viral oncogene homolog (MYC), solute carrier family 2 member 1 (SLC2A1), PKM, plasminogen activator, urokinase (PLAU), peroxisome proliferator activated receptor γ (PPARG), MET proto-oncogene, receptor tyrosine kinase (MET), and integrin subunit α 3 (ITGA3), were identified and found to be up-regulated between PDAC and normal tissues. miR-135b, miR-221, miR-21, miR-27a, miR-199b-5p, miR-143, miR-196a, miR-655, miR-455-3p, miR-744 and hub genes predicted poor OS of PDAC. An integrative bioinformatics analysis identified several hub genes that may serve as potential biomarkers or targets for early diagnosis and precision target treatment of PDAC.

A robust 6-mRNA signature for prognosis prediction of pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is one of the most fatal malignancies worldwide. PDAC prognostic and diagnostic biomarkers are still being explored. The aim of this study is to establish a robust molecular signature that can improve the ability to predict PDAC prognosis. 155 overlapping differentially expressed genes between tumor and non-tumor tissues from three Gene Expression Omnibus (GEO) datasets were explored. A least absolute shrinkage and selection operator method (LASSO) Cox regression model was employed for selecting prognostic genes. We developed a 6-mRNA signature that can distinguish high PDAC risk patients from low risk patients with significant differences in overall survival (OS). We further validated this signature prognostic value in three independent cohorts (GEO batch, P < 0.0001; ICGC, P = 0.0036; Fudan, P = 0.029). Furthermore, we found that our signature remained significant in patients with different histologic grade, TNM stage, locations of tumor entity, age and gender. Multivariate cox regression analysis showed that 6-mRNA signature can be an independent prognostic marker in each of the cohorts. Receiver operating characteristic curve (ROC) analysis also showed that our signature possessed a better predictive role of PDAC prognosis. Moreover, the gene set enrichment analysis (GSEA) analysis showed that several tumorigenesis and metastasis related pathways were indeed associated with higher scores of risk. In conclusion, identifying the 6-mRNA signature could provide a valuable classification method to evaluate clinical prognosis and facilitate personalized treatment for PDAC patients. New therapeutic targets may be developed upon the functional analysis of the classifier genes and their related pathways.

Identification and validation of a potent multi-mRNA signature for the prediction of early relapse in hepatocellular carcinoma

Early recurrence of hepatocellular carcinoma (HCC) is implicated in poor patient survival and is the major obstacle to improving prognosis. The current staging systems are insufficient for accurate prediction of early recurrence, suggesting that additional indicators for early recurrence are needed. Here, by analyzing the gene expression profiles of 12 Gene Expression Omnibus data sets (n = 1533), we identified 257 differentially expressed genes between HCC and non-tumor tissues. Least absolute shrinkage and selection operator regression model was used to identify a 24-messenger RNA (mRNA)-based signature in discovery cohort GSE14520. With specific risk score formula, patients were divided into high- and low-risk groups. Recurrence-free survival within 2 years (early-RFS) was significantly different between these two groups in discovery cohort [hazard ratio (HR): 7.954, 95% confidence interval (CI): 4.596–13.767, P < 0.001], internal validation cohort (HR: 8.693, 95% CI: 4.029–18.754, P < 0.001) and external validation cohort (HR: 5.982, 95% CI: 3.414–10.480, P < 0.001). Multivariable and subgroup analyses revealed that the 24-mRNA-based classifier was an independent prognostic factor for predicting early relapse of patients with HCC. We further developed a nomogram integrating the 24-mRNA-based signature and clinicopathological risk factors to predict the early-RFS. The 24-mRNA-signature-integrated nomogram showed good discrimination (concordance index: 0.883, 95% CI: 0.836–0.929) and calibration. Decision curve analysis demonstrated that the 24-mRNA-signature-integrated nomogram was clinically useful. In conclusion, our 24-mRNA signature is a powerful tool for early-relapse prediction and will facilitate individual management of HCC patients.

Genome-Wide Profiling of Prognostic Alternative Splicing Pattern in Pancreatic Cancer

Alternative splicing (AS) has a critical role in tumor progression and prognosis. Our study aimed to investigate pancreatic cancer-specific AS events using RNA-seq data, gaining systematic insights into potential prognostic predictors. We downloaded 10,623 genes with 45,313 pancreatic cancer-specific AS events from the Cancer Genome Atlas (TCGA) and SpliceSeq database. Cox univariate analyses of overall survival suggested there was a remarkable association between 6,711 AS events and overall survival in pancreatic cancer patients (P < 0.05). The area under the curves (AUC) of the receiver operator characteristic curves (ROC) of risk score was 0.89 for final prognostic predictor. Results indicated that AS events of DAZAP1, RBM4, ESRP1, QKI, and SF1 were significantly associated with overall survival. The results of FunRich showed that transcription factors KLF7, GABPA, and SP1 were the most highly related to survival-associated AS genes. Furthermore, using DriverDBv2, we identified 13 driver genes associated with survival-associated AS events, including TP53 and CDC27. Thus, we concluded that the aberrant AS patterns in pancreatic cancer patients might serve as prognostic predictors.