Identification of liver metastasis-associated genes in human colon carcinoma by mRNA profiling

2019-01-18JianlingLiuDanWangChaoqiZhangZhenZhangXinfengChenJingyaoLianJinboLiuGuixianWangWeitangYuanZhenqiangSunWeijiaWangMengjiaSongYapingWangQianWuLingCaoDongWangYiZhang

Chinese Journal of Cancer Research 2018年6期

Jianling Liu， Dan Wang， Chaoqi Zhang， Zhen Zhang， Xinfeng Chen， Jingyao Lian，Jinbo Liu， Guixian Wang， Weitang Yuan， Zhenqiang Sun， Weijia Wang， Mengjia Song，Yaping Wang， Qian Wu， Ling Cao， Dong Wang， Yi Zhang，3

1Biotherapy Center, the First Affiliated Hospital of Zhengzhou University, Zhengzhou 450052, China; 2Department of Oncology, the First Affiliated Hospital of Zhengzhou University, Zhengzhou 450052, China; 3School of Life Sciences, Zhengzhou University, Zhengzhou 450001,China; 4Department of Anorectal Surgery, the First Affiliated Hospital of Zhengzhou University, Zhengzhou 450052, China; 5Henan Key Laboratory for Tumor Immunology and Biotherapy, Zhengzhou 450052, China

Abstract Objective: Liver metastasis, which contributes substantially to high mortality, is the most common recurrent mode of colon carcinoma. Thus, it is necessary to identify genes implicated in metastatic colonization of the liver in colon carcinoma.Methods: We compared mRNA profiling in 18 normal colon mucosa (N), 20 primary tumors (T) and 19 liver metastases (M) samples from the dataset GSE49355 and GSE62321 of Gene Expression Omnibus (GEO) database.Gene ontology (GO) and pathways of the identified genes were analyzed. Co-expression network and proteinprotein interaction (PPI) network were employed to identify the interaction relationship. Survival analyses based on The Cancer Genome Atlas (TCGA) database were used to further screening. Then, the candidate genes were validated by our data.Results: We identified 22 specific genes related to liver metastasis and they were strongly associated with cell migration, adhesion, proliferation and immune response. Simultaneously, the results showed that C-X-C motif chemokine ligand 14 (CXCL14) might be a favorable prediction factor for survival of patients with colon carcinoma. Importantly, our validated data further suggested that lower CXCL14 represented poorer outcome and contributed to metastasis. Gene set enrichment analysis (GSEA) showed that CXCL14 was negatively related to the regulation of stem cell proliferation and epithelial to mesenchymal transition (EMT).Conclusions: CXCL14 was identified as a crucial anti-metastasis regulator of colon carcinoma for the first time,and might provide novel therapeutic strategies for colon carcinoma patients to improve prognosis and prevent metastasis.

Keywords: Colon carcinoma; liver metastasis; mRNA profiling; functions annotation

Introduction

Colon carcinoma is one of the most common malignant diseases with 945,000 new cases every year and is the fourth cause of cancer-related deaths worldwide (1).Unfortunately, about 70% of colon carcinoma patients develop liver metastases. Curative-intent resections can be performed in only 10%-15% of liver metastases (2). In the majority of metastatic patients, the standard treatment remains palliative chemotherapy. However, most colon cancer patients with active metastasis appear to be resistant,or even non-responsive, to current treatments. A major clinical challenge is to explore possible therapeutic targets that are specifically expressed in liver metastatic settings.

There have been many attempts to determine predictive factors or explain the underlying mechanisms for distant metastasis. MicroRNA 34a, microRNA-34a-5p,microRNA-340 are associated with colon carcinoma cell proliferation and metastasis (3,4). In addition, the CpG island methylator phenotype (CIMP) is concordant between primary colon carcinoma and distant metastases(5). Mitogen-activated protein kinase (MAPK) and phosphatidylinositol 3-kinase (PI3K) signaling pathways inhibit metastasis to the liver (6). Alterations in gene expression, protein expression, posttranslational modification, microRNA and linc-RNA have been reported to act a part of role in tumor progression. However, these have not revealed effective predicted factor which is specific to liver metastasis. Transcriptomic changes inherit from genomic information and take place before protein level.Therefore, we attempt to investigate the malignant features of hepatic metastasis microenvironment by RNAsequencing.

