The extracellular matrix (ECM) is a vital component of the tumor microenvironment, which interplays with stromal and tumor cells to stimulate the capacity of cancer cells to proliferate, migrate, invade, and undergo angiogenesis. Nevertheless, the crucial functions of ECM-related genes (ECMGs) in pancreatic adenocarcinoma (PAAD) have not been systematically evaluated. Hence, a comprehensive evaluation of the ECMGs is required in pan-cancer, especially in PAAD. First, a pan-cancer overview of ECMGs was explored through the integration of expression profiles, prognostic values, mutation information, methylation levels, and pathway-regulation relationships. Seven ECMGs (i.e. LAMB3, LAMA3, ITGB6, ITGB4, ITGA2, LAMC2, and COL11A1) were identified to be hub genes of PAAD, which were obviously up-regulated in PAAD and considerably linked to tumor stage as well as prognosis. Subsequently, patients with PAAD were divided into 3 clusters premised on ECMG expression and ECM scores. Cluster 2 was the subtype with the best prognosis accompanied by the lowest ECM scores, further verifying ECM’s significant contribution to the pathophysiological processes of PAAD. Significant differences were observed for oncogene and tumor suppressor gene expression, immune microenvironment, and chemotherapy sensitivity across three ECM subtypes. After applying a variety of bioinformatics methods, a novel and robust ECM-associated mRNA-lncRNA-based prognostic panel (ECM-APP) was developed and validated for accurately predicting clinical outcomes of patients with PAAD. Patients with PAAD were randomly categorized into the train, internal validation, and external validation cohorts; meanwhile, each patient was allocated into high-risk (unfavorable prognosis) and low-risk (favorable prognosis) populations premised on the expression traits of ECM-related mRNAs and lncRNAs. The discrepancy in the tumor mutation burden and immune microenvironment might be responsible for the difference in prognoses across the high-risk and low-risk populations. Overall, our findings identified and validated seven ECMGs remarkably linked to the onset and progression of PAAD. ECM-based molecular classification and prognostic panel aid in the prognostic assessment and personalized intervention of patients with PAAD.

The worldwide prevalence of pancreatic cancer has been rising in recent decades, and its prognosis has not improved much. The imbalance of substance and energy metabolism in tumour cells is among the primary causes of tumour formation and occurrence, which is often controlled by the neuroendocrine system. We applied Cox and LASSO regression analysis to develop a neuroendocrine regulation- and metabolism-related prognostic risk score model with three genes (GSK3B, IL18 and VEGFA) for pancreatic cancer. TCGA dataset served as the training and internal validation sets, and GSE28735, GSE62452 and GSE57495 were designated as external validation sets. Patients classified as the low-risk population (category, group) exhibited considerably improved survival duration in contrast with those classified as the high-risk population, as determined by the Kaplan-Meier curve. Then, we combined all the samples, and divided them into three clusters using unsupervised clustering analysis. Unsupervised clustering, t-distributed stochastic neighbor embedding (t-SNE), and principal component analysis (PCA) were further utilized to demonstrate the reliability of the prognostic model. Moreover, the risk score was shown to independently function as a predictor of pancreatic cancer in both univariate and multivariate Cox regression analyses. The results of gene set enrichment analysis (GSEA) illustrated that the low-risk population was predominantly enriched in immune-associated pathways. “ESTIMATE” algorithm, single-sample GSEA (ssGSEA) and the Tumor Immune Estimation Resource (TIMER) database showed immune infiltration ratings were enhanced in the low-risk category in contrast with the high-risk group. Tumour immune dysfunction and exclusion (TIDE) database predicted that immunotherapy for pancreatic cancer may be more successful in the high-risk than in the low-risk population. Mutation analysis illustrated a positive link between the tumour mutation burden and risk score. Drug sensitivity analysis identified 44 sensitive drugs in the high- and low-risk population. GSK3B expression was negatively correlated with Oxaliplatin, and IL18 expression was negatively correlated with Paclitaxel. Lastly, we analyzed and verified gene expression at RNA and protein levels based on GENPIA platform, HPA database and quantitative real-time PCR. In short, we developed a neuroendocrine regulation- and metabolism-associated prognostic model for pancreatic cancer that takes into account the immunological microenvironment and drug sensitivity.

Ubiquitination-related genes (URGs) exerted a crucial part in a variety of human disease disorders; however, their association with pancreatic adenocarcinoma (PAAD) had yet to be clearly described. We aimed to comprehensively characterize the contributions of URGs in PAAD through in silico analysis and experimental validation, and then identified a robust mRNA-lncRNA-based molecular prognostic panel for patients with PAAD using bulk RNA-sequencing and single-cell RNA-sequencing data. Initially, we collected the multi-omics data from TCGA platform to depict a comprehensive landscape of URGs in pan-cancer. Furthermore, we were accurate to PAAD for in-depth analysis. Significant differences of the activation of ubiquitination pathways and the expression of URGs were detected between normal and malignant cells. Unsupervised hierarchical clustering determined two PAAD subtypes with distinct clinical outcomes, ubiquitination pathway activities, immune microenvironment, and functional annotation characteristics. The expression profiles of ubiquitination-associated mRNAs and lncRNAs in the training and validation datasets were utilized to develop and verify a novel ubiquitination-related mRNA-lncRNA prognostic panel, which had a satisfied prediction efficiency. Our ubiquitination-associated model could function as an effective prognostic index and outperformed four other recognized panels in evaluating PAAD patients’ survival status. Tumor immune microenvironment, mutation burden, and chemotherapy response were intensively explored to demonstrate the underlying mechanism of prognostic difference according to our panel. Our findings also revealed that FTI-277, a farnesyltransferase inhibitor, had a better curative effect in high-risk patients, while MK-2206, an Akt allosteric inhibitor, had a superior therapeutic effect in low-risk patients. The real-time PCR results uncovered the RNA expression of AC005062.1 in all the three PAAD cell lines was elevated several thousandfold. In conclusion, our URGs-based classification panel could be triumphantly served as a prediction tool for survival evaluation in patients with PAAD, and the genes in this panel could be developed as a potential target in PAAD therapy.