
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Recent Advances and Future Directions in Endocrine Manipulation of Breast Cancer |
1 Department of Genetics, Institute of Cancer Research, Norwegian Radium Hospital; 2 Department of Surgery, Ullevaal Hospital, Oslo, Norway; 3 Department of Oncology, Haukeland Hospital, Bergen, Norway; and 4 Department of Biochemistry, School of Medicine, Fujita Health University, Toyoake, Japan
Requests for reprints: Vessela N. Kristensen, Department of Genetics, Institute of Cancer Research, Norwegian Radium Hospital, Montebello 0310, Oslo, Norway. Phone: 47-22-93-44-17; Fax: 47-22-93-44-40; E-mail: nedelcheva.vessela{at}dnr.uio.no..
| ABSTRACT |
|---|
|
|
|---|
Key Words: estradiol metabolism SNP stage of disease clustering analysis
| INTRODUCTION |
|---|
|
|
|---|
With this in mind, we compared the tumor classes obtained by whole genome gene expression analysis (1) with those obtained by previously published independent quantitation of the estradiol-metabolizing enzymes (2): aromatase (CYP19), HSD1, HSD2, STS, and EST as well as the related signaling members ER
, ERß, cyclin D, and ErbB2 using ß-actin as an internal standard.
| SUMMARY OF MATERIALS AND METHODS |
|---|
|
|
|---|
+ tumor were subsequently treated with tamoxifen for 5 years or to time of relapse.
+ tumors treated with anastrozole (1 mg/d p.o.) for 15 weeks before surgery (6). Microarray expression profiles were studied in the pretreatment biopsies and in samples obtained at final surgery. The primary data tables and the image files are stored in the Stanford Microarray Database (http://genome-www5.stanford.edu/).
Absolute Quantitation of mRNA. Fluorometrical quantitation was done to determine the absolute mRNA contents using fluorescent dyelabeled primers in the presence of an internal standard for each gene: CYP19, HSD1, HSD2, ER
, ERß, STS, EST (SULT1E1), cyclin D, ErbB2, and ß-actin as described in ref. 2.
Statistical Analysis. Clinical variables as well as mRNA levels were collected in a Statistical Package for Science version 11.5 (http://www.spss.org/) file where the initial descriptive statistics, cross-tab (
2 test), and Pearson correlation analysis were done. Average-linkage hierarchical clustering was applied by using the CLUSTER program, and the results were displayed by using TREEVIEW (software available at http://genome-www5.stanford.edu//resources/restech.shtml). Absolute mRNA levels were calculated as relative to median, and clustering analysis of levels above and below the median was done using EPCLUST (http://ep.ebi.ac.uk/EP/EPCLUST/). Gene expression clusters were categorized and the distribution of the genotypes in these categories was studied using
2 statistics.
We applied a recently described analytic method called significance analysis of microarrays to search for genes that correlated with treatment response to anastrozole (7). Class prediction was done by using prediction analysis of microarrays, which is a variant of nearest-centroid classification with an automated gene selection step integrated into the algorithm (8).
| RESULTS |
|---|
|
|
|---|
Estradiol-Metabolizing Enzymes and Correlation with Whole Genome Patterns
Of the 73 tumors from series 2 (see SUMMARY OF MATERIALS AND METHODS) with known whole genome expression pattern, 64 were also characterized for mRNA expression of several estradiol-metabolizing enzymes (2) and classified into four groups: I: ER
+ERß+CYP19+, II: ER
+ERßCYP19, III: ER
ERßCYP19+, and IV: ER
ERßCYP19. When we compared the classification of these tumors by whole genome expression profiling to classification by these four groups, a statistically significant association was observed (Fig. 1). In part, this association could be easily explained through the expression of ER
: The ER
+ERß+CYP19+ and the ER
+ERßCYP19 were most often of the luminal A and B type, and the ER
ERßCYP19 group was basal like. However, there were patients in the normal like and the ERBB2/HER-2 subgroup whose tumor, although ER
+, did not express estradiol-metabolizing enzymes. The normal-like group was particularly marked by tumor samples not expressing the ER but still producing estradiol III: ER
ERßCYP19+ (Fig. 1). These findings point to a possible association between development of a tumor with a particular expression profile and its capacity to synthesize estradiol as measured by the expression of the transcripts for the necessary key enzymes.
|
clustered together and were characterized by a strong basal-like signature highly expressing keratin 5/17, cadherin 3, frizzled, apolipoprotein D, etc. The luminal epithelial tumor cluster, on the other hand, highly expressed ER
, GATA binding protein 3, and N-acetyltransferase. An evident ERBB2 cluster was observed due to the marked overexpression of the ERBB2 gene, GRB7 and PPARBP in this patient material. Further, we wanted to investigate the expression patterns of genes that had been shown previously to be up-regulated and down-regulated by estradiol in MCF-7 cell lines (11) in the tumors of the patients treated with anastrozole, when estradiol synthesis should be inhibited in the tissue. Interestingly, the tumors again segregated in two main clusters, largely representing the ER
+ versus the ER
cases, without actually having the ER
itself among the genes in the cluster (Fig. 2). This classification was mainly due to other luminal class markers such as GATA3 and TFF1 (trefoil factor 1 breast cancer). Using significance analysis of microarrays, we identified 298 genes significantly differently expressed between the partial response and the progressive disease groups. When we clustered all the samples with the 298 gene list, most of the tumor pairs clustered together, indicating some "intrinsic" response signature that does not change with the treatment as such (data not shown).
|
| DISCUSSION |
|---|
|
|
|---|
and the androgen receptor (16). Furthermore, ER+ cell lines, such as MCF-7 and ZR-75-1, showed a rapid and sustained activation of extracellular signal-related kinase in response to estradiol that was substantially prevented by wild-type but not mutant BRCA1, providing a link between BRCA1 mutation status (and the following tumor classification) and estradiol-induced signaling (17).
ER
, whose expression is one of the more important clinical variables of breast cancer, remains a key player behind the gene expression separation of the luminal and basal types of tumor classes. A possibility may be that the observed tumor classes represent separate cell types from which these tumors originate. This hypothesis has been addressed in a study comparing the expression patterns of two epithelial cell types, inner luminal secretory cells and outer contractile myoepithelial cells (18). Comprehensive expression profiles of the two normal cell types, using immunomagnetic cell separation and gene expression microarray analysis, revealed distinct expression profiles of the two cellular types in breast cancer patients. Whereas between 4% and 20% of the luminal epithelial cells are positive for ER and
50% are positive for PR, myoepithelial cells do not express ER. Regulation of the normal mammary gland function is closely related to the preservation of the tissue architecture (19). In vivo, myoepithelial and most luminal cells are in contact with the basement membrane, and communication between the luminal cells either directly or through the myoepithelial cells is essential for mammary differentiation. The proteins that form the intermediate filaments of epithelial cells, the cytokeratins, are heterodimers with very distinct expression patterns. High expression of keratin 5 and 17 was observed in the basal-like tumors characterized by ER negativity and BRCA1 mutations (4, 5). The basement membrane molecules laminin-1 and collagen IV were found to be involved in the maintenance of ER
expression in cultured nonmalignant mammary epithelial cells, and this response could be interfered with by disrupting cell/extracellular matrix adhesion (19).
The relationships between the ER and both its ligand and the enzymes that synthesize it are not well understood. We have reported previously a SNP in the 3' untranslated area of CYP19 mRNA to be associated with higher mRNA levels and increased tumor size in locally advanced breast cancers (stage III; ref. 20). Patients with stage I and II disease with no expression of CYP19, on the other hand, had short disease-free survival (P > 0.0005), possibly attributable to hormone independence and poor response to standard treatment with tamoxifen (2). We have previously reported high correlation of the expression of CYP19 (aromatase) mRNA and a switch from the adipose tissuespecific promoter to the ovary-specific promoter (20). Furthermore, the alternative usage of exon 1 was also associated with shorter disease-free survival; however, the patients with worse prognosis had no expression of CYP19 mRNA at all (2). This is a somewhat paradoxical result, because these and our previous data (20) suggest that in the course of the disease a switch of promoter occurs and results in higher expression levels. The process, however, seems to occur in a given timeframe of the disease on the way to hormone-independent breast cancer. The poor prognosis of patients with null expression of aromatase may also reflect treatment failure of the patients. This was a consecutive series of patients with predominantly stage I/II disease who had been given a standard treatment with tamoxifen. Patients who relapsed after adjuvant tamoxifen then received aromatase inhibitors. Unfortunately, a substantial number of tumors are intrinsically tamoxifen resistant, despite ER positivity, and eventually almost all breast carcinomas acquire resistance to tamoxifen (21). We have no exhaustive data on treatment response to verify whether patients with null aromatase expression in their tumors respond poorly to such standard therapeutic procedure.
Anastrozole has been shown to inhibit both total body and intratumoral aromatization in vivo (6, 22), and to cause tumor repression in the neoadjuvant setting (23). In a recent study (6), we determined the effect of primary medical treatment with anastrozole for 15 weeks on tissue estrogen levels in 12 patients with locally advanced breast cancer. Using a combined high-performance liquid chromatography-RIA method, we found that intratumoral levels of estrone, estradiol, and estrone sulfate decreased by 89% (confidence interval, 73.2-95.5), 83.4% (confidence interval, 63.2-92.5) and 72.9% (confidence interval, 47.2-86.1), respectively (6). Although all patients presented with ER+ tumors in the routine immunohistochemical analysis, in some the receptor expression was low. The estrogen levels dropped in all patients; however, the clinical response varied. The reason for lack of response to hormone manipulation in some breast cancer patients expressing the ER
is poorly understood, and larger studies are needed to validate the significance of these genes for treatment response to anastrozole.
In summary, the analysis of the data sets presented here, studying whole genome expression profiles in both consecutive early-stage breast cancer and selected locally advanced cases, both treated or untreated, suggest common tumor classification groups with different disease-free survival rates. Can this observation cast more light on the origin and on the fate of these breast cancers in terms of prognosis and treatment response? The present tumor classification suggests that there is a more "uneventful" path to develop breast cancer, which is ER+, expresses genes characteristic for the luminal cell type and has a better prognosis in terms of disease-free survival time. On the other hand, there are tumors with more dramatic presentations that are ER, often with amplification of ERBB2, often harboring TP53 mutations, and expressing genes characteristic of the basal cell type. Carriers of germ line mutations in the BRCA1 gene are predominantly of the second type. Do these two types occur independently from different cell entities of the breast (luminal and basal) comprising different types of disease? The existence of intermediate types (e.g., luminal B), which express ER but may also have amplification of ERBB2 and express some basal cell markers, may argue against this hypothesis. It is equally feasible to hypothesize that in the course of tumor development cells express markers that make them resemble either the one or the other cell type. The question would then be whether these expression programs occur simultaneously or in succession during the progression of the disease and whether the results we observe at the time point of our microarray analysis are a product of selection and if so what the driving forces of this selection are.
To be able to use microarray data for prognosis and prediction of treatment response and clinical outcome, we need to further dissect the observed clusters of genes in terms of relevant regulatory networks and incorporate hypotheses regarding functional pathways (24). We are currently performing analysis for common regulatory blocks in the separate gene clusters as well as comparing the patient data with data from cell lines subjected to similar treatments.
| OPEN DISCUSSION |
|---|
|
|
|---|
Dr. Vessela Kristensen: The proliferation genes form a nice cluster by themselves, and we can see that in both the ER+ and mostly in the ER tumors, that is, in the luminal and the basal types. It is exactly the presence of the proliferation cluster that distinguishes the ER-positive tumors with poor prognosis (luminal type B) from those with good prognosis (luminal type A).
Dr. Kent Osborne: You said that ERß was associated with a poor prognosis in tamoxifen-treated patients. ERß has been really difficult to study for lack of an antibody that is reproducible and specific and that works in an immunohistochemistry. We just submitted for publication a study looking at ERß protein expression in a series of adjuvant tamoxifen-treated patients. There, ERß was actually a good factorin fact, independently and statistically significantly so. Some of the data on ERß function are beginning to show that it puts the brake on the system, to put it in a simplistic fashion.
Dr. Donald McDonnell: In general, ER
activates and ERß modulates. What the specific roles of those two proteins are depends on what model systems you've used. There is an interesting possibility here that you're looking at mRNA. There is a negative feedback of the protein to the gene, and if the feedback is disrupted, that is, the protein is not being produced, the mRNA is going to compensate and go up. What could happen is that you're getting a high level of ERß mRNA in those tissues that are actually not expressing very much ERß protein.
Dr. Kristensen: Definitely. We observed the same with aromatase. It was in the patients with successful inhibition of the protein that an increase in mRNA levels was observed. Based on mRNA alone, we would have concluded quite wrongly.
Dr. Russo: In those cases where you compare nontreated tumor versus treated, for example, with tamoxifen, is the genomic profile of the tamoxifen-responsive tumor similar to the normal breast or is it a different profile?
Dr. Kristensen: None of these studies are actually ready to answer the question of the treatment response. There are a lot of statistical issues in that, and it is very difficult to separate the important from the unimportant factors.
Dr. Per Lønning: It is interesting, however, to compare Therese Sorlie's PNAS paper [Proc Natl Acad Sci USA 2003;100:841823] where she applied the same gene analysis to the Dutch van't Veer material. That article reported that the discrimination regarding relapse versus survival was much smaller between the luminal A and B class in the Dutch material compared to the Norwegian material. Now, the difference is that, in the Norwegian material, all of the patients who were ER+ had been exposed to tamoxifen, and in these samples, there was a huge difference between luminal A and luminal B. Whereas in the Dutch material none of these patients received tamoxifen, and the discrimination between luminal A and luminal B regarding relapse was much smaller. So, quite indirectly, that comparison indicates that the luminal B subtype are poor responders to tamoxifen compared to luminal A, but of course this needs to be confirmed.
Dr. Stephen Johnston: The group that would be really interesting to study, based on those data sets, is the patients treated with up-front endocrine therapy in a locally advanced breast cancer setting or a preoperative neoadjuvant setting. Because if you can really differentiate responders from nonresponders, particularly in ER+ patients, then that would help us tremendously.
Dr. Richard Santen: In your data, it looked as if expression of aromatase was a negative survival factor. We did some studies years ago that looked at this. The tumors that were highly differentiated (grade 1) had high aromatase activity. As they became less differentiated, the aromatase activity went down, but then when you got to the grade 3 very aggressive tumors, the aromatase activity went back up again. It appears that as a tumor de-differentiates it begins to secrete proteins or make proteins that weren't there previously. That observation, if one accepts it, offers two ways of thinking about the biology of aromatase. It is making estrogen, but it may also be a marker of a de-differentiated tumor. That might explain why you have poor prognosis with high aromatase activity. But it really says that you can't interpret the biology simply by the over- or underexpression of a gene; there are other regulatory factors.
Dr. Kristensen: That is very interesting to hear because we do studies on genetic polymorphisms and how they influence the observed phenotypes, and we find such patterns as well. We have a polymorphism in the 3' area of aromatase that correlates with high expression levels, which correlate with large tumor. So, you would say if this polymorphism is present, more aromatase will be produced, allowing the tumor to become locally advanced. But in a parallel study from our own lab, we see that if there is no expression of aromatase, the prognosis is very poor.
Dr. McDonnell: In the article from Dennis Sgroi [Proc Natl Acad Sci USA 2003;100:5974-9] using laser capture microdissection to identify factors that predicted for poor prognosis, two or three of the genes that popped up I had never seen before in any of these analyses. I'm wondering is it because the laser capture microdissection looks only at actual epithelial cancerous cells, whereas in these studies here with larger tumors, you are actually looking at other factors from stromal cells or adipocytes?
Dr. Kristensen: Yes, there are pros and cons for both approaches. In this case, we have to use statistical tools to dissect all these genes that are not directly involved in the process. Despite the response differences Prof. Lønning points to, the fact that we see these common classes in both treated and untreated tumors and in all different stages, that means that what we observe is some intrinsic signature. For studies of mechanism, laser capture dissection is more clear-cut because one knows what the experimental settings are. But when it comes to extrapolating these data back to phenotypes like overall survival, then it is more problematic to generalize the observations.
Dr. McDonnell: I was really amazed that when you looked at what emerged as prognostic factors from these two approaches, they were completely different. There was very little overlap between them.
Dr. Lønning: I think we should focus on how can we use microarrays to screen for functional pathways and genes related to functional pathways instead of focusing on prognostic or predictive factors. You don't get a predictive factor straight out of a microarray without incorporating a functional hypothesis.
Dr. Osborne: We are still lacking a clear picture of biology. For all we know, a lot of these genes that we think now are independent of each other might be in a similar overlapping network. We just don't know.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
. Clin Cancer Res 2004;10:4029.
This article has been cited by other articles:
![]() |
L. Lusa, L. M. McShane, J. F. Reid, L. De Cecco, F. Ambrogi, E. Biganzoli, M. Gariboldi, and M. A. Pierotti Challenges in Projecting Clustering Results Across Gene Expression Profiling Datasets J Natl Cancer Inst, November 21, 2007; 99(22): 1715 - 1723. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |