- 1Department of Radiation Oncology, Key Laboratory of Cancer Prevention and Intervention, Ministry of Education, Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
- 2General Surgery, Cancer Center, Department of Breast Surgery, Zhejiang Provincial People’s Hospital (Affiliated People’s Hospital, Hangzhou Medical College), Hangzhou, Zhejiang, China
- 3Department of Pathology, Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
Objective: We aim to explore the clinicopathological features associated with axillary node response and recurrence in breast cancer patients undergoing neoadjuvant treatment (NAT).
Methods: We retrospectively reviewed the medical records of 486 stage I to III breast cancer patients who received NAT and surgery between 2016 and 2021.
Results: A total of 486 cases were reviewed and 154 (31.7%) patients achieved breast pathological complete response (pCR) (ypT0/Tis). Of the 366 cases with initially cN+, 177 (48.4%) cases reach ypN0. Breast pCR is in high accordance to axillary pCR (81.5%). Hormone receptor (HR)-/HER2+ breast cancer patients have the highest axillary pCR rate (78.3%). Patients achieve axillary pCR have a significantly better disease-free survival (DFS) (P=0.0004). Further analysis reveals that the DFS of ypN0 and ypN1 cases are similar (P=0.9049). Moreover, DFS in patients with ypN0 (P<0.0001) and ypN1 (P<0.0001) is significantly better than that in patients with ypN2-3. For post-mastectomy ypN0 cases, radiation could only improve DFS in patients with initially cN+ stage (P=0.0499). Multivariate Cox regression analysis shows that radiation is an independent factor to improve DFS (Hazard ratio (HR): 0.288(0.098-0.841), P=0.0230). Radiation does not improve DFS in pre-cN0/ypN0 patients (P=0.1696).
Conclusion: Axillary pCR rate is higher than breast pCR rate. HR-/HER2+ patients have the highest axillary pCR rate. Axillary pCR is associated with better DFS. Radiation could further improve DFS in ypN0 patients with initially positive nodal disease.
Introduction
Neoadjuvant treatment (NAT) has widely been used in the treatment of early-stage breast cancer. Pathological complete response (pCR), usually defined as ypT0/Tis with or without ypN0, is not the only aim of NAT but undeniably an outstanding marker for predicting prognosis (1). It has been reported that nodal status after NAT but not initially N stage was significantly associated with survival. Axillary pCR is a potential prognostic indicator of better 5-year overall survival (OS) (2). Breast pCR is highly consistent with axillary pCR, especially in HER2+ and triple negative breast cancer (TNBC) (3, 4). The HER2+ breast cancer and TNBC were also reported to have the highest axillary pCR rate (5). In these two subtypes of breast cancer, the risk of residual nodal metastasis in breast pCR and cN0/ycN0 patients is quite low. Therefore, it may be possible to avoid axillary lymph node dissection (ALND) (3, 6). However, axillary pCR is only a moderate predictor of breast pCR, indicating that the response of breast and axillary to NAT is not identical (4).
The pre- and post-NAT nodal stages are important considerations for the use adjuvant radiation treatment (RT). However, it is still various for clinical practice of adjuvant RT in patients with ypN0 (7, 8). In all breast cancer patients after NAT, the use of adjuvant RT is associated with a trend of lower local relapse (8). Data from National Cancer Data Base (NCDB) indicated that adjuvant RT conferred significant overall survival (OS) improvement in ypN1 patients. For ypN0 cases, RT was only associated with OS benefits in post-lumpectomy patients, but not in post-mastectomy ones (7). However, other studies have reported that adjuvant RT could improve OS in all ypN stages subgroups after mastectomy (9). Moreover, the regional nodal irradiation shows no extra survival improvement compared with breast-only RT (7, 9, 10).
We aim to explore the clinicopathological features associated with axillary response and prognosis of post-NAT patients with different ypN stages. We also aim to investigate patients who benefited from post-mastectomy RT (PMRT) most. Here, we present our findings.
Materials and methods
Patients
The data of female breast cancer patients who received NAT and surgery in the Second Affiliated Hospital of Zhejiang University (SAHZU) and Zhejiang Provincial People’s Hospital (ZPPH) from 2016 to 2021 were retrospectively reviewed. The inclusion criteria are as follows: (1) biopsy was performed before treatment, and breast cancer was confirmed by histopathology; (2) patients underwent NAT and post NAT surgery in these two medical centers; (3) follow-up time was longer than one month. Exclusion criteria included: (1) patients initially diagnosed as stage IV; (2) lack of necessary data in medical records; (3) bilateral breast cancer. This study was approved by the Ethics Committee of SAHZU and ZPPH. The study was carried out in accordance with the Declaration of Helsinki.
Statistical analysis
Disease-free survival (DFS) was calculated from the date of surgery to the first local or distant recurrence, death or last follow-up without recurrence. The differences of distribution were analyzed using χ2 or Fisher exact test. A modified Poisson regression model was used to estimate the relative risk between factors and axillary pCR rate. The survival curves were drawn by Kaplan–Meier method. Log-rank P was calculated to compare the differences between survival curves. The hazard ratio (HR) and 95% confidence interval (CI) were calculated using the univariable and multivariable Cox models. The figures were drawn using Graph Prism 9.3.0 and analyses were performed using IBM SPSS for Windows version 26.0. P<0.05 was considered statistically significant.
Results
A total of 486 cases were included in this study. The median age of the cohort was 51 years (24-75 years). The median follow-up time was 25 months (range, 1-70 months). According to immunohistochemistry (IHC) and fluorescence in situ hybridization (FISH) results of breast and lymph node biopsy specimen, 176 (36.2%) cases were categorized as hormone receptor (HR)+/HER2- (Luminal A) type. 140 (28.8%) cases were HR+/HER2- (Luminal B). 80 (16.5%) cases were HER2 amplified and 90 (18.5%) cases were TNBC. Most patients (n=337, 70.4%) were staged as initially cT2 and cT3-4 cases (n=72) accounted for only 15.0%. T stages were unknown in 7 cases. 366 cases (75.3%) were initially cN+ (325 cases were pathologically proven) and 120 cases (24.7%) were initially cN0 (70 cases were pathologically proven). In our cohort, 372 cases (76.5%) underwent mastectomy and 433 cases (89.1%) underwent axillary lymph node dissection (ALND). In total, 154 (31.7%) patients achieved breast pCR (ypT0/Tis) and 269 (55.3%) patients achieved axillary pCR (ypN0). Characteristics of cohort are summarized in Table 1.
Of the 366 pre-cN+ cases, 177 (48.4%) achieved axillary pCR. However, 28 pre-cN0 cases (28/120, 23.3%) turned out to be ypN+. The luminal B (63.2%) and HER2 amplified (78.3%) patients had the highest axillary pCR rate. For all HER2 positive patients, the axillary response rate was 69.1% (vs HER2-: 29.3%, P<0.0001). Estrogen receptor (ER)- (P<0.0001), progesterone receptor (PR) - (P=0.0002) and Ki67>15% (P=0.019) of primary lesions were associated with higher axillary pCR rate (Table 2). The result of Poisson regression analysis indicated that ER positive patients were less likely to achieve axillary pCR (Relative risk (RR) and 95%CI: 1.332(1.077-1.646), P=0.008). HER2 positive patients were more likely to be ypN0 (RR and 95%CI: 0.588(0.468-0.738), P=0.000). Patients with Ki67 10-20%+ had the higher axillary pCR rate than Ki67 0-10%+ group (RR and 95%CI: 0.767(0.611-0.964), P=0.023). The patients with Miller-Payne (MP) grade 5 in breast were significantly likely to be ypN0 (RR and 95%CI: 0.316(0.210-0.477), P=0.000), followed with MP4 patients (RR and 95%CI: 0.745(0.572-0.970), P=0.029) (Figure 1).
Table 2 The correlation between ypN stages after NAT and clinicopathological characteristics in pre-cN+ patients (n=366).
Figure 1 The relative risks between clinicopathological characteristics and axillary response in pre-N+ patients (n=366). ER, estrogen receptor; HER2, human epidermal growth factor receptor2; MP, Miller-Payne; LN, lymph node; CI, confidence interval.
The DFS curves of patients with different ypN stages are shown in Figure 2. The ypN0 patients had significantly better DFS than ypN+ patients in both whole cohort (P=0.0004, Figure 2A) and pre-cN+ group (P=0.0129, Figure 2B). Further analysis revealed that the DFS of ypN0 and ypN1 patients were similar (P=0.9049, Figure 2C). However, DFS was significantly better in patients with ypN0 (P<0.0001) and ypN1 (P<0.0001) than in patients with ypN2-3 (Figure 2C). Adjuvant RT was associated with improved DFS in pre-cN+/ypN0 (P=0.0252, Figure 3C) and post-mastectomy pre-cN+/ypN0 patients (P=0.0499, Figure 3D). Multivariate Cox regression analysis showed that PMRT was an independent factor to improve DFS (HR: 0.288(0.098-0.841), P=0.0230, Table 3) in patients with pre-cN+/ypN0. RT also showed a tendency to improve DFS in ypN+ patients (P=0.0759, Figure 3A). However, adjuvant RT did not benefit DFS in preN0/ypN0 patients (P=0.1696, Figure 3B).
Figure 2 The DFS in patients with different ypN stages from the whole cohort and pre-cN+ groups. (A) The DFS of patients with ypN0 and ypN+ stages from the whole cohort (n=486); (B) The DFS of patients with pre-cN+/ypN0 and pre-cN+/ypN+ stages (n=366); (C) The DFS of patients with ypN0, ypN1 and ypN2-3 stages from pre-cN+ group (n=366).
Figure 3 The DFS in patients with or without adjuvant RT with different pre-cN and ypN stages. (A) The association between adjuvant RT and DFS in ypN+ patients (n=217); (B) The association between adjuvant RT and DFS in pre-cN0/ypN0 patients (n=92); (C) The association between adjuvant RT and DFS in pre-cN+/ypN0 patients (n=177); (D) The association between adjuvant RT and DFS in post-mastectomy pre-cN+/ypN0 patients (n=142). RT, radiation treatment.
Table 3 Results of the univariate and multivariate analysis of DFS and clinicopathological characteristics in post-mastectomy pre-cN+/ypN0 patients (n=142).
Discussion
In this retrospective analysis from two hospitals, we investigated potential clinicopathological characteristics associated with axillary response after NAT in breast cancer patients. We also studied the relationship between ypN stage and DFS. RT is usually considered essential for node-positive patients. Therefore, we explored the role of adjuvant RT in patients with different lymph node status before and after NAT.
NAT has been widely used in the treatment of early breast cancer. Overall, patients administrated with neoadjuvant and adjuvant chemotherapy showed similar long-term survival (11). Patients who achieved pCR were more likely to have a significantly better long-terms outcome (1). The response of the breast and axilla to NAT is usually highly consistent. We also observed that the axillary pCR rate varied with the breast MP grade. However, axillary pCR is only a moderate predictor of breast pCR (3, 5, 12). Therefore, the extent of axillary response and ypN stage may play different roles in predicting prognosis and guiding treatment. The HER2 amplified and triple-negative breast cancer subtypes are highly responsive to NAT and have the highest axillary pCR rate (4, 12). In our study, HER2+ breast cancer also has the highest axillary pCR rate (69.1%) and followed with TNBC (45.2%). In our cohort, about 86% of HER2+ patients received HER2-targeted therapy and achieved high axillary response rate. Anti-HER2 therapy has a vital role in improving the response rate. High Ki67 before treatment is a predictor for more pCR in neoadjuvant settings of breast cancer patients (13). Nevertheless, the cut-off value of Ki67 is still controversial. Here we find that patients with ki67 >15%+ in breast biopsy specimens are more likely to achieve axillary pCR. Further analysis indicates that Ki67 10-20%+ patients may be the most axillary responsive group. Hank Schmidt et al. (3). and Diego Flores-Funes et al. (14). reported that larger tumor size at diagnosis was associated with more residual node disease and more ALND. Our data also show similar but insignificant trend. According to our study, the younger patients show a tendency to have higher axillary pCR rate. Better performance status and compliance of the youngers may be part of the reason.
In our study, patients obtain axillary pCR from the whole cohort and initially cN+ group have significantly better DFS. However, the ypN0 and ypN1 patients have similar DFS which is different from the results of T. J. Nijnatten et al. (15). and G.C. Zhang et al. (16). Further analysis shows that about 52.8% of ypN1 patients in our cohort are categorized as Luminal A (compared with 21.7% of ypN0 cases). This indicates that Luminal A patients with cN+ at diagnosis are more likely to have limited nodal burden after NAT. J. W. Chun et al. (17) reported that the 5-year prognosis of patients with residual N1 disease who underwent sentinel lymph node biopsy (SLNB) only was similar to that of patients after ALND. The Luminal A patients may be suitable to received marker clip placement followed with selected lymph node dissection omission of ALND (18). However, more evidence is needed to certify safety. An ongoing clinical trial aims to compare ALND and axillary radiation in SLNB positive T1-3N1M0 patients treated with NAT, which may give us more message in the future (19).
Adjuvant RT had a vital role in the treatment of breast cancer with nodal metastasis. But it is still various in the clinical practice of adjuvant RT in ypN0 patients (7, 8). Y. Y. Zhang et al. (20). reported that PMRT was an independent factor to improve DFS in ypN1-3 patients. For ypN0 patients, PMRT showed no significant benefit. The results from NCDB varied. O.M. Fayanju et al. (7) indicated that adjuvant RT conferred significant OS improvement in ypN1 patients. For ypN0 cases, RT was only associated with OS benefit in post-lumpectomy patients. However, C.G. Rusthoven et al. (9) reported that PMRT could improve OS in patients with all ypN subgroups regardless of surgery type. The results of a meta-analysis suggested that PMRT might reduce local recurrence in patients with ypN0, but would not improve OS (21). In our cohort study, PMRT is an independent factor in improving DFS in pre-cN+/ypN0 patients. For ypN+ patients, adjuvant RT also shows a trend of benefiting DFS. Most of the pre-cN+/ypN0 patients in our cohort receive PMRT covering chest wall, supraclavicular and infraclavicular nodal regions (± internal mammary node). However, several studies have reported that compared with breast-only RT, regional nodal irradiation showed no extra survival improvement in ypN0 patients (7, 9, 10). Our study also has some limits. The sample size of our study is small and medium follow-up time is short, which may affect the conclusions. Besides, the no-RT post-mastectomy group has a higher fraction of HER2+ patients which is the restriction of retrospective study. The results of some ongoing randomized clinical trials may help to guide clinical practice (22).
Conclusions
ER status, PR status, HER2 status and Ki67 are significantly correlated with ypN stage. Breast pCR is highly consistent with axillary pCR. DFS is significantly better in patients achieve axillary pCR with cN+ at diagnosis. The patients with low nodal burden (1-3 nodes) after NAT are more likely to be less aggressive subtypes and have a similar DFS with ypN0 patients. Adjuvant RT is essential in patients with initially cN+/ypN0 stage.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Ethical Committee of the Second Affiliated Hospital of Zhejiang University. Ethical Committee of the Zhejiang Provincial People’s Hospital. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
XR, JW and QW design the work. LL, WX, RN and SW contribute the collection, analysis, or interpretation of data for the work. YY and WX contribute to the follow-up of the patients. SW helps in reviewing the histopathology slides. XR and YY write the manuscript. All authors finally approved the manuscript version to be published.
Funding
This study was supported by the Public Welfare Technology Application Research Project of Zhejiang Province under Grant No. LGF21H160030, and Medical and Health Science and Technology Project of Zhejiang Province (No. 2021KY061, 2023KY046 and 2023KY543).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1093155/full#supplementary-material
References
1. Cortazar P, Zhang L, Untch M, Mehta K, Costantino JP, Wolmark N, et al. Pathological complete response and long-term clinical benefit in breast cancer: The CTNeoBC pooled analysis. Lancet (2014) 384(9938):164–72. doi: 10.1016/S0140-6736(13)62422-8
2. Zetterlund L, Celebioglu F, Hatschek T, Frisell J, de Boniface J. Long-term prognosis in breast cancer is associated with residual disease after neoadjuvant systemic therapy but not with initial nodal status. Br J Surg (2021) 108(5):583–9. doi: 10.1002/bjs.11963
3. Schmidt H, Zhaveri S, Valente C, Pisapati K, Pickholz E, Weltz S, et al. Response in breast vs axilla after neoadjuvant treatment and implications for nonoperative management of invasive breast cancer. Breast J (2021) 27(2):120–5. doi: 10.1111/tbj.14125
4. Lim DW, Greene BD, Look Hong NJ. Relationship between breast and axillary pathologic complete response in women receiving neoadjuvant chemotherapy for breast cancer. Ann Surg Oncol (2021) 28(10):5495–506. doi: 10.1245/s10434-021-10519-8
5. Samiei S, Simons JM, Engelen SME, Beets-Tan RGH, Classe JM, Smidt ML, et al. Axillary pathologic complete response after neoadjuvant systemic therapy by breast cancer subtype in patients with initially clinically node-positive disease: A systematic review and meta-analysis. JAMA Surg (2021) 156(6):e210891. doi: 10.1001/jamasurg.2021.0891
6. Shi Z, Wang X, Qiu P, Liu Y, Zhao T, Sun X, et al. Predictive factors of pathologically node-negative disease for HER2 positive and triple-negative breast cancer after neoadjuvant therapy. Gland Surg (2021) 10(1):166–74. doi: 10.21037/gs-20-573
7. Fayanju OM, Ren Y, Suneja G, Thomas SM, Greenup RA, Plichta JK, et al. Nodal response to neoadjuvant chemotherapy predicts receipt of radiation therapy after breast cancer diagnosis. Int J Radiat Oncol Biol Phys (2020) 106(2):377–89. doi: 10.1016/j.ijrobp.2019.10.039
8. Haffty BG, McCall LM, Ballman KV, Buchholz TA, Hunt KK, Boughey JC. Impact of radiation on locoregional control in women with node-positive breast cancer treated with neoadjuvant chemotherapy and axillary lymph node dissection: Results from ACOSOG Z1071 clinical trial. Int J Radiat Oncol Biol Phys (2019) 105(1):174–82. doi: 10.1016/j.ijrobp.2019.04.038
9. Rusthoven CG, Rabinovitch RA, Jones BL, Koshy M, Amini A, Yeh N, et al. The impact of postmastectomy and regional nodal radiation after neoadjuvant chemotherapy for clinically lymph node-positive breast cancer: A national cancer database (NCDB) analysis. Ann Oncol (2016) 27(5):818–27. doi: 10.1093/annonc/mdw046
10. Schlafstein A, Liu Y, Goyal S, Kahn S, Godette K, Lin J, et al. Regional nodal irradiation for clinically node-positive breast cancer patients with pathologic negative nodes after neoadjuvant chemotherapy. Clin Breast Cancer (2022) 22(2):127–35. doi: 10.1016/j.clbc.2021.06.003
11. Early Breast Cancer Trialists’ Collaborative G. Long-term outcomes for neoadjuvant versus adjuvant chemotherapy in early breast cancer: Meta-analysis of individual patient data from ten randomised trials. Lancet Oncol (2018) 19(1):27–39. doi: 10.1016/S1470-2045(17)30777-5
12. Vugts G, Van den Heuvel F, Maaskant-Braat AJG, Voogd AC, Van Warmerdam LJC, Nieuwenhuijzen GAP, et al. Predicting breast and axillary response after neoadjuvant treatment for breast cancer: The role of histology vs receptor status. Breast J (2018) 24(6):894–901. doi: 10.1111/tbj.13085
13. Chen X, He C, Han D, Zhou M, Wang Q, Tian J, et al. The predictive value of ki-67 before neoadjuvant chemotherapy for breast cancer: A systematic review and meta-analysis. Future Oncol (2017) 13(9):843–57. doi: 10.2217/fon-2016-0420
14. Flores-Funes D, Aguilar-Jimenez J, Martinez-Galvez M, Ibanez-Ibanez MJ, Carrasco-Gonzalez L, Gil-Izquierdo JI, et al. Development of a predictive score of axillary lymph node dissection based on targeted axillary dissection in patients with breast cancer diagnosis, affected lymph nodes, and neoadjuvant treatment. Surg Oncol (2021) 38:101629. doi: 10.1016/j.suronc.2021.101629
15. van Nijnatten TJ, Simons JM, Moossdorff M, de Munck L, Lobbes MB, van der Pol CC, et al. Prognosis of residual axillary disease after neoadjuvant chemotherapy in clinically node-positive breast cancer patients: Isolated tumor cells and micrometastases carry a better prognosis than macrometastases. Breast Cancer Res Treat (2017) 163(1):159–66. doi: 10.1007/s10549-017-4157-0
16. Zhang GC, Zhang YF, Xu FP, Qian XK, Guo ZB, Ren CY, et al. Axillary lymph node status, adjusted for pathologic complete response in breast and axilla after neoadjuvant chemotherapy, predicts differential disease-free survival in breast cancer. Curr Oncol (2013) 20(3):e180-92. doi: 10.3747/co.20.1294
17. Chun JW, Kim J, Chung IY, Ko BS, Kim HJ, Lee JW, et al. Sentinel node biopsy alone for breast cancer patients with residual nodal disease after neoadjuvant chemotherapy. Sci Rep (2021) 11(1):9056. doi: 10.1038/s41598-021-88442-x
18. Almahariq MF, Levitin R, Quinn TJ, Chen PY, Dekhne N, Kiran S, et al. Omission of axillary lymph node dissection is associated with inferior survival in breast cancer patients with residual N1 nodal disease following neoadjuvant chemotherapy. Ann Surg Oncol (2021) 28(2):930–40. doi: 10.1245/s10434-020-08928-2
19. Khan TM, Rossi AJ, Suman V, Haffty B, Hernandez JM, Boughey JC. Is axillary radiation not inferior to axillary dissection for sentinel lymph node-positive breast cancer after neoadjuvant chemotherapy? Ann Surg Oncol (2022) 29(3):1526–7. doi: 10.1245/s10434-021-10830-4
20. Zhang Y, Zhang Y, Liu Z, Qin Z, Li Y, Zhao J, et al. Impact of postmastectomy radiotherapy on locoregional control and disease-free survival in patients with breast cancer treated with neoadjuvant chemotherapy. J Oncol (2021) 2021:6632635. doi: 10.1155/2021/6632635
21. Wang K, Jin X, Wang W, Yu X, Huang J. The role of postmastectomy radiation in patients with ypN0 breast cancer after neoadjuvant chemotherapy: A meta-analysis. BMC Cancer (2021) 21(1):728. doi: 10.1186/s12885-021-08423-1
22. Krug D, Baumann R, Budach W, Dunst J, Feyer P, Fietkau R, et al. Individualization of post-mastectomy radiotherapy and regional nodal irradiation based on treatment response after neoadjuvant chemotherapy for breast cancer: A systematic review. Strahlenther Onkol (2018) 194(7):607–18. doi: 10.1007/s00066-018-1270-x
Keywords: breast cancer, neoadjuvant treatment, radiation, lymph node, prognosis
Citation: Ren X, Yu Y, Liu L, Xia W, Ni R, Wei S, Wu J and Wei Q (2023) Axillary response and outcome in breast cancer patients after neoadjuvant treatment: The role of radiotherapy in reducing recurrence in ypN0 patients with initially cN+ stage. Front. Oncol. 13:1093155. doi: 10.3389/fonc.2023.1093155
Received: 08 November 2022; Accepted: 23 March 2023;
Published: 03 April 2023.
Edited by:
San-Gang Wu, First Affiliated Hospital of Xiamen University, ChinaReviewed by:
Sheldon Feldman, Montefiore Medical Center, United StatesJun Wang, Sichuan University, China
Copyright © 2023 Ren, Yu, Liu, Xia, Ni, Wei, Wu and Wei. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jun Wu, anVud3VAemp1LmVkdS5jbg==; Qichun Wei, cWljaHVuX3dlaUB6anUuZWR1LmNu