Research Paper Volume 9, Issue 3 pp 914—931

Overexpression of Rab1B and MMP9 predicts poor survival and good response to chemotherapy in patients with colorectal cancer

Xian-Zi Yang1,2,3, *,, Shu-Zhong Cui3, *,, Li-Si Zeng3, , Tian-Tian Cheng3, , Xiao-Xing Li1,2, , Jun Chi1,2,4, , Ren Wang1,2, , X.F. Steven Zheng1,2,5, , Hui-Yun Wang1,2, ,

  • 1 State Key Laboratory of Oncology in South China, Sun Yat-Sen University Cancer Center, Guangzhou, Guangdong, 510060, ; China
  • 2 Collaborative Innovation Center for Cancer Medicine, Sun Yat-Sen University Cancer Center, Guangzhou, Guangdong, 510060, ; China
  • 3 Department of Abdominal Surgery, Affiliated Cancer Hospital & Institute of Guangzhou Medical University of Guangzhou Medical University, Guangzhou Medical University, Guangzhou, Guangdong, 510095, ; China
  • 4 Department of Endoscopy and Laser, Sun Yat-Sen University Cancer Center, Guangzhou, Guangdong, 510060, ; China
  • 5 Rutgers Cancer Institute of New Jersey/Department of Pharmacology, Robert Wood Johnson Medical School, Rutgers, The State University of New Jersey, New Brunswick, NJ 08901, USA
* Equal contribution

Received: January 22, 2017       Accepted: March 6, 2017       Published: March 18, 2017      

https://doi.org/10.18632/aging.101200
How to Cite

Abstract

Rab1B has recently been reported to be involved in human cancer, but the role of Rab1B in colorectal cancer (CRC) remains unclear. In this study, we investigated the expression of Rab1B and MMP9 in CRC by qRT-PCR, immunoblot and immunohistochemistry and analyzed the clinical significance. The results show that Rab1B and MMP9 are increased at both mRNA and protein levels in CRC cell lines and tissues, as measured by qRT-PCR and immunoblotting. The high protein expression of Rab1B and MMP9 in 179 CRC tissues is associated with deep tumor invasion, lymph-node metastasis and advanced TNM stage. Survival analysis indicates that patients with overexpression of Rab1B or MMP9 have significantly worse overall survival and progression-free survival, but better response to chemotherapy than those with low expression of proteins, and that Rab1B is an independent prognostic factor for CRC patients. Furthermore, when Rab1B and MMP9 are combined into a new risk model, it has a remarkably better prediction of prognosis than each protein alone. In conclusion, Rab1B and MMP9 are potential prognostic biomarkers and their combination significantly improves predictive power for survival and chemotherapy response in CRC patients.

Introduction

Colorectal cancer (CRC) is the third most common cancer in men worldwide. Although the 5-year survival rate of CRC patients is more than 60%, it is still the second leading cause of cancer-related death in the developed countries [1,2]. The primary reason for the high mortality of CRC is due to its high recurrence and metastasis in approximately half of all the patients, which are extremely difficult to cure [2,3]. The standard therapy for CRC patients with stage III and high-risk stage II is resection plus postoperative adjuvant chemotherapy [4]. However, a subpopulation of the patients does not benefit from adjuvant chemotherapy, while suffering long-term toxicity [5]. At present, the commonly high-risk features, such as poor differentiation, T4 tumors and margin involvement, cannot precisely distinguish patients with high-risk from low-risk for disease recurrence and metastasis [6]. Because tumor behaviors are determined by oncogenic alterations in molecular and cellular processes, such aberrant events can be explored as predictive biomarkers for the progression and metastasis of human cancer. Moreover, a better understanding of the molecular mechanisms underlying the progression and metastatic process of CRC can lead to new cancer drug target [7]. There is a pressing need for identifying key molecular markers for prediction of recurrence, metastasis and chemotherapeutic outcome to improve the survival of patients with CRC.

Rab1A and Rab1B share highly homologous sequence and function in membrane trafficking between endoplasmic reticulum to Golgi apparatus [8]. They have also been reported to be involved in the progression and metastasis of human cancer by regulating different cell signaling pathways [9]. In the previous studies, we demonstrated that Rab1A is overexpressed and associated with tumor prognosis in patients with CRC and hepatocellular carcinoma (HCC), promoting tumor growth and metastasis through activating mTORC1 signaling [10,11]. It was reported that Rab1A enhances migration of breast cancer cells through promotion of ITGB1 recycling to the plasma membrane [12]. RablA overexpression was also reported in human lung cancer, which is correlated with tumor volume and stage, but the tumorigenic function of Rab1A is not dependent on mTOR or MAPK signaling [13]. These results suggest that in different types of tumor, Rab1A promotes cell migration and invasion through distinct mechanisms. Rab1B, the second Rab1 member, was also found to be increased in cervical cancer [14] and HCC [15], but another study showed that Rab1B protein was down-regulated and inhibited tumor proliferation and migration via regulating TGF-β/Smad pathway in triple-negative breast cancer (TNBC) [16], indicating that Rab1B may have different roles in different cancer types. Recently, Zhai et al reported that Rab1B mRNA was increased in a small sample size (23 cases) of CRC [17]. However, the clinical significance of Rab1B expression in CRC patients has not been evaluated.

Matrix metalloproteinases (MMPs) play an important role in degradation of extracellular matrix (ECM) and basement membranes (BM). In CRC, previous studies indicate that the expression of MMP9 is associated with metastasis [18-20] and poor prognosis [21-23]. Interestingly, a recent study showed that Rab1A knockdown decreases MMP9 expression and inhibites MMP9-mediated invasiveness of human oral squamous cell cancer cells, suggesting that Rab1A regulates MMP9-mediated invasiveness [24]. Given the high homology between Rab1A and Rab1B, it would be interesting to determine whether Rab1B protein promotes invasiveness and metastasis by upregulating the expression of MMP9 in CRC.

In this study, we investigated the expression of Rab1B and MMP9 proteins and their relationship in CRC tissues and cell lines. We further analyzed the correlation between Rab1B and MMP9 expressions and clinicopathological parameters as well as prognosis in CRC patients. Finally, we evaluated the predictive value of Rab1B and/or MMP9 protein expressions in CRC patients who underwent chemotherapy treatment.

Results

The expression of Rab1B and MMP9 is up-regulated in CRC cell lines

To determine the expression of Rab1B and MMP9 in CRC, we first performed immunoblot and RT-PCR in 11 CRC cell lines and one human normal colon cell line CCD112CoN. The result showed that the expression levels of Rab1B and MMP9 proteins were significantly increased in 81.8% (9/11) and 63.6% (7/11) of CRC cell lines, respectively, compared with those in normal cell line (Fig. 1A). As expected, relative expressions of Rab1B and MMP9 mRNA in CRC cell lines were also higher than those in CCD112CoN cells (Fig. 1B). Furthermore, we found that there is a positive correlation trend between Rab1B and MMP9 expression in both mRNA and protein (for mRNA: r = 0.473, P=0.121; for protein: r = 0.537, P = 0.072, Fig. 1C and 1D), although the correlations were not statistically significant. These results indicate that both Rab1B and MMP9 are up-regulated in CRC cells and may have a positive correlation.

Rab1B and MMP9 are overexpressed in CRC cell lines. (A) The expression of Rab1B and MMP9 proteins in a panel of CRC cell lines and an immortalized colon cell line is determined by immunoblot. GAPDH is used as a loading control. (B) Relative expression of Rab1B and MMP9 mRNA (normalized to GAPDH) in the same set of cell lines as in (A) is examined by real-time quantitative PCR. (C) The correlation between Rab1B and MMP9 proteins (normalized to GAPDH) in the CRC cell lines was determined by Spearman correlation assay. (D) Spearman correlation analysis is used to analyze the correlation between Rab1B and MMP9 mRNAs in CRC cell lines.

Figure 1. Rab1B and MMP9 are overexpressed in CRC cell lines. (A) The expression of Rab1B and MMP9 proteins in a panel of CRC cell lines and an immortalized colon cell line is determined by immunoblot. GAPDH is used as a loading control. (B) Relative expression of Rab1B and MMP9 mRNA (normalized to GAPDH) in the same set of cell lines as in (A) is examined by real-time quantitative PCR. (C) The correlation between Rab1B and MMP9 proteins (normalized to GAPDH) in the CRC cell lines was determined by Spearman correlation assay. (D) Spearman correlation analysis is used to analyze the correlation between Rab1B and MMP9 mRNAs in CRC cell lines.

The protein expression of Rab1B and MMP9 is increased and has a positive correlation with each other in CRC tissues

To further investigate the expression and clinical significance of Rab1B and MMP9 protein in CRC, we collected 179 pairs of cancer and corresponding adjacent non-tumor tissues from CRC patients. Their demographic and clinicopathological data are shown in Table 1. The protein expression of Rab1B and MMP9 in the paired CRC and non-tumor tissues was examined by immunohistochemistry (IHC). Both Rab1B and MMP9 proteins are mainly distributed in the cytoplasm and membrane of CRC cells (Fig. 2A) with 81% (145/179) and 71.5% (128/179) of CRC samples displaying higher expression of Rab1B and MMP9 proteins than the matched adjacent non-tumorous tissues, respectively (P <0.0001, Fig. 2B and 2C). To verify the IHC results, we performed immunoblot analysis on another 8 paired CRC and non-tumorous tissues. As expected, high expression of Rab1B and MMP9 was found in 75% (6/8) and 100% CRC tissues (Fig. 2D), respectively, which is consistent with the results of IHC. Furthermore, we analyzed the relationship between Rab1B and MMP9 expressions. In the same cohort of CRC tissue samples, 71.1% (64/90) samples with MMP9 overexpression had high expression of Rab1B protein, while 69.7% (62/89) samples with MMP9 down-regulation also had low expression of Rab1B (P < 0.001, Table 2 and Fig. 2E). Pearson’s correlation analysis showed that there is a positive correlation between Rab1B and MMP9 protein expressions in 179 CRC samples (Fig. 2F, r =0.388, P< 0.001), which is similar to that in CRC cell lines. These data show that Rab1B and MMP9 proteins are significantly co-overexpressed in CRC tissues.

Table 1. Clinical characteristics of patients with colorectal cancer

CharacteristicsN (%)
Gender
  Male95 (53.1%)
  Female84 (46.9%)
Age (years)
  <65117 (65.4%)
  ≥ 6562 (34.6%)
Pathological grade
  I-II147 (82.1%)
  III-IV32 (17.9%)
Tumor location
  Colon94 (52.5%)
  Rectum85 (47.5%)
Tumor size
  < 5 cm97 (54.2%)
  ≥ 5 cm80 (44.7%)
Tumor depth
  Shallow21 (11.7%)
  Deep68 (38.0%)
N stage
  N099 (55.3%)
  N1-280 (44.7%)
TNM stage
  I17 (9.5%)
  II81 (45.3%)
  III81 (45.3%)
Intraoperative chemotherapy
  No126 (70.4%)
  Yes53 (29.6%)
Adjuvant chemotherapy
  No109 (60.9%)
  Yes70 (39.1%)
Preoperative CEA (ng/ml)
  0-5127 (70.9%)
  > 552 (29.1%)
Preoperative CA199 (ng/ml)
  0-35154 (86.0%)
  > 3525 (14.0%)
Shallow: the depth of tumor invasion within mucosa and muscularis; Deep: the depth of tumor invasion beyond serosa; TNM, tumor node metastasis; N, lymph node; CEA, Carcinoembryonic antigen; CA199, carbohydrate antigen 199.
Rab1B and MMP9 expressions are significantly increased in colorectal cancer tissues. (A) Shown are representative immunohistochemistry (IHC) staining of Rab1B and MMP9 in CRC and adjacent non-tumor tissues. Scale bars represent 50 μm. (B) Comparison of Rab1B protein expressions between CRC tissues and matched adjacent non-tumorous tissues. (C) Comparison of MMP9 protein expressions between CRC tissues and matched adjacent non-tumorous tissues. (D) The expression levels of Rab1B and MMP9 proteins in eight pairs of CRC tissues (T) and adjacent non-tumorous tissues (N) are analyzed by immunoblot. (E) Concordance of Rab1B and MMP9 expressions in CRC. Consecutive CRC sections were analyzed for the expression of Rab1B and MMP9 by IHC. Shown are two representative cases. Scale bars represent 50 μm. (F) The correlation between Rab1B and MMP9 protein expressions in CRC tissues as evaluated by Spearman correlation analysis.

Figure 2. Rab1B and MMP9 expressions are significantly increased in colorectal cancer tissues. (A) Shown are representative immunohistochemistry (IHC) staining of Rab1B and MMP9 in CRC and adjacent non-tumor tissues. Scale bars represent 50 μm. (B) Comparison of Rab1B protein expressions between CRC tissues and matched adjacent non-tumorous tissues. (C) Comparison of MMP9 protein expressions between CRC tissues and matched adjacent non-tumorous tissues. (D) The expression levels of Rab1B and MMP9 proteins in eight pairs of CRC tissues (T) and adjacent non-tumorous tissues (N) are analyzed by immunoblot. (E) Concordance of Rab1B and MMP9 expressions in CRC. Consecutive CRC sections were analyzed for the expression of Rab1B and MMP9 by IHC. Shown are two representative cases. Scale bars represent 50 μm. (F) The correlation between Rab1B and MMP9 protein expressions in CRC tissues as evaluated by Spearman correlation analysis.

