Background
Pancreatic ductal adenocarcinoma (PDAC) is an often fatal malignancy with an extremely low survival rate. Liver metastasis, which causes high mortality, is the most common recurring metastasis for PDAC. However, the mechanisms underlying this liver metastasis and associated candidate biomarkers are unknown.
Methods
We performed mRNA profiling comparisons in 8 primary tumors (T) and 12 liver metastases (M) samples using the Gene Expression Omnibus (GEO) database. After determining differentially expressed genes (DEG), gene ontology (GO), pathway enrichment and protein–protein interaction (PPI) network analyses were performed to determine DEG functions. Then, Cytoscape was used to screen out significant hub genes, after which their clinical relevance was investigated using The Cancer Genome Atlas (TCGA) resources. Furthermore, prognosis-associated gene expression was validated using Oncomine and TCGA database. Lastly, associations between prognosis-associated genes, immune cells and immunological checkpoint genes were evaluated using the Tumor Immune Estimation Resource (TIMER).
Results
In total, 102 genes were related to liver metastasis and predominantly involved in cell migration, motility, and adhesion. Using Cytoscape, this number was narrowed down to 16 hub genes. Elevated mRNA expression levels for two of these genes, SPARC (P = 0.019) and TPM1 (P = 0.037) were significantly correlated with poor disease prognosis. For the remaining 14, expression was not related to overall patient survival. SPARC had higher expression in patients with metastatic PDAC than those with non-metastatic PDAC in TCGA dataset. SPARC and TPM1 levels were also positively correlated with the immune infiltration of specific cell types. Additionally, both genes exhibited strong co-expression associations with immune checkpoint genes.
Conclusions
Combined, we suggest SPARC has high potential as biomarker to predict liver metastasis during PDAC. Additionally, both SPARC and TPM1 appeared to recruit and regulate immune-infiltrating cells during these pathophysiological processes.