Purpose:

To validate the clinical performance of the OncoMasTR Risk Score in the biomarker cohort of Austrian Breast and Colorectal Cancer Study Group (ABCSG) Trial 8.

Experimental Design:

We evaluated the OncoMasTR test in 1,200 formalin-fixed, paraffin-embedded (FFPE) surgical specimens from postmenopausal women with estrogen receptor (ER)–positive, human epidermal growth factor receptor 2 (HER2)–negative primary breast cancer with 0 to 3 involved lymph nodes in the prospective, randomized ABCSG Trial 8. Time to distant recurrence (DR) was analyzed by Cox models.

Results:

The OncoMasTR Risk Score categorized 850 of 1,087 (78.2%) evaluable patients as “low risk”. At 10 years, the DR rate for patients in the low-risk group was 5.8% versus 21.1% for patients in the high-risk group (P < 0.0001, absolute risk reduction 15.3%). The OncoMasTR Risk Score was highly prognostic for prediction of DR in years 0 to 10 in all patients [HR 1.91, 95% confidence interval (CI) 1.62–2.26, P < 0.0001; C-index 0.73], in patients that were node negative (HR 1.79, 95% CI, 1.43–2.24, P < 0.0001; C-index 0.72), and in patients with 1 to 3 involved lymph nodes (HR 1.93, 95% CI, 1.44–2.58, P < 0.0001; C-index 0.71). The OncoMasTR Risk Score provided significant additional prognostic information beyond clinical parameters, Ki67, Nottingham Prognostic Index, and Clinical Treatment Score.

Conclusions:

OncoMasTR Risk Score is highly prognostic for DR in postmenopausal women with ER-positive, HER2-negative primary breast cancer with 0 to 3 involved lymph nodes. In combination with prior validation studies, this fully independent validation in ABCSG Trial 8 provides level 1B evidence for the prognostic capability of the OncoMasTR Risk Score.

Translational Relevance

Accurate assessment of the individual risk of breast cancer recurrence after surgery is crucial in making optimal decisions about adjuvant treatment. We evaluated the clinical performance of the OncoMasTR Risk Score in postmenopausal women with estrogen receptor (ER)–positive, human epidermal growth factor receptor 2 (HER2)–negative primary breast cancer with 0 to 3 involved lymph nodes who were included in the prospective, randomized Austrian Breast and Colorectal Cancer Study Group (ABCSG) Trial 8. The OncoMasTR Risk Score is highly prognostic for prediction of distant recurrence in years 0 to 10 in all patients as well as in patients that were node negative and in patients with 1 to 3 involved lymph nodes. In combination with prior validation studies, this fully independent validation in ABCSG Trial 8 provides level 1B evidence for the prognostic capability of the OncoMasTR Risk Score. By gaining additional molecular prognostic information using the OncoMasTR test, patients with breast cancer can be identified with a distant recurrence risk that is low enough to avoid adjuvant chemotherapy or other treatment escalations.

Several prognostic gene signatures have been developed to estimate the risk of breast cancer recurrence after surgery and to guide adjuvant treatment decisions. The MammaPrint Score, the Oncotype DX Recurrence Score, the EndoPredict Score, the Prosigna PAM50 Risk of Recurrence Score, and the Breast Cancer Index have previously been validated as prognostic tests (1–7).

OncoMasTR is a novel prognostic test with the aim of accurately estimating the risk of recurrence for patients with estrogen receptor (ER)–positive, human epidermal growth factor receptor 2 (HER2)–negative, early-stage breast cancer. The signature was developed using a novel bioinformatic approach that identified a shared network of master transcriptional regulators (MTR) from which the combination of only 3 genes, namely Forkhead Box M1 (FOXM1), Pituitary Tumor Transforming Gene 1 (PTTG1), and Zinc Finger Protein 367 (ZNF367), formed the OncoMasTR Molecular Score (8, 9). The final OncoMasTR Risk Score combines the OncoMasTR Molecular Score with nodal status and tumor size. The OncoMasTR test has been analytically validated in accordance with industry standard guidelines from the Clinical Laboratory Standards Institute (CLSI; ref. 10).

The clinical performance of the OncoMasTR test was initially validated in 646 patients with ER-positive, HER2-negative, early-stage breast cancer with 0 to 3 involved lymph nodes enrolled into the ATAC trial (11). These results suggested that the OncoMasTR test was highly prognostic for the overall population and for node-negative disease as well as for both early (0–5 years) and late (5–10 years) distant recurrence (DR). In addition, the OncoMasTR Risk Score added significant prognostic value to Ki67, the Nottingham Prognostic Index (NPI), and the Clinical Treatment Score (CTS; ref. 11).

In the present study, we independently validated the clinical performance of the OncoMasTR Risk Score in a large cohort of patients with early-stage breast cancer who were enrolled into the Austrian Breast & Colorectal Cancer Study Group (ABCSG) Trial 8, a large, prospective, randomized clinical trial with long-term follow-up comparing the efficacy of two adjuvant endocrine treatment regimens (12, 13).

ABCSG trial 8

The objective of ABCSG Trial 8 was to compare the efficacy of the sequence of adjuvant tamoxifen (TAM) for 2 years followed by anastrozole (ANA) for 3 years with tamoxifen alone for 5 years. From January 1996 to June 2004, a total of 3,901 postmenopausal women with hormone receptor–positive, early-stage, grade 1 (G1) or grade 2 (G2) breast cancer who had undergone primary surgery with or without radiotherapy were randomized (12, 13). Patients in the trial did not receive adjuvant chemotherapy and/or anti-HER2 treatment (patients with high-grade tumors were excluded). The results from this study indicate that postmenopausal women already receiving adjuvant tamoxifen for the treatment of hormone receptor–positive early breast cancer benefit from switching to adjuvant anastrozole, rather than continuing with tamoxifen (12, 13). The study was conducted in accordance with the Declaration of Helsinki, approved by the responsible ethics committees, and all patients gave written informed consent.

Specimen collection and OncoMasTR test

All patients (n = 3,901) included in ABCSG Trial 8 were eligible for the translational study and participating centers were requested to provide a tumor block of their patients. Tumor specimens were obtained at the time of surgery prior to the adjuvant therapy. Details with regard to the collection of the samples and the preparation of sections were previously published as part of translational studies in this population (5, 7, 14–18). Formalin-fixed, paraffin-embedded (FFPE) surgical breast cancer specimens from 1,200 patients were retrospectively collected and used for the OncoMasTR test. A flowchart illustrating patient and sample selection is shown in Fig. 1. The collection of test data was performed in the ABCSG central research lab at the Medical University of Vienna, blinded to patient clinical and outcome parameters.

