Comprehensive analysis of metastasis-related genes reveals a gene signature predicting the survival of colon cancer patients

Introduction

Colorectal cancer (CRC) is one of the leading malignant cancers in the world, and colon cancer accounts for a large part of CRC (Siegel et al., 2017; Siegel, Miller & Jemal, 2017). During the last three decades, with the development of new therapies, such as asbevacizumab and anti-epidermal growth factor receptor antibody (Castro et al., 2013; Knijn, Tol & Punt, 2010), survival outcomes for colon cancer patients with localized and regional disease have improved greatly. However, cancer metastasis remains one of the main causes of death in these patients. One studied showed that the prognosis of patients with distant metastasis remained poor compared to the prognosis of patients with early-stage colon cancer (Siegel, Miller & Jemal, 2017). Accumulation of genetic mutations is recognized as one of the important causes for the pathogenesis and progression of colon cancer (Sameer, 2013). Hence, exploring the mechanisms of cancer metastasis and searching suitable predictors are crucial to the diagnosis and treatment of colon cancer.

Previously, the tumor, node, metastasis (TNM) stage and pathological characteristics of colon cancer have been used to predict the prognosis and facilitate treatment for colon cancer patients. However, there are some limitations in using these methods to assess patients (Marzouk & Schofield, 2011). Recently, several novel biomarkers have been tested with the aim of improving the prediction of therapeutic response and prognosis of colon cancer patients (Demirkol et al., 2017; Hu et al., 2014; Xu et al., 2017). These have aided in the diagnosis and treatment of colon cancer, but the results were inconsistent and need to be further studied.

Thus far, metastasis is a major factor for poor prognoses of colon cancer patients, and the liver is the most common organ affected. However, the molecular mechanism underlying distant metastasis remains unclear. Therefore, a comprehensive analysis of the molecular alterations to identify prognostic indicators is key for the management of colon cancer patients with distant metastasis. In this study, we used the colon cancer data from gene expression omnibus (GEO) and The Cancer Genome Atlas (TCGA) and analyzed the data of patients with colon cancer using primary tumor samples and liver metastasis samples to unveil the genes key to the development colon cancer metastasis and the potential prognostic indicators.

Materials and Methods

Results

Overlapped genes of CRC cells and tissues

The DEGs of GSE2509 and GSE40367 were screened on the basis of the selection criteria after preprocessing the raw data. A total of 341 DEGs between colon cancer cells lines SW420 and SW680 cells were identified in GSE2509 dataset, and 7339 DEGs between primary and metastasis tumor specimens were identified in the GSE40367 dataset. By overlapping the DEGs from the two datasets, we obtained 184 overlapped genes that differentially expressed in both CRC cells and tissues. The result is shown in Fig. 1.

GO and KEGG enrichment analysis

The gene functions of the 184 overlapped DEGs were then analyzed by GO and KEGG enrichment analysis. Using the DAVID online tools, we found that the most enriched GO terms of DEGs related to BP were Signal transduction, and the MF was Protein binding, and the CC was Cytoplasm. The KEGG pathway analysis based on the GO results revealed that Thyroid hormone synthesis was the most significant pathway of the overlapped DEGs. The results are shown in Fig. 2.

PPI network and module screening analysis

Using data from the STRING database, a PPI network for the 184 DEGs consisting of 133 nodes and 138 edges was constructed using Cytoscape software. The overall PPI network is shown in Fig. 3A. Following this, the plug-in MCODE was used to detect the modules in the network, and four modules from the PPI network were identified, with the 3-oxoacid CoA-transferase 1 (OXCT1) (sore: 4), actinin alpha 4 (ACTN4) (sore: 3), interleukin 8 (IL-8) (sore: 3), integrin subunit alpha 3 (ITGA3) (sore: 3) as the hub genes of each module. The top three modules are shown in Figs. 3B–3D.

Prognostic value of overlapped DEGs

