Prognostic Utility of the 21-Gene Assay in Hormone Receptor–Positive Operable Breast Cancer Compared With Classical Clinicopathologic Features (original) (raw)

Abstract

Purpose

Adjuvant! is a standardized validated decision aid that projects outcomes in operable breast cancer based on classical clinicopathologic features and therapy. Genomic classifiers offer the potential to more accurately identify individuals who benefit from chemotherapy than clinicopathologic features.

Patients and Methods

A sample of 465 patients with hormone receptor (HR) –positive breast cancer with zero to three positive axillary nodes who did (n = 99) or did not have recurrence after chemohormonal therapy had tumor tissue evaluated using a 21-gene assay. Histologic grade and HR expression were evaluated locally and in a central laboratory.

Results

Recurrence Score (RS) was a highly significant predictor of recurrence, including node-negative and node-positive disease (P < .001 for both) and when adjusted for other clinical variables. RS also predicted recurrence more accurately than clinical variables when integrated by an algorithm modeled after Adjuvant! that was adjusted to 5-year outcomes. The 5-year recurrence rate was only 5% or less for the estimated 46% of patients who have a low RS (< 18).

Conclusion

The 21-gene assay was a more accurate predictor of relapse than standard clinical features for individual patients with HR-positive operable breast cancer treated with chemohormonal therapy and provides information that is complementary to features typically used in anatomic staging, such as tumor size and lymph node involvement. The 21-gene assay may be used to select low-risk patients for abbreviated chemotherapy regimens similar to those used in our study or high-risk patients for more aggressive regimens or clinical trials evaluating novel treatments.

INTRODUCTION

Adjuvant cytotoxic chemotherapy reduces the relative risk of recurrence by approximately 30% in women with operable breast cancer.1 Other adjuvant systemic therapies that reduce recurrence include endocrine therapy for hormone receptor (HR) –positive disease and trastuzumab for disease that overexpresses the human epidermal growth factor receptor 2 (HER-2) protein. With regard to chemotherapy, anthracyclines are now commonly used because they reduce relative recurrence risk by an additional 10% when combined with alkylating agents, such as cyclophosphamide.1 More recently, numerous trials have demonstrated a modest advantage when adding the taxanes paclitaxel or docetaxel to doxorubicin-based chemotherapy (whether used concurrently or sequentially after doxorubicin) and have been referred to as third-generation regimens.2-6 A recent report indicated that the benefits of taxane-containing adjuvant regimens were driven largely by a significant treatment effect in HR-negative disease, with only a modest treatment effect observed in HR-positive disease.5 The lower relapse rates in HR-positive disease observed in contemporary clinical trials have created challenges in evaluating new therapeutic strategies. Extremely large trials that must include thousands of patients are now required to detect improved outcomes for new treatment strategies. Use of standardized validated tools for patient selection could enhance the likelihood of success of such trials and allow them to be performed more efficiently.

Adjuvant! is a standardized validated instrument that projects patient outcomes based on classical clinicopathologic features and therapy, was designed to predict outcome at 10 years, and was validated for this end point.6 Although long-term outcome is generally accepted as fully reflective of the clinical utility of an adjuvant therapy, precise estimation of shorter-term outcomes at 5 years may offer advantages if found to be reflective of longer-term outcomes.

Several studies have demonstrated that multigene molecular markers predict relapse more accurately than clinical criteria in patients with negative or positive axillary nodes and in HR-positive or HR-negative disease.7-11 Of these signatures, only the 21-gene assay (Oncotype DX Recurrence Score [RS]; Genomic Health, Inc, Redwood City, CA) was developed and validated for patients with HR-positive, node-negative disease treated with endocrine therapy and is performed on routinely processed paraffin-embedded tumor tissue. The purpose of this study was to evaluate the prognostic utility of this 21-gene assay in patients with either node-negative or node-positive HR-positive breast cancer treated with doxorubicin-containing chemotherapy and to determine whether it could more reliably predict outcome at 5 years than standard clinicopathologic features. We also compared the prognostic utility of RS compared with clinicopathologic features and treatment information when integrated by an algorithm similar to but not identical to Adjuvant!, with the algorithm adjusted to 5-year outcomes. The purpose of this comparison was not to determine which method for estimating prognosis was superior, but rather to provide a more stringent test of the prognostic utility of the 21-gene assay. The development of such markers offer the potential to more accurately predict patient benefit and outcomes with specific therapies, to make more informed treatment recommendations, and to enrich for high-risk populations to select specific chemotherapy regimens or clinical trials evaluating more novel treatment strategies.

PATIENTS AND METHODS

Study Population and Treatment

The study used tumor specimens and clinical information from patients enrolled on trial E2197 (ClinicalTrials.gov identifier: NCT00003519), coordinated by the Eastern Cooperative Oncology Group (ECOG). The trial included 2,885 eligible and assessable patients with operable breast carcinoma and one to three positive axillary lymph nodes or negative axillary lymph nodes with the primary tumor measuring at least 1.1 cm in size. Patients were randomly assigned to receive four 3-week cycles of doxorubicin 60 mg/m2 plus cyclophosphamide 600 mg/m2 or docetaxel 60 mg/m2, plus endocrine therapy (if HR-positive) for 5 years or longer. Tamoxifen (20 mg daily for 5 years) was recommended beginning after completion of chemotherapy when the trial was initiated, although approximately 40% of patients eventually took an aromatase inhibitor at some point before or after 5 years when it was shown that these agents were more effective than tamoxifen.12 The treatment arms were well balanced with regard to median age (51 years), proportion with lymph node-negative disease (65%), and HR-positive disease (64%). After a median follow-up of 76 months, there was no significant difference between arms in disease-free survival (the primary study end point), relapse-free interval, or overall survival in the entire study group and in the population included in this analysis.

Specimen Analysis

