
[Cancer Research 65, 9155-9158, October 15, 2005]
© 2005 American Association for Cancer Research
Molecular Portraits and 70-Gene Prognosis Signature Are Preserved throughout the Metastatic Process of Breast Cancer
Britta Weigelt1,
Zhiyuan Hu4,5,
Xiaping He4,5,
Chad Livasy6,
Lisa A. Carey7,
Matthew G. Ewend7,
Annuska M. Glas4,
Charles M. Perou4,5,6 and
Laura J. van't Veer1,2,3
Divisions of 1 Experimental Therapy and 2 Diagnostic Oncology, The Netherlands Cancer Institute; 3 Agendia B.V., Amsterdam, the Netherlands; 4 Lineberger Comprehensive Cancer Center; and Departments of 5 Genetics, 6 Pathology and Laboratory Medicine, and 7 Medicine, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina
Requests for reprints: Laura J. van't Veer, Plesmanlaan 121, 1066 CX Amsterdam, the Netherlands. Phone: 31-20-512-2754; Fax: 31-20-512-2759; E-mail: l.vt.veer{at}nki.nl.
 |
Abstract
|
|---|
Microarray analysis has been shown to improve risk stratification of breast cancer. Breast tumors analyzed by hierarchical clustering of expression patterns of "intrinsic" genes have been reported to subdivide into at least four molecular subtypes that are associated with distinct patient outcomes. Using a supervised method, a 70-gene expression profile has been identified that predicts the later appearance or absence of clinical metastasis in young breast cancer patients. Here, we show that distant metastases display both the same molecular breast cancer subtype as well as the 70-gene prognosis signature as their primary tumors. Our results suggest that the capacity to metastasize is an inherent feature of most breast cancers. Furthermore, our data imply that poor prognosis breast carcinomas classified either by the intrinsic gene set or the 70 prognosis genes represent distinct disease entities that seem sustained throughout the metastatic process.
 |
Introduction
|
|---|
DNA microarray technology, which allows the analysis of the expression levels of a thousand of genes in a single experiment, offers great potential to improve our knowledge of tumor molecular biology but also for the discovery of new molecular markers. Over the past years, microarray analysis has been extensively used to improve the diagnosis and risk stratification of many cancers (16). Two major studies have described the use of this technology to assess the molecular classification of human breast cancer and have defined new subgroups that are relevant to patient management. Using hierarchical clustering, Perou et al. showed that breast tumors can be classified into specific subtypes based solely on differences in gene expression patterns named the molecular portraits of breast tumors (3). Three estrogen receptor (ER)negative subtypes of breast carcinomas were identified ("basal like," "HER2+," and "normal breast like"), and later, at least two ER-positive tumor subtypes ("luminal A" and "luminal B") were defined (4). Importantly, the breast cancer subtypes also represent clinically distinct subgroups of patients, as they show differences in metastasis-free and overall survival (4). Sørlie et al. identified the luminal A subgroup of ER-positive tumors to be associated with the best outcome, whereas the basal-like and the HER2+ tumors have the worst outcome (4). The identification of the distinct subtypes was based upon the hierarchical clustering of an "intrinsic" gene set, which comprises genes that show little variance within repeated samplings of the same tumor but which show high variance across different tumors (3). Recently, a new intrinsic gene set of 1,300 genes (i.e., the "Intrinsic/UNC" gene set) was identified using a larger set of breast tumors and genes than in the original report and was validated on a test set of tumors from independent microarray studies.8 The breast tumor subtypes based on the new Intrinsic/UNC gene set were nearly identical to those previously identified and shown to be reproducible across independent data sets, across different microarray platforms, also with respect to clinical outcome.8 In complementary studies, van't Veer et al. used a supervised classification method to identify a gene expression signature in young breast cancer patients with lymph nodenegative tumors that is associated with patient outcomes and prognosis (6). Using this signature consisting of 70 genes, primary breast tumors can be classified as having either a "poor prognosis" signature, which means they are likely to metastasize, or a "good prognosis" signature, meaning that the development of metastases is unlikely. This 70-gene signature was subsequently confirmed in a consecutive series of 295 patients with both lymph nodepositive and lymph nodenegative disease (7). The finding that the metastasis risk of a breast cancer patient can be predicted by the overall gene expression profile of its primary tumor challenged the widely accepted idea that metastatic potential is acquired relatively late during multistep tumorigenesis (8). The concept that the ability to metastasize to distant sites is an early and inherent genetic property of breast cancers seemed supported by our previous finding that the overall molecular program established in a human primary breast carcinoma is preserved in its distant metastasis (9). Here, we test whether the molecular subtypes of breast cancer as well as the 70-gene prognosis profiles are also maintained throughout the metastatic process.
 |