Figure 1.

CONSORT diagram for the process of patient and tumor block selection. 1Includes 1 HER2-positive/≥ 4 positive nodes sample. 2Includes 1 ER-negative/≥ 4 positive nodes sample, 1 HER2-positive/≥ 4 positive nodes sample, and 2 samples with insufficient quality/≥ 4 positive nodes.

Figure 1.

CONSORT diagram for the process of patient and tumor block selection. 1Includes 1 HER2-positive/≥ 4 positive nodes sample. 2Includes 1 ER-negative/≥ 4 positive nodes sample, 1 HER2-positive/≥ 4 positive nodes sample, and 2 samples with insufficient quality/≥ 4 positive nodes.

Close modal

One hematoxylin and eosin (H&E)–stained slide was prepared from each paraffin block and reviewed by an experienced breast pathologist (M. Rudas) to confirm the presence of invasive breast carcinoma. If more than 30% invasive tumor was present, macro-dissection was performed and tumor areas were scraped off the slide. In addition, in cases with extensive intraductal component, ductal carcinoma in situ (DCIS) was removed before performing the assay. Five adjacent sections of 5 μmol/L thickness from the FFPE block, or necessary tissue quantity (if the tumor-cell content was less than the prespecified 30%), was extracted using the Qiagen RNeasy FFPE Kit (Qiagen), followed by an optimized DNase digestion step to remove the genomic DNA background. The extracted RNA was read on NanoDrop for its concentration and purity.

The extracted and quality-controlled RNA was mixed with 1-Step TaqPath Master Mix (OncoMark) and RT-PCR grade water and applied to the OncoMasTR test plate, which is a 96-well plate prespotted with dried oligonucleotides for analyzing the expression levels of individual MTRs and reference genes. The RNA was analyzed in triplicates per gene target. There was an in-line positive control and negative control for each individual run. The OncoMasTR Test was performed on an Applied Biosystems 7500 Fast Dx Real-Time PCR instrument (Applied Biosystems). The OncoMasTR test reagents (including OncoMasTR plate, Reagent Kit, Instructions for Use), Qiagen RNeasy FFPE Kit, and the Applied Biosystems 7500 Fast DX Real-Time PCR instrument were provided by OncoMark for the duration of the analyses. The relative normalized expressions of the prognostic genes were translated into the OncoMasTR Risk Score using a clinically validated algorithm incorporating tumor size and lymph node status (8, 11). The OncoMasTR Molecular Score ranges from 0 to 100 and the OncoMasTR Risk Score ranges from 0 to 10, where a score of <5 corresponds to low risk and a score of ≥5 corresponds to high risk of DR within 10 years of initial diagnosis.

Statistical analysis

All data analyses were carried out according to a predefined statistical analysis plan. The prospectively defined primary endpoint of the study was time to DR defined as time from randomization until the occurrence of distant metastasis (ignoring loco-regional recurrences). All secondary carcinoma including contralateral breast cancer, and all deaths were treated as censoring events. In case of no event, the outcome was censored at the date of last contact. Follow-up was censored at 10 years.

Time to DR was analyzed by univariable and multivariable Cox proportional hazards regression models and described by HRs with 95% confidence intervals (CI), likelihood ratio tests, and by Kaplan–Meier survival curves.

Our primary objectives were to evaluate whether the OncoMasTR Molecular Score as a continuous variable, the OncoMasTR Risk Score as a continuous variable, and the OncoMasTR Risk Category as a categorical variable are prognostic for time to DR in all patients (0–3 involved nodes). The prognostic value of the OncoMasTR Scores (molecular score, risk score, and risk category) was also assessed in patients who were node negative and patients with 1 to 3 involved nodes. Likelihood ratio (LR) tests from uni- and multivariable Cox models were used to evaluate whether additional prognostic information was provided when the OncoMasTR Scores was added to age (as continuous variable), tumor size (as continuous variable), nodal status (pN1 vs. pN0), grading (G2 vs. G1), ER expression (as continuous variable), progesterone receptor (PR) expression (as continuous variable), Ki67 expression (as continuous variable), treatment arm (TAM/ANA vs. TAM), the NPI, or the CTS (19–22). The proportional hazards assumption was checked with weighted Schoenfeld residuals. The continuous relationship between DR risk and numerical risk scores was estimated using Royston/Parmar models (23).

Categorical baseline data according to OncoMasTR Risk Category were compared in univariable analysis using the χ2 or Fisher exact test depending on the expected cell frequencies. Continuous baseline data were compared using Wilcoxon test. All reported P values are from two-sided tests. All results with P values <0.05 were considered statistically significant. Statistical analyses were performed by members of the biostatistics group at ABCSG using statistical analysis system (SAS) software (SAS version 9.4). This study meets the Reporting Recommendations for Tumor Marker Prognostic Studies (REMARK; ref. 24).

Patient characteristics

Patient and sample selection is shown in a consort diagram (Fig. 1). Main clinical and laboratory parameters of all patients included in the present study (also compared with trial patients without tumor blocks) are summarized in Supplementary Table S1.

Mean patient age was 63.3 years. All patients had ER-positive and HER2-negative tumors, 69.5% had tumors smaller than 2 cm, 73.2% had node-negative disease, and 77.6% had G2 grading. Breast conservation was achieved in 82.1%. Data from 1,087 of the 1,200 samples (90.6%) met the prespecified OncoMasTR data quality criteria and were included in the analysis. Within the 1,087 patients, 850 (78.2%) were categorized as low risk and 237 (21.8%) as high risk (Table 1). High-risk patients were significantly older, had larger tumors, more frequently node-positive disease, more tumors with G2 grading (as compared with low grade), more progesterone receptor–negative tumors, and a higher mastectomy rate (Table 1).

Table 1.

Baseline comparison between low-risk and high-risk patients.

Low riskHigh riskTotal
Variable(N = 850)(N = 237)(N = 1,087)P value
Age (years)    <0.001 
N 850 (100.0%) 237 (100.0%) 1,087 (100.0%)  
 Mean (SD) 62.8 (7.8) 65.2 (7.8) 63.3 (7.9)  
 Median (Q1–Q3) 62.0 (57.0–69.0) 65.0 (59.0–71.0) 63.0 (57.0–70.0)  
 Min–max 45.0–79.0 45.0–79.0 45.0–79.0  