Tumor specimens were evaluated for histologic grade using the modified Bloom-Richardson score by a single pathologist (F.L.B.) using 3- or 4-μm tissue sections stained by hematoxylin and eosin13 and for estrogen receptor (ER), progesterone receptor (PR), and HER-2 expression by immunohistochemistry by two pathologists simultaneously (F.L.B., S.B.) using two 1.0-mm tissue microarray cores using previously described methods summarized in the Appendix (online only).14,15 There was good concordance between the local institutional pathology laboratory and central laboratory for HR by immunohistochemistry (90%). All specimens were also analyzed for the Oncotype DX Recurrence Score as previously described.16

Case and Control Selection, End Points, and Statistical Analyses

The analysis included 465 patients with HR-positive disease, of whom 99 patients had experienced relapse (cases) and 366 patients did not (controls). The methods for case and control selection and statistical analyses are described in the Appendix. This includes a description of how the analysis was weighted to compensate for the unbalanced manner in which the cases were sampled relative to the controls. The primary end point was recurrence-free interval (RFI), defined as the time from trial entry to the first evidence of breast cancer recurrence (which included invasive breast cancer in local, regional, or distant sites, including the ipsilateral breast, but excluded new primary breast cancers in the opposite breast).17

RESULTS

Characteristics of Study Population and RS Distribution

There were no significant differences in the characteristics of sample included in this analysis compared with the not-in-sample group with HR-positive disease (Table 1), with the exception of estimated menopause distributions (P = .03) and the proportion with two to three nodes positive (P = .05). The raw numbers of patients in the sampled group and weighted estimates of the percents in the E2197 population with an RS that was low (< 18), intermediate (18 to 30), or high (≥ 31) were 198 patients (46%), 142 patients (30%), and 125 patients (24%), respectively. With regard to HER-2 expression based on central testing, 76 (16% weighted) of 465 patients were positive, including 21 (20% weighted) of the 99 patients who experienced relapse. There was concordance between central and local tumor grade in 249 (57%) of 437 cases in which information regarding local grade was known, a level of concordance that is consistent with that of other reports.16

Table 1.

Distribution of Patient Characteristics in Sample Compared With Not-in-Sample Cohort

Correlation of Recurrence With RS as Categoric or Continuous Variable

Continuous RS, without considering other variables, was a highly significant predictor of recurrence overall and for both the node-negative and node-positive patients (P < .001, linear trend test, in each case). When evaluated as a categoric variable, RS was also predictive of an elevated risk of recurrence regardless of nodal status (Fig 1A). The risk of recurrence was elevated in patients with two to three compared with zero or one positive nodes but was not elevated in patients with one positive versus zero positive nodes. For patients with an RS less than 18, approximately 3.3% (95% CI, 2.2% to 5.0%) of patients experienced recurrence within 5 years if there were zero to one positive nodes, and 7.9% (95% CI, 4.3% to 14.1%) experienced recurrence if there were two to three positive nodes. When modeled as a smooth continuous function using splines, there was a direct correlation between RS and recurrence risk up to an RS of approximately 40 (P < .001; Fig 1B), even if the analysis was restricted to only patients with HER-2–negative disease (P < .001; Fig 1C). Risk of recurrence did not increase with an RS greater than 40 in this group of patients treated with chemotherapy, consistent with previous reports suggesting greater benefit from chemotherapy if the RS was high.18

Fig 1.

(A) Five-year recurrence rates by Recurrence Score (RS) as a categoric variable by nodal status. RS risk categories: low risk is RS less than 18; intermediate risk is RS 18 to 30; high risk is RS ≥ 31. Relationship between RS and recurrence as a continuous variable for (B) all patients with hormone receptor (HR) –positive disease and (C) HR-positive and human epidermal growth factor receptor 2 (HER-2) –negative disease where HER-2–positive disease was excluded. The rug plots above the x-axis show the individual patients with recurrence in black and the individual patients without recurrence in green. The large CIs at higher RS reflect the relatively small number of patients with RS greater than 50. Pts, patients.

RS and Recurrence Risk When Adjusted for Clinical Variables

Proportional hazard models for recurrence were performed that included centrally determined tumor grade and HER-2 expression, age, tumor size, and number of positive axillary lymph nodes. RS was evaluated as a continuous linear variable in a manner that was similar to the B14 validation study; the hazard ratio (HR) for recurrence was calculated relative to an increment of 50 RS units, which is half the range of RS values and thus improves comparability of the HR with the HRs based on the clinical variables. When not including RS, factors associated with a significantly increased recurrence risk included two to three positive axillary nodes (but not one compared with zero positive nodes), young age, and tumor grade. When linear RS was included in the model, only two to three positive nodes and young age were associated with a higher risk, and there was also a strong trend for RS (HR for a 50-point difference in RS = 2.12; 95% CI, 0.97 to 4.65; P = .06, linear trend test; Table 2). If the model included locally determined tumor grade, which is more reflective of actual clinical practice, RS was highly associated with recurrence (HR = 3.13; 95% CI, 1.60 to 6.14; P = .0009). If the model included only 389 patients with HER-2/_neu_–negative tumors, RS was not predictive whether local or central grading was used (P = not significant, linear trend tests).

Table 2.

Results of Proportional Hazards Models for Recurrence

Comparison of Outcomes Predicted by the 21-Gene Assay and an Integrator of Clinicopathologic Information

This analysis was performed in four steps. We used an integrator of clinicopathologic information that was modeled after Adjuvant!, but adjusted to 5-year outcomes rather than 10-year outcomes. We will therefore refer to this tool as an integrator rather than Adjuvant!, because use of the latter term implies the use of the Adjuvant! model that is predictive of 10-year outcomes. First, outcomes projected by the integrator were determined for each patient using both local grade and central grade (Fig 2), which demonstrates some discordance in prediction based on the grade used. However, the 5-year recurrence rate of 10% (using central grade) predicted by the integrator was identical to the 10% (SE of mean ± 0.7%) observed recurrence rate, indicating that the integrator performed well in predicting the average outcome for the entire population.

Fig 2.