Materials and Methods
|
|---|
Detailed information on RNA isolation, amplification, labeling, hybridization, scanning, microarray analysis, and patient information has been described previously (6, 7, 9).8 The microarray data presented here have been deposited into the GEO under the series number of GSE2741, with the 23 samples not described in Hu et al. being GSM52910 to GSM52932.
In brief, seven pairs of matching primary breast tumors and distant metastases of the Netherlands Cancer Institute described previously (ref. 9; pairs 1 and 3-8) were hybridized and analyzed for their 70-gene expression profile (6), and five pairs (pairs 3 and 5-8) analyzed by hierarchical clustering for their molecular breast cancer subtype. The hierarchical clustering analysis was done using the Intrinsic/UNC gene list comprising 1,410 microarray elements (representing 1,300 genes)8 using 156 arrays [Agilent Human 1A Oligo Microarray (V2)] representing 107 patients. The molecular subtypes of the samples were determined by the dendrogram that was associated with characteristic gene expression patterns. In addition, a set of five pairs of primary tumors and lymph node metastases, one pair of primary breast tumor and a brain metastasis, as well as multiple distant metastasis samples of five autopsy patients from the University of North Carolina (UNC) at Chapel Hill were analyzed for their molecular breast cancer subtype.
These studies were approved by the Medical Ethical Committee of the Netherlands Cancer Institute and the Institutional Review Board of UNC.
 |
Results
|
|---|
To test whether the molecular subtype of a primary breast tumor is preserved in its metastasis, five of the seven pairs of primary tumors and matching metastases described in Weigelt et al. (9) were retested and analyzed using hierarchical clustering and 107 additional breast tumors using a new breast "intrinsic gene list."8 All primary breast tumors paired with their matching metastases, even when added to the large data set of breast tumors from 107 patients (Fig. 1, red codes). Two matching pairs were classified into the HER2+ group, one pair into the basal-like group, and one pair was "ER-negative unclassified" and clustered near the basal-like and the HER2+ group. These ER-negative breast cancer subtypes are associated with the shortest relapse-free and overall survival times (4).8 Pair 8 was classified into the luminal group of breast cancer subtypes. The finding that the individual portraits of tumors are maintained in their metastases was further confirmed by the paired clustering of an independent set of five pairs of primary tumors and simultaneous lymph node metastases (Fig. 1, light blue codes) and one primary and metastatic brain tumor pair (Fig. 1, red code) of UNC. In a few additional cases from autopsy patients, we were able to sample and compare multiple metastasis sites from the same individual. The primary tumor and metastases of the spinal cord, liver, adrenal gland, lymph node, and lung of autopsy patient A1 were analyzed for their molecular subtype; of autopsy patient A7, metastases of the liver, kidney, lung, lymph node, diaphragm, and brain were tested. Remarkably, all metastasis samples obtained from one breast cancer patient cluster together and show the same molecular breast cancer subtype (Fig. 1, pink codes).

View larger version (136K):
[in this window]
[in a new window]
|
Figure 1. Determination of the molecular intrinsic subtypes of primary tumors and their metastatic counterparts. A hierarchical clustering analysis was performed using a 1,300-gene intrinsic gene list.8 Row, tumor; column, gene; red codes, six primary tumor-distant metastasis pairs; light blue codes, primary tumor-lymph node metastasis pairs; pink codes, local and distant metastasis samples from autopsy patients (which for patient A1 includes the primary).
|
|
In the next step, the maintenance of the 70-gene prognostic signature throughout the metastatic process was tested for seven pairs of primary tumor and matching metastases (6, 9). Five primary breast carcinomas had a 70-gene signature associated with poor prognosis, as did their distant metastases (Fig. 2). Interestingly, pair 8 that was classified into the good prognosis luminal group of breast tumors also has a good prognosis 70-gene signature in both its primary and recurrence that was separated in time by 15 years. Only the primary tumor of pair 4 has a different 70-gene signature than its metastasis.