Table 2. The relationships of Rab1B and MMP9 expressions with clinicopathological characteristics in patients with colorectal cancer

Clinicopathological characteristicsRab1B ExpressionP valueMMP9 ExpressionP value
LowHighLowHigh
Gender
  Male45 (47.4%)50 (52.6%)0.50348 (50.5%)47 (49.5%)0.698
  Female44 (52.4%)40 (47.6%)40 (47.6%)44 (52.4%)
Age (years)
  <6555 (47.0%)62 (53.0%)0.31956 (47.9%)61 (52.1%)0.633
  ≥ 6534 (54.8%)28 (45.2%)32 (51.6%)30 (48.4%)
Tumor location
  Colon42 (44.7%)52 (55.3%)0.15645 (47.9%)49 (52.1%)0.717
  Rectum47 (55.3%)38 (44.7%)43 (50.6%)42 (49.4%)
Tumor size
  < 5 cm52 (52.5%)47 (47.5%)0.40453 (53.5%)46 (46.5%)0.193
  ≥ 5 cm37 (46.3%)43 (53.8%)35 (43.8%)45 (56.2%)
Pathological grade
  I-II73 (49.7%)74 (50.3%)0.97274 (50.3%)73 (49.7%)0.499
  III-IV16 (50.0%)16 (50.0%)14 (43.8%)18 (56.2%)
Tumor depth
  Shallow21 (91.3%)2 (8.7%)<0.00118 (78.3%)5 (21.7%)0.003
  Deep68 (43.6%)88 (56.4%)70 (44.9%)86 (55.1%)
N stage
  N059 (59.6%)40 (40.4%)0.00367 (67.7%)32 (32.3%)0.000
  N1-230 (37.5%)50 (62.5%)21 (26.2%)59 (73.8%)
TNM stage
  I15 (88.2%)2 (11.8%)0.00116 (94.1%)1 (5.9%)<0.001
  II43 (53.1%)38 (46.9%)50 (61.7%)31 (38.3%)
  III31 (38.3%)50 (61.7%)22 (27.2%)59 (72.8%)
Intraoperative chemotherapy
  No67 (53.2%)59 (46.8%)0.15465 (51.6%)61 (48.4%)0.317
  Yes22 (41.5%)31 (58.5%)23 (43.4%)30 (56.6%)
Adjuvant chemotherapy
  No57 (52.3%)52 (47.7%)0.39056 (51.4%)53 (48.6%)0.460
  Yes32 (45.7%)38 (54.3%)32 (45.7%)38 (54.3%)
Preoperative CEA (ng/ml)
  0-562 (48.8%)65 (51.2%)0.70662 (48.8%)65 (51.2%)0.886
  > 527 (51.9%)25 (48.1%)26 (50.0%)26 (50.0%)
Preoperative CA199 (ng/ml)
  0-3579 (51.3%)75 (48.7%)0.29578 (50.6%)76 (49.4%)0.323
  > 3510 (40.0%)15 (60.0%)10 (40.0%)15 (60.0%)
Rab1B expression
  Low--62 (69.7%)27 (30.3%)<0.001
  High--26 (28.9%)64 (71.1%)
MMP9 expression
  Low62 (70.5%)26 (29.5%)<0.001--
  High27 (29.7%)64 (70.3%)--
Shallow: tumor invasion to mucosa and muscularis; Deep: tumor invasion to or beyond serosa; TNM, tumour node metastasis; N, lymph node; CEA, Carcinoembryonic antigen; CA199, carbohydrate antigen 199.

Rab1B and MMP9 overexpression is correlated with tumor progression and metastasis in CRC patients

To evaluate the clinical relevance of Rab1B and MMP9 protein expression in CRC patients, the median IHC scores of 180 and 156 were defined as the cutoff value for high- and low-expression of Rab1B and MMP9, respectively, which divided CRC patients into high- or low-expression groups. As showed in Table 2, high Rab1B and MMP9 protein expressions in CRC are significantly associated with deep invasion, lymph node metastasis, and advanced TNM stage, suggesting that both Rab1B and MMP9 are involved in the progression and metastasis of CRC.

Rab1B and MMP9 protein expressions are associated with poor prognosis

To explore the prognostic value of Rab1B and MMP9 proteins in this disease, we analyzed overall survival (OS) and progression-free survival (PFS) of CRC patients with Rab1B and/or MMP9 expressions. Patients with Rab1B high-expression have significantly shorter 5-year OS rate (63.5% vs. 92.6%) and 5-year PFS rate (56.3% vs. 88.1%) than those with Rab1B low-expression (All P < 0.001, Fig. 3A). Similarly, patients with MMP9 high-expression have much poorer OS and PFS than those with MMP9 low-expression (Fig. 3B), which is consistent with the conclusion from a meta-analysis of 13 cohort studies on MMP9 expression and prognosis in CRC [25]. These results indicate that the increased Rab1B or MMP9 protein is significantly correlated with poor prognosis of CRC patients.

Overexpression ofRab1B and MMP9 proteins are associated with poor prognosis of CRC patients independent of clinical stage. (A) The overall survival (OS) and progression-free survival (PFS) of CRC patients with high or low Rab1B expression. P value was calculated by Log-rank test. (B) OS and PFS of patients with high or low MMP9 expression. (C) OS and PFS of stage I-II CRC patients with high or low Rab1B expression. (D) OS and PFS of stage III CRC patients with high or low Rab1B expression.

Figure 3. Overexpression ofRab1B and MMP9 proteins are associated with poor prognosis of CRC patients independent of clinical stage. (A) The overall survival (OS) and progression-free survival (PFS) of CRC patients with high or low Rab1B expression. P value was calculated by Log-rank test. (B) OS and PFS of patients with high or low MMP9 expression. (C) OS and PFS of stage I-II CRC patients with high or low Rab1B expression. (D) OS and PFS of stage III CRC patients with high or low Rab1B expression.

