- 1Haematology Department of Shengjing Hospital, China Medical University, Shenyang, China
- 2Neusoft Research Institute, Northeastern University, Shenyang, China
- 3Department of Hematology, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, China
Objective: Light chain amyloidosis (AL) with cardiac involvement is associated with poor prognosis. The existing prognostic assessment system does not consider treatment-related factors, and there is currently no effective system for predicting the response. The purpose of this study was to build an individualized, dynamic assessment model for cardiac response and overall survival (OS) for AL patients with cardiac involvement.
Methods: The records of 737 AL patients with cardiac involvement were collected through cooperation with 18 hospitals in the Chinese Registration Network for Light-chain Amyloidosis (CRENLA). We used univariate and multivariate analyses to evaluate the prognostic factors for OS and cardiac response. Then, two nomogram models were developed to predict OS and cardiac response in AL patients with cardiac involvement.
Results: A nomogram including four independent factors from the multivariate Cox proportional hazards analysis—Mayo staging, courses of treatment, hematologic response, and cardiac response—was constructed to calculate the possibility of achieving survival by adding all the points associated with four variables. The higher the score, the more likely death would occur. The other nomogram model included the courses of treatment, hematological response, and different treatment regimens, and was correlated with cardiac response. The higher the score, the more likely a cardiac response would occur.
Conclusion: In conclusion, based on the large Chinese cohort of patients with AL and cardiac involvement, we identified nomogram models to predict cardiac response and OS. These models are more individualized and dynamic, and therefore, they have important clinical application value.
Introduction
Light chain amyloidosis (AL) is an uncommon systemic disease (1). Abnormal light chain protein deposition and infiltration leads to dysfunction of target organs, including the heart, kidney, liver, and gastrointestinal tract. The incidence of heart involvement is reported to be almost 80% in AL patients (2). Cardiac amyloidosis typically causes restrictive cardiomyopathy which may involve life-threatening manifestations, such as progressive heart failure and fatal arrhythmias (3). Cardiac involvement is a crucial factor in terms of available treatment options and clinical prognosis (4). Despite recent improvements in the outcomes of AL with novel agents (5), the outcomes of advanced cardiac involvement remain dismal with a median survival time of only 6 months from the onset of heart failure (6). Cardiac involvement greatly affects the outcome; therefore, markers of cardiac injury and dysfunction are powerful prognostic factors in AL patients. At present, the most accepted system for estimating survival in patients with amyloidosis and cardiac involvement is the 2012 Mayo stage, which is based on cardiac troponin I (cTnI), N-terminal natriuretic peptide type B (NT-proBNP) and difference between involved and uninvolved free light chain (dFLC) levels (7). However, the 2012 Mayo stage does not take into account therapeutic factors, such as treatment regimen, cycle, and hematological response, that may influence the disease course and response to treatment. Moreover, the current prediction system cannot be used dynamically according to the patient’s individual condition. This study intends to establish an individualized evaluation system for patients with cardiac amyloidosis based on the initial state of the disease and the dynamics of the treatment process, and to evaluate the patient’s response to cardiac treatment and overall survival (OS).
Nomograms have been used as prediction tools to assess survival and prognosis of patients in multiple diseases, especially patients with tumors (8–10). A nomogram can be individualized to predict outcomes and has important clinical application value. To the best of our knowledge, only a few studies about nomograms for AL patients have been reported and no nomogram models have been established to predict cardiac involvement in AL patients. In the present study, novel nomograms for response and survival were developed based on a nationwide survey including 737 cases of AL with cardiac involvement. These models should be helpful for explaining the heterogeneity in the outcomes of AL patients with diverse treatments, thereby providing individual prediction of the probability of outcome events.
Methods
Patients
The study was approved by the Ethics Committee of Shengjing Hospital, China Medical University. The study was performed in accordance with the principles of the Declaration of Helsinki.
Through cooperation with 18 hospitals in the Chinese Registration Network for Light-chain Amyloidosis (CRENLA), a nationwide survey was conducted from 2009 to 2020, and 1269 patients with systemic AL were registered and followed. A total of 737 (58.1%) patients had cardiac involvement. The diagnosis of systemic AL was confirmed by the presence of Congo red-positive fibril deposition. Cardiac involvement was determined either by detection of amyloid deposits on endomyocardial biopsy or on histologic findings of a biopsy from other tissues in patients with abnormal cardiac findings including imaging and laboratory examinations.
We used the 2012 criteria for Mayo Clinic staging (11, 12). Hematologic and cardiac response after 3 cycles of treatment were assessed according to the Consensus Guidelines (13). All patients were followed up at our clinic or by telephone, and detailed follow-up data were recorded.
Construction and Validation of the Nomogram
Of the 737 patients with cardiac involvement, the records of 553 patients were used to develop the nomogram models as training cohort. A validation cohort of 184 patients was used for further evaluation of the new prognostic models. Training cohort and validation cohort were randomly assigned in a 3:1 ratio. To develop the prognostic nomogram, the clinical features that have been shown to be associated with OS and cardiac response were incorporated: age, duration from onset to diagnosis, light chain, cTnI, dFLC, free light chain ratio (rFLC), Mayo classification 2012, New York Heart Association (NYHA) class, treatment regimen, courses of treatment, hematological response, and cardiac response. We used a univariate analysis to evaluate the prognostic factors for OS and cardiac response. Factors with a P value ≤0.05 were used in the multivariate proportional hazards model for OS and multivariate binary logistic regression model for cardiac response. With the nomograms, which were developed based on the multivariate model parameters, survival analysis was verified by C-index and cardiac response was verified by the area under the curve (AUC).
Statistical Analysis
The baseline clinical features were evaluated using a c2 or nonparametric test, and univariate tests were performed using the Kaplan–Meier method with the log rank test for OS and univariate binary logistic regression for cardiac response. The multivariate proportional hazards model was used to assess the hazard ratios and the multivariate binary logistic regression model was used to assess the odds ratios. Two-sided P values <0.05 were considered statistically significant. All statistical analyses were conducted using SPSS version 23.0 (IBM Corp., Armonk, NY, USA). The nomogram and calibration curve were performed using the Hmisc, rms, and ggplot2 package in R software version 4.0.3 (R Project for Statistical Computing, Vienna, Austria). The final model selection for the nomograms was performed by a forward: LR selection process using a threshold of P.
Results
Baseline Characteristics
Table 1 demonstrates the baseline characteristics of 737 patients with cardiac involvement and 532 patients without cardiac involvement. Of the 737 patients with cardiac systemic AL, 473 (64.2%) were male, and the median age was 59 years (range: 33–85 years). The median duration between onset and diagnosis was 9 months. Immunohistochemical staining identified 537 patients (72.86%) with lambda light-chain isotype and 174 patients (23.61%) with kappa light-chain isotype. Mayo (2012) Stage I, II, III, and IV disease was found in 10.99%, 18.32%, 28.22%, and 21.57% of patients, respectively. The median NT-proBNP and dFLC levels were 3773 ng/L (range: 34–103277) and 203.51 mg/L (0.274–12402.78), respectively. The median cTnI was 0.085 ug/L (0.0–1.945). For chemotherapy, bortezomib-containing treatment was administered to 403 (54.68%) patients and melphalan-containing treatment was administered to 64 (8.68%) patients. Sixty-nine patients (9.36%) chose immunomodulatory drugs, including thalidomide- and lenalidomide-based treatments. Only about 20% of all patients completed more than six courses of treatment.
The overall hematological response rate [complete response (CR) + very good partial response (VGPR) + partial response (PR)] after 3 cycles of treatment was 39.08% (288 cases), including 156 patients with CR (21.17%), 78 patients with VGPR (10.58%), and 54 patients with PR (7.33%). All 737 patients were evaluated for cardiac response after 3 cycles of treatment, and 177 (24.02%) achieved a cardiac response characterized by an NT-proBNP>30% and >300ng/L decrease in patients with baseline NT-proBNP ≥650ng/L or NYHA ≥2 class decrease in subjects with baseline NYHA class 3 or 4.
Nomogram Model Developed to Predict Survival
Baseline demographic, laboratory diagnostic variables, and clinical treatment parameters for cardiac response and OS were further determined by univariate and multivariate Cox regression models. The result of univariate analysis showed that age, light chain type, NT-proBNP, cTnl, dFLC, rFLC, NYHA class, Mayo stage 2012, courses of treatment, therapeutic regimen, hematological response, and cardiac response were related to OS (P<0.05), as shown in Table 2. Subsequently, multivariate Cox regression analyses demonstrated that Mayo stage 2012, hematological response, cardiac response and courses of treatment were independently associated with OS.
A nomogram including all four independent factors from the multivariate Cox proportional hazards analysis was constructed to calculate the possibility of achieving survival by adding all the points associated with the four variables (Figure 1). The higher the score, the more likely death would occur. The nomogram indicated that patients with early Mayo stage, administration of 7–9 courses of treatment, better hematological response, and cardiac response were correlated with better OS.
Nomogram Model Developed to Predict Cardiac Response
In addition to survival time, we also hoped to predict the probability of a patient’s cardiac function remission. The result of the univariate analysis showed that light chain type, NT-proBNP, NYHA class, courses of treatment, therapeutic regimen, and hematological response were related to cardiac response (P<0.05), as shown in Table 3. Subsequently, multivariate Cox regression analyses demonstrated that hematological response, courses of treatment, and treatment regimen were independently associated with cardiac response.
A nomogram including all three independent factors from the multivariate Cox proportional hazards analysis was constructed to calculate the possibility of achieving cardiac response by adding all the points associated with the three variables (Figure 2). The higher the score, the more likely a cardiac response would occur. The nomogram indicated that 7–9 courses of treatment, better hematological response, and bortezomib-containing treatment were correlated with better cardiac response.
Validation of the Nomograms
Validation of the nomograms was processed both internally and externally. Analysis of the internal validation cohort (training cohort) showed a C-index value of 0.871 for nomogram-based predictions of OS. Similarly, in the external validation cohort (validation cohort), the C-index value for predicting OS was 0.800. To validate the nomogram for prediction of cardiac response, the AUC values were calculated. Internal validation showed an AUC value of 0.816 for nomogram-based predictions of cardiac response and the external validation AUC was 0.826.
Discussion
AL is a systemic disease which can damage multiple organs throughout the body. Cardiac involvement is the most common organ-related symptom of AL. The severity of cardiac involvement dictates the prognosis in AL. Currently the Mayo scoring system 2004 or 2012 is widely used for assessing amyloidosis and it mainly predicts survival based on indicators of early myocardial damage, without considering the impact of treatment factors on survival. Moreover, there is no model that can predict the effectiveness of cardiac therapy. Therefore, the purpose of this study was to establish a model for predicting survival and a model for predicting cardiac responsiveness using the largest amyloidosis heart involvement database in China, so as to better guide clinical treatment.
Nomograms are pictorial representations of a complex mathematical formula which are widely used as prognostic devices in the field of medicine. By integrating biological and clinical variables, medical nomograms generate a simple graphical representation of a statistical predictive model and create a numerical probability of a clinical event, such as death (14, 15). Therefore, nomograms can contribute to modern medical decision-making by providing estimated prognosis reliably and conveniently. Moreover, through the nomogram model, each patient will have a different, individualized score calculated to predict the probability of clinical outcomes. This result is unique. In the present study, we developed two nomogram models, one for predicting the OS of AL with cardiac involvement, and the other for predicting the cardiac response.
Among the factors used in the nomogram model to predict OS, four indicators were considered. The first was the Mayo staging system that has been widely used for prediction of OS. Our model then adds three treatment-related factors, which are the courses of treatment, whether hematological remission was achieved after treatment, and cardiac remission. Notably, our data further showed that patients achieving cardiac response have better prognosis compared to patients with no cardiac response. It indicated that not only hematologic response, but also heart function remission, finally translate into survival benefits. According to our model, 7–9 courses of treatment are the most suitable. A greater number of courses of treatment will not bring about additional survival benefits. It may be that more courses of treatment may lead to drug-related toxicity. On the contrary, a lower number of courses of treatment cannot achieve the maximum effect of the treatment. The key point is to explore the optimal dose and treatment interval in order to balance the efficacy and treatment-related adverse effects. If the patient quickly obtains hematological remission and cardiac remission after treatment, this indicates a better prognosis. This combines better efficacy indicators with prognosis. At the same time, the nanogram model quantifies and accumulates multiple factors simultaneously to obtain individualized predictive values, which is more clinically meaningful than the current commonly used single staging system.
Currently, very few published series have described the determinants of cardiac response. Our nomogram prediction model is mainly based on the proportion of patients who already have cardiac involvement. It can be seen that the choice of treatment is of great significance to the function of the heart. First, from the point of view of drug selection, the treatment program based on bortezomib is the most effective for obtaining cardiac remission, followed by the program based on melphalan. Bortezomib-based regimens have been the standard therapy for AL and are recommended by both the National Comprehensive Cancer Network and UK guidelines (16) despite the fact that no randomized clinical trial has shown the efficacy and safety of bortezomib-including regimens to date. One European study previously reported the efficacy of bortezomib-cyclophosphamide-dexamethasone (CyBorD) for the initial treatment in 230 patients with AL. The overall hematologic response rate was 60% (23% complete) and cardiac response was 17% (17). In the treatment of AL, the effect of immunomodulatory drugs on cardiac involvement and recovery seems to be limited (18). Both lenalidomide and thalidomide are immunomodulatory derivatives. In a trial of lenalidomide, melphalan, and dexamethasone with 22 of 25 patients having stage II or III cardiac amyloidosis, the overall hematologic response rate was 58% and organ responses were seen in only 8% of patients. Cardiac arrhythmias were seen in 33% (18). Indeed, lenalidomide is frequently observed to be associated with the accumulation of NT-proBNP, and the mechanisms remain unknown (19, 20). Use of thalidomide has been shown to improve response to therapy, but also increases the risk of toxicities (21). The choice of treatment course also shows that 7–9 courses of treatment is most suitable, and more courses of treatment will not bring additional benefits of cardiac response, which also verifies the current view of AL with a fixed course of treatment. Deep relief of hematology brings about benefits in terms of cardiac relief. We integrated a variety of treatment-related factors to predict the success rate of cardiac treatment. This result has important clinical application value.
From our model, different treatment options do not bring about statistically significant OS outcomes, which also suggests that for AL patients with cardiac involvement, more effective treatment options still need to be explored. Daratumumab is an anti-CD38 monoclonal antibody used for multiple myeloma. Daratumumab has been reported to have good effects on AL (22–24). Among patients with newly diagnosed AL, the addition of daratumumab to bortezomib, cyclophosphamide, and dexamethasone was associated with higher rates of hematologic response than a control group (53.3% vs. 18.1%). At 6 months, more cardiac responses occurred in the daratumumab group than in the control group (41.5% vs. 22.2%) (25). In addition, ANDROMEDA study showed Daratumumab subcutaneous(DARA SC)-CyBorD was well tolerated in patients of AL amyloidosis, and importantly, with high rates of deep and durable hematologic and organ responses (24). However, Daratumumab is new antibody drug and the sample size treated with Daratumumab in our database is insufficient to be included in the nomogram model. With the increase in data, our nanogram will continue to incorporate new treatment options for adjustments in the future to better predict clinical outcomes.
In conclusion, based on the largest Chinese cohort of patients with cardiac AL, we identified a nomogram model to predict cardiac response and OS in AL patients with cardiac involvement. This model is more individualized and dynamic, it can accurately predict the clinical outcomes, and it has important clinical application value.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
YL, HW, and JL designed and organized the study. YL and HW wrote the first draft of the manuscript. MZ and YC analyzed the data statistically. WYan, AL, and WYang revised the manuscript critically. All authors contributed to the article and approved the final version of the manuscript.
Funding
This study received technical support from Neusoft Co., Ltd. Neusoft Co., Ltd. was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication. There was no external funding for this study.
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.
This study received technical support from Neusoft Co., Ltd. It had the following involvement with the study: the data was processed and analyzed on their automatic medical analysis platform.
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.
Acknowledgments
The authors would like to thank the organization of Chinese Registration Network for Light-chain Amyloidosis (CRENLA) for providing the database. This work was technically supported by the Neusoft Co., Ltd. We are grateful for their support in terms of data extraction, analysis, and evaluation by leveraging RealMedSci and for allowing us to use their automatic medical analysis platform.
References
1. Rysava R. AL Amyloidosis: Advances in Diagnostics and Treatment. Nephrol Dial Transplant (2019) 34(9):1460–6. doi: 10.1093/ndt/gfy291
2. Milani P, Merlini G, Palladini G. Light Chain Amyloidosis. Mediterr J Hematol Infect Dis (2018) 10(1):e2018022. doi: 10.4084/MJHID.2018.022
3. Tuzovic M, Yang EH, Baas AS, Depasquale EC, Deng MC, Cruz D, et al. Cardiac Amyloidosis: Diagnosis and Treatment Strategies. Curr Oncol Rep (2017) 19(7):46. doi: 10.1007/s11912-017-0607-4
4. Rubin J, Maurer MS. Cardiac Amyloidosis: Overlooked, Underappreciated, and Treatable. Annu Rev Med (2020) 71:203–19. doi: 10.1146/annurev-med-052918-020140
5. Gertz MA, Dispenzieri A. Systemic Amyloidosis Recognition, Prognosis, and Therapy: A Systematic Review. JAMA (2020) 324(1):79–89. doi: 10.1001/jama.2020.5493
6. Stern LK, Kittleson MM. Updates in Cardiac Amyloidosis Diagnosis and Treatment. Curr Oncol Rep (2021) 23(4):47. doi: 10.1007/s11912-021-01028-8
7. Palladini G, Dispenzieri A, Gertz MA, Kumar S, Wechalekar A, Hawkins PN, et al. New Criteria for Response to Treatment in Immunoglobulin Light Chain Amyloidosis Based on Free Light Chain Measurement and Cardiac Biomarkers: Impact on Survival Outcomes. J Clin Oncol (2012) 30(36):4541–9. doi: 10.1200/JCO.2011.37.7614
8. Caulfield S, Menezes G, Marignol L, Poole C. Nomograms Are Key Decision-Making Tools in Prostate Cancer Radiation Therapy. Urol Oncol (2018) 36(6):283–92. doi: 10.1016/j.urolonc.2018.03.017
9. Li W, Zhao K, Wang Z. Prognostic Nomograms Based on Immune Scores for Head-Neck Squamous Cell Carcinoma Patients. Eur Arch Otorhinolaryngol (2021) 278(7):2493–500. doi: 10.1007/s00405-020-06358-0
10. Nieder C, Mehta MP, Geinitz H, Grosu AL. Prognostic and Predictive Factors in Patients With Brain Metastases From Solid Tumors: A Review of Published Nomograms. Crit Rev Oncol Hematol (2018) 126:13–8. doi: 10.1016/j.critrevonc.2018.03.018
11. Dispenzieri A, Gertz MA, Kyle RA, Lacy MQ, Burritt MF, Therneau TM, et al. Serum Cardiac Troponins and N-Terminal Pro-Brain Natriuretic Peptide: A Staging System for Primary Systemic Amyloidosis. J Clin Oncol (2004) 22(18):3751–7. doi: 10.1200/JCO.2004.03.029
12. Kumar S, Dispenzieri A, Lacy MQ, Hayman SR, Buadi FK, Colby C, et al. Revised Prognostic Staging System for Light Chain Amyloidosis Incorporating Cardiac Biomarkers and Serum Free Light Chain Measurements. J Clin Oncol (2012) 30(9):989–95. doi: 10.1200/JCO.2011.38.5724
13. Comenzo RL, Reece D, Palladini G, Seldin D, Sanchorawala V, Landau H, et al. Consensus Guidelines for the Conduct and Reporting of Clinical Trials in Systemic Light-Chain Amyloidosis. Leukemia (2012) 26(11):2317–25. doi: 10.1038/leu.2012.100
14. Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in Oncology: More Than Meets the Eye. Lancet Oncol (2015) 16(4):e173–80. doi: 10.1016/S1470-2045(14)71116-7
15. Zhang Y, Chen XL, Chen WM, Zhou HB. Prognostic Nomogram for the Overall Survival of Patients With Newly Diagnosed Multiple Myeloma. BioMed Res Int (2019) 2019:5652935. doi: 10.1155/2019/5652935
16. Wechalekar AD, Gillmore JD, Bird J, Cavenagh J, Hawkins S, Kazmi M, et al. Guidelines on the Management of AL Amyloidosis. Br J Haematol (2015) 168(2):186–206. doi: 10.1111/bjh.13155
17. Palladini G, Sachchithanantham S, Milani P, Gillmore J, Foli A, Lachmann H, et al. A European Collaborative Study of Cyclophosphamide, Bortezomib, and Dexamethasone in Upfront Treatment of Systemic AL Amyloidosis. Blood (2015) 126(5):612–5. doi: 10.1182/blood-2015-01-620302
18. Dinner S, Witteles W, Afghahi A, Witteles R, Arai S, Lafayette R, et al. Lenalidomide, Melphalan and Dexamethasone in a Population of Patients With Immunoglobulin Light Chain Amyloidosis With High Rates of Advanced Cardiac Involvement. Haematologica (2013) 98(10):1593–9. doi: 10.3324/haematol.2013.084574
19. Dispenzieri A, Dingli D, Kumar SK, Rajkumar SV, Lacy MQ, Hayman S, et al. Discordance Between Serum Cardiac Biomarker and Immunoglobulin-Free Light-Chain Response in Patients With Immunoglobulin Light-Chain Amyloidosis Treated With Immune Modulatory Drugs. Am J Hematol (2010) 85(10):757–9. doi: 10.1002/ajh.21822
20. Tapan U, Seldin DC, Finn KT, Fennessey S, Shelton A, Zeldis JB, et al. Increases in B-Type Natriuretic Peptide (BNP) During Treatment With Lenalidomide in AL Amyloidosis. Blood (2010) 116(23):5071–2. doi: 10.1182/blood-2010-09-305136
21. Palladini G, Russo P, Lavatelli F, Nuvolone M, Albertini R, Bosoni T, et al. Treatment of Patients With Advanced Cardiac AL Amyloidosis With Oral Melphalan, Dexamethasone, and Thalidomide. Ann Hematol (2009) 88(4):347–50. doi: 10.1007/s00277-008-0600-y
22. Sher T, Fenton B, Akhtar A, Gertz MA. First Report of Safety and Efficacy of Daratumumab in 2 Cases of Advanced Immunoglobulin Light Chain Amyloidosis. Blood (2016) 128(15):1987–9. doi: 10.1182/blood-2016-06-722496
23. Abeykoon JP, Zanwar S, Dispenzieri A, Gertz MA, Leung N, Kourelis T, et al. Daratumumab-Based Therapy in Patients With Heavily-Pretreated AL Amyloidosis. Leukemia (2019) 33(2):531–6. doi: 10.1038/s41375-018-0262-2
24. Palladini G, Kastritis E, Maurer MS, Zonder J, Minnema MC, Wechalekar AD, et al. Daratumumab Plus CyBorD for Patients With Newly Diagnosed AL Amyloidosis: Safety Run-in Results of ANDROMEDA. Blood (2020) 136(1):71–80. doi: 10.1182/blood.2019004460
Keywords: light-chain amyloidosis, cardiac involvement, nomogram model, cardiac response, overall survival
Citation: Li Y, Cao Y, Zheng M, Hu J, Yan W, Liu X, Liao A, Yang W, Li J and Wang H (2021) Nomogram Model for Dynamic and Individual Prediction of Cardiac Response and Survival for Light Chain Amyloidosis in 737 Patients With Cardiac Involvement. Front. Oncol. 11:758502. doi: 10.3389/fonc.2021.758502
Received: 14 August 2021; Accepted: 22 November 2021;
Published: 09 December 2021.
Edited by:
María-Victoria Mateos, University Hospital of Salamanca, SpainReviewed by:
Jian Hou, Shanghai JiaoTong University, ChinaMichael Rosenzweig, City of Hope National Medical Center, United States
Copyright © 2021 Li, Cao, Zheng, Hu, Yan, Liu, Liao, Yang, Li and Wang. 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: Huihan Wang, MTM5MTg2MDc2NEBxcS5jb20=