
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Imaging, Diagnosis, Prognosis |
Authors' Affiliations: 1 Institut Jules Bordet, Université Libre de Bruxelles, Brussels, Belgium; 2 International Drug Development Institute, and 3 Machine Learning Group, Université Libre de Bruxelles, Brussels, Belgium; 4 Veridex LLC, a Johnson & Johnson Company, San Diego, California; 5 European Institute of Oncology and University of Milan School of Medicine, Milan, Italy; 6 National Center of Competence in Research Molecular Oncology, Swiss Institute of Experimental Cancer Research, Epalinges and the Swiss Institute of Bioinformatics, Lausanne, Switzerland; 7 Institut Gustave Roussy, Villejuif, France; 8 Karolinska Institute, Stockholm, Sweden; 9 Centre René Huguenin and Institut National de la Santé et de la Recherche Médicale, Saint-Cloud, France; 10 Guy's Hospital, London, United Kingdom; 11 John Radcliffe Hospital, Oxford, United Kingdom; and 12 Erasmus Medical Center/Daniel den Hoed Cancer Center, Rotterdam, the Netherlands
Requests for reprints: Breast International Group/TRANSBIG Secretariat, Institut Jules Bordet, 121 Boulevard de Waterloo, 1000 Brussels, Belgium. Phone: 32-2-5413526; Fax: 32-2-5413199; E-mail: TRANSBIG{at}bordet.be.
| Abstract |
|---|
|
|
|---|
Experimental Design: Gene expression profiling of frozen samples from 198 N systemically untreated patients was done at the Bordet Institute, blinded to clinical data and independent of Veridex. Genomic risk was defined by Veridex, blinded to clinical data. Survival analyses, done by an independent statistician, were done with the genomic risk and adjusted for the clinical risk, defined by Adjuvant! Online.
Results: The actual 5- and 10-year time to distant metastasis were 98% (88-100%) and 94% (83-98%), respectively, for the good profile group and 76% (68-82%) and 73% (65-79%), respectively, for the poor profile group. The actual 5- and 10-year overall survival were 98% (88-100%) and 87% (73-94%), respectively, for the good profile group and 84% (77-89%) and 72% (63-78%), respectively, for the poor profile group. We observed a strong time dependence of this signature, leading to an adjusted hazard ratio of 13.58 (1.85-99.63) and 8.20 (1.10-60.90) at 5 years and 5.11 (1.57-16.67) and 2.55 (1.07-6.10) at 10 years for time to distant metastasis and overall survival, respectively.
Conclusion: This independent validation confirmed the performance of the 76-gene signature and adds to the growing evidence that gene expression signatures are of clinical relevance, especially for identifying patients at high risk of early distant metastases.
A feature common to both signatures is that when they are compared with the classification results with two conventional consensus, the St. Gallen and NIH guidelines (5, 6), both signatures correctly identified the low-risk patients not needing treatment. This means that the number of patients that would receive unnecessary treatment could be significantly reduced, emphasizing the advantage of these gene signatures over the clinical guidelines.
TRANSBIG, a network for translational research established by the Breast International Group (BIG), recently conducted a validation study of the 70-gene signature and showed reproducible prognostic value in a series of 302 patients from five different centers and across different statistical facilities (7). In the present study, we report the validation results of the 76-gene signature using the same approach and patient series, but an independent laboratory. In addition, given the long follow-up period (median, 14 years) in this study, we also investigated the prognostic value of this 76-gene signature over time and intriguingly observed a strong time dependence of the prognostic gene signature.
| Materials and Methods |
|---|
|
|
|---|
5 cm) tumors. Patients in this series had been diagnosed between 1980 and 1998 (median follow-up, 13.6 years) and had been seen at six centers: Institut Gustave Roussy, Villejuif, France (IGR); Karolinska Institute, Stockholm and Uppsala University Hospital, Uppsala, Sweden (KI); Centre René Huguenin, Saint-Cloud, France (CRH); Guy's Hospital, London, United Kingdom (GH); and John Radcliffe Hospital, Oxford, United Kingdom (JRH). Patients with previous malignancies (except basal cell carcinoma) and bilateral synchronous breast tumors were excluded. The corresponding paraffin-embedded tumor samples of these patients were sent to the Department of Pathology at the European Institute of Oncology, Milan, Italy, where ER status (using immunohistochemistry) and histologic grade (using the Elston and Ellis method) were determined by the same pathologist, blinded to the clinical and genomic data. The clinical centers were also visited by two independent auditors who carried out source data verification of all data in the validation series. The validation protocol was finalized in September 2005, and all institutional ethics committees approved the use of the tumor materials for the purposes described in this study. Gene expression analysis. Frozen samples were sent from the clinical centers to the Netherlands Cancer Institute, Amsterdam, where RNA had been extracted for the previous TRANSBIG validation study of the 70-gene profile on the Agilent platform (7), except for the samples from Centre René Huguenin, for which RNA was sent directly to Amsterdam. To carry out the present study, RNAs were then shipped to the Jules Bordet Institute to perform the microarray analyses using the Affymetrix U133a GeneChip blinded to clinical data and independent of Veridex. The quality of the RNA obtained from each tumor sample was assessed via the RNA profile generated by the Agilent bioanalyzer. RNA amplification, hybridization, and image scanning were done according to standard Affymetrix protocols. Expression values for each gene were calculated using Affymetrix GeneChip analysis software MAS 5.0. For chip normalization, probe sets were scaled to a target intensity of 600, and scale mask filters were not selected. Chips with signal to noise ratio < 24 were excluded. Each probe set was considered as a separate gene. Genomic high- and low-risk groups were defined by Veridex blinded to clinical data, as described previously (3, 4), using the array based ER assay results and the 76-gene prognostic signature. The raw and normalized gene expression data, together with the patient's characteristics are publicly available on GEO (http://www.ncbi.nlm.nih.gov/geo), with accession number GSE 7390.
Statistical analyses. All statistical analyses were carried out by the International Drug Development Institute, Brussels, Belgium, using SAS version 9.1 and SPLUS version 7. Patient clinical data and outcomes were blinded to Jules Bordet Institute and Veridex. The end points considered in this study were time from diagnosis to distant metastases (TDM), which was the end point used to identify the gene signature (3), and overall survival, defined as time from diagnosis to death from any cause.
The main analytic approach used to validate the gene signature was to estimate hazard ratios (HR), which quantified the relative risk of an event in the high-risk group compared with the low-risk group. HRs were estimated through Cox's proportional hazard regression models, stratified by clinical center to account for the possible heterogeneity in patient selection or other potential confounders among the various centers. HRs for the risk groups defined by the gene signature were estimated with stratification for clinical risk, using the Adjuvant! Online software13 to reflect the prognostic effect of the gene signature over and above that of clinicopathologic variables ("adjusted HRs"). This clinical risk classification, which calculates a 10-year survival probability based on a patient's age, tumor size and grade, and ER status (and nodal status, which was negative for all patients considered in the present validation; ref. 8) and had been shown to predict the overall survival, breast cancerspecific survival, and event-free survival in 4,083 women diagnosed with early breast cancer in British Columbia between 1989 and 1993 (9), was also used in the TRANSBIG validation of the Amsterdam 70-gene signature (7). It was decided by consensus of the TRANSBIG consortium members14 that the low clinical risk group would be defined as patients with a 10-year survival probability of at least 88% if they had >10% of the cells with positive expression of ERs and of at least 92% if they had not (as calculated by Adjuvant! Online). These two cutoffs were chosen to reflect the fact that nowadays, patients with ER-positive tumors receive adjuvant endocrine therapy (with an estimated absolute 10-year survival benefit of about 4% overall), whereas patients in the validation series were all untreated regardless of their receptor status. The adjusted HRs in the different centers were displayed with their 95% confidence intervals on forest plots and tested for heterogeneity. The Kaplan-Meier product limit estimator was used to display time to event curves. "Sensitivity" was defined as the percentage of the patients with distant metastasis within 5 or 10 years that had been predicted correctly by the 76-gene signature. "Specificity" was defined as the percentage of patients free of distant metastasis for at least 5 or 10 years that had been predicted correctly by the 76-gene signature. The effect of the duration of follow-up on HRs was analyzed by censoring all observations at increasing time points. Cox proportional hazard models were fit on the censored data and stratified by center. HRs for the risk groups defined by the gene signature were adjusted for clinical risk by stratifying by the clinical risk groups. HRs for the risk groups defined by Adjuvant! Online were stratified only by center. To be able to compare the results with those of Adjuvant! Online and those of the recently published validation study of the 76-gene signature (4), the data reported in this study were based on the 10-year follow-up for all the analyses, except for the analysis regarding time dependence.
| Results |
|---|
|
|
|---|
|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
The HRs were equal to 5.78 for TDM and 2.87 for overall survival, and these remained similar after adjustment for the clinical risk, as assessed by Adjuvant! Online. This means that the signature adds prognostic information to that provided by the clinical risk classification. Furthermore, as reported for the initial population (10), the 76-gene signature was able to identify in this cohort a significant percentage of ER-negative tumors in the low genomic risk group (25%), whereas all ER-negative tumors were considered to be at high clinical risk according to Adjuvant! Online. This is in contrast with the Amsterdam 70-gene signature, which identified only 4% of ER-negative patients (5 of 111) in the low genomic risk group but 44% (85 of 191) in the high genetic risk group (7). The independence of the ER status of the 76-gene signature, together with its independence of age and tumor size, as well as its current validation in an independent series of patients support its strong prognostic power and emphasize its broad clinical applicability.
Given the exceptionally long median follow-up time of this validation cohort (i.e., 14 years), as opposed to a median follow-up of about 8 years for the previously reported populations, we decided to further investigate the behavior of the genomic risk over time. Interestingly, we found that the HR for the gene signature adjusted to the clinical risk was the highest when data were censored at 5 years and decreased as censoring time points increased (11). Similar results were recently also reported for the validation of the 70-gene signature in the same validation series (7), although the latter included 104 additional patients. This finding highlights that both signatures are a strong predictor of the development of early distant metastases occurring within 5 years, but that their prognostic ability decreases with increasing follow-up years. Nevertheless, the adjusted HR still remained high (5.11) in our study at 10 years of follow-up. These findings, which were not observed for the clinical risk classification, might not be unexpected because the signatures were built to identify patients with distant metastases within 5 years. When these bell-shaped curves of the HR during long-term follow-up will be confirmed in separate series of patients, it does suggest that different mechanisms may be associated with the development of early and late distant metastases. Klein et al. recently provided new insight into the evolution and progression of breast cancer (12, 13). Their results suggested that during breast cancer tumorigenesis some tumor cells may already disseminate in a much earlier genomic state and remain "dormant" until the microenvironment becomes favorable and advantageous aberrations take place in order for distant metastases to develop. However, once metastases start growing, they will still show the same global genomic phenotype as the primary tumor, as reported by several groups (1418). Indeed, these dormant cells still have the potential to develop more or less the same genetic aberrations and therefore to show a genetic phenotype similar to what is seen in the primary tumor. Thus, according to Klein's findings and to the current clinical observation that early-stage breast cancer patients still relapse despite complete resection of the primary tumor, we assume that cells may disseminate at different stages of breast cancer tumorigenesis and primary tumor development. Cells that disseminated before primary tumor development would still need to go through some genetic changes and receive the right signals from the local microenvironment before becoming fully active, and this phenomenon could explain the occurrence of late metastases. In contrast, cells that disseminated from the primary tumor in a much more advanced genomic state may already have the genomic background needed for metastatic development and could be responsible for the earlier metastases.
This means that by profiling only primary tumor samples, we might essentially be capturing information regarding tumor cells that disseminated from the primary tumor and not necessarily the young disseminated tumor cells that are in a far less progressed genomic state.
Altogether, these findings may have some important clinical implications. Indeed, in this context, adjuvant therapy might help only to eradicate "active" disseminated breast tumor cells, those which have already acquired the full capacity of unrestrained growth and which are therefore responsible for the development of early metastases. Currently administered adjuvant therapies might not be able to kill the "dormant" disseminated cells, suggesting that new targeted therapies should be developed while late adjuvant chemotherapy or switching long-term endocrine therapy (up to more than 10 years) might be considered.
In conclusion, this independent multicenter validation shows the strong reproducibility and robustness of the 76-gene expression signature and adds to the growing evidence that this classifier is of wide clinical relevance, especially to identify patients at high risk of developing early distant metastases.
| Acknowledgments |
|---|
| Footnotes |
|---|
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.
Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).
C. Desmedt and F. Piette are co-first authors. M. Buyse and C. Sotiriou are co-last authors.
All authors have contributed to this work.
These results were presented as a plenary oral presentation at the European Breast Cancer Conference, March 2006, Nice, France.
Conflicts of Interest: Y. Wang and Y. Zhang are employed by Veridex LLC, a Johnson & Johnson company, which is in the business of commercializing diagnostic products.
13 http://www.adjuvantonline.com ![]()
14 http://www.breastinternationalgroup.org ![]()
Received 11/22/06; revised 3/ 2/07; accepted 3/15/07.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. Haibe-Kains, C. Desmedt, C. Sotiriou, and G. Bontempi A comparative study of survival models for breast cancer prognostication based on microarray data: does a single gene beat them all? Bioinformatics, October 1, 2008; 24(19): 2200 - 2208. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Desmedt, B. Haibe-Kains, P. Wirapati, M. Buyse, D. Larsimont, G. Bontempi, M. Delorenzi, M. Piccart, and C. Sotiriou Biological Processes Associated with Breast Cancer Clinical Outcome Depend on the Molecular Subtypes Clin. Cancer Res., August 15, 2008; 14(16): 5158 - 5165. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Schmidt, D. Bohm, C. von Torne, E. Steiner, A. Puhl, H. Pilch, H.-A. Lehr, J. G. Hengstler, H. Kolbl, and M. Gehrmann The Humoral Immune System Has a Key Prognostic Impact in Node-Negative Breast Cancer Cancer Res., July 1, 2008; 68(13): 5405 - 5413. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. V. Glinsky "Stemness" Genomics Law Governs Clinical Behavior of Human Cancer: Implications for Decision Making in Disease Management J. Clin. Oncol., June 10, 2008; 26(17): 2846 - 2853. [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] |
||||
![]() |
L. Pusztai Current Status of Prognostic Profiling in Breast Cancer Oncologist, April 1, 2008; 13(4): 350 - 360. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Cardoso, L. Van't Veer, E. Rutgers, S. Loi, S. Mook, and M. J. Piccart-Gebhart Clinical Application of the 70-Gene Profile: The MINDACT Trial J. Clin. Oncol., February 10, 2008; 26(5): 729 - 735. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Ignatiadis, N. Xenidis, M. Perraki, S. Apostolaki, E. Politaki, M. Kafousi, E. N. Stathopoulos, A. Stathopoulou, E. Lianidou, G. Chlouverakis, et al. Different Prognostic Value of Cytokeratin-19 mRNA Positive Circulating Tumor Cells According to Estrogen Receptor and HER2 Status in Early-Stage Breast Cancer J. Clin. Oncol., November 20, 2007; 25(33): 5194 - 5202. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Harris, H. Fritsche, R. Mennel, L. Norton, P. Ravdin, S. Taube, M. R. Somerfield, D. F. Hayes, and R. C. Bast Jr American Society of Clinical Oncology 2007 Update of Recommendations for the Use of Tumor Markers in Breast Cancer J. Clin. Oncol., November 20, 2007; 25(33): 5287 - 5312. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Loi, B. Haibe-Kains, C. Desmedt, and C. Sotiriou In Reply J. Clin. Oncol., September 20, 2007; 25(27): 4317 - 4318. [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 |