T-stage (%)    <0.001 
 pT1a 1 (0.1%) 1 (0.1%)  
 pT1b 170 (20.0%) 4 (1.7%) 174 (16.0%)  
 pT1c 484 (56.9%) 97 (40.9%) 581 (53.4%)  
 pT2 189 (22.2%) 132 (55.7%) 321 (29.5%)  
 pT3 6 (0.7%) 4 (1.7%) 10 (0.9%)  
Tumor size (cm)    <0.001 
N 848 (99.8%) 235 (99.2%) 1,083 (99.6%)  
 Missing 2 (0.2%) 2 (0.8%) 4 (0.4%)  
 Mean (SD) 1.6 (0.8) 2.2 (0.9) 1.8 (0.9)  
 Median (Q1–Q3) 1.5 (1.1–2.0) 2.2 (1.6–2.6) 1.5 (1.2–2.2)  
 Min–max 0.4–10.0 0.7–7.0 0.4–10.0  
Nodal status (%)    <0.001 
 pN0 680 (80.0%) 116 (48.9%) 796 (73.2%)  
 pN1 170 (20.0%) 121 (51.1%) 291 (26.8%)  
Grading (%)    <0.001 
 G1 181 (21.3%) 18 (7.6%) 199 (18.3%)  
 G2 633 (74.5%) 210 (88.6%) 843 (77.6%)  
 GX 36 (4.2%) 9 (3.8%) 45 (4.1%)  
ER (%)     
 Positive 850 (100.0%) 237 (100.0%) 1,087 (100.0%)  
PR (%)    0.02 
 Negative 53 (6.2%) 25 (10.5%) 78 (7.2%)  
 Positive 797 (93.8%) 212 (89.5%) 1,009 (92.8%)  
Type of surgery (%)    0.003 
 Breast conserving 713 (83.9%) 179 (75.5%) 892 (82.1%)  
 Mastectomy 137 (16.1%) 58 (24.5%) 195 (17.9%)  
Treatment arm (%)    0.55 
 2y TAM + 3y ANA 424 (49.9%) 113 (47.7%) 537 (49.4%)  
 5y TAM 426 (50.1%) 124 (52.3%) 550 (50.6%)  
Low riskHigh riskTotal
Variable(N = 850)(N = 237)(N = 1,087)P value
Age (years)    <0.001 
N 850 (100.0%) 237 (100.0%) 1,087 (100.0%)  
 Mean (SD) 62.8 (7.8) 65.2 (7.8) 63.3 (7.9)  
 Median (Q1–Q3) 62.0 (57.0–69.0) 65.0 (59.0–71.0) 63.0 (57.0–70.0)  
 Min–max 45.0–79.0 45.0–79.0 45.0–79.0  
T-stage (%)    <0.001 
 pT1a 1 (0.1%) 1 (0.1%)  
 pT1b 170 (20.0%) 4 (1.7%) 174 (16.0%)  
 pT1c 484 (56.9%) 97 (40.9%) 581 (53.4%)  
 pT2 189 (22.2%) 132 (55.7%) 321 (29.5%)  
 pT3 6 (0.7%) 4 (1.7%) 10 (0.9%)  
Tumor size (cm)    <0.001 
N 848 (99.8%) 235 (99.2%) 1,083 (99.6%)  
 Missing 2 (0.2%) 2 (0.8%) 4 (0.4%)  
 Mean (SD) 1.6 (0.8) 2.2 (0.9) 1.8 (0.9)  
 Median (Q1–Q3) 1.5 (1.1–2.0) 2.2 (1.6–2.6) 1.5 (1.2–2.2)  
 Min–max 0.4–10.0 0.7–7.0 0.4–10.0  
Nodal status (%)    <0.001 
 pN0 680 (80.0%) 116 (48.9%) 796 (73.2%)  
 pN1 170 (20.0%) 121 (51.1%) 291 (26.8%)  
Grading (%)    <0.001 
 G1 181 (21.3%) 18 (7.6%) 199 (18.3%)  
 G2 633 (74.5%) 210 (88.6%) 843 (77.6%)  
 GX 36 (4.2%) 9 (3.8%) 45 (4.1%)  
ER (%)     
 Positive 850 (100.0%) 237 (100.0%) 1,087 (100.0%)  
PR (%)    0.02 
 Negative 53 (6.2%) 25 (10.5%) 78 (7.2%)  
 Positive 797 (93.8%) 212 (89.5%) 1,009 (92.8%)  
Type of surgery (%)    0.003 
 Breast conserving 713 (83.9%) 179 (75.5%) 892 (82.1%)  
 Mastectomy 137 (16.1%) 58 (24.5%) 195 (17.9%)  
Treatment arm (%)    0.55 
 2y TAM + 3y ANA 424 (49.9%) 113 (47.7%) 537 (49.4%)  
 5y TAM 426 (50.1%) 124 (52.3%) 550 (50.6%)  

Time to DR analyses in the study cohort

Within the 10-year follow-up period, 86 of 1,087 (7.9%) patients of the study population had developed distant metastases, 48 of 796 (6.0%) node-negative patients, and 38 of 291 (13.1%) patients with 1 to 3 positive nodes.

The estimated probability of 10-year DR as continuous functions of the OncoMasTR Molecular Score (A) and the OncoMasTR Risk Score (B) are shown in Fig. 2.

Figure 2.

Estimated DR risk curves for the OncoMasTR Molecular Score (A) and the OncoMasTR Risk Score (B). The blue area along the curves indicates the 95% CI. The blue histogram in the background shows the distribution of scores for the patients and the vertical line represents the threshold of 5 for low and high risk.

Figure 2.

Estimated DR risk curves for the OncoMasTR Molecular Score (A) and the OncoMasTR Risk Score (B). The blue area along the curves indicates the 95% CI. The blue histogram in the background shows the distribution of scores for the patients and the vertical line represents the threshold of 5 for low and high risk.

Close modal

