Treatment response correlation between primary tumor and axillary lymph nodes after neoadjuvant therapy in breast cancer: a retrospective study based on real-world data
Introduction
Breast cancer is the most commonly diagnosed malignancy in women and accounts for more female deaths than any other cancer worldwide (1). So far, a wide variety of trials have been designed and performed for treatments of breast cancer. Notably, neoadjuvant therapy (i.e., preoperative therapy) has been acknowledged for its unique and important position in the comprehensive treatment of breast cancer. The goal of neoadjuvant therapy for breast cancer mainly lies in downstaging the tumor, assessing the sensitivity to systemic therapy regimens (such as chemotherapy) beforehand, improving the possibility of breast-conserving surgery, and predicting the long-term risk of recurrence (2-5).
The absence of residual invasive tumor in the breast and lymph nodes after neoadjuvant therapy, which is widely known as pathologic complete response (pCR), often suggests a better prognosis than residual disease (6,7). The rate of pCR varies among the different subtypes of breast cancer, with approximately 30–50% of human epidermal growth factor receptor 2 (HER2)-positive and triple-negative (defined as hormone receptor-negative and HER2-negative) breast cancer (TNBC) patients achieving pCR after standard neoadjuvant therapy, compared with less than 20% of patients with hormone receptor (HR)-positive breast cancer (8,9). Moreover, the efficacy of neoadjuvant therapy, or the pCR rate, is considered to be related to the characteristics of the patient and tumor before treatment (10,11). For instance, patients with smaller tumors (cT1 or cT2) usually experience higher rates of pCR. In China, however, data regarding the current situation of neoadjuvant therapy are limited, especially the treatment response and corresponding clinicopathological characteristics for breast cancer. Therefore, neoadjuvant therapy for breast cancer patients in China requires further investigation.
At present, in China, the routine management of axillary lymph nodes after neoadjuvant therapy is axillary lymph nodes dissection (ALND), which unavoidably increases the incidence of complications related to ALND, such as lymphedema, shoulder dysfunction, seroma formation, and loss of sensation in the distribution of the intercostobrachial nerve, which negatively impact quality of life (QOL) (12). Thus, identifying exceptional responders after neoadjuvant therapy, who potentially may not require axillary surgery, meets the requirements of precision medicine. Notably, the findings from an early retrospective study also suggested that the invasive tumor content was more likely to be eradicated from the primary tumor if there was no measurable disease in the axillary lymph nodes (13). On the basis of this observation, it is reasonable to hypothesize that an excellent response of the primary tumor after neoadjuvant therapy, or no residual invasive tumor, may indicate a better axillary status (no or few metastases in the axillary lymph nodes) in patients with breast cancer. Interestingly, two studies from MD Anderson Cancer Center suggest that such a correlation most likely exists, and proposed that it is rational to omit axillary surgery in specific populations (cN0 HER2-positive or TNBC diseases with no residual invasive content in the breast after neoadjuvant chemotherapy) treated with neoadjuvant therapy (14,15). In particular, Barron et al. and Tadros et al. reported the rate of pathologic nodal positivity to be extremely low (less than 2%) in HER2-positive or TNBC patients with cN0 diseases (14,15). However, these two studies focused on the population registered in the National Cancer Database (NCDB) of the United States, and the exact characteristics of such patients in the Chinese population remain undetermined.
In this report, we used real-world data from a comprehensive breast cancer institution in western China to investigate the associated factors that may play a role in the efficacy of neoadjuvant therapy for the primary tumor and axillary lymph nodes, respectively. Furthermore, the treatment response correlation between the primary tumor and axillary lymph nodes after neoadjuvant therapy was also analyzed with the aim of guiding the screening of different candidates with corresponding specific characteristics.
We present the following article in accordance with the STROBE reporting checklist (available at http://dx.doi.org/10.21037/gs-20-686).
Methods
Study design and participants
From January 1, 2012 to December 31, 2018, this retrospective study collected and reviewed data from a comprehensive breast cancer institution in western China, the First Affiliated Hospital of Chongqing Medical University. The institutional ethics committee reviewed and approved this study (No. 2020-59). Due to the retrospective nature of the research and the concealment of patient information, the need for informed consent was waived. The study was performed in accordance with the Declaration of Helsinki (as revised in 2013).
Breast cancer patients who received standard neoadjuvant therapy (e.g., anthracycline- and/or taxane-based neoadjuvant chemotherapy, anti-HER2 therapy, neoadjuvant endocrine therapy) followed by breast and axillary surgery were identified. Patients with metastatic disease (M1) were excluded. All patients were diagnosed by core needle biopsy before the initiation of neoadjuvant treatment at our hospital. The suspected axillary lymph nodes were assessed by core needle biopsy to obtain histological evidence before neoadjuvant therapy, or by clinical physical examination or ultrasonography. The following patient demographic and clinical characteristics at baseline were collected: age, menstrual status, lesion side, clinical TNM (cTNM) staging, pathological reports of biopsy [including histology, immunohistochemical analysis, or fluorescence in situ hybridization (FISH)], treatment regimen, and course of chemotherapy.
Definitions
A pCR of the breast primary tumor was defined as the absence of residual invasive carcinoma in the excisional breast specimen, while a pCR of the axillary lymph nodes was defined as the absence of measurable or metastatic disease in the excisional specimens of axillary lymph nodes.
Statistical analysis
Univariate analysis (Pearson’s chi-squared test for categorical variables and non-parametric Wilcoxon rank-sum test for ordinal variables) was used to determine the correlation between baseline characteristics and treatment response of the primary tumor and axillary lymph nodes, respectively. To estimate the effect of these factors on the response of the primary tumor and axillary lymph nodes, in line with the widely used recommendations (16), variables with a P value lower than 0.25 in the univariate analysis and those with clinical importance were selected as independent variables and analyzed by multivariate logistic regression. All tests were two-sided, and P values lower than 0.05 were regarded as being statistically significant (P<0.05).
Association rules analyses were conducted to investigate the association patterns between primary tumor responses and the responses of axillary lymph nodes, as well as to identify the corresponding characteristics of the patients and tumors. Apriori is the most frequently applied machine-learning algorithm for association rules mining (17), which is an effective statistical procedure in which association patterns in data are identified and presented in the form of a series of rules (18). Typically, mined association rules comprise two parts: an antecedent (or left-hand-side, LHS) and a consequent (or right-hand-side, RHS). To describe and represent the interest and significance of a rule, basic metrics (support, confidence, and lift) are utilized. Support is defined as the frequency of one rule occurrence in the total dataset, while confidence is the probability of the RHS in cases of fulfilling the LHS, thus reflecting the reliability of the inference drawn by that rule (LHS → RHS). Additionally, as a measure of significance, lift is a factor defined as the realistic probability of co-occurrence of LHS and RHS in a rule divided by the expected probability of LHS and RHS co-occurring, assuming that both of them are independent. Therefore, lift measures the strength of one rule, which can be interpreted as effective association when it is greater than 1. Moreover, the higher the lift, the more significant the rule is. In the current analyses, association rules mining was performed under the condition of minimum support of 0.05 and minimum confidence of 0.75, according to the recommendations (19).
Univariate analysis and multivariate logistic regression were performed using Statistics Analysis Software (version 9.4, SAS Institute Inc, Cary, NC, USA). Association rules analyses were implemented in the arules R package of R software (version 3.5.2).
Results
A total of 1,312 patients were included in the final analysis. Table 1 summarizes the characteristics of the patients and tumors at baseline. Of the included patients [median age, 49 (IQR 43–56) years], 517 patients were diagnosed during the perimenopausal period, 450 were premenopausal, and 345 were postmenopausal. As far as the lesion side of primary tumor, 695 women had left tumors while 612 patients had right tumors. Most of these patients had cT2 (70.8%) or cN1 (45.8%) disease. Regarding the pathological reports, patients with invasive ductal histology accounted for the largest proportion (97.7%), and most of the tumors were classified as grade II (55.8%). In the immunohistochemical analysis or FISH of the excised specimens, the majority of the patients were estrogen receptor (ER)-positive (60.6%), progesterone receptor (PR)-negative (53.8%), and HER2-negative (50.1%), and had high Ki67 expression (72.6%). After being selected as the candidates for neoadjuvant therapy, most of the included patients were treated with chemotherapy (97.9%) and the courses of treatment were mainly less than five cycles (92.5%).
Full table
The types of surgery and treatment response after neoadjuvant therapy are shown in Table 2. Mastectomy (98.1%) was the most common procedure following neoadjuvant therapy, of which modified radical mastectomy was performed most frequently (96.7%). Furthermore, responses of the primary tumors and axillary lymph nodes were also presented as shown in Table 2.
Full table
Association between clinicopathological factors and the response of primary tumors
In the univariate analysis as shown in Table 3, cT and cN classification, histological grade, ER status, PR status, HER2 and Ki67 expression were selected (P<0.25). These factors and those with known clinical importance (menstrual status, lesion side, course of chemotherapy, and regimen) were then included in the multivariate logistic regression analysis presented in Table 4. The final results showed patients with a larger tumor size [cT2: odds ratio (OR) 0.306; cT3: OR 0.182; P<0.05], higher classification of regional lymph nodes (cN1: OR 0.341; cN2-3: OR 0.357; P<0.05), higher histological grade (Grade II: OR 0.184, P<0.05), and positive PR status (OR 0.375, P<0.05) were less likely to achieve breast pCR. Conversely, tumors with a high expression of Ki67 had a greater likelihood of breast pCR than those with low Ki67 expression (OR 1.907, P<0.05).
Full table
Full table
Association between clinicopathological factors and axilla response
When analyzed in the univariate analysis (Table 5), clinicopathological factors including menstrual status, lesion side, cT and cN classification, histology, histological grade, ER and PR status, and Ki67 expression were identified (P<0.25). With the use of multivariate logistic regression analysis (Table 6), axilla response was less likely to achieve pCR when the clinical classification of regional lymph nodes was more advanced (cN1: OR 0.008; cN2-3: OR 0.006; P<0.05). However, higher histological grade was associated with higher possibility of axillary pCR (Grade II: OR 3.446, P<0.05).
Full table
Full table
Association rules analyses between breast response and axilla response
Table 7 shows the top 20 association rules between breast pCR and axilla pCR, sorted by confidence and lift. Overall, there was a treatment response correlation between breast pCR and axilla pCR. Specifically, when a breast pCR was achieved in patients with cT2N0 and HR-negative (ER-negative and PR-negative) diseases, the probability of axilla pCR was more than 90% [mean confidence, 0.9351 (95% CI, 0.9255–0.9447), SE, 0.00478; mean lift, 2.4321(95% CI, 2.4073–2.4569), SE, 0.01239; Table 7 and Table S1].
Full table
As for unfavorable response (no pCR) in both the breast and the axilla, patients diagnosed with cT3N1-2 disease and treated with fewer than five-cycle chemotherapy were highly likely to achieve no pCR for the primary tumor and axillary lymph nodes at the same time [mean confidence, 0.9637 (95% CI, 0.9516–0.9759), SE, 0.00599; mean lift, 1.3245 (95% CI, 1.2485–1.4005), SE, 0.03757; Tables S2-S5]. Again, this finding demonstrated the consistency between breast response and axillary response under certain specific conditions.
According to the pooled analysis of the rules in Tables S2,S3,S6,S7, compared with advanced regional lymph node classifications (cN1 or cN2), cN0 was shown to be a predictor of axillary pCR [mean confidence, 0.9453 (95% CI, 0.9317–0.9588), SE, 0.00670; mean lift, 1.9898 (95% CI, 1.8458–2.1337), SE, 0.07118].
HR-positive (ER-positive and PR-positive) status, HER2-negative status, and low Ki67 expression were highly likely to be the characteristics responsible for inconsistency between the responses of the breast and the axillary lymph nodes [mean confidence, 0.9563 (95% CI, 0.9460–0.9666), SE, 0.00504; mean lift, 1.5570 (95% CI, 1.2976–1.8164), SE, 0.12683; Tables S6,S8,S9).
Discussion
This study discovered the correlation of some clinicopathological factors with treatment response after neoadjuvant therapy in breast cancer. In terms of the response of the primary tumor, patients with a larger tumor size, higher classification of regional lymph nodes, positive PR status, and higher histological grade were less likely to achieve breast pCR. Furthermore, tumors with high Ki67 expression had a greater likelihood of breast pCR than those with low Ki67 expression. In terms of axilla response, when the clinical classification of regional lymph nodes was more advanced, pCR was less likely. Meanwhile, clinical regional lymph node status (cN) was found to be associated with the treatment response in both the primary tumor and axillary lymph nodes; in other words, a pCR of both the primary tumor and axillary lymph nodes was more likely to be achieved in patients with cN0 disease.
The primary findings of this study demonstrate that a strong correlation exists between the response of the primary tumor and the axilla response. Specifically, once a breast pCR was achieved in patients with cT2N0 disease, invasive ductal histology, and negative HR (ER and PR) status after neoadjuvant chemotherapy, axilla pCR was also highly likely. Moreover, axilla pCR was rarely observed if breast pCR was not achieved among patients with cT3N1-2 disease treated with neoadjuvant chemotherapy (treatment course fewer than 5 cycles). It is worth noting that cN0 classification was mined in almost every rule among the association rules analyses between breast pCR and axilla pCR, which further highlights the considerable importance of this factor in predicting the responses of primary tumors and axillary lymph nodes.
Our findings are generally consistent with those of a previous study reported by Tadros et al. (15) and research performed by Barron et al. (14) For instance, Tadros et al. (15) concluded that breast pCR was strongly associated with pathologic nodal status after neoadjuvant therapy and the risk for missing nodal metastases without axillary surgery in selected population (TNBC and HER2-positive breast cancer) was extremely low. Furthermore, Barron et al. (14) pointed out that the highest rates of breast pCR were found in HER2-positive disease and TNBC, and the pathologic nodal positivity rate was extremely low in patients with cN0 HER2-positive disease or TNBC with breast pCR. Both of these studies strengthen the significance of cN0 disease in the correlation between breast pCR and axilla pCR, and highlight that clinical examination and assessment play a significant role and can indicate the response of axillary lymph nodes. Nevertheless, accurate clinical examination relies on the experience and ability of doctors. Therefore, clinical practitioners should master comprehensive knowledge and improve their personal skills.
Although Kuerer et al. pointed out that the treatment responses of the breast and axilla were inclined to be consistent (13), our findings show that there was inconsistency between these two lesions. The most likely potential factors accounting for such inconsistency were HR-positive (ER-positive, PR-positive), HER2-negative, and low Ki67 expression (i.e., luminal A subtype). More importantly, such a finding is in accordance with those of previous studies which suggested that patients with the luminal A subtype were less likely to achieve pCR (in that case, pCR was defined as both the responses of the primary tumor and axillary lymph nodes were confirmed to be pCR) after neoadjuvant therapy than those with HER2-positive or TNBC disease (8,9). However, to the best of our knowledge, this is the first study showing that the responses of the primary tumor and axillary lymph nodes were obviously different in breast cancer with a luminal A subtype by analyzing these two lesions independently, which can be regarded as an advancement of previous studies (8,9). Combined with results from the univariate analysis and multivariate logistic regression analysis above, such luminal A baseline characteristics were more likely to influence the response of the primary tumor. Instead, axilla response was shown to be associated with clinical regional lymph nodes classification. Therefore, it is necessary to investigate the potential mechanisms behind these phenomena above in further studies.
The research teams from MD Anderson Cancer Center believe that the consideration regarding omission of axillary surgery in exceptional responders after neoadjuvant chemotherapy is reasonable (14,15). Along with the findings from this study, we agree that selecting a proportion of exceptional responders as candidates for omission of axillary surgery on the basis of accurate evaluation of the primary breast tumor is feasible. Further studies and trials are desperately needed in order to build a model for predicting axillary response through accurate evaluation of the response of the breast. Therefore, our study also sets the stage for the selection of potential participants for trials concerning the safe omission of axillary surgery to improve the QOL of Chinese breast cancer patients.
Of note, our study also investigated the current situation in a comprehensive breast cancer institution in western China with different populations and treatment regimens, and can serve as a reference for broader comparisons between different countries and nations.
One of the most valuable strengths of this study is its use of association rules mining to evaluate correlations. To our knowledge, this is the first time this method, which has a positive impact on data mining and can serve as a reference for future studies, has been utilized in this field. Our definition of neoadjuvant therapy, which included chemotherapy and endocrine therapy rather than neoadjuvant chemotherapy only, was another merit of this study. Rigorous statistical methods were used to mitigate potential heterogeneity by excluding groups with a small sample size from analysis, making groups more comparable. However, the current study undoubtedly has limitations, despite it being well designed and employing strict statistical methods for analyses. Firstly, there is no doubt that data loss has affected the interpretation of the results. Secondly, the sample size of this study was not very large. Thirdly, due to the single-centered and retrospective nature of our study, selection bias in the study population and recall bias in data collection might have unavoidably been introduced. Therefore, replication of these findings through large-scale multicenter studies with other comprehensive breast cancer institutions is necessary.
This is the first study to attempt to identify candidates for the omission of axillary surgery in order to eventually improve the QOL of Chinese breast cancer patients. The findings and research methods of this study not only have implications for the treatment of breast cancer in China, but they could potentially serve as a reference for different countries and regions. However, there is still a long way to go to improve the QOL among patients, and further work should concentrate on cultivating better clinical skills to facilitate more accurate assessment of the disease, as well as strengthening international cooperation in this area. Finally, more multi-centered and transnational randomized controlled trials are needed for best evidence to guide practice that conforms to the requirements of precision medicine and brings more and more benefits to patients.
Conclusions
In breast cancer patients, clinicopathological factors including primary tumor size, regional lymph nodes, PR status, Ki67 expression, and histological grade were found to be associated with the treatment response of the primary tumor. In regard to the treatment response of the axilla, clinical classification of lymph nodes and histological grade were identified to be the potential influencing factors.
Importantly, our findings suggest a strong correlation between breast pCR and axilla pCR among patients with specific characteristics (cT2N0 disease, negative HR status). Furthermore, axilla pCR was rarely observed if breast pCR was not achieved among patients with cT3N1-2disease. The factors that were most likely to account for the inconsistency between the response of the breast and the axillary lymph nodes were found to be HR-positive, HER2-negative, and low Ki67 expression. Further studies are needed to build a model to predict axillary response by accurately evaluating the response of the breast first.
Acknowledgments
The authors are grateful to Professor Xiaoni Zhong, from the School of Public Health and Management in Chongqing Medical University, for the critical contributions to the statistical design and analysis in this study. We also thank AME Editing Service for the extensive English language editing of this paper.
Funding: This work was supported by grants from the National Natural Science Foundation of China (No.81772979 and 81472658). The funders of the study had no role in study design, data collection, data analysis, data interpretation, or the writing of the report.
Footnote
Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at http://dx.doi.org/10.21037/gs-20-686
Data Sharing Statement: Available at http://dx.doi.org/10.21037/gs-20-686
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at http://dx.doi.org/10.21037/gs-20-686). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The institutional ethics committee of the First Affiliated Hospital of Chongqing Medical University reviewed and approved this study (No. 2020-59). Due to the retrospective nature of the research and the concealment of patient information, the need for informed consent was waived.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin 2020;70:7-30. [Crossref] [PubMed]
- Mamtani A, Barrio AV, King TA, et al. How Often Does Neoadjuvant Chemotherapy Avoid Axillary Dissection in Patients With Histologicalally Confirmed Nodal Metastases? Results of a Prospective Study. Ann Surg Oncol 2016;23:3467-74. [Crossref] [PubMed]
- Kaufmann M, von Minckwitz G, Bear HD, et al. Recommendations from an international expert panel on the use of neoadjuvant (primary) systemic treatment of operable breast cancer: new perspectives 2006. Ann Oncol 2007;18:1927-34. [Crossref] [PubMed]
- Schwartz GF, Hortobagyi GN. Proceedings of the consensus conference on neoadjuvant chemotherapy in carcinoma of the breast, April 26-28, 2003, Philadelphia, Pennsylvania. Cancer 2004;100:2512-32. [Crossref] [PubMed]
- Shannon C, Smith I. Is there still a role for neoadjuvant therapy in breast cancer? Crit Rev Oncol Hematol 2003;45:77-90. [Crossref] [PubMed]
- Prowell TM, Pazdur R. Pathological complete response and accelerated drug approval in early breast cancer. N Engl J Med 2012;366:2438-41. [Crossref] [PubMed]
- Cortazar P, Zhang L, Untch M, et al. Pathological complete response and long-term clinical benefit in breast cancer: the CTNeoBC pooled analysis. Lancet 2014;384:164-72. [Crossref] [PubMed]
- van la Parra RF, Kuerer HM. Selective elimination of breast cancer surgery in exceptional responders: historical perspective and current trials. Breast Cancer Res 2016;18:28. [Crossref] [PubMed]
- Guarneri V, Broglio K, Kau SW, et al. Prognostic value of pathologic complete response after primary chemotherapy in relation to hormone receptor status and other factors. J Clin Oncol 2006;24:1037-44. [Crossref] [PubMed]
- von Minckwitz G, Untch M, Blohmer JU, 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. [Crossref] [PubMed]
- Houssami N, Macaskill P, von Minckwitz G, et al. Meta-analysis of the association of breast cancer subtype and pathologic complete response to neoadjuvant chemotherapy. Eur J Cancer 2012;48:3342-54. [Crossref] [PubMed]
- Kwan W, Jackson J, Weir LM, et al. Chronic arm morbidity after curative breast cancer treatment: prevalence and impact on quality of life. J Clin Oncol 2002;20:4242-8. [Crossref] [PubMed]
- Kuerer HM, Sahin AA, Hunt KK, et al. Incidence and impact of documented eradication of breast cancer axillary lymph node metastases before surgery in patients treated with neoadjuvant chemotherapy. Ann Surg 1999;230:72-8. [Crossref] [PubMed]
- Barron AU, Hoskin TL, Day CN, et al. Association of Low Nodal Positivity Rate Among Patients With ERBB2-Positive or Triple-Negative Breast Cancer and Breast Pathologic Complete Response to Neoadjuvant Chemotherapy. JAMA Surg 2018;153:1120-6. [Crossref] [PubMed]
- Tadros AB, Yang WT, Krishnamurthy S, et al. Identification of Patients With Documented Pathologic Complete Response in the Breast After Neoadjuvant Chemotherapy for Omission of Axillary Surgery. JAMA Surg 2017;152:665-70. [Crossref] [PubMed]
- Hosmer DW, Lemeshow S. Model-Building Strategies and Methods for Logistic Regression. In: Shewhart WA, Wilks SS, eds. Applied Logistic Regression, 2nd edition: John Wiley & Sons, 2000:91-142.
- Agrawal R, Srikant R. Fast Algorithms for Mining Association Rules in Large Databases. In: Proceedings of the 20th International Conference on Very Large Data Bases, 1994.
- Vougas K, Sakellaropoulos T, Kotsinas A, et al. Machine learning and data mining frameworks for predicting drug response in cancer: An overview and a novel in silico screening process based on association rule mining. Pharmacol Ther 2019;203:107395. [Crossref] [PubMed]
- Hornik K, Grün B, Hahsler M. arules - A Computational Environment for Mining Association Rules and Frequent Item Sets. J Stat Softw 2005;14:1-25.
- Amin MB, Edge S, Greene F, et al. AJCC Cancer Staging Manual, Eight Edition. New York: Springer, 2017.
- Ogston KN, Miller ID, Payne S, et al. A new histological grading system to assess response of breast cancers to primary chemotherapy: prognostic significance and survival. Breast 2003;12:320-7. [Crossref] [PubMed]
- Eisenhauer EA, Therasse P, Bogaerts J, et al. New response evaluation criteria in solid tumours: revised RECIST guideline (version 1.1). Eur J Cancer 2009;45:228-47. [Crossref] [PubMed]