We next analyzed the effect of their expressions and clinicopathological parameters on patient survival by using univariate and multivariate Cox model. Univariate analysis indicates that TNM stage, adjuvant chemotherapy, Rab1B and MMP9 proteins are significant predictors for OS and PFS of CRC patients (All P < 0.05, Table 3), and Tumor depth is a marginally significant predictor for OS and PFS (P =0.065 and P =0.059, respectively, Table 3). Multivariate Cox regression analysis further demonstrates that Rab1B protein is an independent risk predictor for OS (HR: 3.605, 95% CI: 1.481-8.775, P = 0.005) and PFS (HR: 3.394, 95% CI: 1.579-7.297, P = 0.002) of CRC patients (Table 3), and MMP9 protein is a marginally significant independent unfavorable predictor for OS and PFS in CRC patients (Table 3). In addition, TNM stage and adjuvant chemotherapy also are independent prognostic factors for OS and PFS in CRC patients.

Table 3. Cox regression analysis of Rab1B, MMP9 and clinical characteristics associated with survival in patients with colorectal cancer

VariablesUnivariate analysisP valueMultivariate analysisP value
HR95% CIHR95% CI
Overall Survival
Gender (Female vs. Male)0.8620.458-1.6230.645
Age (≥65y vs. <65y)1.1960.621-2.3010.593
Tumor location (Rectum vs. Colon)0.7780.413-1.4650.436
Tumor size (≥5cm vs. <5cm)1.1600.857-3.0230.139
Pathological grade (III-IV vs. I-II)1.5720.766-3.2260.218
Tumor depth (Deep vs. Shallow)6.4770.889-47.1760.0651.5410.187-12.7160.688
TNM stage (III vs. II vs. I)3.9712.028-7.778<0.0012.3561.123-4.9440.023
Intraoperative chemotherapy
  (Yes vs. No)
1.3480.700-2.5930.372
Adjuvant chemotherapy
  (Yes vs. No)
0.3750.172-0.8150.0130.3720.168-0.8240.015
Preoperative serum CEA
  (> 5 ng/ml vs. 0-5 ng/ml)
0.9370.46-1.9220.858
Preoperative serum CA199
  (> 35 ng/ml vs. 0-35 ng/ml)
1.1200.469-2.6740.798
Rab1B expression (High vs. Low)5.2742.327-11.952<0.0013.6051.481-8.7750.005
MMP9 expression (High vs. Low)4.4032.022-9.586<0.0012.0310.872-4.7290.101
Progression-Free Survival
Gender (Female vs. Male)0.9340.530-1.6490.815
Age (≥65y vs. <65y)1.5210.857-2.7010.152
Tumor location (Rectum vs. Colon)0.7390.416-1.3120.302
Tumor size (≥5cm vs. <5cm)1.6200.917-2.8590.096
Pathological grade(III-IV vs. I-II)1.1370.566-2.2810.719
Tumor depth (Deep vs. Shallow)3.9140.950-16.1250.0591.1350.240-5.3650.873
TNM stage (III vs. II vs. I)3.0831.762-5.394<0.0012.0221.080-3.7880.028
Intraoperative chemotherapy
  (Yes vs. No)
1.1710.643-2.1350.605
Adjuvant chemotherapy
  (Yes vs. No)
0.4280.218-0.8390.0130.4150.207-0.8310.014
Preoperative serum CEA
  (>5 ng/ml vs. 0-5 ng/ml)
1.1100.596-2.0690.742
Preoperative serum CA199
  (>35 ng/ml vs. 0-35 ng/ml)
1.2720.595-2.7180.535
Rab1B expression (High vs. Low)4.3562.169-8.747<0.0013.3941.579-7.2970.002
MMP9 expression (High vs. Low)3.0311.602-5.7350.0011.6080.793-3.2590.188
HR, hazard ratio; CI, confidence interval; Shallow: tumor invasion within mucosa and muscularis; Deep: tumor invasion to or beyond serosa; TNM, tumour node metastasis; CEA, Carcinoembryonic antigen; CA199, carbohydrate antigen 199.

Moreover, we found when patients were stratified by TNM stage, stage I-II patients with high level of Rab1B protein has marginally significantly poorer OS and PFS than those with low level (Fig. 3C), and stage III patients with high Rab1B level has significantly poorer survival than those with low level (Fig. 3D). This result reveals that Rab1B can predict survival of CRC patients independent of clinical stage. Therefore, Rab1B may provide additional prognostic information to the current clinical staging system.

Combination of Rab1B and MMP9 expression significantly improves predictive efficiency for the outcome of CRC patients

We found that in patients with low MMP9 expression, those with high Rab1B expression had significantly worse survival than those with low Rab1B expression (Fig. 4A). Similar result was observed in patients with MMP9 high expression (Fig. 4B). These results suggest that Rab1B is a prognostic predictor independent of MMP9. Therefore, we hypothesize that combination of both proteins will improve their predictive efficiency for survival. To this end, a new combined risk score was calculated as the sum of Rab1B score (0 or 1) and MMP9 score (0 or 1) for each case. Survival analysis shows that patients with low- (score 0), intermediate- (score 1) or high-risk (score 2) have significantly different 5-year overall survival rates, 96.7%, 80.2% and 57.5%, respectively. Multivariate Cox regression analysis reveals that the new combined risk score is an independent prognostic factor for OS and PFS in CRC patients (Table S1). Compared with TNM staging system, Rab1B/MMP9 combined risk score has some better predictive accuracy for OS (ROC area: 0.76 vs. 0.71, Fig. S1A) and PFS (ROC area: 0.73 vs. 0.69, Fig. S1B) of CRC patients though there is no statistically significant difference between the two ROC areas. Thus, the new combined risk score is a potentially useful biomarker that provides additional prognostic information for physician to evaluate survival and make decision on adjuvant chemotherapy in post-operative CRC patients.

Combined overexpression of Rab1B/MMP9 further improves predictive efficiency for outcome of CRC patients. (A) OS and PFS of patients with high or low Rab1B expression in those with low-expression of MMP9. (B) OS and PFS of patients with high or low Rab1B expression in those with high-expression of MMP9. (C) OS and PFS of patients who were stratified into three risk groups by the combined risk score of Rab1B and MMP9 proteins. Kaplan-Meier survival was used to predict the outcomes of CRC patients with low, intermediate or high combined risk score.

Figure 4. Combined overexpression of Rab1B/MMP9 further improves predictive efficiency for outcome of CRC patients. (A) OS and PFS of patients with high or low Rab1B expression in those with low-expression of MMP9. (B) OS and PFS of patients with high or low Rab1B expression in those with high-expression of MMP9. (C) OS and PFS of patients who were stratified into three risk groups by the combined risk score of Rab1B and MMP9 proteins. Kaplan-Meier survival was used to predict the outcomes of CRC patients with low, intermediate or high combined risk score.

Rab1B and MMP9 protein expression predicts outcome of adjuvant chemotherapy in CRC patients