At 10 years, the DR rate for all patients in the low-risk group was 5.8% versus 21.1% for patients in the high-risk group (HR 4.33, 95% CI, 2.83–6.61, P < 0.0001; absolute risk reduction 15.3%; Table 2; Fig. 3A). The absolute risk reduction at 10 years was 12.9% in patients that were lymph node negative (HR 4.19, 95% CI, 2.35–7.48, P < 0.0001; Table 2; Fig. 3B) and 14.7% in patients with 1 to 3 positive lymph nodes, respectively (HR 3.15, 95% CI, 1.61–6.15, P = 0.0005; Table 2; Fig. 3C).

Table 2.

Univariable Cox models.

Likelihood ratio test
PopulationCovariateNEventsHR (95% CI)χ2 test statisticP valueC-Index
0–3 nodes and 0–10 years OncoMasTR Molecular Scorea 1,087 86 1.73 (1.45–2.06) 36.7 <0.0001 0.70 
 OncoMasTR Risk Score 1,087 86 1.91 (1.62–2.26) 58.6 <0.0001 0.73 
 OncoMasTR Risk Category 1,087 86 4.33 (2.83–6.61) 42.8 <0.0001 0.66 
 
0–3 nodes and 0–5 years OncoMasTR Molecular Scorea 1,087 53 1.80 (1.44–2.25) 26.2 <0.0001 0.71 
 OncoMasTR Risk Score 1,087 53 1.98 (1.60–2.45) 40.1 <0.0001 0.74 
 OncoMasTR Risk Category 1,087 53 4.34 (2.53–7.45) 26.9 <0.0001 0.66 
 
0–3 nodes and 5–10 years OncoMasTR Molecular Scorea 934 33 1.62 (1.22–2.15) 10.8 0.0010 0.68 
 OncoMasTR Risk Score 934 33 1.81 (1.39–2.38) 18.7 <0.0001 0.71 
 OncoMasTR Risk Category 934 33 4.30 (2.17–8.52) 16.0 <0.0001 0.67 
 
0 nodes and 0–10 years OncoMasTR Molecular Scorea 796 48 1.71 (1.36–2.14) 20.2 <0.0001 0.71 
 OncoMasTR Risk Score 796 48 1.79 (1.43–2.24) 25.2 <0.0001 0.72 
 OncoMasTR Risk Category 796 48 4.19 (2.35–7.48) 20.0 <0.0001 0.63 
 
1–3 nodes and 0–10 years OncoMasTR Molecular Scorea 291 38 1.75 (1.31–2.34) 14.7 0.0001 0.69 
 OncoMasTR Risk Score 291 38 1.93 (1.44–2.58) 20.1 <0.0001 0.71 
 OncoMasTR Risk Category 291 38 3.15 (1.61–6.15) 12.0 0.0005 0.65 
Likelihood ratio test
PopulationCovariateNEventsHR (95% CI)χ2 test statisticP valueC-Index
0–3 nodes and 0–10 years OncoMasTR Molecular Scorea 1,087 86 1.73 (1.45–2.06) 36.7 <0.0001 0.70 
 OncoMasTR Risk Score 1,087 86 1.91 (1.62–2.26) 58.6 <0.0001 0.73 
 OncoMasTR Risk Category 1,087 86 4.33 (2.83–6.61) 42.8 <0.0001 0.66 
 
0–3 nodes and 0–5 years OncoMasTR Molecular Scorea 1,087 53 1.80 (1.44–2.25) 26.2 <0.0001 0.71 
 OncoMasTR Risk Score 1,087 53 1.98 (1.60–2.45) 40.1 <0.0001 0.74 
 OncoMasTR Risk Category 1,087 53 4.34 (2.53–7.45) 26.9 <0.0001 0.66 
 
0–3 nodes and 5–10 years OncoMasTR Molecular Scorea 934 33 1.62 (1.22–2.15) 10.8 0.0010 0.68 
 OncoMasTR Risk Score 934 33 1.81 (1.39–2.38) 18.7 <0.0001 0.71 
 OncoMasTR Risk Category 934 33 4.30 (2.17–8.52) 16.0 <0.0001 0.67 
 
0 nodes and 0–10 years OncoMasTR Molecular Scorea 796 48 1.71 (1.36–2.14) 20.2 <0.0001 0.71 
 OncoMasTR Risk Score 796 48 1.79 (1.43–2.24) 25.2 <0.0001 0.72 
 OncoMasTR Risk Category 796 48 4.19 (2.35–7.48) 20.0 <0.0001 0.63 
 
1–3 nodes and 0–10 years OncoMasTR Molecular Scorea 291 38 1.75 (1.31–2.34) 14.7 0.0001 0.69 
 OncoMasTR Risk Score 291 38 1.93 (1.44–2.58) 20.1 <0.0001 0.71 
 OncoMasTR Risk Category 291 38 3.15 (1.61–6.15) 12.0 0.0005 0.65 

aFor OncoMasTR Molecular Score, the HR depicts the hazard change for a 10-point increase.

Figure 3.

Kaplan–Meier plots for time to DR according to the OncoMasTR Risk Category in the overall population (A), in patients who were node negative (B), and in patients with 1–3 positive nodes (C).

Figure 3.

Kaplan–Meier plots for time to DR according to the OncoMasTR Risk Category in the overall population (A), in patients who were node negative (B), and in patients with 1–3 positive nodes (C).

Close modal

In univariable analyses, the OncoMasTR Risk Score was significantly associated with time to DR in all patients (HR 1.91, 95% CI, 1.62–2.26, P < 0.0001; Table 2), in lymph node–negative patients (HR 1.79, 95% CI, 1.43–2.24, P < 0.0001; Table 2), in patients with 1 to 3 positive nodes (HR 1.93, 95% CI, 1.44–2.58, P < 0.0001; Table 2), in years 0 to 5 (HR 1.98, 95% CI, 1.60–2.45, P < 0.0001; Table 2), and in years 5 to 10 (HR 1.81, 95% CI, 1.39–2.38, P < 0.0001; Table 2). Similar significantly prognostic results were obtained for the OncoMasTR Molecular Score and the OncoMasTR Risk Category (Table 2). In addition, we conducted univariate analyses of the three scores using competing risk models and obtained similar results. We also tested the proportional hazards assumption and none of the assumptions for the OncoMasTR variables were violated.

Next in a series of bivariable models, we evaluated whether the OncoMasTR Scores provided additional prognostic value when added to existing prognostic parameters. In years 0 to 10, OncoMasTR Molecular Score, OncoMasTR Risk Score, and OncoMasTR Risk Category added significant prognostic information to age, tumor size, stage, grading, Ki67, NPI, and CTS in the overall population (Supplementary Fig. S1). Again, we observed similar significantly prognostic results in node-negative and in node-positive patients.

