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.



This article has been 
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.