In clinical practice, physicians determine the use of adjuvant chemotherapy based on the clinical stages and high-risk features of patients with CRC. However, clinical stages and high-risk features do not always accurately predict the outcome of chemotherapy. To investigate the utility of Rab1B and MMP9 proteins in CRC management, we explored the use of Rab1B and/or MMP9 protein expression to predict the response to adjuvant chemotherapy in CRC patients. First, we divided patients into high- and low-Rab1B groups. Survival analysis revealed that patients with low Rab1B protein had a similar outcome (including OS and PFS) regardless whether they had adjuvant chemotherapy (Fig. 5A), indicating that these patients do not benefit from adjuvant chemotherapy. Thus these patients should avoid the unnecessary toxicity and economic burden associated with chemotherapy. We next conducted the same analysis on patients with high Rab1B protein. Surprisingly, patients received adjuvant chemotherapy had markedly better OS and PFS than those not received the chemotherapy (Fig. 5B), clearly indicating that patients with high Rab1B protein level may obtain benefit from adjuvant chemotherapy. We also performed the same analysis on MMP9 and observed similar results (Fig. 5C and 5D). Finally, we combined Rab1B and MMP9 into the aforementioned new risk score model and conducted the same analysis. The results show that in patients with low or intermediate risk, there is no difference in OS and PFS between patients with and without chemotherapy (Fig. 6A and 6B), but in patients with high risk, those with chemotherapy had significantly better survival than those without chemotherapy (Fig. 6C).

Rab1B and MMP9 protein expression predicts outcome of adjuvant chemotherapy in CRC patients. Patients with CRC were stratified into high- or low-expression group by Rab1B or MMP9 expression. (A) Kaplan-Meier survival and Log-rank test were used to compare OS and PFS of CRC patients with or without adjuvant chemotherapy in low Rab1B expression group. (B) OS and PFS of CRC patients with or without adjuvant chemotherapy in the high Rab1B expression group. (C) OS and PFS of CRC patients with or without adjuvant chemotherapy in the low MMP9 expression group. (D) OS and PFS of patients with or without adjuvant chemotherapy in the high MMP9 expression group.

Figure 5. Rab1B and MMP9 protein expression predicts outcome of adjuvant chemotherapy in CRC patients. Patients with CRC were stratified into high- or low-expression group by Rab1B or MMP9 expression. (A) Kaplan-Meier survival and Log-rank test were used to compare OS and PFS of CRC patients with or without adjuvant chemotherapy in low Rab1B expression group. (B) OS and PFS of CRC patients with or without adjuvant chemotherapy in the high Rab1B expression group. (C) OS and PFS of CRC patients with or without adjuvant chemotherapy in the low MMP9 expression group. (D) OS and PFS of patients with or without adjuvant chemotherapy in the high MMP9 expression group.

Combined expression of Rab1B and MMP9 proteins predicts the outcome of adjuvant chemotherapy in CRC patients. CRC Patients were stratified into three risk groups by the combined risk score of Rab1B and MMP9 protein expression. Kaplan-Meier survival was used to compare OS and PFS of CRC patients with or without adjuvant chemotherapy in low risk group (A), intermediate risk group (B), and high risk group (C).

Figure 6. Combined expression of Rab1B and MMP9 proteins predicts the outcome of adjuvant chemotherapy in CRC patients. CRC Patients were stratified into three risk groups by the combined risk score of Rab1B and MMP9 protein expression. Kaplan-Meier survival was used to compare OS and PFS of CRC patients with or without adjuvant chemotherapy in low risk group (A), intermediate risk group (B), and high risk group (C).

According to the current criteria of adjuvant chemotherapy, 109 patients were not administered adjuvant chemotherapy in this study. When these patients were divided into low- and high-expression groups based on Rab1B or MMP9 expression, those with high-expression of Rab1B or MMP9 had significantly worse survival (including OS and PSF) than those with low-expression (Fig. S2A and S2B), indicating that patients with high-expression of Rab1B or MMP9 should be administered adjuvant chemotherapy to improve survival. Next, we combined both proteins to predict outcome of these 109 patients, and the result showed that patients with low combined risk had the best survival and those with high combined risk had the worst survival among the three risk subgroups (Fig. S2C). Finally, we compared single proteins with the combined model using ROC analysis. The result demonstrates that the combination risk score model is remarkably better than single proteins (Fig. S2D), suggesting that the combined risk model can provide more detailed and accurate information to help physicians make chemotherapy decisions.

Discussion

Rab1A and Rab1B are highly homologous and known to share common biological functions in ER to Golgi trafficking and autophagosome formation [26,27]. Our previous studies have demonstrated that Rab1A is an oncogene whose overexpression is correlated with poor prognosis in patients with CRC and HCC [10,11]. Hence, it seems plausible that Rab1B also plays an oncogenic role in CRC. In this study, as expected, we find that the expression of Rab1B mRNA and protein is up-regulated in most of CRC tissues and cell lines, which is consistent with the result recently reported by Zhai et al [17]. Furthermore, another recent study suggested that depletion of Rab1B suppresses tumor growth by inhibiting PI3K/AKT signaling pathway in the lung cancer cell line A549 [28]. However, Rab1B was also reported to be down-regulated in breast cancer and inhibits proliferation and metastasis of breast cancer cells [12]. The reported contradictory functions for Rab1B protein may be due to its distinct roles in a tissue origin- and tumor type-specific manner. Rab1B may participate in distinct signaling pathways, promoting either oncogenic or tumor-suppressing activity in different context.

At present, the most important clinical prognostic predictor of CRC patients is TNM stage. While our data support the clinical utility of the TNM staging system in predicting the prognosis of CRC patients, many CRC patients with same TNM stages have opposite clinical outcomes, suggesting that TNM staging system needs to be improved by additional prognostic factors such as molecular biomarkers for metastasis and relapse. In this study, we explored the prognostic value of Rab1B and MMP9 in CRC patients. Our results demonstrate that Rab1B is a significant prognostic factor independent of the TNM staging system, while MMP9 is a marginally significantly independent predictor for survival of CRC patients, which is consistent with several previous studies [22,29]. However, when Rab1B and MMP9 are combined into a new risk model, it provides much better prediction of survival in CRC patients. Altogether, our result suggests that the combined risk model is a useful biomarker for prognosis and can provide additional prognostic information in CRC patients.