We performed three independent multivariable Cox models combining age, tumor size, nodal status, tumor grade, ER expression, PR expression, Ki67 expression, and treatment arm with either the OncoMasTR Molecular Score, the OncoMasTR Risk Score, or the OncoMasTR Risk Category, respectively (Table 3). Each of the three OncoMasTR scores were significantly prognostic in these three multivariable models. In the first multivariable model, tumor size, nodal status, and OncoMasTR Molecular Score (HR 1.49, 95% CI, 1.20–1.86, P = 0.0004) were independent prognostic factors for time to DR (Table 3). In the second multivariable model, only tumor size and OncoMasTR Risk Score (HR 1.58, 95% CI, 1.26–1.97, P < 0.0001) remained independently associated with time to DR (Table 3). In the third multivariable model, tumor size, nodal status, Ki67 expression, and OncoMasTR Risk Category (HR 2.37, 95% CI, 1.41–3.98, P = 0.001) were highly prognostic for prediction of DR (Table 3).

Table 3.

Multivariable Cox models.

Model with OncoMasTRModel with OncoMasTRModel with OncoMasTR
Molecular ScoreRisk ScoreRisk Category
CovariateHR (95% CI)LR χ2PHR (95% CI)LR χ2PHR (95% CI)LR χ2P
Agea (10 years) 1.04 (0.79–1.38) 0.1 0.79 1.03 (0.78–1.37) 0.0 0.83 1.06 (0.80–1.40) 0.2 0.69 
Tumor sizea (cm) 1.40 (1.19–1.63) 12.7 0.0004 1.31 (1.10–1.56) 7.1 0.008 1.37 (1.16–1.61) 10.1 0.002 
Nodal status (pN1 vs. pN0) 2.06 (1.33–3.19) 10.1 0.002 1.46 (0.92–2.32) 2.5 0.11 1.70 (1.08–2.67) 5.2 0.02 
Grading (G2 vs. G1) 1.15 (0.57–2.34) 0.2 0.69 1.13 (0.56–2.29) 0.1 0.74 1.12 (0.55–2.27) 0.1 0.76 
ERa (10%) 0.96 (0.81–1.14) 0.2 0.65 0.96 (0.81–1.14) 0.2 0.63 0.95 (0.80–1.13) 0.3 0.59 
PRa (10%) 0.98 (0.92–1.05) 0.2 0.64 0.99 (0.92–1.06) 0.2 0.70 0.98 (0.91–1.05) 0.4 0.50 
Ki67a (10%) 1.27 (0.97–1.66) 2.9 0.09 1.23 (0.94–1.60) 2.1 0.14 1.37 (1.06–1.76) 5.5 0.02 
Treatment (TAM/ANA vs. TAM) 0.85 (0.55–1.31) 0.6 0.46 0.84 (0.55–1.30) 0.6 0.44 0.83 (0.54–1.28) 0.7 0.41 
OncoMasTRa 1.49 (1.20–1.86) 12.7 0.0004 1.58 (1.26–1.97) 16.0 <0.0001 2.37 (1.41–3.98) 10.4 0.001 
Model with OncoMasTRModel with OncoMasTRModel with OncoMasTR
Molecular ScoreRisk ScoreRisk Category
CovariateHR (95% CI)LR χ2PHR (95% CI)LR χ2PHR (95% CI)LR χ2P
Agea (10 years) 1.04 (0.79–1.38) 0.1 0.79 1.03 (0.78–1.37) 0.0 0.83 1.06 (0.80–1.40) 0.2 0.69 
Tumor sizea (cm) 1.40 (1.19–1.63) 12.7 0.0004 1.31 (1.10–1.56) 7.1 0.008 1.37 (1.16–1.61) 10.1 0.002 
Nodal status (pN1 vs. pN0) 2.06 (1.33–3.19) 10.1 0.002 1.46 (0.92–2.32) 2.5 0.11 1.70 (1.08–2.67) 5.2 0.02 
Grading (G2 vs. G1) 1.15 (0.57–2.34) 0.2 0.69 1.13 (0.56–2.29) 0.1 0.74 1.12 (0.55–2.27) 0.1 0.76 
ERa (10%) 0.96 (0.81–1.14) 0.2 0.65 0.96 (0.81–1.14) 0.2 0.63 0.95 (0.80–1.13) 0.3 0.59 
PRa (10%) 0.98 (0.92–1.05) 0.2 0.64 0.99 (0.92–1.06) 0.2 0.70 0.98 (0.91–1.05) 0.4 0.50 
Ki67a (10%) 1.27 (0.97–1.66) 2.9 0.09 1.23 (0.94–1.60) 2.1 0.14 1.37 (1.06–1.76) 5.5 0.02 
Treatment (TAM/ANA vs. TAM) 0.85 (0.55–1.31) 0.6 0.46 0.84 (0.55–1.30) 0.6 0.44 0.83 (0.54–1.28) 0.7 0.41 
OncoMasTRa 1.49 (1.20–1.86) 12.7 0.0004 1.58 (1.26–1.97) 16.0 <0.0001 2.37 (1.41–3.98) 10.4 0.001 

Note: 1,083 patients (0–3 nodes and 0–10 years) with 85 DR events for whom all covariates were available were included.

aAge, tumor size, ER, PR, Ki67, and OncoMasTR Molecular Score were included as continuous variables. For age, the HR depicts the hazard change for a 10-year increase. For tumor size, the HR depicts the hazard change for a 1-cm increase. For ER, PR, and Ki67, the HR depicts the hazard change for a 10% increase. For OncoMasTR Molecular Score, the HR depicts the hazard change for a 10-point increase.

The OncoMasTR gene signature was developed with the aim of accurately estimating the risk of DR in the target population: patients with ER-positive, HER2-negative, early-stage breast cancer with 0 to 3 positive lymph nodes. In contrast to currently available gene signatures in this setting, which were discovered using gene-expression profiling, OncoMasTR was trained using a novel bioinformatic approach that identified a shared network of MTRs from two already validated prognostic gene signatures (8, 9). In this approach, FOXM1, PTTG1, and ZNF367 have been identified to play a key role in breast cancer biology and progression and formed the OncoMasTR Molecular Score (8, 11). The OncoMasTR Risk Score combines the molecular score with tumor size and lymph node status and the OncoMasTR Risk Category classifies patients into the two clinically relevant low- or high-risk groups (11).

