Upregulation of SLC2A3 gene and prognosis in colorectal carcinoma: analysis of TCGA data
BMC Cancer volume 19, Article number: 302 (2019)
Upregulation of SLC2A genes that encode glucose transporter (GLUT) protein is associated with poor prognosis in many cancers. In colorectal cancer, studies reporting the association between overexpression of GLUT and poor clinical outcomes were flawed by small sample sizes or subjective interpretation of immunohistochemical staining. Here, we analyzed mRNA expressions in all 14 SLC2A genes and evaluated the association with prognosis in colorectal cancer using data from the Cancer Genome Atlas (TCGA) database.
In the present study, we analyzed the expression of SLC2A genes in colorectal cancer and their association with prognosis using data obtained from the TCGA for the discovery sample, and a dataset from the Gene Expression Omnibus for the validation sample.
SLC2A3 was significantly associated with overall survival (OS) and disease-free survival (DFS) in both the discovery sample (345 patients) and validation sample (501 patients). High SLC2A3 expression resulted in shorter OS and DFS. In multivariate analyses, high SLC2A3 levels predicted unfavorable OS (adjusted HR 1.95, 95% CI 1.22–3.11; P = 0.005) and were associated with poor DFS (adjusted HR 1.85, 95% CI 1.10–3.12; P = 0.02). Similar results were found in the discovery set.
Upregulation of the SLC2A3 genes is associated with decreased OS and DFS in colorectal cancer patients. Therefore, assessment of SLC2A3 gene expression may useful for predicting prognosis in these patients.
Colorectal cancer (CRC) is the third most common cancer and the fourth-leading cause of cancer death in the world [1, 2]. Most CRCs originate from non-cancerous lesions by one or a combination of three different mechanisms: chromosomal instability, CpG island methylator phenotype, and microsatellite instability . Biomarkers of these cytogenetic alterations are of interest for diagnosis, prognostication, and anticancer drug development targeting CRC . Despite research efforts, genetic biomarkers currently have limited value as diagnostic or prognostic markers .
Among the biomarkers, the solute carrier 2A (SLC2A) gene family that encodes glucose transporter (GLUT) proteins has been widely investigated. GLUT proteins facilitate glucose influx into cancer cells which is necessary for cancer cell proliferation. Upregulation of SLC2A genes is associated with poor prognosis in many cancers, including hepatocellular carcinoma, non–small cell lung cancer, and thyroid carcinoma [6,7,8].
An association between overexpression of the subtypes of GLUT proteins and poor clinical outcomes has been reported in CRC . However, these studies were flawed by small sample sizes or subjective interpretation of immunohistochemical (IHC) staining. In this study, we analyzed the mRNA expression of all 14 SLC2A genes (corresponding to 14 GLUT proteins) and evaluated the associations with prognosis in CRC using data from the Cancer Genome Atlas (TCGA) database.
The TCGA CRC data was downloaded from cBioPortal for Cancer Genomics (http://www.cbioportal.org/). The dataset contains survival data with clinical information, somatic mutations, and mRNA expression counts. For validation, we obtained independent microarray datasets (GSE39582) from the Gene Expression Omnibus (GEO). mRNA counts of the validation set were measured by the Affymetrix Human Genome U133 Plus 2.0 Array in a log2 scale. Gene expression of the discovery set was measured by Illumina HiSeq platform and transformed into log2 scale. According to the publication guidelines, the datasets may be used for publication without restriction or limitation (https://cancergenome.nih.gov/publications/publicationguidelines, https://www.ncbi.nlm.nih.gov/geo/info/disclaimer.html).
The putative associations between conventional clinical pathology parameters (age at diagnosis, sex, American Joint Committee on Cancer (AJCC) 7th edition TNM stage, microsatellite instability, and mutational status of KRAS (v-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog) and BRAF (v-Raf murine sarcoma viral oncogene homolog B1 genes) and survival outcome were assessed by Chi-square tests. The Cox proportional-hazards model was used to identify genes associated with survival and to estimate mortality hazard ratios (HRs). The optimal cut-off points for SLC2A3 expression used to divide patients into low-risk and high-risk groups were determined using the MaxStat package of R software (Maximally selected Rank Statistics). Maxstat computes the maximally selected log-rank statistic to identify the cutpoint which provides the best separation (in which the standardized statistics take their maximum) into two groups. Kaplan-Meier analysis was performed to estimate the survival curves of the different subgroups and the log-rank test (Mantel–Cox) was used to compare the curve. Statistical analyses were performed using R statistical software (version 3.4.1) . All P-value were two sided. A P-value less than 0.05 was considered statistically significant.
Our study sample comprised 846 patients. Patient characteristics of the discovery and validation set are shown in Table 1. The TCGA sample of 345 patients were the discovery set, and the GSE39582 sample of 501 patients were the validation set.
Association between clinical parameters and survival outcome
The associations between clinical variables and OS and DFS in the discovery set are summarized in Table 2. Age > 65 was associated with worse OS compared to age ≤ 65 (P < 0.001). TNM stage III and IV was associated with OS (P < 0.001) and DFS (P < 0.001). Other clinicopathological factors (sex, MSI, BRAF, or KRAS status) were not associated with OS or DFS.
Prognostic value of SLC2A3
Table 3 displays the associations between the mRNA expression values of SLC2 family genes and survival outcomes. SLC2A3 was significantly associated with both OS (P = 0.013) and DFS (P = 0.014). There were associations between the expression of SLC2A6 and SLC2A7 with worse OS (P = 0.048 and 0.019); and SLC2A1 with worse DFS (P = 0.018).
Patients were categorized into high and low SLC2A3 expression groups according to the cut-off value determined by Maxstat method. The high expression group had worse OS (P = 0.005) and DFS (P = 0.002) compared to the low expression group (Fig. 1).
We explored the association between SLC2A3 expression and survival outcome in a multivariate context (Table 4). In multivariate Cox regression analysis, after adjusting for clinical factors that were significantly associated with survival outcome (OS: age > 65, AJCC stage III and IV, DFS; AJCC stage III and IV), SLC2A3 remained an independent poor prognostic factor for OS and DFS. SLC2A3 level is associated with shorter OS with an adjusted hazard ratio of 1.98 (95% CI: 1.24–3.17; P = 0.004) and with poor disease-free survival (adjusted HR: 1.85, 95% CI: 1.10–3.12; P = 0.020).
Validation set analysis
SLC2A3 was significantly associated with both OS (P = 0.005) and DFS (P = 0.024). There was associations between the expression of SLC2A1 with worse DFS (P = 0.015), but SLC2A6 was not associated with worse OS (P = 0.940). The expression of SLC2A7 was not provided. Patients in validation set were categorized into high vs. low SLC2A3 expression according to the cut-off point. High expression led to worse OS (P = 0.003) and DFS (P = 0.021) (Fig. 2). In the multivariate Cox regression analysis, SLC2A3 expression is associated with shorter OS (adjusted HR 1.50, 95% CI: 1.11–2.03; P = 0.009) and DFS (adjusted HR 1.38, 95% CI: 1.00–1.91; P = 0.048) after adjusting for clinical factors (age, AJCC stage for OS, AJCC stage for DFS) (Table 5).
Prognostic value of SLC2A1
Patients were categorized into high and low SLC2A1 expression groups according to the cut-off value determined by Maxstat method. The high expression group had worse DFS in both Discovery set and Validation set (P = 0.001 and < 0.001, respectively) compared to the low expression group (Fig. 3). We explored the association between SLC2A1 expression and DFS in a multivariate context of Discovery set and Validation set. In multivariate Cox regression analysis, after adjusting for AJCC TNM stage III and IV that were significantly associated with DFS, SLC2A1 remained an independent poor prognostic factor for DFS in both Discovery set (adjusted HR 1.83, 95% CI: 1.12–3.01; P = 0.017) and Validation set (adjusted HR 1.65, 95% CI: 1.14–2.40; P = 0.009).
Most cancer cells favor glycolytic energy metabolism over mitochondrial metabolism and oxidative phosphorylation chain for energy production, even in the presence of oxygen (Warburg effect). This explains why malignant cells overexpress GLUT family proteins, which is a plasma membrane transport system. Glucose can be translocated into the cell only via GLUT proteins. Expression and subcellular distribution of GLUT proteins are regulated by different signaling molecules and pathways such as PI3K, HIF, p53, Myc, and AMPK LBk1 [11, 12].
In the current study, we used the TCGA database to explore the clinical significance of SLC2 family genes in CRC. IHC staining of GLUT expression in cancer cells can be diverse which is a downside with IHC staining. This may be caused by incorrectly interpreted IHC stain results or different GLUT positive thresholds . Messenger RNA-gene expression analysis from TCGA data is superior to IHC and may best predict cancer prognosis in TCGA data . The current study is the first to report the expression of SLC2A genes in CRC using the TCGA database.
In the current study, age and AJCC TNM stage were associated with survival in CRC patient. As expected, in this study, age > 65 was linked with worse prognosis for CRC patients. This may be due to aging itself, or due to co-morbidities such as cardiovascular disease or postoperative complications . Older patients are less likely to be treated with resectional surgery than younger patients and have poorer survival outcomes . We included age > 65 as a factor for multivariate analysis because age > 65 is useful for prognostication. We found that AJCC stage III and IV, which has nodal metastasis or organ metastasis, was associated with OS and DFS in the TCGA database. AJCC TNM stage was established based on OS, and lymph node involvement in itself is considered to have a strong influence on OS and DFS in CRC . In oral squamous cell carcinoma and papillary thyroid carcinoma, GLUT3 was known to be a prognostic marker for OS, and was associated with advanced cancer stage (AJCC TNM stage III and IV) which has nodal metastasis [8, 17]. In current study, we found high expression of SLC2A3 was a prognostic factor for predicting OS and DFS in CRC and was not associated with AJCC TNM stage. Similar to current findings, GLUT3 is a significant marker of poor prognosis in laryngeal carcinoma, with no significant differences in nodal or distant metastasis between the GLUT3 negative and GLUT3 positive groups . Since SLC2A3 or GLUT3 genes are associated with prognosis in CRC as well as thyroid or laryngeal carcinoma, it would be worthwhile to investigate whether these cancers have similar GLUT-dependent metabolic pathways.
In the current study, we analyzed the mRNA expression values of all 14 SLC2A family genes and found that SLC2A3 is independently associated with both OS and DFS after adjusting for age and AJCC TNM stage in CRC patients. Among the GLUT family, GLUT1 is associated with prognosis in several cancers, including CRC. GLUT3 has high affinity glucose uptake, similar to GLUT1 . GLUT3 is overexpressed in human carcinomas including CRC . Although GLUT1 and GLUT3 have many similarities, including expression pathways, the effects of GLUT3 on outcomes in several cancer varieties are not as well understood as GLUT1. This is the first report to study the effects of SLC2A3 on CRC patients’ survival.
The following reasons may explain why the expression of SLC2A3 affects CRC prognosis. GLUT3 is induced by hypoxia-inducible factor (HIF) formation in response to hypoxia in carcinomas . HIF-1 promotes tumor metastasis into distant and more oxygenated tissue through the transcriptional activation of oncogenic growth factors such as transforming growth factor beta3, epidermal growth factor, and others . Solid tumors with high hypoxia levels are more malignant, more likely to metastasize, and have a worse prognosis . GLUT3 is also induced by Akt involved in the Warburg effect. In cancer cells, Akt increases expression of GLUT1 and GLUT3 by causing degradation of p53.It may reflect the activity of hypoxia independent oncogenic pathways [11, 19].
SLC2A1 expression was associated with poor DFS but not with OS. This corresponds well with a study which reported that GLUT-1 expression is associated with poor DFS but not with OS in rectal cancer patients . To further evaluate the association between the expression of SLC2A1 and rectal cancer, we performed a subgroup analysis for rectal cancer patients using a univariate Cox regression analysis. Although we were unable to confirm our findings in the validation set due to a lack of primary site information, in the discovery set, SLC2A1 expression in rectal cancer patients (n = 72) was significantly associated with DFS (HR 1.57, 95% CI: 1.04–2.38; P = 0.03). We also found no association between SLC2A1 expression and OS (HR 1.54, 95% CI: 0.96–2.46; P = 0.07) in rectal cancer patients. Our findings were consistent with the findings of the previous study.
In addition to SLC2A family genes, the following 12 genes are known to be involved in glucose metabolism: MTOR, RICTOR, HIF1A, MYC, PDK1, PDK2, PDK3, PDK4, PIK3R1, PKM, POU2F-1, and RPTOR . We conducted an analysis to evaluate if the expression of the genes was associated with survival outcomes of CRC. We found no association between the other 12 genes’ mRNA expression and survival outcomes of CRC in both discovery and validation sets (Additional file 1: Table S1). Based on these findings, we concluded that among the glucose metabolism regulating genes, only the SLC2A3 gene is significantly associated with the survival outcomes of CRC.
There are several genes known to be associated with the prognosis of CRC: BRAF, KRAS, HIF, TP53 and thymidylates synthase (TYMS) [22, 23]. The BRAF mutation generates an abnormality in the MEK/ERK signaling pathway in CRC  and has been reported to be associated with poor prognosis by many CRC studies. However, reports vary on its association with survival . Mutation of KRAS, a proto-oncogene, activates RAS signaling pathways, but its association with CRC survival is not clear [24, 26]. In the current study, we analyzed the association of BRAF and KRAS expression level with survival, and found neither to be associated with OS or DFS in the discovery set (Additional file 1: Table S2, Table S3 and Additional file 2: Figure S1). Likewise, in the recent study conducted by authors’ group, the mutational status of BRAF or KRAS was not associated the prognosis of CRC . For HIF, TP53 and TYMS, there are studies which reported the expression of the genes are not related to the prognosis of OS of CRC [28,29,30]. In this study, we also found that HIF, TP53 and TYMS were not associated with OS or DFS (Additional file 1: Table S4). The discrepancies among the studies regarding the prognostic value of these genes may result from differences in patient cohorts, available co-variates, or statistical methods.
The short observation period of the discovery set patients (31.4 months median follow-up period) is a limitation of our study. Secondly, the relationship between the expression of SLC2A mRNA and GLUT protein in CRC has not been confirmed. Although one study reported a close correlation between the expression of SLC2A1/SLC2A3 mRNA and that of GLUT1/GLUT3 proteins in thyroid carcinoma , further studies are needed to investigate whether the expression of SLC2A mRNA correlates with expression of GLUT protein mRNA in CRC. Another limitation of our results is that mRNA gene expression value is not a readily available parameter, especially in clinical settings, due to the high cost of storage and processing of fresh tissue. Its application may become wider in the near future when cost is reduced and stable mRNA expression can be obtained through formalin-fixed paraffin-embedded tissue samples.
In conclusion, upregulation of the SLC2A3 gene is associated with decreased OS and DFS in CRC patients. SLC2A3 gene expression analysis may be useful for predicting prognosis and survival of CRC patients.
American Joint Committee on Cancer
- BRAF :
v-Raf murine sarcoma viral oncogene homolog B1 genes
- KRAS :
v-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog
- SLC2A :
Solute carrier 2A
the cancer genome atlas
Ferlay J, Soerjomataram I, Ervik M, et al. GLOBOCAN 2012 v1.0, Cancer incidence and Mortaligy worldwide: IARC Cancer Base no. 11. Lyon, France: International Agency for Research on. Cancer. 2013.
Arnold M, Sierra MS, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global patterns and trends in colorectal cancer incidence and mortality. Gut. 2017;66(4):683–91.
Tariq K, Chias K. Colorectal cancer carcinogenesis: a review of mechanisms. Cancer Biol Med. 2016;13(1):120–35.
Lech G, Słotwiński R, Słodkowski M, et al. Colorectal cancer tumour markers and biomarkers: recent therapeutic advances. World J Gastroenterol 2016 7;22(5):1745–1755.
Luo HY, Xu RH. Predictive and prognostic biomarkers with therapeutic targets in advanced colorectal cancer. World J Gastroenterol. 2014;20(14):3858–74.
Amann T, Maegdefrau U, Hartmann A, et al. GLUT1 expression is increased in hepatocellular carcinoma and promotes tumorigenesis. Am J Pathol. 2009;174(4):1544–52.
Younes M, Brown RW, Stephenson M, Gondo M, Cagle PT. Overexpression of Glut1 and Glut3 in stage I nonsmall cell lung carcinoma is associated with poor survival. Cancer. 1997;80(6):1046–51.
Chai YJ, Yi JW, Oh SW, et al. Upregulation of SLC2 (GLUT) family genes is related to poor survival outcomes in papillary thyroid carcinoma: analysis of data from the Cancer genome atlas. Surgery. 2017;161(1):188–94.
Yang J, Wen J, Tian T, et al. GLUT-1 overexpression as an unfavorable prognostic biomarker in patients with colorectal cancer. Oncotarget. 2017;8(7):11788–96.
Therneau, R. A package for survival analysis in S. R package version 2.41-3. 2017 Available at: http://CRAN.R-project.org/package=survival.
Barron CC, Bilan PJ, Tsakiridis T, Tsiani E. Facilitative glucose transporters: implications for cancer detection, prognosis and treatment. Metabolism. 2016;65(2):124–39.
Cairns RA, Harris IS, Mak TW. Regulation of cancer cell metabolism. Nat Rev Cancer. 2011;11(2):85–95.
Zhao Q, Shi X, Xie Y, Huang J, Shia B, Ma S. Combining multidimensional genomic measurements for predicting cancer prognosis: observations from TCGA. Brief Bioinform. 2015;16(2):291–303.
Aquina CT, Mohile SG, Tejani MA, et al. The impact of age on complications, survival, and cause of death following colon cancer surgery. Br J Cancer. 2017;116(3):389–97.
Majano SB, Girolamo CD, Rachet B, et al. Surgical treatment and survival from colorectal cancer in Denmark, England, Norway, and Sweden: a population-based study. Lancet Oncol. 2018. https://doi.org/10.1016/S1470-2045(18)30646-6.
Peyravian N, Larki P, Gharib E, et al. The application of gene expression profiling in predictions of occult lymph node metastasis in colorectal Cancer patients. Biomedicines. 2018;6(1):27.
Ayala FR, Rocha RM, Carvalho KC, et al. GLUT1 and GLUT3 as potential prognostic markers for Oral squamous cell carcinoma. Molecules. 2010;15(4):2374–87.
Baer S, Casaubon L, Schwartz MR, Marcogliese A, Younes M. Glut3 expression in biopsy specimens of laryngeal carcinoma is associated with poor survival. Laryngoscope. 2002;112(2):393–6.
Gjedde, A. Glucose metabolism. In Principles and Practice of Nuclear Medicine, 2nd ed.; Wagner. H., Ed.; W B Saunders: Philadelphia, PA, USA, 1995:56–62.
Simpson IA, Dwyer D, Malide D, Moley KH, Travis A, Vannucci SJ. The facilitative glucose transporter GLUT3: 20 years of distinction. Am J Physiol Endocrinol Metab. 2008;295(2):E242–53.
Masoud GN, Li W. HIF-1α pathway: role, regulation and intervention for cancer therapy. Acta Pharm Sin B. 2015;5(5):378–89.
Lech G, Słotwiński R, Słodkowski M, KRASnodębski IW. Colorectal cancer tumour markers and biomarkers: recent therapeutic advances. World J Gastroenterol. 2016;22(5):1745–55.
Aghagolzadeh P, Radpour R. New trends in molecular and cellular biomarker discovery for colorectal cancer. World J Gastroenterol. 2016;22(25):5678–93.
Ahn TS, Jeong D, Son MW, et al. The BRAF mutation is associated with the prognosis in colorectal Cancer. J Cancer Res Clin Oncol. 2014;140(11):1863–71.
Sanz-Garcia E, Argiles G, Elez E, Tabernero J. BRAF mutant colorectal cancer: prognosis, treatment, and new perspectives. Ann Oncol. 2017;28(11):2648–57.
Phipps AI, Buchanan DD, Makar KW, et al. KRAS-mutation status in relation to colorectal cancer survival: the joint impact of correlated tumour markers. Br J Cancer. 2013;108(8):1757–64.
Lee JY, Ahn J, Park WS, et al. Colorectal cancer prognosis is not associated with BRAF and KRAS mutation-a STROBE compliant study. J Clin Med. 2019;8(1):111.
Goto T, Shinmura K, Yokomizo K, et al. Expression levels of thymidylate synthase, dihydropyrimidine dehydrogenase, and thymidine phophorylase in patients with colorectal cancer. Anticancer Res. 2012;32:1757–62.
Morikawa T, Kuchiba A, Liao X, et al. Tumor TP53 expression status, body mass index, and prognosis in colorectal cancer. Int J Cancer. 2012;131(5):1169–78.
da Costa AC, Stevanato Filho PR, Aguiar Júnior S, et al. Prognostic value of factors associated with hypoxia and lipid metabolism in patients with colorectal cancer. Applied Cancer Research. 2017;37:44.
Jóźwiak P, Krześlak A, Pomorski L, Lipińska A. Expression of hypoxia-related glucose transporters GLUT1 and GLUT3 in benign, malignant and non-neoplastic thyroid lesions. Mol Med Rep. 2012;6:601–6.
This research was supported by a grant from the National Research Foundation of Korea (NRF) funded by the Korean government (Ministry of Science, ICT & Future Planning, NRF-2016R1E1A1A01942072, http://www.nrf.re.kr/index) to YJC and the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (grant number: HI17C0048) to KK. The funding bodies had no involvement in the design of the study, collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
Ethics approval and consent to participate
Consent for publication
All authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Univariate cox regression analysis of glycolysis related genes for OS and DFS Table S2. Univariate and multivariate Cox regression analyses of BRAF mutation and the clinicopathological factors in the discovery set (TCGA COADREAD cohort) Table S3. Univariate and multivariate Cox regression analyses of KRAS mutation and the clinicopathological factors in the discovery set (TCGA COADREAD cohort) Table S4. Univariate Cox regression analysis of HIF, TP53, TYMS genes for OS and DFS in the discovery set (TCGA COADREAD cohort) (DOCX 25 kb)
Figure S1. Kaplan-Meier Survival analysis for overall survival and disease-free survival in colorectal cancer patients according to mutational status in the Discovery set. Overall survival (a) and disease-free survival (b) according to BRAF mutation status. Overall survival (c) and disease-free survival (d) according to KRAS mutation status. Figure S2. R code used for statistical analysis. The datasets analyzed during the study are available in the following repositories: TCGA COADREAD : cBioportal for cancer genomics (http://www.cbioportal.org.). GSE39582: Gene Expression Omnibus (GEO)(https://www.ncbi.nlm.nih.gov/geo/.) (PDF 122 kb)
About this article
Cite this article
Kim, E., Jung, S., Park, W.S. et al. Upregulation of SLC2A3 gene and prognosis in colorectal carcinoma: analysis of TCGA data. BMC Cancer 19, 302 (2019). https://doi.org/10.1186/s12885-019-5475-x