- 1Department of Radiation Oncology, Fudan University Shanghai Cancer Center, Shanghai, China
- 2Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China
- 3Shanghai Key Laboratory of Radiation Oncology, Shanghai, China
Purpose: This study aims to demonstrate the feasibility of clinical implementation of automated treatment planning (ATP) using voxel-based dose prediction and post-optimization strategies for rectal cancer on uRT (United Imaging Healthcare, Shanghai, China) treatment planning system.
Methods: A total of 180 previously treated rectal cancer cases were enrolled in this study, including 160 cases for training, 10 for validation and 10 for testing. Using CT image data, planning target volumes (PTVs) and contour delineation of the organs at risk (OARs) as input and three-dimensional (3D) dose distribution as output, a 3D-Uet DL model was developed. Based on the voxel-wise prediction dose distribution, intensity-modulated radiation therapy (IMRT) plans were then generated automatedly using post-optimization strategies, including a complex clinical dose target metrics homogeneity index (HI) and conformation index (CI). To evaluate the performance of the proposed ATP approach, the dose-volume histogram (DVH) parameters of OARs and PTV and the 3D dose distributions of the plan were compared with those of manual plans.
Results: By combining clinical post-optimization strategies, the automatically generated treatment plan can achieve better homogeneous PTV coverage and dose sparing for OARs except the mean dose for femoral-head compared with the use of the mean square error objective function alone. Compared with the manual plan, no statistically significant differences in HI, CI or global maximum dose were found. The manual plans perform slightly better than plans with post-optimization strategies in other dosimetric indexes, but these plans are still within clinical requirements.
Conclusions: With the help of clinical post-optimization strategies, the proposed new ATP solution can generate IMRT plans that are within clinically acceptable levels and comparable to plans manually generated by dosimetrists.
Introduction
Intensity-modulated radiotherapy (IMRT) is an important treatment modality and has been widely used for many types of cancer (1). This radiotherapy technology has the advantage of achieving higher dose coverage of planning target volumes (PTVs) with steep dose gradients at the transition to orangs at risk (OARs). The design of the IMRT plans requires several optimization cycles with manual adjustments of the weights and the addition of new optimization objectives to meet the clinically specific criteria. It is a time-consuming and labor-intensive process and extremely dependent on the experience of the designer. Consequently, the quality of a plan depends drastically on planners (experience, skill and available time) (2). This variability may lead to suboptimal plans, resulting in a negative impact on tumor control in patients (3, 4).
Recently, deep learning (DL) has been used to automate treatment planning and has received considerable attention in the radiotherapy community (5, 6). An approach to reducing user variability and improving the quality and efficiency of treatment plans is the use of so-called knowledge-based planning (KBP). KBP has been implemented in the commercial treatment planning system (TPS) Eclipse as the RapidPlan module (Varian Medical Systems, Palo Alto, CA). This strategy uses a large number of previous acceptable or superior clinical patient databases to estimate specific dose metrics or dose-volume histograms (DVHs) for a new patient (7–11). The main limitation of this method is the lack of spatial information, and planners still need to use hand-crafted features from statistical analysis (9, 12–19). To solve this problem, DL-based 3D dose distribution prediction techniques for automated treatment planning (ATP) have become a major focus of research. Based on popular convolutional neural networks (CNNs), the patient-specific three-dimensional (3D) dose distribution can be achieved and later used as an objective to generate a treatment plan automatically. Accurate spatial dose distribution prediction can eliminate dependence on handcrafted features completely and potentially improve plan quality and consistency (20–22).
However, the clinical implementation of current DL-based ATP solutions remains stymied by two points. First, the performance of the DL method strongly depends on the database used for training. The optimal result of the automated generated plan is to make each voxel dose the same as the predicted dose. Thus, the current voxel-based dose prediction methods cannot obtain a better plan, or the plan quality is limited by the predicted dose results. Second, given the high safety requirements of medical applications, generating executable clinical automated treatment plan based on 3D dose distribution in closed commercial software architecture remains a challenge. To the best of our knowledge, previous studies have generally focused on algorithms to improve the accuracy of 3D dose distribution prediction (21–23) or have used the open-source toolkit matRad to generate radiation treatment planning for educational purposes and research (20, 24).
Consequently, a DL model was developed for predicting a 3D dose distribution in this study. To decrease the quality limits of the training data, clinical post-optimization strategies were used in the process of automated plan generation based on the 3D prediction dose. We extensively collaborated with United Imaging Healthcare (UIH) Co., Ltd. (Shanghai, China) to achieve the clinical implementation of ATP in its TPS uRT.
Materials and Methods
Patient Data Collection
A total of 180 rectal cancer patients undergoing radiotherapy between 2017 and 2019 in our cancer center were enrolled in this study. For analysis, data were randomly divided into 160 training sets, 10 validation sets and 10 testing sets. Simulation CT images (slice thickness 5 mm; 512×512 matrix) were acquired using a Philips Brilliance Big Bore multidetector-row spiral CT scanner (Philips Healthcare, Cleveland, OH). Radiation oncologists delineated the gross tumor volume (GTV), clinical target volume (CTV) and OARs in the planning CT. The prescription of each patient was 50 Gy in 25 fractions (2 Gy per fraction). The IMRT treatment plans were calculated and optimized in Pinnacle 8.0-9.10 TPS (Philips Radiation Oncology Systems, Fitchburg, WI, USA). All these treatment plans were generated with 9 equiangular 6-MV photon beams.
Model Architecture
A 3D U-net was developed in this study and is illustrated in Figure 1. The number of filters for each convolutional layer was 32, 64, 128, 256, and 512, and the feature map size was reduced by half after the max-pooling layer. All convolutional layers applied a 3×3×3 kernel except the output layer with a 1×1×1 kernel. The input data were CT images and contours of regions of interest (ROIs), including the PTV, body bladder, left femoral head and right femoral head. The output data were the 3D predicted dose distribution. The model was implemented in Keras, and the loss function used in the training process was the mean square error. The Adam optimization algorithm was used to minimize the loss function value between the predicted dose and the clinical truth. The network parameters of the mode were initialized using the He_normal initialization method (25).
Automated Planning
To provide a more complete picture, the overall description of the framework of this paper is displayed in Figure 2. First, using previously treated plans generated and optimized in Pinnacle TPS as training data, a DL model was developed for predicting a 3D dose distribution. Second, a new patient CT with contouring information of the target and OARs was input to automatically generate the 3D dose distribution of the current case by our DL model. Third, based on these 3D dose predictions, automated IMRT plans were then generated using two methods in UIH TPS. One approach uses the mean square error (MSE) optimization function only. Another combines MSE and clinical post-optimization strategies, which include the complex clinical dose target metrics homogeneity index (HI) and conformation index (CI). All plans were generated automatically using 9 equally spaced fixed coplanar 6-MV photon beams in UIH TPS.
The function of MSE is described as follows:
where V denotes all voxels, didenotes the dose of ith voxel, and dcalc and dpred represent the input and prediction dose in the iteration process, respectively.
The optimization function of HI is described as follows:
where WHI is the weight of the objective function, and DX is the corresponding dose of X% volume of the target in each iteration.
The optimization function of CI is described as follows:
where WCI is the weight of the objective function, TV is the volume of the target, VRI is the volume of the reference isodose, TVRI is the interaction volume of the reference isodose and TV, and Cindex is constant (typically fixed as 1).
The total normalized objective function is noted as follow:
where Ftotal denotes the total objective function; α and β are the weight of the function of MSE and the clinical optimization index, respectively; and D represents the different index.
Dosimetric Plan Evaluations
To evaluate the performance of the two proposed ATP methods, the DVH parameters of OARs and PTV and the 3D dose distributions of plan were compared between the ATP and original clinical plan. The clinically relevant dosimetric indexes, including the mean dose (Dmean), D2, D5, D95, and D98 for PTV (here, Di means the dose received by i% of PTV volume) and Dmean, V15, V25, V35, and V50 for OARs (here, Vi means volume fraction of OARs irradiated by i Gy), were calculated. HI (26) and CI (27) for PTV were further calculated using the following formulas:
where VPTV and VP are the volume of the PTV and the prescription dose region, respectively, and VR is the irradiated PTV volume of the prescription dose. These results were compared using a paired-samples t-test for the two models with p < 0.05 considered statistically significant.
Results
Isodose comparisons between the two plan-generated methods for a representative example case are shown in Figure 3. As evident in the graph, the MSE optimization plan (Figure 3C) delivered a substantially greater dose (110% of prescription with yellow line) to the PTV compared with the other methods. It is also worth noting that the plan (Figure 3B) using clinical post-optimization strategies resulted in approximately the same homogeneous PTV coverage as the manual plan (Figure 3A). To further illustrate the performance of the three methods. The OAR and PTV dose metrics are shown in Figures 4, 5 as violin plots.
Figure 3 The dose distribution comparison of PTV region (green shaded area) for two methods with manual plan: (A) manual plan, (B) dose prediction based post-optimization plan and (C) MSE optimization plan. The unit of color bar is Gy.
Figure 4 Violin plots representing the mean (white dot), σ (vertical rectangle), 95% percentile (black vertical line) and the probability distribution for the HI (A) and CI (B). The asterisks represent p ≤ 0.05. (A) Manual plan, (B) dose prediction-based post-optimization plan and (C) MSE optimization plan.
Figure 5 This graph displays a column scatter chart with the mean (±SE): (A) minimum, maximum, D98 and D2 of PTV; (B) mean value of PTV, bladder, femoral head-left and femoral head-right; (C) V15, V25, V35, and V50 for bladder; (D) violin plots of global maximum dose for three methods. The asterisks represent p ≤ 0.05. (A) Manual plan, (B) dose prediction-based post-optimization plan and (C) MSE optimization plan.
Figure 4 represents the violin plots for HI and CI. The performance of manual plans and dose prediction-based post-optimization plans were greater than those of MSE optimization plans. Significant differences were observed between MSE optimization plans and the other two approaches. For HI and CI, no significant differences were identified between manual plans and with clinical post-optimization strategies.
More clinically relevant dosimetric indexes and DVH are shown in Figures 5, 6, respectively. For dosimetric indexes of PTV, significant differences were observed between the MSE optimization plan and the other two methods. The MSE optimization plan did not outperform any dose metrics for PTV but had the lowest mean dose for femoral heads compared with the other two approaches. Of the specific dose metric constraints for the bladder, the manual plans achieved the lowest dose volume (V15, V25, V35 and V50) and mean dose followed by the dose prediction based post-optimization plans and then the MSE optimization plans. The global maximum dose doses of these four plans are presented in Figure 5D. Similarly, the MSE optimization approach was inferior to all methods and showed significant differences.
Figure 6 Mean DVHs of the manual plan (dotted line), MSE optimization plan (solid line), and dose prediction based post-optimization plan (dashed line).
Discussion
In this study, new optimization strategies directly related to clinical targets were used to improve the quality of ATPs. To our knowledge, this is the first study that implements clinical post-optimization strategies in the process of automated treatment planning for rectal cancer. We designed a deep learning model with clinical post-optimization strategies, as summarized in Figure 2. With the help of these intelligent strategies, the proposed optimization methods can pull high target dose and spare more OAR dose to obtain a better auto plan. This method can improve the plan quality of automated plans, which are limited by the quality of the training dataset.
There have been some studies on ATP techniques using DL neural networks (20–23, 28–30). Unlike conventional inverse optimization radiotherapy treatment planning with trial and error, ATP can be summarized into two steps: obtaining the predicted dose distribution and generating an automated executable plan (6). It was not easy to predict 3D voxel-wise dose distributions for IMRT plans in previous years due to the complicated relationship between OARs and PTV and the significant variability of PTV shapes. With the rapid advancement of machine learning, the accuracy of 3D dose distribution prediction methods has increased substantially. Compared with the DVH-based prediction algorithms, the prediction models have significant advantage in a way that it could provide spatial dose distribution information. These methods evolve from 2D model (28)to the 3D model (20, 29, 31, 32) and can provide a more robust predicted dose distribution.
However, these existing ATP strategies use the predicted 3D dose distribution as the end without any post-optimization. The optimal result of the automated generated plan is to make each voxel-wise dose the same as the predicted dose. Thus, current voxel-based methods cannot obtain a better plan due to two possible reasons. First, the performance of the deep learning-based model requires a sufficiently large, high-quality database. If the ground truth doses are suboptimal, the predicted doses will also be suboptimal as noted in the “garbage in, garbage out” paradigm. Second, more clinical conditions, such as significant variability in PTV size and complicated spatially neighboring anatomy, also have an impact on dose prediction. Previous studies require only the predicted dose distribution for optimization and did not take into consideration the diverse clinical scenarios for each patient. Therefore, it is challenging to achieve a specific plan for different patients, and the final plan had a worse homogeneity index of the target. Based on the above analysis, the proposed automated treatment planning strategy in our study serves as a step forward and provides a new idea to improve the performance of ATP.
The dose distribution in Figure 3 shows that the method using clinical optimization strategies exhibited more accurate homogeneous PTV coverage than the MSE optimization plan. This conclusion is also strongly supported by more clinical interested dosimetric indexes found on Figures 4–6. These approaches can generate treatment plans that are within clinically acceptable levels and noninferior to plans manually generated by dosimetrists. The quality and consistency of treatment planning for radiotherapy can be largely improved.
However, there are some limitations to address in this study. We investigated the applicability of the clinical optimization of ATP for rectal cancer. The feasibility of the automated plan strategy should be demonstrated for nasopharyngeal cancer patients who have a rather complicated relationship between OARs and PTV and multiple targets, thus requiring different dose prescriptions. These are important priorities of future research.
Conclusion
This paper demonstrates the feasibility of a new automated treatment planning strategy that includes clinical dose metrics post-optimization based on 3D dose distribution prediction. The approach can generate treatment plans that are within clinically acceptable levels and comparable to plans manually generated by dosimetrists. This methodology has great potential to improve the consistency and quality of IMRT planning by minimizing human intervention in the process of plan design.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
Conception, design, and drafting the manuscript were performed by YZ, LY, JW, and WH. Data collection and interpreting were performed by YF, YY, and ZW. All authors contributed to the article and approved the submitted version.
Funding
This work is supported the Shanghai Committee of Science and Technology Fund (19DZ1930902) and Xuhui District Artificial Intelligence Medical Hospital Cooperation Project (2020-009).
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.
Acknowledgments
The authors want to thank United Imaging Healthcare Co., Ltd. (Shanghai, China) for their excellent technical assistance and valuable suggestions.
References
1. Hussein M, Heijmen B, Verellen D, Nisbet A. Automation in Intensity Modulated Radiotherapy Treatment Planning-A Review of Recent Innovations. Br J Radiol (2018) 91(1092):20180270. doi: 10.1259/bjr.20180270
2. Nelms B, Robinson G, Markham J, Velasco K, Boyd S, Narayan S, et al. Variation in External Beam Treatment Plan Quality: An Inter-Institutional Study of Planners and Planning Systems. Pract Radiat Oncol (2012) 2(4):296–305. doi: 10.1016/j.prro.2011.11.012
3. González Ferreira J, Jaén Olasolo J, Azinovic I, Jeremic B. Effect of Radiotherapy Delay in Overall Treatment Time on Local Control and Survival in Head and Neck Cancer: Review of the Literature. Rep Pract Oncol Radiother (2015) 20(5):328–39. doi: 10.1016/j.rpor.2015.05.010
4. Levy A, Blanchard P, Daly-Schveitzer N, Tao Y. Replanning During Intensity Modulated Radiation Therapy Improved Quality of Life in Patients With Nasopharyngeal Carcinoma: In Regard to Yang et al. Int J Radiat Oncol Biol Phys (2013) 86(5):811. doi: 10.1016/j.ijrobp.2013.04.039
5. Kearney V, Chan J, Valdes G, Solberg T, Yom S. The Application of Artificial Intelligence in the IMRT Planning Process for Head and Neck Cancer. Oral Oncol (2018) 87:111–6. doi: 10.1016/j.oraloncology.2018.10.026
6. Wang C, Zhu X, Hong J, Zheng D. Artificial Intelligence in Radiotherapy Treatment Planning: Present and Future. Technol Cancer Res Treat (2019) 18:1533033819873922. doi: 10.1177/1533033819873922
7. Wu B, Ricchetti F, Sanguineti G, Kazhdan M, Simari P, Jacques R, et al. Data-Driven Approach to Generating Achievable Dose-Volume Histogram Objectives in Intensity-Modulated Radiotherapy Planning. Int J Radiat Oncol Biol Phys (2011) 79(4):1241–7. doi: 10.1016/j.ijrobp.2010.05.026
8. Zhu X, Ge Y, Li T, Thongphiew D, Yin F, Wu Q. A Planning Quality Evaluation Tool for Prostate Adaptive IMRT Based on Machine Learning. Med Phys (2011) 38(2):719–26. doi: 10.1118/1.3539749
9. Appenzoller L, Michalski J, Thorstad W, Mutic S, Moore K. Predicting Dose-Volume Histograms for Organs-at-Risk in IMRT Planning. Med Phys (2012) 39(12):7446–61. doi: 10.1118/1.4761864
10. Zarepisheh M, Long T, Li N, Tian Z, Romeijn H, Jia X, et al. A DVH-Guided IMRT Optimization Algorithm for Automatic Treatment Planning and Adaptive Radiotherapy Replanning. Med Phys (2014) 41(6):061711. doi: 10.1118/1.4875700
11. Zhuang Y, Xie Y, Wang L, Huang S, Chen L, Wang Y. Dvh Prediction for VMAT in NPC With GRU-RNN: An Improved Method by Considering Biological Effects. BioMed Res Int (2021) 2021:2043830. doi: 10.1155/2021/2043830
12. Huang Y, Yue H, Wang M, Li S, Zhang J, Liu Z, et al. Fully Automated Searching for the Optimal VMAT Jaw Settings Based on Eclipse Scripting Application Programming Interface (ESAPI) and RapidPlan Knowledge-Based Planning. J Appl Clin Med Phys (2018) 19(3):177–82. doi: 10.1002/acm2.12313
13. Ling C, Han X, Zhai P, Xu H, Chen J, Wang J, et al. A Hybrid Automated Treatment Planning Solution for Esophageal Cancer. Radiat Oncol (Lond Engl) (2019) 14(1):232. doi: 10.1186/s13014-019-1443-5
14. Good D, Lo J, Lee W, Wu Q, Yin F, Das S. A Knowledge-Based Approach to Improving and Homogenizing Intensity Modulated Radiation Therapy Planning Quality Among Treatment Centers: An Example Application to Prostate Cancer Planning. Int J Radiat Oncol Biol Phys (2013) 87(1):176–81. doi: 10.1016/j.ijrobp.2013.03.015
15. Delaney A, Dahele M, Tol J, Kuijper I, Slotman B, Verbakel W. Using a Knowledge-Based Planning Solution to Select Patients for Proton Therapy. Radiother Oncol (2017) 124(2):263–70. doi: 10.1016/j.radonc.2017.03.020
16. Chang A, Hung A, Cheung F, Lee M, Chan O, Philips H, et al. Comparison of Planning Quality and Efficiency Between Conventional and Knowledge-based Algorithms in Nasopharyngeal Cancer Patients Using Intensity Modulated Radiation Therapy. Int J Radiat Oncol Biol Phys (2016) 95(3):981–90. doi: 10.1016/j.ijrobp.2016.02.017
17. Kubo K, Monzen H, Ishii K, Tamura M, Kawamorita R, Sumida I, et al. Dosimetric Comparison of RapidPlan and Manually Optimized Plans in Volumetric Modulated Arc Therapy for Prostate Cancer. Phys Med (2017) 44:199–204. doi: 10.1016/j.ejmp.2017.06.026
18. Faught A, Olsen L, Schubert L, Rusthoven C, Castillo E, Castillo R, et al. Functional-Guided Radiotherapy Using Knowledge-Based Planning. Radiother Oncol (2018) 129(3):494–8. doi: 10.1016/j.radonc.2018.03.025
19. Fusella M, Scaggion A, Pivato N, Rossato M, Zorz A, Paiusco M. Efficiently Train and Validate a RapidPlan Model Through APQM Scoring. Med Phys (2018) 45(6):2611–9. doi: 10.1002/mp.12896
20. Fan J, Wang J, Chen Z, Hu C, Zhang Z, Hu W. Automatic Treatment Planning Based on Three-Dimensional Dose Distribution Predicted From Deep Learning Technique. Med Phys (2019) 46(1):370–81. doi: 10.1002/mp.13271
21. Liu Z, Fan J, Li M, Yan H, Hu Z, Huang P, et al. A Deep Learning Method for Prediction of Three-Dimensional Dose Distribution of Helical Tomotherapy. Med Phys (2019) 46(5):1972–83. doi: 10.1002/mp.13490
22. Barragán-Montero A, Nguyen D, Lu W, Lin M, Norouzi-Kandalan R, Geets X, et al. Three-Dimensional Dose Prediction for Lung IMRT Patients With Deep Neural Networks: Robust Learning From Heterogeneous Beam Configurations. Med Phys (2019) 46(8):3679–91. doi: 10.1002/mp.13597
23. Zhou J, Peng Z, Song Y, Chang Y, Pei X, Sheng L, et al. A Method of Using Deep Learning to Predict Three-Dimensional Dose Distributions for Intensity-Modulated Radiotherapy of Rectal Cancer. J Appl Clin Med Phys (2020) 21(5):26–37. doi: 10.1002/acm2.12849
24. Wieser H, Cisternas E, Wahl N, Ulrich S, Stadler A, Mescher H, et al. Development of the Open-Source Dose Calculation and Optimization Toolkit Matrad. Med Phys (2017) 44(6):2556–68. doi: 10.1002/mp.12251
25. He K, Zhang X, Ren S, Sun J. Delving Deep Into Rectifiers: Surpassing Human-Level Performance on ImageNet Classification. IEEE Conference on Computer Vision and Pattern Recognition (2015) 2015:1026–34. doi: 10.1109/ICCV.2015.123
26. Hodapp N. The ICRU Report 83: Prescribing, Recording and Reporting Photon-Beam Intensity-Modulated Radiation Therapy (IMRT). Strahlenther Onkol (2012) 188(1):97–9. doi: 10.1007/s00066-011-0015-x
27. Paddick I. A Simple Scoring Ratio to Index the Conformity of Radiosurgical Treatment Plans. Technical Note. J Neurosurg (2000) 3:219–22. doi: 10.3171/jns.2000.93.supplement
28. Nguyen D, Long T, Jia X, Lu W, Gu X, Iqbal Z, et al. A Feasibility Study for Predicting Optimal Radiation Therapy Dose Distributions of Prostate Cancer Patients From Patient Anatomy Using Deep Learning. Sci Rep (2019) 9(1):1076. doi: 10.1038/s41598-018-37741-x
29. Nguyen D, Jia X, Sher D, Lin M, Iqbal Z, Liu H, et al. 3D Radiotherapy Dose Prediction on Head and Neck Cancer Patients With a Hierarchically Densely Connected U-Net Deep Learning Architecture. Phys Med Biol (2019) 64(6):065020. doi: 10.1088/1361-6560/ab039b
30. Kajikawa T, Kadoya N, Ito K, Takayama Y, Chiba T, Tomori S, et al. Automated Prediction of Dosimetric Eligibility of Patients With Prostate Cancer Undergoing Intensity-Modulated Radiation Therapy Using a Convolutional Neural Network. Radiol Phys Technol (2018) 11(3):320–7. doi: 10.1007/s12194-018-0472-3
31. Kearney V, Chan J, Haaf S, Descovich M, Solberg T. DoseNet: A Volumetric Dose Prediction Algorithm Using 3D Fully-Convolutional Neural Networks. Phys Med Biol (2018) 63(23):235022. doi: 10.1088/1361-6560/aaef74
Keywords: automated treatment planning, dose distribution prediction, clinical post-optimization strategies, rectal cancer, intensity-modulated radiotherapy
Citation: Zhong Y, Yu L, Zhao J, Fang Y, Yang Y, Wu Z, Wang J and Hu W (2021) Clinical Implementation of Automated Treatment Planning for Rectum Intensity-Modulated Radiotherapy Using Voxel-Based Dose Prediction and Post-Optimization Strategies. Front. Oncol. 11:697995. doi: 10.3389/fonc.2021.697995
Received: 20 April 2021; Accepted: 11 June 2021;
Published: 24 June 2021.
Edited by:
Ning Wen, Henry Ford Health System, United StatesReviewed by:
Yibao Zhang, Peking University Cancer Hospital, ChinaYuenan Wang, Peking University Shenzhen Hospital, China
Copyright © 2021 Zhong, Yu, Zhao, Fang, Yang, Wu, Wang and Hu. 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: Jiazhou Wang, d2ppYXpob3VAZ21haWwuY29t; Weigang Hu, amFja2h1d2dAZ21haWwuY29t
†These authors have contributed equally to this work and share first authorship