The first independent validation of the OncoMasTR Risk Score was carried out in 646 patients from the TransATAC cohort (11). That dataset included postmenopausal women with ER-positive, HER2-negative, breast cancer who received 5 years of endocrine therapy. In the TransATAC study, the OncoMasTR Risk Score was highly prognostic for prediction of DR and provided additional prognostic information beyond NPI and CTS in the overall cohort as well as in node-negative patients. In patients with 1 to 3 positive lymph nodes in that cohort, the OncoMasTR Risk Score was significantly prognostic on its own but did not provide significant additional prognostic value beyond NPI and CTS, possibly due to the relatively small number of patients in this group (11).

The performance of the OncoMasTR test has also been evaluated in a subset of Irish patients enrolled in the TAILORx study (n = 404), showing that OncoMasTR was highly prognostic for DR in this cohort (25). There was only a moderate correlation between OncoMasTR and Oncotype DX: OncoMasTR reclassified a majority of Oncotype DX “intermediate risk” patients to “low risk” and was particularly accurate in reclassifying patients with distant recurrences to “high risk” (25).

In the present study, we report the second independent validation of the clinical performance of the OncoMasTR Risk Score in a large cohort of patients who received endocrine therapy only and who were enrolled into ABCSG Trial 8, a prospective, randomized clinical trial with long-term follow-up. Similar to the TransATAC results, all three OncoMasTR scores (molecular score, risk score, and risk category) significantly predicted DR in year 0 to 10 in all patients, and in patients who were node negative. When compared with the TransATAC results, the OncoMasTR scores show improved performance for the prediction of DR in patients with 1 to 3 affected lymph nodes. Moreover, all three OncoMasTR scores provided significant additional prognostic information when added to standard prognostic variables, NPI, and CTS. As a result, together with the data from TransATAC, there is now level 1B evidence for the OncoMasTR Risk Score in predicting DR.

The strengths of our study are a large, well-documented patient cohort representing the precise clinical target population of the test, the use of analytically validated assays, prospectively defined scores and thresholds, and the blindness of laboratory staff to clinical data. Also, all statistical analyses were carried out according to a predefined statistical analysis plan.

Limitations of this work are the retrospective character and the fact that patients based on the ABCSG Trial 8 inclusion criteria did not receive adjuvant chemotherapy. The latter makes it even theoretically impossible to potentially predict chemotherapy benefit in this cohort. However, if a molecular risk prediction tool accurately defines low risk of distant relapse (as it is the case with the current results: less than 5% DR at 10 years for patients who were node negative, and less than 10% for 1–3 positive nodes), it is highly unlikely that any additional adjuvant therapy (including cytotoxic chemotherapy) could improve those outcomes with reasonable risk/benefit ratio.

The prognostic and/or predictive value of the OncoMasTR Risk Score in other breast cancer subsets such as premenopausal patients with breast cancer will have to be determined in future studies.

Several other gene signatures such as EndoPredict, Prosigna, and MammaPrint were also validated in ABCSG Trial 8 (5, 7, 15, 16, 26, 27). Therefore, we can numerically compare the performance of the OncoMasTR test with other tests. The OncoMasTR test and MammaPrint classified 78% and 71% of ABCSG Trial 8 patients with 0 to 3 involved lymph nodes as “low risk”, respectively (26). In the combined analysis of ABCSG Trial 6 and 8 and in the combined analysis of patients from ABCSG Trial 8 and ATAC, EndoPredict and Prosigna categorized 72% and 58% of the patients with 0 to 3 involved lymph nodes in the low-risk group, respectively (16, 27).

In addition, we performed a cross-study analysis that allows for general comparisons between Prosigna and OncoMasTR in terms of performance. Supplementary Table S2 shows results for patients with ER-positive, HER2-negative breast cancer from ABCSG Trial 8. The patient populations with Prosigna and OncoMasTR test data are not completely identical but, for 1,030 of 1,397 (74%) patients, both Prosigna and OncoMasTR results are available. In contrast to the OncoMasTR subset, which only includes patients with 0 to 3 nodes, the Prosigna subset includes patients with 0 to 9 nodes (7). Since the study endpoints differ in both studies [DRFS (distant recurrences and deaths due to breast cancer) for Prosigna, time to DR for OncoMasTR], we also analyzed DFRS for OncoMasTR. C-Indices are for the continuous scores and are almost identical for both tests. The comparison shown in Supplementary Table S2 demonstrates that both Prosigna and OncoMasTR performed similarly well in ABSCG Trial 8 with regard to distinguishing patients at “high risk” from patients at “low risk”.

However, the present results suggest that the OncoMasTR test classifies a large proportion of patients as “low risk”; the ultimate goal of obtaining additional molecular prognostic information is to identify as many patients with DR risk levels low enough to argue to spare them adjuvant chemotherapy. A formal comparison of the OncoMasTR test with other molecular tests is planned by ABCSG and the next logical step in the development of this signature.

In summary, the OncoMasTR Risk Score is highly prognostic for DR in postmenopausal women with ER-positive, HER2-negative primary breast cancer with 0 to 3 involved lymph nodes. This independent validation in ABCSG Trial 8 provides level 1B evidence for the prognostic role of the OncoMasTR Risk Score.