Five-year recurrence probabilities predicted by the integrator using local and central grade. For cases with the same result for local and central grade, the two values are equal (so the points are on the diagonal). Points above the diagonal have a higher risk grade category for central than for local and patients below the diagonal have a higher risk grade category for local than for central.

Second, we evaluated the concordance between prediction made by RS and the integrator using either risk group or risk percentile classification for comparison. For risk group comparisons, patients were classified as being in low, intermediate, or high integrator risk groups, proportionate to the raw distribution of RS risk groups of low (43%), intermediate (31%), and high (27%). For risk percentile classification, both RS and integrator recurrence probabilities were converted to ranks (that is, the values were used to rank the recurrence risk of the patients), and the ranks were then scaled to the interval (0, 100). Put another way, this method is similar to transformation of the recurrence risks to recurrence percentiles. Although use of the ranks puts the integrator and RS on a similar footing, the estimated effects are only relevant to populations with exactly the same distribution of values as in this study. The overall concordance rates between the integrator and RS risk groups were 166 (36%) of 465 when local grade was used and 179 (38%) of 465 when central grade was used, indicating poor concordance between the RS and integrator risk groups. When evaluated by risk percentile, the correlation coefficient was 0.10 when using local grade (Fig 3A) and 0.19 when using central grade (Fig 3B), again demonstrating poor concordance.

Fig 3.

Recurrence score risk percentiles compared with the integrator risk percentiles calculated using (A) local grade (correlation coefficient = 0.10) or (B) central grade (correlation coefficient = 0.19).

Third, we compared the predictive accuracy of RS with the integrator using three different methods, including a multivariate proportional hazards model, receiver operating characteristic (ROC) curves, and an analysis that evaluated whether RS provided additional information beyond that afforded by the integrator as well as whether the integrator provided additional information that was complementary to RS.

With regard to the proportional hazards model, in a joint linear trend model that included only RS and the integrator risk percentile as variables (and no other factors), the estimated effects for a 50-point difference were HR = 2.64 (95% CI, 1.80 to 3.87; P < .001) for RS and HR = 1.34 (95% CI, 0.94 to 1.91; P = .11) for the integrator using local grade and HR = 2.51 (95% CI, 1.71 to 3.70; P < .001) for RS and HR = 1.51 (95% CI, 1.07 to 2.13; P = .02) for the integrator using central grade. An RS by integrator interaction term added to either model was not significant (P = .94 in the first case, P = .56 in the second), showing that the effect of RS is largely independent of the level of the integrator risk.

With regard to ROC curve analysis, ROCs were developed for RS and the integrator, which facilitates evaluation of whether a marker predicts recurrence by a particular time (Fig 4). The area under the ROC curve (AUC) is a measure of the overall discriminatory power of the marker, with AUC = 0.50 corresponding to no discrimination and AUC = 1.0 corresponding to perfect prediction. The ROC AUC was 0.69 for RS used alone, 0.56 for the integrator using local grade, and 0.61 for the integrator using central grade, indicating that RS performed most accurately in predicting recurrence.

Fig 4.

Area under the receiver operating characteristic curves (AUC) for (A) Recurrence Score (RS), (B) the integrator using local grade, and (C) the integrator using central grade.

Fourth, we evaluated whether RS provided additional information with regard to the relative risk of recurrence for each integrator risk group, using low RS as the referent and central grade when calculating integrator risk (which optimized the performance of the integrator; Fig 5A). For the 43% of those who were in the low integrator risk group, the risk of recurrence was increased 2.6-fold and 4.0-fold for intermediate and high RS, respectively. For the 30% of patients in the intermediate integrator risk group, the risk of recurrence was increased 9.4-fold and 5.8-fold for intermediate and high RS, respectively. For the 24% of patients in the high integrator risk group, only high RS was associated with a significantly increased risk of relapse (2.6-fold). Therefore, RS provided the greatest discriminatory value in patients estimated to be at a low and intermediate risk by the integrator. We also evaluated whether the integrator provided additional information for each RS risk group (Fig 5B). The relative risk of recurrence was increased 3.15-fold for the high-risk integrator group compared with the low risk group for those who had a low RS and 2.44-fold increased for the intermediate integrator risk group compared with the low risk group for those with an intermediate RS, but not for other groups. This demonstrates that the integrator also provides information that is complementary to RS when the RS is low. Finally, we evaluated the absolute risk of recurrence for those with low, intermediate, and high RS in each integrator risk group (Fig 5C). The risk of recurrence was approximately 3% at 5 years for those with a low RS and were in the low risk integrator groups, identifying a group of patients who might do quite well with a short course of chemotherapy as used in this trial. For each integrator risk group, the risk of recurrence was 10% or higher for those with an intermediate or high RS, identifying patients who may be suitable candidates for more aggressive chemotherapy regimens or clinical trials evaluating novel agents.

Fig 5.

(A) Relative risk of recurrence by Recurrence Score (RS) within subsets defined by the integrator predicted risk level. The area of the box is inversely proportional to the variance of the log hazard ratio, and the line indicates the CI. (B) Relative risk of recurrence by integrator risk group within subsets defined by RS risk groups. The area of the box is inversely proportional to the variance of the log hazard ratio, and the line indicates the CI. (C) Absolute 5-year risk of recurrence (and 95% CIs) by RS risk group within subsets defined by integrator-predicted risk level. Intermed, intermediate. Absolute 5-year risk of recurrence (and 95% CIs) by RS risk group within subsets defined by integrator predicted risk level.

DISCUSSION

In the current study of patients with breast cancer with HR-positive disease and zero to three positive nodes who received contemporary chemohormonal therapy, the 21-gene assay was found to predict recurrence more accurately than standard clinicopathologic features individually and when integrated by an algorithm modeled after Adjuvant! that was adjusted to 5-year outcomes. Previous studies have demonstrated that the 21-gene assay predicts distant recurrence more accurately than individual clinical features in patients with HR-positive, axillary lymph node–negative breast cancer treated with tamoxifen.16 In this study, the assay was shown to provide information on distant and local/regional recurrence risk for individual patients treated with standard chemohormonal therapy independent of and beyond that provided by either nodal status or the integrator, which uses a combination of clinical features. Therefore, the 21-gene assay seems to provide complementary information to classical clinicopathologic features in patients with up to three positive axillary nodes.

