- Open Access
Upregulation of the solute carrier family 7 genes is indicative of poor prognosis in papillary thyroid carcinoma
World Journal of Surgical Oncologyvolume 16, Article number: 235 (2018)
The solute carrier (SLC) 7 family genes comprise 14 members and function as cationic amino acid/glycoprotein transporters in many cells, they are essential for the maintenance of amino acid nutrition and survival of tumor cells. This study was conducted to analyze the associations of SLC7 family gene expression with mortality in papillary thyroid carcinoma (PTC).
Clinical features, somatic mutations, and SLC7 family gene expression data were downloaded from The Cancer Genome Atlas database. Linear regression model analysis was performed to analyze the correlations between SLC7 family gene expression and clinicopathologic features. Kaplan-Meier survival and logistic regression analyses were performed to characterize the associations between gene expression and patients’ overall survival.
Patient mortality was negatively associated with age and tumor size but positively increased cancer stage and absence of thyroiditis in PTC patients. Kaplan-Meier survival analysis indicated that patients with high SLC7A3, SLC7A5, and SLC7A11 expression levels exhibited poorer survival than those with low SLC7A3, SLC7A5, and SLC7A11 expression levels (P < 0.05 for all cases). Logistic regression analysis showed that SLC7A3, SLC7A5, and SLC7A11 were associated with increased mortality (odds ratio [OR] 8.61, 95% confidence interval [CI] 2.3–55.91; OR 3.87, 95% CI 1.18–17.31; and OR 3.87, 95% CI 1.18–17.31, respectively.
Upregulation of SLC7A3, SLC7A5, and SLC7A11 expression was associated with poor prognosis in PTC patients, and SLC7 gene expression levels are potentially useful prognostic biomarkers.
PTC has become the fifth most common cancer in women in the USA, and the incidence rate continues to increase worldwide . PTC consists of four main histologic subtypes, including papillary thyroid carcinoma (PTC), follicular, anaplastic, and medullary thyroid cancer. PTC patients typically exhibit a favorable prognosis compared with follicular thyroid carcinoma patients . However, PTC occasionally dedifferentiates and becomes aggressive and lethal. Therefore, the identification of key prognostic biomarkers and effective druggable targets is critical to the survival of PTC patients.
The solute carrier (SLC) 7 family genes comprise 14 genes. The gene family can be divided into two subgroups, the cationic amino acid transporters (SLC7A1-4) and glycoprotein-associated transporters (SLC7 family 5-14) . Tumor cells are highly demanding for nutrients to support their fast growth rate. Therefore, amino acid transporters are essential for the maintenance of amino acid nutrition and survival of tumor cells . To date, two amino acid transporters have been found to be over-expressed in cancers: SLC7A5 and SLC7A11. SLC7A5 is highly expressed in most cancers . The hypoxia-inducible factor HIF2a upregulates SLC7A5 , which plays a critical role in cancer growth and progression. SLC7A5 promoter has canonical binding sites for the oncogene c-Myc, and overexpression of the oncogene causes upregulation of SLC7A5 . Increased expression of SLC7A11 is observed in a variety of cancers where the transporter-assisted promotion of glutathione synthesis reduces oxidative damage and protects the cancer cells from apoptosis .
The associations among mortality, clinicopathologic characteristics, and SLC7 family genes remain unknown in PTC. To address these issues, the present study investigated the relationship between mortality and clinicopathologic characteristics, mortality and SLC7 family gene expression, clinicopathologic characteristics and SLC7 family gene expression by analyzing a large set of PTC patient data from The Cancer Genome Atlas (TCGA) database .
Materials and methods
Data used in the study included normalized SLC7 family gene mRNA expression data and various clinical data from 355 PTC patients. Data were downloaded from the TCGA database  (http://firebrowse.org/?cohort=THCA&download_dialog=true). The clinicopathologic characteristics analyzed in the study included patients’ age, gender, thyroiditis, tumor size, extrathyroidal extension (ETE), American Joint Committee on Cancer (AJCC) stage, survival status, and follow-up data (months).
BRAFV600E mutational status is associated with recurrence and worse prognosis in PTC patients [10,11,12]. RAS mutations may confer a more aggressive phenotype in some PTC cases, increasing a patient’s risk for tumor recurrence, distant metastases, and death [13, 14]. Therefore, to analyze the associations between SLC7 family gene expression and BRAF and RAS mutation status, we downloaded somatic mutations of BRAF and RAS genes from the TCGA database. PTC patients were divided into BRAF- or RAS-mutant and wild-type groups based on their mutation status.
Correlations between clinical characteristics, mortality, and expression of SLC7 family genes
In order to characterize the associations between clinical features and mortality, Student’s t test was applied to compare the age and tumor size between PTC patients who were deceased versus alive. Fisher’s exact test was used to analyze the associations between survival status and patients’ gender, thyroiditis, ETE, AJCC cancer stage, BRAF, and RAS mutation status. Linear regression models were used to evaluate the associations between clinical features and SLC7 family gene expression. For correlation coefficients, t values were extracted. Pearson’s correlation was conducted between different SLC7A mRNA expressions in PTC. All statistical analyses were conducted in R (version 3.2.0), and P < 0.05 was considered statistically significant.
To characterize the association of SLC7 family gene expression with patient survival, Patients were assigned to the “high-expression” group if they exhibited SLC7 family gene expression levels greater than the median values, whereas those patients exhibiting expression levels less than the median values were assigned to the “low-expression” group. Kaplan-Meier survival analysis was performed to plot survival curves, and the log-rank test was utilized to compare the difference in survival rates between the high- and low-expression groups using the R package of survival [15, 16]. Univariate and multivariate survival analyses were performed using logistic regression model. P < 0.05 was considered statistically significant.
General characteristics of 355 PTC patients
The age of 355 PTC patients range from 15 to 88 years old, with the mean of 46.48. The average tumor size was 3 cm. Two hundred fifty-eight patients are female and 97 patients are male. One hundred sixteen and 118 PTC patients had thyroiditis and extrathyroidal extension respectively. One hundred ninety-five cancer samples had BRAF mutations, while 160 samples were wild-type. Twenty-four and 331 cancer samples were RAS-mutant or wild-type respectively. At the last day of follow-up, 14 PTC patients were deceased, and 341 were alive. The mean follow-up was 34.94 months (interquartile range, 14.31–43.87 months).
The association between mortality and clinicopathologic characteristics in PTC
Among the clinicopathologic characteristics, patients with older age or larger tumor sizes exhibited significantly worse mortality than those with younger age or smaller tumor sizes, respectively (P < 0.0001 for all cases, Student’s t test, Fig. 1a, b). Moreover, mortality was significantly associated with increased cancer stage and absence of thyroiditis in PTC patients (P < 0.05 for all cases, Fisher’s exact test, Table 1). However, no significant association was observed between patient mortality and gender, ETE BRAF, or RAS mutation status (P > 0.05 for all cases, Fisher’s exact test, Table 1).
The association between clinicopathologic characteristics and SLC7 family gene expression in PTC
A linear regression model was applied to characterize the associations between clinicopathologic characteristics and mRNA expression values of SLC7 family genes. Tumor size was negatively correlated with SLC7A3 mRNA expression. Thyroiditis was significantly correlated with SLC7A2, SLC7A3, SLC7A7, and SLC7A9 (Table 2). ETE was significantly associated with SLC7A3, SLC7A8, and SLC7A9. AJCC TNM stage was associated with expression of SLC7A3, SLC7A8, and SLC7A11 (Table 2). RAS mutation status was associated with SLC7A3, SLC7A4, SLC7A6, and SLC7A8. BRAF mutation status was associated with SLC7A1, SLC7A2, SLC7A4, SLC7A6, SLC7A8, SLC7A9, and SLC7A14. Mortality was positively associated with greater mRNA expression of SLC7A3 (t = 2.39, P = 0.02) and SLC7A5 (t = 2.02, P = 0.04) expression (Tables 2 and 3). The results demonstrated that SLC7A3 and SLC7A5 expression was significantly correlated with patient mortality in PTC.
The correlations of expression of SLC7 family genes
In order to characterize the correlations between different SLC7A mRNA expressions in PTC, the correlations of expression of SLC7 family genes were analyzed. The expression of SLC7 family genes was highly correlated, with SLC7A1 and SLC7A9 expression most frequently correlated with other SLC7 family members (P < 0.05 for all cases, Pearson’s correlation, Fig. 2). The expression of SLC7A3 was significantly correlated with the expression of SLC7A5, SLC7A7, SLC7A8, and SLC7A9. SLC7A5 expression showed significant correlation with SLC7A1, SLC7A2, SLC7A3, SLC7A4, SLC7A7, SLC7A11, and SLC7A14 genes expression. The expression of SLC7A11 was correlated with SLC7A1, SLC7A2, SLC7A5, SLC7A7, SLC7A9, and SLC7A13 with significant evidence (P < 0.05 for all cases, Pearson’s correlation, Fig. 2).
Survival analyses between patient mortality and SLC7 family gene expression in PTC
To evaluate the association of SLC7 family gene expression with patient survival, the 355 PTC patients were divided into low- and high-expression groups based on the median values. Kaplan-Meier survival analysis indicated that patients with high SLC7A3, SLC7A5, and SLC7A11 expression levels exhibited poorer survival than those with low SLC7A3, SLC7A5, and SLC7A11 expression levels (P < 0.05 for all cases, log rank test) (Fig. 3a–c). Univariate analysis using logistic regression model showed that SLC7A3, SLC7A5, and SLC7A11 were associated with increased mortality (P = 0.005, OR 8.61, 95% CI 2.3–55.91; P = 0.04, OR 3.87, 95% CI 1.18–17.31; and P = 0.04, OR 3.87, 95% CI 1.18–17.31, respectively.
To further validate the association of patients’ survival with SLC7A3, SLC7A5, and SLC7A11 expression levels, multivariate analyses were applied between patients’ survival and the mortality-associated features, including patients’ age, tumor size, cancer stage, thyroiditis, and SLC7A3, SLC7A5, and SLC7A11 expression levels. Multivariate survival analyses confirmed that SLC7A3, SLC7A5, and SLC7A11 were associated with increased mortality (P = 0.03, OR 12.13, 95% CI 1.65–256.97; P = 0.09, OR 4.56, 95% CI 0.88–34.73; and P = 0.03, OR 7.36, 95% CI 1.4–60.47, respectively).
In the present study, we have investigated the associations between mortality and clinicopathologic characteristics, mortality and SLC7 family gene expression, clinicopathologic characteristics, and SLC7 family gene expression in PTC. Many of SLC7 family members are significantly correlated to clinical traits of PTC patients. In particular, we found three genes, SLC7A3, SLC7A5, and SLC7A11, that showed significant associations with mortality, which might have clinical values. The gene SLC7A3 is a sodium-independent cationic amino acid transporter, which is glycosylated and localized to the plasma membrane. SLC7A3 is selective for cationic l-amino acids transportation . So far, there is no report regarding the involvement of SLC7A3 in cancers, for the first time, our study found that SLC7A3 expression was associated with ETE, higher cancer stage, BRAF, RAS mutation, and mortality in PTC. Though the molecular mechanism on how SLC7A3 expression impacts mortality is unknown, SLC7A3 expression was independently associated with the RAS mutation, which is a well-known negative prognostic factor in PTC. This might in part explain the reasons; however, it still needs further studies.
Among the 14 SLC7 family genes, SLC7A5 has been most widely investigated in various cancers. SLC7A5 functions as an l-type amino-acid transporter that transports large neutral amino acids  and is over-expressed in many cancer types, such as prostatic, esophageal, gastric, and pancreatic carcinomas [19,20,21,22]. SLC7A5 is involved in the growth and mortality of a variety of tumor cells [18,19,20,21,22,23]. SLC7A5 expression is indicative of a poor prognosis in pancreas cancer [22, 24], melanoma , bile duct adenocarcinomas , and clear cell renal cell carcinoma . Consistent with previously published results, our study demonstrated that high SLC7A5 expression was related to patient mortality in PTC.
SLC7A11 is a component of the cysteine/glutamate transporter, which plays a key role in glutathione synthesis. It has been evidenced that SLC7A11 has oncogenic functions in glioma. Overexpression of SLC7A11 resulted in decreased endogenous ROS levels as well as decreased migration and invasion in glioblastoma . SLC7A11 expression is associated with accelerated growth and tumor-associated seizures  and predicts poor survival in patients with malignant glioma [29, 30]. In addition, over-expression of SLC7A11 increased resistance to oxidative stress and decreased chemosensitivity to temozolomide in the context of glioma . For the first time, our study identified that SLC7A11 expression was associated with a poor prognosis in PTC, expanding the functions of SLC7A11 in tumorigeneses.
SLC7 family gene expression analysis might have clinical value in the near future. Cytologic or surgical specimens of PTC patients exhibiting high SLC7A3, SLC7A5, and SLC7A11 expression are expected to be associated with a poor prognosis. Therefore, more aggressive treatment or frequent follow-up may be recommended for these patients. Additionally, SLC7A3, SLC7A5, and SLC7A11 may represent novel and potential druggable targets in PTC. Drugs inhibiting SLC7A11 gene could have great therapeutic potential, as demonstrated by in vitro and vivo studies in breast cancer  and resistant head and neck cancer cells  using the SLC7A11 inhibitors.
Though this study provides the evidence that upregulation of SLC7A3, SLC7A5, and SLC7A11 genes is an indicator of poor overall survival in PTC patients, the study has certain limitation. The overexpression of SLC7A5 and SLC7A11 may be partly attributed to macrophages [34, 35], and tumor-associated macrophages are activated and present in thyroid cancer microenvironment . Thus, not the mRNA expression of SLC7A5 and SLC7A11 but the presence of TAM in PTC tissue is the prognostic marker, this possibility needs further studies.
In conclusion, upregulation of SLC7A3, SLC7A5, and SLC7A11 genes is associated with poor overall survival in PTC patients. SLC7 family gene expression represents a potentially prognostic biomarker to predict survival in PTC.
- AJCC :
American Joint Committee on Cancer
Papillary thyroid carcinoma
The Cancer Genome Atlas
Siegel R, Miller K, Jemal A. Cancer statistics , 2015. CA Cancer J Clin. 2015;65(1):29.
Oyer SL, Fritsch VA, Lentsch EJ. Comparison of survival rates between papillary and follicular thyroid carcinomas among 36,725 patients. Ann Otol Rhinol Laryngol. 2014;123(2):94–100.
Closs EI, Boissel JP, Habermeier A, Rotmann A. Structure and function of cationic amino acid transporters (CATs). J Membr Biol. 2006;213(2):67–77.
Bhutia YD, Babu E, Ramachandran S, Ganapathy V. Amino acid transporters in cancer and their relevance to “glutamine addiction”: novel targets for the design of a new class of anticancer drugs. Cancer Res. 2015;75(9):1782–8.
Kaira K, Oriuchi N, Imai H, Shimizu K, Yanagitani N, Sunaga N, et al. L-type amino acid transporter 1 and CD98 expression in primary and metastatic sites of human neoplasms. Cancer Sci. 2008;99(12):2380–6.
Elorza A, Soro-Arnáiz I, Meléndez-Rodríguez F, Rodríguez-Vaello V, Marsboom G, de Cárcer G, et al. HIF2α acts as an mTORC1 activator through the amino acid carrier SLC7A5. Mol Cell. 2012;48(5):681–91.
Gao P, Tchernyshyov I, Chang TC, Lee YS, Kita K, Ochi T, et al. C-Myc suppression of miR-23a/b enhances mitochondrial glutaminase expression and glutamine metabolism. Nature. 2009;458(7239):762–5.
Lewerenz J, Hewett SJ, Huang Y, Lambros M, Gout PW, Kalivas PW, et al. The cystine/glutamate antiporter system x c − in health and disease: from molecular mechanisms to novel therapeutic opportunities. Antioxid Redox Signal. 2013;18(5):522–55.
Network TCGAR, Agrawal N, Akbani R, Aksoy BA, Ally A, Arachchi H, et al. Integrated genomic characterization of papillary thyroid carcinoma. Cell. 2014;159(3):676–90.
Elisei R, Viola D, Torregrossa L, Giannini R, Romei C, Ugolini C, et al. The BRAFV600E mutation is an independent, poor prognostic factor for the outcome of patients with low-risk intrathyroid papillary thyroid carcinoma: single-institution results from a large cohort study. J Clin Endocrinol Metab. 2012;97(12):4390–8.
Prescott JD, Sadow PM, Hodin RA, Le LP, Gaz RD, Randolph GW, et al. BRAF(V600E) status adds incremental value to current risk classification systems in predicting papillary thyroid carcinoma recurrence. Surgery. 2012;152(6):984–90.
Virk RK, Van Dyke AL, Finkelstein A, Prasad A, Gibson J, Hui P, et al. BRAFV600E mutation in papillary thyroid microcarcinoma: a genotype–phenotype correlation. Mod Pathol. 2013;26(1):62–70.
Volante M, Rapa I, Gandhi M, Bussolati G, Giachino D, Papotti M, et al. RAS mutations are the predominant molecular alteration in poorly differentiated thyroid carcinomas and bear prognostic impact. J Clin Endocrinol Metab. 2009;94(12):4735–41.
Fukahori M, Yoshida A, Hayashi H, Yoshihara M, Matsukuma S, Sakuma Y, et al. The associations between RAS mutations and clinical characteristics in follicular thyroid tumors: new insights from a single center and a large patient cohort. Thyroid. 2012;22(7):683–9.
Therneau TM, Grambsch PM. Modeling Survival Data: Extending the Cox Model. New York: Springer; 2000. ISBN 0-387-98784-3.
Fox J. Cox proportional-hazards regression for survival data the cox proportional-hazards model. Most. 2002;2008(June):1–18.
Vékony N, Wolf S, Boissel JP, Gnauert K, Closs EI. Human cationic amino acid transporter hCAT-3 is preferentially expressed in peripheral tissues. Biochemistry. 2001;40(41):12387–94.
Yanagida O, Kanai Y, Chairoungdua A, Kim DK, Segawa H, Nii T, et al. Human L-type amino acid transporter 1 (LAT1): characterization of function and expression in tumor cell lines. Biochim Biophys Acta - Biomembr. 2001;1514(2):291–302.
Sakata T, Ferdous G, Tsuruta T, Satoh T, Baba S, Muto T, et al. L-type amino-acid transporter 1 as a novel biomarker for high-grade malignancy in prostate cancer. Pathol Int. 2009;59(1):7–18.
Kobayashi H, Ishii Y, Takayama T. Expression of L-type amino acid transporter 1 (LAT1) in esophageal carcinoma. J Surg Oncol. 2005;90(4):233–8.
Ichinoe M, Mikami T, Yoshida T, Igawa I, Tsuruta T, Nakada N, et al. High expression of L-type amino-acid transporter 1 (LAT1) in gastric carcinomas: comparison with non-cancerous lesions. Pathol Int. 2011;61(5):281–9.
Yanagisawa N, Ichinoe M, Mikami T, Nakada N, Hana K, Koizumi W, et al. High expression of L-type amino acid transporter 1 (LAT1) predicts poor prognosis in pancreatic ductal adenocarcinomas. J Clin Pathol. 2012;65(11):1019–23.
Kobayashi T, Yan H, Kurahashi Y, Ito Y, Maeda H, Tada T. Role of GalNAc4S-6ST in Astrocytic Tumor Progression. PLoS One. 2013;8(1):1–9.
Kaira K, Sunose Y, Arakawa K, Ogawa T, Sunaga N, Shimizu K, et al. Prognostic significance of L-type amino-acid transporter 1 expression in surgically resected pancreatic cancer. Br J Cancer. 2012;107(4):632–8.
Shimizu A, Kaira K, Kato M, Yasuda M, Takahashi A, Tominaga H, et al. Prognostic significance of L-type amino acid transporter 1 (LAT1) expression in cutaneous melanoma. Melanoma Res. 2015;25(5).
Yanagisawa N, Hana K, Nakada N, Ichinoe M, Koizumi W, Endou H, et al. High expression of L-type amino acid transporter 1 as a prognostic marker in bile duct adenocarcinomas. Cancer Med. 2014;3(5):1246–55.
Betsunoh H, Fukuda T, Anzai N, Nishihara D, Mizuno T, Yuki H, et al. Increased expression of system large amino acid transporter (LAT)-1 mRNA is associated with invasive potential and unfavorable prognosis of human clear cell renal cell carcinoma. BMC Cancer. 2013;13.
Polewski MD, Reveron-Thornton RF, Cherryholmes GA, Marinov GK, Aboody KS. SLC7A11 overexpression in glioblastoma is associated with increased cancer stem cell-like properties. Stem Cells Dev. 2017;26(17):1236–46.
Robert SM, Buckingham SC, Campbell SL, Robel S, Holt KT, Ogunrinu-Babarinde T, et al. SLC7A11 expression is associated with seizures and predicts poor survival in patients with malignant glioma. Sci Transl Med. 2015;7(289):289ra86.
Takeuchi S, Wada K, Toyooka T, Shinomiya N, Shimazaki H, Nakanishi K, et al. Increased xCT expression correlates with tumor invasion and outcome in patients with glioblastomas. Neurosurgery. 2013;72(1):33–41.
Polewski MD, Reveron-Thornton RF, Cherryholmes GA, Marinov GK, Cassady K, Aboody KS. Increased expression of system xc- in glioblastoma confers an altered metabolic state and temozolomide resistance. Mol Cancer Res. 2016;14(12):1229–42.
Timmerman LA, Holton T, Yuneva M, Louie RJ, Padró M, Daemen A, et al. Glutamine sensitivity analysis identifies the xCT antiporter as a common triple-negative breast tumor therapeutic target. Cancer Cell. 2013;24(4):450–65.
Roh JL, Kim EH, Jang H, Shin D. Aspirin plus sorafenib potentiates cisplatin cytotoxicity in resistant head and neck cancer cells through xCT inhibition. Free Radic Biol Med. 2017;104:1–9.
Yoon BR, Oh YJ, Kang SW, Lee EB, Lee WW. Role of SLC7A5 in metabolic reprogramming of human monocyte/macrophage immune responses. Front Immunol. 2018;9:53.
Kandasamy P, Gyimesi G, Kanai Y, Hediger MA. Amino acid transporters revisited: new views in health and disease. Trends Biochem Sci. 2018;43(10):752–89.
Naoum GE, Morkos M, Kim B, Arafat W. Novel targeted therapies and immunotherapy for advanced thyroid cancers. Mol Cancer Mol Cancer. 2018;17(1):1–15.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.