M. Filipits reports grants from OncoMark Ltd. and European Union during the conduct of the study, as well as personal fees from AstraZeneca, Bayer, Biomedica, Boehringer Ingelheim, Eli Lilly, Merck Sharp & Dohme, Myriad Genetics Inc., Pfizer, and Roche outside the submitted work. C.F. Singer reports grants, non-financial support, and other support from Novartis; grants, personal fees, and non-financial support from Amgen; grants and personal fees from AstraZeneca; personal fees and non-financial support from Roche; and grants from Myriad outside the submitted work. F. Fitzal reports personal fees from Novartis, Pfizer, Roche, and AstraZeneca outside the submitted work. Z. Bago-Horvath reports personal fees from Roche, as well as other support from MSD outside the submitted work. R. Greil reports personal fees and other support from AstraZeneca during the conduct of the study, as well as personal fees and other support from Roche, Amgen, Janssen, AstraZeneca, Novartis, MSD, Celgene, Gilead, BMS, AbbVie, Daiichi Sankyo, and Sanofi outside the submitted work. M. Balic reports grants and other support from Medical University of Graz during the conduct of the study. P. Regitnig reports personal fees and non-financial support from Roche and Novartis, as well as personal fees from Diaceutics outside the submitted work. W. Hulla reports personal fees from Biomedica outside the submitted work. D. Egle reports personal fees from AstraZeneca, Daichii Sankyo, MSD, Novartis, and Seagen, as well as personal fees and non-financial support from Pfizer, Roche, and Lilly outside the submitted work. S. Barron was an employee of OncoMark Limited. T. Loughman reports personal fees from OncoMark Limited during the conduct of the study. D. O'Leary reports personal fees from OncoMark Limited outside the submitted work. W.M. Gallagher reports other support from OncoMark Limited during the conduct of the study, as well as other support from Carrick Therapeutics outside the submitted work; in addition, W.M. Gallagher has a patent for EP3194621B1 issued and licensed to OncoMark Limited. D. Hlauschek reports grants from Oncomark Limited during the conduct of the study. M. Gnant reports personal fees from Amgen, Daiichi Sankyo, AstraZeneca, Eli Lilly, LifeBrain, Nanostring, Novartis, and TLC Biopharmaceuticals outside the submitted work; in addition, an immediate family member of M. Gnant is employed by Sandoz. P. Dubsky reports grants from Oncomark Limited during the conduct of the study; P. Dubsky also reports grants from Cepheid/Danaher, Agendia, and Myriad, as well as other support from Myriad, Merck, AstraZeneca, and Roche outside the submitted work. No disclosures were reported by the other authors.

M. Filipits: Conceptualization, resources, data curation, formal analysis, supervision, funding acquisition, validation, investigation, methodology, writing–original draft, project administration, writing–review and editing. M. Rudas: Investigation, methodology, writing–review and editing. V. Kainz: Investigation, methodology, writing–review and editing. C.F. Singer: Resources, investigation, writing–review and editing. F. Fitzal: Resources, investigation, writing–review and editing. Z. Bago-Horvath: Resources, investigation, writing–review and editing. R. Greil: Resources, investigation, writing–review and editing. M. Balic: Resources, investigation, writing–review and editing. P. Regitnig: Resources, investigation, writing–review and editing. S. Halper: Resources, investigation, writing–review and editing. W. Hulla: Resources, investigation, writing–review and editing. D. Egle: Resources, investigation, writing–review and editing. S. Barron: Conceptualization, formal analysis, writing–review and editing. T. Loughman: Conceptualization, resources, investigation, methodology, writing–review and editing. D. O'Leary: Conceptualization, resources, methodology, writing–review and editing. W.M. Gallagher: Conceptualization, investigation, methodology, writing–review and editing. D. Hlauschek: Conceptualization, formal analysis, writing–review and editing. M. Gnant: Resources, investigation, writing–original draft, writing–review and editing. P. Dubsky: Conceptualization, data curation, formal analysis, investigation, methodology, project administration, writing–review and editing.

This work was supported by OncoMark Ltd. The work was supported by European Union's Horizon 2020 research and innovation program awarded to OncoMark under grant agreement number 698630.