For the 46% of patients with a low RS, there was a less than 3% risk of relapse at 5 years if there were zero to one positive nodes and 8% if there were two to three positive nodes. Because there was no arm without chemotherapy treatment, it is not possible here to directly evaluate whether the excellent outcome in those with low RS was related to good prognosis, chemotherapy benefit, or both. The ongoing TAILORx clinical trial (Trial Assigning Individualized Options for treatment) will provide additional information regarding the utility of chemotherapy in patients who have node-negative tumors associated with a midrange RS (11 to 25).19 Although additional studies will be required to determine whether adjuvant chemotherapy may be withheld from patients with one to three positive axillary lymph nodes with low RS, our findings suggest that this may be possible with additional properly designed studies. A recent report suggests that postmenopausal women with positive axillary lymph nodes and a low RS may not benefit from adjuvant anthracycline-based chemotherapy when added to tamoxifen,20 providing additional evidence to support this.

The analyses demonstrated a significant correlation between RS and recurrence when analyzed as continuous variable up to an RS of 40. Over its entire range from 0 to 100, RS was of borderline significance when adjusted for centrally determined tumor grade in the proportional hazards model (Table 2), reflecting the greater treatment effect for chemotherapy previously reported for individuals who have a tumor with a high RS.18 Adjuvant trastuzumab has been shown to reduce the risk of recurrence in HER-2–positive breast cancer,21 and HER-2 expression contributes to the RS.16 When our analysis was restricted to HER-2–negative disease, a significant association between RS and recurrence and pronounced chemotherapy treatment effect for high RS was also observed, indicating that RS provides clinically meaningful information in patients with HER-2–negative disease.

The median follow-up of patients included in this analysis was only 6 years, whereas Adjuvant! was designed for predicting outcomes at 10 years, and both Adjuvant! and the 21-gene profile were validated in populations using data with follow-up exceeding 10 years. The purpose of comparing the 21-gene assay with the integrator modeled after Adjuvant! was to provide a more stringent comparison with classical clinicopathologic factors, as Adjuvant! integrates multiple clinicopathologic and treatment factors. Additional follow-up will be required to determine 10-year outcomes.

Breast cancer mortality rates have declined considerably over the past decade, which has been attributed to improvements in adjuvant therapy and the more widespread use of screening mammography.22,23 Moreover, indications for adjuvant chemotherapy have expanded to include even low-risk patients,24 and more effective endocrine therapies have resulted in lower relapse rates for patients with HR-positive disease.25 More recently, adjuvant trastuzumab has been shown to further reduce the risk of recurrence in patients with HER-2–positive breast cancer (irrespective of HR status).26-28 Taken together, these advances have led to extremely low rates of relapse for patients diagnosed and treated during the past decade5 and have created the need to identify diagnostic tests that may predict outcomes more accurately than clinical criteria. Such diagnostic tests may be used for two purposes: first, to identify individuals with HR-positive disease who may be adequately treated with endocrine therapy alone, or endocrine therapy plus a short course of chemotherapy; second, such tests may be used to select high-risk patients for participation in clinical trials, thereby enriching for individuals most likely to benefit from innovative but potentially more toxic therapies. The TAILORx trial represents the first step toward implementation of this strategy by including patients who meet established clinical criteria for recommending chemotherapy, assigning treatment to hormonal therapy alone if the RS is low (< 11) or chemohormonal therapy if the RS is elevated (> 25), with random assignment between standard chemohormonal therapy to hormonal therapy alone for those with a midrange RS (11 to 25).19 Greater use of selection and enrichment strategies that include molecular markers offers potential for providing more informed treatment recommendations, improving the efficiency of clinical trials, and conserving resources.29,30

AUTHORS’ DISCLOSURES OF POTENTIAL CONFLICTS OF INTEREST

Although all authors completed the disclosure declaration, the following author(s) indicated a financial or other interest that is relevant to the subject matter under consideration in this article. Certain relationships marked with a “U” are those for which no compensation was received; those relationships marked with a “C” were compensated. For a detailed description of the disclosure categories, or for more information about ASCO's conflict of interest policy, please refer to the Author Disclosure Declaration and the Disclosures of Potential Conflicts of Interest section in Information for Contributors.

Employment or Leadership Position: Barrett H. Childs, Sanofi-aventis (C); Carl Yoshizawa, Genomic Health Inc (C); Steve Rowley, Sanofi-aventis (C); Steven Shak, Genomic Health Inc (C); Frederick L. Baehner, Genomic Health Inc (C) Consultant or Advisory Role: George W. Sledge Jr, Genomic Health Inc (C), Sanofi-aventis (C); Edith A. Perez, Sanofi-aventis (C); Joseph A. Sparano, Sanofi-aventis (C) Stock Ownership: Carl Yoshizawa, Genomic Health Inc; Steven Shak, Genomic Health Inc; Frederick L. Baehner, Genomic Health Inc Honoraria: Joseph A. Sparano, Sanofi-aventis Research Funding: Joseph A. Sparano, Sanofi-aventis Expert Testimony: None Other Remuneration: None

AUTHOR CONTRIBUTIONS

Conception and design: Lori J. Goldstein, Robert Gray, Sunil Badve, Steve Shak, Peter M. Ravdin, Nancy E. Davidson, George W. Sledge Jr, Joseph A. Sparano

Financial support: Lori J. Goldstein, Barrett H. Childs, Joseph A. Sparano

Administrative support: Lori J. Goldstein, Joseph A. Sparano

Provision of study materials or patients: Lori J. Goldstein, George W. Sledge Jr, Edith A. Perez, Lawrence N. Shulman, Silvana Martino, Joseph A. Sparano

