Insulin-like growth factor receptor and sphingosine kinase are prognostic and therapeutic targets in breast cancer
BMC Cancer volume 17, Article number: 820 (2017)
Targeting the type 1 insulin-like growth factor receptor (IGF1R) in breast cancer remains an ongoing clinical challenge. Oncogenic IGF1R-signaling occurs via activation of PI3K/AKT/MAPK downstream mediators which regulate cell proliferation and protein synthesis. To further understand IGF1R signaling we have investigated the involvement of the oncogenic IGF1R-related sphingosine kinase (SphK) pathway.
The prognostic (overall survival, OS) and therapeutic (anti-endocrine therapy) co-contribution of IGF1R and SphK1 were investigated using breast cancer patient samples (n = 236) for immunohistochemistry to measure total and phosphorylated IGF1R and SphK1. Kaplan-Meier and correlation analyses were performed to determine the contribution of high versus low IGF1R and/or SphK1 expression to OS in patients treated with anti-endocrine therapy. Cell viability and colony formation in vitro studies were completed using estrogen receptor (ER) positive and negative breast cancer cell-lines to determine the benefit of IGF1R inhibitor (OSI-906) and SphK inhibitor (SKI-II) co-therapy. Repeated measures and 1-way ANOVA were performed to compare drug treatments groups and the Chou-Talalay combination index (CI) was calculated to estimate drug synergism in vitro (CI < 1).
High IGF1R and SphK1 protein co-expression in tumor tissue was associated with improved OS specifically in ER-positive disease and stratified for anti-endocrine therapy. A significant synergistic inhibition of cell viability and/or colony formation following OSI-906 and SKI-II co-treatment in vitro was evident (p < 0.05, CI < 1).
We conclude that high IGF1R and SphK1 co-expression act together as prognostic indicators and are potentially, dual therapeutic targets for the development of a more effective IGF1R-directed combination breast cancer therapy.
Clinically targeting oncogenic signaling pathways in breast cancer, such as those initiated by the estrogen receptor (ER) and the human epidermal growth factor receptor-2 (HER2), has been highly beneficial to the treatment of the disease. However, given the heterogeneity that exists among breast cancer molecular subtypes based on the ER, progesterone receptor (PR) and HER2 status, which modulate many growth factor signaling pathways such as type 1 insulin-like growth factor receptor (IGF1R) signaling [1, 2], it is evident that singularized breast cancer targeted therapies are associated with therapeutic drawbacks such as a propensity to develop therapy resistance [3,4,5].
Specifically, the IGF1R signaling pathway has been shown to play an oncogenic role in both ER-positive and ER-negative breast cancer via the activation of downstream PI3K/AKT/MAPK/FAK signaling mediators to effectively regulate cell proliferation, migration and protein synthesis (i.e. mRNA translation) [6, 7]. However, in conflict with the oncogenic role of IGF1R signaling, lowe IGF1R expression has been reported to be associated with poorer outcomes in ER-negative breast cancer , compared to high IGF1R expression which leads to a better outcome [8, 9]. Moreover, high phosphorylated IGF1R (p-IGF1R) expression in luminal, triple-negative, and HER2 subtypes combined has been shown to be associated with a poorer survival outcome suggesting that IGF1R activation compared to expression may be more important as a prognostic factor . Despite the pre-clinical evidence suggesting that therapeutically targeting the IGF1R-pathway would be clinically effective in some patients, IGF1R monotherapies to date have not shown any improvements in clinical outcome and there is still a need to identify specific IGF1R co-related prognostic factors and therapeutic approaches [6, 11, 12]. Moreover, there is still conflicting prognostic vs. preclinical data in relation to the benefits of IGF1R targeted therapies in breast cancer which highlights the need for a better understanding of IGF1R signaling .
In addition to the ER-signaling pathway, IGF1R is known to regulate the oncogenic lipid kinase, sphingosine kinase 1 (SphK1) pathway which mediates proliferative, migratory and angiogenic effects. These effects are mediated via the intracellular and extracellular actions of the second messenger prosurvival lipid sphingosine 1-phosphate (SIP) and the SIP receptors, S1P1-S1P5 located in the plasma membrane in breast cancer [13,14,15,16,17,18]. SphK1 is known to be expressed in both ER positive and negative breast cancer and is associated with worse disease outcomes in both [19, 20].
Pre-clinical studies using SphK1-targeting therapies have shown that they possess anticancer activity, and recent findings have demonstrated that co-treatment with an epidermal growth factor receptor (EGFR) targeted therapy, gefitinib, and SphK1-targeted therapy has clinical potential in breast cancer [21,22,23,24]. Moreover, expression of IGF1R and SphK1/SIP-receptors has been shown to contribute to tamoxifen resistance in ER-positive breast cancer [16, 25, 26] which further highlights the need to better understand the significance of IGF1R and SphK1 co-expression and their contribution to anti-estrogen therapy resistance in breast cancer.
In order to further understand the prognostic and therapeutic implications of IGF1R and SphK1 co-expression in breast cancer we have analyzed their distribution in human breast cancer formalin-fixed paraffin embedded (FFPE) tissue samples. In addition we have undertaken pre-clinical in vitro studies using the dual IGF1R/insulin receptor (InsR) tyrosine kinase inhibitor, OSI-906 and the SphK inhibitor, SKI-II as a novel IGF1R-directed combination therapy. This study has identified novel relationships between breast cancer patient survival outcome and ER, PR and HER2 status and anti-estrogen therapy, based on IGF1R and SphK1 protein expression. Moreover, our evidence in vitro suggests that therapeutically co-targeting IGF1R and SphK1 has the potential for clinical benefit. In line with the findings of this study, IGF1R and SphK1 expression may have prognostic significance and co-directed combination therapies may be beneficial, specifically for ER-positive breast cancer.
Reagents and drugs
Cell culture reagents were purchased from Trace Biosciences (North Ryde, New South Wales, Australia) and Nunc (Roskilde, Denmark). Bovine insulin, methanol, calcium chloride, magnesium chloride, crystal violet powder and 1-(4,5-Dimethylthiazol-2-yl)-3,5-diphenylformazan were purchased from Sigma-Aldrich. Enhanced chemiluminescence (ECL) reagent was SuperSignal West Pico Chemiluminescent Substrate (Pierce Biotechnology). The dual IGF1R/insulin receptor tyrosine kinase inhibitor (OSI-906; also referred to as linsitinib) was purchased from MedChem Express (Princeton, NJ) and the SphK inhibitor 2-(p-hydroxyanilino)-4-(p-chlorophenyl)thiazole (SKI-II) from Calbiochem . Antibodies raised against phospho-Y1135/1136 IGFR1, IGFR1 beta chain, phospho-Ser473 AKT and total AKT, 4E-BP1 and eIF4E were purchased from Cell Signaling Technology (Beverley, MA). The antibody to detect SphK1 (ab16491) for western blots was purchased from Abcam and for SphK1 immunohistochemistry, from Abgent (AP7237c).
Breast cancer tissues were obtained from the Australian Breast Cancer Tissue Bank (ABCTB), Westmead, NSW, Australia for the purposes of this study. This study was approved by the Human Research Ethics Committee of the Northern Sydney Local Health District (Reference Numbers: RESP/15/125 and LNR/15/HAWKE/182) for the analysis of human breast cancer tissues samples obtained from the Australian Breast Cancer Tissue Bank. All samples obtained from this bank were de-identified and were from donors who had given written informed consent for their banked tumor tissue to be used in future research projects.
A total of 236 FFPE breast tissue samples were approved for use, comprised of five tissue micro-arrays (TMA) in duplicate or triplicate cores (0.6 mm3) (187 patients in total) and 49 whole face tissue sections. All patient samples had molecular subtyping from the ABCTB for ER, PR and HER2 expression by IHC and/or FISH analysis (for HER2) (Table 1). Patient information provided by the ABCTB included gender, disease status (all reported as invasive), pathology notes where applicable, primary histologic diagnosis and histopathological grade (Table 1). Patient follow-up data provided by the ABCTB consisted of diagnosis age, year of first breast event, time of follow-up since diagnosis and follow-up status (median follow up; 61 months (Table 1). In addition the ABCTB provided information relating to the therapy the patients received included the following: 1) anti-endocrine therapy; 2) HER2-therapy and 3) chemotherapy (Table 2). All studies were performed with approval from the Northern Sydney Local Health District (NSLHD) Human Research Ethics Committee (HREC), which assessed it as a low-negligible risk study.
IHC was performed on 4 μm FFPE sections using an automated tissue stainer (Autostainer, DAKO, Glostrup, Denmark) according to standard manufacturer’s operating procedures. Antigen retrieval was performed using a water bath heated to 99.2 °C for 20 min in freshly made 10 mM citric acid monohydrate adjusted to pH 6.0. The sections were quenched in 0.3% hydrogen peroxide for 5 min, blocked with 5% goat serum for 30 min and incubated in primary antibodies: IGF1Rβ antibody (no cross-reaction with the insulin receptor (InsR)) (#3027, Cell Signaling, Danvers, MA, USA 1:100), p-IGF1R (#ab39398, Abcam, Melbourne, VIC, Australia, 1:200) and SphK1 (#AP7237c, Abgent, San Diego, CA, USA, 1:200) for one hour at room temperature. Protein detection was subsequently performed using the DAKO-Envision Dual Link Labelled Polymer (Anti-Rabbit) (#K5007, Dako, Botany, NSW, Australia) for 30 min and the ImmPACT NovaRed Peroxidase Substrate Kit (#SK-4805, Vector Laboratories, Burlingame, CA, USA) for 10 min at room temperature. All antibodies were optimized using a series of dilutions on a TMA comprised of ten ER-positive and negative breast cancer patient tissues in duplicate to determine an optimal dilution for IHC staining. Final dilutions were closely assessed for specific membranous, cytoplasmic and/or nuclear staining in line with the literature. Negative controls were included for all IHC using a rabbit immunoglobulin fraction (#X0936, DAKO) at the final concentrations of the primary antibodies and a tissue sample incubated with the anti-rabbit antibody in the absence of a primary antibody.
Manual scoring was assessed on all samples for subsequent statistical analysis, with examples shown in Fig. 1. IGF1R and p-IGF1R staining expression levels were manually assessed using the HER2 scoring system described in the Hercep Test manual (DAKO) as follows: no staining = 0, faint staining = 1, weak to moderate staining = 2 and strong staining = 3, in line with published studies [27, 28]. Positive staining was defined as membrane/cytoplasmic and/or nuclear staining detectable in ≥10% of tumor epithelial cells. SphK1 staining levels were manually assessed as: no staining (<10% of tumor epithelial cells with cytoplasmic staining) = 0, weak = 1, moderate = 2 and strong = 3, in accordance with previously published data [29,30,31].
Kaplan-Meier (KM) survival analysis was performed using SPSS v.22 (IBM, Armonk, NY, USA) to determine prognostic significance of high protein expression of phospho- and total IGF1R, and SphK1, for overall survival (OS) and disease-free survival (DFS) in non-stratified patient cohorts, and after stratification for high or low ER, PR, and HER2. KM analysis was also performed to determine the relationship between IGF1R and SphK1 protein expression in patients stratified for anti-endocrine therapy treatment. For KM analysis of IGF1R, p-IGF1R, and SphK1, staining scores described above, on a scale of 0 to 3, were converted to a new binary scale where 0 or 1 = low, and 2 or 3 = high. For analysis of IGF1R-SphK1 and p-IGF1R-SphK1 co-expression, a binary co-expression score was used: either or both analytes low = low co-expression, both analytes high = high co-expression. Significance (P < 0.05) was determined by the log-rank test. Multivariate survival analysis was performed in SPSS by the Cox proportional hazards method, using forward stepwise (likelihood ratio) regression, with significance set at P < 0.05. To determine the relationship between IGF1R and SphK1 in non-stratified and stratified patient groups including: ER+/−, PR+/− and HER2+/−, age, tumor grade correlation analysis was performed in SPSS using Spearman correlation analysis with statistical significance set at P < 0.05. All data used in these analyses are included as Additional file 1: Table S1.
Breast cancer cell-lines and culture
The ER-positive MCF7 and T47D and ER-negative HCC1806 and HCC70 breast cancer cell-lines were purchased from the American Tissue Culture Collection (Manassas, VA, USA) and cultured and maintained in phenol-red RPMI-1640 medium supplemented with 5% FBS, 15 mM Hepes and 10 mg/mL bovine insulin at 37 °C in humidified 5% CO2 atmosphere.
MTT-assay and Clonogenic survival assay
The MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) assay and clonogenic survival assay were performed by plating 2 × 103 MCF7, T47D, HCC1806 and HCC70 cells per well of a 96-well plate (MTT-assay) or a 6-well plate (clonogenic assay) in phenol-red RPMI culture media supplemented with 5% FBS and culturing for 24 h. The cells were treated with culture media containing 5% FBS in addition to OSI-906 (0, 0.1, 0.4, 1.6 or 6.4 μM) and/or SKI-II (0, 0.16, 0.8, 4, 10 or 20 μM) in combination for 96 h for MTT-assay and 10-14 days for clonogenic assay. For the MTT-assay 30 μl of 1-(4,5-dimethylthiazol-2-yl)-3,5-diphenylformazan (thiazolyl blue formazan) solution prepared at 2.5 mg/ml, in PBS (containing 0.9 mM calcium chloride plus 0.5 mM magnesium chloride) was added per well and incubated at 37 °C for 4 h. The cells were solubilized in 150 μl of DMSO for 15 min at room temperature (RT). Absorbance was read using a plate reader at 460 nm. For the clonogenic assay, the cells were washed twice in PBS (0.9 mM calcium chloride plus 0.5 mM magnesium chloride) solution at RT, fixed in 0.1% crystal violet prepared in final 20% methanol solution and destained in tap water. The cells were left to air-dry overnight and images obtained using FujiFilm Luminescent Image Analyzer LAS-300 and single colonies were counted using open colony forming unit (CFU) software (http://opencfu.sourceforge.net/).
3 × 105 MCF7 and HCC-1806 breast cancer cell-lines were plated per well of 6-well plates, cultured for 24 h and subsequently treated with OSI-906 (0.1, 0.4 or 1.6 μM) and/or SKI-II (4 μM) for 24 h. Protein lysates were prepared and 20 μg of protein was further run on SDS-polyacrylamide gels, transferred to Hybond C nitrocellulose and probed with antibodies to detect phospho and/or IGF1R, AKT, 4E-BP1, eIF4E (Cell Signaling) and SphK1 (Abgent) steady-state protein levels as previously described . Β-actin antibody was used as a loading control. Subsequent to the addition of the ECL-reagent, a FujiFilm Luminescent Image Analyzer LAS-300 (Stamford, CT) was used for band detection and image production.
Statistical analysis (cell culture)
The effect of drug treatments on MTT and clonogenic assays was first analysed across the full OSI-906 dose-range using 2-way ANOVA for repeated measures (SPSS) with experiment and SKI-II dose as factors, and OSI-906 dose as the repeated measure. Comparison of individual dose combinations was performed using 1-way ANOVA with Tukey’s post-hoc test. (GraphPad Prism v.7, La Jolla, CA, USA). Statistical significance was defined as a p-value <0.05*, p < 0.01**, p < 0.001*** and p < 0.0001****. Drug synergism was determined using the Chou Talalay method to calculate the combination index (CI) <1 .
IGF1R expression is positively associated with overall survival in breast cancer
In FFPE tissue sections, high total IGF1R protein expression, in isolation or combined with SphK1, was associated with a better overall survival (OS) rate when analyzed across all patients suggesting IGF1R, alone or together with SphK1, acts as a positive prognostic indicator (Table 3, p = 0.018, Fig. 2a and p = 0.028,; Fig. 2b, respectively). These highly significant effects on OS were not seen for DFS (data not shown). No significant relationship between SphK1 expression alone and OS was identified (Table 3). In contrast, no significant improvement in OS was detected in patients with high p-IGF1R protein expression alone or combined with SphK1 in our patient analysis (Table 3, p = 0.303, Fig. 2c and p = 0.118, Fig. 2d, respectively). However, pIGF1R and SphK1 high co-expression did lead to lower p-value compared to either protein in isolation (Table 3). After stratification for ER status we found, in line with previous studies, an improved OS in ER-positive breast cancer patients who express high IGF1R protein in isolation (Table 3, p = 0.048, Fig. 3a) or combined with SphK1 (Table 3, p = 0.051, Fig. 3b), compared to ER-negative breast cancers (Table 3, p = 0.445, Fig. 3c and p = 0.582, Fig. 3d, respectively). This effect was further supported by a strong positive correlation between IGF1R and SphK1 in ER-positive (p = 0.001), but not ER-negative (p = 0.936) breast cancer patients (Table 4). In support of the literature reporting that pIGF1R can be detected in the nucleus of the cell, bind DNA and act as a transcription factor, we also found that some patient samples showed positive nuclear IHC staining for pIGF1R (Additional file 2) .
Association between IGF1R and SphK1 expression, and ER, PR, HER2, tumor grade, and age
After stratification of tumors as either HER2 +/− or PR +/−, there was no significant prognostic effect of IGF1R (Additional file 3A-D) or p-IGF1R (Additional file 4A-D) in isolation, or combined with SphK1 expression (data not shown), in our patient cohort. Similarly, no significance was observed for high vs. low SphK1 protein expression alone following stratification for HER2 or PR (Table 3). However, as shown in Table 4, a significant association was observed between both p-IGF1R and IGF1R expression and ER-positivity (p = 0.007 and p < 0.001, respectively), and between IGF1R expression and PR-positivity (p < 0.001). Furthermore, p-IGF1R expression was negatively correlated to HER2-positive expression (Table 4, p = 0.002), yet this was not evident for IGF1R expression. Moreover, SphK1 high expression was inversely associated with both ER-positivity (p < 0.001) and PR-positivity (p = 0.003) in breast tumors (Table 4) which supports literature findings . We were not able to identify any association between IGF1R or SphK1 and tumor grade; however, p-IGF1R was inversely correlated to tumor grade in our analysis (Table 4, p = 0.004). Lastly, we identified an inverse correlation of IGF1R expression to following stratification for age (i.e. ≥51y; average age of menopause) in our analysis (Table 4, p = 0.025).
To determine the combined effect of the IHC measurements and hormone receptor status on patient survival we undertook multivariate Cox regression (proportional hazards) analysis, using the forward stepwise (likelihood ratio) method. Analyzing the influence of IGF1R, p-IGF1R and SphK1 alone on overall patient survival, only IGF1R was significant (Table 5; p = 0.008). None of these variables was significant for DFS (Table 5). When ER and PR status were added to the model, IGF1R remained the only significant variable (Table 5: p = 0.010 for OS, not significant for DFS). However, adding ER, PR and HER status to the model, HER2 became the dominant variable when OS was the survival endpoint (Table 5; p < 0.001), with IGF1R expression approaching significance (Table 5; p = 0.053). For DFS, only HER2 was significant among the six variables (Table 5; p = 0.007).
IGF1R and SphK1 co-expression is associated with improved disease outcome in anti-endocrine therapy treated breast cancer patients
Given that we identified a relationship between high IGF1R protein expression alone and combined with high SphK1 expression, and overall patient survival, that was most evident for ER-positive breast cancers, we undertook further analysis to determine whether a similar relationship existed when samples were stratified for anti-endocrine therapy. KM analysis identified that anti-endocrine therapy is associated with an improved OS (Table 3, p = 0.006, Fig. 4a) in the unstratified cohort. Examining only samples from anti-endocrine therapy treated women, high IGF1R expression (Table 3, p = 0.028, Fig. 4b), and high IGF1R/SphK1 co-expression (Table 3, p = 0.034, Fig. 4c) were both prognostic for improved OS, whereas high SphK1 expression was not (Table 3, p = 0.101, Fig. 4d). This supports our previous findings that high IGF1R/SphK1 co-expression is associated with an improved disease outcome specifically in ER-positive breast cancer patients in relation to anti-endocrine therapy.
Co-targeting IGF1R and SphK1 acts synergistically on breast cancer cell viability and colony formation
The observation that high IGF1R and SphK1 expression are prognostic for improved overall survival is paradoxical, given that the literature suggests IGF1R and SphK1 are oncogenic mediators in breast cancer. Since the positive prognostic effect was only significant for ER-positive cancers, we undertook in vitro studies to compare whether therapeutically co-targeting IGF1R and SphK1 using the dual IGF1R/InR inhibitor OSI-906 and the SphK inhibitor SKI-II might be a viable clinical approach for ER-negative vs. ER-positive cancers. To determine optimal dose ranges to identify a combination effect of the two drugs, initial MTT-assay experiments were performed on each drug alone, i.e. OSI-906 at 0.1, 0.4, 1.6 and 6.4 μM (Fig. 5) and SKI-II at 0.16, 0.8, 4, 10 and 20 μM (Additional file 5), as well as all combinations of these doses (data not shown). Using immunoblot analysis we confirmed inhibition of p-IGF1R phosphorylation and downstream IGF1R signaling (i.e. p-AKT and 4E-BP1 hyper-phosphorylation) following treatment with OSI-906 and SphK1 steady-state protein levels by SphK1 in the MCF7 ER-positive and HCC1806 ER-negative cell-line to demonstrate effective drug-target inhibition by the therapies (Additional file 6). Based on the initial MTT-experiments we identified that a fixed concentration of SKI-II at 4 μM in combination with OSI-906 in the range of 0.1-6.4 μM was the most effective in reducing cell viability (Fig. 5).
OSI-906 showed dose-dependent inhibition of the ER-positive MCF7 breast cancer cell line by both cell-viability and colony-formation assays, with a near-maximal effect seen at 1.6 μM (Fig. 5a and e). The addition of 4 μM SKI-II significantly increased sensitivity to OSI-906 by 2- to 4-fold in both assays (p < 0.001 by repeated measures ANOVA for both assays), with the greatest effect of the drug combination, compared to individual treatments, seen at 4 μM SKI-II and 0.4 μM OSI-906 (Fig. 5a; p < 0.01 and Fig. 5e; p < 0.001, and Additional file 7A). The ER-positive T47D cell-line was somewhat less sensitive to OSI-906, but a similar sensitising effect of 4 μM SKI-II was seen (Fig. 5b and f; p < 0.001 by repeated measures ANOVA for both assays). Greatest combination effects, compared to individual treatments, were seen at 4 μM SKI-II and 1.6 μM OSI-906 co-treatment in the cell-viability assay (Fig. 5f; p = 0.052), and at 4 μM SKI-II and 0.4 μM OSI-906 co-treatment in the colony-formation assay (Fig. 5f; p < 0.01 and Additional file 7B).
The ER-negative HCC1806 cell line showed similar OSI-responsiveness to T47D, and similarly, the addition of 4 μM SKI-II significantly increased responsiveness (Fig. 5c and g; P < 0.001 by repeated measures ANOVA for both assays). The greatest combination effect of SKI-II and OSI-906 co-treatment, compared to individual treatments, was observed on colony formation at the highest dose of 6.4 μM OSI-906 (Fig. 5g; p < 0.05 and Additional file 7C). ER-negative HCC70 cells were highly resistant to inhibition by OSI-906 (Fig. 5d), and even in the presence of 4 μM SKI-II failed to achieve 50% inhibition of cell viability at the highest tested OSI-906 dose of 6.4 μM. HCC70 cells expressed low levels of p-IGF1R and SphK1 protein (data not shown) which is likely a contributing factor to the lack of co-treatment effectiveness.
Calculation of the Chou-Talalay combination index confirmed our findings that co-treatment with OSI-906 and SKI-II led to synergistic effects using OSI-906 at 0.4 μM combined with SKI-II at 4 μM on cell viability and/or colony formation in the MCF7 and T47D ER-positive cell-lines, and at the higher dose of OSI-906 at 6.4 μM in the HCC1806 ER-negative cell-line (i.e. CI < 1; synergism; Table 6). Drug synergy could not be calculated for HCC70 cells owing to their strong resistance to OSI-906.
Given that IGF1R is a well-documented oncogenic factor in many different types of cancers, and therapeutically targeting its activity in isolation has proven to be unsuccessful, a better understanding of its signaling activity is certainly required . Clinically ineffective IGF1R-directed therapies may be associated with the ability of the cancer cell to re-activate the signaling pathway in addition to alternative activation of downstream IGF1R signaling pathways by other growth-promoting receptors in the cell such as HER2/3, EGFR and InsR [12, 36]. This study has provided novel findings demonstrating that the high co-expression of IGF1R and SphK1 may have positive prognostic significance for overall (but not disease free) survival in ER-positive breast cancer, but in line with pre-clinical findings demonstrating IGF1R and SphK1 can promote oncogenic effects in cancer, paradoxically may still be considered for therapeutic targeting .
We have recently reported that high IGF1R and SphK1 gene co-expression, based on data derived from public databases, is associated with a worse survival outcome in ER-positive, but not ER-negative breast cancer . However, our immunohistochemistry findings in this study report the opposite effect for protein expression (i.e. improved OS in ER-positive disease), further highlighting the ongoing discrepancy in findings relating to the relationship between high IGF1R expression and ER-positive breast cancer potentially based on different prognostic parameters (i.e. protein vs gene, patient tumor cohorts, antibodies, cut-off and pIGF1R vs. IGF1R measurements i.e. signalling vs. expression) . In light of the differences between protein versus gene expression data, consideration of whether measurement of gene expression is an accurate indicator of clinical outcomes in women with breast cancer requires further investigation. Importantly, mRNA expression levels have been shown to be a poor predictor of protein levels and in fact translational control is a more accurate predictor of protein expression [7, 37], highlighting that protein abundance rather than mRNA expression is likely to be a more accurate predictor of patient outcomes. Furthermore, expression of IGF1R vs. activity of the signaling pathway (i.e. pIGF1R and activation of IGF1R signaling) play an integral oncogenic role in cancer and likely add to the discrepancy between studies and analysis (i.e. gene vs protein) reporting the prognostic and therapeutic significance of whether blocking IGF1R activity and/or expression is more beneficial in cancer .
In this study p-IGF1R was not a prognostic factor yet Kaplan-Meier analyses revealed that p-IGF1R and SphK1 high co-expression did improve the log-rank p-value compared to either protein in isolation, suggesting that the activity of p-IGF1R signaling may be important and supporting the literature that active IGF1R signaling leads to upregulation of SphK1 expression [38, 39]. Furthermore, we did detect a significant positive correlation of both p-IGF1R and IGF1R to both ER expression and IGF1R to PR expression , and p-IGF1R was additionally negatively correlated to HER2 expression. Nevertheless, high p-IGF1R expression has been reported to be associated with a worse survival outcome in breast cancer  and based on the literature and this study, there are clearly some conflicting data in relation to whether IGF1R and/or p-IGF1R is the more useful breast cancer prognostic marker. Since pIGF1R can be detected also in the nucleus of the cell and act as a transcription factor, the prognostic significance of membranous/cytoplasmic vs. nuclear staining may be of clinical significance and needs to be further investigated  Quantitation of p-IGF1R is potentially confounded by pre-analytical factors such as tissue storage and processing times, since active phosphatases have the ability to cause uncontrolled dephosphorylation. This issue is avoided by the measurement of total IGF1R.
In line with improved therapeutic outcomes associated with ER-directed cancer therapies in ER-positive versus ER-negative breast cancers, our study supports the notion that individuals who receive anti-endocrine targeted-therapies (i.e. diagnosed as ER-positive) have a much better outcome than non-anti-endocrine targeted-therapies . Moreover, a better prognosis has been shown to be associated with non-triple-negative versus triple-negative breast cancer , as expected since ER-positive disease at present has better clinical management. IGF1R mRNA has been reported as a good prognostic factor specifically in luminal breast cancer subtypes which was also correlated to IGF1R protein expression [9, 40, 43]. However, within the luminal A subtype, high IGF1R mRNA expression is associated with a worse patient outcome than low expression and has been suggested to be a contributing factor to anti-estrogen therapy resistance .
Based on the literature linking both IGF1R and SphK1 to ER-positive breast cancer, we also demonstrated that IGF1R and SphK1 are positively correlated in ER-positive but not ER-negative disease. Our findings support the literature reporting that higher SphK1 expression is evident in ER-negative breast cancer (i.e. inverse correlation of SphK1 to ER status), compared to ER-positive disease . This finding was also extended to PR status, which supports the biology of PR as an ER-regulated gene . Interestingly we further detected an inverse correlation of IGF1R to age (i.e. average age of menopausal status) and p-IGF1R to tumor grade and HER2 status. Moreover, IGF1R expression is reported to be significantly associated to HER2-positivity and poorer disease-free survival in premenopausal women, suggesting a link between IGF1R and HER2 specifically in pre-menopausal breast cancer . In this context it is notable that multivariate survival analysis on our predominantly post-menopausal cohort found that high IGF1R expression was the only significant variable associated with overall survival if HER2 status was excluded from the model, but that HER2 status assumed the dominant effect when added to the analysis, with high IGF1R expression as the only other measured factor approaching significance. This may indicate that the variable findings in the literature in regard to the prognostic value of IGF1R measurement could be influenced by concomitant HER2 status, and suggests that further examination of the effects of IGF1R and HER2 co-expression is warranted. Serum IGF-1 has also been reported to be associated with mammographic density in pre-menopausal women and this may be linked to IGF1R activity .
It is well accepted that IGF1R can regulate ER, and SphK1 is regulated by estrogen, in breast cancer [13, 47, 48], hence the therapeutic implication that co-targeting them may be clinically effective in ER-positive breast cancer. Despite the lack of patient survival significance in relation to high IGF1R and Sphk1 co-expression in our ER-negative patient cohort, IGF1R and SphK1 are expressed in ER-negative cells and may have therapeutic significance [35, 49]. Our findings show that, despite high IGF1R and/or SphK1 being associated with a positive clinical outcome in ER-positive breast cancer while having no effect in ER-negative breast cancer, IGF1R and SphK1 co-targeting is worth investigating for clinical benefit in these women. The molecular basis to this finding may be associated with the improvement in anti-endocrine therapy in women who express IGF1R and SphK1. Further prognostic studies using a recurrent disease population vs. non-recurrent ER-positive breast cancer population would be helpful since both IGF1R and SphK1 are reported to contribute to tamoxifen resistance [16, 25, 26]. Moreover, since both IGF1R and SphK1 are reported to be associated with anti-estrogen therapy resistance [16, 25, 26], it is important that future studies using IGF1R-targeted therapies which have not been successful as monotherapies also investigate the benefit of including other IGF1R co-related factors such as SphK1 as targets.
We conducted in vitro studies in ER-positive and -negative cell lines in an attempt to dissect the different prognostic potential of IGF1R expression in ER-positive and –negative cancers. The two ER-negative cell lines tested were both more resistant to OSI-906 than the ER-positive lines, but the synergistic co-inhibition of cell viability by IGF1R and SphK1 inhibition was at least as effective in the ER-negative line HCC1806 as in the two ER-positive lines, providing a possible rationale for evaluating a combination of IGF1R and SphK1 inhibitors in vivo, particularly in ER-negative tumors. IGF1R signaling can occur in response to the formation of either a homodimer or a hybrid receptor (i.e. heterodimer) with the InsR. Since OSI-906 is a dual IGF1R/InsR inhibitor and both receptors are potential therapeutic targets in breast cancer  it is important to consider the inhibitory contribution of both receptors to OSI-906 effects. However, since we did not use a totally selective InsR inhibitor or measure changes to p-InsR it is difficult to ascertain the contribution of OSI-906 mediated inhibition of InsR vs IGF1R. Nevertheless, we have shown that OSI-906 reduces both p-IGF1R steady-state protein and IGF1R signaling (i.e. p-AKT and 4E-BP-1 hyper-phosphorylation) in the MCF7 ER-positive and HCC1806 ER-negative breast cancer cell-lines, consistent with the assumption that the therapeutic inhibition of OSI-906 is mediated at least in part via IGF1R.
In this study we have provided new evidence demonstrating a prognostic relationship between IGF1R and the co-related SphK1 signaling pathway in relation to: 1) ER, PR and HER2 status and 2) anti-endocrine targeted therapies. Moreover, our preclinical findings suggest that co-targeting IGF1R and SphK1 may have benefit in some women with breast cancer and requires further investigation.
In conclusion the findings in this study support the proposal that high tissue IGF1R levels, alone and in combination with high SphK1 levels, are prognostic indicators in ER-positive breast cancer and in response to anti-endocrine therapy. Furthermore, since IGF1R is an established oncogenic factor in breast cancer and has been the focus of numerous clinical trials, it is potentially of therapeutic importance that we have shown that co-targeting both IGF1R and SphK1 in vitro has a synergistic benefit. There is an ongoing need to understand how to effectively target IGF1R in breast cancer and to identify new IGF1R-related co-factors for clinical intervention. This study may provide a better understanding of IGF1R-signaling and help to explain the lack of clinical outcome benefit of IGF1R-directed monotherapies in breast cancer.
Perou CM. Molecular stratification of triple-negative breast cancers. Oncologist. 2011;15(Suppl 5):39–48.
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.
Andre F, Slimane K, Bachelot T, Dunant A, Namer M, Barrelier A, Kabbaj O, Spano JP, Marsiglia H, Rouzier R, et al. Breast cancer with synchronous metastases: trends in survival during a 14-year period. J Clin Oncol. 2004;22(16):3302–8.
Dawood S, Broglio K, Buzdar AU, Hortobagyi GN, Giordano SH. Prognosis of women with metastatic breast cancer by HER2 status and trastuzumab treatment: an institutional-based review. J Clin Oncol. 2010;28(1):92–8.
Pal SK, Childs BH, Pegram M. Triple negative breast cancer: unmet medical needs. Breast Cancer Res Treat. 2011;125(3):627–36.
Ochnik AM, Baxter RC. Combination therapy approaches to target insulin-like growth factor receptor signaling in breast cancer. Endocr Relat Cancer. 2016;23(11):R513–36.
Ochnik AM, Peterson MS, Avdulov SV, Oh AS, Bitterman PB, Yee D. Amplified in Breast Cancer Regulates Transcription and Translation in Breast Cancer Cells. Neoplasia. 2016;18(2):100–10.
Aaltonen KE, Rosendahl AH, Olsson H, Malmstrom P, Hartman L, Ferno M. Association between insulin-like growth factor-1 receptor (IGF1R) negativity and poor prognosis in a cohort of women with primary breast cancer. BMC Cancer. 2014;14:794.
Fu P, Ibusuki M, Yamamoto Y, Hayashi M, Murakami K, Zheng S, Iwase H. Insulin-like growth factor-1 receptor gene expression is associated with survival in breast cancer: a comprehensive analysis of gene copy number, mRNA and protein expression. Breast Cancer Res Treat. 2011;130(1):307–17.
Law JH, Habibi G, Hu K, Masoudi H, Wang MY, Stratford AL, Park E, Gee JM, Finlay P, Jones HE, et al. Phosphorylated insulin-like growth factor-i/insulin receptor is present in all breast cancer subtypes and is related to poor survival. Cancer Res. 2008;68(24):10238–46.
Yang Y, Yee D. Targeting insulin and insulin-like growth factor signaling in breast cancer. J Mammary Gland Biol Neoplasia. 2012;17(3-4):251–61.
Farabaugh SM, Boone DN, Lee AV. Role of IGF1R in breast cancer subtypes, Stemness, and lineage differentiation. Front Endocrinol. 2015;6:59.
Fagan DH, Yee D. Crosstalk between IGF1R and estrogen receptor signaling in breast cancer. J Mammary Gland Biol Neoplasia. 2008;13(4):423–9.
Pitson SM, Moretti PA, Zebol JR, Lynn HE, Xia P, Vadas MA, Wattenberg BW. Activation of sphingosine kinase 1 by ERK1/2-mediated phosphorylation. EMBO J. 2003;22(20):5491–500.
Maiti A, Takabe K, Hait NC. Metastatic triple-negative breast cancer is dependent on SphKs/S1P signaling for growth and survival. Cell Signal. 2017;32:85–92.
Maczis M, Milstien S, Spiegel S. Sphingosine-1-phosphate and estrogen signaling in breast cancer. Adv Biol Regul. 2016;60:160–5.
Antoon JW, Meacham WD, Bratton MR, Slaughter EM, Rhodes LV, Ashe HB, Wiese TE, Burow ME, Beckman BS. Pharmacological inhibition of sphingosine kinase isoforms alters estrogen receptor signaling in human breast cancer. J Mol Endocrinol. 2011;46(3):205–16.
Takabe K, Kim RH, Allegood JC, Mitra P, Ramachandran S, Nagahashi M, Harikumar KB, Hait NC, Milstien S, Spiegel S. Estradiol induces export of sphingosine 1-phosphate from breast cancer cells via ABCC1 and ABCG2. J Biol Chem. 2010;285(14):10477–86.
Ohotski J, Long JS, Orange C, Elsberger B, Mallon E, Doughty J, Pyne S, Pyne NJ, Edwards J. Expression of sphingosine 1-phosphate receptor 4 and sphingosine kinase 1 is associated with outcome in oestrogen receptor-negative breast cancer. Br J Cancer. 2012;106(8):1453–9.
Ruckhaberle E, Rody A, Engels K, Gaetje R, von Minckwitz G, Schiffmann S, Grosch S, Geisslinger G, Holtrich U, Karn T, et al. Microarray analysis of altered sphingolipid metabolism reveals prognostic significance of sphingosine kinase 1 in breast cancer. Breast Cancer Res Treat. 2008;112(1):41–52.
Martin JL, de Silva HC, Lin MZ, Scott CD, Baxter RC. Inhibition of insulin-like growth factor-binding protein-3 signaling through sphingosine kinase-1 sensitizes triple-negative breast cancer cells to EGF receptor blockade. Mol Cancer Ther. 2014;13(2):316–28.
Pitman MR, Powell JA, Coolen C, Moretti PA, Zebol JR, Pham DH, Finnie JW, Don AS, Ebert LM, Bonder CS, et al. A selective ATP-competitive sphingosine kinase inhibitor demonstrates anti-cancer properties. Oncotarget. 2015;6(9):7065–83.
Pitman MR, Woodcock JM, Lopez AF, Pitson SM. Molecular targets of FTY720 (fingolimod). Curr Mol Med. 2012;12(10):1207–19.
Martin JL, Julovi SM, Lin MZ, de Silva HC, Boyle FM, Baxter RC. Inhibition of basal-like breast cancer growth by FTY720 in combination with epidermal growth factor receptor kinase blockade. Breast Cancer Res. 2017;19(1):90.
Fagan DH, Uselman RR, Sachdev D, Yee D. Acquired resistance to tamoxifen is associated with loss of the type I insulin-like growth factor receptor: implications for breast cancer treatment. Cancer Res. 2012;72(13):3372–80.
Watson C, Long JS, Orange C, Tannahill CL, Mallon E, McGlynn LM, Pyne S, Pyne NJ, Edwards J. High expression of sphingosine 1-phosphate receptors, S1P1 and S1P3, sphingosine kinase 1, and extracellular signal-regulated kinase-1/2 is associated with development of tamoxifen resistance in estrogen receptor-positive breast cancer patients. Am J Pathol. 2010;177(5):2205–15.
Jaillardon L, Abadie J, Godard T, Campone M, Loussouarn D, Siliart B, Nguyen F. The dog as a naturally-occurring model for insulin-like growth factor type 1 receptor-overexpressing breast cancer: an observational cohort study. BMC Cancer. 2015;15:664.
Shin SJ, Gong G, Lee HJ, Kang J, Bae YK, Lee A, Cho EY, Lee JS, Suh KS, Lee DW, et al. Positive expression of insulin-like growth factor-1 receptor is associated with a positive hormone receptor status and a favorable prognosis in breast cancer. J Breast Cancer. 2014;17(2):113–20.
Liu G, Zheng H, Zhang Z, Wu Z, Xiong H, Li J, Song L. Overexpression of sphingosine kinase 1 is associated with salivary gland carcinoma progression and might be a novel predictive marker for adjuvant therapy. BMC Cancer. 2010;10:495.
Li J, Guan HY, Gong LY, Song LB, Zhang N, Wu J, Yuan J, Zheng YJ, Huang ZS, Li M. Clinical significance of sphingosine kinase-1 expression in human astrocytomas progression and overall patient survival. Clin Cancer Res. 2008;14(21):6996–7003.
Matula K, Collie-Duguid E, Murray G, Parikh K, Grabsch H, Tan P, Lalwani S, Garau R, Ong Y, Bain G, et al. Regulation of cellular sphingosine-1-phosphate by sphingosine kinase 1 and sphingosine-1-phopshate lyase determines chemotherapy resistance in gastroesophageal cancer. BMC Cancer. 2015;15:762.
Martin JL, Lin MZ, McGowan EM, Baxter RC. Potentiation of growth factor signaling by insulin-like growth factor-binding protein-3 in breast epithelial cells requires sphingosine kinase activity. J Biol Chem. 2009;284(38):25542–52.
Chou TC. Drug combination studies and their synergy quantification using the Chou-Talalay method. Cancer Res. 2010;70(2):440–6.
Sarfstein R, Werner H. Minireview: nuclear insulin and insulin-like growth factor-1 receptors: a novel paradigm in signal transduction. Endocrinology. 2013;154(5):1672–9.
Datta A, Loo SY, Huang B, Wong L, Tan SS, Tan TZ, Lee SC, Thiery JP, Lim YC, Yong WP, et al. SPHK1 regulates proliferation and survival responses in triple-negative breast cancer. Oncotarget. 2014;5(15):5920–33.
Beckwith H, Yee D. Insulin-like growth factors, insulin, and growth hormone signaling in breast cancer: implications for targeted therapy. Endocr Pract. 2014;20(11):1214–21.
Schwanhausser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, Chen W, Selbach M. Global quantification of mammalian gene expression control. Nature. 2011;473(7347):337–42.
Granata R, Trovato L, Lupia E, Sala G, Settanni F, Camussi G, Ghidoni R, Ghigo E. Insulin-like growth factor binding protein-3 induces angiogenesis through IGF-I- and SphK1-dependent mechanisms. J Thromb Haemost. 2007;5(4):835–45.
Dai L, Qi Y, Chen J, Kaczorowski D, Di W, Wang W, Xia P. Sphingosine kinase (SphK) 1 and SphK2 play equivalent roles in mediating insulin's mitogenic action. Mol Endocrinol. 2014;28(2):197–207.
Kim JH, Cho YH, Park YL, Sohn JH, Kim HS. Prognostic significance of insulin growth factor-I receptor and insulin growth factor binding protein-3 expression in primary breast cancer. Oncol Rep. 2010;23(4):989–95.
Peiro G, Adrover E, Sanchez-Tejada L, Lerma E, Planelles M, Sanchez-Paya J, Aranda FI, Giner D, Gutierrez-Avino FJ. Increased insulin-like growth factor-1 receptor mRNA expression predicts poor survival in immunophenotypes of early breast carcinoma. Mod Pathol. 2011;24(2):201–8.
Bahhnassy A, Mohanad M, Shaarawy S, Ismail MF, El-Bastawisy A, Ashmawy AM, Zekri AR. Transforming growth factor-beta, insulin-like growth factor I/insulin-like growth factor I receptor and vascular endothelial growth factor-a: prognostic and predictive markers in triple-negative and non-triple-negative breast cancer. Mol Med Rep. 2015;12(1):851–64.
Fu P, Ibusuki M, Yamamoto Y, Yamamoto S, Fujiwara S, Murakami K, Zheng S, Iwase H. Quantitative determination of insulin-like growth factor 1 receptor mRNA in formalin-fixed paraffin-embedded tissues of invasive breast cancer. Breast Cancer. 2012;19(4):321–8.
Nardulli AM, Greene GL, O'Malley BW, Katzenellenbogen BS. Regulation of progesterone receptor messenger ribonucleic acid and protein levels in MCF-7 cells by estradiol: analysis of estrogen's effect on progesterone receptor synthesis and degradation. Endocrinology. 1988;122(3):935–44.
Sun WY, Yun HY, Song YJ, Kim H, Lee OJ, Nam SJ, Koo JS. Insulin-like growth factor 1 receptor expression in breast cancer tissue and mammographic density. Mol Clin Oncol. 2015;3(3):572–80.
Diorio C, Pollak M, Byrne C, Masse B, Hebert-Croteau N, Yaffe M, Cote G, Berube S, Morin C, Brisson J. Insulin-like growth factor-I, IGF-binding protein-3, and mammographic breast density. Cancer Epidemiol Biomark Prev. 2005;14(5):1065–73.
Sukocheva O, Wadham C. Role of sphingolipids in oestrogen signalling in breast cancer cells: an update. J Endocrinol. 2014;220(3):R25–35.
Sukocheva O, Wadham C, Gamble J, Xia P. Sphingosine-1-phosphate receptor 1 transmits estrogens' effects in endothelial cells. Steroids. 2015;104:237–45.
Davison Z, de Blacquiere GE, Westley BR, May FE. Insulin-like growth factor-dependent proliferation and survival of triple-negative breast cancer cells: implications for therapy. Neoplasia. 2011;13(6):504–15.
We thank Professor Anthony Gill for assessment and advice associated with the IHC staining and quantification.
Aleksandra Ochnik was supported by research fellowship funding from Cancer Institute NSW through the Sydney Vital Translational Cancer Research Centre. The funding bodies had no role in the design of the study, collection, analysis, and interpretation of data, or in writing the manuscript.
Availability of data and materials
All data generated or analysed during this study are included in this published article.
Ethics approval and consent to participate
This study was approved by the Human Research Ethics Committee of the Northern Sydney Local Health District (Reference Numbers: RESP/15/125 and LNR/15/HAWKE/182) for the analysis of human breast cancer tissues samples obtained from the Australian Breast Cancer Tissue Bank. All samples obtained from this bank were de-identified and were from donors who had given written informed consent for their banked tumor tissue to be used in future research projects.
Consent for publication
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Full dataset used in survival analyses. (PDF 86 kb)
Immunohistochemistry and manual scoring analysis of Australian Breast Cancer Tissue Bank patient samples. Refer to Fig. 1 for experimental details. (DOCX 332 kb)
Prognostic survival outcomes in relation to IGF1R protein expression in HER2 and PR-positive and negative breast cancer tissues. Kaplan-Meier analysis was performed to measure the overall survival (OS) following stratification for high vs. low IGF1R protein expression, stratified for HER2 and PR expression: A. IGF1R (HER2-negative); B. IGF1R (HER2-positive); C. IGF1R (PR-negative) and D. IGF1R (PR-negative) as described under Methods. (DOCX 264 kb)
Survival outcomes in relation to p-IGF1R protein expression in HER2 and PR-positive and negative breast cancer tissues. Kaplan-Meier analysis was performed to measure the overall survival (OS) following stratification for high vs. low p-IGF1R protein expression, stratified for HER2 and PR expression: A. p-IGF1R (HER2-negative); B. p-IGF1R (HER2-positive); C. p-IGF1R (PR-negative) and D. p-IGF1R (PR-negative) as described under Methods. (DOCX 279 kb)
Breast cancer cell viability dose curves in response to SKI-II. Refer to Fig. 5 for experimental details. A. MCF7, B. T47D, C. HCC-1806 and D. HCC70 were treated with single agents of the SphK1 inhibitor (SKI-II; 0.16, 0.8, 4, 10 and 20 μM) for 96 h. Graphs depict experimental data normalized to zero treatment vehicle control and 1-way ANOVA followed by Tukey’s test was performed to determine significance between treatment groups and significance accepted p-values *p < 0.05, **p < 0.01, ***p < 0.001 and ****p < 0.0001. (DOCX 186 kb)
OSI-906 inhibits p-IGF1R and IGF1R signaling factors in breast cancer cell-lines. 3 × 105 ER-positive A. MCF7 and ER-negative B. HCC-1806 breast cancer cells were plated per well of a 6-well plate and cultured for 24 h and subsequently treated with the dual IGF1R/InsR dual tyrosine kinase inhibitor (OSI-906; 0.1, 0.4 and 1.6 μM) and/or SphK1 inhibitor (SKI-II; 4 μM) for 24 h. Protein lysates were collected and 20 μg of protein was used for immunoblot analysis to measure changes to IGF1R signaling (i.e. AKT and eIF4E-BP1 (protein translation)) and SphK1 steady-state protein expression levels. B-actin was used as a loading control. Abbreviations: Hyper-P = hyper-phosphorylation and Hypo-P = hypophosphorylated. (DOCX 255 kb)
Effects of co-targeting IGF1R and SphK1 on colony formation in breast cancer cell-lines. Refer to Fig. 5 for experimental details. Images represent individual wells of colonies following drug treatments for 10-14 d of ER-positive A. MCF7 and B. T47D and ER-negative C. HCC1806 breast cancer cell-lines. Images were captured using FujiFilm Luminescent Image Analyzer LAS-300 (Stamford, CT) and single colonies were counted using open colony forming unit (CFU) software (http://opencfu.sourceforge.net/). (DOCX 545 kb)
About this article
Cite this article
Ochnik, A.M., Baxter, R.C. Insulin-like growth factor receptor and sphingosine kinase are prognostic and therapeutic targets in breast cancer. BMC Cancer 17, 820 (2017). https://doi.org/10.1186/s12885-017-3809-0