In clinical practice, TNM staging system also is the main method for making decision on adjuvant chemotherapy. However, our data shows that the patients who are not administered chemotherapy based on the TNM stages have remarkably different survivals when they are stratified by Rab1B or MMP9 expression, indicating that the TNM stage is not adequate for making decision on adjuvant chemotherapy. In contrast, patients with low Rab1B or MMP9 expression show the similar survival rates no matter whether they received adjuvant chemotherapy or not, suggesting that these patients do not need chemotherapy. In patients with high Rab1B or MMP9 expression, those received chemotherapy have significantly better survival than those did not, indicating that patients with high expression of Rab1B or MMP9 should receive chemotherapy. Furthermore, when combination of Rab1B and MMP9 expression is used to predict survival of patients who did not receive adjuvant chemotherapy based on TNM stage, the combined risk model performs significantly better than Rab1B or MMP9 protein alone. Therefore, our study demonstrates for the first time that Rab1B and MMP9 individually or in combination are useful biomarkers for making decision on adjuvant chemotherapy.

In this study, we reveal that there is a significantly positive correlation between Rab1B and MMP9 expressions in CRC tissues, suggesting that Rab1B and MMP9 have an interaction. Interestingly, it was reported that MMP9 secretion is controlled by Rab1A mediated membrane trafficking [24], suggesting that Rab1B has a similar role on MMP9 secretion. Another recent study reported that PITPNC1 drives metastasis of multiple prevalent cancer types by enhancing Rab1B-mediated vesicular secretion of several pro-metastatic genes, which include MMP1 [30]. It is tempting to speculate that Rab1B regulates MMP9 in such fashion, which warrants further investigation.

In conclusion, our study reveals that the elevated expression of Rab1B and MMP9 is common event in CRC tissues and cell lines. Overexpression of Rab1B and MMP9 alone or together is significantly associated with poor prognosis, suggesting that they are useful biomarkers for prediction of outcome in CRC patients. More importantly, our study demonstrates for the first time that they can accurately predict outcome of adjuvant chemotherapy, which is valuable in guiding the chemotherapy for the post-operative CRC patients. Further studies are required to decipher the molecular mechanism by which Rab1B interacts with MMP9 to promote tumor progression and metastasis in CRC.

Materials and Methods

Patients and tissue samples

The tumor and matched colorectal tissues were obtained from 179 consecutive patients with stage I-III CRC who underwent radical resection at Sun Yat-Sen University Cancer Center (Guangzhou, China) between January 2009 and December 2010. None of patients received preoperative chemotherapy and/or radiotherapy. After radical resection, patients with stage III and high-risk stage II CRC further received postoperative adjuvant chemotherapy. Fluorouracil-containing regimens are standard, including FOLFOX, XELOX and Xeloda. All the samples were pathologically diagnosed by two experienced pathologists. Histological classification and tumor differentiation were determined according to the criteria of the World Health Organization. Clinical stage was defined according to the 7th edition of tumor-node metastasis (TNM) classification of the American Joint Committee on Cancer Staging (AJCC) [31]. This cohort of CRC patients included 95 males (52.8%) and 85 females (47.2%), with a mean age of 59.5 years old. During the follow-up, tumor assessment was done by colonoscopy, abdominal ultrasound or computed tomography scanning at 3-month intervals for the first 2 years, then at 6-month intervals for 3 years for a total of 5 years. The average follow-up time was 54.2 months, ranged from 1 to 80 months. Overall survival (OS) is defined as the time from the date of surgery to the date of death from any cause or last date of follow-up; Progression-free survival (PFS) is defined as the time from the date of surgery to the date of relapse or metastasis of CRC or death from CRC or last date of follow-up.

Another randomly selected eight pairs of fresh CRC tissues and matched adjacent non-tumorous colorectal tissues from patients undergoing surgical resection in 2011 were collected for immunoblot analysis. This study was reviewed and approved by the Ethical Committees of Sun Yat-Sen University Cancer Center. Written informed consent was obtained from all patients before surgery.

Cell lines and cell culture

Eleven human CRC cell lines (SW480, SW620, HT29, CACO2, HCT8, HCT116, CW2, DLD-1, RKO, LS174T and LoVo) and human immortalized colon cell lines (CCD-112CoN) were obtained from the American Type Culture Collection (ATCC). All cell lines were cultured in the conditions specified by the manufacturer, and authenticated by short tandem repeat DNA fingerprinting and tested for mycoplasma before use at Medicine Laboratory of the Department of Forensic Medicine, Sun Yat-Sen University (Guangzhou, China).

Immunohistochemistry (IHC)

IHC staining was performed using standard streptavidin-peroxidase complex method as described previously [32]. Briefly, fresh surgical tissue specimens were fixed in 10% formaldehyde and routinely processed for paraffin embedding. Then these blocks were cut into 4μm thick sections. Endogenous peroxidase activity was blocked by 3% hydrogen peroxide in methanol. After antigen retrieval, the sections were incubated with rabbit polyclonal Rab1B antibody (1:200, 17824-1-AP,Proteintech Group, USA) or mouse monoclonal MMP9 antibody (1:50, sc-21733,Santa Cruz, USA) overnight at 4°C and a negative control was set up by replacing the primary antibody with phosphate buffered solution (PBS). Subsequently, the sections were incubated with horseradish peroxidase (HRP)-conjugated anti-rabbit or anti-mouse IgG secondary antibody (Dako, Denmark) for 30 min at room temperature, followed by developing using 3, 5-diaminobenzidine (DAB, Dako, Denmark) substrate and counterstaining with Mayer’s hematoxylin for the nuclei.

Scoring of immunostaining

The scoring of immunostaining was independently performed by two experienced pathologists who were blinded to the clinical information. The immunostaining intensity was scored as 0 - 3 (0, no staining; 1, weak staining; 2, moderate staining; 3, strong staining), and percentage of immunostained tumor cells was scored as 0 - 100. The immunostaining intensity was multiplied by the percentage of stained tumor cells, which resulted in a semi-quantitative immunostaining score (ISS) between 0 and 300.

Immunoblotting

