- 1The Second Affiliated Hospital, Jiangxi Medical College, Nanchang University, Nanchang, China
- 2Department of Urology, Second Affiliated Hospital of Nanchang University, Nanchang, China
Background and purpose: Distant metastasis in bladder cancer is linked to poor prognosis and significant mortality. Machine learning (ML), a key area of artificial intelligence, has shown promise in the diagnosis, staging, and treatment of bladder cancer. This study aimed to employ various ML techniques to predict distant metastasis in patients with bladder cancer.
Patients and methods: Patients diagnosed with bladder cancer in the Surveillance, Epidemiology, and End Results (SEER) database from 2000 to 2021 were included in this study. After a rigorous screening process, a total of 4,108 patients were selected for further analysis, divided in a 7:3 ratio into a training cohort and an internal validation cohort. In addition, 118 patients treated at the Second Affiliated Hospital of Nanchang University were included as an external validation cohort. Features were filtered using the least absolute shrinkage and selection operator (LASSO) regression algorithm. Based on the significant features identified, three ML algorithms were utilized to develop prediction models: logistic regression, support vector machine (SVM), and linear discriminant analysis (LDA). The predictive performance of the three models was evaluated by obtaining the area under the receiver operating characteristic (ROC) curve (AUC), the precision, the accuracy, and the F1 score.
Results: According to the statistical results, the final probability of distant metastasis in the population was 12.0% (n = 495). LASSO regression analysis revealed that age, chemotherapy, tumor size, the examination of non-regional lymph nodes, and regional lymph node evaluation were significantly associated with distant metastasis of bladder cancer. In the internal validation cohort, the prediction accuracy rates for logistic regression, SVM, and LDA were 0.874, 0.877, and 0.845, respectively. The precision rates were 0.805, 0.769, and 0.827, respectively, and the F1 scores were 0.821, 0.819, and 0.835, respectively. The ROC curve demonstrated that the AUC for all models was greater than 0.7. In the external validation cohort, the prediction accuracy rates for logistic regression, SVM, and LDA were 0.856, 0.848, and 0.797, respectively, with the ROC curve indicating that the AUC also exceeded 0.7. The precision rates were 0.877, 0.718, and 0.736, respectively, and the F1 scores were 0.797, 0.778, and 0.762, respectively. Among the algorithms used, logistic regression demonstrated better predictive efficiency than the other two methods. The top three variables with the highest importance scores in the logistic regression were non-regional lymph nodes, age, and chemotherapy.
Conclusion: The prediction model developed using three ML algorithms demonstrated strong accuracy and discriminative capability in predicting distant metastasis in patients with bladder cancer. This might help clinicians in understanding patient prognosis and in formulating personalized treatment strategies, ultimately improving the overall prognosis of patients with bladder cancer.
Introduction
Bladder cancer (Bca) is one of the most common malignant tumors of the urinary tract. As the ninth most prevalent malignant tumor globally, its incidence and prevalence are increasing year by year. Approximately 550,000 patients are newly diagnosed with Bca each year (1, 2). Men tend to have a higher incidence, and smoking is believed to be a contributing factor (3). Based on the ability of the tumor to invade, Bca can be classified into non-muscle-invasive bladder cancer (NMIBC), muscle-invasive bladder cancer (MIBC), and metastatic forms of the disease (4). Approximately 70% of newly diagnosed cases are NMIBC, while approximately 30% are diagnosed as MIBC, frequently with metastatic characteristics. Half of the patients with muscle-invasive disease will die from metastases within 2 years. Compared with the 5-year survival rate of 77% for all stages of Bca, the 5-year survival rate for metastatic Bca is only 5% (5–7). Therefore, it is critical to develop models for the prediction of distant metastasis of Bca.
As an essential branch of artificial intelligence (AI), machine learning (ML) develops predictive models by automatically learning from large datasets to improve prediction algorithms. This advancement assists clinicians in identifying high-risk patients and in evaluating the prognosis of various diseases (8, 9). Diagnostic and prognostic models built using ML algorithms based on pathomic data have demonstrated remarkable efficiency in distinguishing patients with Bca from those with glandular cystitis, as well as in predicting the survival outcomes of BCa (10). Network analysis methods rooted in ML frameworks can effectively identify biomarkers linked to immunotherapy responses, resulting in robust predictions for precision oncology (11). Furthermore, a ML model utilizing full-sequence MRI can accurately forecast the depth of invasion of Bca prior to surgery, helping clinicians in recognizing pathological features associated with tumor invasion ahead of invasive procedures (12). In addition, ML has been widely applied in prostate cancer and kidney cancer research, and it also holds significant potential in various benign conditions, such as urinary tract stones (12, 13).
In this study, data on patients with Bca and their clinical and pathological characteristics were retrieved from the Surveillance, Epidemiology, and End Results (SEER) database from 2000 to 2021. Our aim was to utilize ML algorithms to develop a reliable model for the prediction of distant metastasis of Bca, thereby providing clinical support for treatment decisions and individualized prognosis assessment.
Patients and methods
Patients
Detailed data of patients with Bca from 17 registration centers between 2000 and 2021 were collected from SEER*Stat (version 8.4.3), including the demographic, clinical, and pathological characteristics. The inclusion criterion comprised patients diagnosed with Bca during this time frame. The exclusion criteria were cases with unknown race, unknown pathological grade, unknown radiotherapy status, unknown tumor size or blank records, unknown non-regional lymph nodes, and unknown T, N, or M stage, as well as unknown regional lymph node examination results. Figure 1 illustrates the complete screening process. As SEER is a public database and all records are de-identified, no additional ethical approval was required.
Figure 1. Surveillance, Epidemiology, and End Results (SEER) database screening and research flowchart.
The external validation cohort comprised 118 eligible patients diagnosed with Bca between 2016 and 2020. Upon review and analysis, a total of 18 patients were found to have ultimately developed distant metastasis. The final follow-up was scheduled for December 2023. This study was approved by the Ethics Committee of the Second Affiliated Hospital of Nanchang University, with a waiver of informed consent.
Screening for clinical and pathological features
The variables selected included demographic characteristics (i.e., age, sex, and race), tumor characteristics (i.e., tumor size, grade, histology, T stage, and N stage), treatment information (e.g., chemotherapy and radiation therapy), and other variables (e.g., regional lymph nodes examined and non-regional lymph nodes). The primary endpoint of this study was the development of distant metastases in Bca. For better analysis, some of the variables were processed in the SEER database. A number of continuous variables, i.e., age, tumor size, regional lymph nodes examined, and non-regional lymph nodes, were transformed into categorical variables. All clinical and pathological features were used in the training set to screen the features most related to the distant metastasis of Bca using least absolute shrinkage and selection operator (LASSO) regression.
Variables were selected not only for their individual significance but also for their potential interactions and contributions to the overall predictive model. The LASSO regression technique applies a penalty proportional to the absolute value of the coefficient sizes, effectively shrinking some coefficients to zero. This process facilitates variable selection while preventing overfitting, which is particularly important when dealing with high-dimensional datasets. The tuning parameter (λ) regulates the strength of the penalty: a larger λ leads to more coefficients being set to zero, thereby highlighting the most important predictors.
The clinical rationale for selecting these variables is grounded in their established correlations with outcomes in previous studies and clinical guidelines. For instance, age and gender have been shown to significantly influence the prognosis of cancer and the response to treatment. Tumor characteristics such as size, grade, and stage (T and N) are critical in the assessment of the aggressiveness of a cancer and the likelihood of metastasis. In addition, treatment variables, particularly the effects of chemotherapy and radiation therapy on cancer progression, have been well documented and are essential for our analysis. Lastly, examination of regional and non-regional lymph nodes is crucial for the assessment of the metastatic spread, thus helping in the development of predictive models for distant metastasis.
Construction of the machine learning model
As a common linear classifier, logistic regression establishes a regression formula for the classification boundary based on existing data, enabling the accurate identification of binary classification problems. It primarily implements the algorithm through the logistic function, specifically the sigmoid function. Support vector machine (SVM), as a binary classification model, is fundamentally defined as a linear classifier that maximizes the margin in the feature space. By maximizing this margin, SVM can ultimately be transformed into a solution for a convex quadratic programming problem. Linear discriminant analysis (LDA), a classic method in statistics and ML, aims to maximize the separation between samples of different categories while minimizing the scatter within the same category through a projection method, thus achieving linear classification of data.
The clinical and pathological characteristics identified above were utilized in constructing the ML models. Three types of ML algorithms were trained using a binary classification approach. The training set (70%) was used to train the model, while the validation set (30%) was employed to assess the predictive performance of the model, followed by external validation. A weighted average of the precision, accuracy, and F1 score was employed to minimize the impact of sample imbalance on the evaluation results. The predictive values of logistic regression, SVM, and LDA were compared by calculating the prediction accuracy, precision, F1 score, and the area under the receiver operating characteristic (ROC) curve (AUC).
Statistical analysis
Data analysis was conducted using SPSS 27.0, R language software (version 4.3.1; http://www.r-project.org/), and Python language software (https://www.python.org/downloads/release/python-380/). For the analysis and processing, all categorical variables were expressed as numbers with percentages, and LASSO regression was employed to identify key predictors. LASSO regression was carried out using the glmnet package in R or the scikit-learn package in Python, and the ROC curve was plotted using the matplotlib package in Python.
Results
Patient characteristics and metastasis
The characteristics of the training cohort, the internal validation cohort, and the external validation cohort of patients are presented in Table 1. The probabilities of distant metastasis for Bca in the three cohorts were 11.9% (n = 343), 12.3% (n = 152), and 19.4% (n = 14), respectively.
Identification of significant features
All clinical and pathological characteristics were analyzed using LASSO regression. The results indicated that age, chemotherapy, tumor size, non-regional lymph nodes, and regional lymph node examination were the significant variables associated with the distant metastasis of Bca, as detailed in Figure 2.
Figure 2. (A) Least absolute shrinkage and selection operator (LASSO) coefficient distribution of the clinical and pathological characteristics. (B) LASSO regression analysis was performed using the minimum criterion and a 10-fold cross-validation method. By introducing a penalty adjustment parameter (λ), the coefficients of the features are compressed towards zero, enabling automatic feature selection.
Construction and evaluation of machine models
Three ML models were developed using the significant variables identified earlier. A weighted average of the precision, accuracy, and the F1 score was utilized to minimize the effects of sample imbalance on the evaluation results. In the internal validation cohort, the prediction accuracy rates for logistic regression, SVM, and LDA were 0.874, 0.877, and 0.845, respectively. The accuracy rates were 0.805, 0.769, and 0.827, respectively, and the F1 scores were 0.821, 0.819, and 0.835, respectively. The AUC (95% CI) values were 0.701 (0.634–0.767), 0.701 (0.641–0.766), and 0.701 (0.635–0.780), respectively. In the external validation cohort, the prediction accuracy rates for logistic regression, SVM, and LDA were 0.856, 0.848, and 0.797, respectively. The precision rates were 0.877, 0.718, and 0.736, respectively, and the F1 scores were 0.797, 0.778, and 0.762, respectively. The AUC (95% CI) values were 0.740 (0.682–0.795), 0.740 (0.674–0.792), and 0.740 (0.685–0.796), respectively. Details can be found in Figure 3. These results indicate that all three models can effectively predict distant metastasis in Bca. However, after comprehensive comparison, logistic regression was found to have the best predictive value and identification ability. The top three variables with logistic regression importance scores were non-regional lymph nodes, age, and chemotherapy (Figure 4).
Figure 3. Receiver operating characteristic (ROC) curves of the internal validation cohort (A) and the external validation cohort (B).
Discussion
Distant metastasis of Bca significantly impacts the survival rates and treatment prospects of patients. The common sites of metastasis include the liver, kidneys, and bones. Patients with Bca face an increased risk of death due to distant metastasis (14, 15). Those with localized Bca often experience an excellent quality of life and progression-free survival with aggressive treatment. Despite advancements in various targeted therapies and immunotherapies, distant metastasis remains highly lethal when it occurs (16–18). The metastatic spread in patients with Bca is believed to be random and unpredictable. Although early systemic treatment for Bca has made progress and is somewhat effective in preventing tumor progression, some individuals still inevitably develop recurrent or metastatic disease (19, 20). Therefore, it is crucial to develop new models for predicting the metastasis of Bca.
ML, a subfield of AI, emerged from the intersection of statistics and computer science. Among the algorithms that demonstrate superior performance, LDA and SVM are the most commonly recognized (21, 22). The continuous advancement of computing systems and recognition software has contributed to the popularity of ML-based systems. These systems are capable of performing complex tasks in bioinformatics and medical imaging, assisting healthcare professionals in improving the diagnostic accuracy, predicting disease progression, and even exploring new therapeutic avenues (23–25). A retrospective study by Denget al. showed that a ML model combining radiomic features and clinical variables could accurately predict the pathological grade of Bca (26). In addition, the quantitative assessment of tumor-infiltrating lymphocytes using ML mitigates the inter- and intra-observer variability, effectively predicting survival in MIBC and identifying patients who may benefit from immunotherapy or treatment adjustments (27). Liu et al. predicted lung metastasis in Bca by developing a nomogram model; however, this model did not include all distant metastases and highlighted the need for improved algorithmic approaches (28). In the study by Chen et al., multiple clinical variables were utilized to develop a nomogram for predicting distant metastasis in patients with urothelial Bca (29). In this research, a more advanced ML algorithm was employed, incorporating AI methods to successfully predict distant metastasis in Bca. ML utilizes advanced algorithms to analyze various clinical variables and integrate multiple data sources, resulting in a more comprehensive assessment of the condition of a patient. This approach improves the diagnostic accuracy and reduces the likelihood of misdiagnosis (9).
Our research utilized three ML models: logistic regression, SVM, and LDA. After screening the clinical variables and pathological characteristics using LASSO regression, we incorporated the identified significant variables and features into the construction of our ML models. In the internal validation cohort, we evaluated the predictive capability of the three models by calculating the prediction accuracy, precision, F1 score, and the AUC. The results indicated that all three models effectively predicted distant metastasis in Bca. This was further confirmed with the external validation cohort. By utilizing the model described above, the stage and the pathological type of patients with Bca can be predicted in real time, allowing for targeted treatment and follow-up for those at high risk of distant metastasis. However, this study has several limitations. Firstly, as a retrospective analysis based on the SEER database, it may have lacked certain potentially important factors, such as vascular invasion and preoperative laboratory results. In addition, patients with unknown variable information were excluded, which means that the results were inevitably impacted by selection bias. Secondly, variations in the sequencing of radiotherapy and chemotherapy, the types of radiation used, and the chemotherapy agents administered may have affected the likelihood of Bca metastasizing to distant sites. Finally, while we have validated our model using an external cohort, it is essential to include additional clinical and molecular biological features in the future to enhance its predictive power in complex clinical settings. Furthermore, multicenter prospective studies are essential for the development of more robust algorithms, such as deep learning and ensemble models. These studies should also focus on comparing model performance across different groups (such as gender and age) in multicenter external validations to further verify our findings.
Conclusion
In summary, we successfully developed and validated three ML models for the prediction of distant metastasis in Bca. In particular, the logistic regression algorithm achieved better prediction results. The application of these models can assist clinicians in making more accurate treatment decisions and in providing personalized follow-up management, ultimately improving patient prognosis.
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 humans were approved by Biomedical Research Ethics Committee of the Second Affiliated Hospital of Nanchang University. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because it’s a retrospective study, and is difficult to obtain informed consent, and no loss to the patient.
Author contributions
XZ: Data curation, Formal Analysis, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing. XR: Supervision, Writing – review & editing. JH: Writing – review & editing. JZ: Supervision, Validation, Writing – review & editing. TZ: Funding acquisition, Supervision, Writing – review & editing. HC: Methodology, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by the National Natural Science Foundation of China (no. 82260598) and the Jiangxi Provincial Academic and Technical Leader Training Program in Major Disciplines (no. 20225BCJ22009).
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.
References
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2021) 71(3):209–49. doi: 10.3322/caac.21660
2. Richters A, Aben KKH, Kiemeney LALM. The global burden of urinary bladder cancer: an update. World J Urol. (2020) 38(8):1895–904. doi: 10.1007/s00345-019-02984-4
3. Lobo N, Afferi L, Moschini M, Mostafid H, Porten S, Psutka SP, et al. Epidemiology, screening, and prevention of bladder cancer. Eur Urol Oncol. (2022) 5(6):628–39. doi: 10.1016/j.euo.2022.10.003
4. Compérat E, Amin MB, Cathomas R, Choudhury A, De Santis M, Kamat A, et al. Current best practice for bladder cancer:a narrative review of diagnostics and treatments. Lancet. (2022) 400(10364):1712–21. doi: 10.1016/S0140-6736(22)01188-6
5. Patel VG, Oh WK, Galsky MD. Treatment of muscle-invasive and advanced bladder cancer in 2020. CA Cancer J Clin. (2020) 70(5):404–23. doi: 10.3322/caac.21631
6. Ghatalia P, Kaur J, Sonpavde G. Muscle invasive bladder cancer: where is the field headed? Expert Opin Biol Ther. (2023) 23(9):913–27. doi: 10.1080/14712598.2023.2238607
7. Dobruch J, Oszczudłowski M. Bladder cancer : current challenges and future directions. Medicina(Kaunas). (2021) 57(8):749. doi: 10.3390/medicina57080749
8. Moore NS, McWilliam A, Aneja S. Bladder cancer radiation oncology of the future : prognostic modelling, radiomics, and treatment planning with artificial intelligence. Semin Radiat Oncol. (2023) 33(1):70–5. doi: 10.1016/j.semradonc.2022.10.009
9. Borhani S, Borhani R, Kajdacsy-Balla A. Artificial intelligence: A promising frontier in bladder cancer diagnosis and outcome prediction. Crit Rev Oncol Hematol. (2022) 171:103601. doi: 10.1016/j.critrevonc.2022.103601
10. Chen S, Jiang L, Zheng X, Shao J, Wang T, Zhang E, et al. Clinical use of machine learning-based pathomics signature for diagnosis and survival prediction of bladder cancer. Cancer Sci. (2021) 112(7):2905–14. doi: 10.1111/cas.14927
11. Kong J, Ha D, Lee J, Kim I, Park M, Im SH, et al. Network-based machine learning approach to predict immunotherapy response in cancer patients. Nat Commun. (2022) 13(1):3703. doi: 10.1038/s41467-022-31535-6
12. Hamad J, McCloskey H, Milowsky MI, Royce T, Smith A. Bladder preservation in muscle-invasive bladder cancer: a comprehensive review. Int Braz J Urol. (2020) 46(2):169–84. doi: 10.1590/S1677-5538.IBJU.2020.99.01
13. Chen H, Yang W, Ji Z. Machine learning-based identification of tumor-infiltrating immune cell-associated model with appealing implications in improving prognosis and immunotherapy response in bladder cancer patients. Front Immunol. (2023) 14:1171420. doi: 10.3389/fimmu.2023.1171420
14. Alfred Witjes J, Max Bruins H, Carrión A, Cathomas R, Compérat E, Efstathiou JA, et al. European association of urology guidelines on muscle-invasive and metastatic bladder cancer : summary of the 2023 guidelines. Eur Urol. (2024) 85(1):17–31. doi: 10.1016/j.eururo.2023.08.016
15. Safiri S, Kolahi AA, Naghavi M, Meraz-Munoz A, North SA, Alimohamed NS, et al. Global, regional and national burden of bladder cancer and its attributable risk factors in 204 countries and territories,1990-2019:a systematic analysis for the Global Burden of Disease study 2019. BMJ Glob Health. (2021) 6(11):e004128. doi: 10.1136/bmjgh-2020-004128
16. Lopez-Beltran A, Cookson MS, Guercio BJ, Cheng L. Advances in diagnosis and treatment of bladder cancer. BMJ. (2024) 384:e076743. doi: 10.1136/bmj-2023-076743
17. Jiang DM, Gupta S, Kitchlu A, Meraz-Munoz A, North SA, Alimohamed NS, et al. Defining cisplatin eligibility in patients with muscle-invasive bladder cancer. Nat Rev Urol. (2021) 18(2):104–14. doi: 10.1038/s41585-020-00404-6
18. Kwon WA, Seo HK. Emerging agents for the treatment of metastatic urothelial cancer. Investig Clin Urol. (2021) 62(3):243–55. doi: 10.4111/icu.20200597
19. Charpentier M, Gutierrez C, Guillaudeux T, Verhoest G, Pedeux . Noninvasive urine-based tests to diagnose or detect recurrence of bladder cancer. Cancers (Basel). (2021) 13(7):1650. doi: 10.3390/cancers13071650
20. Ruan W, Chen X, Huang M, Wang H, Chen J, Liang Z, et al. A urine-based DNA methylation assay to facilitate early detection and risk stratification of bladder cancer. Clin Epigenet. (2021) 13(1):91. doi: 10.1186/s13148-021-01073-x
21. Ferro M, Falagario UG, Barone B, Maggi M, Crocetto F, Busetto GM, et al. Artificial intelligence in the advanced diagnosis of bladder cancer-comprehensive literature review and future advancement. Diagnostics (Basel). (2023) 13(13):2308. doi: 10.3390/diagnostics13132308
22. Barberis E, Khoso S, Sica A, Falasca M, Gennari A, Dondero F, et al. Precision medicine approaches with metabolomics and artificial intelligence. Int J Mol Sci. (2022) 23(19):11269. doi: 10.3390/ijms231911269
23. Benning L, Peintner A, Peintner L. Advances in and the applicability of machine learning-based screening and early detection approaches for cancer: A primer. Cancers (Basel). (2022) 14(3):623. doi: 10.3390/cancers14030623
24. Choi RY, Coyner AS, Kalpathy-Cramer J, Chiang MF, Campbell JP. Introduction to machine learning, neural networks, and deep learning. Transl Vis Sci Technol. (2020) 9(2):14. doi: 10.1167/tvst.9.2.14
25. Castiglioni I, Rundo L, Codari M, Di Leo G, Salvatore C, Interlenghi M, et al. AI applications to medical images: From machine learning to deep learning. Phys Med. (2021) 83:9–24. doi: 10.1016/j.ejmp.2021.02.006
26. Deng Z, Dong W, Xiong S, Jin D, Zhou H, Zhang L, et al. Machine learning models combining computed tomography semantic features and selected clinical variables for accurate prediction of the pathological grade of bladder cancer. Front Oncol. (2023) 13:1166245. doi: 10.3389/fonc.2023.1166245
27. Zheng Q, Yang R, Ni X, Yang S, Jiao P, Wu J, et al. Quantitative assessment of tumor-infiltrating lymphocytes using machine learning predicts survival in muscle-invasive bladder cancer. J Clin Med. (2022) 11(23):7081. doi: 10.3390/jcm11237081
28. Liu L, Sun FZ, Zhang PY, Xiao Y, Ni HX. Development and validation a model for predicting overall survival of bladder cancer with lung metastasis: a population-based study. Eur J Med Res. (2023) 28(1):279. doi: 10.1186/s40001-023-01261-w
Keywords: machine learning, bladder cancer, SEER database, distant metastasis, predictive value
Citation: Zou XC, Rao XP, Huang JB, Zhou J, Chao HC and Zeng T (2024) Predicting distant metastasis of bladder cancer using multiple machine learning models: a study based on the SEER database with external validation. Front. Oncol. 14:1477166. doi: 10.3389/fonc.2024.1477166
Received: 07 August 2024; Accepted: 19 November 2024;
Published: 13 December 2024.
Edited by:
Jacopo Troisi, University of Salerno, ItalyReviewed by:
Bernhard Ralla, Charité University Medicine Berlin, GermanyYanyang Jin, First Affiliated Hospital of Jinzhou Medical University, China
Copyright © 2024 Zou, Rao, Huang, Zhou, Chao and Zeng. 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: Tao Zeng, dGFvemVuZzQwNzA5QHNpbmEuY29t
†These authors have contributed equally to this work