Collection and assembly of data: Lori J. Goldstein, Robert Gray, Joseph A. Sparano

Data analysis and interpretation: Lori J. Goldstein, Robert Gray, Sunil Badve, Carl Yoshizawa, Steve Rowley, Steven Shak, Frederick L. Baehner, Peter M. Ravdin, Joseph A. Sparano

Manuscript writing: Lori J. Goldstein, Robert Gray, Sunil Badve, Steven Shak, Peter M. Ravdin, Nancy E. Davidson, George W. Sledge Jr, Edith A. Perez, Joseph A. Sparano

Final approval of manuscript: Lori J. Goldstein, Robert Gray, Sunil Badve, Barrett H. Childs, Carl Yoshizawa, Steve Rowley, Steven Shak, Frederick L. Baehner, Peter M. Ravdin, Nancy E. Davidson, George W. Sledge Jr, Edith A. Perez, Lawrence N. Shulman, Silvana Martino, Joseph A. Sparano

Glossary Terms

Histologic grade:

Histologic grade provides prognostic information in many tumors, including ovarian cancer. It is based on a combination of cellular features (nuclear, cytologic, and architectural). The more nuclear atypia and mitotic figures, the higher the grade.

Immunohistochemistry:

The application of antigen-antibody interactions to histochemical techniques. Typically, a tissue section is mounted on a slide and is incubated with antibodies (polyclonal or monoclonal) specific to the antigen (primary reaction). The antigen-antibody signal is then amplified using a second antibody conjugated to a complex of peroxidase-antiperoxidase (PAP), avidin-biotin-peroxidase (ABC) or avidin-biotin alkaline phosphatase. In the presence of substrate and chromogen, the enzyme forms a colored deposit at the sites of antibody-antigen binding. Immunofluorescence is an alternate approach to visualize antigens. In this technique, the primary antigen-antibody signal is amplified using a second antibody conjugated to a fluorochrome. On UV light absorption, the fluorochrome emits its own light at a longer wavelength (fluorescence), thus allowing localization of antibody-antigen complexes.

ROC (receiver operating characteristic) curves:

ROC curves plot the true-positive rate (sensitivity) against the false-positive rate (1-specificity) for different cut-off levels of a test. The area under the curve is a measure of the accuracy of the test. An area of 1.0 represents a perfect test (all true positives), whereas an area of 0.5 represents a worthless test.

Tissue microarray:

Used to analyze the expression of genes of interest simultaneously in multiple tissue samples, tissue microarrays consist of hundreds of individual tissue samples placed on slides ranging from 2 to 3 mm in diameter. Using conventional histochemical and molecular detection techniques, tissue microarrays are powerful tools to evaluate the expression of genes of interest in tissue samples. In cancer research, tissue microarrays are used to analyze the frequency of a molecular alteration in different tumor type, to evaluate prognostic markers, and to test potential diagnostic markers.

Multivariate proportional hazards model:

Proportional hazards or Cox regression modeling is a general method in medical statistics to analyze the influence of several (patient-specific) covariates on time-to-event end points. No assumption is made concerning the form of the underlying time-to-event curve. The only assumption made is that the effect of the covariates on the hazard rate in the study population is multiplicative and does not change over time.

ER (estrogen receptor):

Belonging to the class of nuclear receptors, estrogen receptors are ligand-activated nuclear proteins present in many breast cancer cells that are important in the progression of hormone-dependent cancers. After binding, the receptor-ligand complex activates gene transcription. There are two types of estrogen receptors (á and â). ERá is one of the most important proteins controlling breast cancer function. ERâ is present in much lower levels in breast cancer and its function is uncertain. Estrogen-receptor status guides therapeutic decisions in breast cancer.

PgR (progesterone receptor):

Like estrogen receptors, progesterone receptors are also nuclear proteins that are activated by the hormone progesterone in breast cancer cells that are hormone-dependent.

HER-2/neu (human epithelial growth factor receptor-2):

Also called ErbB2, HER-2/neu belongs to the EGFR family and is overexpressed in several solid tumors. Like EGFR, it is a tyrosine kinase receptor whose activation leads to proliferative signals within the cells. On activation, the HER family of receptors are known to form homodimers and heterodimers, each with a distinct signaling activity. Because HER-2 is the preferred dimerization partner when heterodimers are formed, it is important for signaling through ligands specific for any members of the family. It is typically overexpressed in several epithelial tumors.

Table A1.

Reasons for Patients Inclusion and Exclusion From This Analysis

Acknowledgments

We thank Adekunle Raji and other staff members at the Eastern Cooperative Oncology Group Pathology Coordinating Office at the Robert H. Lurie Comprehensive Cancer Center, Chicago, IL.

Appendix

CASE-CONTROL SELECTION