The prognostic value of the 184 overlapping DEGs was analyzed using the COAD dataset of TCGA by multiple Cox regression analysis after adjusting the data for age, sex, and TNM stage. The results showed that only aldehyde dehydrogenase 2 (ALDH2), neural precursor cell expressed, developmentally down-regulated 9 (NEDD9), filamin A (FLNA), lamin B receptor (LBR), twinfilin actin binding protein 1 (TWF1), and serine and arginine rich splicing factor 1 (SRSF1) were independent genes that associated with the OS of colon cancer patients, with β values of −1.343, −0.051, 0.492, −0.020, −0.181, and −1.938, respectively. We developed a six-gene signature by calculating the risk score of each gene, and divided the patients into a high-risk group and a low-risk group on the basis of the median risk score (Fig. 4A), the survival status, and the genes expression level shown in Figs. 4B–4C. The survival analysis revealed that this six-gene signature within the high-risk group predicted poor OS of colon cancer patients compared to that in the low-risk group, with p < 0.001 (Fig. 5A). Using survival ROC analysis, we found that the risk score of this six-gene signature could moderately predict the 1-, 3-, 5-year OS of colon cancer patients, as the value of the area under ROC curve was 0.686, 634, and 618, respectively (Fig. 5B).

Discussion

Metastasis accounts for 90% of the mortalities of colon cancer patients and thus has become the most lethal characteristic of colon cancer (Li et al., 2017). Colon cancer patients with localized and regional disease have a good prognosis (the 5-year survival rate is up to 91.1%), but patients with distant metastasis have a much worse prognosis (the 5-year survival rate has dropped to 13.3%) (Siegel et al., 2017). Furthermore, the failure of treatment is mostly caused by the metastatic dissemination of primary tumors (Deliu, Georgescu & Bezna, 2014; Stein & Schlag, 2007). At a molecular level, distinct metastasis of colon cancer is molecularly and clinically distinct from the primary site of origin (Zarour et al., 2017). Thus, analyzing the molecular alteration of colon cancer with distinct metastasis is beneficial for the identification of candidate targets of early diagnosis and treatment of advanced stage colon cancer.

To date, some biomarkers have been identified as candidate targets for early diagnosis and treatment of colon cancer and rectal cancer, including genes, miRNA, lncRNA, and the related signatures. Some gene signatures related to the metastasis potential of each tumor have been described, with promising results. Vellinga et al. (2017) designed a lymphangiogenic gene set, applied it to large datasets of CRC, and found that this lymphangiogenic gene set was related to a worse prognosis. Rokavec et al. (2017) reported a single gene, RNA binding motif protein 47, down-regulated during CRC progression may promote epithelial–mesenchymal transition and metastasis. Furthermore, proteomic studies used exosomes that from cancer cell lines and identified four candidate genes MET proto-oncogene, receptor tyrosine kinase, S100A8, S100A9, tenascin C associated with CRC metastasis (Ji et al., 2013). Other biomarkers, such as a four-miRNA signature (let-7i, miR-10b, miR-221, and miR-320a) (Hur et al., 2015), and a six-lncRNA signature were reported to be promising biomarkers for the metastasis and prognosis of CRC (Hu et al., 2014).

In this study, we identified four hub genes from the subnetwork of the PPI network: OXCT1, ACTN4, IL-8, and ITGA3. OXCT1 is a key enzyme involved in the process of ketone body metabolism, and catalyzes the first and rate-determining step of ketolysis. OXCT1 can be converted into acetyl-CoA during the process of metabolism and takes part in the tricarboxylic acid cycle for the oxidation and the production of adenosine triphosphate (ATP) (Zhang & Xie, 2017). The role of OXCT1 has been implicated in several cancers, included CRC, and OXCT1 was overexpressed in the metastatic CRC cell line CC-M3 (Lee et al., 2016). ACTN4 is a non-muscle-type alpha-actinin, which plays an important role in regulating cytoskeleton organization and is involved in transcriptional regulation of gene expression. ACTN4 encodes a nonmuscle, alpha-actinin isoform that is concentrated in the cytoplasm; participates in metastatic processes; and facilitates the motility, invasion, and metastasis of cancer cells. In CRC, ACTN4 was reported to promote CRC cell line invasion by suppressing focal adhesion maturation (Fukumoto et al., 2015). ITGA3 belongs to the integrin family, joins a β₁ subunit to form an intact integrin and interacts with several extracellular matrix proteins (Nagata et al., 2013). ITGA3 was found to be associated with lymphatic dissemination and local invasiveness in cancers. One study showed that ITGA3 was overexpressed in stages III vs I of CRC patients and is related to OS and disease-free survival (Linhares et al., 2015). IL-8 is an important pro-inflammatory chemokine and plays a role in the recruitment of leukocytes to the sites of infection or tissue injury. Growing evidence suggests that paracrine signaling by tumor-derived IL-8 promotes the trafficking of neutrophils and myeloid-derived suppressor cells into the tumor microenvironment, which is associated with dampened anti-tumor immune responses (David et al., 2016). Lambrechts et al. (2015) observed that IL-8 plasma levels at baseline and subsequent increases in IL-8 were associated with worse progression-free survival in metastatic CRC patients. These studies confirmed our results that the hub genes of the PPI network were crucial to the metastasis of colon cancer, but further studies are warranted to determine the underlying mechanism.

