- Research article
- Open Access
Prediction and clinical utility of a contralateral breast cancer risk model
Breast Cancer Research volume 21, Article number: 144 (2019)
Breast cancer survivors are at risk for contralateral breast cancer (CBC), with the consequent burden of further treatment and potentially less favorable prognosis. We aimed to develop and validate a CBC risk prediction model and evaluate its applicability for clinical decision-making.
We included data of 132,756 invasive non-metastatic breast cancer patients from 20 studies with 4682 CBC events and a median follow-up of 8.8 years. We developed a multivariable Fine and Gray prediction model (PredictCBC-1A) including patient, primary tumor, and treatment characteristics and BRCA1/2 germline mutation status, accounting for the competing risks of death and distant metastasis. We also developed a model without BRCA1/2 mutation status (PredictCBC-1B) since this information was available for only 6% of patients and is routinely unavailable in the general breast cancer population. Prediction performance was evaluated using calibration and discrimination, calculated by a time-dependent area under the curve (AUC) at 5 and 10 years after diagnosis of primary breast cancer, and an internal-external cross-validation procedure. Decision curve analysis was performed to evaluate the net benefit of the model to quantify clinical utility.
In the multivariable model, BRCA1/2 germline mutation status, family history, and systemic adjuvant treatment showed the strongest associations with CBC risk. The AUC of PredictCBC-1A was 0.63 (95% prediction interval (PI) at 5 years, 0.52–0.74; at 10 years, 0.53–0.72). Calibration-in-the-large was -0.13 (95% PI: -1.62–1.37), and the calibration slope was 0.90 (95% PI: 0.73–1.08). The AUC of Predict-1B at 10 years was 0.59 (95% PI: 0.52–0.66); calibration was slightly lower. Decision curve analysis for preventive contralateral mastectomy showed potential clinical utility of PredictCBC-1A between thresholds of 4–10% 10-year CBC risk for BRCA1/2 mutation carriers and non-carriers.
We developed a reasonably calibrated model to predict the risk of CBC in women of European-descent; however, prediction accuracy was moderate. Our model shows potential for improved risk counseling, but decision-making regarding contralateral preventive mastectomy, especially in the general breast cancer population where limited information of the mutation status in BRCA1/2 is available, remains challenging.
Breast cancer (BC) is a major burden for women’s health . Survival has improved substantially over the past half century due to earlier detection and advanced treatment modalities, for example in the Netherlands, 10-year survival of a first primary BC improved from 40% in 1961–1970 to 79% in 2006–2010 . Consequently, increasing numbers of BC survivors are at risk to develop a new primary tumor in the opposite (contralateral) breast, with subsequent treatment and potentially less favorable prognosis . BC survivors are more likely to develop contralateral breast cancer (CBC) compared to healthy women to develop a first primary BC .
Women at elevated CBC risk have been identified to be BRCA1/2 and CHEK2 c.1100del mutation carriers and to have a BC family history, particularly a family history of bilateral BC [5,6,7,8,9,10]. For BRCA1/2 mutation carriers, in whom CBC risk is high, contralateral preventive mastectomy (CPM) is often offered . However, the average risk of CBC among all first BC survivors is still relatively low, with an incidence of ~ 0.4% per year [12,13,14]. Despite this, in recent years, CPM frequency has increased among women in whom CBC risk is low . For these reasons, there is an urgent need for improved individualized prediction of CBC risk, both to facilitate shared decision-making of physicians and women regarding treatment and prevention strategies for those at high CBC risk and to avoid unnecessary CPM or surveillance mammography after first primary BC when CBC risk is low.
To our knowledge, only one specific CBC risk prediction model (CBCrisk) has been developed to date. CBCrisk used data on 1921 CBC cases and 5763 matched controls with validation in two independent US studies containing a mix of invasive and in situ BC [16, 17]. Moreover, the level of prediction performance measures such as calibration and discrimination needed for a CBC risk prediction to be clinically useful have not yet been addressed .
Our aim was twofold: first, to develop and validate a CBC risk prediction model using a large international series of individual patient data including 132,756 patients with a first primary invasive BC between 1990 and 2013 from multiple studies in Europe, USA, and Australia with 4682 incident CBCs, and second, to evaluate the potential clinical utility of the model to support decision-making.
Material and methods
We used data from five main sources: three studies from the Netherlands, 16 studies from the Breast Cancer Association Consortium (BCAC), and a cohort from the Netherlands Cancer Registry [19,20,21,22]. For details regarding data collection and patient inclusion, see Additional file 1: Data and patient selection and Table S1, and Additional file 1: Table S2. We included female patients with invasive non-metastatic first primary BC with no prior history of cancer (except for non-melanoma skin cancer). The studies were either population- or hospital-based series; most women were of European-descent. We only included women diagnosed after 1990 to have a population with diagnostic and treatment procedures likely close to modern practice and at the same time sufficient follow-up to study CBC incidence; in total 132,756 women from 20 studies were included. All studies were approved by the appropriate ethics and scientific review boards. All women provided written informed consent or did not object to secondary use of clinical data in accordance with Dutch legislation and codes of conduct [23, 24].
Available data and variable selection
Several factors have been shown or suggested to be associated with CBC risk, including age at first BC, family history for BC, BRCA1/2 and CHEK2 c.1100del mutations, body mass index (BMI), breast density change, (neo)adjuvant chemotherapy, endocrine therapy, CPM, and characteristics of the first BC such as histology (lobular vs ductal), estrogen receptor (ER) status, lymph node status, tumor size, and TNM stage [5, 9, 12, 25,26,27,28,29,30,31,32,33,34,35,36]. The choice of factors to include in the analyses was determined by evidence from literature, availability of data in the cohorts, and current availability in clinical practice. We extracted the following information: BRCA1/2 germline mutation, (first degree) family history of primary BC, and regarding primary BC diagnosis: age, nodal status, size, grade, morphology, ER status, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) status, administration of adjuvant and/or neoadjuvant chemotherapy, adjuvant endocrine therapy, adjuvant trastuzumab therapy, radiotherapy. We excluded PR status and TNM stage of the primary BC due to collinearity with ER status and the size of the primary tumor, respectively. In the current clinical practice, only patients with ER-positive tumors receive endocrine therapy and only patients with HER2-positive tumors receive trastuzumab; these co-occurrences were considered in the model by using composite categorical variables. More information is available online about the factors included in the analyses (Additional file 1: Data patient selection and Additional file 2: Figure S1), follow-up per dataset, and study design (Additional file 1: Table S2).
All analyses were performed using SAS (SAS Institute Inc., Cary, NC, USA) and R software .
Primary endpoint, follow-up, and predictors
The primary endpoint in the analyses was in situ or invasive metachronous CBC. Follow-up started 3 months after invasive first primary BC diagnosis, in order to exclude synchronous CBCs, and ended at date of CBC, distant metastasis (but not at loco-regional relapse), CPM, or last date of follow-up (due to death, being lost to follow-up, or end of study), whichever occurred first. The follow-up of 27,155 (20.4%) women from the BCAC studies, recruited more than 3 months after diagnosis of the first primary BC (prevalent cases), started at recruitment (left truncation). Distant metastasis and death due to any cause were considered as competing events. Patients who underwent CPM during the follow-up were censored because the CBC risk was almost zero after a CPM . Missing data were multiply imputed by chained equations (MICE) to avoid loss of information due to case-wise deletion [39, 40]. Details about the imputation model, strategy used, and the complete case analysis are provided in Additional file 1: Multiple Imputation of missing values, complete case analysis, and model diagnostics and baseline recalibration and Additional file 1: Tables S3 and S4.
Model development and validation
For model development, we used a multivariable Fine and Gray model regression to account for death and distant metastases as competing events [41, 42]. Heterogeneity of baseline risks between studies was taken into account using the study as a stratification term. A stratified model allows the baseline subdistribution hazard to be different across the studies, and parameter estimation is performed by maximization of the partial likelihood per study. A Breslow-type estimator was used to estimate the baseline cumulative subdistribution hazard per study. The assumption of proportional subdistribution hazards was graphically checked using Schoenfield residuals . The resulting subdistributional hazard ratios (sHRs) and corresponding 95% confidence intervals (CI) were pooled from the 10 imputed data sets using Rubin’s rules . We built a nomogram for estimating the 5- and 10-year cumulative incidence of CBC as a graphical representation of the multivariable risk prediction model .
The validity of the model was investigated by leave-one-study-out cross-validation, i.e., in each validation step, all studies are used except one in which the validity of the model is evaluated [46, 47]. Since the ABCS study and some studies from BCAC had insufficient CBC events required for reliable validation, we used the geographic area as unit of splitting. We had 20 studies in five main sources: 17 out of 20 studies that were combined in 4 geographic areas. In total, 3 studies and 4 geographic areas were used to assess the prediction performance of the model (see Additional file 1: Leave-one-study-out cross-validation and Additional file 1: Table S5, [47, 48].
The performance of the model was assessed by discrimination ability to differentiate between patients who experienced CBC and those who did not, and by calibration, which measures the agreement between observed and predicted CBC risk. Discrimination was quantified by time-dependent area under the ROC curves (AUCs) based on Inverse Censoring Probability Weighting at 5 and 10 years [49, 50]. In the presence of competing risks, the R package timeROC provides two types of AUC according to a different definition of time-dependent cases and controls. AUCs were calculated considering a patient who developed a CBC as a case and a patient free of any event as a control at 5 and 10 years . Values of AUCs close to 1 indicate good discriminative ability, while values close to 0.5 indicated poor discriminative ability. Calibration was assessed by the calibration-in-the-large and slope statistic . Calibration-in-the-large lower or higher than 0 indicates that prediction is systematically too high or low, respectively. A calibration slope of 1.0 indicates good overall calibration; slopes below (above) 1.0 indicate over (under) estimation of risk by the model.
To allow for heterogeneity among studies, a random-effect meta-analysis was performed to provide summaries of discrimination and calibration performance. The 95% prediction intervals (PI) indicated the likely range for the prediction performances of the model in a new dataset. Further details about the validation process are provided in Additional file 1: Leave-one-study-out cross-validation.
The clinical utility of the prediction model was evaluated using decision curve analysis (DCA) [52, 53]. Such a decision may apply to more or less intensive screening and follow-up or to decision of a CPM. The key part of the DCA is the net benefit, which is the number of true-positive classifications (in this example: the benefit of CPM to a patient who would have developed a CBC) minus the number of false-positive classifications (in this example: the harm of unnecessary CPM in a patient who would not have developed a CBC). The false positives are weighted by a factor related to the relative harm of a missed CBC versus an unnecessary CPM. The weighting is derived from the threshold probability to develop a CBC using a defined landmark time point (e.g., CBC risk at 5 or 10 years) . For example, a threshold of 10% implies that CPM in 10 patients, of whom one would develop CBC if untreated, is acceptable (thus performing 9 unnecessary CPMs). The net benefit of a prediction model is traditionally compared with the strategies of treat all or treat none. Since the use of CPM is generally only suggested among BRCA1/2 mutation carriers, for a more realistic illustration, the decision curve analysis was reported among BRCA1/2 mutation carriers and non-carriers . See Additional file 1: Clinical utility for details.
A total of 132,756 invasive primary BC women diagnosed between 1990 and 2013, with 4682 CBC events, from 20 studies, were used to derive the model for CBC risk (Additional file 1: Table S2). Median follow-up time was 8.8 years, and CBC cumulative incidences at 5 and 10 years were 2.1% and 4.1%, respectively. Details of the studies and patient, tumor, and treatment characteristics are provided in Additional file 1: Table S6. The multivariable model with estimates for all included factors is shown in Table 1 and Additional file 3. BRCA1/2 germline mutation status, family history, and systemic adjuvant treatment showed the strongest associations with CBC risk.
The prediction performance of the main model (PredictCBC, version 1A) based on the leave-one-study-out cross-validation method is shown in Fig. 1. The AUC at 5 years was 0.63 (95% confidence interval (CI): 0.58–0.67; 95% prediction interval (PI): 0.52–0.74)); the AUC at 10 years was also 0.63 (95% CI: 0.59–0.66; 95% PI: 0.53–0.72). Calibrations showed some indications of overestimation of risk. The calibration-in-the-large was − 0.13 (95% CI: -0.66–0.40; 95% PI: -1.62–1.37). The calibration slope was 0.90 (95% CI: 0.79–1.02; 95% PI: 0.73–1.08) in the cross-validation. Calibration plots are provided in Additional file 2: Figure S2 and S3.
The nomogram representing a graphical tool for estimating the CBC cumulative incidence at 5 and 10 years based on our model and the estimated baseline of the Dutch Cancer Registry is shown in Fig. 2. In the nomogram, the categories of each factor are assigned a score using the topmost “Points” scale, then all scores are summed up to obtain the “Total points”, which relate to the cumulative incidence of CBC. The formulae of the models (PredictCBC-1A and 1B) providing the predicted cumulative incidence are given in Additional file 1: Formula to estimate the CBC risk and formula to estimate CBC risk in patients not tested for BRCA.
The DCAs for preventive contralateral mastectomy showed the potential clinical utility of PredictCBC-1A between thresholds of 4–10% 10-year CBC risk for BRCA1/2 mutation carriers and non-carriers (Table 2 and Additional file 3). For example, if we find it acceptable that one in 10 patients for whom a CPM is recommended develops a CBC, a risk threshold of 10% may be used to define high and low risk BRCA1/2 mutation carriers based on the absolute 10-year CBC risk prediction estimated by the model. Compared with a strategy recommending CPM to all carriers of a mutation in BRCA1/2, this strategy avoids 161 CPMs per 1000 patients. In contrast, almost no non BRCA1/2 mutation carriers reach the 10% threshold (the general BC population, Fig. 3). The decision curves provide a comprehensive overview of the net benefit for a range of harm-benefit thresholds at 10-year CBC risk (Fig. 4).
We also derived a risk prediction model (PredictCBC, version 1B) omitting BRCA status to provide CBC risk estimates for first BC patients not tested for BRCA1/2 mutations. This model has slightly lower prediction performance; AUC at 5 and 10 years was both 0.59 (at 5 years: 95% CI: 0.54–0.63, 95% PI: 0.46–0.71; at 10 years: 95% CI: 0.56–0.62, 95% PI: 0.52–0.66), calibration-in-the-large was − 0.17 (95% CI: -0.72–0.38; 95% PI: -1.70–1.36), and calibration slope was 0.81 (95% CI 0.63–0.99; 95% PI: 0.50–1.12) (Additional file 1: Results of the prediction model without BRCA mutation). Details of development, validation, and clinical utility are provided in Additional file 1: Tables S8–S10 and Figure S5–S10.
In a sensitivity analysis (see Additional file 1: Assessment of limited information of CPM), we studied the impact of CPM on our results using two studies, in which CPM information was (almost) completely available. The lack of CPM information on cumulative incidence estimation hardly affected the results of our analyses (Additional file 2: Figure S11).
Using established risk factors for CBC which are currently available in clinical practice, we developed PredictCBC, which can be used to calculate 5- and 10-year absolute CBC risk. The risk prediction model includes carriership of BRCA1/2 mutations, an important determinant of CBC risk in the decision-making process .
The calibration of the model was reasonable and discrimination moderate within the range of other tools commonly used for routing counseling and decision-making in clinical oncology for primary BC risk [56,57,58,59]. As expected, the prediction accuracy was lower when we omitted the BRCA mutation carrier status although the prevalence of BRCA mutations among BC patients is quite low (2–4%) [60, 61].
In the breast cancer population, CBC is a relatively uncommon event (~ 0.4% per year) and difficult to predict. Therefore, physicians should carefully consider which patients should consider CPM using a prediction model . The current clinical recommendations of CPM are essentially based on the presence of a mutation in the BRCA1/2 genes. Based on the risk distribution defined by the current model (Fig. 3), this is a reasonable approach: essentially no non-carrier women reach a 10% risk 10-year threshold. However, more than 50% of carriers do not reach this threshold either, suggesting that a significant proportion of BRCA1/2 carriers might be spared CPM. Contralateral surveillance mammography may also be avoided although detection and knowledge of recurrences may be necessary for better defined individualized follow-up and patient-tailored treatment strategies [63, 64].
CBC risk patterns and factors were identified previously in a large population-based study with 10,944 CBC of 212,630 patients from the Surveillance, Epidemiology and End Results (SEER) database diagnosed from 1990 to 2013 . However, SEER does not include details of endocrine treatment and chemotherapy, therapies administrated to reduce recurrences and CBCs [13, 66]. Furthermore, in this study, the model was not validated or evaluated based on prediction accuracy, nor was a tool provided. Another study provided general guidelines for CPM by calculating the lifetime risk of CBC based on a published systematic review of age at first BC, BRCA1/2 gene mutation, family history of BC, ER status, ductal carcinoma in situ, and oophorectomy [34, 67]. However, the authors specified that the calculation of the CBC lifetime risk should be considered only as a guide for helping clinicians to stratify patients into risk categories rather than a precise tool for the objective assessment of the risk.
Only one other prediction model (CBCrisk) has been developed and validated using data of 1921 CBC cases and 5763 matched controls . External validation of CBCrisk of two independent datasets using 5185 and 6035 patients with 111 and 117 CBC assessed a discrimination between 0.61 and 0.65 . The discrimination of our PredictCBC model at 5 and 10 years was similar; however, the geographic diversity of the studies gave a more complete overview of external validity . Moreover, we showed the net benefit of our model using decision curve analysis since standard performance metrics of discrimination, calibration, sensitivity, and specificity alone are insufficient to assess the clinical utility [18, 53].
Some limitations of our study must be recognized. First, reporting of CBC was not entirely complete in all studies and information about CPM was limited in most datasets, which may have underestimated the cumulative incidence, although the overall 10-year cumulative incidence of 4.1% is in line with other data [5, 34]. Second, some women included in the Dutch studies (providing specific information on family history, BRCA mutation, or CPM) were also present in our selection of the Netherlands Cancer Registry population. Privacy and coding issues prevented linkage at the individual patient level, but based on the hospitals from which the studies recruited, and the age and period criteria used, we calculated a maximum potential overlap of 3.4%. Third, in the US and Australian datasets, the prediction performance was uncertain due to the limited sample size and missing values. Moreover, some important predictors such as family history and especially BRCA mutation status were only available in a subset of the women (from familial- and unselected hospital-based studies) and patients with data on BRCA mutation status might have been insufficiently represented for tested populations and further development and validation of PredictCBC-1A will be necessary. However, although BRCA1/2 mutation information was unavailable in 94% of our data, the approach of the imputation led to consistently good performing models [68,69,70]. The remaining factors were quite complete: ~ 79% of patients had at most one missing factor, which provided good imputation diagnostic performances. Since most BC patients are not currently tested in the clinical practice for BRCA1/2 mutations, we assessed the clinical utility of PredictCBC version 1B to provide individualized CBC risk estimates for first BC patients not tested for BRCA1/2 germline mutations [60, 71]. Our PredictCBC version 1B model provides less precise estimates, but may be useful in providing general CBC risk estimates, which could steer women away from CPM or trigger BRCA testing.
Last but not the least, adequate presentation of the risk estimates from the PredictCBC-1A and PredictCBC-1B is crucial for effective communication about CBC risk during doctor-patient consultations [72, 73]. A nomogram is an important component to communicate the risk of modern medical decision-making, although it may be difficult to use and might potentially make it more difficult to interpret the risks for laymen  An online tool is being implemented, and a pilot study will be conducted among patients and clinicians to assess how the risk estimates from PredictCBC-1A and 1B can best be visualized to facilitate communication with patients. Other factors, which were not available in our study, predict breast cancer risk and their inclusion may further improve the discrimination and clinical utility of our CBC risk model: these factors include CHEK2 c.1100del mutation carriers, polygenic risk scores based on common genetic variants, breast density, and reproductive and lifestyle factors such as BMI and age at menarche . Additional data with complete information of BRCA1/2 mutation should be also considered in the model upgrade to reduce uncertainty of CBC risk estimates. External validation in other studies, including patients of other ethnicities, will also be important. In the meantime, our model provides a reliable basis for CBC risk counseling.
In conclusion, we have developed and cross-validated risk prediction models for CBC (PredictCBC) based on different European-descent population and hospital-based studies. The model is reasonably calibrated and prediction accuracy is moderate. The clinical utility assessment of PredictCBC showed potential for improved risk counseling, although the decision regarding CPM in the general breast cancer population remains challenging. Similar results have been found for PredictCBC version 1B, a CBC risk prediction model that calculates individualized CBC risk for first BC patients not tested for BRCA1/2 germline mutation.
Availability of data and materials
All data relevant to this report are included in this published article and its supplementary information files. The datasets analyzed during the current study are not publicly available due to protection of participant privacy and confidentiality, and ownership of the contributing institutions, but may be made available in an anonymized form via the corresponding author on reasonable request and after approval of the involved institutions.
Area under the ROC curve
Breast Cancer Association Consortium
Body mass index
Contralateral breast cancer
Contralateral preventive mastectomy
Decision curve analysis
Human epidermal growth receptor 2
Multiple imputation by chained equations
Surveillance, Epidemiology and End Results
TNM Classification of Malignant Tumors
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Survival and prevalence of cancer. https://www.cijfersoverkanker.nl. Accessed Mar 2019.
Schaapveld M, Visser O, Louwman WJ, Willemse PH, de Vries EG, van der Graaf WT, Otter R, Coebergh JW, van Leeuwen FE. The impact of adjuvant therapy on contralateral breast cancer risk and the prognostic significance of contralateral breast cancer: a population based study in the Netherlands. Breast Cancer Res Treat. 2008;110(1):189–97.
Brenner DJ. Contralateral second breast cancers: prediction and prevention. J Natl Cancer Inst. 2010;102(7):444–5.
van den Broek AJ, van ‘t Veer LJ, Hooning MJ, Cornelissen S, Broeks A, Rutgers EJ, Smit VT, Cornelisse CJ, van Beek M, Janssen-Heijnen ML, et al. Impact of age at primary breast cancer on contralateral breast cancer risk in BRCA1/2 mutation carriers. J Clin Oncol. 2016;34(5):409–18.
Malone KE, Begg CB, Haile RW, Borg A, Concannon P, Tellhed L, Xue S, Teraoka S, Bernstein L, Capanu M, et al. Population-based study of the risk of second primary contralateral breast cancer associated with carrying a mutation in BRCA1 or BRCA2. J Clin Oncol. 2010;28(14):2404–10.
Evans DG, Ingham SL, Baildam A, Ross GL, Lalloo F, Buchan I, Howell A. Contralateral mastectomy improves survival in women with BRCA1/2-associated breast cancer. Breast Cancer Res Treat. 2013;140(1):135–42.
Graeser MK, Engel C, Rhiem K, Gadzicki D, Bick U, Kast K, Froster UG, Schlehe B, Bechtold A, Arnold N, et al. Contralateral breast cancer risk in BRCA1 and BRCA2 mutation carriers. J Clin Oncol. 2009;27(35):5887–92.
Weischer M, Nordestgaard BG, Pharoah P, Bolla MK, Nevanlinna H, Van't Veer LJ, Garcia-Closas M, Hopper JL, Hall P, Andrulis IL, et al. CHEK2*1100delC heterozygosity in women with breast cancer associated with early death, breast cancer-specific death, and increased risk of a second breast cancer. J Clin Oncol. 2012;30(35):4308–16.
Kuchenbaecker KB, Hopper JL, Barnes DR, Phillips KA, Mooij TM, Roos-Blom MJ, Jervis S, van Leeuwen FE, Milne RL, Andrieu N, et al. Risks of breast, ovarian, and contralateral breast cancer for BRCA1 and BRCA2 mutation carriers. JAMA. 2017;317(23):2402–16.
Domchek SM. Risk-reducing mastectomy in BRCA1 and BRCA2 mutation carriers: a complex discussion. JAMA. 2019;321(1):27.
Chen Y, Thompson W, Semenciw R, Mao Y. Epidemiology of contralateral breast cancer. Cancer Epidemiol Biomark Prev. 1999;8(10):855–61.
Kramer I, Schaapveld M, Oldenburg HSA, Sonke GS, McCool D, van Leeuwen FE, Van de Vijver KK, Russell NS, Linn SC, Siesling S, et al. The influence of adjuvant systemic regimens on contralateral breast cancer risk and receptor subtype. J Natl Cancer Inst. 2019;111(7):709–18.
Portschy PR, Abbott AM, Burke EE, Nzara R, Marmor S, Kuntz KM, Tuttle TM. Perceptions of contralateral breast cancer risk: a prospective, longitudinal study. Ann Surg Oncol. 2015;22(12):3846–52.
Murphy JA, Milner TD, O'Donoghue JM. Contralateral risk-reducing mastectomy in sporadic breast cancer. Lancet Oncol. 2013;14(7):e262–9.
Chowdhury M, Euhus D, Onega T, Biswas S, Choudhary PK. A model for individualized risk prediction of contralateral breast cancer. Breast Cancer Res Treat. 2017;161(1):153–60.
Chowdhury M, Euhus D, Arun B, Umbricht C, Biswas S, Choudhary P. Validation of a personalized risk prediction model for contralateral breast cancer. Breast Cancer Res Treat. 2018;170:415. https://link.springer.com/article/10.1007/s10549-018-4763-5#citeas.
Vickers AJ, Van Calster B, Steyerberg EW. Net benefit approaches to the evaluation of prediction models, molecular markers, and diagnostic tests. BMJ. 2016;352:i6.
Michailidou K, Lindstrom S, Dennis J, Beesley J, Hui S, Kar S, Lemacon A, Soucy P, Glubb D, Rostamianfar A, et al. Association analysis identifies 65 new breast cancer risk loci. Nature. 2017;551(7678):92–4.
Schmidt MK, Tollenaar RA, de Kemp SR, Broeks A, Cornelisse CJ, Smit VT, Peterse JL, van Leeuwen FE, Van't Veer LJ. Breast cancer survival and tumor characteristics in premenopausal women carrying the CHEK2*1100delC germline mutation. J Clin Oncol. 2007;25(1):64–9.
Schmidt MK, van den Broek AJ, Tollenaar RA, Smit VT, Westenend PJ, Brinkhuis M, Oosterhuis WJ, Wesseling J, Janssen-Heijnen ML, Jobsen JJ, et al. Breast cancer survival of BRCA1/BRCA2 mutation carriers in a hospital-based cohort of young women. J Natl Cancer Inst. 2017;109(8). https://academic.oup.com/jnci/article/109/8/djw329/3064570.
Font-Gonzalez A, Liu L, Voogd AC, Schmidt MK, Roukema JA, Coebergh JW, de Vries E, Soerjomataram I. Inferior survival for young patients with contralateral compared to unilateral breast cancer: a nationwide population-based study in the Netherlands. Breast Cancer Res Treat. 2013;139(3):811–9.
Riegman PH, van Veen EB. Biobanking residual tissues. Hum Genet. 2011;130(3):357–68.
Foundation Federation of Dutch Medical Scientific Societies: Human tissue and medical research: code of conduct for responsible use. 2011.
Vichapat V, Garmo H, Holmqvist M, Liljegren G, Warnberg F, Lambe M, Fornander T, Adolfsson J, Luchtenborg M, Holmberg L. Tumor stage affects risk and prognosis of contralateral breast cancer: results from a large Swedish-population-based study. J Clin Oncol. 2012;30(28):3478–85.
Vichapat V, Gillett C, Fentiman IS, Tutt A, Holmberg L, Luchtenborg M. Risk factors for metachronous contralateral breast cancer suggest two aetiological pathways. Eur J Cancer. 2011;47(13):1919–27.
Mariani L, Coradini D, Biganzoli E, Boracchi P, Marubini E, Pilotti S, Salvadori B, Silvestrini R, Veronesi U, Zucali R, et al. Prognostic factors for metachronous contralateral breast cancer: a comparison of the linear Cox regression model and its artificial neural network extension. Breast Cancer Res Treat. 1997;44(2):167–78.
Reiner AS, Lynch CF, Sisti JS, John EM, Brooks JD, Bernstein L, Knight JA, Hsu L, Concannon P, Mellemkjaer L, et al. Hormone receptor status of a first primary breast cancer predicts contralateral breast cancer risk in the WECARE study population. Breast Cancer Res. 2017;19(1):83.
Sisti JS, Bernstein JL, Lynch CF, Reiner AS, Mellemkjaer L, Brooks JD, Knight JA, Bernstein L, Malone KE, Woods M, et al. Reproductive factors, tumor estrogen receptor status and contralateral breast cancer risk: results from the WECARE study. Springerplus. 2015;4:825.
Healey EA, Cook EF, Orav EJ, Schnitt SJ, Connolly JL, Harris JR. Contralateral breast cancer: clinical characteristics and impact on prognosis. J Clin Oncol. 1993;11(8):1545–52.
Gao X, Fisher SG, Emami B. Risk of second primary cancer in the contralateral breast in women treated for early-stage breast cancer: a population-based study. Int J Radiat Oncol Biol Phys. 2003;56(4):1038–45.
Brooks JD, John EM, Mellemkjaer L, Lynch CF, Knight JA, Malone KE, Reiner AS, Bernstein L, Liang X, Shore RE, et al. Body mass index, weight change, and risk of second primary breast cancer in the WECARE study: influence of estrogen receptor status of the first breast cancer. Cancer Med. 2016;5(11):3282–91.
Knight JA, Blackmore KM, Fan J, Malone KE, John EM, Lynch CF, Vachon CM, Bernstein L, Brooks JD, Reiner AS, et al. The association of mammographic density with risk of contralateral breast cancer and change in density with treatment in the WECARE study. Breast Cancer Res. 2018;20(1):23.
Basu NN, Barr L, Ross GL, Evans DG. Contralateral risk-reducing mastectomy: review of risk factors and risk-reducing strategies. Int J Surg Oncol. 2015;2015:901046.
Akdeniz D, Schmidt MK, Seynaeve CM, McCool D, Giardiello D, van den Broek AJ, Hauptmann M, Steyerberg EW, Hooning MJ. Risk factors for metachronous contralateral breast cancer: a systematic review and meta-analysis. Breast. 2018;44:1–14.
Edge SB, Compton CC. The American Joint Committee on Cancer: the 7th edition of the AJCC cancer staging manual and the future of TNM. Ann Surg Oncol. 2010;17(6):1471–4.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2017. http://www.R-project.org/, https://cran.r-project.org/doc/FAQ/R-FAQ.html#Citing-R.
van den Broek AJ, Schmidt MK, van ‘t Veer LJ, HSA O, Rutgers EJ, Russell NS, Smit V, Voogd AC, Koppert LB, Siesling S, et al. Prognostic impact of breast-conserving therapy versus mastectomy of BRCA1/2 mutation carriers compared with noncarriers in a consecutive series of young breast cancer patients. Ann Surg. 2019;270(2):364–72.
Resche-Rigon M, White IR, Bartlett JW, Peters SA, Thompson SG, Group P-IS. Multiple imputation for handling systematically missing confounders in meta-analysis of individual participant data. Stat Med. 2013;32(28):4890–905.
Sv B. Flexible imputation of missing data. Boca Raton: CRC Press; 2012.
Geskus RB. Cause-specific cumulative incidence estimation and the fine and gray model under both left truncation and right censoring. Biometrics. 2011;67(1):39–49.
Fine JP, Gray RJ. A proportional hazards model for the subdistribution of a competing risk. J Am Stat Assoc. 1999;94(446):496–509.
Schoenfeld DA. Sample-size formula for the proportional-hazards regression model. Biometrics. 1983;39(2):499–503.
Little RJA, Rubin DB. Statistical analysis with missing data. New York: Wiley; 1987.
Zhang Z, Geskus RB, Kattan MW, Zhang H, Liu T. Nomogram for survival analysis in the presence of competing risks. Ann Transl Med. 2017;5(20):403.
Steyerberg EW, Harrell FE Jr. Prediction models need appropriate internal, internal-external, and external validation. J Clin Epidemiol. 2016;69:245–7.
Austin PC, van Klaveren D, Vergouwe Y, Nieboer D, Lee DS, Steyerberg EW. Geographic and temporal validity of prediction models: different approaches were useful to examine model performance. J Clin Epidemiol. 2016;79:76–85.
Collins GS, Ogundimu EO, Altman DG. Sample size considerations for the external validation of a multivariable prognostic model: a resampling study. Stat Med. 2016;35(2):214–26.
Steyerberg EW. Clinical prediction models: a practical approach to development, validation and updating. New York: Springer; 2010.
Blanche P, Dartigues JF, Jacqmin-Gadda H. Estimating and comparing time-dependent areas under receiver operating characteristic curves for censored event times with competing risks. Stat Med. 2013;32(30):5381–97.
Snell KI, Hua H, Debray TP, Ensor J, Look MP, Moons KG, Riley RD. Multivariate meta-analysis of individual participant data helped externally validate the performance and implementation of a prediction model. J Clin Epidemiol. 2016;69:40–50.
Vickers AJ, Elkin EB. Decision curve analysis: a novel method for evaluating prediction models. Med Decis Mak. 2006;26(6):565–74.
Kerr KF, Brown MD, Zhu K, Janes H. Assessing the clinical impact of risk prediction models with decision curves: guidance for correct interpretation and appropriate use. J Clin Oncol. 2016;34(21):2534–40.
Vickers AJ, Cronin AM, Elkin EB, Gonen M. Extensions to decision curve analysis, a novel method for evaluating diagnostic tests, prediction models and molecular markers. BMC Med Inform Decis Mak. 2008;8:53.
Heemskerk-Gerritsen BA, Rookus MA, Aalfs CM, Ausems MG, Collee JM, Jansen L, Kets CM, Keymeulen KB, Koppert LB, Meijers-Heijboer HE, et al. Improved overall survival after contralateral risk-reducing mastectomy in BRCA1/2 mutation carriers with a history of unilateral breast cancer: a prospective analysis. Int J Cancer. 2015;136(3):668–77.
Rockhill B, Spiegelman D, Byrne C, Hunter DJ, Colditz GA. Validation of the Gail et al. model of breast cancer risk prediction and implications for chemoprevention. J Natl Cancer Inst. 2001;93(5):358–66.
Elmore JG, Fletcher SW. The risk of cancer risk prediction: “what is my risk of getting breast cancer”? J Natl Cancer Inst. 2006;98(23):1673–5.
Wishart GC, Azzato EM, Greenberg DC, Rashbass J, Kearins O, Lawrence G, Caldas C, Pharoah PD. PREDICT: a new UK prognostic model that predicts survival following surgery for invasive breast cancer. Breast Cancer Res. 2010;12(1):R1.
Goldstein LJ, Gray R, Badve S, Childs BH, Yoshizawa C, Rowley S, Shak S, Baehner FL, Ravdin PM, Davidson NE, et al. Prognostic utility of the 21-gene assay in hormone receptor-positive operable breast cancer compared with classical clinicopathologic features. J Clin Oncol. 2008;26(25):4063–71.
van den Broek AJ, de Ruiter K, van ‘t Veer LJ, Tollenaar RA, van Leeuwen FE, Verhoef S, Schmidt MK. Evaluation of the Dutch BRCA1/2 clinical genetic center referral criteria in an unselected early breast cancer population. Eur J Hum Genet. 2015;23(5):588–95.
Gail MH, Pfeiffer RM. Breast cancer risk model requirements for counseling, prevention, and screening. J Natl Cancer Inst. 2018;110(9):994–1002.
O'Donnell M. Estimating contralateral breast cancer risk. Current Breast Cancer Rep. 2018;10(2):91–7.
van Maaren MC, de Munck L, Strobbe LJA, Sonke GS, Westenend PJ, Smidt ML, Poortmans PMP, Siesling S. Ten-year recurrence rates for breast cancer subtypes in the Netherlands: a large population-based study. Int J Cancer. 2019;144(2):263–72.
Lu W, Schaapveld M, Jansen L, Bagherzadegan E, Sahinovic MM, Baas PC, Hanssen LM, van der Mijle HC, Brandenburg JD, Wiggers T, et al. The value of surveillance mammography of the contralateral breast in patients with a history of breast cancer. Eur J Cancer. 2009;45(17):3000–7.
Xiong Z, Yang L, Deng G, Huang X, Li X, Xie X, Wang J, Shuang Z, Wang X. Patterns of occurrence and outcomes of contralateral breast cancer: analysis of SEER data. J Clin Med. 2018;7(6):133. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6025574/.
Langballe R, Mellemkjaer L, Malone KE, Lynch CF, John EM, Knight JA, Bernstein L, Brooks J, Andersson M, Reiner AS, et al. Systemic therapy for breast cancer and risk of subsequent contralateral breast cancer in the WECARE study. Breast Cancer Res. 2016;18(1):65.
Basu NN, Ross GL, Evans DG, Barr L. The Manchester guidelines for contralateral risk-reducing mastectomy. World J Surg Oncol. 2015;13:237.
Nieboer D, Vergouwe Y, Ankerst DP, Roobol MJ, Steyerberg EW. Improving prediction models with new markers: a comparison of updating strategies. BMC Med Res Methodol. 2016;16(1):128.
Collins GS, Altman DG. Predicting the 10 year risk of cardiovascular disease in the United Kingdom: independent and external validation of an updated version of QRISK2. BMJ. 2012;344:e4181.
Madley-Dowd P, Hughes R, Tilling K, Heron J. The proportion of missing data should not be used to guide decisions on multiple imputation. J Clin Epidemiol. 2019;110:63–73.
Childers CP, Childers KK, Maggard-Gibbons M, Macinko J. National estimates of genetic testing in women with a history of breast or ovarian cancer. J Clin Oncol. 2017;35(34):3800–6.
Bonnett LJ, Snell KIE, Collins GS, Riley RD. Guide to presenting clinical prediction models for use in clinical settings. BMJ. 2019;365:l737.
Van Belle V, Van Calster B. Visualizing risk prediction models. PLoS One. 2015;10(7):e0132614.
Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol. 2015;16(4):e173–80.
Mavaddat N, Michailidou K, Dennis J, Lush M, Fachal L, Lee A, Tyrer JP, Chen TH, Wang Q, Bolla MK, et al. Polygenic risk scores for prediction of breast cancer and breast cancer subtypes. Am J Hum Genet. 2019;104(1):21–34.
We thank all individuals who took part in these studies and all researchers, clinicians, technicians, and administrative staff who have enabled this work to be carried out.
ABCFS thank Maggie Angelakos, Judi Maskiell, Gillian Dite. ABCS and BOSOM thanks all the collaborating hospitals and pathology departments and many individual that made this study possible; specifically, we wish to acknowledge Annegien Broeks, Sten Cornelissen, Frans Hogervorst, Laura van ‘t Veer, Floor van Leeuwen, Emiel Rutgers. EMC thanks J.C. Blom-Leenheer, P.J. Bos, C.M.G. Crepin, and M. van Vliet for the data management. CGPS thanks staff and participants of the Copenhagen General Population Study. For the excellent technical assistance: Dorthe Uldall Andersen, Maria Birna Arnadottir, Anne Bank, Dorthe Kjeldgård Hansen. HEBCS thanks Taru A. Muranen, Kristiina Aittomäki, Karl von Smitten, and Irja Erkkilä. KARMA thanks the Swedish Medical Research Counsel. LMBC thanks Gilian Peuteman, Thomas Van Brussel, EvyVanderheyden, and Kathleen Corthouts. MARIE thanks Petra Seibold, Dieter Flesch-Janys, Judith Heinz, Nadia Obi, Alina Vrieling, Sabine Behrens, Ursula Eilber, Muhabbet Celik, Til Olchers, and Stefan Nickels. ORIGO thanks E. Krol-Warmerdam, and J. Blom for patient accrual, administering questionnaires, and managing clinical information. The authors thank the registration team of the Netherlands Comprehensive Cancer Organisation (IKNL) for the collection of data for the Netherlands Cancer Registry as well as IKNL staff for scientific advice. PBCS thanks Louise Brinton, Mark Sherman, Neonila Szeszenia-Dabrowska, Beata Peplonska, Witold Zatonski, Pei Chao, and Michael Stagner. The ethical approval for the POSH study is MREC /00/6/69, UKCRN ID: 1137. We thank the SEARCH team.
This work is supported by the Alpe d’HuZes/Dutch Cancer Society (KWF Kankerbestrijding) project 6253.
BCAC is funded by Cancer Research UK [C1287/A16563, C1287/A10118], the European Union’s Horizon 2020 Research and Innovation Programme (grant numbers 634935 and 633784 for BRIDGES and B-CAST respectively), and the European Community’s Seventh Framework Programme under grant agreement number 223175 (grant number HEALTH-F2-2009-223175) (COGS). The EU Horizon 2020 Research and Innovation Programme funding source had no role in the study design, data collection, data analysis, data interpretation, or writing of the report.
The Australian Breast Cancer Family Study (ABCFS) was supported by grant UM1 CA164920 from the National Cancer Institute (USA). The ABCFS was also supported by the National Health and Medical Research Council of Australia, the New South Wales Cancer Council, the Victorian Health Promotion Foundation (Australia), and the Victorian Breast Cancer Research Consortium. J.L.H. is a National Health and Medical Research Council (NHMRC) Senior Principal Research Fellow. M.C.S. is a NHMRC Senior Research Fellow. The ABCS study was supported by the Dutch Cancer Society [grants NKI 2007-3839; 2009 4363]. The work of the BBCC was partly funded by ELAN-Fond of the University Hospital of Erlangen. BOSOM was supported by the Dutch Cancer Society grant numbers DCS-NKI 2001-2423, DCS-NKI 2007-3839, and DCSNKI 2009-4363; the Cancer Genomics Initiative; and notary office Spier & Hazenberg for the coding procedure. The EMC was supported by grants from Alpe d’HuZes/Dutch Cancer Society NKI2013-6253 and from Pink Ribbon 2012.WO39.C143. The HEBCS was financially supported by the Helsinki University Hospital Research Fund, the Finnish Cancer Society, and the Sigrid Juselius Foundation.
Financial support for KARBAC was provided through the regional agreement on medical training and clinical research (ALF) between Stockholm County Council and Karolinska Institutet, the Swedish Cancer Society, The Gustav V Jubilee foundation and Bert von Kantzows foundation. The KARMA study was supported by Märit and Hans Rausings Initiative Against Breast Cancer. LMBC is supported by the ‘Stichting tegen Kanker’. The MARIE study was supported by the Deutsche Krebshilfe e.V. [70-2892-BR I, 106332, 108253, 108419, 110826, 110828], the Hamburg Cancer Society, the German Cancer Research Center (DKFZ) and the Federal Ministry of Education and Research (BMBF) Germany [01KH0402]. MEC was supported by NIH grants CA63464, CA54281, CA098758, CA132839, and CA164973. The ORIGO study was supported by the Dutch Cancer Society (RUL 1997-1505) and the Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL CP16). The PBCS was funded by Intramural Research Funds of the National Cancer Institute, Department of Health and Human Services, USA. Genotyping for PLCO was supported by the Intramural Research Program of the National Institutes of Health, NCI, Division of Cancer Epidemiology and Genetics. The POSH study is funded by Cancer Research UK (grants C1275/A11699, C1275/C22524, C1275/A19187, C1275/A15956) and Breast Cancer Campaign 2010PR62, 2013PR044. PROCAS is funded from NIHR grant PGfAR 0707-10031. SEARCH is funded by Cancer Research UK [C490/A10124, C490/A16561] and supported by the UK National Institute for Health Research Biomedical Research Centre at the University of Cambridge. SKKDKFZS is supported by the DKFZ. The SZBCS (Szczecin Breast Cancer Study) was supported by Grant PBZ_KBN_122/P05/2004 and The National Centre for Research and Development (NCBR) within the framework of the international ERA-NET TRANSAN JTC 2012 application no. Cancer 12-054 (Contract No. ERA-NET-TRANSCAN / 07/2014).
Ethics approval and consent to participate
Each study was approved by its institutional ethical review board.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1 Table S1. Data source flowchart. Table S2. Description of the studies included in the analyses. Table S3. Patients and first primary breast cancer characteristics used in the contralateral breast cancer risk prediction model in the complete case and all case analyses. Table S4. Results of multivariable subdistributional hazard model using the complete case dataset. Table S5. List of BCAC studies (including ABCS source) with the corresponding country and geographic area. Table S6. Main patient and disease characteristics. Table S7. Clinical utility of the 5-year contralateral breast cancer risk prediction model. Table S8. Results of multivariable subdistributional hazard model for breast cancer patients without BRCA mutations. Table S9. Clinical utility of the 5-year contralateral breast cancer risk prediction model in non-BRCA tested patients. Table S10. Clinical utility of the 10-year contralateral breast cancer risk prediction model in non-BRCA tested patients.
Additional file 2 Figure S1. Graphical assessment of non-linear relationship of age with contralateral breast cancer risk. Figure S2. Visual assessment of calibration through calibration plots in the internal-external cross-validation at 5 years for the contralateral breast cancer risk model with BRCA mutation information. Figure S3. Visual assessment of calibration through calibration plots in the internal-external cross-validation at 10 years for the contralateral breast cancer risk model with BRCA mutation information. Figure S4. Decision curve analysis at 5 years for the contralateral breast cancer risk model including BRCA1/2 mutation information. Figure S5. Results of the leave-one-study-out cross-validation for the contralateral breast cancer risk model at 5 and 10 years without BRCA mutation information. Figure S6. Visual assessment of calibration through calibration plots in the internal-external cross-validation at 5 years for the contralateral breast cancer risk model without BRCA mutation information. Figure S7. Visual assessment of calibration through calibration plots in the internal-external cross-validation at 10 years for the contralateral breast cancer risk model without BRCA mutation information. Figure S8. Density distribution of 10-year predicted absolute risk in patients with no family history and patients with a family history. Figure S9. Decision curve analysis at 5 years for the contralateral breast cancer risk model without BRCA mutation information. Figure S10. Decision curve analysis at 10 years for the contralateral breast cancer risk model without BRCA mutation information. Figure S11. Assessment of inclusion of information of contralateral preventive mastectomy (CPM).
About this article
Cite this article
Giardiello, D., Steyerberg, E.W., Hauptmann, M. et al. Prediction and clinical utility of a contralateral breast cancer risk model. Breast Cancer Res 21, 144 (2019). https://doi.org/10.1186/s13058-019-1221-1
- Contralateral breast cancer
- Risk prediction model
- Clinical decision-making
- BRCA mutation carriers