In model 1, when adjusted for tumor size, grade, LVI, ER, PR, and HER2, larger tumor size, higher grade, LVI, PR+, and HER2+ were risk factors of axillary lymph nodes metastasis. This evidence concerns the gene ESR1 and neoplasm.