Lymph node positivity in different early breast carcinoma phenotypes: a predictive model
BMC Cancer volume 19, Article number: 45 (2019)
A strong correlation between breast cancer (BC) molecular subtypes and axillary status has been shown. It would be useful to predict the probability of lymph node (LN) positivity. Objective: To develop the performance of multivariable models to predict LN metastases, including nomograms derived from logistic regression with clinical, pathologic variables provided by tumor surgical results or only by biopsy.
A retrospective cohort was randomly divided into two separate patient sets: a training set and a validation set. In the training set, we used multivariable logistic regression techniques to build different predictive nomograms for the risk of developing LN metastases. The discrimination ability and calibration accuracy of the resulting nomograms were evaluated on the training and validation set.
Consecutive sample of 12,572 early BC patients with sentinel node biopsies and no neoadjuvant therapy. In our predictive macro metastases LN model, the areas under curve (AUC) values were 0.780 and 0.717 respectively for pathologic and pre-operative model, with a good calibration, and results with validation data set were similar: AUC respectively of 0.796 and 0.725.
Among the list of candidate’s regression variables, on the training set we identified age, tumor size, LVI, and molecular subtype as statistically significant factors for predicting the risk of LN metastases.
Several nomograms were reported to predict risk of SLN involvement and NSN involvement. We propose a new calculation model to assess this risk of positive LN with similar performance which could be useful to choose management strategies, to avoid axillary LN staging or to propose ALND for patients with high level probability of major axillary LN involvement but also to propose immediate breast reconstruction when post mastectomy radiotherapy is not required for patients without LN macro metastasis.
A retrospective cohort of 12.572 early BC patients with SN biopsies was randomly divided into two separate patient sets to develop, validate and compare different predictive nomograms for the risk of developing LN metastases from clinical and pathologic variables provided by tumor surgical results or by biopsy.
In breast cancer (BC), nodal status is a major prognostic factor that determines therapeutic decisions to a large extent. Sentinel lymph node biopsy (SLNB) provides a reliable assessment of the axilla status in early clinically node-negative BC . Since it also causes less morbidity than axillary lymph node dissection (ALND), it is now considered as a standard of care procedure. The omission of completion ALND in patients with negative sentinel lymph nodes (SLN) has been recognized as a reasonable attitude since the publication of the NSABP B-32 results . Moreover, it is likely that it can be safely expanded to patients with minimal SLN involvement (isolated tumor cells and micro metastases), with regard to survival outcomes [3, 4]. Indeed, 40 to 70% of these patients do not have metastatic non-sentinel lymph nodes (NSLN) . Main predictors of LN metastases are tumor size, grade, lymphovascular invasion (LVI), age at diagnosis, extracapsular extension of the positive SLN, and hormonal and HER2 receptor status [6,7,8,9,10]. In addition, a strong correlation between BC molecular subtypes and /or tumor phenotypes on the one hand (determined by hormonal receptor and HER2 status) and axillary status on the other hand has been shown in numerous studies [11,12,13,14,15,16].
The determination of the risk of positive axillary LN can significantly contribute to therapeutic decisions. However, this risk cannot be immediately induced from the results of multivariate analyses that provide broad statistical information. Only an appropriate prediction tool, using a nomogram, can indicate the individual risk of a given patient. These nomograms can also be used to compare populations from different studies. A large cohort is necessary to reliably determine the probability of positive SN, particularly for less frequent tumor phenotypes. Reyal et al. published such a nomogram predictive of the risk of developing SN metastases in 2011 , built on a training set made of 1543 early-stage BC patients, and validated on two cohorts of 615 and 496 patients respectively. This model was further validated in a cohort of 755 consecutive patients treated at Institut Curie in 2009 .
The aim of our study was to develop and compare the performance of multivariable models to predict LN metastases, including nomograms derived from logistic regression with clinical, pathologic variables provided by tumor surgical results or only provided by biopsy as explanatory variables.
Our cohort consisted of 12,572 consecutive patients with small (≤ 30 mm based on clinical and radiologic findings), clinically node-negative invasive BC, who did not receive neoadjuvant therapy, and underwent SLNB between 1999 and 2012 at 13 French centers. HER2 status was determined for all patients. During the first years of the study, ALND was systematically performed in some sites; thereafter, ALND was performed only in case of SN involvement, this attitude being homogeneous within all the participating sites.
The following data were retrieved: characteristics of patients (age at the time of SLNB), and tumors [size, clinical stage, histological type, estrogen (ER), progesterone (PR) and HER2 status, LVI, Scarff-Bloom-Richardson (SBR) grade], description of ALND (number of LN sampled and involved), and results of the pathological examination of surgical resection specimens. Tumor size was determined on the results of pathological examination but could be evaluated pre operatively by mammography, sonography and in selected cases by MRI (clinical T stage). LVI was detected on surgical specimen.
Tumor phenotype was defined by the combination of ER, PR and HER2 status, evaluated by immuno-histochemistry (IHC) and confirmed by FISH in case of IHC-HER2 2+. Positivity for ER and PR was determined according to French guidelines (≥ 10% of cancer cells expressing ER/PR). Five molecular subtypes were defined according to clinico-pathological criteria . Because information on Ki-67 was not available, we used grade to capture cell proliferation, as described by von Minckwitz et al  The following definitions were used: triple-negative (basal-like, HER2-/HR-), HER2 positive (non-luminal, HER2+/HR-), and luminal (HR+), divided into luminal A (HR+/HER2−/grade1 or 2), luminal B-HER2-negative like (HR+/HER2−/grade 3), and luminal B-HER2-positive like (HR+/HER2+ all grades).
Although the methods used for histological examination were not standardized in the protocol, all sites proceeded similarly: serial sections were performed every 200 μm and stained with standard hematoxylin and eosin. The number of sections was six to ten, or pursued until node exhaustion in case of large SN. Additional IHC analysis was done in case of negative results at standard examination. For additional nodes identified by completion ALND, routine HE analysis was performed.
Five categories of LN status were defined: negative LN (pN0i-), isolated tumor cells (pN0(i+): < 0.2 mm), detected either by hematoxylin and eosin (HE) staining or by cytokeratin IHC, micro metastases (pN1mi: > 0.2 mm and < 2 mm), and macro metastases (> 2 mm), divided into single and multiple macro metastases .
Our main objective was to create prediction models for the risk of LN positivity and the risk of LN macroscopic metastases from clinical and pathologic variables provided by tumor surgical results or by biopsy, and evaluate their performance with respect to three main features: discrimination (i.e. whether the relative ranking of individual predictions is in the correct order), calibration (i.e. agreement between observed outcomes and predictions) and clinical utility defined as proportions of patients classified into risk categories using predefined cutoff values (< 10%, between 10 and 20%, between 20 and 30%, between 30 and 40%, and > = 40%). Our main evaluation criteria were based on the final status of LN metastases (pN0(i+), pN1mi or pN1ma) as the result of SLNB alone or the final result of both SLNB and ALND. LN positivity was defined as the presence of isolated tumor cells, micro or macro LN metastases. We used logistic regression models  including age (<=40, 41–75,> 75), tumor size (<=20, 20–30, > = 30 mm) or clinical T stage (T0-T1, T2, T3-T4), tumor grade, histology type, LVI, and molecular subtypes as predictor factors to predict each individual risks. The list of predictor factors was set beforehand, based on the investigator’s experience and some reference papers [6,7,8,9,10,11, 13,14,15, 17]. No additional procedure was used in regression analysis to reduce the list of only 5 or 6 predictor factors identified beforehand. Prior to analysis, we randomly divided our initial cohort (N = 12,572) in two separate sub-cohorts: a large training cohort (N = 8381) to create prediction models and a confirmatory cohort (N = 4191) to evaluate their individual’s prediction performance. A split-sample approach was adopted in order to estimate unbiasedly the model performance, as these estimates are known to be biased upwards when regression parameters are estimated on the same dataset . First we performed a descriptive analysis using the following criteria: patient’s age at SN biopsy, clinical and pathological tumor size, tumor grade and histology type, lymphovascular invasion or not (LVI), presence of estrogen (ER), progesterone (PR) and hormonal receptors (RH), Her2 positivity, tumor subtype, number of SN removed and final LN status. The evaluation of each model was assessed in the training sample and the confirmatory sample. Differences in patient’s and tumor’s characteristics were compared using Chi Square or exact Fisher test, Student or Wilcoxon rank sum tests as appropriate. The discrimination ability was evaluated by the area under the ROC (Receiver Operating Characteristic) curve (AUC). We used the functions roc and pROC implemented in R to estimate AUC with 95%CI and test for difference in AUCs along the Delong’s method in the confirmatory sample . Empirical distributions of AUC observed after re-fitting a model on bootstrap replicates (B = 2000) were used to estimate AUC and difference in AUCs with 95% Ci in the training sample. Model calibration was evaluated using Hosmer goodness-of-fit test . All statistical analyses were conducted in the R Language and Environment for Statistical Computing version 3.2.5 (The R Foundation, Vienna, Austria).
Patients’ main characteristics are summarized in Table 1. SBR grade was 1, 2 and 3 in 34, 46 and 20% of cases respectively. Hormone receptor-positive tumors (ER+ and/or PR+) accounted for 88% of cases (11,013 patients). Final LN status, taking into account ALND results when performed, was: pN0(i-) in 8253 patients (66%), pN0(i+) in 355 (3%), pN1mi in 970 (8%) and macro metastasis in 2994 (24%). The comparison between patients with positive and negative final LN status, and between patients with LN macro metastases versus pN0 or pNo(i+) or pN1mi showed statistically significant differences with regard to age, pathologic tumor size, SBR grade, LVI, histological type and distribution of molecular subtypes (Tables 2 and 3).
We first predicted the individual probabilities of final LN positivity and of detecting LN macro metastases from selected clinico-pathologic predictor factors provided by tumor surgical results. The model AUCs with 95% CIs for confirmatory and training samples were respectively 0.767 [0.750–0.783] and 0.755 [0.744–0.767]. Calibration plot and Hosmer-Lemeshow test revealed that the calibration is adequate (p = 0.332 in confirmatory sample, p = 0.158 in training sample). With respect to clinical utility in confirmatory and training samples, the probability of positive LN were respectively below 10% for 7 patients (< 1%) and 19 patients (< 1%), between 10 and 20% for 1096 (31%) and 2255 (32%), and ≥ 20% for 2409 (68.6%) and 4859 (68.1%) patients (Table 4) (Fig. 1A, Additional file 1: Figure S1A and Additional file 2: Figure S2A). The second pathological model estimated the probability of detecting LN macro metastases only. The AUC values for confirmatory and training samples were respectively 0.798 (0.780–0.815) and 0.780 [0.767–0.790]. Clinical utility measures, estimated the probability of LN macro metastases respectively in confirmatory and training samples below 10% for 1004 patients (29%) and 2029 patients (28%), between 10 and 20% for 1075 patients (31%) and 2289 patients (32%), and > 20% for 1433 patients (41%) and 2815 patients (39.4%). The Hosmer-Lemeshow test revealed a poor calibration of the model (p = 0.024 in confirmatory sample, p = 0.427 in training sample) (Table 4 and Additional file 1: Table S1) (Fig. 1B, Additional file 2: Figure S1B, Additional file 3: Figure S2B).
We evaluated the loss in discrimination ability in pre-operative prediction models omitting the information about LVI and substituting pathological tumor size information by clinical T stage. For the overall probability of LN positivity, the AUC values for confirmatory and training samples were respectively 0.687 [0.669–0.705] and 0.682 (0.669–0.694). For the probability of detecting LN macro metastases, the observed AUC results for confirmatory and training samples were respectively 0.727 [0.707–0.746] and 0.717 (0.703–0.732). The calibration of both pre-operative models was found satisfactory. (Table 4) (Fig. 1 C-D, Additional file 2: Figure S1C-D, Additional file 3: Figure S2C-D).
The change in AUCs between pathological and per-operative model were found statistically significantly decreased (p < 0.001). We also evaluated in the confirmatory sample the discrimination ability of the prediction models obtained when treating the variable age and tumor size as continuous. The AUC values for predicting LN positivity and the presence of LN metastases were respectively 0.774 [0.758, 0.79] and 0.805 [0.789–0.823]. The observed increases were significantly (p = 0.041 and p = 0.026), but the results in terms of calibration were judged inadequate (Hosmer-Lemeshow p value < 0.001).
The aim of this study was to better understand the relationships between tumor characteristics and the probability of axillary LN positivity. The large cohort used in our study is appropriate for less frequent tumor phenotypes (namely Her2+ and HR-Her2-). We distinguished between various histological tumor types, showing a lower LN positivity rate in tumors other than ductal, lobular or mixt, as previously reported for BC with favorable histology (tubular, mucinous, papillary, medullary, adenoid cystic and secretory) that are associated with a very low LN positivity rate .
In our model, we used the same independent variables as Reyal et al. , namely age, tumor size, molecular subtypes and LVI, and we added grade and histological type. However, age intervals were different, as well as tumor phenotype definitions (ER only in the Reyal model) and tumor size description (continuous variable in the Reyal model). We obtained different odds ratios for the same variables and clinical utility results were different and higher for low probability of positive lymph node, particularly for macro metastases in our population for both models. Clinical utility results for low probability of positive lymph node could be contributive to avoid surgical axillary staging by sentinel lymph node biopsy or axillary lymph node dissection.
The models were less reliable when information about LVI was missing. LVI could be detected on pre-operative biopsies but the difference in accuracy is obviously large in comparison with surgical specimen analysis.
The HER2 status was unknown in old studies  and others studies were based on small number of patients. We found that HER2 negative tumors were associated with LN positivity less frequently than HER2 positive tumors (22.9% vs. 31.9%). Lu et al. published that the lowest probability of node metastasis was for ER- / HER2- tumors . Similarly in our study, triple negative tumors had the lowest probability of node metastasis, while HR- / Her2+ tumors had the highest probability. Reyal et al. hypothesized that the axillary LN metastatic process is predominantly related to intrinsic biological properties in ER-negative and HER2-negative BC, while tumor size, proliferation rate and LVI are the main determinants in the ER positive or HER2 positive breast cancers. However, positive axillary lymph nodes in triple negative BC were pejorative prognostic factors for sentinel node macro-metastases but also for occult sentinel node involvement (pN0(i+) and pN1mi) .
A reliable predictive model of LN positivity, based on pathologic parameters, can be used to compare populations from different studies, particularly for trials with or without axillary surgical procedure. Above all, it might allow avoiding SN biopsy when the probability of positivity is very low (< 10%). Some authors already suggested that SN biopsy could be omitted in tumors with good-prognosis subtypes  or that axillary dissection is useless in older patients . We believe that these criteria lack accuracy and we prefer a decision-making approach, based on molecular subtypes. However, we must be aware of the risk of insufficient treatment in small tumors with favorable prognostic factors, in which LN status is a major determinant of adjuvant chemotherapy and regional radiotherapy. Moreover, the model is less reliable when LVI is not documented, which is usually the case before surgery. Ultra-sonography of the axilla and percutaneous biopsy is a growing practice. These clinical predictive tools may be helpful relative to the use of axillary ultra-sonography with percutaneous LN biopsy for patients with high level risk of axillary LN involvement.
These models can also be contributive in order to determined indications of post mastectomy radiotherapy for patients with axillary lymph nodes macro-metastases , particularly when immediate breast reconstruction can be proposed.
We reported a reliable predictive model of LN positivity according to different early breast carcinoma phenotypes in a large cohort. The determination of the risk of positive axillary LN can significantly contribute to therapeutic decisions. These models, with or without LVI results, can also be used to determine the risk of positive axillary LN or the risk of LN macro-metastasis. Before surgery, clinical models can be used to propose SLNB or not according to LN involvement probability. After surgery, in case of SLNB omission, if LN involvement probability is high, with eventually modifications of adjuvant treatment indications according to LN status, a re-operation can be proposed (SLNB or cALND). Thus clinical and pathologic models should be helpful in surgical planning, in the setting of a clinical trial and in clinical practice to avoid SLNB for very low risk of LN involvement and to avoid re-operation in case of SLNB omission or to propose ALND for patients with high level probability of major axillary LN involvement but also to propose immediate breast reconstruction when PMRT is not required for.
Lyman GH, Giuliano AE, Somerfield MR, Benson AB 3rd, Bodurka DC, Burstein HJ, et al. American Society of Clinical Oncology guideline recommendations for sentinel lymph node biopsy in early-stage breast cancer. J Clin Oncol. 2005;23:7703–20.
Krag DN, Anderson SJ, Julian TB, Brown AM, Harlow SP, Costantino JP, et al. Sentinel-lymph-node resection compared with conventional axillary-lymph-node dissection in clinically node-negative patients with breast cancer: overall survival findings from the NSABP B-32 randomised phase 3 trial. Lancet Oncol. 2010;11:927–33.
Giuliano AE, Hunt KK, Ballman KV, Beitsch PD, Whitworth PW. Blumencranz PWet al. Axillary dissection vs no axillary dissection in women with invasive breast cancer and sentinel node metastasis: a randomized clinical trial. JAMA. 2011;305:569–75.
Galimberti V, Cole BF, Zurrida S, Viale G, Luini A, Veronesi P, et al. Axillary dissection versus no axillary dissection in patients with sentinel-node micrometastases (IBCSG 23-01): a phase 3 randomised controlled trial. Lancet Oncol. 2013;14:297–305.
Friedman D, Gipponi M, Murelli F, Meszaros P, Solari N, Massa M, et al. Predictive factors of non-sentinel lymph node involvement in patients with invasive breast cancer and sentinel node micrometastases. Anticancer Res. 2013;33:4509–14.
Rivadeneira DE, Simmons RM, Christos PJ, Hanna K, Daly JM, Osborne MP. Predictive factors associated with axillary lymph node metastases in T1a and T1b breast carcinomas: analysis in more than 900 patients. J Am Coll Surg. 2000;191:1–6 discussion 6-8.
González-Vela MC, Garijo MF, Fernández FA, Buelta L, Val-Bernal JF. Predictors of axillary lymph node metastases in patients with invasive breast carcinoma by a combination of classical and biological prognostic factors. Pathol Res Pract. 1999;195:611–8.
Fehm T, Maul H, Gebauer S, Scharf A, Baier P, Sohn C, et al. Prediction of axillary lymph node status of breast cancer patients by tumor biological factors of the primary tumor. Strahlenther Onkol. 2005;181:580–6.
Yip CH, Taib NA, Tan GH, Ng KL, Yoong BK, Choo WY. Predictors of axillary lymph node metastases in breast cancer: is there a role for minimal axillary surgery? World J Surg. 2009;33:54–7.
Bader AA, Tio J, Petru E, Bühner M, Pfahlberg A, Volkholz H, et al. T1 breast cancer: identification of patients at low risk of axillary lymph node metastases. Breast Cancer Res Treat. 2002;76:11–7.
Reyal F, Rouzier R, Depont-Hazelzet B, Bollet MA, Pierga JY, Alran S, et al. The molecular subtype classification is a determinant of sentinel node positivity in early breast carcinoma. PLoS One. 2011;6:e20297.
Lu X, Lu X, Wang ZC, Iglehart JD, Zhang X, Richardson AL. Predicting features of breast cancer with gene expression patterns. Breast Cancer Res Treat. 2008;108:191–201.
Crabb SJ, Cheang MC, Leung S, Immonen T, Nielsen TO, Huntsman DD, et al. Basal breast cancer molecular subtype predicts for lower incidence of axillary lymph node metastases in primary breast cancer. Clinical Breast Cancer. 2008;8:249–56.
Van Calster B, Vanden Bempt I, Drijkoningen M, Pochet N, Cheng J, Van Huffel S, et al. Axillary lymph node status of operable breast cancers by combined steroid receptor and HER-2 status: triple positive tumours are more likely lymph node positive. Breast Cancer Res Treat. 2009;113:181–7.
Lee JH, Kim SH, Suh YJ, Shim BY, Kim HK. Predictors of axillary lymph node metastases (ALNM) in a Korean population with T1-2 breast carcinoma: triple negative breast cancer has a high incidence of ALNM irrespective of the tumor size. Cancer Res Treat. 2010;42:30–6.
Kim MJ, Ro JY, Ahn SH, Kim HH, Kim SB, Gong G. Clinicopathologic significance of the basal-like subtype of breast cancer: a comparison with hormone receptor and Her2/neu-overexpressing phenotypes. Hum Pathol. 2006;37:1217–26.
Ngo C, Mouttet D, De Rycke Y, Reyal F, Fourchotte V, Hugonnet F, et al. Validation over time of a nomogram including HER2 status to predict the sentinel node positivity in early breast carcinoma. Eur J Surg Oncol. 2012;38:1211–7.
Goldhirsch A, Wood WC, Coates AS, Gelber RD, Thürlimann B, Senn HJ, et al. Strategies for subtypes--dealing with the diversity of breast cancer: highlights of the St. Gallen international expert consensus on the primary therapy of early breast Cancer 2011. Ann Oncol. 2011;22:1736–47.
von Minckwitz G, Untch M, Blohmer JU, Costa SD, Eidtmann H, Fasching PA, et al. Definition and impact of pathologic complete response on prognosis after neoadjuvant chemotherapy in various intrinsic breast cancer subtypes. J Clin Oncol. 2012;30:1796–804.
Houvenaeghel G, Goncalves A, Classe JM, Garbay JR, Giard S, Charytensky H, et al. Characteristics and clinical outcome of T1 breast cancer: a multicenter retrospective cohort study. Ann Oncol. 2014.
Harrell FE. Regression Modeling Strategies with Applications to Linear Models, Logistic Regression and Survival Analysis. New York: Springer-Verlag; 2001. p. 568. ISBN 0–387–95232-2.
Gerds TA, Cai T, Schumacher M. The performance of risk prediction models. Biom J. 2008;(4):457–79.
Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez J-C, Müller M. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.
Hosmer DW Jr., Lemeshow S. Applied logistic regression, Second Edition. New York: Wiley; 2000.
Mendez JE, Fey JV, Cody H, Borgen PI, Sclafani LM, et al. Can sentinel lymph node biopsy be omitted in patients with favorable breast cancer histology? Ann Surg Oncol. 2005;12:24–8.
Houvenaeghel G, Sabatier R, Reyal F, Classe JM, Giard S, Charitansky H, Rouzier R, Faure C, Garbay JR, Daraï E, Hudry D, Gimbergues P, Villet R, Lambaudie E. Axillary lymph node micrometastases decrease triple-negative early breast cancer survival. Br J Cancer. 2016;115(9):1024–31.
Martelli G, Miceli R, Daidone MG, Vetrella G, Cerrotta AM, Piromalli D, et al. Axillary dissection versus no axillary dissection in elderly patients with breast cancer and no palpable axillary nodes: results after 15 years of follow-up. Ann Surg Oncol. 2011;18:125–33.
Forissier V, Tallet A, Cohen M, Classe JM, Reyal F, Chopin N, Mazouni C, Gimbergues P, Daraï E, Colombo PE, Azuar P, Lambaudie E, Houvenaeghel G. Is post-mastectomy radiation therapy contributive in pN0-1mi breast cancer patients? Results of a French multi-centric cohort. Eur J Cancer. 2017;87:47–57.
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Availability of data and materials
The datasets generated and analyzed during the current study are available in the ClinicalTrials.gov (Identifier: NCT02869607) repository, https://clinicaltrials.gov/ct2/show/NCT02869607?cond=Breast+Cancer&cntry=FR&city=Marseille&rank=2
Ethics approval and consent to participate
This work was approved by our institutional review board (IPC - Comité d’Orientation Stratégique). All procedures performed in this study involving human participants were done in accordance with the French ethical standards and with the 2008 Helsinki declaration. All included patients provided written informed consent.
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.
Table S1. Logistic regression results. (DOCX 22 kb)
Figure S1. Calibration plots of our models. 1A: Calibration of predictive LN Involvement for validation set– Pathologic model. 1B: Calibration of predictive LN macro metastases for validation set – Pathologic model. 1C: Calibration of predictive LN Involvement for validation set – Clinical model. 1D: Calibration of predictive LN macro metastases for validation set – Clinical model. (DOCX 45 kb)
Figure S2. ROC curves of our models. 2A: ROC curves of predictive LN Involvement – Pathologic model. 2B: ROC curves of predictive LN macro metastases – Pathologic model. 2C: ROC curves of predictive LN Involvement – Clinical model. 2D: ROC curves of predictive LN macro metastases – Clinical model. (DOCX 44 kb)
About this article
Cite this article
Houvenaeghel, G., Lambaudie, E., Classe, JM. et al. Lymph node positivity in different early breast carcinoma phenotypes: a predictive model. BMC Cancer 19, 45 (2019). https://doi.org/10.1186/s12885-018-5227-3
- Breast cancer
- Sentinel node
- Risk prediction
- Molecular subtype