Mathematical Biosciences and Engineering (Apr 2021)

Transcriptome profiling reveals liver metastasis-associated genes in pancreatic ductal adenocarcinoma

  • Yang Yu,
  • Zhe Wang ,
  • Dai hai Mo,
  • Zhen Wang,
  • Gang Li

DOI
https://doi.org/10.3934/mbe.2021088
Journal volume & issue
Vol. 18, no. 2
pp. 1708 – 1721

Abstract

Read online

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies worldwide. However, the mechanisms underlying the acquisition of the metastatic potential in PDAC has not been well understood. In this study, we identified a total of 154 genes upregulated in primary tissues of PDAC with liver metastasis using the Genome Cancer Atlas (TCGA) and GSE151580 cohorts. The epithelial-mesenchymal transition and glycolysis were enriched by the liver metastasis-related genes, indicating that the liver metastasis-related genes might be functionally relevant to liver metastasis in PDAC. Moreover, we also found that the liver metastasis-related genes were primarily regulated at epigenetic level. Particularly, SFN, a cell cycle checkpoint protein, and KRT19, a marker gene for ductal cells, were predicted to be regulated by multiple methylation sites at the promoter. Clinically, we for the first time defined a liver metastasis score (LMS), which was derived from liver metastasis-related genes, and closely associated with clinical characteristics such as disease type and tumor grade, in PDAC. Furthermore, we also divided the samples into high and low LMS groups using three cohorts with long-term follow-up (TCGA, GSE71729, and E-MTAB-6134), which exhibited significantly different prognostic outcomes across three PDAC cohorts, suggesting that the LMS might be a good indicator for risk stratification in PDAC. Furthermore, we also found that the liver metastasis-related genes were primarily expressed in malignant ductal cells by integrative analysis of the bulk and single-cell gene expression data. Moreover, the malignant ductal cells and M0 macrophages were highly correlated with LMS, indicating that the two cell types might function as tumor-promoting cells in PDAC. In summary, the systematic analysis for the first time defined an LMS score to evaluate the risk of liver metastasis in PDAC and revealed that malignant ductal cells might promote PDAC liver metastasis, which greatly improves our understanding of the liver metastasis-related genes, their underlying mechanisms, and association with prognosis in PDAC.

Keywords