View larger version (58K):
[in this window]
[in a new window]
|
Figure 2. Expression data matrix of 70 prognostic marker genes from 78 primary breast tumors plus seven pairs of matching pairs of primary tumors and distant metastases (6, 9). Row, tumor; column, gene. Genes are ordered according to their correlation coefficient with the two prognostic groups; tumors are ordered by the correlation to the average profile of the good prognosis group (right). Above the yellow line, patients have a good prognosis signature, below the line, a poor prognosis signature.
|
|
 |
Discussion
|
|---|
Gene expression profiling confirms that breast cancer is a heterogeneous disease, both biologically and clinically (3, 4, 6, 10, 11).8 Nevertheless, the genome-wide expression analysis of breast cancers has made it possible to identify signature patterns in primary breast carcinomas that are associated with patient outcomes and prognoses. Our data show that distant metastases mirror the specific prognostic profiles, the molecular breast cancer subtypes and 70-gene prognostic signatures, of their primary breast tumors. Remarkably, multiple metastases from one patient all display the same molecular breast cancer subtype independent of the organ in which they developed and still maintain the unique molecular identity of the primary that they arose from. Our findings support the hypothesis that the molecular subtypes might originate from different cell types within the breast and therefore reflect different biological entities (12), which are maintained throughout the multistep metastatic process.
Almost all pairs are either of the poor prognosis molecular subtype HER2+ or basal-like or have a poor prognosis signature of 70 genes, which is what might be expected because all tumors metastasized. Primary tumor 8 and its distant metastasis sample in the ovary, however, had a good prognosis signature of 70 genes and a luminal breast cancer subtype, both of which are associated with a low metastasis risk. Of note, the metastasis of this primary tumor developed 15 years after primary diagnosis, compared with a mean of 4.5 years (range, 1.6-6.3) in the other six patients. Only one primary breast carcinoma did not maintain the 70-gene signature in its distant metastasis (pair 4), although these two samples were relatively close as observed by their Pearson correlation coefficient. This pair was due to limited amounts of RNA not tested for its breast cancer subtype. For the samples analyzed from UNC, all primary tumors and their associated metastases were collected at the same time.
The analysis of the 70-gene signature revealed that in six of seven pairs, the metastasis has a higher correlation with the poor prognosis signature than its primary tumor (see Fig. 2). Genes correlated with the good prognosis signature of 70 genes remain virtually unchanged between primary and metastatic tumors within a patient. A small number of genes correlated with the poor prognosis signature is, however, up-regulated in the metastases compared with their matching primary breast carcinomas (data not shown), which leads to the "poorer" 70-gene signature of the metastases. The genes up-regulated in the metastatic tumors are among others involved in DNA replication (RFC4 and ORC6L) and signal transduction (IGFBP5 and PRC1; data not shown). Interestingly, only matrix metalloproteinase 9, which plays an important role in the proteolysis of the extracellular matrix (13), is down-regulated in five of the six metastases showing a poorer profile than their primary tumors (data not shown).
Our results presented here emphasize that the metastatic nature of poor prognosis breast carcinomas, which are depicted by the poor prognosis 70-gene profile or the luminal B, HER2+, or basal-like molecular subtype, is an inherent feature of breast cancers that remains stable with time and across distinct tumor outgrowth locations within the same individual. Because both the molecular breast cancer subtype and prognostic expression profile of a primary breast tumor are maintained throughout the metastatic process, future treatment decisions based on the expression profile of a primary tumor is a rational approach towards preventing the outgrowth of metastases.
 |
Acknowledgments
|
|---|
Grant support: Dutch Cancer Society (B. Weigelt and L.J. van't Veer), Cancer Genomics Initiative (B. Weigelt and L.J. van't Veer), National Cancer Institute Breast Specialized Programs of Research Excellence/UNC at Chapel Hill grant P50-CA58223-09A1 (C.M. Perou), National Institute of Environmental Health Sciences grant U19-ES11391-03 (C.M. Perou), National Cancer Institute grant RO1-CA-101227-01 (C.M. Perou), and NIH General Clinical Research Centers grant M01RR00046 (L.A. Carey).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank Drs. Phil Bernard, Juan Palazzo, and Olufunmilayo I. Olopade for their contributions towards the 107-patient sample data set.
 |
Footnotes
|
|---|
8 Z. Hu et al. A new breast tumor intrinsic gene list identifies novel characteristics that are conserved across microarray platforms, submitted for publication. 
Received 7/20/05.
Revised 8/23/05.
Accepted 8/26/05.
 |