The fresh CRC tissue samples and CRC cell lines were directly homogenized in RIPA Lysis buffer. 30μg of lysate protein were run on a 10% SDS-PAGE gel and then transferred onto polyvinylidene fluoridemembranes (PVDF, Millipore, USA). After blocking nonspecific binding site with TBST containing 5% non-fat milk, the membranes were incubated with a rabbit polyclonal anit-Rab1B antibody (1:300, sc-599, Santa Cruz, USA), a mouse monoclonal anti-MMP9 antibody (1:200, sc-21733, Santa Cruz, USA) or GAPDH antibody (1:6000, #5174, Cell Signaling Technology, USA) at 4°C overnight. Then the membranes were incubated with HRP-conjugated secondary antibody (1:5000, Jackson Immunoresearch Inc, USA) for about 60 min at room temperature. The blots were scanned and the intensities of protein bands were quantitated by the Bio-Rad software Quantity One (Bio-Rad Laboratories Inc., USA).

Isolation of total RNA and Real-time quantitative PCR (RT-PCR)

Total RNA from cancer cells was isolated using TRIzol reagent (Invitrogen) as per the manufacturer’s instruction. RNA samples were purified and extracted with phenol and chloroform. The quantity of RNA samples was measured by a NanoDropTM 2000 spectrophotometer (ThermoFisher Scientific, Waltham, MA, USA). One microgram of total RNA was used to synthesize cDNA via reverse transcription reaction according to the protocol of GoScriptTM Reverse Transcription System kit (Promega, A5001, Madison, WI, USA). After 1:20 dilution, 2μl of cDNA products was employed for PCR, which was performed on a Roche Lightcycler 96 real time PCR machine (Roche Diagnostics, Indianapolis, IN, USA) according to a standard method as described previously [33]. All samples were amplified in triplicate and GAPDH was detected as an internal control. The relative quantification of target genes was calculated using the comparative 2-ΔΔCT method [34]. Primer sequences for this experiment are as follows:Rab1B forward primer: 5’-GGACTTCAAGATCCGAACCAT-3’, reverse primer: 5’-ATACACCACGATGATGCCA-3’, and the amplicon length: 135bp; MMP9 forward primer: 5’-GGGACGCAGACATCGTCATC-3’, reverse primer: 5’-TCGTCATCGTCGAAATGGGC-3’, and the amplicon length: 139bp; GAPDH forward primer: 5’-CTCCTCCTGTTCGACAGTCAGC-3’, reverse primer: 5’-CCCAATACGACCAAATCCGTT-3’, and the amplicon length: 204bp.

Statistical analysis

SPSS 17.0 for windows (SPSS, Inc., Chicago, IL, USA) and GraphPad PrismV6 (GraphPad Prism, Inc., USA) were used for statistical analyses. The results were expressed as mean ± SD or SEM. The correlation between Rab1B expression and clinicopathological parameters was analyzed by Chi-square test or Fisher’s exact test. The Student’s t-test was used for the analysis of comparisons. The relationship between Rab1B and MMP9 expressions was tested by Pearson’s correlation analysis. The prognostic variables were evaluated by univariate and multivariate Cox proportional hazards regression model. Kaplan-Meier plots and log-rank test were used to estimate cancer specific survival curves.

Supplementary Materials

Supplementary File

Conflicts of Interest

There is no conflict of interest disclosures from any authors.

Funding

This work was supported by the Research Fund of State Key Laboratory of Oncology in South China (No. 010811).

References

  • 1. Brenner H, Kloor M, Pox CP. Colorectal cancer. Lancet. 2014; 383:1490–502. https://doi.org/10.1016/S0140-6736(13)61649-9 [PubMed]
  • 2. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015; 65:87–108. https://doi.org/10.3322/caac.21262 [PubMed]
  • 3. Lin GN, Liu PP, Liu DY, Peng JW, Xiao JJ, Xia ZJ. Prognostic significance of the pre-chemotherapy lymphocyte-to-monocyte ratio in patients with previously untreated metastatic colorectal cancer receiving FOLFOX chemotherapy. Chin J Cancer. 2016; 35:5. https://doi.org/10.1186/s40880-015-0063-1 [PubMed]
  • 4. Barton MK. Oxaliplatin in the adjuvant treatment of colon cancer. CA Cancer J Clin. 2012; 62:3–4. https://doi.org/10.3322/caac.21131 [PubMed]
  • 5. Li P, Fang YJ, Li F, Ou QJ, Chen G, Ma G. ERCC1, defective mismatch repair status as predictive biomarkers of survival for stage III colon cancer patients receiving oxaliplatin-based adjuvant chemotherapy. Br J Cancer. 2013; 108:1238–44. https://doi.org/10.1038/bjc.2013.83 [PubMed]
  • 6. Oue N, Anami K, Schetter AJ, Moehler M, Okayama H, Khan MA, Bowman ED, Mueller A, Schad A, Shimomura M, Hinoi T, Aoyagi K, Sasaki H, et al. High miR-21 expression from FFPE tissues is associated with poor survival and response to adjuvant chemotherapy in colon cancer. Int J Cancer. 2014; 134:1926–34. https://doi.org/10.1002/ijc.28522 [PubMed]
  • 7. Chaffer CL, Weinberg RA. A perspective on cancer cell metastasis. Science. 2011; 331:1559–64. https://doi.org/10.1126/science.1203543 [PubMed]
  • 8. Hutagalung AH, Novick PJ. Role of Rab GTPases in membrane traffic and cell physiology. Physiol Rev. 2011; 91:119–49. https://doi.org/10.1152/physrev.00059.2009 [PubMed]
  • 9. Yang XZ, Li XX, Zhang YJ, Rodriguez-Rodriguez L, Xiang MQ, Wang HY, Zheng XF. Rab1 in cell signaling, cancer and other diseases. Oncogene. 2016; 35:5699–704. https://doi.org/10.1038/onc.2016.81 [PubMed]
  • 10. Thomas JD, Zhang YJ, Wei YH, Cho JH, Morris LE, Wang HY, Zheng XF. Rab1A is an mTORC1 activator and a colorectal oncogene. Cancer Cell. 2014; 26:754–69. https://doi.org/10.1016/j.ccell.2014.09.008 [PubMed]
  • 11. Xu BH, Li XX, Yang Y, Zhang MY, Rao HL, Wang HY, Zheng XF. Aberrant amino acid signaling promotes growth and metastasis of hepatocellular carcinomas through Rab1A-dependent activation of mTORC1 by Rab1A. Oncotarget. 2015; 6:20813–28. https://doi.org/10.18632/oncotarget.5175 [PubMed]
  • 12. Wang F, Zheng Z, Guo J, Ding X. Correlation and quantitation of microRNA aberrant expression in tissues and sera from patients with breast tumor. Gynecol Oncol. 2010; 119:586–93. https://doi.org/10.1016/j.ygyno.2010.07.021 [PubMed]
  • 13. Wang X, Liu F, Qin X, Huang T, Huang B, Zhang Y, Jiang B. Expression of Rab1A is upregulated in human lung cancer and associated with tumor size and T stage. Aging (Albany NY). 2016; 8:2790–98. https://doi.org/10.18632/aging.101087 [PubMed]
  • 14. Nikoshkov A, Broliden K, Attarha S, Sviatoha V, Hellström AC, Mints M, Andersson S. Expression pattern of the PRDX2, RAB1A, RAB1B, RAB5A and RAB25 genes in normal and cancer cervical tissues. Int J Oncol. 2015; 46:107–12. [PubMed]
  • 15. Nalesnik MA, Tseng G, Ding Y, Xiang GS, Zheng ZL, Yu Y, Marsh JW, Michalopoulos GK, Luo JH. Gene deletions and amplifications in human hepatocellular carcinomas: correlation with hepatocyte growth regulation. Am J Pathol. 2012; 180:1495–508. https://doi.org/10.1016/j.ajpath.2011.12.021 [PubMed]
  • 16. Jiang HL, Sun HF, Gao SP, Li LD, Hu X, Wu J, Jin W. Loss of RAB1B promotes triple-negative breast cancer metastasis by activating TGF-β/SMAD signaling. Oncotarget. 2015; 6:16352–65. https://doi.org/10.18632/oncotarget.3877 [PubMed]
  • 17. Zhai H, Song B, Xu X, Zhu W, Ju J. Inhibition of autophagy and tumor growth in colon cancer by miR-502. Oncogene. 2013; 32:1570–79. https://doi.org/10.1038/onc.2012.167 [PubMed]
  • 18. Matsuyama Y, Takao S, Aikou T. Comparison of matrix metalloproteinase expression between primary tumors with or without liver metastasis in pancreatic and colorectal carcinomas. J Surg Oncol. 2002; 80:105–10. https://doi.org/10.1002/jso.10106 [PubMed]
  • 19. Lee MA, Park JH, Rhyu SY, Oh ST, Kang WK, Kim HN. Wnt3a expression is associated with MMP-9 expression in primary tumor and metastatic site in recurrent or stage IV colorectal cancer. BMC Cancer. 2014; 14:125. https://doi.org/10.1186/1471-2407-14-125 [PubMed]
  • 20. Kim HC, Kim YS, Oh HW, Kim K, Oh SS, Kim JT, Kim BY, Lee SJ, Choe YK, Kim DH, Kim SH, Chae SW, Kim KD, Lee HG. Collagen triple helix repeat containing 1 (CTHRC1) acts via ERK-dependent induction of MMP9 to promote invasion of colorectal cancer cells. Oncotarget. 2014; 5:519–29. https://doi.org/10.18632/oncotarget.1714 [PubMed]
  • 21. Bendardaf R, Buhmeida A, Hilska M, Laato M, Syrjänen S, Syrjänen K, Collan Y, Pyrhönen S. MMP-9 (gelatinase B) expression is associated with disease-free survival and disease-specific survival in colorectal cancer patients. Cancer Invest. 2010; 28:38–43. https://doi.org/10.3109/07357900802672761 [PubMed]
  • 22. Zhang Y, Guan XY, Dong B, Zhao M, Wu JH, Tian XY, Hao CY. Expression of MMP-9 and WAVE3 in colorectal cancer and its relationship to clinicopathological features. J Cancer Res Clin Oncol. 2012; 138:2035–44. https://doi.org/10.1007/s00432-012-1274-3 [PubMed]
  • 23. Goos JA, Coupé VM, van de Wiel MA, Diosdado B, Delis-Van Diemen PM, Hiemstra AC, de Cuba EM, Beliën JA, Menke-van der Houven van Oordt CW, Geldof AA, Meijer GA, Hoekstra OS, Fijneman RJ, and DeCoDe PET Group. A prognostic classifier for patients with colorectal cancer liver metastasis, based on AURKA, PTGS2 and MMP9. Oncotarget. 2016; 7:2123–34. [PubMed]
  • 24. Hwang YS, Park KK, Chung WY. Kalopanaxsaponin A inhibits the invasion of human oral squamous cell carcinoma by reducing metalloproteinase-9 mRNA stability and protein trafficking. Biol Pharm Bull. 2012; 35:289–300. https://doi.org/10.1248/bpb.35.289 [PubMed]
  • 25. Li CY, Yuan P, Lin SS, Song CF, Guan WY, Yuan L, Lai RB, Gao Y, Wang Y. Matrix metalloproteinase 9 expression and prognosis in colorectal cancer: a meta-analysis. Tumour Biol. 2013; 34:735–41. https://doi.org/10.1007/s13277-012-0601-2 [PubMed]
  • 26. Szatmári Z, Sass M. The autophagic roles of Rab small GTPases and their upstream regulators: a review. Autophagy. 2014; 10:1154–66. https://doi.org/10.4161/auto.29395 [PubMed]
  • 27. Martinez H, García IA, Sampieri L, Alvarez C. Spatial-Temporal Study of Rab1b Dynamics and Function at the ER-Golgi Interface. PLoS One. 2016; 11:e0160838. https://doi.org/10.1371/journal.pone.0160838 [PubMed]
  • 28. Wheeler DB, Zoncu R, Root DE, Sabatini DM, Sawyers CL. Identification of an oncogenic RAB protein. Science. 2015; 350:211–17. https://doi.org/10.1126/science.aaa4903 [PubMed]
  • 29. Araújo RFJr, Lira GA, Vilaça JA, Guedes HG, Leitão MC, Lucena HF, Ramos CC. Prognostic and diagnostic implications of MMP-2, MMP-9, and VEGF-α expressions in colorectal cancer. Pathol Res Pract. 2015; 211:71–77. https://doi.org/10.1016/j.prp.2014.09.007 [PubMed]
  • 30. Halberg N, Sengelaub CA, Navrazhina K, Molina H, Uryu K, Tavazoie SF. PITPNC1 Recruits RAB1B to the Golgi Network to Drive Malignant Secretion. Cancer Cell. 2016; 29:339–53. https://doi.org/10.1016/j.ccell.2016.02.013 [PubMed]
  • 31. Ueno H, Mochizuki H, Akagi Y, Kusumi T, Yamada K, Ikegami M, Kawachi H, Kameoka S, Ohkura Y, Masaki T, Kushima R, Takahashi K, Ajioka Y, et al. Optimal colorectal cancer staging criteria in TNM classification. J Clin Oncol. 2012; 30:1519–26. https://doi.org/10.1200/JCO.2011.39.4692 [PubMed]
  • 32. Huang Y, Zhang J, Cui ZM, Zhao J, Zheng Y. Expression of the CXCL12/CXCR4 and CXCL16/CXCR6 axes in cervical intraepithelial neoplasia and cervical cancer. Chin J Cancer. 2013; 32:289–96. https://doi.org/10.5732/cjc.012.10063 [PubMed]
  • 33. Cha YL, Li PD, Yuan LJ, Zhang MY, Zhang YJ, Rao HL, Zhang HZ, Zheng XF, Wang HY. EIF4EBP1 overexpression is associated with poor survival and disease progression in patients with hepatocellular carcinoma. PLoS One. 2015; 10:e0117493. https://doi.org/10.1371/journal.pone.0117493 [PubMed]
  • 34. Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001; 29:e45. https://doi.org/10.1093/nar/29.9.e45 [PubMed]