- Open Access
Fully automated treatment planning of spinal metastases – A comparison to manual planning of Volumetric Modulated Arc Therapy for conventionally fractionated irradiation
Radiation Oncology volume 12, Article number: 33 (2017)
Planning for Volumetric Modulated Arc Therapy (VMAT) may be time consuming and its use is limited by available staff resources. Automated multicriterial treatment planning can eliminate this bottleneck. We compared automatically created (auto) VMAT plans generated by Erasmus-iCycle to manually created VMAT plans for treatment of spinal metastases.
Forty-two targets in 32 patients were analyzed. Lungs and kidneys were defined as organs at risk (OARs). Twenty-two patients received radiotherapy on kidney levels, 17 on lung levels, and 3 on both levels.
All Erasmus-iCycle plans were clinically acceptable. When compared to manual plans, planning target volume (PTV) coverage of auto plans was significantly better. The Homogeneity Index did not differ significantly between the groups. Mean dose to OARs was lower in auto plans concerning both kidneys and the left lung. One hotspot (>110% of D50%) occurred in the spinal cord of one auto plan (33.2 Gy, D50%: 30 Gy). Treatment time was 7% longer in auto plans.
Erasmus-iCycle plans showed better target coverage and sparing of OARs at the expense of minimally longer treatment times (for which no constraint was set).
Fractionated or single-dose radiotherapy to spinal metastases of solid tumors is one of the most frequently performed radiation treatments . Highly conformal radiation techniques such as Stereotactic Body Radiation Therapy (SBRT), Intensity Modulated Radiation Therapy (IMRT) or VMAT are increasingly being used as their application becomes easier and dose distributions confer theoretical advantages , albeit so far without proof of superior clinical outcome. On the other hand, generation of an optimal IMRT or VMAT plan is often time consuming and it has been shown that staff limitations are correlated with restricted use of new techniques such as IMRT even in the developed world [3, 4]. Automated planning of IMRT and VMAT may reduce the workload which is associated with manual “trial-and-error” approaches by around 50% .
Erasmus-iCycle, developed at the Erasmus MC-Cancer Institute, is an optimizer for multicriterial beam profile optimization and beam angle selection for coplanar and non-coplanar IMRT [5–10]. Other solutions for multicriterial beam angle optimization have been proposed by Schreibmann et al. , Craft and Monz presented a full multibeam space Pareto navigation tool . However, both are a posteriori methods, i.e., the algorithm generates different sets of beam angles and intensity profiles from which the user selects the plan afterward. The Erasmus-iCycle solution is an a priori approach which enables the user to define a site specific set of criteria (wish-list) which may not be violated (constraints) or have to be met as well as possible, or better (objectives). Objectives have assigned priorities to steer the multicriterial planning towards favourable trade-offs between the various treatment goals. If the priority of an objective is higher, the probability that the corresponding objective is met increases. Hard constraints are always respected in Erasmus-iCycle plans.
Beam directions are selected from candidate directions which can be restricted, e.g., in case of coplanar treatments [6, 7, 13]. For fully automated generation of plans that are clinically delivered, Erasmus-iCycle auto plans are automatically reconstructed and segmented in the clinical treatment planning system (TPS) . This study intends to validate VMAT plans for treatment of spinal metastases, generated with this approach, and to compare the quality of automatically generated plans with plans that were manually created by experienced treatment planners. Different spine regions pose different optimization problems as a consequence of different OARs being relevant for the treatment. Therefore cervical, thoracic, and lumbar targets were included in the study design.
Manual plan generation
Forty-two clinical target volumes (CTVs) in the spinal column of 32 patients were manually delineated for clinical routine treatments. CTVs were anisotropically expanded to PTVs that were the basis for all further analyses. All PTVs were reviewed by expert radiation oncologists (FW, FL).
Manual plans were created by expert treatment planners and also reviewed by radiation oncologists. All manually generated plans were calculated using the Monaco® treatment planning system (Elekta Ltd, Crawly, UK) version 3.2 or later which supports static IMRT, dynamic IMRT, and VMAT. Informed consent was obtained from all patients for anonymized processing of their clinical data. The study was approved by the ethics committee of Heidelberg University, Medical Faculty Mannheim (2016-806R-MA). For this study we assumed lungs or kidneys to be (potentially pre-irradiated) OARs, therefore we included only spinal regions which were at the level of kidneys (n = 22), lungs (n = 17), or both (n = 3), resulting in 90 (45*2) OARs and 42 PTVs in total. Spinal reirradiation plans were not tested in this study; nevertheless, we contoured the spinal cord to identify any hotspots in this area within the target volume. Most plans were at the thoracic, thoracolumbar or lumbar level. In 4 cases, soft tissue metastases or rib metastases were included in the PTV. Prescribed doses ranged from 30 to 40 Gy in 10–20 fractions which are commonly applied treatment regimens in patients with sufficient life expectancy and good general condition (Eastern Cooperative Oncology Group, ECOG 0/1; Karnofsky Performance Status, KPS > 80%). Generally, 40 Gy regimens were applied to patients who received postoperative radiotherapy. We did not include patients who had to be treated with spinal cord sparing plans after prior full dose irradiation of the spinal cord. While some patients had in fact received prior radiotherapy, this previous therapy did not require sparing of spinal cord to meet dose constraints . Details on irradiation sites are summarized in Table 1. Further information such as primary tumor site of each patient is provided in Additional file 1: Table S1.
Automated VMAT plan generation with Erasmus-iCycle (auto)
Target volumes were identical to those used in the manual planning approach. General principles of iCycle plan generation are described above. Further details on auto-planning with the Erasmus-iCycle/Monaco system have been previously provided by Voet et al. . A site-specific wish-list for our patient cohort is shown in Additional file 2: Table S2. Apart from objectives for the PTV, the kidneys and the lungs, shells around the PTV are used to steer on conformality. Two cases were included which required beam restrictions, both because patients were unable to lift their arms.
Plan quality and statistical analysis
Plan quality was estimated by calculating the dose to 98% (D98%) and to 2% (D2%) of the PTV. Both parameters were used to compute the homogeneity index (HI, ) which was defined as the ratio between the difference of D2% and D98% and the median dose (D50%), i.e.,:
In addition, we calculated the volume receiving 95% of the prescribed dose (V95% i.e., V28.5Gy, or V38Gy, depending if the prescribed dose was 30 Gy or 40 Gy), and based on V95% a conformity index (CI, reviewed in ) was calculated as follows:
An optimal CI would be considered to be 1; however, values can be above or below 1. To consider both over- and underdosage for statistical comparison of CI, we calculated the difference to 1 as follows:
For OARs, we calculated the mean dose (Dmean), and in case of the spinal cord, the maximum dose to any hotspot (Dmax). All parameters were compared directly between manually generated plans and auto plans. Differences between variables were computed using the two-sample paired Wilcoxon test. Statistical significance was defined as p ≤ 0.05 (two-sided testing). Clinical acceptability of the plans was evaluated by an expert radiation oncologist (FL).
All automatically created plans were clinically deliverable and acceptable. Figures 1 and 2 show dose distributions, and dose-volume histogram comparisons for a kidney-level target, and a lung-level target, respectively.
PTV coverage, defined as V38Gy or V28.5Gy was higher in auto plans when compared to manual plans when the whole patient population was analyzed, with an average difference of 1.47% (SD = 3.17%, p = 0.008, n = 42), see Fig. 3. When subgroups were analyzed, the difference was only significant in lung-level plans (p = 0.004, n = 20) but not at kidney levels (p = 0.110, n = 25). This pattern occurred irrespective whether 3 patients with both kidneys and lungs as OARs in one plan were added to the lung or to the kidney group or to both. Dose parameters were as follows: D98% was numerically higher in auto plans when compared to manual plans but this was not statistically significant in the whole population (p = 0.479), nor in any subgroup (p > 0.05 for both lung and kidney levels). D2% was significantly higher in manual plans over the whole patient population (p < 0.001), however the difference was clinically irrelevant (mean: 33.4 Gy vs. 33.2 Gy; median 31.1 Gy vs. 31.0 Gy). Subgroup analysis showed that D2% was significantly lower in auto plans at kidney level (p = 0.003), however the difference was marginal and not clinically relevant. In lung-level plans, D2% did not differ significantly between auto and manual plans.
HI as defined above was not significantly different between auto and manual plans in the whole group or in any subgroup (p > 0.05 for all comparisons). Considering the whole population, CI was slightly but statistically significant higher in auto plans than in manual plans: 1.014 ± 0.12 (1 SD) vs. 1.044 ± 0.11. For targets at the kidney level, the mean CIs for manual plans and auto plans were 0.99 ± 0.071 and 1.022 ± 0.067, respectively, and for lung-level targets 1.036 ± 0.157 and 1.065 ± 0.144. The differences between auto and manual plans can be explained by a slightly more pronounced underdosage in manual plans (minimum CI in manual plans: 0.88 [kidney level] and 0.81 [lung level]; minimum CI in auto plans: 0.92, and 0.95, respectively). Cdiff did not differ significantly between auto and manual plans.
Organs at risk
In the pooled analysis of all plans, mean dose to OARs (kidneys and lungs pooled) was significantly lower in auto plans when compared to manual plans on both the right side (p = 0.001, n = 45, relative difference -10.9% [mean value auto vs. manual]) and on the left side (p < 0.001, n = 45, -19.5%). Subgroup analyses in patients with kidneys as OARs, including 3 patients with lungs and kidneys, showed that auto plans were associated with lower dose to both kidneys (p < 0.001, n = 25, -18.6% right kidney, and p < 0.001, n = 25, -23.1%, left kidney), see Fig. 4a for details. In lung-level plans, again including 3 plans in which both lungs and kidneys were included as OARs, only the left lung showed a significantly lower mean dose when auto-planning was used (p = 0.004, n = 20, -13.7%). Mean dose to the right lung did not differ significantly between manual and auto plans (p = 0.827, n = 20, -1.3%), see Fig. 4b for details. There was no significant difference between the maximum doses (Dmax) to the spinal cord between auto and manual groups in the whole population or in any subgroup. However, one patient in the auto group had a hotspot in the spinal cord (>110% of prescribed dose: 33.2 Gy, D50%: 30 Gy).
Estimated treatment time and required monitor units
Auto plans required slightly more (estimated) treatment time in the whole population (p < 0.001, relative difference 7%) and in the subgroups of kidney-level (p < 0.001, relative difference 9.7%), and lung-level plans (p = 0.046, relative difference 3.6%) as no hard constraint was placed on treatment time. This was a consequence of auto plans requiring more monitor units in the whole population and in any subgroup (p-values and relative differences as follows: whole population, p < 0.001, 17.9%; kidney levels: p < 0.001, 15.4%; lung levels: p = 0.001, 20.6%). All target volume information, monitor units, OAR dose values, and treatment times for each patient are also shown in Additional file 1: Table S1.
Though highly conformal radiotherapy is not mandatory in the treatment of painful bone metastases , intensity modulation with static beams or VMAT allows to deliver a highly focused dose distribution in every clinical situation and may be beneficial for patients with oligometastatic disease  who require dose escalation or patients who had received prior courses of radiotherapy limiting doses to OARs in case of reirradiation [18, 19]. Highly conformal radiotherapy is, however, not always applied when appropriate  for its perceived resource intensity, both in personnel and planning time. Unexpectedly, the reason for not performing advanced radiotherapy techniques despite a clear indication seems to be a shortage of sufficiently trained personnel rather than machine shortage [3, 4, 20]. AlDuhaiby et al. described in a Canadian survey in 2012 that limitations in CT simulator or linear accelerator configuration impaired IMRT implementation in only 10% of cases . The most relevant factors hindering IMRT implementation was the need to train existent treatment planners (50%) as well as the need to hire more planners (30%) . Personal shortage is even more relevant in VMAT settings as these are more time consuming as compared to fixed field IMRT ; “increased planning time required for generating the VMAT plans” has been described by Rao et al. as the main disadvantage of VMAT . Therefore, there is a need for tools that reduce planning time and ease the process of developing a radiotherapy plan without compromising plan quality.
Various auto-planning strategies are currently being developed , these include approaches for beam angle selection (e.g., ), and integrated beam weight optimization algorithms. The latter can be based on global optimization , or on sequential beam selection . As discussed above, multicriterial optimization systems have been proposed [11, 12], that enable the user to select the plan out of different sets of beam angles and intensity profiles after optimization (a posteriori setting). Template based auto-planning solutions are already commercially available and have been applied to various target paradigms such as head and neck , breast  and esophageal targets .
Erasmus-iCycle plans, which provide a template-free approach have been successfully applied to head and neck , prostate , and cervical cancer  target volumes to create clinically deliverable plans.
Our results show that fully automated VMAT treatment planning with the Erasmus-iCycle/Monaco system for spinal metastases was non-inferior to conventional treatment planning by expert dosimetrists or medical physicists. Automatically generated plans with this system outperformed manual plans in terms of sparing of OARs in 3 of 4 predefined organs and also PTV coverage was favorable when compared to manual plans. We observed one hotspot in an iCycle plan in the spinal cord (33.2 Gy, D50%: 30 Gy). This problem can be addressed by modifying the wish list to include a maximum dose cost function on the spinal cord. In our analysis, however, we applied a quadratic overdose cost function which lead to maximum doses that were comparable to maximum doses observed in manual plans in all but this patient and allowed for more degrees of freedom regarding dose reduction to other OAR. If the maximum dose has to be strictly controlled, e.g., in case of reirradiation, a maximum dose cost function should be applied. Otherwise, the situation would be evaluated during the manual approval process which might prompt occasional replanning with a maximum dose cost function. Monitor units and estimated treatment time were higher in auto plans; however, this difference was not considered clinically relevant in the treated patient population and resulted in a negligible prolongation of treatment times. The difference in treatment times was expected as the aim of this study was to generate high-quality plans with clinically acceptable delivery times. Monitor units and delivery times could be reduced in the auto-plans if required. In contrast to a study in a more complex setting such as head and neck cancer , the observed dosimetric advantages of Erasmus-iCycle for treatment of spinal metastases were on average of low clinical relevance although advantages across the population were also observed.
External audits have shown that experienced treatment centers may yield superior IMRT plans . Therefore, as postulated most recently by Fogliata et al. , quality improvements of automated planning vs. manual plan generation may even be more pronounced in target paradigms that are not frequently treated in a particular center.
Clinical target volume
Eastern Cooperative Oncology Group
Intensity modulated radiation therapy
Karnofsky performance status
Organ at risk
Planning target volume
Stereotactic body radiation therapy
Treatment planning system
Volumetric modulated arc therapy
Tiwana MS, Barnes M, Kiraly A, Olson RA. Utilization of palliative radiotherapy for bone metastases near end of life in a population-based cohort. BMC Palliat Care. 2016;15:2.
Stieler F, Wolff D, Bauer L, Wertz HJ, Wenz F, Lohr F. Reirradiation of spinal column metastases: comparison of several treatment techniques and dosimetric validation for the use of VMAT. Strahlenther Onkol. 2011;187:406–15.
Mayles WP, Radiotherapy Development B. Survey of the availability and use of advanced radiotherapy technology in the UK. Clin Oncol (R Coll Radiol). 2010;22:636–42.
Shikama N, Tsujino K, Nakamura K, Ishikura S. Survey of advanced radiation technologies used at designated cancer care hospitals in Japan. Jpn J Clin Oncol. 2014;44:72–7.
Voet PW, Dirkx ML, Breedveld S, Al-Mamgani A, Incrocci L, Heijmen BJ. Fully automated volumetric modulated arc therapy plan generation for prostate cancer patients. Int J Radiat Oncol Biol Phys. 2014;88:1175–9.
Breedveld S, Storchi PR, Keijzer M, Heemink AW, Heijmen BJ. A novel approach to multi-criteria inverse planning for IMRT. Phys Med Biol. 2007;52:6339–53.
Breedveld S, Storchi PR, Heijmen BJ. The equivalence of multi-criteria methods for radiotherapy plan optimization. Phys Med Biol. 2009;54:7199–209.
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.
Rossi L, Breedveld S, Heijmen BJ, Voet PW, Lanconelli N, Aluwini S. On the beam direction search space in computerized non-coplanar beam angle optimization for IMRT-prostate SBRT. Phys Med Biol. 2012;57:5441–58.
Voet PW, Breedveld S, Dirkx ML, Levendag PC, Heijmen BJ. Integrated multicriterial optimization of beam angles and intensity profiles for coplanar and noncoplanar head and neck IMRT and implications for VMAT. Med Phys. 2012;39:4858–65.
Schreibmann E, Lahanas M, Xing L, Baltas D. Multiobjective evolutionary optimization of the number of beams, their orientations and weights for intensity-modulated radiation therapy. Phys Med Biol. 2004;49:747–70.
Craft D, Monz M. Simultaneous navigation of multiple Pareto surfaces, with an application to multicriteria IMRT planning with multiple beam angle configurations. Med Phys. 2010;37:736–41.
Sharfo AW, Voet PW, Breedveld S, Mens JW, Hoogeman MS, Heijmen BJ. Comparison of VMAT and IMRT strategies for cervical cancer patients using automated planning. Radiother Oncol. 2015;114:395–401.
Marks LB, Yorke ED, Jackson A, Ten Haken RK, Constine LS, Eisbruch A, Bentzen SM, Nam J, Deasy JO. Use of normal tissue complication probability models in the clinic. Int J Radiat Oncol Biol Phys. 2010;76:S10–9.
Ohtakara K, Hayashi S, Hoshi H. The relation between various conformity indices and the influence of the target coverage difference in prescription isodose surface on these values in intracranial stereotactic radiosurgery. Br J Radiol. 2012;85:e223–8.
Lutz S, Berk L, Chang E, Chow E, Hahn C, Hoskin P, Howell D, Konski A, Kachnic L, Lo S, et al. Palliative radiotherapy for bone metastases: an ASTRO evidence-based guideline. Int J Radiat Oncol Biol Phys. 2011;79:965–76.
Corbin KS, Hellman S, Weichselbaum RR. Extracranial oligometastases: a subset of metastases curable with stereotactic radiotherapy. J Clin Oncol. 2013;31:1384–90.
Kirkpatrick JP, van der Kogel AJ, Schultheiss TE. Radiation dose-volume effects in the spinal cord. Int J Radiat Oncol Biol Phys. 2010;76:S42–9.
Damast S, Wright J, Bilsky M, Hsu M, Zhang Z, Lovelock M, Cox B, Zatcky J, Yamada Y. Impact of dose on local failure rates after image-guided reirradiation of recurrent paraspinal metastases. Int J Radiat Oncol Biol Phys. 2011;81:819–26.
AlDuhaiby EZ, Breen S, Bissonnette JP, Sharpe M, Mayhew L, Tyldesley S, Wilke DR, Hodgson DC. A national survey of the availability of intensity-modulated radiation therapy and stereotactic radiosurgery in Canada. Radiat Oncol. 2012;7:18.
Oliver M, Ansbacher W, Beckham WA. Comparing planning time, delivery time and plan quality for IMRT, RapidArc and Tomotherapy. J Appl Clin Med Phys. 2009;10:3068.
Rao M, Yang W, Chen F, Sheng K, Ye J, Mehta V, Shepard D, Cao D. Comparison of Elekta VMAT with helical tomotherapy and fixed field IMRT: plan quality, delivery efficiency and accuracy. Med Phys. 2010;37:1350–9.
Potrebko PS, McCurdy BM, Butler JB, El-Gubtan AS. Improving intensity-modulated radiation therapy using the anatomic beam orientation optimization algorithm. Med Phys. 2008;35:2170–9.
Lee EK, Fox T, Crocker I. Simultaneous beam geometry and intensity map optimization in intensity-modulated radiation therapy. Int J Radiat Oncol Biol Phys. 2006;64:301–20.
de Pooter JA, Mendez Romero A, Wunderink W, Storchi PR, Heijmen BJ. Automated non-coplanar beam direction optimization improves IMRT in SBRT of liver metastasis. Radiother Oncol. 2008;88:376–81.
Hazell I, Bzdusek K, Kumar P, Hansen CR, Bertelsen A, Eriksen JG, Johansen J, Brink C. Automatic planning of head and neck treatment plans. J Appl Clin Med Phys. 2016;17:5901.
Fogliata A, Nicolini G, Bourgier C, Clivio A, De Rose F, Fenoglietto P, Lobefalo F, Mancosu P, Tomatis S, Vanetti E, et al. Performance of a knowledge-based model for optimization of volumetric modulated Arc therapy plans for single and bilateral breast irradiation. PLoS One. 2015;10:e0145137.
Fogliata A, Nicolini G, Clivio A, Vanetti E, Laksar S, Tozzi A, Scorsetti M, Cozzi L. A broad scope knowledge based model for optimization of VMAT in esophageal cancer: validation and assessment of plan quality among different treatment centers. Radiat Oncol. 2015;10:220.
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.
Nelms BE, Robinson G, Markham J, Velasco K, Boyd S, Narayan S, Wheeler J, Sobczak ML. Variation in external beam treatment plan quality: an inter-institutional study of planners and planning systems. Pract Radiat Oncol. 2012;2:296–305.
Availability of data and materials
All data generated or analyzed during the current study are available from the corresponding author on reasonable request.
DB analyzed patient data, wrote the article and performed statistical analysis; AWMS performed iCycle planning procedures and statistical analysis; BJMH, PWJV, and SB provided assistance with study design and iCycle procedures. FW and FL reviewed PTVs and provided support with study design. FS reviewed irradiation plans and performed manual planning procedures. All authors read and approved the final manuscript.
DB reports personal fees from Siemens AG, personal fees from NB Capital ApS, outside the submitted work; AWMS reports grants from Elekta AB, Stockholm, Sweden, during the conduct of the study; grants from Accuray Inc, Sunnyvale, USA, outside the submitted work; BJMH reports grants from Elekta AB, Stockholm, Sweden, during the conduct of the study; grants from Accuray Inc, Sunnyvale, USA, outside the submitted work; PWJV reports other from COS Research Clinical Liaison, outside the submitted work; and Working as an employee at Elekta B.V., The Netherlands; SB reports grants from Elekta AB, Stockholm, Sweden, during the conduct of the study; grants from Accuray Inc, Sunnyvale, USA, outside the submitted work; FW reports grants and personal fees from Elekta, during the conduct of the study; FL reports grants and personal fees from Elekta AB, grants and personal fees from IBA, personal fees from C-RAD, during the conduct of the study. FS has nothing to disclose.
Consent for publication
All patient data were anonymized. Approval for publication of anonymized data and images were obtained from all patients.
Ethics approval and consent to participate
Informed consent was obtained from all patients for anonymized processing of their clinical data. The study was approved by the ethics committee of Heidelberg University, Medical Faculty Mannheim (2016-806R-MA).