We thank our patients who contributed to ABCSG Trial 8, the ABCSG investigators, ABCSG centers, and the ABCSG central office for ongoing support. We thank Jana Link and Magdalena Schwarz for assistance in project coordination, Victoria Kügler and Katharina Hofbauer for technical assistance, and Chan-Ju Angel Wang for administrative assistance.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
van de Vijver
MJ
,
He
YD
,
van't Veer
LJ
,
Dai
H
,
Hart
AAM
,
Voskuil
DW
, et al
A gene-expression signature as a predictor of survival in breast cancer
.
N Engl J Med
2002
;
347
:
1999
2009
.
2.
van't Veer
LJ
,
Dai
H
,
van de Vijver
MJ
,
He
YD
,
Hart
AAM
,
Mao
M
, et al
Gene expression profiling predicts clinical outcome of breast cancer
.
Nature
2002
;
415
:
530
6
.
3.
Paik
S
,
Shak
S
,
Tang
G
,
Kim
C
,
Baker
J
,
Cronin
M
, et al
A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer
.
N Engl J Med
2004
;
351
:
2817
26
.
4.
Ma
X-J
,
Salunga
R
,
Dahiya
S
,
Wang
W
,
Carney
E
,
Durbecq
V
, et al
A five-gene molecular grade index and HOXB13:IL17BR are complementary prognostic factors in early stage breast cancer
.
Clin Cancer Res
2008
;
14
:
2601
8
.
5.
Filipits
M
,
Rudas
M
,
Jakesz
R
,
Dubsky
P
,
Fitzal
F
,
Singer
CF
, et al
A new molecular predictor of distant recurrence in ER-positive, HER2-negative breast cancer adds independent information to conventional clinical risk factors
.
Clin Cancer Res
2011
;
17
:
6012
20
.
6.
Jankowitz
RC
,
Cooper
K
,
Erlander
MG
,
Ma
X-J
,
Kesty
NC
,
Li
H
, et al
Prognostic utility of the breast cancer index and comparison to Adjuvant! Online in a clinical case series of early breast cancer
.
Breast Cancer Res
2011
;
13
:
R98
.
7.
Gnant
M
,
Filipits
M
,
Greil
R
,
Stoeger
H
,
Rudas
M
,
Bago-Horvath
Z
, et al
Predicting distant recurrence in receptor-positive breast cancer patients with limited clinicopathological risk: using the PAM50 Risk of Recurrence score in 1478 postmenopausal patients of the ABCSG-8 trial treated with adjuvant endocrine therapy alone
.
Ann Oncol
2014
;
25
:
339
45
.
8.
Lanigan
F
,
Brien
GL
,
Fan
Y
,
Madden
SF
,
Jerman
E
,
Maratha
A
, et al
Delineating transcriptional networks of prognostic gene signatures refines treatment recommendations for lymph node-negative breast cancer patients
.
FEBS J
2015
;
282
:
3455
73
.
9.
Moran
B
,
Rahman
A
,
Palonen
K
,
Lanigan
FT
,
Gallagher
WM
. 
Master transcriptional regulators in cancer: discovery via reverse engineering approaches and subsequent validation
.
Cancer Res
2017
;
77
:
2186
90
.
10.
Loughman
T
,
Chan-Ju Wang
A
,
Dynoodt
P
,
Fender
B
,
Lopez Ruiz
C
,
Barron
S
, et al
Analytical validation of OncoMasTR, a multigene test for predicting risk of distant recurrence in hormone receptor-positive early stage breast cancer
.
Ann Oncol
2018
;
29
:
viii65
.
11.
Buus
R
,
Sestak
I
,
Barron
S
,
Loughman
T
,
Fender
B
,
Ruiz
CL
, et al
Validation of the OncoMasTR risk score in estrogen receptor-positive/HER2-negative patients: a TransATAC study
.
Clin Cancer Res
2020
;
26
:
623
31
.
12.
Jakesz
R
,
Jonat
W
,
Gnant
M
,
Mittlboeck
M
,
Greil
R
,
Tausch
C
, et al
Switching of postmenopausal women with endocrine-responsive early breast cancer to anastrozole after 2 years' adjuvant tamoxifen: combined results of ABCSG trial 8 and ARNO 95 trial
.
Lancet
2005
;
366
:
455
62
.
13.
Dubsky
PC
,
Jakesz
R
,
Mlineritsch
B
,
Pöstlberger
S
,
Samonigg
H
,
Kwasny
W
, et al
Tamoxifen and anastrozole as a sequencing strategy: a randomized controlled trial in postmenopausal patients with endocrine-responsive early breast cancer from the Austrian Breast and Colorectal Cancer Study Group
.
J Clin Oncol
2012
;
30
:
722
8
.
14.
Bago-Horvath
Z
,
Rudas
M
,
Dubsky
P
,
Jakesz
R
,
Singer
CF
,
Kemmerling
R
, et al
Adjuvant sequencing of tamoxifen and anastrozole is superior to tamoxifen alone in postmenopausal women with low proliferating breast cancer
.
Clin Cancer Res
2011
;
17
:
7828
34
.
15.
Filipits
M
,
Nielsen
TO
,
Rudas
M
,
Greil
R
,
Stöger
H
,
Jakesz
R
, et al
The PAM50 risk-of-recurrence score predicts risk for late distant recurrence after endocrine therapy in postmenopausal women with endocrine-responsive early breast cancer
.
Clin Cancer Res
2014
;
20
:
1298
305
.
16.
Filipits
M
,
Dubsky
P
,
Rudas
M
,
Greil
R
,
Balic
M
,
Bago-Horvath
Z
, et al
Prediction of distant recurrence using endopredict among women with ER(+), HER2(-) node-positive and node-negative breast cancer treated with endocrine therapy only
.
Clin Cancer Res
2019
;
25
:
3865
72
.
17.
Nielsen
TO
,
Leung
SCY
,
Rimm
DL
,
Dodson
A
,
Acs
B
,
Badve
S
, et al
Assessment of Ki67 in breast cancer: updated recommendations from the international Ki67 in breast cancer working group
.
J Natl Cancer Inst
2020
;
113
:
808
19
.
18.
Sestak
I
,
Filipits
M
,
Buus
R
,
Rudas
M
,
Balic
M
,
Knauer
M
, et al
Prognostic value of EndoPredict in women with hormone receptor-positive, HER2-negative invasive lobular breast cancer
.
Clin Cancer Res
2020
;
26
:
4682
7
.
19.
Haybittle
JL
,
Blamey
RW
,
Elston
CW
,
Johnson
J
,
Doyle
PJ
,
Campbell
FC
, et al
A prognostic index in primary breast cancer
.
Br J Cancer
1982
;
45
:
361
6
.
20.
Galea
MH
,
Blamey
RW
,
Elston
CE
,
Ellis
IO
. 
The nottingham prognostic index in primary breast cancer
.
Breast Cancer Res Treat
1992
;
22
:
207
19
.
21.
Cuzick
J
,
Dowsett
M
,
Pineda
S
,
Wale
C
,
Salter
J
,
Quinn
E
, et al
Prognostic value of a combined estrogen receptor, progesterone receptor, Ki-67, and human epidermal growth factor receptor 2 immunohistochemical score and comparison with the Genomic Health recurrence score in early breast cancer
.
J Clin Oncol
2011
;
29
:
4273
8
.
22.
Dowsett
M
,
Sestak
I
,
Lopez-Knowles
E
,
Sidhu
K
,
Dunbier
AK
,
Cowens
JW
, et al
Comparison of PAM50 risk of recurrence score with oncotype DX and IHC4 for predicting risk of distant recurrence after endocrine therapy
.
J Clin Oncol
2013
;
31
:
2783
90
.
23.
Royston
P
,
Parmar
MK
. 
Flexible parametric proportional-hazards and proportional-odds models for censored survival data, with application to prognostic modelling and estimation of treatment effects
.
Stat Med
2002
;
21
:
2175
97
.
24.
McShane
LM
,
Altman
DG
,
Sauerbrei
W
,
Taube
SE
,
Gion
M
,
Clark
GM
. 
Reporting recommendations for tumor marker prognostic studies
.
J Clin Oncol
2005
;
23
:
9067
72
.
25.
O'Connor
D
,
Kelly
CM
,
Crown
J
,
Russell
N
,
Barron
S
,
Loughman
T
, et al
Additional prognostic value of OncoMasTR multigene prognostic signature to clinicopathological information in patients with HR-positive, HER2-negative, lymph node-negative breast cancer from the TAILORx Tissue Bank, Ireland
.
J Clin Oncol
2019
;
37
:
535
.
26.
Dubsky
P
,
Van't Veer
L
,
Gnant
M
,
Rudas
M
,
Bago-Horvath
Z
,
Greil
R
, et al
A clinical validation study of MammaPrint in hormone receptor-positive breast cancer from the Austrian Breast and Colorectal Cancer Study Group 8 (ABCSG-8) biomarker cohort
.
ESMO Open
2021
;
6
:
100006
.
27.
Sestak
I
,
Cuzick
J
,
Dowsett
M
,
Lopez-Knowles
E
,
Filipits
M
,
Dubsky
P
, et al
Prediction of late distant recurrence after 5 years of endocrine treatment: a combined analysis of patients from the Austrian breast and colorectal cancer study group 8 and arimidex, tamoxifen alone or in combination randomized trials using the PAM50 risk of recurrence score
.
J Clin Oncol
2015
;
33
:
916
22
.