- Open Access
Knowledge-based IMRT planning for individual liver cancer patients using a novel specific model
Radiation Oncology volume 13, Article number: 52 (2018)
The purpose of this work is to benchmark RapidPlan against clinical plans for liver Intensity-modulated radiotherapy (IMRT) treatment of patients with special anatomical characteristics, and to investigate the prediction capability of the general model (Model-G) versus our specific model (Model-S).
A library consisting of 60 liver cancer patients with IMRT planning was used to set up two models (Model-S, Model-G), using the RapidPlan knowledge-based planning system. Model-S consisted of 30 patients with special anatomical characteristics where the distance from planning target volume (PTV) to the right kidney was less than three centimeters and Model-G was configurated using all 60 patients in this library. Knowledge-based IMRT plans were created for the evaluation group formed of 13 patients similar to those included in Model-S by Model-G, Model-S and manually (M), named RPG-plans, RPS-plans and M-plans, respectively. The differences in the dose-volume histograms (DVHs) were compared, not only between RP-plans and their respective M-plans, but also between RPG-plans and RPS-plans.
For all 13 patients, RapidPlan could automatically produce clinically acceptable plans. Comparing RP-plans to M-plans, RP-plans improved V95% of PTV and had greater dose sparing in the right kidney. For the normal liver, RPG-plans delivered similar doses, while RPS-plans delivered a higher dose than M-plans. With respect to RapidPlan models, RPS-plans had better conformity index (CI) values and delivered lower doses to the right kidney V20Gy and maximizing point doses to spinal cord, while delivering higher doses to the normal liver.
The study shows that RapidPlan can create high-quality plans, and our specific model can improve the CI of PTV, resulting in more sparing of OAR in IMRT for individual liver cancer patients.
Liver cancer is the fifth most common type of cancer and the third leading cause of cancer-related death worldwide. Consequently, liver cancer is an issue that needs to be urgently addressed . Intensity-modulated radiotherapy (IMRT), based on computerized treatment plan optimization, permits the delivery of higher therapeutic doses to the target volume, while reducing the impact on adjacent normal tissue . However, each step in the clinical workflow, from contouring to delivery, contains variability and uncertainties that ultimately translate into inconsistencies and inefficiency [3, 4].
A variety of solutions were studied, to reduce the variability present the back and forth between plan creation and approval [5,6,7,8,9,10]. Among these methods, knowledge-based planning is the most commonly used. It can predict achievable planning target volume (PTV) and DVHs for organs-at-risk (OARs) for prospective patients, utilizing models built using a library consisting of previous plans [11,12,13,14,15,16,17,18,19]. Chanyavanich et al.  used an algorithm based on mutual information to identify similar patients, and then generate new plans for the target cases. Zhu et al.  developed an evaluation tool for quantification, that generates DVHs based on organ volumes, as well as distance-to-target histograms (DTH), using a machine learning approach. Wu et al.  predicted DVHs of target cases by establishing the overlap volume histogram (OVH). Based on that, Varian developed the Geometry-based Expected Dose (GED) algorithm and build a commercially available knowledge-based planning solution (RapidPlan, Varian Medical Systems, Palo Alto, CA), which semi-automates the treatment planning process. It uses a library of previous plans to build a model that can achieve OAR DVHs prediction range for a new patient, and subsequently guides the volumetric modulated arc therapy (VMAT) or IMRT optimization process using the Eclipse treatment planning system, along the inferior boundary of the DVH-prediction range. Previous work suggested that RapidPlan achieved clinically acceptable plans for different treatment sites [21,22,23,24,25,26,27]. To verify whether a model is suboptimal, Hussein et al.  suggested that there should be an insignificant effect on resulting plan quality when removing dosimetric outliers from the model training set. Similarly, Delaney et al.  found that statistical outliers removed from or added to models (5-10 outliers) had only a marginal impact on plan results. In a study by Tol et al. , a model created using 30 plans generated plans that were similar to those in a model based on 60 plans, when their plans were selected arbitrarily from all of 90 plans.
In this study, two models were developed, one general model and one specific model. All the studies mentioned above are based on the general model, which suits a wide variety of patient cases. However, the prediction ability of specific models tailored for specific patient cases, with specific anatomical features, has not been investigated. The aim of this study was (1) to evaluate the accuracy of RapidPlan prediction capability in IMRT (Eclipse, Varian) plans for individual liver cancer patients by using model libraries consisting of different total number of plans, with different similarity; and (2) to investigate the prediction capability of the general model vs. our specific model.
Liver cancer patients were treated with IMRT from 2015 to 2016, planned based on the Eclipse treatment planning system. For all patients, the IMRT plans were created using 5 non-uniformly distributed coplanar fields (0°, 200°, 240°, 280°, 320°) with the same photon beam setting of DVO, including number and beams energy. The prescribed dose was set as 50 Gy in 25 fractions. All plans were normalized to a mean dose of PTV (and isodose 95% was set to the prescribed dose) in order to make plan comparisons valid. OARs planning goals included maximizing point doses to the spinal cord and their planned at-risk volumes (3-mm expansion) below their tolerance dose levels, while lowering the endpoint dose to normal liver tissue and right kidney as much as possible. The PTV dose-volume constraints and OARs dose constraints are shown in Table 1. All clinical plans were manually optimized by an expert liver dosimetrist and each IMRT plan met the clinical protocol. Patients exhibited high variability in tumor position and size, as well as OAR exposure. All optimization and dose calculations were performed the dose volume optimization (DVO) version 13.5.35 and the anisotropic analytical algorithm (AAA) version 13.5.35 with a calculation grid of 2.5 mm. IMRT planning in Eclipse creates highly conformal dose distributions for liver cancer by continuously optimizing the beam intensity modulation to satisfy the institutional optimization protocol. Eclipse IMRT planning combines intensity modulation and inverse planning to accomplish this goal .
Model library and DVH estimation model configuration
The model library consisted of 60 liver cancer patients treated as above. From this library, all 60 patients were selected for Model-G. The average target volume was 147.2 ± 83.6 cm3 (range: 15.7-298.3 cm3). No other specific criteria were applied. For Model-S, 30 patients were selected, with the distance from PTV to the right kidney of less than 3 cm. The average target volume was 155.69 ± 77.5 cm3 (range: 30.4-306 cm3). Table 2 shows the volumes of the PTV, liver, L-PTV (normal liver) and right kidney, included in the two models.
RapidPlan was used to create a knowledge-based model that predicts achievable ranges of DVHs for individual OARs of prospective patients. The model libraries contain all planning CTs, structure sets, and dose distributions of previously treated patients. The model configuration should conceptually consist of multiple kinds of PTV and OAR geometries and alterations in OAR dosimetry, resulting exclusively from geometric alterations.
During the training, the system analyzed the patient anatomy and DVHs in the plans using principal component analysis (PCA), and created the final mathematical DVH estimation model. Then, the results of the model training were verified using statistical presentations of the training set. Regression, residual and DVH-plots help in estimating the quality of the model and finding potential outlier values that differ from the average in the training set . The outliers must be processed, and after that, the plan data are re-extracted and the model is retrained iteratively until the results are acceptable.
An OAR is designated as an outlier when one or more of these metrics lie outside the range of values found in the model. Using this strategy, 14 of 60 patients in Model-G and 6 of 30 patients in Model-S were identified as containing 1 outlier OAR. RapidPlan requires each OAR to be present in at least 20 plans included in the model. We removed statistical outliers from the training set, rather than deleting the whole plan. Table 3 shows the number of structures included in the two models.
An evaluation group consisting of 13 previous patients, treated in 2016, was used to test the RapidPlan results. Patients in the evaluation group were not included in the RapidPlan model libraries. The evaluation group was similar to that included in Model-S, where the distance from the right kidney to PTV was less than 3 cm. In the evaluation group, the average target volume was 163.1 ± 73.9 cm3 (range: 24-287.5 cm3). OARs typically included the liver, L-PTV (normal liver), the spinal-cord and their planned at-risk volumes (3-mm expansion). The volumes of the PTV, liver, L-PTV and right kidney, included in the evaluation group, are shown in Table 2.
The two models were used to generate optimization objectives and automatically optimize treatments for patients in the evaluation group, using the Eclipse treatment planning system. The PTV and OAR objectives are shown in Table 4. Optimization and dose calculation were performed using the Photon Optimization (PO) version 13.5.35 and AAA version 13.5.35. Knowledge-based IMRT plans were created for the evaluation group by Model-G and Model-S, named RPG-plans and RPS-plans, respectively, collectively termed RP-plans. Furthermore, these plans were manually optimized by an expert liver physicist, defined as M-plans.
Evaluating the performance of RapidPlan
Comparisons of the differences in the DVHs, not only between RP-plans and their respective M-plans, but also between RPG-plans and RPS-plans were performed. RapidPlan results were compared based on target dose coverage and normal tissue sparing. The target dose coverage includes: (1) the homogeneity index (HI) calculated for PTV using [HI = ((D2% - D98%)/Dp) × 100%], where Dp = prescribed dose; a lower HI value indicates that the dose coverage is more homogeneous ; (2) the conformity index (CI) proposed by Nakamura et al. : CI = TV × PIV/TVPIV2, where TV = target volume, PIV = prescribed isodose volume, and TVPIV = target volume receiving the prescribed dose. The dose coverage and conformity are better when the value of CI is closer to 1. The normal tissue sparing is based on statistical average doses to normal liver (Dmean, V30Gy V30Gy, V40Gy), spinal-cord (Dmax), and right kidney (Dmean, V5Gy, V10Gy, V15Gy). Paired t-tests were performed to determine significant differences (p < 0.05) between RP-plans and their respective M-plans. The target and normal tissue constraints shown in Table 1 were used to compare all patients in the evaluation group.
All knowledge-based plans were deemed to conform with the liver IMRT clinical protocol used at our institution, regarding dose-volume constraints. Table 5 and Fig. 1 summarize the RapidPlan results for the evaluation group, averaged over all patients, whereas Fig. 2 shows results for individual patients.
The comparison of RP-plans to M-plans revealed that RP-plans significantly improved V95% of PTV (RPG-plans vs. RPS-plans vs. M-plans: 98.4% vs. 98.4% vs. 97.9%). The D98%, D2%, HI and CI values of PTV were improved, although most differences were not significant. For the OAR, RP-plans had significantly more sparing of the right kidney (Dmean: RPG-plans vs. RPS-plans vs. M-plans, 3.9 Gy vs. 3.8 Gy vs. 4.5 Gy). For the normal liver, RPG-plans delivered similar doses, while RPS-plans delivered higher doses than those delivered by M-plans (Dmean: RPG-plans vs. RPS-plans vs. M-plans, 10.8 Gy vs. 11.1 Gy vs. 10.8 Gy). No difference was found between RP-plans and M-plans in the maximum dose to the spinal cord.
With respect to RapidPlan models, there are different. For the PTV coverage, the V95%, D98%, D2%, and HI values were very similar, but RPS-plans had better CI values (1.117 vs. 1.102, p < 0.05). For OAR doses paring, RPS-plans delivered a lower dose than the maximizing point doses to spinal cord (19.0 Gy vs. 20.7 Gy). Regarding the dose sparing for the right kidney, RPS-plans were better than RPG-plans, especially with V20Gy (4.5% vs. 5.6%, p < 0.05), although others were very similar. Figure 1 shows the advantages of RPS-plans in the DVH curves for the right kidney and spinal cord. For normal liver, RPS-plans delivered a higher dose (Dmean: 10.8 Gy vs. 11.1 Gy, p < 0.05; V20Gy: 20.1% vs. 22.6%, p < 0.05), although all RPS-plans were deemed conform with the liver IMRT clinical protocol. However, for certain individual patients, such as patient 11 (Fig. 2), RPG-plans are better for normal liver sparing.
Fogliata et al.  tested RapidPlan for the optimization of RapidArc plans, and generated clinically acceptable plans for hepatocellular cancer radiotherapy. This shows that the model is reliable when no special selection criteria are applied to generate the training, i.e. including all cases, for which the only requirement is to be clinically acceptable. However, liver cancer has high variability in tumor size and position, and the general model (Model-G) may have limited accuracy for special patients. Thus, a specific model (Model-S) was established in this study, using plans with specific anatomical features for individual patients. For Model-G, Jol et al.  demonstrated that 30 plans were sufficient for building a general model. However, in this study, Model-G consisted of 60 patients, to ensure the variety of the data, whereas Model-S consisted of 30 patients, to guarantee high similarity in the geometry of the region of interest. The scope of the minimum reasonable sample size will be addressed in further studies. Nevertheless, Jim et al.  observed that more OAR outliers did not necessarily translate into a worse OAR dose. Hussein et al.  showed that there were insignificant effects on resulting plan quality when removing dosimetric outliers from the model training set.
In this study, we evaluated the prediction capability of two models (Model-G & Model-S) in knowledge-based IMRT planning for individual liver cancer patients. Pooled results show that, generally, RapidPlan can improve planning quality and efficiency for liver IMRT, and the prediction ability of the two models with different configurations have a remarkable difference. RP-plans significantly improved the target coverage and the sparing of the right kidney compared to the M-plans. The advantages of some RP-plans compared to the M-plans may be due to challenges in performing an interactive planning of plans that contain many OARs optimally and consistently in a limited number of iterations. Comparison of the two models revealed that Model-S improved the CI and delivered lower dose to the right kidney (V20Gy) and spinal cord, while Model-G delivered lower dose to normal liver tissue. This indicates that the RP-plans are sensitive to the configuration of the model library and the anatomical characteristics of the patient that knowledge-based planning is performed on. The degree of similarity of the cases that make up the model library has a significant effect on the predictive capability of the model.
The selection criteria for establishing the specific model in this study was that the distance from PTV to the right kidney should measure less than 3 cm. We expected to better protect the right kidney, while improving the target coverage, and the experimental results show that Model-S achieved substantial gains compared to Model-G. The pooled data (Table 5 and Fig. 2) shows that RPS-plans slightly improved the right kidney sparing. However, for certain individual patients, RPS-plans had a significant advantage when compared to RPG-plans. For example, in patient 3 (Fig. 1), high variability in right kidney sparing was obtained, RPS-plans decreased V10Gy, V15Gy, V20Gy and Dmean to the right kidney compared to RPG-plans. In patient 11 (Fig. 2), RPG-plans are better for normal liver sparing. From an anatomical point of view, the reason may be that the PTV was small (24 cm3), and located in the inferior segment of the left side of the liver. In addition, it is encouraging that RPS-plans significantly reduced the maximum dose of the spinal cord compared to RPG-plans. This may be due to the fact that the spinal cord has a relatively fixed geometry, adjacent to the kidneys. Therefore, the specific model has great potential in clinical applications for individual patients and we will focus on this.
Some of the findings indicate the different tradeoffs in knowledge-based planning results. According to the pooled results of dosimetry, Model-G is better than the Model-S for normal liver sparing. The following factors may contribute to this: (1) the large size of the liver, only considering the cases with short distance between PTV and the right kidney may not meet the requirements for livers located far away from right kidney; (2) Model-G consisted of more cases, which increases geometric heterogeneity. It is conceivable that further improvements could be made and/or some complex tradeoffs should be addressed. Therefore, some parameters may be adjusted during the optimization process. It is worth noting that a review of plans used for RapidPlan models is still required.
RapidPlan results depend on a variety of factors and the following two are foremost, the geometry of the region of interest and the quality of the plans contained in the model libraries. To ensure the unity of the variables, we built a specific model, which only considers the geometric distance from the PTV to the right kidney, with no restrictions on the other OAR. In addition, one limitation of the study is the small sample size of the evaluation group (n = 13). Our study can provide some guidance for clinical applications. Future research will focus on providing optimal allocation of model libraries for individual patients.
This study shows that RapidPlan can create high quality plans and significantly improve the planning efficiency of IMRT for individual liver cancer patients. Furthermore, these findings demonstrate that the specific model can result in more sparing of OAR, while increasing the conformity index of PTV for liver cancer. Although more systematic studies are needed before a broad clinical application of the proposed methodology, this specific model might be considered as a way to improve the planning quality. Further studies are needed to determine the optimal composition of model libraries, including the relationship between model composition and dosimetry of subsequent plans.
Anisotropic analytical algorithm
Dose volume optimization
Geometry-based expected dose
Plans were created by manually
Overlap volume histogram
Principal component analysis
Planning target volume
Plans were created by Model-G
Plans were created by Model-S
Volumetric modulated arc therapy
Kudo M. Surveillance, diagnosis, treatment, and outcome of liver cancer in Japan. Liver Cancer. 2015;4:39–50.
Djajaputra D, Wu Q, Wu Y, Mohan R. Algorithm and performance of a clinical IMRT beam-angle optimization system. Phys Med Biol. 2003;48:3191.
Nelms BE, 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:296–305.
Das IJ, Cheng C-W, Chopra KL, Mitra RK, Srivastava SP, Glatstein E. Intensity-modulated radiation therapy dose prescription, recording, and delivery: patterns of variability among institutions and treatment planning systems. J Natl Cancer I. 2008;100:300–7.
Voet PW, Dirkx ML, Breedveld S, Fransen D, Levendag PC, Heijmen BJ. Toward fully automated multicriterial plan generation: a prospective clinical study. Int J Radiat Oncol Biol Phys. 2013;85:866–72.
Hong TS, Craft DL, Carlsson F, Bortfeld TR. Multicriteria optimization in intensity-modulated radiation therapy treatment planning for locally advanced cancer of the pancreatic head. Int J Radiat Oncol Biol Phys. 2008;72:1208–14.
Craft DL, Hong TS, Shih HA, Bortfeld TR. Improved planning time and plan quality through multicriteria optimization for intensity-modulated radiotherapy. Int J Radiat Oncol Biol Phys. 2012;82:e83–90.
Chen H, Craft DL, Gierga DP. Multicriteria optimization informed VMAT planning. Med Dosim. 2014;39:64–73.
Wu B, Pang D, Simari P, Taylor R, Sanguineti G, McNutt T. Using overlap volume histogram and IMRT plan data to guide and automate VMAT planning: a head-and-neck case study. Med Phys. 2013;40:021714.
Breedveld S, Storchi PR, Voet PW, Heijmen BJ. iCycle: integrated, multicriterial beam angle, and profile optimization for generation of coplanar and noncoplanar IMRT plans. Med Phys. 2012;39:951–63.
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:011709.
Wu B, McNutt T, Zahurak M, Simari P, Pang D, Taylor R, et al. Fully automated simultaneous integrated boosted–intensity modulated radiation therapy treatment planning is feasible for head-and-neck cancer: a prospective clinical study. Int J Radiat Oncol Biol Phys. 2012;84:e647–e53.
Lian J, Yuan L, Ge Y, Chera BS, Yoo DP, Chang S, et al. Modeling the dosimetry of organ-at-risk in head and neck IMRT planning: an intertechnique and interinstitutional study. Med Phys. 2013;40:121704.
Liu BJ. A knowledge-based imaging informatics approach for managing proton beam therapy of cancer patients. Technol CANCER Res T. 2007;6:77–84.
Lee KJ, Barber DC, Walton L. Automated gamma knife radiosurgery treatment planning with image registration, data-mining, and Nelder-mead simplex optimization. Med Phys. 2006;33:2532–40.
Wu B, Ricchetti F, Sanguineti G, Kazhdan M, Simari P, Chuang M, et al. Patient geometry-driven information retrieval for IMRT treatment plan quality control. Med Phys. 2009;36:5497–505.
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:1241–7.
Chanyavanich V, Das SK, Lee WR, Lo JY. Knowledge-based IMRT treatment planning for prostate cancer. Med Phys. 2011;38:2515–22.
Good D, Lo J, Lee WR, Wu QJ, Yin F-F, Das SK. 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:176–81.
Zhu X, Ge Y, Li T, Thongphiew D, Yin F-F, Wu QJ. A planning quality evaluation tool for prostate adaptive IMRT based on machine learning. Med Phys. 2011;38:719–26.
Fogliata A, Wang P-M, Belosi F, Clivio A, Nicolini G, Vanetti E, et al. Assessment of a model based optimization engine for volumetric modulated arc therapy for patients with advanced hepatocellular cancer. Radiat Oncol. 2014;9:236–48.
Fogliata A, Belosi F, Clivio A, Navarria P, Nicolini G, Scorsetti M, et al. On the pre-clinical validation of a commercial model-based optimisation engine: application to volumetric modulated arc therapy for patients with lung or prostate cancer. Radiother Oncol. 2014;113:385–91.
Wu H, Jiang F, Yue H, Zhang H, Wang K, Zhang Y. Applying a RapidPlan model trained on a technique and orientation to another: a feasibility and dosimetric evaluation. Radiat Oncol. 2016;11:108.
Hussein M, South CP, Barry MA, Adams EJ, Jordan TJ, Stewart AJ, et al. Clinical validation and benchmarking of knowledge-based IMRT and VMAT treatment planning in pelvic anatomy. Radiother Oncol. 2016;120:473–9.
Chang AT, Hung AW, Cheung FW, Lee MC, Chan OS, 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:981–90.
Berry SL, Ma R, Boczkowski A, Jackson A, Zhang P, Hunt M. Evaluating inter-campus plan consistency using a knowledge based planning model. Radiother Oncol. 2016;120:349–55.
Tol JP, Delaney AR, Dahele M, Slotman BJ, Verbakel WF. Evaluation of a knowledge-based planning solution for head and neck cancer. Int J Radiat Oncol Biol Phys. 2015;91:612–20.
Delaney AR, Tol JP, Dahele M, Cuijpers J, Slotman BJ, Verbakel WF. Effect of dosimetric outliers on the performance of a commercial knowledge-based planning solution. Int J Radiat Oncol Biol Phys. 2016;94:469–77.
Nakamura JL, Verhey LJ, Smith V, Petti PL, Lamborn KR, Larson DA, et al. Dose conformity of gamma knife radiosurgery and risk factors for complications. Int J Radiat Oncol Biol Phys. 2001;51:1313–9.
We want to thank Dr. Lei Xing and Dr. Yong Yang (Department of Radiation Oncology, Stanford University School of Medicine, Stanford, CA, USA) for the fruitful discussions of knowledge-based models and optimization.
This work was jointly supported by National Natural Science Foundation of China (61471226, 81530060), Natural Science Foundation for Distinguished Young Scholars of Shandong Province (JQ201516), Natural Science Foundation of Shandong Province (ZR2017BF005). We also thanks to the supporting of Taishan scholar project of Shandong Province (NO.tsqn20161023).
Availability of data and materials
All data generated or analysed during this study are included in this published article.
Ethics approval and consent to participate
Consent for publication
Written informed consent was obtained from the patient for the publication of this report.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yu, G., Li, Y., Feng, Z. et al. Knowledge-based IMRT planning for individual liver cancer patients using a novel specific model. Radiat Oncol 13, 52 (2018) doi:10.1186/s13014-018-0996-z
- Specific model
- Knowledge-based planning
- Liver cancer