- 1Department of Radiology, The Affiliated Hospital of Southwest Medical University, Luzhou, China
- 2Department of Oncology, The Affiliated Hospital of Southwest Medical University, Luzhou, China
- 3Department of Nursing College, Southwest Medical University, Luzhou, China
Introduction: This study aimed to establish a support vector machine (SVM) model to predict the dose for organs at risk (OARs) in intracavitary brachytherapy planning for cervical cancer with tandem and ovoid treatments.
Methods: Fifty patients with loco-regionally advanced cervical cancer treated with 200 CT-based tandem and ovoid brachytherapy plans were included. The brachytherapy plans were randomly divided into the training (N = 160) and verification groups (N = 40). The bladder, rectum, sigmoid colon, and small intestine were divided into sub-OARs. The SVM model was established using MATLAB software based on the sub-OAR volume to predict the bladder, rectum, sigmoid colon, and small intestine . Model performance was quantified by mean squared error (MSE) and δ . The goodness of fit of the model was quantified by the coefficient of determination (R2). The accuracy and validity of the SVM model were verified using the validation group.
Results: The value of the bladder, rectum, sigmoid colon, and small intestine correlated with the volume of the corresponding sub-OARs in the training group. The mean squared error (MSE) in the SVM model training group was <0.05; the R2 of each OAR was >0.9. There was no significant difference between the -predicted and actual values in the validation group (all P > 0.05): bladder δ = 0.024 ± 0.022, rectum δ = 0.026 ± 0.014, sigmoid colon δ = 0.035 ± 0.023, and small intestine δ = 0.032 ± 0.025.
Conclusion: The SVM model established in this study can effectively predict the for the bladder, rectum, sigmoid colon, and small intestine in cervical cancer brachytherapy.
Introduction
Cervical cancer is the most common malignancy among women in developing countries (1). Depending on the stage of diagnosis, the treatment strategies for cervical cancer mainly include surgery, along with radiotherapy and chemotherapy (2). For patients with locally advanced cervical cancer, brachytherapy combined with external-beam radiotherapy is the prevalent standard treatment (3). Three-dimensional brachytherapy is widely applied in clinical practice, and computed tomography (CT)- or magnetic resonance imaging (MRI)-based treatment planning systems (TPS) provide accurate tumor and organs at risk (OARs) dose information. However, the experience of brachytherapy planners and knowledge of the Radiation Therapy Oncology Group guidelines, as well as clinical expertise and intuition, have a significant effect on the quality of a brachytherapy plan (4). If a planner can predict the OAR dose before designing a brachytherapy plan, the quality of the brachytherapy plan can be controlled, and the interfering factors can be minimized. Previous reports on cervical cancer brachytherapy have described the effects of the volume of the OARs on the dose to the bladder, rectum, sigmoid colon, and small intestine (5). Although there is a correlation between the dose to the OARs and their volumes, information to predict the dose to the OARs is limited. In recent years, methods for predicting the dose to the OARs have been widely introduced in external irradiation intensity-modulated radiotherapy (6–10). These approaches typically use libraries of existing patient plans to create models that predict the extent of OAR sparing that can be achieved in a new patient based on, for example, the planning target volume (PTV)-OAR distance and overlap (11). In this study, we examined factors relevant for the dose to the bladder, rectum, sigmoid colon, and small intestine in cervical cancer brachytherapy based on the Fletcher applicator. The bladder, rectum, sigmoid colon, and small intestine were divided into sub-OARs. We analyzed the correlation between the sub-organ volume and of each OAR, and the SVM prediction model based on the correlation was established to predict the dose of each OAR before brachytherapy; the model can be used as an evaluation standard for brachytherapy plans to minimize the effects of confounding factors on the quality of the plans. To our knowledge, this study is the first to apply the SVM model to OAR dosimetric prediction based only on the contours of the organs and targets. This approach has been granted a Chinese invention patent (patent no.: 201610529290.8).
Materials and Methods
Patients
We retrospectively selected 50 patients with loco-regionally advanced cervical cancer treated with 200 CT-based tandem and ovoid brachytherapy plans between 2016 and 2018 in the Affiliated Hospital of Southwest Medical University. The patients treated with brachytherapy were randomly divided into the training (N = 160) and verification groups (N = 40). The cervical cancer stages ranged from ІІB to IVA, according to the International Federation of Gynecology and Obstetrics system.
Targets and Delineation of the OARs
The high-risk clinical target volume (HR-CTV) contours were generated for each treatment based on the Gynaecological European Society for Radiotherapy and Oncology Working Group I (Gyn GEC-ESTRO WG I) recommendations (12). The HR-CTV covered the entire cervix and macroscopic extent of the disease, based on clinical examinations and as depicted in CT images. The OARs included the bladder, rectum, sigmoid colon, and small intestine. The same radiation oncologist performed the target and delineation of the OARs.
Prescription Dose and Limiting Requirements for the OARs
After receiving 45 Gy intensity-modulated radiation therapy (IMRT), the per fraction prescription dose (Dprescription) for the HR-CTV was defined as 7 Gy with a total of four fractions for brachytherapy. A prescription dose delivered to 90% of the HR-CTV was considered. Combined with the IMRT dose, the total EQD2 (equivalent dose in 2 Gy, α/β = 10) for HR-CTV and IR-CTV was 85 and 60 Gy, respectively. We applied dose constraints for the OARs according to the following principles: combined IMRT dose, of EQD2 of ≤90 Gy (α/β = 3) for bladder, ≤75 Gy (α/β = 3) for rectum, ≤75 Gy (α/β = 3) for sigmoid colon, and ≤75 Gy (α/β = 3) for small intestine. These dose constraints were primarily based on the Gynaecological European Society for Radiotherapy and Oncology Working Group II (Gyn GEC-ESTRO WG II) recommendations (13). The 192Ir-source was delivered using the Fletcher applicator. To avoid bladder and rectum volume variations, the bladder of all patients was emptied and subsequently filled with 50 ml of saline solution; they accepted an enema to empty the rectum before brachytherapy.
Brachytherapy Plans
The Oncentra 4.3 treatment planning system (Elekta Brachytherapy, Veenendaal, the Netherlands) was used for the brachytherapy plans. All brachytherapy plans in this study were developed using a manual and/or graphical optimization approach to repeatedly optimize the plan and thus ensure that the dose administered to 90% of the HR-CTV reached the prescribed dose (Dprescription), whereas the dose to the OARs was lower. For the optimization of the single brachytherapy plan, the prescription dose (7 and 4.2 Gy) was administered to 90% of HR-CTV and IR-CTV; of the bladder < 5.2 Gy, of the rectum, sigmoid colon, and small intestine < 4.7 Gy.
Deriving Sub-OARs From the OARs
The HR-CTV was externally expanded to a plurality of rings (ring1–ringn) with a width of 0.5 cm using the Oncentra 4.3 treatment planning system. Ring1–ringn and different OAR intersection regions (ring1–ringn∩OAR) were used as independent sub-OARs, with ring1∩OAR defined as the sub-OAR1, and so on; ringn∩OAR was defined as sub-OARn. The total sub-OARs are controlled within 10 and the statistics of the volume of each sub-OARs. The intersecting regions for ring1–ring9 and the bladder in patient 15 are shown in Figure 1.
Figure 1 The red line indicated HR-CTV, The green shadow indicated the intersection of ring1 and bladder, the blue shadow indicated the intersection of ring2 and bladder. The yellow shadow indicated the intersection of ring3 and bladder. The purple shadow indicated the intersection of ring4 and bladder. The sky blue shadow indicated the intersection of ring5 and bladder. The lavender shadow indicated the intersection of ring6 and bladder. The orange shadow indicated the intersection of ring7 and bladder. The forest shadow indicated the intersection of ring8 and bladder. The slate blue shadow indicated the intersection of ring9 and bladder.
SVM Model Development
In machine learning, support vector machine (SVM) are supervised learning models with the associated learning algorithms used to analyze data for classification and regression analysis. In our study, we applied a radial basis function kernel for binary classification. We used MATLAB (R2017a, MathWorks, Inc., Natick, MA, USA) software to read, prepare, process, and output the predicted value. The SVM models were trained, validated, and tested for prediction accuracy using a self-written algorithm in MATLAB. A common radial basis function kernel was used:
where xi and xj are two data points, and γ is the shape parameter that represents the equivalent to the standard deviation in Gaussian distribution. To deal with the problem of regularization for noisy data, a user-specified cost parameter C is introduced, which acts to soften the margin. The cost parameter C controls the trade-off between allowing transgression of data points across the margin edges toward the other class and a more complex boundary, which might lead to overfitting. The evaluation and choice of C and γ were conducted using a grid search. The optimal parameters were estimated using the training and validation sets. We analyze the correlation between the sub-organ volume and of each OAR and establish the SVM prediction model based on the correlation. The volumes of the sub-OARs were used as the independent variable in the SVM model, and the ratios were used as the dependent variable.
For the verification group, the performance of the SVM model was investigated to predict per fraction in the bladder, rectum, sigmoid, and small intestine using the volumes of the corresponding sub-OARs. The volumes of the sub-OARs were used as the input values for the SVM model, and the ratios were used as the output values. The performance of the model can be characterized by mean squared error (MSE) and δ . The goodness of fit of the model was quantified by the coefficient of determination (R2 = 1 − the ratio of the sum of squares regressed to the total sum of squares). R2 indicates the proportionate amount of variation in the response variable explained by the independent variables in the model. They measure the fitting performance of a model from different perspectives. The closer the δ is to 0, the closer the actual and prescription values are to each other. Furthermore, the closer the R2 is to 1, the higher the fitting degree.
Statistical Analysis
Significant differences were determined using a two-sided paired t-test with SPSS 19.0 software (SPSS, Inc., Chicago, IL, USA). Correlations were tested by performing the Pearson correlation coefficient analysis. P <0.05 indicates that there is a correlation between the two variables, and P <0.01 indicates that there is a significant correlation between the two variables.
Results
The volume of each sub-OAR (Vsub-OAR) was correlated with the of the respective OAR. The volume of the HR-CTV (VHR-CTV) was correlated with the of the bladder, rectum, and sigmoid colon (all correlations, P < 0.05). The volume of the bladder (Vbladder) and the of the small intestine were correlated. The correlation coefficient (r, a statistical index used to describe the degree of linear correlation between two variables), and P values are shown in Tables 1–4. Therefore, these data can be used to predict the of each OAR using the SVM model. The MSE and the R2 of each OAR in the SVM model prediction group are shown in Table 5.
The predicted and actual values for the bladder, rectum, sigmoid colon, and small intestine in the validation group are shown in Figure 2. There was no statistically significant difference between the predicted and actual values for the bladder (P = 0.68), rectum (P = 0.16), sigmoid colon (P = 0.14), and small intestine (P = 0.77) in the validation group. The δ value for the bladdera of the verification group was 0.024 ± 0.022, the corresponding rectum δ value was 0.026 ± 0.014, the sigmoid colon δ value was 0.035 ± 0.023, and the small intestine δ value was 0.032 ± 0.025.
Figure 2 The predicted and actual values for the bladder, rectum, sigmoid colon, and small intestine.
Discussion
The quality control of radiotherapy plan has always been a research hotspot in the field of radiotherapy (14–17). The most critical aspect is the prediction of the dose to the OAR before designing the radiotherapy plan. It has been reported that the OAR dose in the brachytherapy plan could be predicted by the overlapping volume of the OAR with the targeted area and knowledge-based tool (18, 19). Ours is a relatively simple mathematical model that uses prescription dose and Vsub-OAR to predict the bladder, rectum, and sigmoid for brachytherapy; this does not require buying new modules of TPS or extracting the distance of each sampling point of the OAR with the dose. We also divided the OARs into multiple sub-OARs to predict the OAR dose in the external IMRT plan (20, 21). In contrast to previous studies, the focus of our study is to determine the correlation between the Vsub-OAR and of each OAR in brachytherapy, therefore, this method has been granted the Chinese invention patent. Owing to this correlation, we could fit the data of the training group using the SVM model approach. To rule out the effects of different prescription doses on the of each OAR, we divided each by 90% of the HR-CTV that reached the dose (Dprescription).
As shown in Figure 2, our SVM estimation system predicted that the of the OARs is very close to the actual value. There was no significant difference between the predicted and actual values for each OAR. The δ values of the bladder, rectum, sigmoid colon, and small intestine were 0.024 ± 0.022, 0.026 ± 0.014, 0.035 ± 0.023, and 0.032 ± 0.025, respectively. The abovementioned statistics, MSE, and R2 of the SVM prediction model indicated that the prediction model was reliable. We used a relatively simple mathematical model, which does not require the acquisition of new modules of TPS software. The process model for the acquisition of the sub-OAR can be edited into scripts to improve efficiency and effectiveness.
Our model can be used as a component of a quality assurance tool to detect suboptimal treatment plans in OAR sparing. A properly trained model will provide an estimate of the OAR doses required for appropriate planning and will detect outlines that require further review. Specifically, considering δ, a δ value closer to 0 indicated a closer relationship between the planned and predicted values of . A standard δ threshold can be set for the of each OAR, and the value above the threshold should be further optimized or the position of the applicator should be re-adjusted, until a satisfactory δ value is obtained. Predictions using the SVM model can be conducted for the quality control of the brachytherapy plan and for minimizing the effect of subjective factors (22).
Our study has some limitations. It was restricted to a single institution and considered only standard tandem and ovoid cases. Further research is needed comprising multiple centers and more cervical cancer brachytherapy plan data sets for analysis. If the data set is large enough, a neural network model can be developed, which will generate predictions with higher accuracy of the OAR dose for cervical cancer brachytherapy plans. The SVM models discussed herein may be applied beyond gynecologic brachytherapy. The application of our models to prostate brachytherapy as well can be considered after validation.
Conclusion
The SVM model can be applied to not only predict the dose to the OARs for the high-dose rate brachytherapy of cervical cancer but also develop quality assurance tools for designing brachytherapy plans.
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 The Affiliated Hospital of Southwest Medical University Ethics Committee. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
Guarantors of integrity of entire study, HP. Study concepts/study design or data acquisition or data analysis/interpretation, all authors. Manuscript drafting or manuscript revision for important intellectual content, all authors. Resolution of any questions related to the work, all authors. Literature research, PZ, SL, and HP. Statistical analysis, PZ, XL, and HP. Manuscript editing, PZ, SL, and HP. All authors contributed to the article and approved the submitted version.
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.
References
1. Small W, Bacon MA, Bajaj A, Chuang LT, Fisher BJ, Harkenrider MM, et al. Cervical Cancer: A Global Health Crisis. Cancer (2017) 123(13):2404–12. doi: 10.1002/cncr.30667
2. Cibula D, Pötter R, Planchamp F, Avall-Lundqvist E, Fischerova D, Haie Meder C, et al. The European Society of Gynaecological Oncology/European Society for Radiotherapy and Oncology/European Society of Pathology Guidelines for the Management of Patients With Cervical Cancer. Int J Gynecol Cancer (2018) 28(4):641–55. doi: 10.1097/IGC.0000000000001216
3. Mayadev J, Viswanathan A, Liu Y, Li CS, Albuquerque K, Damato AL, et al. American Brachytherapy Task Group Report. American Brachytherapy Task Group Report: A Pooled Analysis of Clinical Outcomes for High-Dose-Rate Brachytherapy for Cervical Cancer. Brachytherapy (2017) 16(1):22–43. doi: 10.1016/j.brachy.2016.03.008
4. Wu B, Kusters M, Kunze-Busch M, Dijkema T, McNutt T, Sanguineti G, et al. Cross-Institutional Knowledge-Based Planning (KBP) Implementation and its Performance Comparison to Auto-Planning Engine (APE). Radiother Oncol (2017) 123(1):57–62. doi: 10.1016/j.radonc.2017.01.012
5. Siavashpour Z, Aghamiri MR, Jaberi R, Manshadi HR, Ghaderi R, Kirisits C. Optimum Organ Volume Ranges for Organs at Risk Dose in Cervical Cancer Intracavitary Brachytherapy. J Contemp Brachytherapy (2016) 8(2):135–42. doi: 10.5114/jcb.2016.59687
6. Moore KL, Brame RS, Low DA, Mutic S. Experience-Based Quality Control of Clinical Intensity-Modulated Radiotherapy Planning. Int J Radiat Oncol Biol Phys (2011) 81(2):545–51. doi: 10.1016/j.ijrobp.2010.11.030
7. Zawadzka A, Nesteruk M, Brzozowska B, Kukołowicz PF. Method of Predicting the Mean Lung Dose Based on a Patient’s Anatomy and Dose-Volume Histograms. Med Dosim (2017) 42(1):57–62. doi: 10.1016/j.meddos.2016.12.001
8. Tol JP, Dahele M, Delaney AR, Slotman BJ, Verbakel WF. Can Knowledge-Based DVH Predictions be Used for Automated, Individualized Quality Assurance of Radiotherapy Treatment Plans? Radiat Oncol (2015) 10(1):234. doi: 10.1186/s13014-015-0542-1
9. Ge Y, Wu QJ. Knowledge-Based Planning for Intensity-Modulated Radiation Therapy: A Review of Data-Driven Approaches. Med Phys (2019) 46(6):2760–75. doi: 10.1002/mp.13526
10. Jiao SX, Chen LX, Zhu JH, Wang ML, Liu XW. Prediction of Dose-Volume Histograms in Nasopharyngeal Cancer IMRT Using Geometric and Dosimetric Information. Phys Med Biol (2019) 64(23):23NT04. doi: 10.1088/1361-6560/ab50eb
11. Yang Y, Ford EC, Wu B, Pinkawa M, van Triest B, Campbell P, et al. An Overlap-Volume-Histogram Based Method for Rectal Dose Prediction and Automated Treatment Planning in the External Beam Prostate Radiotherapy Following Hydrogel Injection. Med Phys (2013) 40(1):011709. doi: 10.1118/1.4769424
12. Haie-Meder C, Pötter R, Van Limbergen E, Briot E, De Brabandere M, Dimopoulos J, et al. Recommendations From Gynaecological (GYN) GEC-ESTRO Working Group (I): Concepts and Terms in 3D Image Based 3D Treatment Planning in Cervix Cancer Brachytherapy With Emphasis on MRI Assessment of GTV and CTV. Radiother Oncol (2005) 74(3):235–45. doi: 10.1016/j.radonc.2004.12.015
13. Pötter R, Haie-Meder C, Van Limbergen E, Barillot I, De Brabandere M, Dimopoulos J, et al. Recommendations From Gynaecological (GYN) GEC ESTRO Working Group (II): Concepts and Terms in 3D Image-Based Treatment Planning in Cervix Cancer Brachytherapy-3D Dose Volume Parameters and Aspects of 3D Image-Based Anatomy, Radiation Physics, Radiobiology. Radiother Oncol (2006) 78(1):67–77. doi: 10.1016/j.radonc.2005.11.014
14. Yusufaly TI, Meyers SM, Mell LK, Moore KL. Knowledge-Based Planning for Intact Cervical Cancer. Semin Radiat Oncol (2020) 30(4):328–39. doi: 10.1016/j.semradonc.2020.05.009
15. Krayenbuehl J, Di Martino M, Guckenberger M, Andratschke N. Improved Plan Quality With Automated Radiotherapy Planning for Whole Brain With Hippocampus Sparing: A Comparison to the RTOG 0933 Trial. Radiat Oncol (2017) 12(1):161. doi: 10.1186/s13014-017-0896-7
16. Court LE, Kisling K, Mccarroll R, Zhang L, Yang J, Simonds H, et al. Radiation Planning Assistant - A Streamlined, Fully Automated Radiotherapy Treatment Planning System. J Vis Exp (2018) 134(134):57411. doi: 10.3791/57411
17. Buschmann M, Sharfo AWM, Penninkhof J, Seppenwoolde Y, Goldner G, Georg D, et al. Automated Volumetric Modulated Arc Therapy Planning for Whole Pelvic Prostate Radiotherapy. Strahlenther Onkol (2018) 194(4):333–42. doi: 10.1007/s00066-017-1246-2
18. Damato AL, Viswanathan AN, Cormack RA. Validation of Mathematical Models for the Prediction of Organs-at-Risk Dosimetric Metrics in High-Dose-Rate Gynecologic Interstitial Brachytherapy. Med Phys (2013) 40(10):101711. doi: 10.1118/1.4819946
19. Yusufaly TI, Kallis K, Simon A, Mayadev J, Yashar CM, Einck JP, et al. A Knowledge-Based Organ Dose Prediction Tool for Brachytherapy Treatment Planning of Patients With Cervical Cancer. Brachytherapy (2020) 19(5):624–34. doi: 10.1016/j.brachy.2020.04.008
20. Pang H, Sun X, Yang B, Wu J. A Quality Control Method for Intensity-Modulated Radiation Therapy Planning Based on Generalized Equivalent Uniform Dose. J Appl Clin Med Phys (2018) 19(3):276–82. doi: 10.1002/acm2.12331
21. Pang H, Sun X, Yang B, Wu J. Predicting the Dose Absorbed by Organs at Risk During Intensity Modulated Radiation Therapy for Nasopharyngeal Carcinoma. Br J Radiol (2018) 91(1092):20170289. doi: 10.1259/bjr.20170289
22. Klement RJ, Allgäuer M, Appold S, Dieckmann K, Ernst I, Ganswindt U, et al. Support Vector Machine-Based Prediction of Local Tumor Control After Stereotactic Body Radiation Therapy for Early-Stage Non-Small Cell Lung Cancer. Int J Radiat Oncol Biol Phys (2014) 88(3):732–8. doi: 10.1016/j.ijrobp.2013.11.216
Keywords: brachytherapy, cervical cancer, organs at risk, support vector machine, dose prediction
Citation: Zhou P, Li X, Zhou H, Fu X, Liu B, Zhang Y, Lin S and Pang H (2021) Support Vector Machine Model Predicts Dose for Organs at Risk in High-Dose Rate Brachytherapy of Cervical Cancer. Front. Oncol. 11:619384. doi: 10.3389/fonc.2021.619384
Received: 20 October 2020; Accepted: 25 June 2021;
Published: 15 July 2021.
Edited by:
Gene A. Cardarelli, Warren Alpert Medical School of Brown University, United StatesReviewed by:
Aaron Howard Wolfson, University of Miami, United StatesMark Rivard, Brown University, United States
Copyright © 2021 Zhou, Li, Zhou, Fu, Liu, Zhang, Lin and Pang. 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: Sheng Lin, lslinsheng@163.com; Haowen Pang, haowenpang@foxmail.com
†These authors have contributed equally to this work