References
|
|---|
- Golub TR, Slonim DK, Tamayo P, et al. Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science 1999;286:5317.[Abstract/Free Full Text]
- Alizadeh AA, Eisen MB, Davis RE, et al. Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature 2000;403:50311.[CrossRef][Medline]
- Perou CM, Sorlie T, Eisen MB, et al. Molecular portraits of human breast tumours. Nature 2000;406:74752.[CrossRef][Medline]
- Sorlie T, Perou CM, Tibshirani R, et al. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci U S A 2001;98:1086974.[Abstract/Free Full Text]
- Pomeroy SL, Tamayo P, Gaasenbeek M, et al. Prediction of central nervous system embryonal tumour outcome based on gene expression. Nature 2002;415:43642.[CrossRef][Medline]
- van 't Veer LJ, Dai H, van de Vijver MJ, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature 2002;415:5306.[CrossRef][Medline]
- van de Vijver MJ, He YD, van't Veer LJ, et al. A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med 2002;347:19992009.[Abstract/Free Full Text]
- Fidler IJ, Kripke ML. Metastasis results from preexisting variant cells within a malignant tumor. Science 1977;197:8935.[Abstract/Free Full Text]
- Weigelt B, Glas AM, Wessels LF, Witteveen AT, Peterse JL, van't Veer LJ. Gene expression profiles of primary breast tumors maintained in distant metastases. Proc Natl Acad Sci U S A 2003;100:159015.[Abstract/Free Full Text]
- Gruvberger S, Ringner M, Chen Y, et al. Estrogen receptor status in breast cancer is associated with remarkably distinct gene expression patterns. Cancer Res 2001;61:597984.[Abstract/Free Full Text]
- West M, Blanchette C, Dressman H, et al. Predicting the clinical status of human breast cancer by using gene expression profiles. Proc Natl Acad Sci U S A 2001;98:114627.[Abstract/Free Full Text]
- Sorlie T. Molecular portraits of breast cancer: tumour subtypes as distinct disease entities. Eur J Cancer 2004;40:266775.
- Bodey B, Bodey B, Jr., Siegel SE, Kaiser HE. Matrix metalloproteinases in neoplasm-induced extracellular matrix remodeling in breast carcinomas. Anticancer Res 2001;21:20218.[Medline]
This article has been cited by other articles:

|
 |

|
 |
 
P. S. Steeg
Heterogeneity of Drug Target Expression Among Metastatic Lesions: Lessons from a Breast Cancer Autopsy Program
Clin. Cancer Res.,
June 15, 2008;
14(12):
3643 - 3645.
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
L. Xu, S. S. Shen, Y. Hoshida, A. Subramanian, K. Ross, J.-P. Brunet, S. N. Wagner, S. Ramaswamy, J. P. Mesirov, and R. O. Hynes
Gene Expression Changes in an Animal Melanoma Model Correlate with Aggressiveness of Human Melanoma Metastases
Mol. Cancer Res.,
May 1, 2008;
6(5):
760 - 769.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
J. S. Ross, C. Hatzis, W. F. Symmans, L. Pusztai, and G. N. Hortobagyi
Commercialized Multigene Predictors of Clinical Outcome for Breast Cancer
Oncologist,
May 1, 2008;
13(5):
477 - 493.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
C. E. Caldon, A. Swarbrick, C. S.L. Lee, R. L. Sutherland, and E. A. Musgrove
The Helix-Loop-Helix Protein Id1 Requires Cyclin D1 to Promote the Proliferation of Mammary Epithelial Cell Acini
Cancer Res.,
April 15, 2008;
68(8):
3026 - 3036.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
I. G. Yulug and B. Gur-Dedeoglu
Functional genomics in translational cancer research: focus on breast cancer
Brief Funct Genomic Proteomic,
March 7, 2008;
(2008)
eln009v1.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
D. J. Wong, D. S.A. Nuyten, A. Regev, M. Lin, A. S. Adler, E. Segal, M. J. van de Vijver, and H. Y. Chang
Revealing Targeted Therapy for Human Cancer by Gene Module Maps
Cancer Res.,
January 15, 2008;
68(2):
369 - 378.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
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]
|
 |
|

|
 |

|
 |
 