Similar to the previous studies, in this study, a signature constructed by six genes was shown to be a good predictor of the OS of patients with colon cancer. Among these six genes, the ALDH2, NEDD9, SRSF1, and FLNA were reported to be associated with the CRC in several studies. ALDH2 is essential for the metabolism and detoxification of a wide range of endogenous and exogenous aldehyde substrates. ALDH2 is the rate-limiting enzyme in the ethanol metabolism, oxidizing acetaldehyde to acetic acid both in the liver and other tissues (Chen, Joshi & Mochly-Rosen, 2016). As a novel biological marker, ALDH2 is an attractive prospect in the screening, diagnosis, and evaluation of the prognosis of many diseases, and the genetic polymorphism of ALDH2 significantly correlated with the susceptibility to CRC (Li et al., 2016). NEDD9 is a non-catalytic scaffolding protein that assembles complexes involving oncogenic kinases and regulates the magnitude and duration of cell signaling cascades that control multiple processes important to the development and progression of tumors (Shagisultanova et al., 2015). Studies have shown that downregulation of NEDD9 by apigenin can suppresses migration, invasion, and metastasis of CRC cells (Dai et al., 2016). With regard to SRSF1, one study reported that phosphorylation of SRSF1 regulated alternative splicing of tumor-related Rac1b in CRC cells (Goncalves et al., 2014). SRSF1 is classified as an exonic splicing enhancer and recognizes degenerate purine-rich sequence motifs. In addition, it can promote the recognition of both constitutive and alternative exons during the process of spliceosomal assembly (Sanford et al., 2009). FLNA is an actin-binding protein expressed ubiquitously in the body, involved in many cell-signaling pathways, and it is important in the maintenance of cell shape and motility. The mutation of FLNA has been shown to be associated with the neuronal migration, vascular function, connective tissue integrity, and skeletal development (Shelmerdine et al., 2017); FLNA showed low expression in CRC patients and was significantly correlated with the incidence and development of CRC (Tian et al., 2015). LBR, a transmembrane protein of the inner nuclear membrane, interacts with heterochromatin and B-type lamins through its nucleoplasmic amino-terminal domain and is phosphorylated throughout the cell cycle (Duband-Goulet, Courvalin & Buendia, 1998). TWF1 belongs to the ADF-H family and is a conserved actin-binding protein. TWF1 regulates diverse morphological processes through sequestering ADP-actin monomers or capping filament barbed ends (Paavilainen et al., 2007). However, no study has yet reported the role of LBR and TWF1 in CRC and therefore their role needs be investigated further in colon cancers.

Compared with previous studies, this study screened the metastasis-related genes by overlapping the DEGs from cancer tissues and cell lines. These tissues included primary colon cancer tissues and liver metastasis tissues, and the cell lines included primary (SW480) and metastatic (SW620) human isogenic CRC cell lines. Thus, the overlapped DEGs could more reliably reflect the gene alteration in metastatic colon cancer. Some limitations of this study need to be noted. First, although metastasis-related genes of colon cancer were identified and their prognostic value was validated in our study, the results were calculated from microarray or RNA-sequencing technique datasets and thus lack of functional validation of the target genes is one of the major limitations of this study. Therefore, thorough functional experiments for these genes and corresponding downstream events to reveal novel diagnostic and therapeutic targets for colon cancer is necessary. Second, the development of colon cancer metastasis can be caused by many factors, such as KRAS, BRAF mutation, and microsatellite instability, which have been proven to be closely related to colon cancer, but because of the limitations of the datasets, we did not perform a stratified analysis on the basis of these factors. Future studies should be conducted to explore the difference under different conditions. Third, the mean time of follow-up in the validation cohort was 65.3 months; thus, a study including a longer follow-up time is warranted to validate our results in the future. Fourth, this study only used a GSE dataset to perform bioinformatics to search for biomarkers. Future studies should implement Partek Genomics software and Illumina Correlation Engine for bioset analyses.

Conclusions

In this study, we screened gene profiles involved in the metastasis of colon cancer and identified four hub genes from the profiles. We also identified and validated a six-gene signature that can serve as a prognostic indicator of colon cancer. Some genes that have not yet been proven to be associated with colon cancer metastasis may represent new therapeutic targets.