The amino acid transporter SLC7A5 confers a poor prognosis in the highly proliferative breast cancer subtypes and is a key therapeutic target in luminal B tumours
Breast Cancer Research volume 20, Article number: 21 (2018)
Breast cancer (BC) is a heterogeneous disease characterised by variant biology and patient outcome. The amino acid transporter, SLC7A5, plays a role in BC although its impact on patient outcome in different BC subtypes remains to be validated. This study aimed to determine whether the clinicopathological and prognostic value of SLC7A5 is different within the molecular classes of BC.
SLC7A5 was assessed at the genomic level, using Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) data (n = 1980), and proteomic level, using immunohistochemical analysis and tissue microarray (TMA) (n = 2664; 1110 training and 1554 validation sets) in well-characterised primary BC cohorts. SLC7A5 expression correlated with clinicopathological and biological parameters, molecular subtypes and patient outcome.
SLC7A5 mRNA and protein expression were strongly correlated with larger tumour size and higher grade. High expression was observed in triple negative (TN), human epidermal growth factor receptor 2 (HER2)+, and luminal B subtypes. SLC7A5 mRNA and protein expression was significantly associated with the expression of the key regulator of tumour cell metabolism, c-MYC, specifically in luminal B tumours only (p = 0.001). High expression of SLC7A5 mRNA and protein was associated with poor patient outcome (p < 0.001) but only in the highly proliferative oestrogen receptor (ER)+/ luminal B (p = 0.007) and HER2+ classes of BC (p = 0.03). In multivariate analysis, SLC7A5 protein was an independent risk factor for shorter breast-cancer-specific survival only in ER+ high-proliferation tumours (p = 0.02).
SLC7A5 appears to play a role in the aggressive highly proliferative ER+ subtype driven by MYC and could act as a potential therapeutic target. Functional assessment is necessary to reveal the specific role played by this transporter in the ER+ highly proliferative subclass and HER2+ subclass of BC.
Altered metabolic pathways have been readily accepted as part of the revised hallmarks of cancer where cancer cells are able to regulate their metabolism to provide energy and cellular building blocks required for their unremitting proliferation . Many cancer cells are highly reliant on amino acids for their growth, not only because they are a nitrogen donor for the synthesis of nucleotides and other amino acids, but also because they activate mammalian target of rapamycin complex1 (mTORC1) through nutrient signalling pathways which in turn regulates protein translation and cell growth [2, 3]. There is also increasing evidence that oncogenes and/or tumour-suppressor genes can reprogramme tumour cell metabolism including the direct regulation of the amino acid transporter, solute carrier family 7 member 5 (SLC7A5), by the oncogene MYC [4, 5] and the regulation of expression of the glutamine transporter, SLC1A5, by the tumour suppressor retinoblastoma (Rb) .
SLC7A5 is a sodium-independent transporter and acts as an amino acid exchanger by transporting large neutral amino acids such as leucine, phenylalanine and tryptophan by exchange with intracellular glutamine . It therefore functions in supplying amino acids to cancer cells and maintaining intra-cellular leucine, which is considered a master regulator of the mTORC1 signalling pathway [8,9,10]. For functional expression on the plasma membrane, SLC7A5 must heterodimerise with the heavy chain of SLC3A2 [7, 11].
It has been reported that SLC7A5 is highly expressed in a variety of cancers including oesophageal carcinoma , oral cancer  and lung adenocarcinoma . SLC7A5 is co-expressed with the glutamine transporter, SLC1A5, in many cancers suggesting a functional coupling of these transporters in supporting tumour progression . In this study, we aimed to assess SLC7A5 gene copy number and mRNA expression, alongside protein expression in large and well-characterised annotated cohorts of BC to determine its biological, clinicopathological and prognostic value in the different BC molecular classes with particular interest in the highly proliferative aggressive subgroups.
SLC7A5 copy number and gene expression
A cohort of 1980 BC tumours in the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC)  was used to evaluate SLC7A5 gene copy number aberrations (CNA) and gene expression. DNA/RNA was isolated from fresh frozen samples and genomic and transcriptional profiling was obtained using the Affymetrix SNP 6.0 and Illumina HT-12v3 platforms respectively. CNA were considered at the gene level by segments and the Šidák correction  was applied for multiple testing. Gene expression data were pre-processed and normalised as described previously . In this cohort, patients included were oestrogen receptor (ER)-positive (ER+) and/or lymph node (LN)-negative (LN-) and did not receive adjuvant chemotherapy, whereas ER- and LN+ patients received adjuvant treatment. X-tile (version 3.6.1, Yale University, USA) was applied to dichotomise SLC7A5 mRNA expression, based on prediction of breast-cancer-specific survival (BCSS). The association between the SLC7A5 mRNA expression and clinicopathological parameters, molecular subtypes, and patient outcome was investigated. The online dataset, Breast Cancer Gene Expression Miner v4.0 (http://bcgenex.centregauducheau.fr) and breast cancer data from The Cancer Genome Atlas (TCGA)  were used for external validation of SLC7A5 copy number/or mRNA expression.
Patients and tumours
This study evaluated well-characterised cohorts of patients with early-stage primary operable invasive BC, who presented aged ≤70 years. Patients in the discovery set (n = 1110) presented at Nottingham City Hospital between 1989 and 1998, while the validation set (n = 1554) includes patients who were presented between 1998 and 2006. Patient management was uniform and based on tumour characteristics by Nottingham Prognostic Index (NPI) and hormone receptor status. Patients within the NPI excellent prognostic group (score ≤3.4) received no adjuvant therapy, but those patients with NPI >3.4 received tamoxifen if ER-positive (± goserelin (Zoladex) in case the patients were premenopausal). Conversely, classical cyclophosphamide, methotrexate and 5-flurouracil (CMF) were used if the patients were ER-negative and fit enough to receive chemotherapy. None of the patients in this study received neoadjuvant therapy. Clinical history, tumour characteristics and information on therapy and outcomes are prospectively maintained. Outcome data included development and time to distant metastasis (DM) and breast-cancer-specific survival (BCSS). There was no difference in the distribution of clinicopathological parameters between the discovery and validation cohorts or the METABRIC series of patients (all correlation coefficients ≥0.91, all p < 0.0001) (Additional file 1: Table S1).
The antibody specificity of anti-SLC7A5 (EPR17573, Abcam, UK) was validated using western blotting in human embryonic kidney (HEK) 293 T over expression lysate (Origene Technologies, Rockville, MD, USA) and MDA-MB-175 (luminal B-like), T47D and MCF7 (luminal A)  breast cancer lysate (American Type Culture Collection; Rockville, MD, USA). A dilution of 1:200 of the primary antibody and 1:1000 HRP-conjugated (Dako) secondary antibodies were applied: 5% milk /PBS-Tween (0.1%) (Marvel Original Dried Skimmed Milk, Premier Food Groups Ltd., UK) was used for blocking. Mouse monoclonal anti-β-actin primary antibody was used as a marker of endogenously expressed control. SLC7A5 bands were visualised using enhanced chemiluminescence (ECL) showing a single specific band at the correct predicted size (40 kDa) for the SLC7A5 protein.
Tissue arrays and immunohistochemical analysis
The discovery set (n = 1110) were arrayed as previously described using a single 0.6-mm core sampled from the periphery of each invasive tumour . The validation set (n = 1554) were similarly arrayed using a tissue microarray (TMA) GrandMaster (3D Histech). Immunohistochemical (IHC) staining was performed on 4-μm TMA sections using the Novolink polymer detection system (Leica Biosystems, RE7150-K). Briefly, tissue slides were deparaffinised with xylene and rehydrated through three changes of alcohol. Heat-induced antigen epitope retrieval was performed in citrate buffer (pH 6.0) for 20 min using a microwave oven (Whirpool JT359 Jet Chef 1000 W). Endogenous peroxidase activity was blocked by peroxidase block for 5 min. Slides were washed with Tris-buffered saline (TBS, pH 7.6), followed by application of protein block for 5 min. Following another TBS wash, sections were incubated with the primary SLC7A5 antibody diluted at 1:50 in Leica antibody diluent (RE AR9352, Lieca, Biosysytems, UK) overnight at 4 °C. Slides were washed with TBS followed by incubation with post primary block for 30 min followed by a TBS wash. Novolink polymer was applied for 30 min: 3,3′-diaminobenzidine (DAB) chromogen was applied for 5 min. Slides were counterstained with Novolink haematoxylin for 6 min, dehydrated and coverslipped.
Stained TMA sections were scored using high resolution digital images (NanoZoomer; Hamamatsu Photonics, Welwyn Garden City, UK), at × 20 magnification. Evaluation of staining for SLC7A5 was based on a semi-quantitative assessment of digital images of the cores using a modified histochemical score (H-score) which includes an assessment of both the intensity and the percentage of stained cells . Staining intensity was assessed as follows: 0, negative; 1, weak; 2, medium; 3, strong, and the percentage of the positively stained tumour cells was estimated subjectively. The final H-score was calculated by multiplying the percentage of positive cells (0–100) by the intensity (0–3), producing a total range of 0–300. Dichotomisation of protein expression in predicting BCSS was determined using x-tile software.
Immunhistochemical staining and dichotomisation of the other biomarkers included in this study were as per previous publications [20, 22,23,24,25,26,27,28,29,30]. ER and progesterone receptor (PgR) positivity was defined as ≥ 1% staining. Immunoreactivity of HER2 in TMA cores was scored using standard HercepTest guidelines (Dako). Chromogenic in situ hybridisation (CISH) was used to quantify HER2 gene amplification in borderline cases using the HER2 FISH pharmDx™ plus HER2 CISH pharmDx™ kit (Dako) and was assessed according to the American Society of Clinical Oncology guidelines. BC molecular subtypes were defined based on tumour IHC profile and the Elston-Ellis  mitotic score as: ER+/HER2- low proliferation (mitotic score 1), ER+/HER2- high proliferation (mitotic score 2 and 3); HER2-positive class: HER2+ regardless of ER status; triple negative: ER-, PgR- and HER2- . Basal-like phenotype was defined as tumours expressing cytokeratin (Ck) 5/6, and/or Ck14 and/or Ck17.
Statistical analysis was performed using SPSS 22.0 statistical software (SPSS Inc., Chicago, IL, USA). Spearman’s correlation coefficient was calcualted to examine the association between continuous variables. The chi-square test was performed to analyse relationships between categorical variables. For the continuous variables, differences between three or more groups were assessed using one-way analysis of variance (ANOVA) with the post-hoc Tukey multiple comparison test (for normally distributed data) or Kruskal-Wallis test (for non-normal distribution). Differences between two groups were assessed using the t test (normally distributed data) or Mann-Whitney test (non-normal distribution). Survival curves were analysed by the Kaplan-Meier and log rank test. Cox’s proportional hazard method was performed for multivariate analysis to identify the independent prognostic/predictive factors. The p values were adjusted for multiple testing using the Bonferroni correction. A p value ˂0.05 was considered significant. The study endpoints were 10-year BCSS or distant metastasis-free survival (DMFS).
SLC7A5 genomic profiling
SLC7A5 was amplified in 0.3% and 0.6% of BC cases in the METABRIC and TCGA datasets, respectively, while deletion (deep and shallow) was detected in 56% and 68% of cases in the same cohorts respectively. Point mutations in SLC7A5 were extremely rare, where TCGA data reported just one case with a missense mutation (Additional file 2: Figure S4A) [33, 34]. SLC7A5 is situated on chromosome 16 (16q24.2); all the annotated genes, which were located on 16q  were selected to determine their CNV in relation to SLC7A5 and assess whether these aberrations were locus-specific or involved large chromosomal segments. There was significant positive correlation between SLC7A5 deletion and the deletion of all genes (p < 0.001, Additional file 3: Table S2) in both the METABRIC and TCGA data. However, amplification of three genes (FANCA, CBFA2T3 and CDT1) showed significant association with the amplified SLC7A5 (p ≤ 0.03, data not shown) in the afore-mentioned datasets together.
SLC7A5 expression in breast cancer
SLC7A5 protein expression was observed, predominantly in the membrane of invasive breast cancer cells, with expression levels varying from absent to high (Fig. 1b and c). The distribution of the SLC7A5 protein expression was unimodal and left-skewed. The SLC7A5 mRNA expression had a normal distribution. Expression of SLC7A5 mRNA and protein were dichotomised using cut points derived from prediction of patient survival using X-tile (https://medicine.yale.edu/lab/rimm/research/software.aspx; Yale University). Positive SLC7A5 expression (>15 H-score) was observed in 191/1110 (17%) and 268/1554 (17%) of cases in the discovery and validation sets, respectively, while high SLC7A5 mRNA expression (log2 intensity >8) was observed in 1019/1923 (53%) of the METABRIC breast cancer cases. A total of 49/1980 (2.4%) of cases had a copy number (CN) gain of SLC7A5 and 530/1980 (26.7%) a CN loss. Significant association was observed between SLC7A5 copy number variation (CNV) and SLC7A5 mRNA expression (p < 0.001, Fig. 2).
SLC7A5 and clinicopathological parameters
Table 1 summarises the associations between SLC7A5 protein expression including larger tumour size, high tumour grade, and poor Nottingham Prognostic Index (NPI) (all p < 0.001). Regarding BC metastatic sites, high SLC7A5 protein levels were associated with the development of distant metastases (DM) in the brain (p < 0.001) and lung (p = 0.04), while there was no association with development of DM in the bone or liver.
High SLC7A5 mRNA expression was significantly associated with larger tumour size (Fig. 3a, p < 0.001), higher grade (Fig. 3b, p < 0.001), positive nodal metastasis (Fig. 3c, p< 0.001) and poor NPI (Fig. 3d, p < 0.001). Both SLC7A5 mRNA and SLC7A5 protein were associated with medullary-like tumours. Where data were available, these associations were confirmed using the Breast Cancer Gene-Expression Miner v4.0 (Additional file 4: Figure S1A, B) and the TCGA data (Additional file 2: Figure S4B). In addition SLC7A5 copy number loss was significantly associated with good prognostic parameters including, lower grade and good NPI (Table 2, p < 0.001). There was positive association between SLC7A5 copy number gain and MYC gain (p < 0.001, Table 2).
SLC7A5 expression in molecular BC subtypes
SLC7A5 protein expression was associated with negative hormone receptor status and HER2+ tumours (all p ≤ 0.002, Table 3) and it was highly expressed in triple negative (TN) and basal-like phenotype malignancies compared to non-TN and non-basal-like tumours (p < 0.001, Table 3). Similarly, high expression of SLC7A5 mRNA was significantly associated with hormone receptor negative (ER- and PgR-) and HER2+ tumours (all p < 0.001, Table 4). These results were in concordance with the findings of the Breast Cancer Gene-Expression Miner v4.0 (Additional file 4: Figure S1C-F) and TCGA data analysis (Additional file 2: Figure S4C-E).
When comparing the levels of SLC7A5 mRNA expression in the intrinsic (prediction analysis of microarray 50 (PAM50)) subtypes , high expression was observed in basal-like, HER2+ and lLuminal B tumours (Fig. 3f, p < 0.001). Similarly, within the METABRIC integrative clusters, high SLC7A5 mRNA expression was associated with clusters 5 (ERBB2 amplified), 9 (luminal B subgroup) and 10 (TN/basal-like) (p < 0.001, Fig. 3e). In the SCMGENE subtypes there was higher expression of SLC7A5 mRNA in the ER+/HER2- high proliferation class (luminal B) compared with the ER+/HER2- low proliferation class (luminal A) (p < 0.001, Fig. 3g). Association between SLC7A5 mRNA and PAM50 subtypes was confirmed using the Breast Cancer Gene-Expression Miner v4.0 (Additional file 4: Figure S1G). There was lower expression of SLC7A5 protein in the low-proliferation tumours than in the other defined molecular subtypes subtypes (p < 0.001, Table 1).
At the gene level, there was a greater copy number gain of SLC7A5 (p < 0.001, Table 2 in the basal-like subtype, while SCL7A5 copy number loss was mainly observed in the luminal A subtype (p < 0.001, Table 2).
SLC7A5 expression and other associated markers
Correlation between SLC7A5 mRNA and associated genes was investigated in the METABRIC dataset (Table 5). The genes were selected based on previous publications, and were either regulatory genes or others that share or support the SLC7A5 biological function, which focused mainly on glutamine transport and glutamine metabolism [2, 5, 15, 37,38,39,40,41]. There was positive correlation between SLC7A5 mRNA expression and the expression of regulatory genes, several amino acid transporters and genes involved in the glutamine-proline regulatory axis. There was a positive relationship between SLC7A5 and MYC, mTOR and ATF4 (p < 0.001) and the positive relationship between MYC, HIF2A and SLC7A5 was only observed in luminal B tumours (p = 0.01 and p < 0.001, respectively).
High SLC7A5 mRNA expression was specifically associated with the enzymes involved with conversion of glutamine (Gln) to proline, where PYCR1 and ALDH18A1 showed a positive relationship with SLC7A5 in almost all subtypes (p < 0.02).
The majority of glutamine transporters were significantly associated with SLC7A5 expression primarily in triple negative tumours and to a lesser extent luminal and HER2+ subtypes. SLC7A5 was significantly correlated with SLC1A5 in all subtypes (p < 0.02).
TP53 mutations were also highly prevalent in breast tumours where there was high SLC7A5 mRNA expression (p < 0.001, Tables 3 and 4). Moreover, high SLC7A5 protein was positively associated with high p53 protein (p < 0.001).
SLC7A5 protein expression was significantly expressed in breast tumours with high Ki67, and the upstream effector MYC (p < 0.001, Table 6). SLC1A5, GLS, PYCR1 and PIK3CA were significantly expressed in breast tumours with high expression of SLC7A5 (p < 0.001), while the low expression of SLC7A5 was associated with high levels of p-mTORC1 (p < 0.001) (Table 6).
SLC7A5 expression and patient outcome
Both high SLC7A5 mRNA (p < 0.001) (Fig. 4a) and protein (p < 0.001) expression were associated with poor BCSS in the discovery and validation sets (Fig. 5a, b). This association was anticipated as the cutoff was based on the prediction of BCSS.
While SLC7A5 mRNA expression was not predictive of BCSS in any specific molecular class (Fig. 4b-e), high expression of SLC7A5 protein was only predictive of shorter BCSS in ER+ high proliferation (p = 0.007, Fig. 5d) and HER2+ tumours (p = 0.03, Fig. 5f). There was no association between SLC7A5 protein and outcome in ER+ low proliferation (Fig. 5c) or in TNBC (Fig. 5e). In multivariate Cox regression analysis, SLC7A5 mRNA was a predictor of shorter BCSS independent of tumour size, grade or lymph node stage (p = 0.006, Additional file 5: Table S3) but not in any specific subtype. However, SLC7A5 protein was significant only in the ER+ high-proliferation tumours (p = 0.02, Table 7) and not in any other subtypes (data not shown).
Likewise, high SLC7A5 protein expression was associated with shorter distant metastases-free survival (DMFS) (p < 0.001; Additional file 6: Figure S2A, B) within the ER+ high-proliferation class (p = 0.03, Additional file 6: Figure S2D) but not in other subtypes (Additional file 6: Figure S2C, E, F). This association was identified in the discovery set and validated in the validation set. The relationship between high SLC7A5 mRNA expression and poor patient outcome in ER+ disease, but not ER- disease, was confirmed using Breast Cancer Gene-Expression Miner (Additional file 7: Figure S3A, B, C).
Breast cancer is a heterogeneous disease with various subtypes  differing in terms of morphology, molecular and biological profiles, response to therapy and clinical behaviour. In addition, different subtypes exhibit disparity in their metabolic pathways and their nutritional needs. The most common form of BC (~ 55–80%) is the ER+/luminal tumour [43, 44], and tumours that belong to this class are also variable in terms of recurrence, mortality rates and disease prognosis . Therefore, understanding the biology of BC and exploring the metabolic pathways could help to identify potential novel therapeutic targets.
Cancer cells must alter their metabolism in order to satisfy the demands of necessary energy and cellular building blocks. It is widely known that amino acid transport systems play a principal role in sustaining the proliferation of cancer cells by supplying the required amino acids for protein synthesis and by activation of nutrient signalling through the mTORC1 complex. This study has revealed for the first time that SLC7A5 is a key amino acid transporter in the more aggressive and highly proliferative ER+ tumours.
SLC7A5 is located in 16q24.2. This study showed that SLC7A5 deletion, but not amplification, was significantly correlated with all the annotated genes located in the same chromosomal region, indicating that the deletion was not locus-specific. Interestingly, E-cadherin (CDH1), which was located in 16q22.1, was also implicated. It is widely known that most lobular tumours harbour loss of heterozygosity (LOH) at chromosome 16 and are missing the wild type CDH1 allele . In this study, approximately 40% of METABRIC cases with SLC7A5 loss were associated with invasive lobular histology. In addition, SLC7A5 protein expression in lobular carcinoma has a relatively lower mean rank value compared to the other histological subtypes, confirming that deletions involve large segments of q16, which can reflect the BC phenotype.
SLC7A5 is widely expressed in many human cancers and various cancer cell lines . The current study included two large discovery and validation cohorts to confirm the significant association between the high SLC7A5 protein expression and the poor prognostic clinico-pathological parameters, including larger tumour size, higher grade and poor NPI. Furthermore, high SLC7A5 expression was significantly associated with higher expression of the proliferative marker (Ki67). This supports the results of previous studies which demonstrated that these two biomarkers are significantly correlated in tongue cancer , neuroendocrine carcinoma of the lung , thymic carcinoma  and breast cancer , indicating that SLC7A5 is critical for proliferation in cancer cells.
With respect to the breast cancer ER+ subtypes, SLC7A5 expression was lower in ER+ tumours that have low proliferation (luminal A subtype) compared with the highly proliferative ER+ (luminal B) malignancies, and it was primarily associated with poor patient outcome and shorter DMFS in the latter class only. This is most likely due to their heavier energy and nutrient requirements for cell survival, proliferation and metastasis. This was anticipated, as it has been shown that over expression of SLC7A5 is actively involved in the proliferation of vascular smooth muscle cells  and it is co-expressed with vascular endothelial growth factor (VEGF) in the primary and metastatic sites of many cancers , which may be implicated for the metastatic process. In this study the most significant positive correlation between mRNA expression of SLC7A5 and VEGFA was identified in the luminal B subtype. In this regard, Bartlett et al. included SLC7A5 as a part of the five-gene Mammostrat® immunohistochemistry panel, where the higher expression is used to predict recurrence-free survival (RFS), DMFS and overall survival (OS) in ER+ breast cancer during endocrine therapy . However, they did not consider the different molecular subtypes of BC.
SLC7A5 mRNA and protein was also highly expressed in TNBC and HER2+ BC, in concordance with Furuya et al. . However, in these subtypes the significant association between SLC7A5 protein expression and patient outcome was only observed in the HER2+ tumours. Among all BC subtypes, SLC7A5 protein expression was an independent predictor of short BCSS in ER+ high-proliferation tumours only. In this regard, the larger sample of ER+ high-proliferation cases might reflect the stronger significance compared with the smaller sample of HER2+ and TNBC tumours. We therefore suggest that further confirmation in larger cohorts of HER2+ and TN tumours is warranted.
Previous studies have shown regulation of SLC7A5 by other proteins including the tumour oncogene Myc, which induces SLC7A5 [4, 5]. In the current study, the relationship between SLC7A5 and other regulatory proteins in both mRNA and protein expression was investigated. A positive relationship was observed between SLC7A5 and Myc in both protein and mRNA levels, and this correlation remained significant only in luminal B subtype, when different subtypes were investigated. ATF4-dependent transcripts also encode for SLC7A5 upon amino acid deprivation  and in this study there was positive correlation between ATF4 and SLC7A5 gene expression, in line with expectations. A previous study showed that activation of the HIF2α pathway increases mTORC1 activity by upregulating expression of the amino acid carrier SLC7A5  and the current study confirmed the positive correlation between HIF2α and SLC7A5, which was only observed in luminal B tumours. SLC7A5 functions by importing essential amino acids to cancer cells and research has detailed the role of amino acids, particularly leucine, in activating mTORC1, which in turn controls protein translation and cell proliferation, and prevents apoptosis in malignant cells [2, 3]. This study showed positive correlation between SLC7A5 and mTOR at the mRNA level. However, there were conflicting results in the analysis of protein levels of SLC7A5, whereby high SLC7A5 expression was associated with lower expression of the mTORC1 phosphorylated at ser (2448), which was included in this study. This was unsurprising as Cheng et al. confirmed that phosphorylation of mTORC1 at ser (2448), which is stimulated by growth factors, was mutually exclusive with mTORC1 phosphorylated at thr (2446), which is regulated by amino acids . These observations may explain why SLC7A5 over expression is primarily associated with poor outcome only in the high proliferation ER+ tumours.
This study further investigated the association of SLC7A5 expression with other glutamine transporters, in which some variability in the expression of amino acid transporters across molecular subtypes was observed. For example, the TN subtype was the main class associated with the transporters required for glutamine influx, perhaps because it depends on delivery of glutamine instead of synthesis. In contrast, positive correlation between SLC7A5 and the glutamine synthase enzyme GLUL was observed in luminal A tumours, suggesting that this subtype might rely on glutamine neosynthesis rather than uptake. SLC1A5 functionally couples with SLC7A5 to allow the cellular influx and efflux of glutamine, as SLC1A5 mediates uptake of glutamine, while SLC7A5 uses intracellular glutamine concentrations to adjust the essential amino acid cytoplasmic pool for metabolic demands and signalling to mTORC1 . Here we observed that SLC7A5 and SLC1A5 are positively correlated in all the BC subtypes.
Previous studies have raised awareness and revealed the importance of the proline-glutamine (Pro-Gln) regulatory axis in BC. SLC7A5 appears to have a pivotal role in this regulatory axis, as its expression was highly associated with the enzymes that mediate glutamate degradation to form the amino acid proline, which has been shown to play a role in assisting tumour growth by different mechanisms .
Blocking of SLC7A5 using its inhibitor, 2-aminobicyclo-(2,2,1)-heptane-2-carboxylic acid (BCH), efficiently decreased colony formation of MDA-MB-231 TNBC cells . Even though the consequences of blocking SLC7A5 in the highly proliferative ER+ tumours remain undetermined, this study suggests that SLC7A5 can potentially be used as a therapeutic target for luminal B BC.
This study revealed and confirmed that the glutamine transporter SLC7A5 is associated with poor prognostic characteristics and poor survival outcome. Over expression of SLC7A5 appears to play a role in the proliferation and progression of the aggressive ER+ subtype of breast cancer, thus it could act as a potential therapeutic target. Functional assessment is necessary to reveal the specific role played by this amino acid transporter in the highly proliferative subclass and HER2+ BC.
Breast cancer specific survival
Chromogenic in situ hybridisation
Copy number aberration
Copy number variation
Distant metastasis-free survival
Human epidermal growth factor receptor 2
Loss of heterozygosity
Molecular Taxonomy of Breast Cancer International Consortium
Mammalian target of rapamycin complex1
Nottingham Prognostic Index
Prediction analysis of microarray 50
Single nucleotide polymorphism
The Cancer Genome Atlas
Triple negative breast cancer
Vascular endothelial growth factor
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.
Bond P. Regulation of mTORC1 by growth factors, energy status, amino acids and mechanical stimuli at a glance. J Int Soc Sports Nutr. 2016;13:8.
Bar-Peled L, Sabatini DM. Regulation of mTORC1 by amino acids. Trends Cell Biol. 2014;24(7):400–6.
Gao P, Tchernyshyov I, Chang TC, Lee YS, Kita K, Ochi T, Zeller KI, De Marzo AM, Van Eyk JE, Mendell JT, et al. c-Myc suppression of miR-23a/b enhances mitochondrial glutaminase expression and glutamine metabolism. Nature. 2009;458(7239):762–5.
Hayashi K, Jutabha P, Endou H, Anzai N. c-Myc is crucial for the expression of LAT1 in MIA Paca-2 human pancreatic cancer cells. Oncol Rep. 2012;28(3):862–6.
Reynolds MR, Lane AN, Robertson B, Kemp S, Liu Y, Hill BG, Dean DC, Clem BF. Control of glutamine metabolism by the tumor suppressor Rb. Oncogene. 2014;33(5):556–66.
Kanai Y, Segawa H, Miyamoto K, Uchino H, Takeda E, Endou H. Expression cloning and characterization of a transporter for large neutral amino acids activated by the heavy chain of 4F2 antigen (CD98). J Biol Chem. 1998;273(37):23629–32.
Jewell JL, Kim YC, Russell RC, Yu FX, Park HW, Plouffe SW, Tagliabracci VS, Guan KL. Metabolism. Differential regulation of mTORC1 by leucine and glutamine. Science (New York, NY). 2015;347(6218):194–8.
Kimball SR. Regulation of translation initiation by amino acids in eukaryotic cells. Prog Mol Subcell Biol. 2001;26:155–84.
Bonfils G, Jaquenoud M, Bontron S, Ostrowicz C, Ungermann C, De Virgilio C. Leucyl-tRNA synthetase controls TORC1 via the EGO complex. Mol Cell. 2012;46(1):105–10.
Yanagida O, Kanai Y, Chairoungdua A, Kim DK, Segawa H, Nii T, Cha SH, Matsuo H, Fukushima J, Fukasawa Y, et al. Human L-type amino acid transporter 1 (LAT1): characterization of function and expression in tumor cell lines. Biochim Biophys Acta. 2001;1514(2):291–302.
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.
Yoon JH, Kim IJ, Kim H, Kim HJ, Jeong MJ, Ahn SG, Kim SA, Lee CH, Choi BK, Kim JK, et al. Amino acid transport system L is differently expressed in human normal oral keratinocytes and human oral cancer cells. Cancer Lett. 2005;222(2):237–45.
Nakanishi K, Matsuo H, Kanai Y, Endou H, Hiroi S, Tominaga S, Mukai M, Ikeda E, Ozeki Y, Aida S, et al. LAT1 expression in normal lung and in atypical adenomatous hyperplasia and adenocarcinoma of the lung. Virchows Arch. 2006;448(2):142–50.
Fuchs BC, Bode BP. Amino acid transporters ASCT2 and LAT1 in cancer: partners in crime? Semin Cancer Biol. 2005;15(4):254–66.
Curtis C, Shah SP, Chin SF, Turashvili G, Rueda OM, Dunning MJ, Speed D, Lynch AG, Samarajiwa S, Yuan Y, et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature. 2012;486(7403):346–52.
Šidák Z. Rectangular confidence regions for the means of multivariate normal distributions. J Am Stat Assoc. 1967;62(318):626–33.
Ciriello G, Gatza ML, Beck AH, Wilkerson MD, Rhie SK, Pastore A, Zhang H, McLellan M, Yau C, Kandoth C, et al. Comprehensive molecular portraits of invasive lobular breast cancer. Cell. 2015;163(2):506–19.
Kao J, Salari K, Bocanegra M, Choi YL, Girard L, Gandhi J, Kwei KA, Hernandez-Boussard T, Wang P, Gazdar AF, et al. Molecular profiling of breast cancer cell lines defines relevant tumor models and provides a resource for cancer gene discovery. PLoS One. 2009;4(7):e6146.
Abd El-Rehim DM, Ball G, Pinder SE, Rakha E, Paish C, Robertson JF, Macmillan D, Blamey RW, Ellis IO. High-throughput protein expression analysis using tissue microarray technology of a large well-characterised series identifies biologically distinct classes of breast cancer confirming recent cDNA expression analyses. Int J Cancer. 2005;116(3):340–50.
McCarty KS Jr, McCarty KS Sr. Histochemical approaches to steroid receptor analyses. Semin Diagn Pathol. 1984;1(4):297–308.
Green AR, Aleskandarany MA, Agarwal D, Elsheikh S, Nolan CC, Diez-Rodriguez M, Macmillan RD, Ball GR, Caldas C, Madhusudan S, et al. MYC functions are specific in biological subtypes of breast cancer and confers resistance to endocrine therapy in luminal tumours. Br J Cancer. 2016;114(8):917–28.
Aleskandarany MA, Rakha EA, Ahmed MA, Powe DG, Ellis IO, Green AR. Clinicopathologic and molecular significance of phospho-Akt expression in early invasive breast cancer. Breast Cancer Res Treat. 2011;127(2):407–16.
Abdel-Fatah TM, Powe DG, Ball G, Lopez-Garcia MA, Habashy HO, Green AR, Reis-Filho JS, Ellis IO. Proposal for a modified grading system based on mitotic index and Bcl2 provides objective determination of clinical outcome for patients with breast cancer. J Pathol. 2010;222(4):388–99.
Elsheikh S, Green AR, Aleskandarany MA, Grainge M, Paish CE, Lambros MB, Reis-Filho JS, Ellis IO. CCND1 amplification and cyclin D1 expression in breast cancer and their relation with proteomic subgroups and patient outcome. Breast Cancer Res Treat. 2008;109(2):325–35.
Aleskandarany MA, Rakha EA, Macmillan RD, Powe DG, Ellis IO, Green AR. MIB1/Ki-67 labelling index can classify grade 2 breast cancer into two clinically distinct subgroups. Breast Cancer Res Treat. 2011;127(3):591–9.
Jerjees DA, Negm OH, Alabdullah ML, Mirza S, Alkaabi M, Hameed MR, Abduljabbar R, Muftah A, Nolan CC, Green AR, et al. The mammalian target of rapamycin complex 1 (mTORC1) in breast cancer: the impact of oestrogen receptor and HER2 pathways. Breast Cancer Res Treat. 2015;150(1):91–103.
Aleskandarany MA, Rakha EA, Ahmed MA, Powe DG, Paish EC, Macmillan RD, Ellis IO, Green AR. PIK3CA expression in invasive breast cancer: a biomarker of poor prognosis. Breast Cancer Res Treat. 2010;122(1):45–53.
Lancashire LJ, Powe DG, Reis-Filho JS, Rakha E, Lemetre C, Weigelt B, Abdel-Fatah TM, Green AR, Mukta R, Blamey R, et al. A validated gene expression profile for detecting clinical outcome in breast cancer using artificial neural networks. Breast Cancer Res Treat. 2010;120(1):83–93.
Green AR, Caracappa D, Benhasouna AA, Alshareeda A, Nolan CC, Macmillan RD, Madhusudan S, Ellis IO, Rakha EA. Biological and clinical significance of PARP1 protein expression in breast cancer. Breast Cancer Res Treat. 2015;149(2):353–62.
Elston CW, Ellis IO. Pathological prognostic factors in breast cancer. I. The value of histological grade in breast cancer: experience from a large study with long-term follow-up. C. W. Elston & I. O. Ellis. Histopathology 1991; 19; 403-410. Histopathology. 2002;41(3a):151–2. discussion 152–153
Senkus E, Kyriakides S, Ohno S, Penault-Llorca F, Poortmans P, Rutgers E, Zackrisson S, Cardoso F. Primary breast cancer: ESMO clinical practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2015;26(Suppl 5):v8–30.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.
Huret JL, Ahmad M, Arsaban M, Bernheim A, Cigna J, Desangles F, Guignard JC, Jacquemot-Perbal MC, Labarussias M, Leberre V, et al. Atlas of genetics and cytogenetics in oncology and haematology in 2013. Nucleic Acids Res. 2013;41(Database issue):D920–4.
Parker JS, Mullins M, Cheang MC, Leung S, Voduc D, Vickery T, Davies S, Fauron C, He X, Hu Z, et al. Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol. 2009;27(8):1160–7.
Kaira K, Oriuchi N, Imai H, Shimizu K, Yanagitani N, Sunaga N, Hisada T, Tanaka S, Ishizuka T, Kanai Y, 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-Arnaiz I, Melendez-Rodriguez F, Rodriguez-Vaello V, Marsboom G, de Carcer G, Acosta-Iborra B, Albacete-Albacete L, Ordonez A, Serrano-Oviedo L, et al. HIF2alpha acts as an mTORC1 activator through the amino acid carrier SLC7A5. Mol Cell. 2012;48(5):681–91.
Adams CM. Role of the transcription factor ATF4 in the anabolic actions of insulin and the anti-anabolic actions of glucocorticoids. J Biol Chem. 2007;282(23):16744–53.
Craze ML, Cheung H, Jewa N, Coimbra ND, Soria D, El-Ansari R, Aleskandarany MA, Cheng KW, Diez-Rodriguez M, Nolan CC. MYC regulation of glutamine-proline regulatory axis is key in luminal B breast cancer. Br J Cancer. 2018;118(2):258-265.
El Ansari R, McIntyre A, Craze ML, Ellis IO, Rakha EA, Green AR. Altered glutamine metabolism in breast cancer; subtype dependencies and alternative adaptations. Histopathology. 2018;72(2):183–90.
Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, et al. Molecular portraits of human breast tumours. Nature. 2000;406(6797):747–52.
Dawson SJ, Rueda OM, Aparicio S, Caldas C. A new genome-driven integrated classification of breast cancer and its implications. EMBO J. 2013;32(5):617–28.
Rakha EA, El-Sayed ME, Green AR, Paish EC, Powe DG, Gee J, Nicholson RI, Lee AH, Robertson JF, Ellis IO. Biologic and clinical characteristics of breast cancer with single hormone receptor positive phenotype. J Clin Oncol. 2007;25(30):4772–8.
Berx G, Cleton-Jansen AM, Strumane K, de Leeuw WJ, Nollet F, van Roy F, Cornelisse C. E-cadherin is inactivated in a majority of invasive human lobular breast cancers by truncation mutations throughout its extracellular domain. Oncogene. 1996;13(9):1919–25.
Wang Q, Holst J. L-type amino acid transport and cancer: targeting the mTORC1 pathway to inhibit neoplasia. Am J Cancer Res. 2015;5(4):1281–94.
Toyoda M, Kaira K, Ohshima Y, Ishioka NS, Shino M, Sakakura K, Takayasu Y, Takahashi K, Tominaga H, Oriuchi N, et al. Prognostic significance of amino-acid transporter expression (LAT1, ASCT2, and xCT) in surgically resected tongue cancer. Br J Cancer. 2014;110(10):2506–13.
Kaira K, Oriuchi N, Imai H, Shimizu K, Yanagitani N, Sunaga N, Hisada T, Kawashima O, Iijima H, Ishizuka T, et al. Expression of L-type amino acid transporter 1 (LAT1) in neuroendocrine tumors of the lung. Pathol Res Pract. 2008;204(8):553–61.
Kaira K, Oriuchi N, Imai H, Shimizu K, Yanagitani N, Sunaga N, Hisada T, Ishizuka T, Kanai Y, Endou H, et al. L-type amino acid transporter 1 (LAT1) is frequently expressed in thymic carcinomas but is absent in thymomas. J Surg Oncol. 2009;99(7):433–8.
Furuya M, Horiguchi J, Nakajima H, Kanai Y, Oyama T. Correlation of L-type amino acid transporter 1 and CD98 expression with triple negative breast cancer prognosis. Cancer Sci. 2012;103(2):382–9.
Liu XM, Reyna SV, Ensenat D, Peyton KJ, Wang H, Schafer AI, Durante W. Platelet-derived growth factor stimulates LAT1 gene expression in vascular smooth muscle: role in cell growth. FASEB J. 2004;18(6):768–70.
Bartlett JM, Thomas J, Ross DT, Seitz RS, Ring BZ, Beck RA, Pedersen HC, Munro A, Kunkler IH, Campbell FM, et al. Mammostrat as a tool to stratify breast cancer patients at risk of recurrence during endocrine therapy. Breast Cancer Res. 2010;12(4):R47.
Cheng SW, Fryer LG, Carling D, Shepherd PR. Thr2446 is a novel mammalian target of rapamycin (mTOR) phosphorylation site regulated by nutrient status. J Biol Chem. 2004;279(16):15719–22.
Phang JM, Liu W, Hancock CN, Fischer JW. Proline metabolism and cancer: emerging links to glutamine and collagen. Curr Opin Clin Nutr Metab Care. 2015;18(1):71–7.
Liang Z, Cho HT, Williams L, Zhu A, Liang K, Huang K, Wu H, Jiang C, Hong S, Crowe R, et al. Potential biomarker of L-type amino acid transporter 1 in breast cancer progression. Nucl Med Mol Imaging. 2011;45(2):93–102.
We thank the Nottingham Health Science Biobank and Breast Cancer Now Tissue Bank for the provision of tissue samples.
Availability of data and materials
The datasets generated during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
This study was approved by the Nottingham Research Ethics Committee 2 under the title “Development of a molecular genetic classification of breast cancer”.
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.
Table S1. Clinicopathological parameters of the METABRIC and Nottingham discovery and validation series. (DOCX 16 kb)
Figure S4. SLC7A5 mRNA expression, in the TCGA data, and its association with copy number alteration (A), staging system (B), ER status (C), PR status (D) and HER2 (E) status. (PPTX 176 kb)
Table S2. List of genes with copy number loss that was significantly associated with SLC7A5 deletion in the METABRIC and TCGA data. (DOCX 16 kb)
Figure S1. SLC7A5 gene expression and its association, using Breast Cancer Gene-Expression Miner v4.0, with: tumour grade (A), NPI (B), ER status (C), PR status (D), HER2 status (E), Triple Negative status (F)and PAM50 subtypes (G). (PPTX 117 kb)
Table S3. SLC7A5 mRNA and patient outcome. (DOCX 14 kb)
Figure S2. SLC7A5 vs DMFS in all cases in the discovery set (A), all cases in the validation set (B), ER + low proliferation tumours in the combined discovery and validation set cases (C), ER + high proliferation tumours in the combined discovery and validation set cases (D), triple negative tumours in the combined discovery and validation set cases (E) and HER2+ tumours in the combined discovery and validation set (F). (PPTX 407 kb)
Figure S3. SLC7A5 mRNA and breast cancer patient outcome using Breast Cancer Gene-Expression Miner in unselected cases (K), ER+ disease (L) and ER- disease (M). (PPTX 80 kb)
About this article
Cite this article
El Ansari, R., Craze, M.L., Miligy, I. et al. The amino acid transporter SLC7A5 confers a poor prognosis in the highly proliferative breast cancer subtypes and is a key therapeutic target in luminal B tumours. Breast Cancer Res 20, 21 (2018). https://doi.org/10.1186/s13058-018-0946-6
- Breast cancer