Research article Special Issues

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

  • 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.

    Citation: Yang Yu, Zhe Wang, Dai hai Mo, Zhen Wang, Gang Li. Transcriptome profiling reveals liver metastasis-associated genes in pancreatic ductal adenocarcinoma[J]. Mathematical Biosciences and Engineering, 2021, 18(2): 1708-1721. doi: 10.3934/mbe.2021088

    Related Papers:

    [1] Yan Li, Yuzhang Zhu, Guiping Dai, Dongjuan Wu, Zhenzhen Gao, Lei Zhang, Yaohua Fan . Screening and validating the core biomarkers in patients with pancreatic ductal adenocarcinoma. Mathematical Biosciences and Engineering, 2020, 17(1): 910-927. doi: 10.3934/mbe.2020048
    [2] Jinqi He, Wenjing Zhang, Faxiang Li, Yan Yu . Development of metastasis-associated seven gene signature for predicting lung adenocarcinoma prognosis using single-cell RNA sequencing data. Mathematical Biosciences and Engineering, 2021, 18(5): 5959-5977. doi: 10.3934/mbe.2021298
    [3] Xinwang Yan, Xiaowen Zhao, Qing Yan, Ye Wang, Chunling Zhang . Analysis of the role of METTL5 as a hub gene in lung adenocarcinoma based on a weighted gene co-expression network. Mathematical Biosciences and Engineering, 2021, 18(5): 6608-6619. doi: 10.3934/mbe.2021327
    [4] Marcello Delitala, Tommaso Lorenzi . Emergence of spatial patterns in a mathematical model for the co-culture dynamics of epithelial-like and mesenchymal-like cells. Mathematical Biosciences and Engineering, 2017, 14(1): 79-93. doi: 10.3934/mbe.2017006
    [5] Yuexian Li, Zhou Su, Biwei Wei, Mengbin Qin, Zhihai Liang . Bioinformatics analysis identified MMP14 and COL12A1 as immune-related biomarkers associated with pancreatic adenocarcinoma prognosis. Mathematical Biosciences and Engineering, 2021, 18(5): 5921-5942. doi: 10.3934/mbe.2021296
    [6] Yong Ding, Jian-Hong Liu . The signature lncRNAs associated with the lung adenocarcinoma patients prognosis. Mathematical Biosciences and Engineering, 2020, 17(2): 1593-1603. doi: 10.3934/mbe.2020083
    [7] Ji-Ming Wu, Wang-Ren Qiu, Zi Liu, Zhao-Chun Xu, Shou-Hua Zhang . Integrative approach for classifying male tumors based on DNA methylation 450K data. Mathematical Biosciences and Engineering, 2023, 20(11): 19133-19151. doi: 10.3934/mbe.2023845
    [8] Huili Yang, Wangren Qiu, Zi Liu . Anoikis-related mRNA-lncRNA and DNA methylation profiles for overall survival prediction in breast cancer patients. Mathematical Biosciences and Engineering, 2024, 21(1): 1590-1609. doi: 10.3934/mbe.2024069
    [9] Yi Shi, Xiaoqian Huang, Zhaolan Du, Jianjun Tan . Analysis of single-cell RNA-sequencing data identifies a hypoxic tumor subpopulation associated with poor prognosis in triple-negative breast cancer. Mathematical Biosciences and Engineering, 2022, 19(6): 5793-5812. doi: 10.3934/mbe.2022271
    [10] Mengyang Han, Xiaoli Wang, Yaqi Li, Jianjun Tan, Chunhua Li, Wang Sheng . Identification of coagulation-associated subtypes of lung adenocarcinoma and establishment of prognostic models. Mathematical Biosciences and Engineering, 2023, 20(6): 10626-10658. doi: 10.3934/mbe.2023470
  • 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.



    Pancreatic cancer is a malignancy with a very poor prognosis, as the number of its new cases almost equals to related deaths [1]. The major type of pancreatic tumors is pancreatic ductal adenocarcinoma (PDAC), where metastasis is present in most cases, leading to an extremely low 5-year survival rate [2]. PDAC samples can be classified into four major molecular subtypes including squamous, pancreatic progenitor, immunogenic, and aberrantly differentiated endocrine exocrine (ADEX) subtypes, and squamous subtype had the worse prognosis [3]. Meanwhile, it has been demonstrated that PDAC patient with isolated liver metastases exhibited worse overall survival when compared to those with other sites of metastases [4]. Thus, for patients with PDAC, it is of great importance to investigate molecular mechanism behind liver-metastatic behaviors and identify markers to assess the risks of metastasis.

    So far, several researches have focused on demonstrating the genomic changes behind metastatic PDAC. Unique gene expression profiles were uncovered between liver and peritoneal metastatic PDAC cell lines, suggesting that different mechanisms were involved when site-specific metastases were present [5]. Aberrant expression of Zeb1, EphA2, DPC4 and MKK4 was found to be associated with metastatic pancreatic cancers [6,7], and dysregulated PRDM14 expression was confirmed to affect liver metastasis in mice [8]. Notably, Zeb1 is a key transcription factor promoting epithelial-mesenchymal transition (EMT) in pancreatic cancer [9]. Moreover, metastatic PDAC is often accompanied with an immunosuppressive tumor microenvironment (TME), and evidences suggest that stemness and differentiation of pancreatic ductal epithelial cells were associated with inflammatory cytokines in the hepatic microenvironment [10]. For example, interleukin (IL)-1α, a cytokine that mostly secreted by macrophages, were found to induce metastasis in PDAC, while a recent study has reported an association between macrophage-derived granulin and CD8+ T cell infiltration in PDAC metastatic tumors [11]. Pharmacological depletion of macrophages was deemed effective in reducing metastatic behavior in mice with pancreatic cancer [12]. However, molecular mechanisms possibly underlying altered macrophage secretion and infiltration are seldom addressed, thus it is of great necessity to explore significantly dysregulated genes in varied types of immune cells in the metastatic TME of PDAC patients [13]. In the present study, we focused on exploring the RNA signatures of liver metastasis in PDAC, anticipated to identify a set of liver metastasis-related genes, evaluate their prognostic values, and reveal key cell types involved in liver metastasis of PDAC.

    We collected four gene expression datasets including The Cancer Genome Atlas [14] (TCGA), GSE151580, GSE71729 [15] and E-MTAB-6134 [16] from UCSC Xena (https://xena.ucsc.edu/), Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/), and ArrayExpress (https://www.ebi.ac.uk/arrayexpress/), for the systematic data analysis, respectively. The detailed clinical characteristics of the PDAC samples were summarized in Table 1. The first two gene expression datasets were used to identify liver metastasis-related genes. The RNA-seq data were normalized to FPKM or TPM, and the microarray data were pre-normalized by previous studies. Moreover, we also collected the cell types and their specific marker genes from previous study [17].

    Table 1.  The clinical characteristics of the PDAC samples used for data analysis.
    Clinical variables TCGA (n = 182) GSE71729 (n = 145) E-MTAB-6134 (n = 309) GSE151580 (n = 13)
    TNM stage
    I 11.8% (21) 6.8% (21) 0% (0)
    II 84.3% (150) 93.2% (288) 0% (0)
    III 1.7% (3) 0% (0) 0% (0)
    IV 2.2% (4) 0% (0) 100% (13)
    Gender
    Male 54.7% (99) 57.9% (179)
    Female 45.3% (82) 42.1% (130)
    Overall survival
    Dead 51.9% (94) 67.2% (84) 62.8% (181)
    Alive 48.1% (87) 32.8% (41) 37.2% (107)
    Liver metastasis
    Yes 11.5% (21) 100% (13)
    No 88.5% (161) 0% (0)

     | Show Table
    DownLoad: CSV

    The two-sample comparison was conducted to identify differentially expressed genes, and the differential gene expression levels were test by t test and fold changes. The thresholds of P-value and log2 (fold change) were selected at 0.01 and 0.25.

    The Gene set overrepresentation enrichment analysis was used to test the degree of overlapping genes between two gene sets. The gene set enrichment analysis was employed to test whether the gene set of interest was expressed at the lowest or highest order of all the genes ranked by a specified statistic. The two analyses were implemented in R clusterProfiler [18] package.

    The liver metastasis score was calculated using single-sample gene set enrichment analysis (ssGSEA) [19], which defines an enrichment score that represents the degree of absolute enrichment of a gene set in each sample within a given data set. The enrichment score, a variant of Kolmogorov-Smirnov statistic, was proposed by the previous study [20]. The ssGSEA was implemented in R GSVA [21] package.

    The overall survival time and status, and the liver metastasis score were used as the response and predictor variables for the Cox proportional hazard regression model. The Cox model was fitted using R survival [22] package (https://cran.r-project.org/web/packages/survival/index.html), and visualized by R survminer package (https://cran.r-project.org/web/packages/survminer/index.html). The log-rank test was employed to test the difference of survival time between the two groups.

    The statistical analysis was implemented in R-4.0.2. The two-sample comparisons were performed by Wilcoxon rank-sum test or student-t test. Multiple-sample comparisons were performed by analysis of variance (ANOVA) or Kruskal-Wallis test. P-value of 0.05 was indicated as statistical significance.

    To identify the liver metastasis-related genes in PDAC, we conducted a systematic analysis of two gene expression datasets. Specifically, we compared the primary tissues of liver-metastatic PDAC patients with those of PDAC without liver metastasis using TCGA cohort [14] (n = 182). Moreover, the primary tissues of liver-metastatic PDAC patients were also compared with their adjacent normal tissues using GSE151580 cohort (n = 19). Totally, we identified 154 genes upregulated in primary tissues of liver-metastatic PDAC patients as compared with both primary tissues of PDAC patients without liver metastasis and adjacent normal tissues (P-value < 0.01 and log2 (fold change) > 0.25, Supplementary Table S1).

    As shown in Figure 1A and 1B, the 154 genes had significantly different expression profiles between those tissues, and were termed as liver metastasis-related genes. Furthermore, the gene set enrichment analysis was conducted to characterize the biological function of those liver metastasis-related genes, and epithelia mesenchymal transition (EMT) was significantly enriched by those genes (Figure 1C, Supplementary Table S2), indicating that EMT was significantly associated with liver metastasis in PDAC. In addition, estrogen response, p53, and glycolysis pathways were also enriched by those genes, further suggesting their potential roles in the liver metastasis of PDAC.

    Figure 1.  The genes associated with liver metastasis in pancreatic ductal adenocarcinoma (PDAC). The upregulated genes in primary PDAC tissues with liver metastasis, which was compared with adjacent normal tissues (A) and primary tissues of PDAC patients with liver metastasis (B). (C) The pathways enriched by the liver metastasis-associated genes in PDAC. The node size and color represent the number of genes within the pathway and adjusted P-value for enrichment. The gene ratio was calculated by dividing the total number of genes by the number of upregulated genes within the pathway.

    To further investigate whether the liver metastasis-related genes were regulated at epigenetic level, we conducted correlation analysis between the DNA methylation level within promoter regions and their corresponding target genes using TCGA cohort. Among the liver metastasis-related genes, 18 were predicted to be epigenetically regulated by their promoter DNA methylation at a stringent threshold (Figure 2A, correlation test, FDR < 0.05, n = 182). Notably, SFN (Stratifin), a gene encoding a cell cycle checkpoint protein [23], was reversely correlated with the DNA methylation levels of 8 sites, including cg06720467, cg07786675, cg11348165, cg13374701, cg13466284, cg14825555, cg17330303, and cg21950166, within the promoter (Figure 2B). Moreover, KRT19 (Keratin-19), a well-known marker for ductal cell [24], was also negatively correlated with its DNA methylation levels of 6 sites. These results suggested that DNA methylation acted as a pivotal role in the regulation of the liver metastasis-related genes.

    Figure 2.  The regulation of genes associated with liver metastasis by DNA methylation. (A) The liver metastasis-associated genes significantly regulated by DNA methylation. The color bars on the top represent the genes regulated by DNA methylation, and the width of the bars represent the number of CpG sites. The green and red colors in the heatmap represent the negative and positive correlation, respectively. The DNA methylation and gene expression of two well-recognized metastatic genes, SFN (B) and KRT19 (C). The x and y axis represent the beta value of methylation (probe-level) and log2-transformed FPKM (Fragment Per Kilo Million).

    To reveal the association of the liver metastasis-related genes with the clinical characteristics of PDAC, we derived a liver metastasis score (LMS) based on the liver metastasis-related genes for TCGA PDAC patients. Specifically, the LMS of each PDAC patient was estimated by single-sample gene set enrichment analysis (ssGSEA). The PDAC patients with liver metastasis had significantly higher LMS than those without liver metastasis (Figure 3A, t test, P-value < 0.0001). Moreover, the LMS was observed lower in patients with stage I as compared with those with stage II (Figure 3B, P-value < 0.05). Furthermore, the comparative analysis of the disease types and tumor grades revealed that the ductal and lobular type had a higher liver metastasis score than adenomas and adenocarcinoma (Figure 3C, P-value < 0.0001), and higher LMS was observed in PDAC patients with higher grade (Figure 3D, P-value < 0.05), suggesting that patients with ductal and lobular type or higher grade might have a higher risk of liver metastasis. In addition, there were no significant associations of LMS with tumor sites and gender (Figure 3E, F). These results demonstrated that the LMS derived from liver metastasis-related genes was closely associated with the clinical characteristics of disease type and tumor grade in PDAC.

    Figure 3.  The association of liver metastasis score (LMS) with clinical factors. The differential LMS levels of PDAC patients with different metastatic status, TNM stages, pathological classifications, grades, anatomical locations, and genders are displayed in (A–F).

    To evaluate the prognostic significance of the liver metastasis-related genes in PDAC, we built a Cox proportional hazard regression model based on the liver metastasis score (LMS). The samples in The Cancer Genome Atlas (TCGA) cohort were divided into high and low LMS groups by the LMS median. Specifically, the high LMS group showed significantly shorter overall survival than low LMS group (Figure 4A, log-rank test, P-value < 0.05). To further confirm the prognostic value of LMS, we also divided the samples from two public datasets, termed as GSE71729 [15] (n = 145) and E-MTAB-6134 [16] (n = 309), into high and low LMS groups, respectively. Consistently, the prognostic differences were also observed in the two datasets (Figure 4B, C, Table 2, P-value < 0.05). Collectively, these results indicated that the LMS might be a prognostic predictor for PDAC.

    Figure 4.  The association of liver metastasis score (LMS) with PDAC prognosis. The association of liver metastasis score (LMS) with PDAC prognosis is evaluated in three independent cohorts, TCGA (A), E-MTAB-6134 (B), and GSE71729 (C). The horizontal and vertical dashed lines represent the 50% survival rate and corresponding survival time. The number at risk means the number of samples survived before that time point. The number of censoring means the number of samples whose follow-up is censored at that time point.
    Table 2.  The hazard ratio and coefficient for LMS in the Cox model.
    Cohort log2 HR HR log2 HR standard error Z P-value
    TCGA −0.56 0.57 0.21 −2.66 7.92 × 10−3
    GSE71729 −0.60 0.55 0.22 −2.65 8.01 × 10−3
    E-MTAB-6134 −0.66 0.52 0.15 −4.37 1.27 × 10−5

     | Show Table
    DownLoad: CSV

    As the tumor metastasis was regulated by the cell types infiltrating into the tumor microenvironment and some specific tumor cells, we then investigated whether some specific cell types were responsible for liver metastasis in PDAC. Firstly, we collected the marker genes of cell types in PDAC tissues from earlier study [17]. The overrepresentation test of the liver metastasis-related genes revealed that those genes were highly enriched in the marker genes of malignant ductal cell (Figure 5A). Specifically, some metastasis-related genes such as MET (MET Proto-Oncogene), KRT19 (Keratin-19), KRT7 (Keratin-7), AGRN (Agrin), SDC4 (Syndecan-4), and SFN (Stratifin) were found to be specifically expressed in malignant ductal cell, indicating that the liver metastasis-related genes might be expressed in malignant ductal cells, which might be responsible for the liver metastasis in PDAC. Moreover, we also conducted correlation analysis between LMS and gene expression levels. Gene set enrichment analysis (GSEA) revealed that most of the marker genes in malignant ductal cell were positively correlated with LMS (Figure 5B, FDR < 0.05), further suggesting that malignant ductal cells were responsible for liver metastasis via expressing the liver metastasis-related genes. Furthermore, marker genes of M0 macrophage were also highly positively correlated with LMS (Figure 5C, FDR < 0.05). In addition, as motivated by previous study [25], we also estimated the relative abundance of macrophage M0 in PDAC tissues using gene expression and macrophage M0-specific genes by ssGSEA method. Consistently, the LMS was highly positively correlated with the relative activity of macrophage M0 (Figure 5D, Spearman correlation: 0.784). Therefore, the positive correlation between LMS and malignant ductal cell or macrophage M0 specific genes indicated that the two cell types might have the potential to promote liver metastasis of PDAC.

    Figure 5.  The association of metastasis-associated genes with malignant ductal cells and macrophage M0. (A) The shared number of genes between metastasis-associated genes and malignant ductal cell-specific genes. (B) Positive correlation of malignant ductal cell-specific genes with liver metastasis score (LMS). (C) Positive correlation of macrophage M0-specific genes with liver metastasis score (LMS). The genes were ranked by the correlation coefficients by decreasing order. The black lines in the middle represent the cell type specific genes. These cell-type specific genes were mainly clustered within the gene sets with high positive correlation with LMS (on the left). (D) The correlation between LMS and the relative abundance of macrophage M0.

    Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies worldwide [26]. Liver metastasis accounts for the high mortality rate and confers poor prognosis in PDAC [27]. Therefore, a better understanding of the mechanisms underlying the acquisition of the metastatic potential in PDAC is highly desirable.

    In this study, we conducted a systematic analysis of gene expression data to identify liver metastasis-related genes in PDAC. Specifically, a total of 154 genes were identified to be upregulated in primary tissues of PDAC with liver metastasis. These liver metastasis-related genes were involved in signaling pathways such as EMT, estrogen response, p53, and glycolysis. EMT, a developmental program that enables stationary epithelial cells to gain the ability to migrate and invade, has been implicated in carcinogenesis and cancer metastasis by enhancing mobility and invasion [28]. Consistently, EMT was also found as the key players in PDAC metastasis [29]. Particularly, glycolysis was also identified as a key pathway involved in liver metastasis of PDAC. The metabolic reprogramming from oxidative phosphorylation to glycolysis in pancreatic cancer cells, has been found to enhance the invasion-metastasis cascade by promoting EMT, angiogenesis and metastatic colonization of distant organs [30]. Combined with the result in this study, the liver metastasis-related genes were functionally relevant to liver metastasis in PDAC.

    The integrative analysis of DNA methylation and gene expression data revealed 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. It should be noted that the two genes have been found as novel biomarkers of poor prognosis in PDAC [31,32]. We thus speculated that the overexpression of SFN might be associated with uncontrolled cell cycle progression and stemness maintenance, while the high expression of KRT19 gave us a hint that ductal cells might play key roles in liver metastasis in PDAC. In addition, some other DNA methylation-associated genes were also reported to be associated with tumor initiation or progression. For example, TACSTD2 (tumor associated calcium signal transducer 2) encoded a carcinoma-associated antigen [33], suggesting that it might be a potential diagnostic marker in PDAC.

    Clinically, we found that liver metastasis score (LMS), derived from liver metastasis-related genes, was found to be closely associated with clinical characteristics associated with prognostic outcomes [34], such as disease type and tumor grade, in PDAC. Furthermore, we also divided the samples from TCGA, GSE71729 and E-MTAB-6134 cohorts into high and low LMS groups, and the two groups exhibited significantly different prognostic outcomes across the three cohorts, suggesting that the LMS might be a promising biomarker for risk stratification in PDAC. However, as the liver metastasis-related genes used to derive LMS were mainly enriched in glycolysis and EMT pathway, which were altered in many other cancers, the LMS might not be the PDAC-specific prognostic marker.

    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. An earlier study [17] reported that the malignant ductal cells might be transited from ductal cells with abnormal gene expression profiles to malignant ductal cells by trajectory analysis. The malignant ductal cells were characterized to have highly proliferative and migratory subpopulations [17]. Moreover, the high correlation between LMS and M0 macrophage indicated that M0 macrophage might promote liver metastasis of PDAC. The macrophage cell subpopulations in tumors commonly refer to M2 macrophage, exhibiting anti-inflammatory and pro-tumoral effects [35], however, the role of unpolarized macrophage (M0) in tumor tissues was rarely reported. In accordance with our finding, low expression of M0 macrophage is associated with better clinical prognosis in bladder cancer patients [36]. Collectively, the malignant ductal cells and M0 macrophage might function as tumor-promoting cells in PDAC.

    However, this study still has some limitations. First, the biological function of some novel liver metastasis-related genes has not been validated by experiments. Second, the clinical significance of malignant ductal cells needs to be evaluated by clinical specimens. In addition, the protein expression levels of some important liver metastasis-related genes need to be evaluated. In summary, the systematic analysis greatly improves our understanding of the critical cell types and genes expressed by these cell types in the process of liver metastasis in PDAC.

    This project is supported by the Lead project of Western medicine of the Shanghai Science Committee (No. 16411967200), the PLA's youth training program of the medical science and technology (No. 16QNP095), and National Natural Science Foundation of China (No. 81871992).

    The authors declare that they have no conflicts of interest.



    [1] R. L. Siegel, K. D. Miller, A. Jemal, Cancer statistics, 2020, CA Cancer J. Clin., 70 (2020), 7–30. doi: 10.3322/caac.21590
    [2] E. Giovannetti, C. L. van der Borden, A. E. Frampton, A. Ali, O. Firuzi, G. J. Peters, Never let it go: Stopping key mechanisms underlying metastasis to fight pancreatic cancer, Semin. Cancer Biol., 44 (2017), 43–59. doi: 10.1016/j.semcancer.2017.04.006
    [3] C. Torres, P. J. Grippo, Pancreatic cancer subtypes: a roadmap for precision medicine, Ann. Med., 50 (2018), 277–287. doi: 10.1080/07853890.2018.1453168
    [4] H. Oweira, U. Petrausch, D. Helbling, J. Schmidt, M. Mannhart, A. Mehrabi, et al., Prognostic value of site-specific metastases in pancreatic adenocarcinoma: A Surveillance Epidemiology and End Results database analysis, World J. Gastroenterol., 23 (2017), 1872–1880. doi: 10.3748/wjg.v23.i10.1872
    [5] K. Ohno, F. Hata, H. Nishimori, T. Yasoshima, Y. Yanai, K. Sogahata, et al., Metastatic-associated biological properties and differential gene expression profiles in established highly liver and peritoneal metastatic cell lines of human pancreatic cancer, J. Exp. Clin. Cancer Res., 22 (2003), 623–631.
    [6] S. Yachida, C. A. Iacobuzio-Donahue, The pathology and genetics of metastatic pancreatic cancer, Arch. Pathol. Lab. Med., 133 (2009), 413–422.
    [7] A. M. Krebs, J. Mitschke, M. L. Losada, O. Schmalhofer, M. Boerries, H. Busch, et al., The EMT-activator Zeb1 is a key factor for cell plasticity and promotes metastasis in pancreatic cancer, Nat. Cell Biol., 19 (2017), 518–529. doi: 10.1038/ncb3513
    [8] C. Moriya, H. Taniguchi, K. Miyata, N. Nishiyama, K. Kataoka, K. Imai, Inhibition of PRDM14 expression in pancreatic cancer suppresses cancer stem-like properties and liver metastasis in mice, Carcinogenesis, 38 (2017), 638–648. doi: 10.1093/carcin/bgx040
    [9] J. A. Galvan, I. Zlobec, M. Wartenberg, A. Lugli, B. Gloor, A. Perren, et al., Expression of E-cadherin repressors SNAIL, ZEB1 and ZEB2 by tumour and stromal cells influences tumour-budding phenotype and suggests heterogeneity of stromal cells in pancreatic cancer, Br. J. Cancer, 112 (2015), 1944–1950. doi: 10.1038/bjc.2015.177
    [10] H. Knaack, L. Lenk, L. M. Philipp, L. Miarka, S. Rahn, F. Viol, et al., Liver metastasis of pancreatic cancer: the hepatic microenvironment impacts differentiation and self-renewal capacity of pancreatic ductal epithelial cells, Oncotarget, 9 (2018), 31771–31786. doi: 10.18632/oncotarget.25884
    [11] V. Quaranta, C. Rainer, S. R. Nielsen, M. L. Raymant, M. S. Ahmed, D. D. Engle, et al., Macrophage-Derived Granulin Drives Resistance to Immune Checkpoint Inhibition in Metastatic Pancreatic Cancer, Cancer Res., 78 (2018), 4253–4269. doi: 10.1158/0008-5472.CAN-17-3876
    [12] H. Griesmann, C. Drexel, N. Milosevic, B. Sipos, J. Rosendahl, T. M. Gress, et al., Pharmacological macrophage inhibition decreases metastasis formation in a genetic model of pancreatic cancer, Gut, 66 (2017), 1278–1285. doi: 10.1136/gutjnl-2015-310049
    [13] F. Runa, S. Hamalian, K. Meade, P. Shisgal, P. C. Gray, J. A. Kelber, Tumor microenvironment heterogeneity: challenges and opportunities, Curr. Mol. Biol. Rep., 3 (2017), 218–229. doi: 10.1007/s40610-017-0073-7
    [14] Cancer Genome Atlas Research Network, Integrated Genomic Characterization of Pancreatic Ductal Adenocarcinoma, Cancer Cell, 32 (2017), 185–203.
    [15] R. A. Moffitt, R. Marayati, E. L. Flate, K. E. Volmar, S. G. Loeza, K. A. Hoadley, et al., Virtual microdissection identifies distinct tumor- and stroma-specific subtypes of pancreatic ductal adenocarcinoma, Nat. Genet., 47 (2015), 1168–1178. doi: 10.1038/ng.3398
    [16] F. Puleo, R. Nicolle, Y. Blum, J. Cros, L. Marisa, P. Demetter, et al., Stratification of Pancreatic Ductal Adenocarcinomas Based on Tumor and Microenvironment Features, Gastroenterology, 155 (2018), 1999–2013. doi: 10.1053/j.gastro.2018.08.033
    [17] J. Peng, B. F. Sun, C. Y. Chen, J. Y. Zhou, Y. S. Chen, H. Chen, et al., Single-cell RNA-seq highlights intra-tumoral heterogeneity and malignant progression in pancreatic ductal adenocarcinoma, Cell Res., 29 (2019), 725–738. doi: 10.1038/s41422-019-0195-y
    [18] G. Yu, L. G. Wang, Y. Han, Q. Y. He, ClusterProfiler: an R package for comparing biological themes among gene clusters, OMICS, 16 (2012), 284–287. doi: 10.1089/omi.2011.0118
    [19] D. A. Barbie, P. Tamayo, J. S. Boehm, S. Y. Kim, S. E. Moody, I. F. Dunn, et al., Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1, Nature, 462 (2009), 108–112. doi: 10.1038/nature08460
    [20] A. Subramanian, P. Tamayo, V. K. Mootha, S. Mukherjee, B. L. Ebert, M. A. Gillette, et al., Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles, Proc. Natl. Acad. Sci. U. S. A., 102 (2005), 15545–15550. doi: 10.1073/pnas.0506580102
    [21] S. Hanzelmann, R. Castelo, J. Guinney, GSVA: gene set variation analysis for microarray and RNA-seq data, BMC Bioinf., 14 (2013), 7. doi: 10.1186/1471-2105-14-7
    [22] T. M. Therneau, T. Lumley, A. Elizabeth, C. Cynthia, Survival: Survival Analysis, 2020. Available from: https://cran.r-project.org/web/packages/survival/index.html.
    [23] L. Phan, P. C. Chou, G. Velazquez-Torres, I. Samudio, K. Parreno. Y. Huang, et al., The cell cycle regulator 14-3-3sigma opposes and reverses cancer metabolic reprogramming, Nat. Commun., 6 (2015), 7530. doi: 10.1038/ncomms8530
    [24] C. A. Rimland, S. G. Tilson, C. M. Morell, R. A. Tomaz, W. Y. Lu, S. E. Adams, et al., Regional differences in human biliary tissues and corresponding in vitro derived organoids, Hepatology, 2020 (2020).
    [25] Y. Senbabaoglu, R. S. Gejman, A. G. Winer, M. Liu, E. M. Van Allen, G. de Velasco, et al., Tumor immune microenvironment characterization in clear cell renal cell carcinoma identifies prognostic and immunotherapeutically relevant messenger RNA signatures, Genome Biol., 17 (2016), 231. doi: 10.1186/s13059-016-1092-z
    [26] F. Fu, S. Liu, S. Zeng, H. Shen, The critical roles of activated stellate cells-mediated paracrine signaling, metabolism and onco-immunology in pancreatic ductal adenocarcinoma, Mol. Cancer, 17 (2018), 62. doi: 10.1186/s12943-018-0815-z
    [27] T. Y. S. Le Large, M. F. Bijlsma, G. Kazemier, H. W. M. van Laarhoven, E. Giovannetti, C. R. Jimenez, Key biological processes driving metastatic spread of pancreatic cancer as identified by multi-omics studies, Semin. Cancer Biol., 44 (2017), 153–169. doi: 10.1016/j.semcancer.2017.03.008
    [28] K. T. Yeung, J. Yang, Epithelial-mesenchymal transition in tumor metastasis, Mol. Oncol., 11 (2017), 28–39. doi: 10.1002/1878-0261.12017
    [29] E. Rodriguez-Aznar, L. Wiesmuller, B. Sainz, P. C. Hermann, EMT and Stemness-Key Players in Pancreatic Cancer Stem Cells, Cancers (Basel), 11 (2019), 1136. doi: 10.3390/cancers11081136
    [30] J. Yang, B. Ren, G. Yang, H. Wang, G. Chen, L. You, et al., The enhancement of glycolysis regulates pancreatic cancer metastasis, Cell. Mol. Life Sci., 77 (2020), 305–321. doi: 10.1007/s00018-019-03278-z
    [31] F. Robin, G. Angenard, L. Cano, L. Courtin-Tanguy, E. Gaignard, Z. E. Khene, et al., Molecular profiling of stroma highlights stratifin as a novel biomarker of poor prognosis in pancreatic ductal adenocarcinoma, Br. J. Cancer, 123 (2020), 72–80. doi: 10.1038/s41416-020-0863-1
    [32] H. Yao, Z. Yang, Z. Liu, X. Miao, L. Yang, D. Li, et al., Glypican-3 and KRT19 are markers associating with metastasis and poor prognosis of pancreatic ductal adenocarcinoma, Cancer Biomark, 17 (2016), 397–404.
    [33] M. Fornaro, R. Dell'Arciprete, M. Stella, C. Bucci, M. Nutini, M. G. Capri, et al., Cloning of the gene encoding Trop-2, a cell-surface glycoprotein expressed by human carcinomas, Int. J. Cancer, 62 (1995), 610–618. doi: 10.1002/ijc.2910620520
    [34] M. Tanaka, A. L. Mihaljevic, P. Probst, M. Heckler, U. Klaiber, U. Heger, et al., Meta-analysis of recurrence pattern after resection for pancreatic cancer, Br. J. Surg., 106 (2019), 1590–1601. doi: 10.1002/bjs.11295
    [35] M. Najafi, N. Hashemi Goradel, B. Farhood, E. Salehi, M. S. Nashtaei, N. Khanlarkhani, et al., Macrophage polarity in cancer: A review, J. Cell. Biochem., 120 (2019), 2756–2765. doi: 10.1002/jcb.27646
    [36] W. Li, J. Zeng, B. Luo, Y. Mao, Y. Liang, W. Zhao, et al., High expression of activated CD4(+) memory T cells and CD8(+) T cells and low expression of M0 macrophage are associated with better clinical prognosis in bladder cancer patients, Xi Bao Yu Fen Zi Mian Yi Xue Za Zhi, 36 (2020), 97–103.
  • mbe-18-02-088-s001.xlsx
    mbe-18-02-088-s002.xlsx
  • This article has been cited by:

    1. Tina Draškovič, Nina Zidar, Nina Hauptman, Circulating Tumor DNA Methylation Biomarkers for Characterization and Determination of the Cancer Origin in Malignant Liver Tumors, 2023, 15, 2072-6694, 859, 10.3390/cancers15030859
    2. Tao Liu, Jian Chen, An-an Liu, Long Chen, Xing Liang, Jun-Feng Peng, Ming-Hui Zheng, Ju-Dong Li, Yong-Bing Cao, Cheng-Hao Shao, Identifying Liver Metastasis-Related Genes Through a Coexpression Network to Construct a 5-Gene Model for Predicting Pancreatic Ductal Adenocarcinoma Patient Prognosis, 2023, 52, 1536-4828, e151, 10.1097/MPA.0000000000002229
    3. Yuting Jiang, Chengyu Liao, Jianlin Lai, Yunyi Peng, Qilin Chen, Xiaoling Zheng, KRT7 promotes pancreatic cancer metastasis by remodeling the extracellular matrix niche through FGF2-fibroblast crosstalk, 2025, 15, 2045-2322, 10.1038/s41598-024-84129-1
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4348) PDF downloads(272) Cited by(3)

Figures and Tables

Figures(5)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog