Given the higher proportion of axillary lymph node involvement in patients with HER2-positive BC, we performed Firth logistic regressions with nodal status as covariate and confirmed the associations independently of the nodal status (Supplementary Fig. 1 and Supplementary Table 2). This evidence concerns the gene ERBB2 and breast cancer.