Gene expression profiling has become a strategy to identify genes involved in the progression and the prognosis of different cancers. Few attentions were focused on the gene signatures associated with metastatic disease(7). Two studies presented gene signatures associated with metastatic disease containing more than 400 genes. Such long lists of genes are difficult to be used for the development of new therapies (8,9). Pairs of primary and metastatic tumors were analyzed and the samples clustered by patients but not the tissue origin (10,11). The identified genes are specific to colon carcinoma and hepatic metastases, but the precise target is still unknown (12).Comparative profiling of primary colon carcinomas and liver metastases identifies lymphoid enhancer factor-1(LEF1) as a prognostic biomarker (13). However, this research only focused on the development of diagnostic and prognostic markers without trying to identify gene signatures able to distinguish metastatic from primary cancer tissues (13). Therefore, it is most important for us to investigate effective targets for the treatment of liver metastasis.

To identify genes implicated in metastatic colonization of the liver in colon carcinoma, we compared mRNA expression between groups of normal colon mucosa (N),primary tumors (T) and liver metastases (M) samples which from Gene Expression Omnibus (GEO) database. The expression of the differential genes was processed by gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology and Signal network, which are all effective bioinformatics analytical methods. We then verified the clinical significance of identified genes using clinical samples. Our data provide novel information and help further understanding of the liver metastasis cascade of colon carcinoma.

Materials and methods

Microarray data

The transcriptional expression data (GSE49355 and GSE62321) of human colon tumor were downloaded from the GEO database. They were from the same set of patients. It contained 18 normal colon mucosa (N), 20 primary tumors (T) and 19 liver metastases (M) samples.Platforms information were GPL96 [HG-U133A] and GPL97 [HG-U133B] Affymetrix Human Genome U133A/B Array and the datasets were already normalized.

Investigating of differential expression genes (DEGs)

Genes were standardized and interpreted functionally before comparison. Using random variance model (RVM)t-test (14) and the normal colon mucosa group as the control group, the P value and the false discovery rate(FDR) were calculated for each DEG. FDR was calculated to correct the P-value, which controls type I errors. With a threshold of P＜0.05, FDR＜0.05 and fold change (FC) ＞2,DEGs were picked out.

Hierarchical cluster analysis

Hierarchical cluster analysis was performed to ensure good characterizations of screened DEGs between different groups (15). In hierarchical cluster analysis, Pearson correlation was used to calculate the correlation between the genes and samples.

Venn analysis