Because of the low overall recurrence rate in E2197, a stratified sampling design (similar to the stratified case-cohort design of Borgan et al [Lifetime Data Anal 6:39-58, 2000] was used, with recurrences sampled more heavily than nonrecurrences. Patients were sampled separately within groups defined by recurrence status, hormone receptor (HR) status (as determined by local laboratories), axillary nodal status (positive v negative), and treatment arm (doxorubicin plus cyclophosphamide v doxorubicin plus docetaxel), giving eight sampling strata with separate sampling from recurrences and nonrecurrences in each stratum.

E2197 included 2,952 patients with operable breast cancer, of whom 1,579 patients were potentially eligible for inclusion in this analysis, with reasons for exclusion listed in Table 1 and described herein. Although this report includes only patients with HR-positive disease (defined as being estrogen receptor [ER] and/or progesterone receptor [PR] positive), the sampling strategy included HR-positive and HR-negative disease.

Of the 1,579 potential patients eligible for the sampling and analysis, 191 patients (12%) had a recurrence and 1,388 (88%) did not have a recurrence. All patients with recurrence were included in the sampling (defined as case sample), plus a randomly selected sample of 641 patients without recurrence (defined as the control sample, based on the planned ratio of 1:3.5 for the case-control sampling), yielding a total of 832 patients for the analysis. Although the sampling stratification was based on HR expression determined in local laboratories, the final classification of HR status in this analysis is based on central HR expression testing (as described in the Methods section of the text).

Of the 832 patients identified, samples from each patient were sent from the Eastern Cooperative Oncology Group (ECOG) Pathology Coordinating Office (PCO) to Genomic Health. The case and control sample was selected by the coordinating statistician (R.G.), and all specimens were processed by the ECOG PCO and Genomic Health without knowledge of the recurrence status. Of the 832 samples sent to Genomic Health, 56 patients (7%) were excluded for the reasons shown in Appendix Table A1, leaving 776 patients (93%) with genomic data. Of these 776 patients, there were 465 patients that had HR-positive disease (as determined centrally), of whom 99 patients (21%) had a recurrence and 366 patients (79%) did not have a recurrence.

CENTRAL PATHOLOGY REVIEW

A representative tumor block from the primary tumor specimen was submitted to the ECOG PCO at the Robert H. Lurie Comprehensive Cancer Center (Chicago, IL), and all specimens underwent routine quality control evaluation to assure that there was adequate tumor. Central immunohistochemistry (IHC) for ER (antibodies 1D5 and ER-2-123) and PR (PcR 1294) was performed using the DakoCytomation EnVision+ System (Dako, Carpinteria, CA) in a two-step technique and were done by designating a proportion score (PS; range = 0 to 5), intensity score (IS; range = 0 to 3), and Allred score (AS = PS + IS; range = 0 to 8) for each case; an AS of greater than 2 was defined as positive, as previously described. HER-2 expression was defined as positive (3+) if there was intense membrane staining in at least 30% of cells using the DAKO HercepTest.15 Seven cases in the total sample were not assessable on central IHC. ER and PR status for these cases was determined using local results and HER-2 status was determined using genomic results.

REGULATORY APPROVALS AND MANUSCRIPT DEVELOPMENT PROCESS

The E2197 protocol was approved by the institutional review boards of all participating institutions and was carried out in accordance with the Declaration of Helsinki, current United States Food and Drug Administration Good Clinical Practices, and local ethical and legal requirements. ECOG designed and coordinated the study and was responsible for all aspects of the data collection and analysis. Other members of the North American Breast Cancer Intergroup participated and contributed patients to the study, including the Southwest Oncology Group, Cancer and Leukemia Group B, and the North Central Cancer Treatment Group. Only patients who gave consent for future research of their tumor specimen were included in the analysis. The use of specimens for this project was approved by the North American Intergroup Correlative Science Committee (http://ctep.cancer.gov/resources/tbci/correlative_studies.htmlwww.tbci.gov) and by the Northwestern University institutional review board (which oversees the ECOG PCO, where the specimens were banked and evaluated).

The authors assume responsibility for the overall content and integrity of the manuscript and vouch for the accuracy and completeness of the reported data; their views do not necessarily represent the official views of the National Cancer Institute (NCI). The authors made the decision to publish the analysis. The data analysis was performed by R.G., and the manuscript was written by J.A.S. and L.J.G.

DESCRIPTION OF STATISTICAL ANALYSIS COMPARING RECURRENCE SCORE WITH CLINICAL VARIABLES

The analyses described herein refer to the 465 patients with HR-positive disease. The primary test for the effect of Oncotype DX Recurrence Score (RS; Genomic Health, Inc, Redwood City, CA) on recurrence risk was prespecified as the weighted partial likelihood Wald test (see Description of Weighted Analysis Methods). The hypothesis was to be tested in the model with only RS, with the log hazard modeled as a linear function of RS (the Wald test for the coefficient of RS as a linear covariate is denoted as the linear trend test in the text of the article). If this test was significant, then the tests adjusting for other clinical and pathologic factors were to be done. The factors prespecified for inclusion in the model were tumor size, number of positive nodes, histologic grade (as determined by the local institution pathologist), and HER-2 as determined centrally, with age and menopause to be included if significant. Planned secondary analyses included evaluating the prognostic utility of grouped RS values with splits based on tertiles and on the standard cutoffs of 18 and 31. Because the final availability of tissue and assessability rates for the assays were unknown at the time the project was planned, power calculations were based on a case-cohort sample with 250 recurrences and 500 nonrecurrence patients. Simulations indicated that the primary comparison should have more than 80% power for an effect size corresponding to a 12% increase in recurrence risk for each 10-point increase in RS. Simulations with only 92 recurrences in the HR-positive group, somewhat less than the final sample, showed that there was still at least 80% power for an effect size corresponding to a 14% increase in recurrence risk for each 10-point increase in RS.

DESCRIPTION OF WEIGHTED ANALYSIS METHODS

Unlike the standard case-cohort design, only a subset of the recurrences from the E2197 study are included in the sample here, and the analysis methods are based on the general sampling theory of Lin (Biometrika 87:37-47, 2000). To estimate the magnitude of effects in the full E2197 population, sampling weights for each of the 16 groups in the sample are defined by the number of patients in the E2197 study in that group divided by the number in the sample. In the weighted analyses, contributions to estimators and other quantities, such as partial likelihoods, are multiplied by these weights. If the patients included in the case-control sample are randomly selected from the recurrences and nonrecurrences within each stratum, then the weighted estimators give consistent estimates of the corresponding quantities from the full E2197 sample. Because availability and analyzability of tissue samples was a factor in the selection, the possibility of systematic bias in the selection cannot be completely ruled out, but the comparisons in Table 1 of the text suggest that for many purposes, weighted analysis of the case-control sample should be representative of the full E2197 study. The weighted partial likelihood computed in this fashion is used for estimating hazard ratios and testing effects. This essentially gives the weighted pseudolikelihood estimator β̂3 of Chen and Lo (Biometrika 86:755-764, 1999), generalized to allow subsampling of cases, which was also discussed in the context of stratified case-cohort sampling as Estimator II in Borgan et al (Lifetime Data Anal 6:39-58, 2000). The variance of the partial likelihood estimators is estimated using the general approach of Lin (Biometrika 87:37-47, 2000), which leads to a generalization of the variance estimator from Borgan et al (Lifetime Data Anal 6:39-58, 2000) to allow subsampling of cases. Natural splines with 3 df (as computed by the R function ‘ns’) were used to obtain a flexible smooth model for the effect of RS. Tests for an RS effect based on this model use the 3 df Wald test. Weighted Kaplan-Meier estimators are used to estimate event-free rates. SEs of event-free rates are obtained by using general finite population sampling theory to estimate the variance of the corresponding weighted Nelson-Aalen cumulative hazard estimator L(t), and using the delta method to obtain the large sample variance of exp{-L(t)}, which is asymptotically equivalent to the weighted Kaplan-Meier estimator. CIs on event (or event-free) rates were computed using the normal approximation to the log cumulative hazard estimates and transformed to the event scale. Estimated event-free rates from Cox models were calculated using the generalization of Breslow's underlying cumulative hazard estimator given by Lin (Biometrika 87:37-47, 2000), with the variance of this estimator determined using the general theory given in that article. Weighted averages, with proportions estimated using weighted averages of indicator variables, are also used for estimating the distribution of factors and for comparing the distributions between the overall E2197 study population and the genomic sample. Tests comparing factor distributions are based on asymptotic normality of the difference in weighted averages.

DESCRIPTION OF ANALYTIC METHODS FOR COMPARISON OF RS WITH ADJUVANT!

The predicted 5-year recurrence risk was computed using a batch processor provided by P.M.R., which he had adapted for estimation of 5-year recurrence risk. In calculating the recurrence risk using this processor, HR status (HR positive = ER or PR positive v HR negative = ER and PR negative) was based on central IHC if available and local determination if not (seven cases). Predicted recurrence risk was computed using both local grade and central grade. Because information on comorbidities was not available, all cases were coded as Minor Problems. Chemo was coded as CA*4, CMF, FE(50)C, and hormonal therapy was coded as tamoxifen. Although some patients have received aromatase inhibitor therapy, most have only begun aromatase inhibitors fairly recently, so the tamoxifen classification was thought to be reasonably reflective of the hormonal therapy benefit here, even for these cases. HER-2 status (based on central IHC) was included as an additional risk factor. Most patients who were classified locally as ER- and PR-negative and centrally as ER- or PR-positive were not given hormonal therapy. Because Adjuvant! adjusts recurrence risk for HR-positive patients based on hormonal therapy, it seemed appropriate to use the chemo only predicted recurrence risk for this group instead of the chemo plus hormones risk used for the other patients, unless the patients were reported to have received hormonal therapy (this increases the recurrence risk for 39 patients). In all cases, it seemed that this change gave slightly better prediction of recurrence.

Receiver operating characteristic curves for 5-year recurrence risk are computed using the naïve approach described in Section 2.2 of Heagerty et al (Biometrics 56:337-344, 2000), with weighted Kaplan-Meier estimators substituted for the ordinary estimators in the formulas there. Because these naive receiver operating characteristic estimators are generally not monotone, isotonic regression was applied to give a monotone curve.

published online ahead of print at www.jco.org August 4, 2008.

Supported in part by the United States Department of Health and Human Services and the National Institutes of Health (Grants No. CA23318 to the Eastern Cooperative Oncology Group [ECOG] statistical center, CA66636 to the ECOG data management center, CA21115 to the ECOG coordinating center, CA25224 to North Central Cancer Treatment Group, CA32291 to Cancer and Leukemia Group B, and CA32012 to Southwest Oncology Group), and by a grant from Sanofi-aventis.

Presented in part at the 43rd Annual Meeting of the American Society of Clinical Oncology, June 1-5, 2007, Chicago, IL; and at the San Antonio Breast Cancer Symposium, San Antonio, TX, December 13-16, 2007.

Authors’ disclosures of potential conflicts of interest and author contributions are found at the end of this article.

Clinical trial information can be found for the following: NCT00003519.

REFERENCES

• 1.Early Breast Cancer Trialists Collaborative Group: Effects of chemotherapy and hormonal therapy for early breast cancer on recurrence and 15-year survival: An overview of the randomised trials. Lancet 365:1687-1717, 2005 [DOI] [PubMed] [Google Scholar]
• 2.Henderson IC, Berry DA, Demetri GD, et al: Improved outcomes from adding sequential paclitaxel but not from escalating doxorubicin dose in an adjuvant chemotherapy regimen for patients with node-positive primary breast cancer. J Clin Oncol 21:976-983, 2003 [DOI] [PubMed] [Google Scholar]
• 3.Mamounas EP, Bryant J, Lembersky B, et al: Paclitaxel after doxorubicin plus cyclophosphamide as adjuvant chemotherapy for node-positive breast cancer: Results from NSABP B-28. J Clin Oncol 23:3686-3696, 2005 [DOI] [PubMed] [Google Scholar]
• 4.Martin M, Pienkowski T, Mackey J, et al: Adjuvant docetaxel for node-positive breast cancer. N Engl J Med 352:2302-2313, 2005 [DOI] [PubMed] [Google Scholar]
• 5.Berry DA, Cirrincione C, Henderson IC, et al: Estrogen-receptor status and outcomes of modern chemotherapy for patients with node-positive breast cancer. JAMA 295:1658-1667, 2006 [DOI] [PMC free article] [PubMed] [Google Scholar]
• 6.Olivotto IA, Bajdik CD, Ravdin PM, et al: Population-based validation of the prognostic model ADJUVANT! for early breast cancer. J Clin Oncol 23:2716-2725, 2005 [DOI] [PubMed] [Google Scholar]
• 7.van 't Veer LJ, Dai H, van de Vijver MJ, et al: Gene expression profiling predicts clinical outcome of breast cancer. Nature 415:530-536, 2002 [DOI] [PubMed] [Google Scholar]
• 8.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 347:1999-2009, 2002 [DOI] [PubMed] [Google Scholar]
• 9.Foekens JA, Atkins D, Zhang Y, et al: Multicenter validation of a gene expression-based prognostic signature in lymph node-negative primary breast cancer. J Clin Oncol 24:1665-1671, 2006 [DOI] [PubMed] [Google Scholar]
• 10.Chang HY, Nuyten DS, Sneddon JB, et al: Robustness, scalability, and integration of a wound-response gene expression signature in predicting breast cancer survival. Proc Natl Acad Sci U S A 102:3738-3743, 2005 [DOI] [PMC free article] [PubMed] [Google Scholar]
• 11.Ma XJ, Hilsenbeck SG, Wang W, et al: The HOXB13:IL17BR expression index is a prognostic factor in early-stage breast cancer. J Clin Oncol 24:4611-4619, 2006 [DOI] [PubMed] [Google Scholar]
• 12.Baum M, Buzdar A, Cuzick J, et al: Anastrozole alone or in combination with tamoxifen versus tamoxifen alone for adjuvant treatment of postmenopausal women with early-stage breast cancer: Results of the ATAC (Arimidex, Tamoxifen Alone or in Combination) trial efficacy and safety update analyses. Cancer 98:1802-1810, 2003 [DOI] [PubMed] [Google Scholar]
• 13.Dalton LW, Page DL, Dupont WD: Histologic grading of breast carcinoma: A reproducibility study. Cancer 73:2765-2770, 1994 [DOI] [PubMed] [Google Scholar]
• 14.Harvey JM, Clark GM, Osborne CK, et al: Estrogen receptor status by immunohistochemistry is superior to the ligand-binding assay for predicting response to adjuvant endocrine therapy in breast cancer. J Clin Oncol 17:1474-1481, 1999 [DOI] [PubMed] [Google Scholar]
• 15.Wolff AC, Hammond ME, Schwartz JN, et al: American Society of Clinical Oncology/College of American Pathologists Guideline Recommendations for Human Epidermal Growth Factor Receptor 2 Testing in Breast Cancer. Arch Pathol Lab Med 131:18, 2007 [DOI] [PubMed] [Google Scholar]
• 16.Paik S, Shak S, Tang G, et al: A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N Engl J Med 351:2817-2826, 2004 [DOI] [PubMed] [Google Scholar]
• 17.Hudis CA, Barlow WE, Costantino JP, et al: Proposal for standardized definitions for efficacy end points in adjuvant breast cancer trials: The STEEP system. J Clin Oncol 25:2127-2132, 2007 [DOI] [PubMed] [Google Scholar]
• 18.Paik S, Tang G, Shak S, et al: Gene expression and benefit of chemotherapy in women with node-negative, estrogen receptor-positive breast cancer. J Clin Oncol 24:3726-3734, 2006 [DOI] [PubMed] [Google Scholar]
• 19.Sparano JA: TAILORx: Trial assigning individualized options for treatment (Rx). Clin Breast Cancer 7:347-350, 2006 [DOI] [PubMed] [Google Scholar]
• 20.Albain K, Barlow W, Shak S, et al: Prognostic and predictive value of the 21-gene recurrence score assay in postmenopausal, node-positive, ER-positive breast cancer (S8814,INT0100). Breast Cancer Res Treat 106, 2007. (suppl 1; abstr 10)
• 21.Viani GA, Afonso SL, Stefano EJ, et al: Adjuvant trastuzumab in the treatment of HER-2-positive early breast cancer: A meta-analysis of published randomized trials. BMC Cancer 7:153, 2007 [DOI] [PMC free article] [PubMed] [Google Scholar]
• 22.Berry DA, Cronin KA, Plevritis SK, et al: Effect of screening and adjuvant therapy on mortality from breast cancer. N Engl J Med 353:1784-1792, 2005 [DOI] [PubMed] [Google Scholar]
• 23.Shen Y, Yang Y, Inoue LY, et al: Role of detection method in predicting breast cancer survival: Analysis of randomized screening trials. J Natl Cancer Inst 97:1195-1203, 2005 [DOI] [PubMed] [Google Scholar]
• 24.Carlson RW, Brown E, Burstein HJ, et al: NCCN Task Force Report: Adjuvant Therapy for Breast Cancer. J Natl Compr Canc Netw 4:S1-S26, 2006. (suppl) [PubMed] [Google Scholar]
• 25.Ryan PD, Goss PE: The evolving role of aromatase inhibitors in the treatment of early-stage breast cancer. Nat Clin Pract Oncol 2:596-597, 2005 [DOI] [PubMed] [Google Scholar]
• 26.Romond EH, Perez EA, Bryant J, et al: Trastuzumab plus adjuvant chemotherapy for operable HER2-positive breast cancer. N Engl J Med 353:1673-1684, 2005 [DOI] [PubMed] [Google Scholar]
• 27.Joensuu H, Kellokumpu-Lehtinen PL, Bono P, et al: Adjuvant docetaxel or vinorelbine with or without trastuzumab for breast cancer. N Engl J Med 354:809-820, 2006 [DOI] [PubMed] [Google Scholar]
• 28.Piccart-Gebhart MJ, Procter M, Leyland-Jones B, et al: Trastuzumab after adjuvant chemotherapy in HER2-positive breast cancer. N Engl J Med 353:1659-1672, 2005 [DOI] [PubMed] [Google Scholar]
• 29.Sargent DJ, Conley BA, Allegra C, et al: Clinical trial designs for predictive marker validation in cancer treatment trials. J Clin Oncol 23:2020-2027, 2005 [DOI] [PubMed] [Google Scholar]
• 30.Jiang W, Freidlin B, Simon R: Biomarker-adaptive threshold design: A procedure for evaluating treatment with possible biomarker-defined subset effect. J Natl Cancer Inst 99:1036-1043, 2007 [DOI] [PubMed] [Google Scholar]