Proteomics profiling identify CAPS as a potential predictive marker of tamoxifen resistance in estrogen receptor positive breast cancer
Clinical Proteomics volume 12, Article number: 8 (2015)
Despite the success of tamoxifen since its introduction, about one-third of patients with estrogen (ER) and/or progesterone receptor (PgR) - positive breast cancer (BC) do not benefit from therapy. Here, we aim to identify molecular mechanisms and protein biomarkers involved in tamoxifen resistance.
Using iTRAQ and Immobilized pH gradient-isoelectric focusing (IPG-IEF) mass spectrometry based proteomics we compared tumors from 12 patients with early relapses (<2 years) and 12 responsive to therapy (relapse-free > 7 years). A panel of 13 proteins (TCEAL4, AZGP1, S100A10, ALDH6A1, AHNAK, FBP1, S100A4, HSP90AB1, PDXK, GFPT1, RAB21, MX1, CAPS) from the 3101 identified proteins, potentially separate relapse from non-relapse BC patients. The proteins in the panel are involved in processes such as calcium (Ca2+) signaling, metabolism, epithelial mesenchymal transition (EMT), metastasis and invasion. Validation of the highest expressed proteins in the relapse group identify high tumor levels of CAPS as predictive of tamoxifen response in a patient cohort receiving tamoxifen as only adjuvant therapy.
This data implicate CAPS in tamoxifen resistance and as a potential predictive marker.
Consensus guidelines for adjuvant breast cancer (BC) therapy advise different treatment modalities to diminish the risk of recurrence after surgery for primary breast cancer . Factors taken into account beside stage of the disease are histopathological parameters, expression of steroid receptors, overexpression or amplification of HER2, and proliferation. Adjuvant endocrine therapy for a minimum of five years postoperatively is advised to patients with BC expressing estrogen and/or progesterone receptors (ER and PgR) and half the recurrence rate in this group. However, about a third of the eligible patients will relapse during or after tamoxifen therapy and even more so patients with advanced BC . A lot of effort has been made to find markers to endocrine therapy, see Musgroove and Sutherland for a review .
Following the work by Sørlie and colleagues that presented the intrinsic molecular subgroups of BC based on gene expression patterns, a substantial amount of information has elucidated the complexity in pathways driving the different BC subgroups. The intrinsic subgroups differ molecularly, in prognosis as well as relapse rates after different therapy modalities . The two luminal subgroups originating from the well differentiated luminal layer are exclusively ER positive and have in general a favorable prognosis. However, there is a substantial heterogeneity within the ER positive groups where luminal B has been characterized by a more aggressive disease course compared to luminal A in terms of recurrence rate which could partly be explained by the proportion of HER2 positive patients in this group . Genomic studies have great impact on BC classification enabling the identification of subgroups within this heterogeneous disease which is consequently taking us a step closer to personalizing therapy.
Proteomics is a mean of complementing the genomics information since mRNA and protein levels don’t always correlate . Several different proteins and signaling pathways have been suggested to be part of the tamoxifen resistance mechanism, for example kinase expression levels and activity, transcription factors and their coregulators, as well as downstream intracellular events as the PIK3/AKT/mTOR pathway [3,6]. In addition, other nuclear receptors, as the retinoic acid receptor alpha are involved in tamoxifen resistance .
Apart from binding and inhibiting ER, tamoxifen also bind and inhibit the calcium binding protein calmodulin (CALM) [8,9]. CALM regulates many cellular protein kinases, phosphatases and transmembrane ion transporters, mainly in a calcium dependent manner. CALM interacts and modulate ER activity . Another member of the EF hand motif family is Calcyphosine (CAPS) - for calcium binding and regulated by cyclic AMP through phosphorylation protein. CAPS has been suggested as an alternative calcium signaling route to CALM . CAPS have high levels in endometrial tumors, whose proliferation is known to be induced by tamoxifen, compared to normal proliferative tissue .
Here we use mass spectrometry (MS)-based proteomics to discover potential predictive biomarkers for adjuvant tamoxifen therapy in a patient population that received adjuvant tamoxifen as the only systemic adjuvant therapy. We identified 13 proteins showing significant differential expression in relapsing patients compared to our defined control group. These were involved in processes such as calcium (Ca2+) signaling, metabolism, epithelial mesenchymal transition, metastasis and invasion [12-14]. Validation of calcyphosine (CAPS) in the entire clinical cohort suggests that high levels predict relapse, and that CAPS is a potential predictor of tamoxifen response.
Experimental design and mass spectrometry based proteomic
To do an unbiased search for tamoxifen predictive markers, we performed quantitative proteomics on tumor homogenates from BC patients. We selected tumors from 12 patients who relapsed within 2 years of tamoxifen treatment (referred to as relapse) and 12 patients with a disease-free follow up time of more than 7 years (referred to as control). Patients were matched into 12 pairs defined by age, tumor size, and node status. All patients were ductal and ER positive cancers (Table 1). Quantitative mass spectrometry based proteomics on these patient samples was performed by nanoLC-MS/MS using LTQ-Orbitrap Velos mass spectrometer on fractions from peptide isoelectric focusing, pH 3.4–4.8. See Figure 1A for workflow. MS based proteomics yielded a total of 3101 identified proteins, of which 550 overlapped between all the 4 iTRAQ sets and used in the analysis (Additional file 1).
Statistical analysis of proteomics data to identify tamoxifen-predictive markers
Uni- and multivariate analysis by SAM and OPLS were used to identify potential tamoxifen predictive markers [15,16]. These statistical analyses revealed a 13 protein signature, which could separate relapse vs. control groups, P-value 2.2e-005 (Figure 1B). Principal component analysis (PCA) displayed no bias between and within the relapse and control groups (Additional file 2A-C). Protein identities and iTRAQ protein quantities of the potential 13 protein signature are shown in Figure 1C. Many of the proteins in the 13 protein panel have been connected to BC and other cancer types before (Table 2). There is connection to EMT via S100A4, an EMT marker, as well as S100A10 and AHNAK, who are protein complex components together with E-cadherin at the plasma membrane.
Validation of potential predictive protein biomarkers in tumor homogenates
A positive selection marker is generally preferred over a negative marker and since our aim was to identify patients relapsing on tamoxifen, we choose to verify expression of the 2 proteins with the highest expression ratio in the relapse group, CAPS and MX1, from the proteomics data (Figure 1C). An initial verification was performed by western blot (WB) on four randomly selected matched pairs of cytosols from all 24 patients included in the study. The WB showed that both CAPS and MX1 had overall higher protein levels in the relapse group compared to control (Figure 2A). Correlation between MS and WB data were 0.8 R 2 (p = 0.0037) for CAPS and 0.6 R 2 (p = 0.024) for MX1 (Additional file 3).
Based on this small verification we performed protein quantification by ELISA for CAPS and MX1 on 79 and 89 breast tumor homogenates respectively. See Table 3 for clinical characteristics. Relapse-free survival (RFS) and breast cancer-specific survival (BCS) were used as primary end points, based on the time from diagnosis to the first event of loco-regional or distant recurrence, and time from diagnosis to breast cancer death, respectively.
High protein levels of CAPS was associated with increased RFS (p = 0.049, 85% power) and BCS (p = 0.11, 80% power) (Figure 2B, C). Selection of lower tertile vs. the 2 higher tertiles for CAPS was done based on difference in HR to continuous data (Additional file 4). MX1 protein levels were not correlated to RFS or BCS in this clinical cohort (Additional file 5).
CAPS remained an independent prognostic marker in a Cox proportional multivariate analysis for RFS (HR = 3.6; p = 0.011), while nodal status (node-negative versus node-positive) (HR = 1.7; p = 0.22), tumor size (<20 mm versus ≥ 20 mm) (HR = 1.4; p = 0.43), age (HR = 1.0; p = 0.87), ER (HR = 0.88; p = 0.073) and PgR (HR = 0.96; p = 0.48) were not (Table 4). CAPS was also statistical significant factor in multivariate analysis for BCS (HR = 4.0 p = 0.043). Despite a total of 29 relapses and 17 breast cancer deaths were registered among the 79 patients with data on CAPS.
Receiver operator characteristic (ROC) for CAPS yielded an area of 0.62 (p = 0.08) and 0.58 (p = 0.29) for RFS and BCS, respectively (Additional file 6).
The vast majority of BC patients are ER positive, making them eligible for adjuvant endocrine treatment. Tamoxifen has been the corner stone in breast cancer treatment since 40 years, although in part replaced with the aromatase inhibitors (AI) during the last decade. However, due to the marginal benefit and the severe side effects accomplished with aromatase inhibitor therapy, tamoxifen still plays an important role in breast cancer treatment. Tamoxifen alone is not sufficient in about 30% of these patients .
Genomics together with proteomics are on their way to decipher the complex molecular events that are responsible for the heterogeneity within subgroups of this disease . In this study we compared the protein differences between patients with a relapse-free survival of more than 7 years and those exhibiting a relapse within 2 years of tamoxifen treatment. We identified a potential signature of 13 proteins that was able to significantly differentiate relapsing patients from our defined control group (Figure 1) Expression analysis of the 2 proteins with highest levels in the relapse group confirmed that high levels of CAPS, but not MX1, to be related to tamoxifen non-responsive patients.
High levels of CAPS were statistically significantly correlated to lower RFS both in uni- and multivariate analyses, including important breast cancer prognostic markers as nodal status, tumor size, age and ER and PgR levels. Several gene array based tests are validated and commercially available for breast cancer care. Some are developed to gain information about endocrine therapy. The Oncotype dx separates patients where endocrine therapy solely is enough versus patients that are in need of additional chemotherapy . The majority of markers included in these profiles have been linked to proliferation. Endopredict can identify ER+, HER2 negative patients with an increased likelihood of development of late distant metastasis . Interestingly, Endopredict include CALM2, a Ca2+ binding protein. Our ROC analysis also indicates that CAPS protein levels needs to be used together with other markers to obtain good sensitivity and specificity to predict outcome. However, we can not exclude CAPS as a prognostic factor with the clinical cohort used in this study.
CAPS is a Ca2+ binding protein whose synthesis and activation is induced by the cAMP cascade (PKA) and connected to signaling for cellular proliferation and differentiation . Interestingly, tamoxifen also has effects on Ca2+ signaling, aside from its antiestrogenic properties by binding to and inhibiting CALM activity . Our data of CAPS overexpression in non-responsive patients raises the hypothesis that tamoxifen fails to inhibit the calcium signaling response in the cell, and that high CAPS expression can function as an alternative compensatory mechanism when CALM is inhibited. CAPS has previously been suggested to be an alternative signaling step to CALM . Moreover, recent studies found higher CAPS levels in endometrial cancer compared to normal proliferative tissue [12,21]. It is of note that tamoxifen has been shown to increase the risk for endometrial cancer due to proliferative effects of the drug in this tissue .
The MX1 results are surprising since MX1 showed increased mRNA and protein levels in tamoxifen and fulvestrant resistance models [23,24]. The number of psms for quantification differs between the 4 iTRAQ sets, with 1, 3, 8, and 27 psms used for iTRAQ quantification, which could have introduced quantification bias. However, the negative result for MX1 with ELISA is probably due to antibody effects, i.e. different or mixed epitope recognition by the antibody compared to proteomics data. In comparison, CAPS had 28, 5, 21 and 14 psms used for quantification in the 4 different iTRAQ 8plex sets. Hence, CAPS had a more robust MS quantification that is less sensitive to noise and more likely to identify the major protein variant that will ease ELISA validation.
Previous studies using DNA microarrays and proteomics have shown upregulation of MX1 in vivo and in vitro in cells resistant to tamoxifen (in vivo) and to the ER downregulator fulvestrant (in vitro) [23,24]. A recent in vitro study implicated the proliferative PIK3/AKT pathway as part of the regulatory cascade inducing MX1 expression in response to IFNα . The upregulation of MX1 in the relapsing patients may be a consequence of induction of growth signaling through various pathways.
Some draw backs with the present study must be discussed. So far, proteomics have been very labor intensive making discovery work in a large number of tumors difficult. The small sample set for the discovery phase potentially increase putative markers without any value. Larger test sets of at least 100 patients in each group are probably needed for identification of more robust markers. Another drawback is that lobular patients were excluded. These patients are often ER positive and subjected to adjuvant endocrine treatment, thus the role of CAPS has to be studied in all ER positive breast cancer types. However, an advantage of proteomics is the direct identification of a particular protein, or a set of proteins, which are measurable by more time-effective, simpler and less costly routine analysis methods such as IHC and ELISA in the clinic.
The potential 13 protein panel shows connection to breast and other forms and cancer (Table 2). Among the proteins, there is connection to invasion, motility and epithelial mesenchymal transition (EMT) by RAB21, MX1, S100A4, S100A10, AHNAK, AZGP1. S100A10, AHNAK and E-cadherin form protein complexes on the plasma membrane [26,27]. S100A4 (also known as FSP1 and MTS-1) is an EMT marker . There is also a connection to metabolic changes by GFPT1, PDXK, FBP1, ALDH6A1 and AZGP1. GFPT1 is the rate limiting enzyme into the hexosamine synthesis pathway and increased downstream GlcNAcylation is increased in breast tissue compared to normal tissue [29,30]. PDXK catalyze the synthesis of the prosthetic group pyridoxal-5-phosphate from vitamin B6. The Snail-G9a-DNMT1 protein complex has been shown not only to decrease E-cadherin levels but also FBP1, increasing glucose uptake, inducing glycolysis and cancer stem cell like characteristics .
We have performed MS based proteomics on 12 tumors from relapse and 12 tumors from control patients and associated a 13 protein panel with tamoxifen resistance. Validation of the highest expressed protein in the relapse group, CAPS, in the whole cohort implicate CAPS as a potential predictive factor in ER positive breast cancer receiving adjuvant tamoxifen.
Patient characteristics and selection
The study investigate differences in protein expression profile between patients exhibiting relapses within two years of tamoxifen therapy (referred to as relapse group) and patients with a relapse-free survival of more than seven years (referred to as control group) included twenty-four patients (12 patients in each group). The reason for choosing a relapse-free period of more than 7 years was to include the possible hang-over effect by tamoxifen. In order to avoid bias from biological differences, the two groups were matched for tumor size, nodal status, and age. Only ductal cancers were included. Data on histopathological grade was not available. These patients originate from a BC material previously described in detail . In short, a total of 402 patients with ER positive BC subjected to adjuvant tamoxifen as the only adjuvant systemic therapy (no chemotherapy allowed) were identified from a larger consecutive cohort of 711 patients diagnosed with a primary breast stage I-III from 1991 to 1996. Patients with a locally advanced BC, displaying distant metastases at diagnosis, or having received neoadjuvant therapy, were excluded from the cohort of 402 patients. The median age of patients included in the present study was 61.6 years for the control group and 65.1 for the relapse group. A description of these patients’ characteristics is shown in Table 1. Being unable to retrieve S-phase data for all 24 patients, this factor was excluded. The study was approved by the research ethical boards of Linköping University Hospital, Linköping, Sweden. During the study period, there was no requirement to receive an informed consent from each patient for storage and usage of the remaining tumor homogenates according to Swedish law.
Preparation of tumor homogenates
Representative tumor tissue was homogenized in a microdismembrator (Braun, Melsungen, Germany) and suspended in cold potassium phosphate buffer (5 mM, pH 7.4, 10% glycerol v/v, 1 mM dithiothreitol). Supernatants containing the cytosolic fractions were collected after refrigerated centrifugation at 20,000 g, used for steroid receptor content analysis and stored at –70°C. The pellet fractions were analyzed by the method of Burton, in order to evaluate DNA concentration.
Sample preparation for mass spectrometry
Tumor homogenates (cytosols) were acetone precipitated and pellets were dissolved in 1% SDS, which were diluted to 0.1% before protein concentration determination with Bio-rad DC protein assay. Triethylammoniumbicarbonat was added to each sample to give a final concentration of 0.5 M. Proteins were reduced by adding tris-(2-carboxyethyl) phosphine (TCEP) and alkylated by methyl methanethiosulfonate (MMTS). Trypsin (modified sequence grade, Promega, Madision WI, USA) was added (1:20, trypsin:protein) and the samples were incubated overnight at 37°C. iTRAQ labelling of the peptides were done according to the manufacturer’s protocol (Applied Biosystems) and cleaned by a strata-X-C-cartridge (Phenomenex). Four different iTRAQ 8plex sets was used to include all patients, and the 4 sets were connected by 2 iTRAQ channels with pooled samples, one for the 12 controls (PC) and one for the 12 relapse (PR) samples. Three pairs of matched control (115, 116, 117) and relapse (118, 119, 121) samples were included together with the pooled samples (113, 114) in each iTRAQ 8plex set.
The iTRAQ labelled peptides, 360 μg for each of the 4 iTRAQ sets, were separated by immobilized pH gradient - isoelectric focusing (IPG-IEF) on a narrow range pH 3.4–4.8 strip as described by Branca et al. . Peptides were extracted from the strips by a prototype liquid handling robot, kindly supplied by GE Healthcare Bio-Sciences AB. A plastic device with 72 wells was put onto each strip and 50 μl of MQ was added to each well. After 30 minutes incubation, the liquid was transferred to a 96 well plate and the extraction was repeated 2 more times. The extracted peptides were dried in speedvac and dissolved in 3% acetronitrile (ACN), 0.1% formic acid.
Before analysis on the LTQ-Orbitrap Velos (Thermo Fischer Scientific, San Jose, CA, USA), peptides were separated using an Agilent 1200 nano-LC system. Samples were trapped on a Zorbax 300SB-C18, and separated on a NTCC-360/100-5-153 (Nikkyo Technos., Ltd) column using a gradient of A (3% ACN, 0.1% FA) and B (95% ACN, 0.1% FA), ranging from 3% to 40% B in 45 min with a flow of 0.4 μl/min. The LTQ-Orbitrap Velos was operated in a data-dependent manner, selecting 5 precursors for sequential fragmentation by CID and HCD, and analyzed by the linear iontrap and orbitrap, respectively. The survey scan was performed in the Orbitrap at 30.000 resolution (profile mode) from 300–2000 m/z, using lock mass at m/z 445.120025, with a max injection time of 500 ms and AGC set to 1 × 106 ions. For generation of HCD fragmentation spectra, a max ion injection time of 500 ms and AGC of 2 × 105 were used before fragmentation at 50% normalized collision energy. For FTMS MS2 spectra, normal mass range was used, centroiding the data at 7500 resolution. Peptides for CID were accumulated for a max ion injection time of 200 ms and AGC of 3 × 104, fragmented with 35% collision energy, wideband activation on, activation q 0.25, activation time 10 ms before analysis at normal scan rate and mass range in the linear iontrap. Precursors were isolated with a width of 2 m/z and put on the exclusion list for 60 s. Single and unassigned charge states were rejected from precursor selection.
Peptide and protein identification
All Orbitrap data was searched by Mascot 2.2 (Matrix Science Limited, London, UK) under the software platform Proteome Discoverer 1.1 (Thermo) against the human swissprot database (build 57.13) and results were limited to a false discovery rate of <1%. Precursor mass tolerance was set to 10 ppm, and product mass tolerances of 0.015 Da for HCD-FTMS and 0.7 Da for CID-ITMS were used. Oxidized methionine was set as dynamic modification and methylthio, N-terminal 8plex iTRAQ, and lysyl 8plex iTRAQ as fixed modifications. Quantification of iTRAQ 8plex reporter ions was done by Proteome Discoverer on HCD-FTMS tandem mass spectra using an integration window tolerance of 20 ppm. Proteins with 1 or more of 1% FDR confident peptides were used in the following data analysis.
Identification of potential biomarkers
Potential biomarkers were selected by multivariate analysis of the BC patient samples with orthogonal partial least square (OPLS) analysis using Simca software, using the S-plot and loading plot described by Wiklund et al. 2008 and univariate analysis corrected for multiple testing using SAM [15,16].
Western blot and ELISA
The randomly selected 8 patients included in the verification represent 4 control and 4 relapse patients from all 24 patients included in the study. 30 μg of tumor homogenate per sample were run on SDS-PAGE gels (NuPage Bis Tris 4–12%, Invitrogen), blotted and subsequently incubated overnight with CAPS (sc-134298), (Santa Cruz Biotech) and MX1 (HPA030917) (Sigma-Aldrich). After incubation with anti-rabbit or anti-mouse IgG-HRP Conjugate (BioRad), bands were detected using ECL™ Western Blotting Detection Reagents (Amersham). Commercially available enzyme-linked immune sorbent assays (ELISA) for quantification of CAPS (E92359Hu) and MX1 (E80763Hu) were performed according to the protocols of the manufacturer (USCN, Cologne, Germany).
Statistical analysis of validation cohort
Differences in CAPS levels in relation to patient’s age and tumor characteristics were analyzed with Mann-Whitney’s test. Relapse-free survival (RFS) and breast cancer-specific survival (BCS) were chosen as primary end points, based on the time from diagnosis to the first event of loco-regional or distant recurrence, and time from diagnosis to breast cancer death, respectively. Survival curves and probabilities of RFS and BCS were estimated using the Kaplan-Meier method. Hazard ratios were calculated using Cox hazard regression analysis. Multivariate Cox models included the variables age, nodal status, tumor size, ER and PgR levels and CAPS levels. The software STATISTICA 10 (StatSoft Inc., Tulsa, OK, USA) was used for statistical calculations. P-values less than 0.05 in two-sided tests were considered significant. Graphpad prism 6 was used to visualize survival curves and perform ROC and correlation analysis. Power calculations were done using Graphpad StateMate 2.00.
Goldhirsch A, Wood WC, Gelber RD, Coates AS, Thurlimann B, Senn HJ. Progress and promise: highlights of the international expert consensus on the primary therapy of early breast cancer 2007. Ann Oncol. 2007;18:1133–44.
Early Breast Cancer Trialists’ Collaborative Group (EBCTCG). Effects of chemotherapy and hormonal therapy for early breast cancer on recurrence and 15-year survival: an overview of the randomised trials. Lancet. 2005;365:1687–717.
Musgrove EA, Sutherland RL. Biological determinants of endocrine resistance in breast cancer. Nat Rev Cancer. 2009;9:631–43.
Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, et al. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci U S A. 2001;98:10869–74.
Schwanhausser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, et al. Global quantification of mammalian gene expression control. Nature. 2011;473:337–42.
Grunt TW, Mariani GL. Novel approaches for molecular targeted therapy of breast cancer: interfering with PI3K/AKT/mTOR signaling. Curr Cancer Drug Targets. 2013;13:188–204.
Johansson HJ, Sanchez BC, Mundt F, Forshed J, Kovacs A, Panizza E, et al. Retinoic acid receptor alpha is associated with tamoxifen resistance in breast cancer. Nat Commun. 2013;4:2175.
Lam HY. Tamoxifen is a calmodulin antagonist in the activation of cAMP phosphodiesterase. Biochem Biophys Res Commun. 1984;118:27–32.
Lopes MC, Vale MG, Carvalho AP. Ca2(+)-dependent binding of tamoxifen to calmodulin isolated from bovine brain. Cancer Res. 1990;50:2753–8.
Garcia Pedrero JM, Del Rio B, Martinez-Campa C, Muramatsu M, Lazo PS, Ramos S. Calmodulin is a selective modulator of estrogen receptors. Mol Endocrinol. 2002;16:947–60.
Lecocq R, Lamy F, Erneux C, Dumont JE. Rapid purification and identification of calcyphosine, a Ca(2+)-binding protein phosphorylated by protein kinase A. Biochem J. 1995;306(Pt 1):147–51.
Voisin SN, Krakovska O, Matta A, DeSouza LV, Romaschin AD, Colgan TJ, et al. Identification of novel molecular targets for endometrial cancer using a drill-down LC-MS/MS approach with iTRAQ. PLoS One. 2011;6:e16352.
Zeisberg M, Neilson EG. Biomarkers for epithelial-mesenchymal transitions. J Clin Invest. 2009;119:1429–37.
Oxelmark E, Roth JM, Brooks PC, Braunstein SE, Schneider RJ, Garabedian MJ. The cochaperone p23 differentially regulates estrogen receptor target genes and promotes tumor cell adhesion and invasion. Mol Cell Biol. 2006;26:5205–13.
Wiklund S, Johansson E, Sjostrom L, Mellerowicz EJ, Edlund U, Shockcor JP, et al. Visualization of GC/TOF-MS-based metabolomics data for identification of biochemically interesting compounds using OPLS class models. Anal Chem. 2008;80:115–22.
Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98:5116–21.
Galvao ER, Martins LM, Ibiapina JO, Andrade HM, Monte SJ. Breast cancer proteomics: a review for clinicians. J Cancer Res Clin Oncol. 2011;137:915–25.
Paik S, Tang G, Shak S, Kim C, Baker J, Kim W, et al. Gene expression and benefit of chemotherapy in women with node-negative, estrogen receptor-positive breast cancer. J Clin Oncol. 2006;24:3726–34.
Dubsky P, Brase JC, Jakesz R, Rudas M, Singer CF, Greil R, et al. The EndoPredict score provides prognostic information on late distant metastases in ER+/HER2- breast cancer patients. Br J Cancer. 2013;109:2959–64.
Morley P, Whitfield JF. Effect of tamoxifen on carbachol-triggered intracellular calcium responses in chicken granulosa cells. Cancer Res. 1994;54:69–74.
Li Z, Huang C, Bai S, Pan X, Zhou R, Wei Y, et al. Prognostic evaluation of epidermal fatty acid-binding protein and calcyphosine, two proteins implicated in endometrial cancer using a proteomic approach. Int J Cancer. 2008;123:2377–83.
Yasue A, Hasegawa K, Udagawa Y. Effects of tamoxifen on the endometrium and its mechanism of carcinogenicity. Hum Cell. 2011;24:65–73.
Becker M, Sommer A, Kratzschmar JR, Seidel H, Pohlenz HD, Fichtner I. Distinct gene expression patterns in a tamoxifen-sensitive human mammary carcinoma xenograft and its tamoxifen-resistant subline MaCa 3366/TAM. Mol Cancer Ther. 2005;4:151–68.
Huber M, Bahr I, Kratzschmar JR, Becker A, Muller EC, Donner P, et al. Comparison of proteomic and genomic analyses of the human breast cancer cell line T47D and the antiestrogen-resistant derivative T47D-r. Mol Cell Proteomics. 2004;3:43–55.
Chai Y, Huang HL, Hu DJ, Luo X, Tao QS, Zhang XL, et al. IL-29 and IFN-alpha regulate the expression of MxA, 2′,5′-OAS and PKR genes in association with the activation of Raf-MEK-ERK and PI3K-AKT signal pathways in HepG2.2.15 cells. Mol Biol Rep. 2011;38:139–43.
Benaud C, Gentil BJ, Assard N, Court M, Garin J, Delphin C, et al. AHNAK interaction with the annexin 2/S100A10 complex regulates cell membrane cytoarchitecture. J Cell Biol. 2004;164:133–44.
Yamada A, Irie K, Hirota T, Ooshio T, Fukuhara A, Takai Y. Involvement of the annexin II-S100A10 complex in the formation of E-cadherin-based adherens junctions in Madin-Darby canine kidney cells. J Biol Chem. 2005;280:6016–27.
Kalluri R, Weinberg RA. The basics of epithelial-mesenchymal transition. J Clin Invest. 2009;119:1420–8.
Hanover JA, Krause MW, Love DC. The hexosamine signaling pathway: O-GlcNAc cycling in feast or famine. Biochim Biophys Acta. 1800;2010:80–95.
Gu Y, Mi W, Ge Y, Liu H, Fan Q, Han C, et al. GlcNAcylation plays an essential role in breast cancer metastasis. Cancer Res. 2010;70:6344–51.
Dong C, Yuan T, Wu Y, Wang Y, Fan TW, Miriyala S, et al. Loss of FBP1 by Snail-mediated repression provides metabolic advantages in basal-like breast cancer. Cancer Cell. 2013;23:316–31.
Sanchez BC, Sundqvist M, Fohlin H, Spyratos F, Nordenskjold B, Stal O, et al. Prolonged tamoxifen treatment increases relapse-free survival for patients with primary breast cancer expressing high levels of VEGF. Eur J Cancer. 2010;46:1580–7.
Branca RM, Orre LM, Johansson HJ, Granholm V, Huss M, Perez-Bercoff A, et al. HiRIEF LC-MS enables deep proteome coverage and unbiased proteogenomics. Nat Methods. 2014;11:59–62.
Akaishi J, Onda M, Okamoto J, Miyamoto S, Nagahama M, Ito K, et al. Down-regulation of transcription elogation factor A (SII) like 4 (TCEAL4) in anaplastic thyroid cancer. BMC Cancer. 2006;6:260.
Kong B, Michalski CW, Hong X, Valkovskaya N, Rieder S, Abiatari I, et al. AZGP1 is a tumor suppressor in pancreatic cancer inducing mesenchymal-to-epithelial transdifferentiation by inhibiting TGF-beta-mediated ERK signaling. Oncogene. 2010;29:5146–58.
Parris TZ, Danielsson A, Nemes S, Kovacs A, Delle U, Fallenius G, et al. Clinical implications of gene dosage and gene expression patterns in diploid breast carcinoma. Clin Cancer Res. 2010;16:3860–74.
Henshall SM, Horvath LG, Quinn DI, Eggleton SA, Grygiel JJ, Stricker PD, et al. Zinc-alpha2-glycoprotein expression as a predictor of metastatic prostate cancer following radical prostatectomy. J Natl Cancer Inst. 2006;98:1420–4.
Carlsson H, Petersson S, Enerback C. Cluster analysis of S100 gene expression and genes correlating to psoriasin (S100A7) expression at different stages of breast cancer development. Int J Oncol. 2005;27:1473–81.
Phipps KD, Surette AP, O’Connell PA, Waisman DM. Plasminogen receptor S100A10 is essential for the migration of tumor-promoting macrophages into tumor sites. Cancer Res. 2011;71:6676–83.
Perroud B, Ishimaru T, Borowsky AD, Weiss RH. Grade-dependent proteomics characterization of kidney cancer. Mol Cell Proteomics. 2009;8:971–85.
Liu X, Wang X, Zhang J, Lam EK, Shin VY, Cheng AS, et al. Warburg effect revisited: an epigenetic link between glycolysis and gastric carcinogenesis. Oncogene. 2010;29:442–50.
Boye K, Maelandsmo GM. S100A4 and metastasis: a small actor playing many roles. Am J Pathol. 2010;176:528–35.
Helfman DM, Kim EJ, Lukanidin E, Grigorian M. The metastasis associated protein S100A4: role in tumour progression and metastasis. Br J Cancer. 2005;92:1955–8.
Rudland PS, Platt-Higgins A, Renshaw C, West CR, Winstanley JH, Robertson L, et al. Prognostic significance of the metastasis-inducing protein S100A4 (p9Ka) in human breast cancer. Cancer Res. 2000;60:1595–603.
Cheng Q, Chang JT, Geradts J, Neckers LM, Haystead T, Spector NL, et al. Amplification and high-level expression of heat shock protein 90 marks aggressive phenotypes of human epidermal growth factor receptor 2 negative breast cancer. BCR. 2012;14:R62.
Caldwell SA, Jackson SR, Shahriari KS, Lynch TP, Sethi G, Walker S, et al. Nutrient sensor O-GlcNAc transferase regulates breast cancer tumorigenesis through targeting of the oncogenic transcription factor FoxM1. Oncogene. 2010;29:2831–42.
Pellinen T, Arjonen A, Vuoriluoto K, Kallio K, Fransen JA, Ivaska J. Small GTPase Rab21 regulates cell adhesion and controls endosomal traffic of beta1-integrins. J Cell Biol. 2006;173:767–80.
Hooper S, Gaggioli C, Sahai E. A chemical biology screen reveals a role for Rab21-mediated control of actomyosin contractility in fibroblast-driven cancer invasion. Br J Cancer. 2010;102:392–402.
Mushinski JF, Nguyen P, Stevens LM, Khanna C, Lee S, Chung EJ, et al. Inhibition of tumor cell motility by the interferon-inducible GTPase MxA. J Biol Chem. 2009;284:15206–14.
Skubitz AP, Pambuccian SE, Argenta PA, Skubitz KM. Differential gene expression identifies subgroups of ovarian carcinoma. Transl Res. 2006;148:223–48.
Wang D, Jensen RH, Williams KE, Pallavicini MG. Differential protein expression in MCF7 breast cancer cells transfected with ErbB2, neomycin resistance and luciferase plus yellow fluorescent protein. Proteomics. 2004;4:2175–83.
Wang S, Chen JZ, Zhang Z, Huang Q, Gu S, Ying K, et al. Cloning, characterization, and expression of calcyphosine 2, a novel human gene encoding an EF-hand Ca(2+)-binding protein. Biochem Biophys Res Commun. 2002;291:414–20.
Funding from the Swedish Research Council, Swedish Cancer Society, Swedish State under the LUA-agreement (Sahlgrenska University Hospital, Gothenburg), the King Gustav V Jubilee Clinic Cancer Research Foundation, Gothenburg.
The authors declare that they have no competing interests.
HJJ, BCS, JL, BKL designed the study and wrote the manuscript. HF, OS, BKL provided clinical material and selected tumor samples. HJJ performed proteomics experiments, HJJ, JF performed statistical analysis of proteomics data, BCS performed western blot and ELISA, OS, HJJ performed statistical analysis of ELISA data. RL, JB, PH provided financial support and input on experimental design. All authors read and approved the final manuscript.
Proteomics data. Protein identifications in the discovery cohort.
Principal component analysis (PCA) of quantitative data from overlapping protein identifications between the iTRAQ sets. (A) Both control and relapse patients, denoted C and R, respectively. (B) Control patients (C) Relapse patients.
Correlation between MS and WB data. Correlation between MS and WB data for (A) CAPS and (B) MX1.
Continuous vs. tertile analysis of CAPS. Table showing continuous vs. tertile analysis of CAPS ELISA data in validation cohort.
Survial analysis for MX1. Relapse free survival (A) and breast cancer survival (B) for MX1. MX1 was not associated with relapse free or overall survival.
ROC analysis for CAPS. ROC analysis using CAPS measurements by ELISA for (A) relapse and (B) overall survival.
About this article
Cite this article
Johansson, H.J., Sanchez, B.C., Forshed, J. et al. Proteomics profiling identify CAPS as a potential predictive marker of tamoxifen resistance in estrogen receptor positive breast cancer. Clin Proteom 12, 8 (2015). https://doi.org/10.1186/s12014-015-9080-y