- Research article
- Open Access
Dynamic changes in gene expression in vivo predict prognosis of tamoxifen-treated patients with breast cancer
Breast Cancer Research volume 12, Article number: R39 (2010)
Tamoxifen is the most widely prescribed anti-estrogen treatment for patients with estrogen receptor (ER)-positive breast cancer. However, there is still a need for biomarkers that reliably predict endocrine sensitivity in breast cancers and these may well be expressed in a dynamic manner.
In this study we assessed gene expression changes at multiple time points (days 1, 2, 4, 7, 14) after tamoxifen treatment in the ER-positive ZR-75-1 xenograft model that displays significant changes in apoptosis, proliferation and angiogenesis within 2 days of therapy.
Hierarchical clustering identified six time-related gene expression patterns, which separated into three groups: two with early/transient responses, two with continuous/late responses and two with variable response patterns. The early/transient response represented reductions in many genes that are involved in cell cycle and proliferation (e.g. BUB1B, CCNA2, CDKN3, MKI67, UBE2C), whereas the continuous/late changed genes represented the more classical estrogen response genes (e.g. TFF1, TFF3, IGFBP5). Genes and the proteins they encode were confirmed to have similar temporal patterns of expression in vitro and in vivo and correlated with reduction in tumour volume in primary breast cancer. The profiles of genes that were most differentially expressed on days 2, 4 and 7 following treatment were able to predict prognosis, whereas those most changed on days 1 and 14 were not, in four tamoxifen treated datasets representing a total of 404 patients.
Both early/transient/proliferation response genes and continuous/late/estrogen-response genes are able to predict prognosis of primary breast tumours in a dynamic manner. Temporal expression of therapy-response genes is clearly an important factor in characterising the response to endocrine therapy in breast tumours which has significant implications for the timing of biopsies in neoadjuvant biomarker studies.
The majority of human breast cancers express estrogen receptor alpha (ERα) and are estrogen responsive . Tamoxifen is still the most widely prescribed anti-estrogen for patients with ER-positive breast cancer and has improved survival in women initially receiving this drug as adjuvant therapy . However, although the majority of women respond to this agent, not all patients benefit and there is a need to identify with greater precision which tumors are sensitive and responding to this therapy. Dynamic changes in specific marker genes in biopsy material at early treatment points could be informative and might indicate whether a tumor is likely to regress or progress.
Although many in vitro studies have explored estrogen- and tamoxifen-regulated changes on gene expression [3–7], we are unaware of any xenograft studies that have investigated the temporal regulation of expression changes produced by tamoxifen in an ER-positive model in vivo. Previous attempts to characterize the gene expression response to tamoxifen in breast tumors in vivo have been limited to single time points [8, 9]. A recent time course experiment demonstrated dynamic gene expression changes in response to estradiol (E2) in ZR-75-1 cell lines in vitro . Xenograft models allow assessment of dynamic changes in tissue gene expression at multiple time points from tissue, which is not feasible in the clinical setting. Furthermore, an in vivo model allows the effect of stromal elements and matrix elements to contribute to expression, which cannot be easily reproduced in vitro.
A number of studies have investigated whether differences in gene expression in primary tumors (prior to treatment) are associated with or can predict the response to tamoxifen [11–13]. Vendrell and colleagues recently described a candidate molecular signature associated with tamoxifen failure in primary breast cancer by examining gene expression in tumors following tamoxifen treatment . An alternative to measuring gene expression differences in the primary static situation is to compare tumor biopsies matched before and after treatment in so-called neoadjuvant 'window of opportunity studies' ; these are likely to generate interesting results in the near future.
We have previously used the ER-positive ZR-75-1 breast cancer xenograft model to demonstrate that tamoxifen causes significant changes in apoptosis, proliferation and angiogenesis within two days of initiating therapy, which both antedated any evidence of growth response and persisted for up to 14 days [3, 16]. Here we present the first study to look at the dynamic changes in gene expression using multiple time points following treatment with tamoxifen in vivo in order to better understand the temporal response to therapy.
Materials and methods
The ZR-75-1, MCF-7 and MDA-MB-231 breast cancer cell lines were obtained from the American Type Culture Collection. Cells were maintained in DMEM (Life Technologies, Paisley, Scotland) containing 10% heat-inactivated FCS, penicillin (100 units/mL), and streptomycin (100 μg/mL). Cells were maintained routinely at 37°C in a humidified atmosphere of 5% carbon dioxide in air. Forty-eight hours before treatment, medium was changed to phenol-red-free DMEM containing 5% double charcoal stripped FCS, glutamine (2 mM), penicillin (100 U/ml) and streptomycin (100 μg/ml). For temporal analysis of gene expression, ZR-75-1, MCF-7 and MDA-MB-231 were incubated with 0.1 nM 17 β-E2 and/or tamoxifen (1 μM) (Sigma-Aldrich Chemical Co St. Louis, MO, USA) or in serum-free medium alone for 0, 6 and 24 hours.
All mouse experiments were performed in accordance with Home Office guidelines. For the xenograft studies, the ZR-75-1 cell line was first implanted into female nu/nu mice. Animals received a subcutaneous slow-release E2 pellet (0.72 mg released over 60 days, Innovative Research of America, Sarasota, FL, USA) on the day of tumor implant. The tumor was then maintained subcutaneously in the flanks of recipient animals by passaging 1 mm3 fragments as required, approximately every eight weeks. For microarray experiments, ZR-75-1 fragments were implanted subcutaneously into animals and allowed to grow to a mean size of 0.25 cm3. All animals received E2. On day 0, animals were randomly allocated to tamoxifen (2.5 mg released over 60 days, Innovative Research of America, Sarasota, FL, USA) or E2-only control groups. There were 20 mice, with tumors in each flank, in both the control and treatment groups of this experiment. Tumor volumes were measured using vernier callipers. Bidimensional tumor diameters were recorded and volumes calculated as vol = πDd2/6, where D is the larger of the two diameters. Four mice from each group were sacrificed at each time point.
Tumor xenografts treated with E2-only or E2 and tamoxifen were obtained from animals sacrificed on days 0, 1, 2, 4, 7 and 14. These were homogenized in lysis buffer and total RNA was extracted using the Qiagen RNeasy® kit (Qiagen, Valencia, CA, USA), according to the manufacturer's instructions. The concentration and purity of RNA were determined by measuring spectrophotometric absorption at 260 to 280 nm. To verify the integrity of the total RNA, the samples were electrophoresed on a 1% agarose gel in RNA loading buffer. A pool of total RNA from xenografts collected on days 0, 1, 2 and 4, treated with E2 only, was used as the reference population for all cDNA microarray hybridisations. This provided an internal standard when compared with each experimental sample.
Probe preparation, labelling, hybridisation and scanning of microarrays
Total RNA (100 ug), spiked with a bacterial-RNA mixture for control, was used to prepare direct Cy3- and Cy5-labelled first-strand cDNA probes using a single-base anchored oligo dT17 primer (Sigma, St. Louis, MO, USA) and Superscript II reverse transcriptase (Invitrogen, Carlsbad, CA, USA). Unincorporated nucleotides were removed using QIAquick PCR purification kit (QIAGEN, Valencia, CA, USA) and Cy3- and Cy5-labelled probes were coprecipitated with 16 μg human Cot 1 DNA (Invitrogen, Carlsbad, CA, USA) and 8 μg polyA (Sigma, St. Louis, MO, USA). The pellets were resuspended in 8 μl of H2O and 40 μl of hybridization buffer (5 × Saline-Sodium Citrate (SSC), 6 × Denhardt's solution, 60 mM Tris HCl pH 7.6, 0.12% sarkosyl, 48% formamide) boiled for five minutes and cooled at room temperature for 10 minutes. The mix was overlaid with a coverslip and hybridized at 47°C for 12 to 24 hours in a humidified atmosphere to Sanger Hver 1.3.1 cDNA microarrays (Sanger, Cambridge, UK) as part of the CRUK/LICR Microarray Consortium, contain 9,930 sequence-validated cDNA clones representing approximately 6,000 unique sequences. Microarrays were washed sequentially with 2 × SSC, 0.1 × SSC/0.1% SDS, and 0.1 × SSC and were air-dried by briefly spinning in a centrifuge to remove excess liquid. Fluorescent images of hybridized microarrays were captured using a ScanArray Express 3.0 scanner and ScanArray software (both from Perkin Elmer Waltham, MA, USA).
Analysis of microarray data
Comparisons were made between pooled E2 only-treated controls and E2 plus tamoxifen-treated samples across the following time points - days 1, 2, 4, 7 and 14 - and included reciprocal dye labelling to exclude gene-specific dye bias. Expression ratios (Cy5/Cy3) were calculated following background correction using the R programming language  and the BioConductor  package limma  to account for dye bias. Intensity dependant (Loess) and quantile normalization were also performed. Fold changes were calculated as the relative mean difference between treated and untreated dye-swap replicates. Normalized data and raw gene expression files are publicly available from in NCBI's Gene Expression Omnibus (GEO)  and are accessible through GEO Series accession number GSE22386. Clustering was performed using the Cluster and TreeView  programs. Kaplan Meier analysis was performed using SPSS version 14 (an IBM Company, Chicago, IL, USA). Estrogen-response elements (EREs) were identified using the Dragon program . Genes with the greatest prognostic power were identified using supervized principle components analysis  using version 3.5 of BRB ArrayTools  as previously described . The Database for Annotation, Visualization and Integrated Discovery (DAVID)  was used to identify KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways and Gene Ontology terms that were significantly over-represented in gene lists above the level expected by chance.
Validation of targets by quantitative RT-PCR
The expression of putative tamoxifen-regulated genes in the ER-positive cell lines ZR-75-1 and MCF7 and the ER-negative cell line MDA-MB-231 was performed by quantitative RT-PCR. Cells were maintained as outlined above. A specific set of primers was designed for each target [see Additional File 1]. Total RNA was extracted from log-phase cells using TRI reagent (Sigma, Poole, UK) following the manufacturer's instructions and treated with DNAse I (Roche, Indianapolis, IN, USA). RNA was analyzed by real-time RT-PCR using Rotorgene (Corbett Research, San Francisco, CA, USA) and the QuantiTect SYBR Green system (QIAGEN, Valencia, CA, USA) according to the manufacturers' instructions. Thermal cycling conditions were as follows: RT at 50°C for 30 minutes; PCR: polymerase activation 95°C for 15 minutes, followed by 40 cycles of denaturation at 94°C for 15 seconds, annealing at 57°C for 30 seconds and extension at 72°C for 45 seconds. After a final extension at 72°C for 5 minutes, the melt profile was obtained across a 65°C to 99°C ramp, with 5 second ramps of 1°C. All reactions were performed in triplicate for standard curve samples and in quadruplicate for experimental and negative (no template) samples. Analysis and quantification was performed using Rotorgene v6 software. Relative quantification was calculated by extrapolation of the standard curve and calculation of ratio levels compared with β-Actin.
All experiments involving human tissues were conducted with the permission of the local medical ethics advisory board. A series of women over the age of 70 years with large operable or locally advanced primary breast cancer without metastatic disease presenting to the Edinburgh Breast Unit between October 1991 and October 1993 have previously been described . All had tumors greater than 2 cm in maximum diameter, confirmed as ER-positive invasive breast cancer. All patients received 20 mg tamoxifen daily for three months. Tumor size was monitored by ultrasound measurements, and clinical response defined as the percentage volume reduction between the initial and final tumor volumes at three months.
Formalin fixed paraffin-embedded blocks from the initial wedge biopsy and at definitive loco-regional surgery three months later were available for 28 of these patients and 3 μm tissue sections were cut. Formalin-fixed paraffin embedded (FFPE) blocks were also available from the original parallel xenograft study, which analyzed proliferation and apoptosis changes after tamoxifen treatment in the ZR-75-1 xenograft . Sections were deparaffinized and rehydrated by standard methods and endogenous peroxidase activity blocked by incubation in 3% H2O2 for 30 minutes. Sections were immersed in citrate buffer (0.005 M, pH 6.0) and microwaved three times for five minutes and then allowed to stand for 20 minutes. Slides were washed in 0.05 M Tris/NaCl buffer (pH 7.6) and then incubated in 20% FCS for 10 minutes prior to the addition of the primary antibodies at room temperature for 90 minutes in a humidified container. Optimal conditions for antigen retrieval, and primary antibody dilutions were previously determined for each antibody, as follows: trefoil factor (TFF) 3 (1/5, Calbiochem, Nottingham, UK), PDZK1 (1/10, Abnova, Taipei City Taiwan), insulin growth factor receptor binding protein (IGFBP) 4 (17661, US Biological, Swampscott, MA, USA; 1:3 dilution) and IGFBP5 (Ab4255, Abcam Cambridge, MA, USA; 1:300 dilution). After primary antibody incubation, sections were washed in Tris/NaCl buffer for 10 minutes. A Streptavidin-biotin multilink method (StrAviGen Multilink kit; Biogenex, San Ramon, CA, USA) was used for detection. The sections were incubated with secondary multilink antibody (1:20 dilution for 30 minutes) followed by a horseradish-peroxidase-labelled streptavidin complex (1:20 dilution for 30 minutes) at room temperature. Diaminobenzidine tetrachloride was applied for five minutes prior to washing in water for two minutes. Slides were then counterstained in hematoxylin, dehydrated and mounted. Expression was measured using a scoring system consisting of the product of the percentage of positive cells and their intensity of staining (0 to 3) producing a Histoscore ranging from 0 to 300. All tumor cells in the section were counted in the scoring system. Sections were scored by three independent readers and mean values obtained. Where initial scoring produced a value divergent by more than 10%, these sections were rescored until agreement was reached.
Dynamic changes in gene expression produced by tamoxifen
The effect of tamoxifen on tumor volume growth and gene expression were studied on days 1, 2, 4, 7 and 14 after initiation of tamoxifen treatment and compared with tumors grown in the absence of tamoxifen. Tumor volumes were expressed relative to the initial tumor volume (Figure 1a). A reduction in tumor volume was clearly evident at day 7 and by day 14 the curves had significantly diverged. The graphs are significantly different at day 14 and are diverging by day 7 (P < 0.05; Student's t-test).
Across the five time points, 333 probes representing 253 genes showed evidence of at least a 1.5-fold change in level of expression (using a P ≤ 0.05) [full list in the Additional File 1]. There was good agreement between the expression levels of xenograft replicates at most time points and the pattern of expression of these genes over the five time points was most consistently separated into six sets using hierarchical clustering (Figure 1b). These six sets of differentially expressed genes can be divided into three general groups: early/transient response (sets 1 and 2), variable response (sets 3 and 4) and continuous/late response (sets 5 and 6), relative to untreated samples. The early/transient-response genes were repressed relative to untreated samples, the variable-response genes were initially induced and then repressed and the continuous/late-response genes were both repressed (set 5) and induced (set 6). A large percentage of the genes in set 1 were very strongly associated with cell cycle regulation, among them AURKA, BUB1B, CCNA2, CDC25B, CDKN3, CENPF, CDC28, CKS2, DLG7, MKI67, NEK2, PRC1, STMN1, TACC3, UBE2C and ZWINT. BUB1, CKS2, PRC1, UBE2C and ZWINT have previously been shown to be estrogen-regulated in model systems . Set 2 genes included mini-chromosome maintenance (MCM) 2 and MCM6, components of the replication fork , which may account for a primary response soon after treatment reducing DNA replication and regulation. Another member of set 2 was IGFBP4, which has been widely detected in breast tumors and cell lines, and previously correlated with ER expression .
Most of the variable-response genes in sets 3 and 4 responded rapidly to tamoxifen treatment, although they were both up- and down-regulated with some variation between replicates. The genes in set 3 were predominantly involved in cell proliferation, adhesion, and apoptosis including BTG2, MYB, MYBL2 and CELSR1, whereas genes in set 4, such as IRS and IGFR1 are involved in insulin receptor signalling. Set 5 and set 6 represent genes with the greatest down- and up-regulation at day 14, respectively (Figure 1b). Set 5 contained many classical ER-response genes including TFF1, TFF3 and MYC. Serpins A1, A4 and A6 were also strongly down-regulated. These genes play a key role in the control of tissue homeostasis and have previously been shown to be up-regulated in response to E2 in normal human breast tissue . The cluster of up-regulated genes in set 6 was the largest cluster representing a wide variety of signalling pathways and processes. EREs were found in the promoter regions of a similar proportion (34 to 42%) of all six clusters of genes [genes shown in bold in the Additional File 1]. Studies by Carroll and colleagues have shown that ERs only sometimes regulate genes using EREs from proximal promoter regions and generally use distal enhancers and other binding sequences, such as Forkhead binding sites . The observation of 34 to 42% of genes containing EREs in their promoter regions is consistent with these studies.
Tamoxifen response compared with the response to estradiol over time
Many of the genes identified as changing in response to tamoxifen have also been identified in previous single time-point experiments, either in the opposite direction in response to 17β-E2 or in the same direction with tamoxifen in both in vivo and in vitro studies [8, 30]. In order to establish whether the dynamic changes observed in this study reflected the reverse of the response to E2 over time, we compared our results with those from an in vitro time-course experiment, which also utilized the ZR-75-1 cell line . Although that study had 12 time-points, with the final one being at 32 hours following addition of E2, the vast majority of genes showed the expected reciprocal changes in expression to those seen in the six clusters for the initial time points following treatment with tamoxifen in the present study (Figure 1c).
Gene expression changes in vitro
To obtain further confirmation that the expression changes observed in response to tamoxifen were valid, 15 genes were selected for in vitro validation. These were analyzed in ZR-75-1 cells treated with either 0.1 nM E2 or 1 μM tamoxifen or both agents together to assess whether the genes were not only tamoxifen-regulated but also estrogen-regulated and whether tamoxifen was antagonising the estrogen-modulation or working via some other mechanism. A second ERα-responsive cell line, the MCF-7 line, was also used to assess whether the expression changes could be observed in an independent genotype. The ERα-negative cell line, the MDA-MB-231 line, was used to assess the specificity of these changes to involvement of ERα. Expression changes were measured at both 6 and 24 hours [Additional File 2]. The gene expression changes that were observed in vivo were also observed in these in vitro experiments and the changes seen in ZR-75-1 cells were mirrored in MCF-7 cells (Figure 2). Rather than reversing the expression change produced by E2, IER3 produced a greater change in the same direction. In contrast, there were no significant changes for any of these genes in the MDA-MB-231 cell line (Figure 2).
Dynamic changes in protein expression within the ZR-75-1 xenograft model
Four candidate genes were selected to evaluate whether expression changes at the protein level over time within this xenograft model were consistent with those seen at the gene expression level (Figure 3a). MCM2 and CKS2 were chosen as examples of early/transiently changing cell cycle-associated genes (sets 1 and 2). IGFBP5 and TFF3 were selected as examples of late/continuously up- and down-regulated genes, respectively (sets 5 and 6). Sections of the xenografts were assessed by semi-quantitative immunohistochemistry and histoscores related to the initial values. The protein expression of MCM2 and CKS2 had a similar profile to that seen at the gene expression level, although protein expression was higher at day 21 than at 0 and 14 days with gene expression. TFF3 expression decreased with time whereas IGFBP5 expression increased (Figure 3a).
Protein expression changes in breast cancers treated with tamoxifen
To establish whether the change in expression of identified proteins correlated with changes seen in vivo, a series of primary breast cancers for which material was available pre- and post-tamoxifen treatment and in which the response at three months had been measured were analyzed by immunohistochemistry. Four proteins were selected (IGFBP5, TFF3, IGFBP4 and PDZK1) to represent the early/transient and later/continuous patterns of expression seen at the transcript. Breast cancer samples pre- and post-tamoxifen treatment were available for 28 patients and information on the percentage change in tumor volume was known. Change in histoscores in pre- and post-treatment paired samples were compared with the change in tumor volume (Figure 3b). The change in IGFBP5, TFF3 or both was significantly associated with change in tumor volume (P = 0.0135, 0.018 and 0.0002, respectively; Spearman rank test). This contrasted with data for IGFBP4 and PDZK1 where there was no significant association. PDZK was selected as a known estrogen-regulated gene that has been identified within a number of clinical data sets [32, 33].
Are the dynamically changing genes able to predict prognosis?
To evaluate whether the genes identified as dynamically changing in response to tamoxifen are associated with long-term follow up we downloaded four Affymetrix primary breast tumor datasets [11, 34, 35] from the NCBI GEO for patients that had all been treated with tamoxifen and for whom corresponding outcome data were available (Table 1). Affymetrix probesets representing the genes in the six gene sets with similar temporal profiles of expression were identified and clustered to separate tumors with high or low expression of each set of representative probesets [Additional File 3], as described previously . Kaplan Meier plots were generated and log rank (Mantel-Cox) statistics calculated to see if the level of these sets of genes could discriminate between patients with good or poor outcomes. The set 1 cluster of genes was highly prognostic with all four datasets. Set 2, set 4 and all six gene sets combined also had some predictive power, although this was not consistent across the four datasets (Table 1). The genes driving this prognostic separation appear to be those involved with cell cycle and proliferation, patients with high levels of these genes at presentation are known to be at high risk of recurrence . Additional File 3 illustrates the level of expression of the tamoxifen-response genes in primary tumors at presentation. Of the genes in set 1, one-third (11 out of 32) are represented in the 97 gene Genomic Grade Index that is undoubtedly associated with prognosis ; set 1 genes were also able to clearly distinguish between the histological grade of the tumors.
Do the genes that are most changed at independent time points following treatment with tamoxifen predict prognosis?
Different numbers of probes were identified to be significantly differentially expressed at each time point. In order to compare the relative prognostic value of a profile of genes that are most changed at an individual time point following treatment, we identified lists of 50 probes with the greatest fold changes (up or down) at each time point among the list of 333 most changed probes as described above. Some probes were most changed at more than one time point [see Additional File 1]. Patients whose gene expression profile at presentation was more like that of xenograft tumors following treatment had a poorer prognosis [see Additional File 3]. Figure 4 demonstrates that profiles of the most differentially expressed genes at two, four and seven days following tamoxifen treatment were able to predict prognosis, while lists of genes most changed initially (day 1) or later (day 14) cannot. None of these most changed gene lists were significantly prognostic in two datasets [37, 38] of ER-positive tumors that did not receive adjuvant therapy (Figure 4). Supervized principle components analysis  was also used to identify which genes within the profiles have the greatest prognostic power. The genes changed at each of the five independent time points were dominated by the late/continuously up- and down-regulated genes (sets 1 and 2), with early changes (sets 5 and 6) less represented and the transient changes hardly represented at all [see Additional File 1]. Known estrogen-response genes including IGFBP5, TFF3, TFF1, PDZK1 and SERPINA genes appear to dominate in prognostic performance. IGFBP5 expression is higher in patients with poor prognosis, as noted previously [29, 39]; however, it is not significantly changed at day one, but is at subsequent time points, as seen at the protein level (Figure 3a). The heatmaps in Additional File 3 also indicate that IGFBP5 may be a good biomarker of outcome on tamoxifen.
This study is, to the best of our knowledge, the first to define tamoxifen-regulated gene expression profiles at multiple time points after long-term anti-estrogen treatment in an in vivo model of ER-positive breast cancer. The use of multiple time points over the 14-day period allowed analysis of the temporal patterns of gene expression profiles. Three basic patterns of change were observed; early/transient changes, continuous/late changes and more variable changes. The pattern of expression of representatives of these sets of genes was confirmed by quantitative RT-PCR and at the protein level by semi-quantitative immunohistochemistry. The changes observed in the expression of IGFBP5 and TFF3 correlated with reductions in tumor volume in primary tumors treated with tamoxifen. Two different approaches were used to evaluate whether those genes for which there was clear evidence of tamoxifen-induced changes in expression level would themselves be prognostic for patients treated with adjuvant tamoxifen. The early/transient pattern of gene expression associated with a reduction in cell cycle/proliferation genes and genes that were most differentially expressed on days two, four and seven were able to predict prognosis of primary breast tumors treated with tamoxifen. The IMPACT (Immediate Preoperative Arimidex Compared to Tamoxifen) trial demonstrated that Ki67 level two weeks after treatment was predictive of long-term outcome . The timing of measurement of gene expression changes appears to be critical for certain groups of genes. A number of neoadjuvant 'window of opportunity studies' are underway to characterize changes in gene expression in response to treatment and establish whether clinical response after a couple of weeks or several months is predictive of long-term outcome. Studies such as that described here may provide insights as to when significant changes are detectable and which genes may represent good markers of response. It also highlights the risk that clinical snapshots of treatment could miss informative changes in expression.
There have been multiple short-term studies of in vitro profiling after estrogen treatment, predominantly in MCF-7 cells [41, 42] and also in ZR-75-1 cells [4, 10]. In a study using the T47D model of ER-positive breast cancer, Harvell and colleagues stated that E2 regulates different genes in human breast tumor xenografts compared with the identical cells in culture . However, Creighton and colleagues found that genes regulated by estrogen in breast tumor cells in vitro are similarly regulated in vivo in tumor xenografts and human tumors . Disparities between approaches may be the result of differences in time points as well as the differences in microenvironment. Our in vitro study allowed exploration of whether tamoxifen's effects were antagonistic to estrogen or not. Fifteen genes were selected and all were modulated by estrogen in vitro. Of these genes, tamoxifen reversed the estrogen modulation in 14 cases but not for IER3. This gene was of particular interest in that while it was estrogen up-regulated, tamoxifen produced a greater up-regulation, it was also continuously up-regulated at all five time points and among the genes of set 6.
Our previous study  demonstrated early changes in apoptotic and mitotic indices (days 2 and 4) predated tumor volume changes, we speculate that the earlier/transient expression changes observed are more likely to be causative and primary events for tumor volume inhibition whereas later/continuous expression changes are possibly only consequential and secondary to the volume changes. These may represent changes in stromal elements and infiltrating cell populations. Ongoing studies are seeking to develop a putative model of how the early/transient changes interact with the later/continuous changes.
Many of the breast cancer gene expression signatures that have previously been developed highlight a number of genes associated with cell cycle and proliferation [35, 36, 43–46], which has been suggested to be largely a reflection of tumor grade. These genes appear to have most prognostic value for ER-positive breast tumors, generally differentiating between luminal A and luminal B subtypes, both prior to or following treatment with tamoxifen  or chemotherapy [47, 48]. Our results are consistent with this idea and the suggestion that a lower response to E2 or growth factor signaling , also a feature of luminal B tumors, may also be prognostic. It is not clear from our results the extent to which prognosis or prediction of response to therapy is an intrinsic property of tumors.
Relatively high levels of the set 1 pattern of dynamically changing cell cycle/proliferation genes at presentation in primary tumors was associated with poor prognosis; however, it was relatively low levels of the classically up-regulated E2-response genes (down-regulated by tamoxifen), such as TFF1, TFF3, AREG and IGFBP4, at presentation that were among the genes most changed at day 4 and associated with poor prognosis. Conversely, a reduction in TFF3 (or an increase in IGFBP5) following tamoxifen treatment in vivo and the protein levels in primary tumors correlated with a reduction in tumor volume in the 28 patients treated with tamoxifen for three months. This apparent contradiction between the direction of change in genes upon treatment and their relative level in primary tumors as long-term predictors of outcome may be due to the complexity of estrogen signaling, the agonistic and antagonistic roles of estrogen and tamoxifen on the ER and/or a difference between short- and long-term effects on both tumors and normal tissues. We also recently demonstrated that proliferation genes were strongly down-regulated following treatment with the mTOR inhibitor, everolimus, despite these often being considered markers of prognosis .
The aim of this study was to assess the dynamic response to tamoxifen, not to find the definitive tamoxifen-response signature or biomarker. A better test of the tamoxifen-response genes in primary tumors would be a dataset from a neoadjuvant 'window study'  of gene expression before and after tamoxifen with both clinical or pathological endpoints and long-term follow-up. It would be interesting to measure gene expression at multiple time points in a number of different cell line xenograft models or primary tumors in order to fully investigate patient-patient variation in temporal response to tamoxifen. This approach would also benefit from single-color microarrays in order to evaluate the relative merits of pre- and post-treatment samples avoiding the limitation of using comparative two-colour cDNA arrays, as in this study. We did examine gene expression of the different response patterns (sets 1 to 6) and individual time points in matched before and after breast biopsies from patients treated with 14 days of neoadjuvant letrozole  and found largely consistent changes for most genes with those of the in vivo study in the majority of cases (data not shown). Further work is required to fully assess how the response to different hormonal therapies and short-term molecular changes correlate with long-term outcome. We have previously demonstrated that estrogen-regulated gene expression predicts response to endocrine therapy in patients with ovarian cancer , and in this study we demonstrate for the first time that tamoxifen-response genes identified from a xenograft breast cancer model with different profiles of expression can predict prognosis in primary tumors treated with tamoxifen.
The genes highlighted in this study are now being explored in clinical material collected by biopsy from patients pre- and post-treatment with tamoxifen and who are known to have either responded to or progressed on treatment. This will help determine which of the genes identified in this study have the potential to be predictive markers of response. This study also suggests that future studies searching for genes predictive of outcome on therapy could perhaps be informed by studies that identify which genes demonstrate early dynamic response to therapy, rather than those with sustained changes. This is reminiscent of data from early positron emission tomography (PET) scans that suggest the patients with the best outcome on therapy are those with pronounced early reduction in PET signal [52, 53].
Both early/transient/proliferation-response genes and continuous/late/estrogen-response genes are able to predict prognosis of primary breast tumors in a dynamic manner. Temporal expression of therapy-response genes is clearly an important factor in the response to endocrine therapy in breast tumors which has significant implications for the timing of biopsies in neoadjuvant biomarker studies.
CDC28 protein kinase regulatory subunit
Database for Annotation, Visualization and Integrated Discovery
Dulbecco's modified Eagle's medium
estrogen receptor alpha
fetal calf serum
insulin growth factor receptor binding protein
- NCBI GEO:
National Centre for Biotechnology Information Gene Expression Omnibus
positron emission tomography
reverse-transcription polymerase chain reaction
Ali S, Coombes RC: Endocrine-responsive breast cancer and strategies for combating resistance. Nat Rev Cancer. 2002, 2: 101-112. 10.1038/nrc721.
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-1717. 10.1016/S0140-6736(05)66544-0.
Cameron DA, Ritchie AA, Langdon S, Anderson TJ, Miller WR: Tamoxifen induced apoptosis in ZR-75 breast cancer xenografts antedates tumour regression. Breast Cancer Res Treat. 1997, 45: 99-107. 10.1023/A:1005850827825.
Soulez M, Parker MG: Identification of novel oestrogen receptor target genes in human ZR75-1 breast cancer cells by expression profiling. J Mol Endocrinol. 2001, 27: 259-274. 10.1677/jme.0.0270259.
Charpentier AH, Bednarek AK, Daniel RL, Hawkins KA, Laflin KJ, Gaddis S, MacLeod MC, Aldaz CM: Effects of estrogen on global gene expression: identification of novel targets of estrogen action. Cancer Res. 2000, 60: 5977-5983.
Cicatiello L, Scafoglio C, Altucci L, Cancemi M, Natoli G, Facchiano A, Iazzetti G, Calogero R, Biglia N, De Bortoli M, Sfiligoi C, Sismondi P, Bresciani F, Weisz A: A genomic view of estrogen actions in human breast cancer cells by expression profiling of the hormone-responsive transcriptome. J Mol Endocrinol. 2004, 32: 719-775. 10.1677/jme.0.0320719.
Inoue A, Yoshida N, Omoto Y, Oguchi S, Yamori T, Kiyama R, Hayashi S: Development of cDNA microarray for expression profiling of estrogen-responsive genes. J Mol Endocrinol. 2002, 29: 175-192. 10.1677/jme.0.0290175.
Creighton CJ, Cordero KE, Larios JM, Miller RS, Johnson MD, Chinnaiyan AM, Lippman ME, Rae JM: Genes regulated by estrogen in breast tumor cells in vitro are similarly regulated in vivo in tumor xenografts and human breast tumors. Genome Biol. 2006, 7: R28-10.1186/gb-2006-7-4-r28.
Harvell DM, Richer JK, Allred DC, Sartorius CA, Horwitz KB: Estradiol regulates different genes in human breast tumor xenografts compared to the identical cells in culture. Endocrinology. 2006, 147: 700-713. 10.1210/en.2005-0617.
Mutarelli M, Cicatiello L, Ferraro L, Grober OM, Ravo M, Facchiano AM, Angelini C, Weisz A: Time-course analysis of genome-wide gene expression data from hormone-responsive human breast cancer cells. BMC Bioinformatics. 2008, 9 (Suppl 2): S12-10.1186/1471-2105-9-S2-S12.
Zhang Y, Sieuwerts AM, McGreevy M, Casey G, Cufer T, Paradiso A, Harbeck N, Span PN, Hicks DG, Crowe J, Tubbs RR, Budd GT, Lyons J, Sweep FC, Schmitt M, Schittulli F, Golouh R, Talantov D, Wang Y, Foekens JA: The 76-gene signature defines high-risk patients that benefit from adjuvant tamoxifen therapy. Breast Cancer Res Treat. 2008, 116: 303-309. 10.1007/s10549-008-0183-2.
Musgrove EA, Sergio CM, Loi S, Inman CK, Anderson LR, Alles MC, Pinese M, Caldon CE, Schutte J, Gardiner-Garden M, Ormandy CJ, McArthur G, Butt AJ, Sutherland RL: Identification of functional networks of estrogen- and c-Myc-responsive genes and their relationship to response to tamoxifen therapy in breast cancer. PLoS ONE. 2008, 3: e2987-10.1371/journal.pone.0002987.
Loi S, Haibe-Kains B, Desmedt C, Wirapati P, Lallemand F, Tutt AM, Gillet C, Ellis P, Ryder K, Reid JF, Daidone MG, Pierotti MA, Berns EM, Jansen MP, Foekens JA, Delorenzi M, Bontempi G, Piccart MJ, Sotiriou C: Predicting prognosis using molecular profiling in estrogen receptor-positive breast cancer treated with tamoxifen. BMC Genomics. 2008, 9: 239-10.1186/1471-2164-9-239.
Vendrell JA, Robertson KE, Ravel P, Bray SE, Bajard A, Purdie CA, Nguyen C, Hadad SM, Bieche I, Chabaud S, Bachelot T, Thompson AM, Cohen PA: A candidate molecular signature associated with tamoxifen failure in primary breast cancer. Breast Cancer Res. 2008, 10: R88-10.1186/bcr2158.
Sims AH, Bartlett JM: Approaches towards expression profiling the response to treatment. Breast Cancer Res. 2008, 10: 115-10.1186/bcr2196.
Cameron DA, Ritchie AA, Miller WR: The relative importance of proliferation and cell death in breast cancer growth and response to tamoxifen. Eur J Cancer. 2001, 37: 1545-1553. 10.1016/S0959-8049(01)00166-6.
Ihaka R, Gentleman R: R: a language for data analysis and graphics. Journal of Computational and Graphical Statistics. 1996, 5: 299-314. 10.2307/1390807.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.
Smyth GK, Speed T: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.
Barrett T, Suzek TO, Troup DB, Wilhite SE, Ngau WC, Ledoux P, Rudnev D, Lash AE, Fujibuchi W, Edgar R: NCBI GEO: mining millions of expression profiles--database and tools. Nucleic Acids Res. 2005, 33: D562-566. 10.1093/nar/gki022.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.
Bajic VB, Tan SL, Chong A, Tang S, Strom A, Gustafsson JA, Lin CY, Liu ET: Dragon ERE Finder version 2: A tool for accurate detection and analysis of estrogen response elements in vertebrate genomes. Nucleic Acids Res. 2003, 31: 3605-3607. 10.1093/nar/gkg517.
Bair E, Tibshirani R: Semi-supervised methods to predict patient survival from gene expression data. PLoS Biol. 2004, 2: E108-10.1371/journal.pbio.0020108.
BRB ArrayTools. [http://linus.nci.nih.gov/BRB-ArrayTools.html]
Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57. 10.1038/nprot.2008.211.
Cameron DA, Keen JC, Dixon JM, Bellamy C, Hanby A, Anderson TJ, Miller WR: Effective tamoxifen therapy of breast cancer involves both antiproliferative and pro-apoptotic changes. Eur J Cancer. 2000, 36: 845-851. 10.1016/S0959-8049(00)00013-7.
Walker G, MacLeod K, Williams AR, Cameron DA, Smyth JF, Langdon SP: Estrogen-regulated gene expression predicts response to endocrine therapy in patients with ovarian cancer. Gynecol Oncol. 2007, 106: 461-468. 10.1016/j.ygyno.2007.05.009.
Labib K, Tercero JA, Diffley JF: Uninterrupted MCM2-7 function required for DNA replication fork progression. Science. 2000, 288: 1643-1647. 10.1126/science.288.5471.1643.
Walker G, MacLeod K, Williams AR, Cameron DA, Smyth JF, Langdon SP: Insulin-like growth factor binding proteins IGFBP3, IGFBP4, and IGFBP5 predict endocrine responsiveness in patients with ovarian cancer. Clin Cancer Res. 2007, 13: 1438-1444. 10.1158/1078-0432.CCR-06-2245.
Wilson CL, Sims AH, Howell A, Miller CJ, Clarke RB: Effects of oestrogen on gene expression in epithelium and stroma of normal human breast tissue. Endocr Relat Cancer. 2006, 13: 617-628. 10.1677/erc.1.01165.
Carroll JS, Liu XS, Brodsky AS, Li W, Meyer CA, Szary AJ, Eeckhoute J, Shao W, Hestermann EV, Geistlinger TR, Fox EA, Silver PA, Brown M: Chromosome-wide mapping of estrogen receptor binding reveals long-range regulation requiring the forkhead protein FoxA1. Cell. 2005, 122: 33-43. 10.1016/j.cell.2005.05.008.
Ghosh MG, Thompson DA, Weigel RJ: PDZK1 and GREB1 are estrogen-regulated genes expressed in hormone-responsive breast cancer. Cancer Res. 2000, 60: 6367-6375.
Mackay A, Urruticoechea A, Dixon JM, Dexter T, Fenwick K, Ashworth A, Drury S, Larionov A, Young O, White S, Miller WR, Evans DB, Dowsett M: Molecular response to aromatase inhibitor treatment in primary breast cancer. Breast Cancer Res. 2007, 9: R37-10.1186/bcr1732.
Loi S, Haibe-Kains B, Desmedt C, Lallemand F, Tutt AM, Gillet C, Ellis P, Harris A, Bergh J, Foekens JA, Klijn JG, Larsimont D, Buyse M, Bontempi G, Delorenzi M, Piccart MJ, Sotiriou C: Definition of clinically distinct molecular subtypes in estrogen receptor-positive breast carcinomas through genomic grade. J Clin Oncol. 2007, 25: 1239-1246. 10.1200/JCO.2006.07.1522.
Sotiriou C, Wirapati P, Loi S, Harris A, Fox S, Smeds J, Nordgren H, Farmer P, Praz V, Haibe-Kains B, Desmedt C, Larsimont D, Cardoso F, Peterse H, Nuyten D, Buyse M, Van de Vijver MJ, Bergh J, Piccart M, Delorenzi M: Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis. J Natl Cancer Inst. 2006, 98: 262-272. 10.1093/jnci/djj052.
Liu J, Campen A, Huang S, Peng SB, Ye X, Palakal M, Dunker AK, Xia Y, Li S: Identification of a gene signature in cell cycle pathway for breast cancer prognosis using gene expression profiling data. BMC Med Genomics. 2008, 1: 39-10.1186/1755-8794-1-39.
Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, Meijer-van Gelder ME, Yu J, Jatkoe T, Berns EM, Atkins D, Foekens JA: Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365: 671-679.
Desmedt C, Piette F, Loi S, Wang Y, Lallemand F, Haibe-Kains B, Viale G, Delorenzi M, Zhang Y, d'Assignies MS, Bergh J, Lidereau R, Ellis P, Harris AL, Klijn JG, Foekens JA, Cardoso F, Piccart MJ, Buyse M, Sotiriou C: Strong time dependence of the 76-gene prognostic signature for node-negative breast cancer patients in the TRANSBIG multicenter independent validation series. Clin Cancer Res. 2007, 13: 3207-3214. 10.1158/1078-0432.CCR-06-2765.
Mita K, Zhang Z, Ando Y, Toyama T, Hamaguchi M, Kobayashi S, Hayashi S, Fujii Y, Iwase H, Yamashita H: Prognostic significance of insulin-like growth factor binding protein (IGFBP)-4 and IGFBP-5 expression in breast cancer. Jpn J Clin Oncol. 2007, 37: 575-582. 10.1093/jjco/hym066.
Ellis MJ, Tao Y, Luo J, A'Hern R, Evans DB, Bhatnagar AS, Chaudri Ross HA, von Kameke A, Miller WR, Smith I, Eiermann W, Dowsett M: Outcome prediction for estrogen receptor-positive breast cancer based on postneoadjuvant endocrine therapy tumor characteristics. J Natl Cancer Inst. 2008, 100: 1380-1388. 10.1093/jnci/djn309.
Cunliffe HE, Ringner M, Bilke S, Walker RL, Cheung JM, Chen Y, Meltzer PS: The gene expression response of breast cancer to growth regulators: patterns and correlation with tumor expression profiles. Cancer Res. 2003, 63: 7158-7166.
Coser KR, Chesnes J, Hur J, Ray S, Isselbacher KJ, Shioda T: Global analysis of ligand sensitivity of estrogen inducible and suppressible genes in MCF7/BUS breast cancer cells by DNA microarray. Proc Natl Acad Sci USA. 2003, 100: 13994-13999. 10.1073/pnas.2235866100.
Sotiriou C, Neo SY, McShane LM, Korn EL, Long PM, Jazaeri A, Martiat P, Fox SB, Harris AL, Liu ET: Breast cancer classification and prognosis based on gene expression profiles from a population-based study. Proc Natl Acad Sci USA. 2003, 100: 10393-10398. 10.1073/pnas.1732912100.
Sotiriou C, Pusztai L: Gene-expression signatures in breast cancer. N Engl J Med. 2009, 360: 790-800. 10.1056/NEJMra0801289.
Colozza M, Azambuja E, Cardoso F, Sotiriou C, Larsimont D, Piccart MJ: Proliferative markers as prognostic and predictive tools in early breast cancer: where are we now?. Ann Oncol. 2005, 16: 1723-1739. 10.1093/annonc/mdi352.
Dai H, van't Veer L, Lamb J, He YD, Mao M, Fine BM, Bernards R, van de Vijver M, Deutsch P, Sachs A, Stoughton R, Friend S: A cell proliferation signature is a marker of extremely poor outcome in a subpopulation of breast cancer patients. Cancer Res. 2005, 65: 4059-4066. 10.1158/0008-5472.CAN-04-3953.
Liedtke C, Hatzis C, Symmans WF, Desmedt C, Haibe-Kains B, Valero V, Kuerer H, Hortobagyi GN, Piccart-Gebhart M, Sotiriou C, Pusztai L: Genomic grade index is associated with response to chemotherapy in patients with breast cancer. J Clin Oncol. 2009, 27: 3185-3191. 10.1200/JCO.2008.18.5934.
Jones RL, Salter J, A'Hern R, Nerurkar A, Parton M, Reis-Filho JS, Smith IE, Dowsett M: The prognostic significance of Ki67 before and after neoadjuvant chemotherapy in breast cancer. Breast Cancer Res Treat. 2009, 116: 53-68. 10.1007/s10549-008-0081-7.
Loi S, Sotiriou C, Haibe-Kains B, Lallemand F, Conus NM, Piccart MJ, Speed TP, McArthur GA: Gene expression profiling identifies activated growth factor signaling in poor prognosis (Luminal-B) estrogen receptor positive breast cancer. BMC Med Genomics. 2009, 2: 37-10.1186/1755-8794-2-37.
Sabine VS, Sims AH, Macaskill EJ, Renshaw L, Thomas JS, Dixon JM, Bartlett JM: Gene expression profiling of response to mTOR inhibitor everolimus in pre-operatively treated post-menopausal women with oestrogen receptor-positive breast cancer. Breast Cancer Res Treat. 2010, 122: 419-28. 10.1007/s10549-010-0928-6.
Miller WR, Larionov AA, Renshaw L, Anderson TJ, White S, Murray J, Murray E, Hampton G, Walker JR, Ho S, Krause A, Evans DB, Dixon JM: Changes in breast cancer transcriptional profiles after treatment with the aromatase inhibitor, letrozole. Pharmacogenet Genomics. 2007, 17: 813-826. 10.1097/FPC.0b013e32820b853a.
Li D, Chen JH, Wang J, Ling R, Yao Q, Wang L: Value of fused 18F-FDG PET/CT images in predicting efficacy of neoadjuvant chemotherapy on breast cancer. Ai Zheng. 2007, 26: 900-904.
Smith IC, Welch AE, Hutcheon AW, Miller ID, Payne S, Chilcott F, Waikar S, Whitaker T, Ah-See AK, Eremin O, Heys SD, Gilbert FJ, Sharp PF: Positron emission tomography using (18)F-fluorodeoxy-D-glucose to predict the pathologic response of breast cancer to primary chemotherapy. J Clin Oncol. 2000, 18: 1676-1688.
We are grateful to Cancer Research UK, Breakthrough Breast Cancer and the Charon Fund for supporting this work. We are grateful to Xu Tang for assistance with certain of the xenograft immunohistochemistry studies.
The authors declare that they have no competing interests.
SPL, DAC and DH conceived and directed the study. KT undertook the microarray and RT-PCR studies. GW, BK and SPL performed the immunohistochemistry on the clinical breast cancer cases. DAC obtained tissue blocks. AHS and LL analyzed the microarray data. MM carried out the xenograft work. DF undertook the immunohistochemistry on the xenograft material. JMD collected the clinical material. KT, AHS and SL drafted the manuscript. All authors read and approved the final manuscript.
Karen J Taylor, Andrew H Sims contributed equally to this work.
Electronic supplementary material
Additional file 1: Lists of differentially expressed genes. Microsoft Excel Workbook containing probe and gene names, plus Ensembl identifiers and mean fold changes for the 333 significantly differentially expressed probes. Also provided are lists of the 50 most changed genes at the five timepoints following treatment with tamoxifen. (XLS 108 KB)
Additional file 2: Quantitative PCR results. Gene expression in vitro measured by quantitative RT-PCR for ZR75 (royal blue), MCF7 (dark blue) and MDA-MB-231 (red) before (0), 6 and 24 hours following no treatment (C), addition of estradiol (E), tamoxifen (T) and estradiol plus tamoxifen (ET). (PDF 27 KB)
Additional file 3: Examples of heatmap clustering. Heatmaps showing the level of expression of the (a) set 1 and (b) day 4 tamoxifen-response genes in primary tumors at presentation. Patients whose expression of set 1 genes correlate with post-treatment xenograft samples have a good prognosis (blue). However, patients whose expression of genes at presentation is more like those that were differentially expressed at day 4 following tamoxifen treatment tend to have a poor prognosis (green). See Table 1 and Figure 4 for survival analysis results. (PDF 1 MB)
About this article
Cite this article
Taylor, K.J., Sims, A.H., Liang, L. et al. Dynamic changes in gene expression in vivo predict prognosis of tamoxifen-treated patients with breast cancer. Breast Cancer Res 12, R39 (2010) doi:10.1186/bcr2593
- Breast Cancer
- Gene Expression Change
- Tamoxifen Treatment
- Multiple Time Point