Open AccessTranscriptome profiling reveals liver metastasis-associated genes in pancreatic ductal adenocarcinoma
Author(s) -
Yang Yu,
Zhe Wang,
Dai hai Mo,
Zhen Wang,
Gang Li
Publication year - 2021
Publication title -
mathematical biosciences and engineering
Language(s) - Uncategorized
Resource type - Journals
SCImago Journal Rank - 0.451
H-Index - 45
eISSN - 1551-0018
pISSN - 1547-1063
DOI - 10.3934/mbe.2021088
Subject(s) - metastasis , transcriptome , cancer research , biology , pancreatic cancer , epigenetics , gene expression profiling , epithelial–mesenchymal transition , gene , oncology , cancer , medicine , gene expression , genetics
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.