- Methodology article
Improvement of the clinical applicability of the Genomic Grade Index through a qRT-PCR test performed on frozen and formalin-fixed paraffin-embedded tissues
BMC Genomicsvolume 10, Article number: 424 (2009)
Proliferation and tumor differentiation captured by the genomic grade index (GGI) are important prognostic indicators in breast cancer (BC) especially for the estrogen receptor positive (ER+) disease. The aims of this study were to convert this microarray index to a qRT-PCR assay (PCR-GGI), which could be realized on formalin fixed paraffin embedded samples (FFPE), and to assess its prognostic performance and predictive value of clinical benefit in early and advanced ER+ BC patients treated with tamoxifen.
The accuracy and concordance of the PCR-GGI with the GGI was assessed using BC patients for which frozen and FFPE tissues as well as microarray data were available (n = 19). The evaluation of the prognostic value of the PCR-GGI was assessed on FFPE material using a consecutive series of 212 systemically treated early BC patients. The predictive performance for tamoxifen benefit was assessed using two ER+ BC populations treated either with adjuvant tamoxifen only (n = 77+139) or first-line tamoxifen for advanced disease (n = 270).
The PCR-GGI is based on the expression of 8 genes (4 representative of the GGI and 4 reference genes). A significant correlation was observed between the microarray-derived GGI and the qRT-PCR assay using frozen (ρ = 0.95, p < 10E-06) and FFPE material (ρ = 0.89, p < 10E-06). The prognostic performance of the PCR-GGI was confirmed on FFPE samples (HRunivar. = 1.89; [95CI:1.01-3.54], p = 0.05). The PCR-GGI further identified two subgroups of patients with statistically different time to distant metastasis free survival (DMFS) across the two cohorts of ER+ BC patients treated with adjuvant tamoxifen. Additionally, the PCR-GGI was associated with response to tamoxifen in the advanced setting (HRunivar. = 1.98; [95CI:1.51-2.59], p = 6.9E-07).
PCR-GGI recapitulates in an accurate and reproducible manner the performances of the GGI using frozen and FFPE samples.
Gene expression profiling has resulted in paradigm shift in the way that researchers and clinical investigators view breast cancer (BC) biology. During the last years, several groups have evaluated the potential of gene expression profiles to improve BC prognostication. This has resulted in the identification of several gene signatures, most of them outperforming current clinico-pathological parameters in predicting clinical outcome (reviewed in ).
One example is the genomic grade index (GGI) developed by our group . The hypothesis beyond the development of this index was to improve BC grading and its prognostic value. This index includes 97 genes that are consistently differentially expressed between low and high histological grade breast tumors. Of interest, the majority of them are related to cell cycle and proliferation. One of the major findings of that study was that GGI could assign intermediate histological grade tumors into two distinct subgroups whose gene expression profiles and clinical outcome was similar to the ones of low and high histological grade tumors respectively. In addition, GGI could identify two clinically relevant ER+ subtypes with very distinct clinical outcomes in both systemically untreated and tamoxifen only treated BC patients .
Finally, comparative studies of the GGI with other prognostic gene signatures in the context of a large meta-analysis involving ~3000 BC patients suggest that most of these signatures have a similar prognostic performance which is limited to the ER+ disease, and that the proliferation and cell-cycle genes represent the driving force of these signatures [4, 5].
In this study we aimed to transpose the GGI onto a qRT-PCR assay based on a minimal set of genes that could recapitulate in an accurate and reproducible manner the prognostic performance of GGI using both frozen and paraffin-embedded (FFPE) tumor samples, to facilitate its use in clinical practice. Furthermore, we investigated the performance of this assay to predict benefit to adjuvant tamoxifen in early BC patients or to first line tamoxifen in advanced BC patients.
Development of the reduced genomic qRT-PCR grade index (PCR-GGI)
The first step in the development of the PCR-GGI was to identify a minimum number of genes that could recapitulate the performance of the original GGI.
To this end, we first tested the correlation of the expression values of several random combinations of minimum 4 genes, covering all phases of the cell cycle, with the original GGI. This correlation was performed using gene expression data from two publicly available microarray data sets namely the NKI  and VDX . We selected a set of four genes (MYBL2, KPNA2, CDC2 and CDC20) based on their high significant correlation with the original GGI (0.95 and 0.94 (all P < 10-6) in the NKI and VDX data sets respectively). We also selected 4 reference genes (GUS, TBP, RPLPO and TFRC). Although all housekeeping genes gave similar results for normalizing the data, the variance of the estimation of the housekeeping Ct value was reduced by using the mean of these 4 reference genes (data not shown).
The prognostic performance of this reduced gene set (GGIreduced) was then compared to the performance of the original GGI in the NKI and VDX data sets. This comparison revealed similar prognostic performances (NKI data set: HRGGI = 1.73; [95%CI: 1.40-2.14], p = 4.2E-07 and HRGGIreduced = 1.49; [95%CI: 1.21-1.82], p = 1.2E-04; VDX data set: HRGGI = 1.41; [95%CI: 1.13-1.76], p = 2.0E-03 and HRGGreduced = 1.36; [95%CI:1.09-1.68], p = 5.8E-03). Of interest, the correlation values with the original GGI and the prognostic performances were not improved when additional genes were taken into consideration (data not shown).
The second step in developing the PCR-GGI was to convert these four genes into a qRT-PCR index. To this end, we designed primers for the four selected genes as well as for four reference genes considering the SYBR-based real-time PCR assay platform. We assessed the concordance between the microarray-derived GGI and the PCR-GGI measured on frozen and FFPE tissues using a series of 19 primary BC from which frozen, FFPE tissues and their corresponding microarray data were available, further referred as "IJBtest". A statistically significant correlation was observed between the expression levels of the microarray-derived GGI and the PCR-GGI assay using frozen samples (ρ = 0.95 [95%CI: 0.86-0.98], p = 3.6E-09) as well as FFPE samples (ρ = 0.89 [95%CI: 0.72-0.96], p = 8.26E-07) (see Table 1). Similarly, a statistically significant correlation was found between the expression levels of the PCR-GGI derived from frozen and FFPE samples respectively (ρ = 0.85 [95%CI: 0.64-0.94], p = 7.7E-06).
Interestingly, the correlations between the microarray-derived GGI and the PCR-GGI assay on frozen and FFPE samples were higher when the 4 genes were combined together compared to the individual genes (see table 1).
Prognostic performance of the PCR-GGI
To evaluate the prognostic performance of the PCR-GGI, we considered an independent population of 212 primary BC FFPE samples originated from patients consecutively diagnosed at the Jules Bordet Institute between 1995 and 1996 (referred as IJB95/96). The clinico-pathological characteristics are summarized in Table 2.
The prognostic performances of the individual four genes (included in the PCR-GGI) were compared to the one of the PCR-GGI and summarized in Table 3. The group of patients with high PCR-GGI scores was associated with a higher risk of recurrence than the group of patients with low PCR-GGI (HR = 1.89; [95%CI: 1.01-3.54], p = 0.05; see Figure 1a). However, this association was not significant but showed a trend towards significance in the multivariate analysis. Indeed, when selecting the age of the patient, the size of the tumor, the centrally reviewed histological grade, the nodal status and the PCR-GGI as variables for the multivariate analysis using a backward stepwise selection for the identification of the covariates to be kept in the final Cox regression model for distant metastasis free survival (DMFS), only the nodal status (HR = 3.10; [95%CI: 1.53-6.27], p = 0.002) and PCR-GGI (HR = 1.89; [95%CI: 0.94-3.77], p = 0.07) were still present in the last step of the analysis.
As expected, histological grade 2 tumors (n = 89) spanned the expression levels of PCR-GGI of histological grade 1 and 3 (data not shown). Similarly to GGI, PCR-GGI divided histological grade 2 tumors into two subgroups with distinct clinical outcome although of borderline significance (HR = 2.27; [95%CI: 0.94-5.49], p = 0.07; see Figure 1b).
Since we previously demonstrated that the GGI was particularly performant at identifying high and low risk patients within the ER-positive population , we restricted our analysis to the ER+ subpopulation of the "IJB95/96" cohort. In order to identify these ER-positive patients, we evaluated ER by qRT-PCR and defined a cut-off for positivity based on the binomial repartition of ER score (see Material and methods). Sixty-four percent (136/212) of the samples were considered as ER-positive. Overall concordance with IHC data assessed on 162 FFPE samples (76%) for which both IHC and qRT-PCR results were available was of 75% with a correlation of 0.38 (p = 9E-07). When PCR-GGI was applied to these ER+ samples, Kaplan-Meier survival curves and Cox analyses revealed a significant association between a high risk of recurrence and a high PCR-GGI score (HR = 2.26; [95%CI: 1.08-4.75], p = 0.03; see Figure 1c).
Several studies have demonstrated that the widely used proliferation marker Ki-67 evaluated by IHC predicts clinical outcome . We therefore evaluated the prognostic value of Ki-67 by IHC in the IJB95/96 series (N = 160). There was no significant association between Ki-67 and the risk of recurrence both when considering Ki-67 as a continuous variable (HR = 1.01; [95%CI: 0.99-1.03], p = 0.19) or as a binary variable using the cut-off of >15%;(HR = 1.03; [95%CI: 0.51-2.07], p = 0.93).
Assessment of the performance of the PCR-GGI in tamoxifen-treated BC patients
In order to assess a potential predictive value of the PCR-GGI for clinical benefit to tamoxifen, we evaluated its performance in three independent ER-positive BC populations treated either with adjuvant tamoxifen only or with first-line tamoxifen for advanced disease.
i) ER-positive BC patients treated with adjuvant tamoxifen only
The potential value of the PCR-GGI to predict DMFS in tamoxifen-treated patients was first evaluated in a patient's series on which the original GGI was previously computed  (n = 77), offering the possibility to compare the performance of GGI and PCR-GGI (referred as OXFD, see Table 2). Of note, 62% (45/73) of the samples originated from node-negative patients.
We observed a significant association between a high PCR-GGI score and a higher risk of recurrence (HR = 3.34; [95%CI: 1.19-9.37], p = 0.01). This performance appeared similar to the one of the original GGI (HRunivar. = 4.89; [95%CI: 1.73-13.83]; p = 0.003) and superior to the one of the histological grade (HRunivar. = 2.01; 95%CI: [0.87-4.67]; p = 0.10). The Kaplan-Meier curves for the histological grade, the PCR-GGI and the original GGI are illustrated in Figure 2. When selecting the age of the patient, the size of the tumor, the histological grade, the nodal status and the PCR-GGI as variables for the multivariate analysis using a backward stepwise selection for the identification of the covariates to be kept in the final Cox regression model for DMFS, only the tumor size (HR = 6.22; [95%CI: 1.35-28.75], p = 0.02), the age of the patient (HR = 0.27; [95%CI: 0.07-1.03], p = 0.06), and PCR-GGI (HR = 2.53; [95%CI: 0.80-7.99], p = 0.12) were still present in the last step of the analysis.
The PCR-GGI was further applied to an independent population of 139 ER-positive BC samples originated from node positive patients that received adjuvant tamoxifen for primary disease (referred as JNIadj; see Table 2). Again, we observed a statistically significant association between a high PCR-GGI score and a higher risk of recurrence (HR = 1.87; [95%CI: 1.13-3.09], p = 0.015). The Kaplan-Meier curves are illustrated in Figure 3. Interestingly, a 3-years delay of recurrence was observed between the patients with low and high PCR-GGI.
The PCR-GGI remained significantly associated with DMFS in the multivariate analysis (HR = 1.97; [95%CI: 1.09-3.55], p = 0.024) together with tumor size (HR = 2.65; [95%CI: 1.25-5.61], p = 0.011), but not with the age of the patient. Of note, the nodal status and the histological grade were not included in the analysis due to the node-positivity of all patients and the high number of missing values for the histological grade.
ii) ER+ BC patients treated with first-line tamoxifen for metastatic disease (JNImeta)
The PCR-GGI was applied to an independent population of 270 ER+ BC samples originated from patients that received first-line tamoxifen for advanced disease (referred as JNImeta; see Table 2).
A statistically significant association was observed between a high PCR-GGI score and progression-free survival (PFS) (HR = 1.98; [95CI: 1.51-2.59], p = 6.9E-07); disease progression after start of first-line tamoxifen being observed with a delay of 7.5 months in the patients with low PCR-GGI compared to the patients with high PCR-GGI at 50% PFS (see Figure 4). The PCR-GGI remained significantly associated with PFS in the multivariate analysis (HR = 1.98; [95CI: 1.51-2.59], p = 6.9E-07). We also observed a significant association between the PCR-GGI and response to tamoxifen as the patients with a low PCR-GGI score had a higher probability of response to tamoxifen (74.5% versus 53% for the patients with high PCR-GGI scores, p = 5E-04). The PCR-GGI correctly classified 81 of the 106 non-responders patients in the high grade subgroup of samples (76.5% sensitivity, 95%CI: 67-84%) and 73 of the 164 patients with objective response (OR) or stable disease lasting over six months in the low grade subgroup of samples (44.5% specificity, 95%CI: 37-52%) with an odds ratio of 2.59 (95%CI: 1.51-4.48).
In this study, we have developed a qRT-PCR grade index (PCR-GGI) composed of four genes associated with cell cycle progression and proliferation initially included in the GGI . This gene set includes MYBL2, KPNA2, CDC2 and CDC20, which together cover all phases of the cell cycle. Indeed, MYBL2 encodes a nuclear protein involved in cell cycle progression; KPNA2 is involved in the import of proteins to the nuclear envelope and acts as a regulator of cell cycle checkpoint mediators; CDC2 encodes a catalytic subunit of the highly conserved protein kinase complex known as M-phase promoting factor (MPF) which is essential for G1/S and G2/M phase transitions; and, CDC20 encodes a protein acting as a regulatory protein interacting with several other proteins at multiple points in the cell cycle.
We demonstrated that a quantitative assessment of a very small number of genes is sufficient to recapitulate the performance of the original GGI, especially in ER-positive breast cancer. This could be explained by the fact that the expression levels of most genes of the GGI signature are highly correlated with one another. Two out of the four genes included in our PCR-GGI signature were also present in a number of previously reported prognostic signatures; MYBL2 was part of the Recurrence Score developed by Paik et al.  as well as the microarray signatures defined by Naderi et al.  and Miller et al. , whereas the KPNA2 gene was present in the 76-gene signature of Wang et al. , which has been applied successfully on a cohort of node-negative estrogen receptor positive breast cancer patients to predict benefit of tamoxifen .
We further showed that the averaged expression of the four selected genes evaluated by qRT-PCR using FFPE samples accurately recapitulates the prognostic performance of the GGI. Moreover the PCR-GGI was conducted with success from 10 years old FFPE samples. Fourteen percent of the samples have been excluded from the study and this number is expected to be lower using recent FFPE samples.
Additionally, as illustrated in this manuscript, the PCR-GGI has several advantages which highlight its clinical relevance.
First, as illustrated previously for the GGI, the PCR-GGI is able to identify high and low-risk patients within the subgroup of patients with intermediate histological grade tumors, possibly improving treatment decision making for these patients.
Second, as the GGI, the PCR-GGI is able to identify a higher proportion of early breast cancer patients at low risk of recurrence than do the clinical guidelines. This means that the number of patients that would receive unnecessary treatment could be reduced by applying this molecular assay.
Third, the PCR-GGI is not subject to reproducibility and inter-observer variability problems such as the histological grading. Moreover, the PCR-GGI's procedure developed at Jules Bordet was easily transferred at the Josephine Nefkens Institute, Rotterdam following some inter-laboratory tests which gave great confidence in the reproducibility of the assessment of the PCR-GGI. Also, the PCR-GGI cut-off for positivity identified at the Jules Bordet Institute was easily applied on the Josephine Nefkens Institute assays.
Fourth, the PCR-GGI can classify ER-positive BC patients treated either with adjuvant tamoxifen only or first-line tamoxifen for advanced disease into clinically relevant subgroups. We first observed a statistically significant association between a high PCR-GGI score and a higher risk of recurrence in ER-positive BC samples from patients that received adjuvant tamoxifen for primary disease across two different BC populations. We also observed a statistically significant association between a high PCR-GGI score and a higher risk of progression in ER+ BC samples originated from patients that received first-line tamoxifen for advanced disease. However, the determination of ER positivity might also play a role in the assessment of the performance of the PCR-GGI since we report here a discordance of 25% between IHC and RT-PCR assessment of ER.
Finally, the PCR-GGI has two practical advantages: 1/ it requires only very small amounts of routinely available FFPE samples, whereas the original GGI requires fresh-frozen material, and 2/ it allows the use of a SYBR-based technology, instead of the specific Taqman probes (Applied Biosystems), providing the advantages of being easy to use and cheaper than specific probes. Therefore, we might envisage that this test could be carried out in local certified pathology laboratories.
Additionally, several studies have demonstrated that tumors with characteristics associated with poor histological grade and high proliferating index tended to respond better to chemotherapy [13, 14]. As our signature was derived from the histological grade we might then hypothesize that the PCR-GGI not only quantifies the likelihood of BC recurrence in women with node-negative ER+ BC, it might also predict the magnitude of chemotherapy benefit.
As for the GGI, the clinical potential of the PCR-GGI is obvious since: 1/ the three category of histological grading is replaced by a two-category grading clinically more relevant, 2/ it classifies ER-positive BC patients into two clinically relevant subgroups, and 3/ the genes expression score is not subject to inter-observer variability as is the histological grade. Additionally, the PCR-GGI presents the advantage that it can be evaluated on FFPE samples, which are more widely available than frozen samples.
Patients and Samples Collection
A small series of 19 primary BC from which frozen, FFPE tissues and their corresponding microarray data were available, further referred as "IJBtest", was used to compare the accuracy and the concordance of the microarray-derived GGI to the qRT-PCR genomic grade index (PCR-GGI).
The prognostic value of the PCR-GGI was evaluated in a consecutive series of 212 FFPE primary BC samples collected at the Jules Bordet Institute from 1995 to 1996, referred as "IJB95/96". Median follow-up was of 7 years.
The clinical value of the PCR-GGI for tamoxifen benefit was evaluated in three independent BC populations derived from two independent institutions:
1) A series of 77 frozen ER+, node-negative and -positive adjuvant tamoxifen-only treated primary BC collected at the John Radcliff Hospital (Oxford, UK) for which microarray data were available, referred as "OXFD". Median follow-up was of 5 years.
2) A series of 139 frozen ER+, node positive, adjuvant tamoxifen only treated primary BC collected from 1979 to 2001 at the Josephine Nefkens Institute of the Erasmus MC Rotterdam (NL), referred as "JNIadj". Median follow-up was of 3.5 years.
3) A series of 270 frozen ER+ BC samples from hormone-naïve patients that received first-line tamoxifen for metastatic disease [of which 52 patients received adjuvant chemotherapy], collected from 1979 to 2001 at the Josephine Nefkens Institute of the Erasmus MC Rotterdam, (NL), referred as "JNImeta". Median follow-up was of 5 years.
All frozen as well as FFPE samples contained >30% of invasive tumor cells. This study was approved by the ethics committee of the Institut Jules Bordet, the John Radcliffe Hospital and the Josephine Nefkens Institute of the Erasmus MC Rotterdam (MEC 02.953). The JNIadj and JNImeta samples involved coded tumor tissues and execution of this part of the study was performed in accordance with the Code of Conduct of the Federation of Medical Scientific Societies in the Netherlands http://www.fmwv.nl/.
RNA extraction and purification
RNA from frozen samples was extracted for the OXFD and IJBtest cohorts from four 10-μm sections using Trizol reagent according to the supplier's instructions (Invitrogen Life Technologies, Carlsbad, CA). For the JNIadj and JNImeta cohorts RNA was extracted from 30-μm sections using RNABee (Campro Scientific) as previously described .
RNA from FFPE samples was extracted from three 10-μm sections using the MasterPure Purification kit (Epicentre, Madison, WI) after paraffin removal with xylene. A DNase I treatment step was included.
To evaluate RNA quality and concentration, RNA from frozen samples were loaded onto an Agilent RNA 6000 NanoLabChip (RNA LabChip, Agilent Technology, Santa Clara, CA) in combination with the Bioanalyser 2100. Samples with a total area under the 28S and 18S bands of less than 15% of the total RNA band area, as well as a 28S/18S ratio of less than 1.1, were considered to be degraded. The RNA concentration from FFPE samples was analyzed with the NanoDrop® ND-1000 UV-Vis Spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA).
For the "OXFD" and "IJB95/96" cohorts reverse transcription (RT) was performed using a Super-Script First-Strand Synthesis kit for RT-PCR (Invitrogen Corp., Carlsbad, CA). Total RNA (300 ng) was reverse transcribed in a final volume of 21 μl with 50 ng of random hexamers. For the "JNIadj" and "JNImeta" cohorts cDNA was synthesized from 2 μg total RNA using random hexamers and oligo dTs, as previously described . A RNase H treatment step was included in all protocols.
Quantitative PCR (qPCR) reactions were performed in 96-well plates using Applied Biosystems Prism 7900 HT (TaqMan instruments) and a Stratagene Mx3000P QPCR System (Agilent Technologies, Waldbronn, Germany). Gene expressions were measured in duplicate using 5 ng equivalent cDNA per reaction well. Amplifications were performed in 25 μl PCR mixture containing 300-600 nM of each primer and 12,5 μl 2× SYBR Green PCR Master Mix (Applied Biosystems). After 2 min at 50°C and 10 min at 95°C, cDNA was subjected to 40 cycles of PCR with a denaturation step at 95°C for 30 sec followed by a combined annealing/extension step at 60°C for 1 min.
Primers of the selected genes were designed using the Primer Express software (PE Applied Biosystems-see Table 4 for primer sequences).
Four housekeeping genes were selected on the basis of the literature as reference genes for data normalization: TFRC, GUS, RPLPO and TBP . All PCR assays were run in duplicate. For all our samples the average value for the housekeeping gene is highly similar therefore all samples with no expression or a lower Ct value of the housekeeping genes compared to the value of the global population were excluded from the analyses. Six percent of frozen and 14% of FFPE samples were then excluded from the analyses.
PCR-GGI assay was developed at the Jules Bordet Institute on Jules Bordet populations and validated both at the Jules Bordet (IJB95/96, OXFD) and the Josephine Nefkens Institutes (JNIadj, JNImeta) on independent populations.
qRT-PCR data normalization and PCR-GGI computation
The average of the housekeeper genes values was used as reference and a Ct value of GGI genes was defined for each sample by taking off this average. To calculate the reduced gene expression grade index, we averaged the normalized values of the four GGI genes. This is illustrated by the following formula:
where s j is the PCR-GGI score for patient j, G i is the raw gene expression of the ith gene in the PCR-GGI signature for patient j, and is the mean of the housekeeping genes used for normalization.
In order to classify patients into low- and high-risk groups, a cut-off was defined as the middle point between the averages of the values of histological grade 1 and histological grade 3 tumors. This cut-off was defined on frozen samples (OXFD). To correct for the difference in scale of PCR-GGI scores between frozen and FFPE samples, we transformed the qRT-PCR scores derived from FFPE ER-positive patients such that the mean and standard deviation were equal to those from frozen samples:
where s j is the PCR-GGI score for patient j from a FFPE sample, μ f and μ p are the mean value of the PCR-GGI scores for ER-positive frozen and FFPE tumor samples respectively, and σ f and σ p are the standard deviation of the PCR-GGI scores of ER-positive frozen and FFPE tumor samples respectively.
Evaluation of estrogen receptor (ER) and Ki-67 expression by immunohistochemistry (IHC)
Paraffin-embedded blocks routinely fixed in neutral buffered formalin were cut on poly-L-Lysine-coated slides and stained with antibodies to ER clone 6F11 (dilution 1/40, Novocastra, Newcastle, UK) and to Ki-67 clone MIB-1 (dilution 1/50; DAKO, Carpinteria, CA). Antigen retrieval was performed in citrate buffer pH6 as previously described . The Ventana Nexes automated immunostainer (Ventana Medical Systems, Tucson, Az) was used with standard reagents. Tumors were defined as ER negative if < 10% of tumor cells had positive immunostaining and as highly proliferating samples if >15% of tumors cells had positive Ki-67 immunostaining.
ER evaluation by qRT-PCR
In order to consistently identify the ER status of BC patients, we clustered the tumors using the ESR1 qRT-PCR expressions by fitting Gaussian mixture models  with equal and diagonal variance for all clusters. We used the Bayesian Information Criterion  to test the number of components. Each tumor was then automatically classified as ER-positive or ER-negative using the maximum posterior probability of membership to these two clusters.
We considered DMFS of BC patients as endpoint for the "IJB95/96" and "JNIadj" data sets and PFS as endpoint for the "JNImeta" data set. Criteria for follow up, type of response and response to therapy were defined by standard International Union Against Cancer criteria of tumor response . PFS results were described previously . Survival curves were computed using Kaplan-Meier product limit estimator. Hazard ratios for continuous and discrete variables were estimated through Cox's proportional hazard regression models. A backward stepwise selection based on likelihood ratio tests was used for the identification of the covariates to be kept in the final Cox regression models.
All p-values were two-tailed and p-values < 0.05 were considered statistically significant. All statistical analyses were carried out using R version 2.5.1 and SPSS version 15 (SPSS Inc. 1999, Chicago IL.).
- BC :
- cDNA :
complementary deoxyribonucleic acid
- Ct :
- DMFS :
distant metastasis free survival
- ER :
- FFPE :
formalin fixed paraffin embedded
- GGI :
genomic grade index
- HR :
- IHC :
- IJB95/96 :
series of 212 FFPE primary BC sample collected at Jules Bordet Institute from 1995 to 1996
- IJBtest :
series of 19 primary BC from which frozen, FFPE tissues and corresponding microarray data were available
- JNIadj :
series of 139 frozen ER+, node positive, adjuvant tamoxifen treated primary BC collected from 1979 to 2001 at the Josephine Nefkens Institute of the Erasmus MC Rotterdam
- JNImeta :
series of 270 frozen ER+ BC from hormone-naïve patients that received first-line tamoxifen for metastatic disease, collected from 1979 to 2001 at the Josephine Nefkens Institute of the Erasmus MC Rotterdam
- MPF :
M-phase promoting factor
- OR :
- OXFD :
series of 77 frozen ER+ adjuvant tamoxifen treated primary BC collected at the John Radcliff Hospital for which microarray data were available
- PCR-GGI :
genomic grade index converts to a qRT-PCR assay
- PFS :
progression free survival
- qRT-PCR :
real-time reverse transcription polymerase chain reaction
- RNA :
Sotiriou C, Piccart MJ: Taking gene-expression profiling to the clinic: when will molecular signatures become relevant to patient care?. Nat Rev Cancer. 2007, 7 (7): 545-553. 10.1038/nrc2173. Review
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, Vijver Van de 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.
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 (10): 1239-1246. 10.1200/JCO.2006.07.1522. Erratum in: J Clin Oncol 2007, 25(24):3790
Wirapati P, Sotiriou C, Kunkel S, Farmer P, Pradervand S, Haibe-Kains B, Desmedt C, Ignatiadis M, Sengstag T, Schütz F, Goldstein DR, Piccart M, Delorenzi M: Meta-analysis of gene expression profiles in breast cancer: toward a unified understanding of breast cancer subtyping and prognosis signatures. Breast Cancer Res. 2008, 10 (4): R65-10.1186/bcr2124.
Desmedt C, Haibe-Kains B, Wirapati P, Buyse M, Larsimont D, Bontempi G, Delorenzi M, Piccart M, Sotiriou C: Biological processes associated with breast cancer clinical outcome depend on the molecular subtypes. Clin Cancer Res. 2008, 14 (16): 5158-5165. 10.1158/1078-0432.CCR-07-4756.
Vijver van de MJ, He YD, van't Veer L, Dai H, Hart AM, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ, Parrish M, Atsma D, Witteveen A, Glas A, Delahaye L, Velde van der T, Bartelink H, Rodenhuis S, Rutgers ET, H Friend S, Bernards R: A gene expression signature as a predictor of survival in breast cancer. The new England Journal of Medecine. 2002, 347 (25): 1999-2009. 10.1056/NEJMoa021967.
Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, van Gelder MEM, Yu J, Jatkoe T, Berns EM, Atkins D, Foekens J: A Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365: 671-679.
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 (11): 1723-1739. 10.1093/annonc/mdi352.
Paik S, Tang G, Shak S, Kim C, Baker J, Kim W, Cronin M, Baehner FL, Watson D, Bryant J, Costantino JP, Geyer CE, Wickerham DL, Wolmark N: Gene expression and benefit of chemotherapy in women with node-negative, estrogen receptor-positive breast cancer. J ClinOncol. 2006, 24 (23): 3726-3734. 10.1200/JCO.2005.04.7985.
Naderi A, Teschendorff AE, Barbosa-Morais NL, Pinder SE, Green AR, Powe DG, Roberston JF, Aparicio S, Ellis IO, Brenton JD, Caldas C: A Gene-Expression Signature as a Predictor of Survival in Breast Cancer across independent data sets. Oncogene. 2007, 26 (10): 1507-1516. 10.1038/sj.onc.1209920.
Miller LD, Smeds J, George J, Vega VB, Vergara L, Ploner A, Pawitan Y, Hall P, Klaar S, Liu ET, Bergh J: An expression signature for p53 status in human breast cancer predicts mutation status, transcriptional effects, and patient survival. Proc NatlAcadSci USA. 2005, 102: 13550-13555. 10.1073/pnas.0506230102.
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 (2): 303-9. 10.1007/s10549-008-0183-2.
Sjostrom J, Krajewski S, Franssila K, Niskanen E, Wasenius VM, Nordling S, Reed JC, Blomqvist C: A multivariate analysis of tumor biological factors predicting response to cytotoxic treatment in advanced breast cancer. Br J Cancer. 1998, 78: 812-815.
C Liedtke, C Hatzis, Symmans WF, C Desmedt, B Haibe-Kains, V Valero, H Kuerer, Hortobagyi GN, M Piccart-Gebhart, C Sotiriou, L Pusztai: 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.
Sieuwerts AM, Meijer-van Gelder ME, Timmermans M, Trapman AM, Garcia RR, Arnold M, Goedheer AJ, Portengen H, Klijn JG, Foekens JA: How ADAM-9 and ADAM-11 differentially from estrogen receptor predict response to tamoxifen treatment in patients with recurrent breast cancer: a retrospective study. Clin Cancer Res. 2005, 11 (20): 7311-7321. 10.1158/1078-0432.CCR-05-0560.
de Cremoux P, Bieche I, Tran-Perennou C, Vignaud S, Boudou E, Asselain B, Lidereau R, Magdelénat H, Becette V, Sigal-Zafrani B, Spyratos F: Inter-laboratory quality control for hormone-dependent gene expression in human breast tumors using real-time reverse transcription-polymerase chain reaction. Endocr Relat Cancer. 2004, 11 (3): 489-495. 10.1677/erc.1.00808.
Leake R, Barnes D, Pinder S, Ellis I, Anderson L, Anderson T, Adamson R, Rhodes T, Miller K, Walker R: Immunohistochemical detection of steroid receptors in breast cancer: a working protocol. J Clin Pathol. 2000, 53: 634-635. 10.1136/jcp.53.8.634.
McLachlan G, Peel D: Finite Mixture Models. 2000, John Wiley & Sons. New York, 419-
Schwarz G: Estimating the dimension of a model. Annals of Statistics. 1978, 6: 461-464. 10.1214/aos/1176344136.
Hayward JL, Carbone PP, Heuson JC, Kumaoka S, Segaloff A, Rubens RD: Assessment of response to therapy in advanced breast cancer: A project of the Programme on Clinical Oncology of the International Union Against Cancer, Geneva, Switzerland. Cancer. 1977, 39: 1289-1294. 10.1002/1097-0142(197703)39:3<1289::AID-CNCR2820390340>3.0.CO;2-F.
Foekens JA, Peters HA, Grebenchtchikov N, Look MP, Meijer-van Gelder ME, Geurts-Moespot A, Kwast van der TH, Sweep CG, Klijn JG: High tumour levels of vascular endothelial growth factor predict poor response to systemic therapy in advanced breast cancer. Cancer Res. 2001, 61: 5407-5414.
The authors thank the Breast Cancer Research Foundation (BCRF), the Swiss MEDIC Foundation, the Fonds National de la Recherche Scientifique (FNRS), the Fonds Ithier and the Netherlands Genomics Initiative (NGI)/Netherlands Organization for Scientific Research (NWO) for their research support. J Toussaint was supported by a grant from the FNRS (FRIA).
CS is named inventor on a patent application for the Genomic Grade signature used in this study. There are no other conflicts of interest.
JT has made contributions to method conception and design, acquisition of expression data, analysis and interpretation of data. He has been involved in drafting the manuscript and given final approval of the version to be published. AMS has made contributions to acquisition of expression data, analysis, interpretation of data and in drafting the manuscript. BHK and CD have contributed to statistical analyses, contributed to method design and in drafting the manuscript.GR and DL have been involved in the acquisition of data. VD has been involved in method design, acquisition of data, statistical analyses and in drafting the manuscript. ALH, MP, JAF and CS: have been involved in method design and in drafting the manuscript. All authors read and approved the final manuscript.
Jérôme Toussaint, Anieta M Sieuwerts contributed equally to this work.