To adjust for potential confounding variables (tumor grade (GIII and GI-II), molecular subtype (HR+, HER2+, HER2+HR+, and TNBC), and age (young and old)), the 555 identified gene were subjected to multivariate Cox proportional hazard model, which identified 176 genes associated with worse prognosis and 218 genes associated with better prognosis (Table S15). Here, ERBB2 is linked to neoplasm.