Proteomics profiling identify CAPS as a potential predictive marker of tamoxifen resistance in estrogen receptor positive breast cancer
© Johansson et al.; licensee BioMed Central. 2015
Received: 20 October 2014
Accepted: 28 February 2015
Published: 21 March 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
Clinicopathological characteristics of patients included in the discovery proteomics
Control n =12
Relapse n = 12
Planned tamoxifen regimen
ER (fmol/μg DNA)
PgR (fmol/μg DNA)
Statistical analysis of proteomics data to identify tamoxifen-predictive markers
Connection of the 13 protein panel to cancer
Connection to cancer
Transcription elongation factor A (SII) like 4 (TCEAL4) is down-regulated in anaplastic thyroid cancer.
Stimulates lipid degradation. AZGP1 is a tumor suppressor in pancreatic cancer inducing mesenchymal-to-epithelial (MET) transdifferentiation by inhibiting TGF-β-mediated ERK signaling. The percentage of IHC positive stained cells in (Pancreatic Intraepithelial Neoplasia) PanIN lesions, primary and metastatic PDAC, gradually decreases from 48 to 26 and 5%, respectively.
More malignant breast tumors showed downregulated AZGP1 mRNA.
Low AZGP1 expression, by IHC, was associated with clinical recurrence in prostate cancer.
S100A10 was down-regulated in breast cancer, irrespective of pathological stage.
Complex with AnnexinA2 and AHNAK and E-cadherin at the plasma membrane.
S100A10 is required for recruitment of macrophages to tumor sites and tumor growth.
Catalyzes the irreversible oxidative decarboxylation of malonate and methylmalonate semialdehydes to acetyl- and propionyl-CoA. Decreased expression with increasing grade in kidney cancer.
Associates with S100A10 and E-cadherin at the plasma membrane.
Converts fructose-1,6-bisphosphate to fructose 6-phosphate in gluconeogenesis. Antagonize glycolysis. Downregulated through NF-kappaB pathway in Ras-transformed NIH3T3 cells. Restoration of FBP1 expression suppressed anchorage-independent growth.
Loss of FBP1 by Snail-G9a-Dnmt1 complex increase glucose uptake, glycolysis and induce a cancer stem cell like characteristics.
Involved in regulation of angiogenesis, cell survival, motility, and invasion. EMT marker.
Staining for S100A4 is associated with poorer survival in BC.
High levels of HSP90AB1 correlates with poor prognosis in HER2-/ER+ BC.
Rate-limiting enzyme of hexosamine biosynthetic pathway (HBP).
IHC analysis indicated elevated GlcNAcylation levels in breast tumor tissue as compared to adjacent tissue. GlcNAcylation was significantly enhanced in metastatic lymph nodes compared with their corresponding primary tumor tissues.
BC cells upregulate HBP, including increased O-GlcNAcation and elevated expression of O-GlcNAc transferase (OGT), which is the enzyme catalyzing the addition of O-GlcNAc to proteins.
Required for synthesis of pyridoxal-5-phosphate from vitamin B6. pyridoxal-5-phosphate is a prosthetic group of some enzymes.
Overexpression stimulates cell migration. Small GTPase Rab21 regulates cell adhesion and controls endosomal traffic of β1-integrins.
Rab21 is required for Carcinoma-associated fibroblasts (CAFs) to promote the invasion of cancer cells. RAB21 enables the accumulation of integrin a5 at the plasma membrane and subsequent force-mediated matrix remodelling.
Stable exogenous MX1 expression inhibited in vitro motility and invasiveness of human prostate carcinoma cell line PC-3 M.
Upregulated in a mammary carcinoma xenograft model of tamoxifen resistance.
Upregulated in a fulvestrant-resistant cell line T47D-r on the mRNA and protein level compared to T47D.
High expression gives bad prognosis in endometrial cancer.
Over-expressed in ovarian cancer.
Upregulated 30× in MCF7 ErbB2 compared to MCF7.
Phosphorylation substrate for protein kinase A.
Validation of potential predictive protein biomarkers in tumor homogenates
Clinicopathological characteristics of patients included in the validation cohort
PgR (fmol/μg DNA)
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).
Univariate and Cox proportional multivariate analyses for recurrence (n = 29) and breast cancer death (n = 17) in the 79 patients validation cohort with ER positive breast cancer, treated with tamoxifen only as systemic adjuvant therapy
Tumor size (T1 vs. T2–3)
Node status (pos vs neg)
CAPS (high vs low)
Age (<50 vs ≥50 years)
ER cont (fmol/μg DNA)
PgR cont (fmol/μg DNA)
CAPS (high vs low)
Tumor size (T1 vs T2–3)
Node status (pos vs neg)
CAPS (high vs low)
Age (<50 vs ≥50 years)
ER cont (fmol/μg DNA)
PgR cont (fmol/μg DNA)
CAPS (high vs low)
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.
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.
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Musgrove EA, Sutherland RL. Biological determinants of endocrine resistance in breast cancer. Nat Rev Cancer. 2009;9:631–43.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- Lam HY. Tamoxifen is a calmodulin antagonist in the activation of cAMP phosphodiesterase. Biochem Biophys Res Commun. 1984;118:27–32.View ArticlePubMedGoogle Scholar
- Lopes MC, Vale MG, Carvalho AP. Ca2(+)-dependent binding of tamoxifen to calmodulin isolated from bovine brain. Cancer Res. 1990;50:2753–8.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- Zeisberg M, Neilson EG. Biomarkers for epithelial-mesenchymal transitions. J Clin Invest. 2009;119:1429–37.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- Morley P, Whitfield JF. Effect of tamoxifen on carbachol-triggered intracellular calcium responses in chicken granulosa cells. Cancer Res. 1994;54:69–74.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Yasue A, Hasegawa K, Udagawa Y. Effects of tamoxifen on the endometrium and its mechanism of carcinogenicity. Hum Cell. 2011;24:65–73.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Kalluri R, Weinberg RA. The basics of epithelial-mesenchymal transition. J Clin Invest. 2009;119:1420–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Hanover JA, Krause MW, Love DC. The hexosamine signaling pathway: O-GlcNAc cycling in feast or famine. Biochim Biophys Acta. 1800;2010:80–95.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Perroud B, Ishimaru T, Borowsky AD, Weiss RH. Grade-dependent proteomics characterization of kidney cancer. Mol Cell Proteomics. 2009;8:971–85.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Boye K, Maelandsmo GM. S100A4 and metastasis: a small actor playing many roles. Am J Pathol. 2010;176:528–35.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- Skubitz AP, Pambuccian SE, Argenta PA, Skubitz KM. Differential gene expression identifies subgroups of ovarian carcinoma. Transl Res. 2006;148:223–48.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.