To identify specific genes of liver metastasis, genes expression in each tissue were input to the web tool Venn Diagrams (http://bioinformatics.psb.ugent.be/webtools/Venn).

GO annotation analysis

Functional analysis of differentially expressed genes was carried out by the GO project (http://www.geneontology.org) on the basis of biological process (16).

Pathway annotation analysis

Pathway analysis was used to identify significant pathways involving DEGs, according to KEGG, BioCarta, and Reactome.

Co-expression network analysis

For each pair of genes, the Pearson correlation coefficient was calculated, and 0.8 was defined as the threshold to construct the network. Within the network analysis, degree of the association is an important factor to determine the relative importance of a gene. We have employed different colors and sizes of node to discriminate the degree of the associations for one gene with the surrounding nodes.The co-expression networks were constructed by Cytoscape (17).

PPI network construction

In order to reveal functional associations between proteins in a genome-wide scale, STRING online tool (18,19) was used to construct a PPI network. In the PPI network, each node represents a protein, and each edge represents an interaction of pairwise proteins. The nodes with a relatively large number of edges were defined as hub proteins.

Gene set enrichment analysis (GSEA)

GSEA was performed by the GSEA software and gene sets used in this work were downloaded from the Molecular Signatures Database. The MSigDB collects various types of gene set and the online pathway database included 1,320 Canonical pathways derived from the pathway databases of BioCarta, KEGG, PID, Reactome and others databases.The data for GSEA analysis is from The Cancer Genome Atlas (TCGA).

TCGA database analysis

TGCA database was derived from UCSC Cancer Browser(https://genome-cancer.ucsc.edu). Overall survival (OS)analysis of colon cancer patients with high and low levels of different genes was shown by using a Kaplan-Meier survival plot. The cut-off values for the genes were the median respectively. We used Kaplan-Meier curves to present the prognosis of the high and low groups. The Wilcoxon logrank test was then conducted on the Kaplan-Meier curves to detect the survival difference between these two groups.All survival analysis was conducted using the R software.

Clinical specimens

Specimens were from colon carcinoma patients who were diagnosed and received operation in the Department of Anus and Intestine Surgery of the First Affiliated Hospital of Zhengzhou University (Zhengzhou, China) from 2011 to 2013. Pre- and post-operative clinical data and other survival-related data were perfected by reviewing the medical records and following-up the patients by telephone. All postoperative specimens were examined by one pathologist and reviewed by another pathologist. Of them, all patients were used as the basis of the present study. The clinical data of the patients are shown in Table 1. Collection of samples in this study was approved by Institutional Ethics Committee of the First Affiliated Hospital of Zhengzhou University (Ethics approval number: Science-2010-LW-1213), and informed consent was obtained from each patient with available follow-up information.

Quantitative real-time polymerase chain reaction (qRTPCR)

Tumor or marginal tissues were cut into 20 mm of pieces and mechanically grinded. Then, total RNA was extracted using Trizol solution (Invitrogen, Waltham, MA, USA).qRT-PCR was performed using specific primers and SYBR Green qPCR Master Mix (Takara, Japan). Listed primers were used: 5’-GGAGCCAAAAGGGTCATCATCTC-3’sense primer and 5’-GAGGGGCCATCCACAGTCTTC T-3’ antisense primer for GAPDH, 5’- CGCTACAGCG ACGTGAAGAA-3’ sense primer and 5’-GTTC CAGGCGTTGTACCAC-3’ antisense primer for CXC chemokine ligand 14 (CXCL14). GAPDH was used as an internal control. With the 2-ΔΔCtmethod, we compared the expression level of clinical samples (20). For each sample,the expression of CXCL14 as well as GAPDH was examined, the relative expression of CXCL14 was calculated by using the 2-ΔCtvalue of CXCL14 dividing the 2-ΔCtvalue of GAPDH (20).

Immunohistochemistry

Paraffin-embedded tissues of 45 colon cancer samples wereexamined for the expression of CXCl14 protein (Abcam,Cambridge, UK; 1:200). Sections were treated with 3%H2O2and 5% bull serum albumin (BSA) and incubated with primary antibodies overnight at 4 °C. After incubation with horseradish peroxidase (HRP)-conjugated secondary antibody for 1 h at 37 °C, sections were washed and counterstained with hematoxylin, and visualized under a microscope (Olympus, Shinjuku, Japan) (21).

Table 1 Characteristics of patients with colon carcinoma

Statistical analysis

Clinicopathologic factors were compared by using the χ2test and continuous variables were compared by using the Student t test or one-way analysis of variance (ANOVA)analysis. Kaplan-Meier analysis and the log-rank test were used for survival analysis. Univariate and multivariate logistic regression models identified the association between CXCL14 expression and clinical characteristics.P＜0.05 was considered statistically difference. All statistics associated with clinical samples were performed using Prism 7 (GraphPad Software Inc., La Jolla, USA).Statistical analysis of significance was calculated by ANOVA followed by Tukey’s post hoc test with SPSS 16.0 for Windows (SPSS Inc., Chicago, IL, USA). The bioinformatics analysis was used by using R software(Version 3.4; R Foundation for Statistical Computing,Vienna, Austria).

Results

Gene expression analysis

We used the public transcriptome sequencing dataset(GSE49355 and GSE62321) from GEO database,including 18 normal colon mucosa (N), 20 primary tumors(T) and 19 liver metastases (M) samples. Detailed sample information could be found in Supplementary Table S1.Expression profiling of the 57 samples was conducted on Affymetrix human U133A/B chips. Expression profiling of the 57 samples was conducted on Affymetrix human U133A chips containing 22,200 probes corresponding to about 12,700 genes. These gene expression data have been performed normalization and log2 transformation.Hierarchical cluster analysis showed that normal samples clustered together and were relatively well separated from T and M samples in GSE49355 and GSE62321 (Figures 1A-D).

Identification of specific gene signatures

To identify molecular signatures that regulate distant metastasis in colon carcinoma, we compared mRNA expression levels in T vs. N and M vs. N. After analyzing the transcriptomic changes of T vs. N, a total of 1,646 DEGs including 861 up-regulated and 785 down-regulated transcription factors were screened out from GSE49355,and a total of 868 DEGs including 477 up-regulated and 391 down-regulated transcription factors were also identified in GSE62321. Of 1,809 DEGs, 869 were downregulated and 940 overexpressed in M vs. N in GSE49355.The volcano plot of DEGs distribution was also presented 934 DEGs including 468 up-regulated and 466 downregulated when comparing M with N in GSE62321 (Figure 2A, B) (P＜0.05, FDR＜0.05, FC＞2, respectively). Based on the fact that the GSE49355 and GSE62321 were from the same panel of patients but different platform, union analysis was first performed and 719 specific genes related to liver metastasis of colon carcinoma were identified(Supplementary Table S2). However, taken into account that some of the 719 genes might be due to a single platform error, we took the intersection analysis here for obtaining higher accurate genes. The results showed that 179 genes might play an important role in the metastasis of cancer and were altered in M vs. N. Excluding 157 genes associated with tumor development, 22 genes were specific for liver metastasis (Figure 2C). Subsequently, unsupervised hierarchical cluster analysis was performed on selected 22 genes expression data using Pearson correlation-based distance and average clustering. Considerable patients’non-pairing of N and M samples was observed in the dendrogram. Most of the specific genes showed a significantly differential expression between N and M samples (Figure 2D). Details were shown in Supplementary Table S3.

Significant GOs and pathways

Figure 2 Identification of specific genes associated with liver metastasis in colon carcinoma. (A, B) With a threshold of P＜0.05, false discovery rate (FDR)＜0.05 and fold change ＞2, differential expression genes (DEGs) were picked out by volcano plot when comparing 20 primary tumors (T) samples with 18 normal colon mucosa (N) samples and 19 liver metastases (M) samples with 18 N samples from GSE49355 and GSE62321; (C) Venn diagram of commonly DEGs in comparison groups; (D) Hierarchical clustering analysis of specific genes associated with liver metastasis in the two datasets.

All of the identified genes were used to predict the functional categories with GO annotation. They were involved in different biological processes, molecular functions and cellular components. It was found that the differential expression of the 22 genes mainly participated in 153 significant GOs (Supplementary Table S4). It was concluded that the specific genes were mainly involved in immune response, metabolic process and cell adhesion.Blood coagulation, platelet activation and degranulation,acute-phase responses, negative regulation of endopeptidase activity, and complement activation may take part in liver metastasis of colon carcinoma.

In order to identify the key pathways the specific genes were involved in, we performed pathway analysis. Fifty-six KEGG biological pathways were annotated (Supplementary Table S5). The major regulated biological pathways include complement and coagulation cascades, metabolic pathways,PI3K-protein kinase B (AKT) signaling pathway, pathways in cancer, focal adhesion, Staphylococcus aureus infection,carbon metabolism, chemokine signaling pathway, and biosynthesis of amino acids. The results revealed the genes play an important role in pathways related to cancer cell migration, such as PI3K-AKT signaling pathway, focal adhesion and chemokine signaling pathway.

Dynamic gene network analysis

All the screened 22 DEGs were then subjected to a gene co-expression analysis network with k-core algorithm to determine which genes may play a potential role in the colon carcinoma metastasis. The gene-gene interaction network was constructed as shown in Figure 3A. The degree of a node describes the number of links of one gene with others, which had shown in the gene network. The node with larger diameter in the network means more important values. Importantly, six genes (SERPINA1, UPB1,FTCD, F5, EMCN, GIMAP6) belonged to the most significant genes, which involved in acute-phase response,metabolic process, angiogenesis, endothelial cell migration and proliferation, cell adhesion. The SERPINA1, UPB1,FTCD and F5 genes were up-regulated, but EMCN and GIMAP6 genes were down-regulated (Figure 3A).Furthermore, it was obvious that CXCL14, which was associated with cell migration and immune response, was also down-regulated.

Figure 3 Functional associations between screened genes. (A) From the total differential genes, 22 specific genes about liver metastasis of colon carcinoma were constructed a gene co-expression network with k-core algorithm. Red cycle nodes represent up-regulated genes, blue cycle nodes represent down-regulated genes; (B) Protein-protein interaction (PPI) network of screened genes. Each node represents one gene; edges indicate the interaction relationship.

A PPI network of genes

To further define the interaction between the screened 22 DEGs, we used STRING database to construct the PPI network. The PPI network consisted of 6 nodes interacting by 29 edges, the remaining 16 DEGs failed to form the PPI pairs. It was concluded that FTCD, APOB, APOH, PLG, F5,SERPINA1 were closely linked (Figure 3B).

Prognostic values of highlighted DEGs

To evaluate the prognostic values of the 22 DEGs, we further investigated the associations of the DEGs with OS of patients by Kaplan-Meier and log-rank analysis. Because neither ACSM2A nor FTCD’s positive expression rate, the percentage of sample numbers with gene expression accounting for all sample numbers, was less than 50% in TCGA database, survival analysis was used to estimate the prognosis value of the other 20 genes. We found that patients with lower CXCL14, SERPINA1 expression demonstrated poorer survival than patients with higher expression (P=0.0388; P=0.0109; Figure 4). However, it was contradictory that SERPINA1 expression up-regulated in liver metastasis tissues indicated benefit prognosis. We therefore further researched the gene CXCL14 which were specifically involved in anti-liver metastasis process of colon carcinoma and predicted beneficial prognosis.

GSEA analysis of CXCL14

Based on above results, we have found that CXCL14 play a key role in liver metastasis of colon carcinoma. Then, it was quite necessary to predict biological functions of this gene.Analysis of GSEA, a powerful tool to infer the biological function, was performed. The results showed that genes associated with cell aging, negative regulation of stem cell proliferation and epithelial to mesenchymal transition(EMT), which were closely related to cancer metastasis(22-24) were significantly enriched in CXCL14-high samples of colon carcinoma (Figure 5). These observations suggested that CXCL14 may be a predicted indicator of patients with colon carcinoma liver metastasis.

Validation of CXCL14 expression and its clinical relevance with clinical samples

To further demonstrate the clinical significance of CXCL14 expression in patients with colon carcinoma, the association between CXCL14 expression and various clinicopathological variables was investigated by real-time quantitative PCR in 103 colon carcinoma patients. The clinicopathological data of the patients are detailed in Table 1. CXCL14 expression showed a high level in colon carcinoma patients with early stage, non-liver metastasis,middle histological differentiation (Figure 6A, B, C). At the protein level, the results also showed that CXCL14 expression was lower in patients with liver metastasis(Figure 6D, E). Then 103 colon carcinoma samples were stratified into “high” and “low” according to the median 0.045 127 of CXCL14 level. We found that low expression of CXCL14 was strongly correlated with advanced liver metastasis (P=0.01), overall stage (P=0.0001), abnormal CA72-4 value (P=0.0001), tumor size (P=0.001) and site of lesion (P=0.006) (Table 2).

Figure 4 Association of expression of C-X-C motif chemokine ligand 14 (CXCL14) and SERPINA1 with overall survival (OS) of 250 patients from The Cancer Genome Atlas (TCGA) data. Kaplan-Meier survival analysis of OS based on expression status provided associations of differential expression genes (DEGs) with OS of 250 patients from TCGA data. Cut-off values for genes were the median respectively. (A) CXCL14 [hazard rate (HR)=1.551; P=0.0388]; (B) SERPINA1 (HR=1.703; P=0.0109). With x-axis from left to right, the expression of CXCL14 was from high to low.

Figure 5 Gene set enrichment analysis (GSEA) analysis of C-X-C motif chemokine ligand 14 (CXCL14). GSEA showed that CXCL14 was associated with (A) Cell aging; (B) Stem cell proliferation; and (C) Epithelial to mesenchymal transition (EMT).

To examine the potential of CXCL14 to predict liver metastasis, logistic regression analysis was used. Univariate analyses revealed that low CXCL14 level [odds ratio(OR)=2.13; P=0.04], high CA 72-4 level (OR=6.9; P=0.01), and advanced overall stage (OR=6.0; P=0.02)were associated with liver metastasis. In multivariate analyses, CXCL14 (OR=1.24; P=0.03) and CA 72-4 levels(OR=2.35; P=0.04) were independent predictor of liver metastasis (Table 3).

Figure 6 Association of C-X-C motif chemokine ligand 14 (CXCL14) expression with clinical characteristics and overall survival (OS) of patients with colon carcinoma. The mRNA expression level of CXCL14 in different groups of (A) Colon carcinoma patients with liver metastasis (Yes) and without liver metastasis (No); (B) TNM stage; and (C) Histological differentiation. L, low differentiation; M, moderate differentiation (n=103). (D) Kaplan-Meier curves show the association between mRNA expression level of CXCL14 and OS (n=103); (E)Immunohistochemical staining results of tumor tissue in colon cancer with liver metastasis, and without liver metastasis (×200); (F) Immune responsive score (IRS) of CXCL14 in colon cancer with liver metastasis (Yes) and without liver metastasis (No) (n=45); (G) Kaplan-Meier curves show the association between expression of CXCL14 and OS according to the immunohistochemical results (n=45)(＊, P＜0.05).

Table 2 Association between CXCL14 expression and clinicopathological features of patients with colon carcinoma (N=99)

Moreover, Kaplan-Meier survival analysis revealed that higher CXCL14 expression was significantly associated with better survival in patients with colon carcinoma(Figure 6D, G). Simultaneously, as shown in Table 4,univariate analysis revealed that low CXCL14 level [hazard ratio (HR)=0.348; P=0.001] and liver metastasis(HR=2.742; P=0.037) were significantly associated with poor prognosis. Importantly, a multivariate Cox’s regression analysis revealed that CXCL14 level (HR=0.388;P=0.0001) and liver metastasis (HR=1.174; P=0.045) were independent prognostic factors for the OS of patients with colon carcinoma (Table 4). Collectively, these results suggest that CXCL14 expression status plays an important role in predicting prognosis and liver metastasis in patients with colon carcinoma.

Table 3 Logistic regression model analysis of liver metastasis predictors in patients with colon carcinoma

Discussion

In relation to disease relapse, liver metastasis is the most major recurrent mode of colon carcinoma. When patients were firstly diagnosed, some of them were found to have distant metastasis, which might result in unfavorable prognosis. Thus, it is critical to identify an effective indicator that predicts the liver metastasis of colon carcinoma to provide new methods for therapy. It is notable that RNA-sequencing data and microarray-based expression profiling data provide a more comprehensive and accurate understanding of carcinogenesis and cancer progression at the molecular level. In this study, mRNA profiling by microarray from GEO was used to identify a number of novel genes related to colon carcinoma liver metastasis.

To identify new predictors regulating liver metastasis in colon carcinoma, we compared mRNA expression levels in M vs. N and T vs. N. Based on the fact that the GSE49355 and GSE62321 were from the same panel of patients but different platform, we took the intersection analysis ratherthan the union analysis for obtaining more accurate genes.The results showed that 22 genes were specifically related to liver metastasis. To further demonstrate their function and signaling pathway, we performed annotation analysis and verified that the genes were strongly associated with: 1)cell migration, adhesion, proliferation (cell adhesion/focal adhesion/ chemokine signaling pathway/PI3K-AKT signaling pathway/APOH/F5/CXCL14); and 2) immune response (innate immune response/complement activation/acute-phase response/SERPINA1/CXCL14).

Table 4 Cox’s proportional hazard model analysis of prognostic factors in patients with colon carcinoma

It is well known that metastasis is closely related to colon cancer patients’ survival, and almost 80% of metastases occurred in liver. Therefore, we analyzed the prognosis value of the screened 22 specific liver metastasis genes though TCGA database. Because neither ACSM2A nor FTCD’s positive expression rate was less than 50%,survival analysis was used to focus on the other 20 specific genes. The results suggest that CXCL14 and SERPINA1 may be favorable prediction factors for colon carcinoma patients’ survival. However, our present study found that SERPINA1 expressed higher level in metastatic liver tissues when comparing with normal tissues. Recent studies have been reported that SERPINA1, a protease inhibitor that can act on a variety of targets such as serine proteases,has been proposed as a poor prognosis biomarker for various diseases, including papillary thyroid carcinoma (25),lung cancer (26) and breast carcinoma (27). As for our inconsistent results of SERPINA1, we decided to focus on the CXCL14 in proceeding research.

CXCL14, is an orphan member of the CXC chemokine subfamily. CXCL14 mRNA and protein are ubiquitously expressed in normal tissues, but are absent in tumor cell lines and in primary tumors (28,29). CXCL14 level in colon carcinoma tissues with lymphoid metastasis was significantly lower than that in tumor tissues without lymphoid metastasis (30). However, the effect of CXCL14 on colon carcinoma liver metastasis remains unclear. In our study, we provided evidences to show that the expression of CXCL14 was down-regulated in M vs. N and it was closely correlated with a beneficial survival outcome. Combined with the similar results that CXCL14 mediated suppression of tumor metastasis in lung cancer and Ewing sarcoma (31),we could conclude that it may play an important role in regulating colon carcinoma metastasis. On the contrary,Liu and colleagues previously described that CXCL14 induced metastasis (32). CXCL14-positive cancer associated fibroblast involved in ovarian cancer metastatic progression. This inconsistency may be caused by the intrinsic characteristic differences of different subtypes of human cancer.

Here, we firstly reported the role of CXCL14 in colon carcinoma liver metastasis. However, its underlying mechanism remains to be elucidated. Combined with the evidence that: 1) absence of CXCL14 expression in many malignant tissues is in agreement with the deficiency of effective antitumor immune responses in cancer patients.CXCL14 may act as chemo-attractant for monocytes,dendritic cells (DC) and (natural killer) NK cells; 2)CXCL14 may also influ¬ence the proliferation, invasion and migration of tumor cells via auto/paracrine pathways;3) CXCL14 may suppress tumor vasculature by inhibiting the chemotaxis of vascular smooth muscle cells and the formation of microvas¬cular systems (33,34), and therefore suppresses the metabolism and growth of a tumor. Thus, it might indicate that CXCL14 plays an important role in regulating liver metastasis of colon carcinoma through suppression of cancer cells and migration of leukocytes.

Our data suggest that modulating CXCL14 expression could exert tumor suppression effects on colon carcinoma.CXCL14 expression is suppressed by epidermal growth factor (EGF) and can be restored by treatment with an EGF receptor (EGFR) tyrosine kinase inhibitor in head and neck squamous cell carcinoma (HNSCC) cells (35).Reversing the promoter hypermethylation of CXCL14 could be a feasible approach to restore anti-tumor immune responses to treat oral cancers (36). In summary, upregulating the CXCL14 level may be a valuable adjuvant treatment to improve the outcomes of patients.

Conclusions

Taken together, our data firstly indicated that the CXCL14 expression level was down-regulated in metastatic liver tissues compared to non-tumor tissues. The absence of CXCL14 contributed to the cancer metastasis that then causes poor outcomes of patients. This is the first report that CXCL14 exerts an anti-metastasis effect on colon carcinoma via the screening of bioinformatics and the further validation of clinical samples. The expression level of CXCL14 may be a valuable adjuvant parameter to predict the liver metastasis and prognosis of patients with colon carcinoma and provides a potential future therapeutic strategy.