J. C. Harrell, W. W. Dye, D. M.E. Harvell, M. Pinto, P. Jedlicka, C. A. Sartorius, and K. B. Horwitz
Estrogen Insensitivity in a Model of Estrogen Receptor Positive Breast Cancer Lymph Node Metastasis
Cancer Res.,
November 1, 2007;
67(21):
10582 - 10591.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
B. Schade, S. H.L. Lam, D. Cernea, V. Sanguin-Gendreau, R. D. Cardiff, B. L. Jung, M. Hallett, and W. J. Muller
Distinct ErbB-2 Coupled Signaling Pathways Promote Mammary Tumors with Unique Pathologic and Transcriptional Profiles
Cancer Res.,
August 15, 2007;
67(16):
7579 - 7588.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
R. D. Loberg, D. A. Bradley, S. A. Tomlins, A. M. Chinnaiyan, and K. J. Pienta
The Lethal Phenotype of Cancer: The Molecular Basis of Death Due to Malignancy
CA Cancer J Clin,
July 1, 2007;
57(4):
225 - 241.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
W. W.-L. Wong, J. W. Clendening, A. Martirosyan, P. C. Boutros, C. Bros, F. Khosravi, I. Jurisica, A. K. Stewart, P. L. Bergsagel, and L. Z. Penn
Determinants of sensitivity to lovastatin-induced apoptosis in multiple myeloma
Mol. Cancer Ther.,
June 1, 2007;
6(6):
1886 - 1897.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
C. Desmedt, F. Piette, S. Loi, Y. Wang, F. Lallemand, B. Haibe-Kains, G. Viale, M. Delorenzi, Y. Zhang, M. S. d'Assignies, et al.
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.,
June 1, 2007;
13(11):
3207 - 3214.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
P. Stafford and M. Brun
Three methods for optimization of cross-laboratory and cross-platform microarray expression data
Nucleic Acids Res.,
May 11, 2007;
35(10):
e72 - e72.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. Nicolau, R. Tibshirani, A.-L. Borresen-Dale, and S. S. Jeffrey
Disease-specific genomic analysis: identifying the signature of pathologic biology
Bioinformatics,
April 15, 2007;
23(8):
957 - 965.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
E. I. Chen, J. Hewel, J. S. Krueger, C. Tiraby, M. R. Weber, A. Kralli, K. Becker, J. R. Yates III, and B. Felding-Habermann
Adaptation of Energy Metabolism in Breast Cancer Brain Metastases
Cancer Res.,
February 15, 2007;
67(4):
1472 - 1486.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
R. J. Pelham, L. Rodgers, I. Hall, R. Lucito, K. C. Q. Nguyen, N. Navin, J. Hicks, D. Mu, S. Powers, M. Wigler, et al.
Identification of alterations in DNA copy number in host stromal cells during tumor progression
PNAS,
December 26, 2006;
103(52):
19848 - 19853.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
B. G. Haffty, Q. Yang, M. Reiss, T. Kearney, S. A. Higgins, J. Weidhaas, L. Harris, W. Hait, and D. Toppmeyer
Locoregional Relapse and Distant Metastasis in Conservatively Managed Triple Negative Early-Stage Breast Cancer
J. Clin. Oncol.,
December 20, 2006;
24(36):
5652 - 5657.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
S. K Gruvberger-Saal, H. E Cunliffe, K. M Carr, and I. A Hedenfalk
Microarrays in breast cancer research and clinical practice - the future lies ahead
Endocr. Relat. Cancer,
December 1, 2006;
13(4):
1017 - 1031.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. Lacroix
Significance, detection and markers of disseminated breast cancer cells
Endocr. Relat. Cancer,
December 1, 2006;
13(4):
1033 - 1067.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
P. Roepman, A. de Jager, M. J.A. Groot Koerkamp, J. A. Kummer, P. J. Slootweg, and F. C.P. Holstege
Maintenance of Head and Neck Tumor Gene Expression Profiles upon Lymph Node Metastasis
Cancer Res.,
December 1, 2006;
66(23):
11110 - 11114.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. S. Wicha
Cancer stem cells and metastasis: lethal seeds.
Clin. Cancer Res.,
October 1, 2006;
12(19):
5606 - 5607.
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. V. Fournier, K. J. Martin, P. A. Kenny, K. Xhaja, I. Bosch, P. Yaswen, and M. J. Bissell
Gene expression signature in organized and growth-arrested mammary acini predicts good outcome in breast cancer.
Cancer Res.,
July 15, 2006;
66(14):
7095 - 7102.
[Abstract]
[Full Text]
[PDF]
|
 |
|