Dosimetric robustness against setup errors in charged particle radiotherapy of skull base tumors
Radiation Oncology volume 9, Article number: 279 (2014)
It is expected that physical dose deposition properties render charged particle dose distributions sensitive to targeting uncertainties. Purpose of this work was to investigate the robustness of scanned-beam particle therapy plans against setup errors for different optimization modalities, beam setups and ion species.
Material and methods
For 15 patients with skull base tumors, localized in regions of severe tissue density heterogeneity, scanned lateral-opposed-beam treatment plans were prepared with the treatment planning system TRiP98, employing different optimization settings (single- and multiple-field modulation) and ion species (carbon ions and protons). For 10 of the patients, additional plans were prepared with individually selected beam setups, aiming at avoiding severe tissue heterogeneities. Subsequently, multiple rigid positioning errors of magnitude 1–2 mm (i.e. within planning target expansion) were simulated by introducing a shift of the irradiation fields with respect to the computed tomography (CT) data and recomputing the plans.
In presence of shifts, in carbon ion plans using a lateral-opposed beam setup and fulfilling clinical healthy tissue dose constraints, the median reduction in CTV V95% was up to 0.7 percentage points (pp) and 3.5 pp, for shifts of magnitude 1 mm and 2 mm respectively, however, in individual cases, the reduction reached 5.1 pp and 9.7 pp. In the corresponding proton plans similar median CTV V95% reductions of up to 0.9 pp (1 mm error) and 3.4 pp (2 mm error) were observed, with respective individual-case reductions of at most 3.2 pp and 11.7 pp. Unconstrained plans offered slightly higher coverage values, while no relevant differences were observed between different field modulation methods. Individually selected beam setups had a visible dosimetric advantage over lateral-opposed beams, for both particle species. While carbons provided more conformal plans and generally more advantageous absolute dose values, in presence of setup errors, protons showed greater dosimetric stability, in most of the investigated scenarios.
Residual patient setup errors may lead to substantial dose perturbation in scanned-beam particle therapy of skull base tumors, which cannot be dealt with by planning target expansion alone. Choice of irradiation directions avoiding extreme density heterogeneities can improve plan stability against such delivery-time uncertainties.
The so-called inverse depth-dose profile and the sharp lateral penumbra make charged particles, like protons and carbon ions, apt for high-precision dose-escalated radiotherapy, e.g. in the treatment of intracranial tumors, due to the direct vicinity of critical structures ,. Carbon ions offer additionally an enhanced relative biological effectiveness (RBE) .
Yet, based on physical properties of dose deposition by charged particle beams, it is expected that these advantages come to the price of greater sensitivity to delivery-time targeting uncertainties as both, the geometrical miss due to sharp lateral gradients and the displacement of the density heterogeneity interfaces within the treatment field, may result in delivered dose deterioration, which may not be completely dealt with by planning target expansion . For scanned-beam intensity-modulated treatments, which can produce steep dose gradients and often rely on them to achieve their superior conformity, the problem of dosimetric robustness is likely of greater concern . In this respect, a role may also be played by the physical dose modulation typically present in biologically-optimized carbon ion treatment plans, on account of the variable RBE, depending on multiple factors, including local dose .
The importance of the selection of an optimal plan, not only in terms of clinical dosimetric quality (target dose coverage versus normal tissue sparing), but also taking into account robustness, i.e. degree of sensitivity of the plan to the uncertainties involved in the treatment process, has been recognized by the ICRU  and emphasized, more recently, by Cantone et al.. One method to assess the relevance of the above-mentioned dosimetric effects and to support clinical decision making is recomputation of treatment plans in geometries simulating the expected errors. Jäkel et al. reported this procedure as part of patient-specific quality assurance for scanned carbon ion therapy, while a method for the evaluation of such tests was demonstrated by Lomax .
While the dosimetric robustness of particle therapy treatment plans has been extensively studied for other indications, e.g. prostate tumors ,, the published literature on cranial tumors is sparse. Recently, Hopfgartner et al. investigated the robustness of proton treatment plans with respect to interfractional setup uncertainties, while the only previous work addressing this issue for carbon ions, a conference communication by Ellerbrock et al., focused primarily on the evaluation of the multiple-field optimization approach, i.e. full intensity modulation (typically referred to as IMPT), as opposed to the single-field uniform-dose (SFUD) approach.
In this study we apply extensive plan recomputation, on a cohort of patients with skull base tumors, in order to systematically assess dosimetric robustness of scanned beam charged particle treatment plans in presence of setup errors and critical heterogeneities for (1) various plan optimization modalities, (2) different beam setups and (3) two particle species (carbon ions and protons).
Material and methods
The data of 15 patients with skull base tumors (P1–P15) treated at our institution with intensity-modulated radiotherapy (IMRT) (P1, P13) or with stereotactic radiotherapy (SRT) (remaining) were chosen, with tumor localization typical for particle therapy indications (chordoma, chondrosarcoma, adenocystic carcinoma), in areas of high density heterogeneity and close to dose limiting structures.
Planning computed tomographies (CT) were acquired, without contrast medium, with in-slice resolution and slice thickness of 0.98 mm/pixel and 3 mm (P1, P13), 0.59 mm/pixel and 3 mm (P2, P11–P12) and 0.59 mm/pixel and 2.5 mm (remaining). The original clinical target volume (CTV) and organs at risk (OAR) contours, retained for the study, were delineated in the Pinnacle3 (Philips Healthcare, Best, The Netherlands) treatment planning system (TPS) for IMRT cases and in the Virtuos (Stryker-Leibinger, Freiburg, Germany) system for SRT cases, excluding unnecessary air cavities from target volumes. The median (range) CTV volume was 45.2 (15.6–90.7) cm3. The planning target volume (PTV) was defined, similarly to Orecchia et al., by applying a margin of 2 mm in anteroposterior (AP) and mediolateral (LR) direction and of 2.5–3 mm in superoinferior (SI) direction, depending on CT slice thickness. The CTV-to-PTV expansion was performed on a single TPS to avoid potential inter-system algorithm variability. A contour separating the skin surface from the immobilization mask was added and densities outside it overridden with air, to exclude mask-specific heterogeneities from the investigation.
Following our institutional ethics procedures, all patients were informed of potential retrospective research use of their data and given the choice to opt out. The data presented in this manuscript do not reveal, either directly or indirectly, the patients’ identities.
Carbon ion and proton treatment plans for raster-scanning delivery  were prepared using the TRiP98 planning system (GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany) ,. The conversion program dcm2trip, developed at our institution and distributed with TRiP98, was used to convert the DICOM data (CT images and structure sets) to the VOXELPLAN format used by TRiP98 .
The local effect model LEM I was used for calculation of the RBE of carbon ions (α/β=2 Gy, threshold dose Dt=30 Gy, nuclear radius r=4.5μ m), while a fixed value of 1.1 was employed for proton RBE.
The prescription dose was set to 60–63 Gy(RBE) delivered in 20–21 fractions for carbon ion plans and 64–74 Gy(RBE) in 32–37 fractions for proton plans similarly to Nikoghosyan et al., and to Ares et al.. The planning objective was to deliver at least 95% of the prescription dose to 95% of the PTV (V95%≥95%) and at least 95% of the prescription dose to 98% of the CTV volume (V95%≥98%), if possible without compromising organs at risk.
Also OAR dose constraints were derived from the protocols adopted in the above-mentioned studies. A maximum dose (Dmax) constraint at 54 Gy(RBE) was defined for optic nerves, optic chiasm and brainstem, however, the brainstem surface (1% of volume) abutting the PTV was allowed to receive a dose of up to 60 Gy(RBE) ,. In proton plans, the corresponding limits were 60 Gy(RBE) and 63 Gy(RBE) .
For each patient two carbon ion treatment plans and a proton plan were prepared using two lateral-opposed beams (LR), reflecting the geometry typically available at fixed-nozzle carbon ion facilities.
The two carbon ion plans had, respectively, OAR constraints requiring per-field modulation (single-field optimization, denoted as sc) and full plan modulation (multiple field optimization, denoted as mc). The latter modality was also used in proton planning. Additionally, for comparison, plans without OAR constraints (unconstrained mode, denoted as uc) were prepared for both ion species.
Finally, for selected patients (n=10), carbon ion and proton plans using full modulation (mc) were prepared using two beams with individually selected directions. By visually analyzing potential beam directions, those were selected that appeared to minimize the presence of density interfaces along the entrance channels, that could introduce range variations in case of setup errors. These beam setups, referred to as robust (ROB) in the remainder of the text, employ a range of isocentric couch rotations (yaw) towards cranial direction, combined, in the majority of cases, with beam inclination about the longitudinal axis. Such gantry-rotation-like setups are enabled, at our facility, by combined use of 45-degree-inclined beam lines, 6-degree-of-freedom robotic couches (enabling patient roll and pitch up to 15 degrees with proper immobilization) and, when required, individualized solutions (immobilization with an advantageous head rotation or treatment in prone position).
At synchrotron based facilities, the available pencil beam transversal widths (spot sizes) in air at the treatment room isocenter are a function of the beam energy . For carbon ion plans pencil beam full-width at half-maximum varied from 7.5 mm to 5.0 mm, as specified for the required energy range in the TPS synchrotron library, compatible with the GSI facility during the German ion beam therapy pilot project and representative of modern combined-beam facilities, like the Particle Therapy Center in Marburg. For proton plans spot sizes between 11.0 mm and 5.5 mm were enabled by treatment planning and delivery with shorter nozzle-to-patient distance than for carbon ions . These planning settings enabled the highest available comparability of treatment plans, as they resulted in proton spot sizes, measured at patient surface, as close as possible to the corresponding ones of carbon ions.
In the irradiation raster, the transversal distance between neighboring spots (pitch) was set to 2 mm for carbon ion plans and 3 mm for proton plans. A tolerance of 0.4–1.0 times the spot size was set, allowing the TPS to place additional raster spots outside the target volume projection, transversely to the beam, to ensure PTV coverage without the need for high fluence spots at the target border. For carbon ion plans, a ripple filter was used to broaden the narrow pristine Bragg peaks and enable an in-depth peak positioning step of 3 mm . For proton planning an in-depth peak distance of 2 mm was used. In the plans of one patient (P03) a bolus of 20 mm of water-equivalent material was used to ensure coverage of a tumor section shallower than the penetration depth of the lowest energy available. The dose grid resolution was the same of the planning CT.
Dosimetric robustness tests
Rigid positioning setup errors were simulated by introducing shifts of the irradiation fields with respect to the CT data and recomputing the dose without plan reoptimization.
All plans were tested against 52 shifts, divided in two groups respectively of 3D magnitude 1 mm and 2 mm. For each magnitude the group represents a complete enumeration of shifts with components, in the patient’s coordinate system (LR, AP, SI), either zero or equal in absolute value. In the figures, individual shifts are expressed through their magnitude and a triplet, with the + and - signs indicating the direction of the (equal) non-zero components. The following formula generates the set S of all 52 tested shifts, expressed in mm, in the patient’s coordinate system:
The treatment plans were compared in terms of dose distribution and dose-volume histograms (DVH). For a quantitative assessment of the dose distribution in the PTV and CTV the 95% isodose coverage index (V95%) and the homogeneity index (HI), calculated as the difference (D2% – D98%), were used , along with the PTV conformity index (CI) . To evaluate OAR dose, the maximum dose (Dmax) was reported for the optic structures and the near-maximum dose at 1% volume (Dnear-max) for the brainstem. Visual and quantitative dose distribution analysis was entirely performed with trip2png, developed at our institution.
To analyze dosimetric plan robustness, the differences of the index values between the recomputed (shifted) and the initial plan were employed and denoted in the following as Δ V95%, Δ HI and Δ D(near-)max respectively.
For selected comparisons and indices, statistical significance was assessed through a non-parametric, paired-sample sign test with a significance level of 0.05, using the R environment for statistical computing .
Initial dosimetric quality of treatment plans
The mean (range) prescribed dose was 61.8 (60–63) Gy(RBE) for carbon ion plans and 66.9 (64–74) Gy(RBE) for proton plans. The initial (no shift) plans showed comparable target coverage for all planning approaches. The median PTV and CTV V95%, PTV and CTV HI, PTV CI and OAR D(near-)max values for all planning scenarios are shown in Table 1. In all carbon ion treatment plans, the planning objective of PTV V95%≥95.0% was fulfilled, except in one of the sc LR plans, where it was compromised to 93.1% in order to meet OAR sparing constraints. For the same reason, for protons the planning objective (PTV V95%≥95.0%) had to be compromised twice in mc LR plans (to 94.4% and 90.4%) and once in mc ROB plans (to 93.1%).
Figure 1 presents carbon ion and proton axial dose distributions of an example patient, as planned and as resulting from exposure to setup errors of magnitude 1 mm and 2 mm, for various planning scenarios: a plan with multiple-field optimization and lateral-opposed beam directions (mc LR) and a plan with multiple-field optimization and individually chosen, robust beam directions (mc ROB). The lateral-opposed-beam dose distributions for both particle species show cold spots reaching into the CTV volume, despite the 2 mm CTV-to-PTV margin, a typical consequence of particle beam displacement with respect to a strong tissue density interface. The underdosage is not present in dose distributions from robust beam setups.
The corresponding dose-volume histograms of CTV and chiasm are presented in Figure 2, together with the DVHs of all other simulated shifts, illustrating the range of target coverage reduction and OAR involvement variability.
The changes of carbon ion and proton CTV V95% and HI induced by selected shifts, across the studied patient cohort, are shown in Figure 3. Since a shift along the beam does not introduce any density variation, i.e. any dose perturbation, for each pair of setup errors differing only in the sign of the LR component, only one is shown (positive LR component), for the sake of figure readability.
In general, setup errors caused a reduction of target coverage and homogeneity in the dose distributions. For lateral-opposed carbon ion plans making use of beam modulation (mc LR), the median reduction in CTV V95% was up to 0.7 pp and 3.5 pp, for the worst-case shifts of magnitude 1 mm and 2 mm respectively. In individual cases, notably stronger coverage deterioration has been observed, with CTV V95% decreasing by up to 5.1 pp and 9.7 pp respectively. The corresponding values for the proton plans of the same type were comparable, with median CTV V95% reduction of at most 0.9 pp (1 mm error) and 3.4 pp (2 mm error) and respective individual case reduction of at most 3.2 pp and 11.7 pp.
Comparison of optimization modalities
Carbon ion plans
Carbon ion plans introducing OAR sparing constraints through full intensity-modulation (mc LR) exhibited, in presence of setup errors, higher losses in target coverage, in comparison to the unconstrained counterparts. As visible in Figure 3, these differences in CTV Δ V95% were only slight in median value (within 1 pp), for all simulated shifts, while in individual cases they reached 3.7 pp. This reduction of robustness, induced by the modulation of the irradiation fields, is also reflected in the histograms in Figure 4, with the predominance of smaller CTV V95% changes in the unconstrained plans (about 70% vs. 63% of displacements resulting in Δ V95% smaller than 1 pp in absolute value). Additionally, in about half of the simulated shifts, the reduction of CTV V95%, induced by the additional field modulation, was found to be statistically significant, although the generally low differences in median value suggest a lack of clinical relevance.
The effect of setup errors on target dose homogeneity appeared to be virtually unaffected by the difference in field modulation. For all simulated shifts, median CTV Δ HI of mc LR plans oscillated about the corresponding uc LR values without a clear trend and always below 1 pp in absolute value. Also the CTV HI histograms in Figure 4 show no distinct effect.
Increased incidence of larger Dmax changes in critical structures was also observed for constrained plans (Figure 4). In particular, for the ipsilateral optic nerve in 16% more of the cases the Dmax change exceeded 2 Gy(RBE) in absolute value. For the contralateral optic nerve this difference was 6% and for optic chiasm 10%. The increase in beam modulation almost did not affect Dnear-max stability of the brainstem: the number of observations with an absolute deviation larger than 2 Gy(RBE) varied by less than 4% of the total, between the two optimization approaches. In the majority of cases under investigation, the brainstem was abutting the PTV, with its surface located in the strong dose gradient region for all planning modalities, rendering Dnear-max subject to significant changes for all plans.
Regarding sc LR treatment plans, even smaller differences between them and the corresponding mc LR were observed. Throughout all simulated shifts the variation in CTV Δ V95% was at most 0.4 pp in median value and 1.3 pp in single cases. Also the smaller deviations in the histograms (Figure 4) suggest that in practice no difference in robustness against setup errors exists between the two constrained planning approaches, at least at the levels of field modulation required for the indications under investigation.
For proton plans, the effect of modulation on CTV coverage, in presence of setup errors, was less pronounced and did not present a clear trend (Figure 3). Also for protons, median CTV V95% was slightly better preserved in unconstrained plans, throughout the simulated shifts, but this was not always true in individual cases. Together with the small differences in median value, a paired-sample sign test reported, for almost all investigated shifts, no statistically significant difference between the two proton planning approaches. The similar behavior of mu LR and mc LR plans is also confirmed by the respective CTV Δ V95% histograms in Figure 4.
Median values throughout all simulated setup errors seem to suggest a better HI preservation in the case of constrained plans (Figure 3). This is also somewhat visible in the histogram in Figure 4, where smaller CTV HI changes appear to be a little more common (4% of cases more) in constrained than in unconstrained plans.
With the introduction of modulation, also an incidence of larger D(near-)max changes for critical structures was observed, which was mostly pronounced for the ipsilateral optic nerve (8% cases more exceeding 2 Gy(RBE) in absolute value) and less marked for the other optic nerve, for the optic chiasm and for the brainstem (respectively 4%, 5% and 5%).
Comparison of beam setups
Carbon ion plans
Carbon ion treatment plans with beam setups selected to avoid severe density inhomogeneities (mc ROB) exhibited, in presence of setup errors, greater dosimetric stability than the corresponding lateral-opposed plans (mc LR), in terms of preservation of both target coverage and dose homogeneity. Among all simulated shifts, mc ROB plans were affected by CTV V95% loss of at most 2.5 pp in median value, reaching 6.3 pp in individual cases, which represents a clear improvement in comparison with the mc LR plans of the same patient cohort, with respective values of 4 pp and 9.7 pp. Similarly, for the CTV HI, mc ROB plans presented better median value preservation (by up to 1.6 pp) and particularly strong worst case reduction, almost halving the difference.
While a direct shift-by-shift comparison between the plans is not applicable, as the influence of each setup error is specific to irradiation direction, the advantage of robust over lateral-opposed beam setups is visible in the histograms in Figure 5 as strong dominance of smaller changes (below 1 pp in absolute value) in CTV V95% (16% of cases more) and CTV HI (18% of cases more) for mc ROB plans and in the prolonged tails toward larger variations for mc LR plans.
Also organs at risk benefitted from the greater dosimetric stability, afforded by robust beam selection, in presence of setup errors. Reduced incidence of larger Dmax changes, in comparison with the lateral-opposed plans, was observed for most of the critical structures: the index deterioration remained below 2 Gy(RBE) in absolute value in 25% of the shifts more for the chiasm, 11% more for the contralateral optic nerve and 7% more for the ipsilateral optic nerve. Only for the brainstem, which for nearly all patients was in contact with the target dose gradient region, no considerable improvement in Dnear-max stability was observed.
In presence of setup errors, proton mc ROB treatment plans achieved, similarly to the equivalent carbon ion ones, better preservation of target coverage and dose homogeneity in comparison to plans with a lateral-opposed beam setup. In mc ROB plans, median and worst-case CTV V95% losses were 2.8 pp and 5.8 pp respectively, while in the corresponding mc LR plan they were 3.6 pp and 11.7 pp. HI median and worst-case variations were instead 3.4 pp and 6 pp for mc ROB plans and 3.9 pp and 14.2 pp for mc LR plans. The dosimetric advantage in target dose conservation through robust beam configurations is visible in the V95% and HI deterioration histograms in Figure 5, where mc ROB exhibit a larger number of cases (respectively 7% and 8%) affected by smaller absolute changes.
Proton mc ROB treatment plans also offered improved OAR dose stability. In comparison to the corresponding lateral-opposed plans, in a higher number of cases D(near-)max remained in absolute value within 2 Gy(RBE) for the contralateral optic nerve (16% of cases more), ipsilateral optic nerve (9%), optic chiasm (29%) and brainstem (7%).
Comparison of ion species
The preceding subsections already report data on the investigated scenarios, for both carbon ions and protons. Some results of potential interest are better described through direct comparison of the two ion species, as presented in Figure 6, for mc LR and mc ROB plans, which both represent clinically valid treatment planning approaches.
With a lateral-opposed beam setup protons presented greater dosimetric stability in the target, a result that is consistent with the lower initial conformity of proton plans used in this study (Table 1). Use of robust beam setups reversed this situation, affording greater CTV V95% reproducibility to carbon ions and making the two particle species practically equal in terms of HI preservation. This result shows that applying a planning technique aiming at reducing physical effects of setup errors, like selecting less inhomogeneous beam paths, brought greater benefit to the particle species with the most conformal dose distributions. For both mc LR and mc ROB plans, protons showed greater reproducibility of planned OAR dose, in presence of setup errors, with the greatest advantage for Dnear-max in the brainstem, that is the OAR closest to the target dose gradient (Figure 6).
Finally, it should be emphasized that, while in terms of relative changes, between optimized plans and simulated delivery with setup errors, protons demonstrated greater dosimetric stability (i.e. a larger number of smaller changes) for the majority of investigated indices and plan types, in terms of absolute dose values carbon ion plans used for this study remained superior in terms of both median target coverage and OAR sparing (Table 1).
Systematic plan robustness tests, like those presented in this manuscript, are necessary to understand the magnitude of delivery-time dose uncertainties, in order to optimize treatment planning methodologies and also to ensure the selection of the best treatment for individual patients.
In fact, despite its heavy requirements in terms of time and computing resources, a procedure based on multiple plan recomputations, similar to the one used in our study, was reported by Jäckel et al. as part of the clinical patient-specific quality assurance process, for scanned carbon ion therapy. An alternative measurement-based method was proposed by Albertini et al., to assess plan robustness by recording the effects of setup (or range) uncertainties on a plan delivered to an anthropomorphic phantom, thus excluding the inaccuracies of analytical dose computation algorithms, typically employed in clinical planning systems. While this is certainly a considerable advantage, this approach is of limited applicability to carbon ions, for which only the underlying absorbed dose distribution deterioration can be measured. Furthermore, relying on the standardized phantom anatomy appears more suitable a method to compare the influence of optimization approaches or planning settings for typical situations, approximated by target volume definition in the phantom. Instead, as also shown in our study, significant variability between the patients may be present, stemming from the differences between individual target volumes and anatomical features, as well as their spatial relation, which emphasizes the importance of conducting robustness assessments for specific patients (or well-defined populations) in order to support clinical decision making.
Alternatively, accurate dose computation on real patient geometries can be achieved through Monte Carlo methods, a particularly interesting option for proton therapy, for which a dedicated high-performance implementation was recently demonstrated . For carbon ions, instead, while novel integrated tools can simplify clinical applications ,, the prolonged computational times must be weighed against the expected precision gain, especially in large scale computations. For one of the cases presented in this manuscript, for instance, Monte Carlo was employed on a single treatment plan and a single recomputation, demonstrating slight discrepancies in highly-heterogeneous regions and absorbed-dose DVH differences in the order of 2% , which warrant more extensive future investigations.
A critical point, when investigating the dosimetric effects of setup errors, is selecting the range of errors to simulate. Two approaches have been proposed in the literature, simulating either a small number of extreme displacements, representing the worst-case scenarios, or calculating numerous dose distributions, reflecting the entire possible range (and probability) of displacements. Data obtained in our study suggests that the former may lead to false conclusions, as the worst case is not known a priori and not always strictly correlated to the displacements magnitude. More importantly, the worst-case displacement is not the same for different beam setups, rendering direct comparisons of limited value. An example of this is the comparison of a lateral-opposed beam setup against other beam configurations in presence of mediolateral shifts, by which the former remains unaffected.
Therefore, the study presented in this manuscript was based on massive plan recomputation, with real patient data and a comprehensive, fine-grained set of simulated setup errors, and is, to our best knowledge, the first of this kind investigating multiple aspects, that affect carbon ion and proton dosimetric robustness, in skull base radiotherapy, on a multi-patient cohort.
The only previous study addressing the issue of robustness against setup errors for carbon ion plans, a conference communication by Ellerbrock et al., reported an absolute reduction of the CTV V90% by 0–1 pp for 1 mm displacements and 2–4 pp for 2 mm displacements in one patient. While this work, like ours, tested a broad range of errors, it focused primarily on usability of IMPT and employed plans with systematically varying weight of constraints, rather than clinically valid ones. It should also be noted that, to counterbalance the computation-time explosion of multiple optimizations, only absorbed dose computations were performed, neglecting both the additional field modulation and the potential compensation effects, introduced by variable carbon ion RBE. More recently, Hopfgartner et al. systematically investigated the robustness of proton treatment plans with respect to interfraction setup uncertainties in 7 skull base (and 7 paranasal sinus) tumor cases for 6 extreme orthogonal displacements, reporting median CTV V95% reductions of up to about 3 pp. In contrast to this study, which aimed at quantifying robustness in presence of extreme displacements for simple and more complex beam setups, selected with sole consideration of critical structure sparing, the goals of our work included extensive investigation of the robustness loss due to field modulation as well as its improvement through avoidance of critical heterogeneities in the beam channels.
Our results demonstrate that residual repositioning errors may, in some cases, lead to a significant perturbation of the delivered dose in scanned-beam carbon ion and proton therapy for skull base tumors (by up to 10 pp for carbon ions and up to 12 pp in proton plans, in individual cases), which cannot be completely dealt with by planning target expansion alone. Only rigid setup errors were simulated, as they represent the predominant source of targeting uncertainty in skull base treatments, with magnitudes (1–2 mm) consistent with the precision reported for modern non-invasive cranial immobilization devices and patient positioning systems, in terms of residual setup errors as well as intrafraction motion ,-. No rotational errors were simulated, as in comparison to translational displacements, angular deviations cause smaller effects as established in the preliminary work for this study  and also observed for other indications .
The clinical plans (fulfilling OAR constraints) displayed a slightly reduced robustness, understood as the spectrum of variation of dosimetric indexes induced by presence of simulated setup errors, as compared to plans prepared without OAR constraints, illustrating the isolated influence of the dose gradients within the individual, modulated fields. For individual shifts these differences however seldom reached statistical significance, confirming the observation by Ellerbrock et al., that the generally stronger dose gradients at the field edge have more significant interactions with positioning uncertainties, while internal dose gradients, caused by inhomogeneous dose delivery by individual fields, are weaker and thus have minor influence on the robustness. Additionally, in carbon ion plans, only minimal differences were observed between the two modulation approaches (single- vs. multiple-field optimization). The strength of the planning constraints, which determines the magnitude of field modulation, is likely a factor influencing this result. In our study, typical dose modulation (defined here as the ratio between the maximum OAR dose and the prescription dose) was 0.85. The robustness of treatment plans where the application of stronger OAR constraints might be necessary, e.g. for re-irradiations, should be investigated for each individual case. Further complexity lies in the potentially large degeneracy of the optimization result. It has been demonstrated that dosimetrically equivalent IMPT plans employing different modulation strategies may be characterized by differing degrees of robustness .
One of the most interesting applications of robustness tests against an extensive set of setup errors is the direct comparison, in terms of dosimetric stability, of different beam-setups. In our study, for a subset of patients, improved plan robustness, with respect to both CTV and OAR dose, was demonstrated when using beam setups avoiding strongly heterogeneous tissue. It should be noted that the beam directions for this set of treatment plans were selected with sole consideration of the density homogeneity in the beam channels, i.e. not considering the position of critical structures, which were however spared through plan optimization (in compliance with the planning protocol).
Currently, clinical robustness testing of irradiation setups must be usually carried out by expert personnel, on a case-by-case basis, evaluating possible tradeoffs with planning constraints, like OAR sparing. This underlines the need for automatized tools to support such multi-criteria, optimal beam setup selection. Verification of the entrance channel homogeneity through pattern matching techniques, that mimic human analysis, has been demonstrated a promising approach, which can be easily combined, through an optimization process, with clinical constraints, like patient positioning limitations or OAR involvement ,. Methods have also been proposed to include robustness in the plan optimization process in proton therapy  or in margin design ,, but they suffer from the limitation of being based on assumed distributions of positioning errors and do not have, up to now, widespread clinical application.
Finally, our study presents the first attempt to compare the dosimetric robustness of treatment plans employing different ion species. Between protons and carbon ions, used in our study, differences have been demonstrated in the basic dose deposition properties, that may be relevant in their application as therapeutic radiation types . In our study, we strove to maintain proton and carbon ion plans as comparable as possible, in both technical (e.g. using similar raster spot sizes) and clinical terms (e.g. using comparable planning objectives). Also, for both particle species plans with two beams were used, as typical in carbon ion therapy of the skull base and not uncommon in proton therapy, especially at fixed-nozzle combined facilities ,,. Yet, in order to achieve clinically acceptable plans for both particle species, not all planning settings could be kept strictly equal, which should not be neglected when interpreting our results. For instance the significant beam broadening in matter undergone by protons required a larger target projection expansion, to achieve acceptable PTV coverage, and resulted in a lower initial conformity of the proton plans. While these were comparable, in terms of CTV coverage, to those investigated in a similar study by Hopfgartner et al., it should be noted that, in our study, realistic proton spot sizes attainable at synchrotron-based combined ion therapy facilities were used, which were inferior to the idealized ones used by Hopfgartner et al.. Additional proton beams, as often used when a rotating gantry is available, would have likely improved the initial conformity, but they would have also altered the robustness, as observed e.g. by Hopfgartner et al..
Most notably, from a strictly numerical interpretation of our results, in most of the scenarios considered, in presence of setup errors protons appear to deviate less from the planned values, both in CTV coverage and OAR involvement, with a particular stability advantage for the brainstem, the organ closest to the target and often in direct contact with it. This observation could be of importance in clinical decision-making, e.g. when defining guidelines and thresholds for image-guided radiotherapy (IGRT) protocols, and is likely related to the more pronounced scattering protons undergo, resulting in typically shallower gradients. Another important observation, related to potential clinical applications, is that carbon ion plans, with their higher initial conformity and more advantageous absolute dosimetric values, received greater advantage, in terms of target dose reproducibility, through use of robust beam setups.
Ultimately, purpose of the comparison was not to identify the best particle species, but rather to demonstrate how the differences in the physical properties of ions, and in the clinical applications they enable, may affect our understanding of robustness issues. In this sense, a certainly valid general conclusion is that per-indication studies should be carried out for each ion species and some care should be used when translating results from one particle type to another. This is particularly important in the current context of combined facilities, where protons and carbon ions are used for different indications, and as newer ion species, e.g. oxygen or helium, are being considered for introduction in clinical particle therapy . Similarly, where a radiobiological model is used to compute effective dose, results should not be considered valid across different models, model versions or model input parameters ,.
When assessing the impact of these results on patient treatments, one aspect, crucial for the interpretation of this study, should not be neglected. The application of the total therapeutic dose is typically fractionated, causing delivery uncertainties to be statistically smoothed, provided that they are not systematic. It is therefore advisable that, for a thorough clinical assessment, a similar study be conducted in the future, employing e.g. the results of daily positioning verification data, to assess cumulative dose. Finally, the positioning errors are only one of many potential delivery-time uncertainties that can affect the final dose distribution, like range uncertainties , anatomy changes , beam fluctuations , etc. For novel techniques, like carbon ion scanned-beam radiotherapy, the relevance of such uncertainties and of their interplay remains mostly to be quantified.
Schulz-Ertner D, Tsujii H: Particle radiation therapy using proton and heavier ion beams. J Clin Oncol 2007, 25: 953-964. 10.1200/JCO.2006.09.7816
Orecchia R, Vitolo V, Fiore M. R, Fossati P, Iannalfi A, Vischioni B, Srivastava A, Tuan J, Ciocca M, Molinelli S, Mirandola A, Vilches G, Mairani A, Tagaste B, Riboldi M, Fontana G, Baroni G, Rossi S, Krengli M: Proton beam radiotherapy: report of the first ten patients treated at the “Centro Nazionale di Adroterapia Oncologica (CNAO)” for skull base and spine tumours. Radiol med 2014, 119: 277-282. 10.1007/s11547-013-0345-0
Lomax AJ: Intensity modulated proton therapy and its sensitivity to treatment uncertainties 2: the potential effects of inter-fraction and inter-field motions. Phys Med Biol 2008, 53: 1043-1056. 10.1088/0031-9155/53/4/015
Krämer M, Scholz M: Treatment planning for heavy-ion radiotherapy: calculation and optimization of biologically effective dose. Phys Med Biol 2000, 45: 3319-3330. 10.1088/0031-9155/45/11/314
International Commission on Radiation Units and Measurements (ICRU): ICRU Report 78: Prescribing, Recording, and Reporting Proton-Beam Therapy, Bethesda, MD: ICRU; 2007.
Cantone MC, Ciocca M, Dionisi F, Fossati P, Lorentini S, Krengli M, Molinelli S, Orecchia R, Schwarz M, Veronese I, Vitolo V: Application of failure mode and effects analysis to treatment planning in scanned proton beam radiotherapy. Radiat Oncol 2013, 8: 127. 10.1186/1748-717X-8-127
Jäkel O, Krämer M, Karger CP, Debus J: Treatment planning for heavy ion radiotherapy: clinical implementation and application. Phys Med Biol 2001, 46: 1101-1116. 10.1088/0031-9155/46/4/314
Meyer J, Bluett J, Amos R, Levy L, Choi S, Nguyen QN, Zhu XR, Gillin M, Lee A: Spot scanning proton beam therapy for prostate cancer: treatment planning technique and analysis of consequences of rotational and translational alignment errors. Int J Radiat Oncol Biol Phys 2010, 78: 428-434. 10.1016/j.ijrobp.2009.07.1696
Jelen U, Ammazzalorso F, Chanrion MA, Gräf S, Zink K, Engenhart-Cabillic R, Wittig A: Robustness against interfraction prostate movement in scanned ion beam radiation therapy. Int J Radiat Oncol Biol Phys 2012, 84: 257-262. 10.1016/j.ijrobp.2012.03.058
Hopfgartner J, Stock M, Knäusl B, Georg D: Robustness of IMPT treatment plans with respect to inter-fractional set-up uncertainties: Impact of various beam arrangements for cranial targets. Acta Oncol 2013, 52: 570-579. 10.3109/0284186X.2012.744874
Ellerbrock M, Jäkel O, Krämer M, Nikoghosyan A, Schulz-Ertner D, Karger CP, Ackermann B, Heeg P, Debus J: Clinical implementation of intensity modulated radiotherapy using carbon ions [abstract]. Int J Radiat Oncol Biol Phys 2008,72(Suppl 1):592. 10.1016/j.ijrobp.2008.06.199
Haberer T, Becher W, Schardt D, Kraft G: Magnetic scanning system for heavy ion therapy. Nucl Instrum Methods A 1993, 330: 296-305. 10.1016/0168-9002(93)91335-K
Krämer M, Jäkel O, Haberer T, Kraft G, Schardt D, Weber U: Treatment planning for heavy-ion radiotherapy: physical beam modeling and dose optimization. Phys Med Biol 2000, 45: 3299-3317. 10.1088/0031-9155/45/11/313
TRiP98 Online Documentation [http://bio.gsi.de/DOCS/TRiP98/DOCS/trip98.html] TRiP98 Online Documentation
Nikoghosyan AV, Karapanagiotou-Schenkel I, Münter MW, Jensen AD, Combs SE, Debus J: Randomised trial of proton vs and carbon ion radiation therapy in patients with chordoma of the skull base, clinical phase III study HIT-1-Study. BMC Cancer 2010, 10: 607. 10.1186/1471-2407-10-607
Nikoghosyan AV, Rauch G, Münter MW, Jensen AD, Combs SE, Kieser M, Debus J: Randomised trial of proton vs and carbon ion radiation therapy in patients with low and intermediate grade chondrosarcoma of the skull base, clinical phase III study. BMC Cancer 2010, 10: 606. 10.1186/1471-2407-10-606
Ares C, Hug EB, Lomax AJ, Bolsi A, Timmermann B, Rut HP, Schuller JC, Pedroni E, Goitein G: Effectiveness and safety of spot scanning proton radiation therapy for chordomas and chondrosarcomas of the skull base: first long-term report. Int J Radiat Oncol Biol Phys 2009, 75: 1111-1118. 10.1016/j.ijrobp.2008.12.055
Jelen U, Bubula ME, Ammazzalorso F, Engenhart-Cabillic R, Weber U, Wittig A: Dosimetric impact of reduced nozzle-to-isocenter distance in intensity-modulated proton therapy of intracranial tumors in combined proton-carbon fixed-nozzle treatment facilities. Radiat Oncol 2013, 8: 218. 10.1186/1748-717X-8-218
Paddick I: A simple scoring ratio to index the conformity of radiosurgical treatment plans. J Neurosurg 2000,92(Suppl 3):219-222.
Ammazzalorso F, Chanrion MA, Graef S, Jelen U: A free software display and analysis tool for photon and particle radiotherapy dose distributions [abstract]. Proceedings to the 52nd Annual Meeting for the Particle Therapy Cooperative Group (PTCOG). Int J Particle Ther 2014, 1: 312–313.
R Development Core Team: R: A Language and Environment for Statistical Computing, Vienna: The R Foundation for Statistical Computing; 2012. ISBN 3-900051-07-0. . [http://www.R-project.org/] R Development Core Team: R: A Language and Environment for Statistical Computing, Vienna: The R Foundation for Statistical Computing; 2012. ISBN 3-900051-07-0. .
Albertini F, Casiraghi M, Lorentini S, Rombi B, Lomax A: Experimental verification of IMPT treatment plans in an anthropomorphic phantom in the presence of delivery uncertainties. Phys Med Biol 2011, 56: 4415-4431. 10.1088/0031-9155/56/14/012
Jia X, Schümann J, Paganetti H, Jiang SB: GPU-based fast Monte Carlo dose calculation for proton therapy. Phys Med Biol 2012, 57: 7783-7797. 10.1088/0031-9155/57/23/7783
Tessonnier T, Mairani A, Cappucci F, Mirandola A, Freixas GV, Molinelli S, Donetti M, Ciocca M: Development and application of tools for Monte Carlo based simulations in a particle beam radiotherapy facility. Appl Radiat Isot 2014, 83: 155-158. 10.1016/j.apradiso.2012.12.019
Jelen U, Radon M, Santiago A, Wittig A, Ammazzalorso F: A Monte Carlo tool for raster-scanning particle therapy dose computation. J Phys Conf Ser 2014, 489: 012013. 10.1088/1742-6596/489/1/012013
Schulte RW, Fargo RA, Meinass HJ, Slater JD, Slater JM: Analysis of head motion prior to and during proton beam therapy. Int J Radiat Oncol Biol Phys 2000, 47: 1105-1110. 10.1016/S0360-3016(00)00551-4
Bolsi A, Lomax AJ, Pedroni E, Goitein G, Hug E: Experiences at the Paul Scherrer Institute with a remote patient positioning procedure for high-throughput proton radiation therapy. Int J Radiat Oncol Biol Phys 2008, 71: 1581-1590. 10.1016/j.ijrobp.2008.02.079
Jelen U: Recommendations for organ depending optimized fixation systems. Deliverable D.JRA 5.1. EU-Project ULICE (GA no. 228436), 2011. . [http://ulice.web.cern.ch/ULICE/] Jelen U: Recommendations for organ depending optimized fixation systems. Deliverable D.JRA 5.1. EU-Project ULICE (GA no. 228436), 2011. .
Jensen AD, Winter M, Kuhn SP, Debus J, Nairz O, Münter MW: Robotic-based carbon ion therapy and patient positioning in 6 degrees of freedom: setup accuracy of two standard immobilization devices used in carbon ion therapy and IMRT. Radiat Oncol 2012, 7: 51. 10.1186/1748-717X-7-51
Nairz O, Winter M, Heeg P, Jäkel O: Accuracy of robotic patient positioners used in ion beam therapy. Radiat Oncol 2013, 8: 124. 10.1186/1748-717X-8-124
Ammazzalorso F, Jelen U, Strassmann G, Engenhart-Cabillic R: Influence of the choice of beam directions on the dosimetric quality and robustness of heavy ion radiotherapy plans [abstract]. In Proceedings of the Heavy Ions in Therapy and Space Symposium: 6–10 July 2009; Cologne, 2009:74.
Lomax AJ: Intensity modulated proton therapy and its sensitivity to treatment uncertainties: 1. the potential effects of calculational uncertainties. Phys Med Biol 2008, 53: 1027-1042. 10.1088/0031-9155/53/4/014
Ammazzalorso F, Jelen U, Krämer M, Strassmann G, Engenhart-Cabillic R: Validation of a homogeneity index for the optimal selection of robust beam configurations in heavy ion radiotherapy planning [abstract]. Radiother Oncol 2009,92(Suppl 1):109. 10.1016/S0167-8140(12)72875-X
Ammazzalorso F, Bednarz T, Jelen U: GPU-accelerated automatic identification of robust beam setups for proton and carbon-ion radiotherapy. J Phys Conf Ser 2014, 489: 012043. 10.1088/1742-6596/489/1/012043
Unkelbach J, Bortfeld T, Martin B. C, Soukup M: Reducing the sensitivity of IMPT treatment plans to setup errors and range uncertainties via probabilistic treatment planning. Med Phys 2009, 39: 149-163. 10.1118/1.3021139
Park PC, Zhu XR, Lee AK, Sahoo N, Melancon AD, Zhang L, Dong L: A beam-specific planning target volume (PTV) design for proton therapy to account for setup and range uncertainties. Int J Radiat Oncol Biol Phys 2012, 82: 329-336. 10.1016/j.ijrobp.2011.05.011
Cabal GA, Jäkel O: Dynamic target definition: a novel approach for PTV definition in ion beam therapy. Radiother Oncol 2013, 107: 227-233. 10.1016/j.radonc.2013.03.010
Weber U, Kraft G: Comparison of carbon ions versus protons. Cancer J 2009, 15: 325-332. 10.1097/PPO.0b013e3181b01935
Fuchs H, Strobele J, Schreiner T, Hirtl A, Georg D: A pencil beam algorithm for helium ion beam therapy. Med Phys 2012, 39: 6726-6737. 10.1118/1.4757578
Grün R, Friedrich T, Elsässer T, Krämer M, Zink K, Karger CP, Durante M, Engenhart-Cabillic R, Scholz M: Impact of enhancements in the local effect model (LEM) on the predicted RBE-weighted target dose distribution in carbon ion therapy. Phys Med Biol 2012, 57: 7261-7274. 10.1088/0031-9155/57/22/7261
Chanrion M-A, Sauerwein W, Jelen U, Wittig A, Engenhart-Cabillic R, Beuve M: The influence of the local effect model parameters on the prediction of the tumor control probability for prostate cancer. Phys Med Biol 2014,59(12):3019-3040. 10.1088/0031-9155/59/12/3019
Lassen-Ramshad Y, Vestergaard A, Muren LP, Høyer M, Petersen JB: Plan robustness in proton beam therapy of a childhood brain tumour. Acta Oncol 2011, 50: 791-796. 10.3109/0284186X.2011.590149
Chanrion MA, Ammazzalorso F, Wittig A, Engenhart-Cabillic R, Jelen U: Dosimetric consequences of pencil beam width variations in scanned beam particle therapy. Phys Med Biol 2013, 58: 3979-3993. 10.1088/0031-9155/58/12/3979
This study was supported by a research grant of the University Medical Center Giessen and Marburg UKGM (Project Number 35/2010 MR) and by the ULICE project co-funded by the European Commission under FP7 (Grant Agreement Number 228436). The authors would like to acknowledge Dr. M. Krämer for enabling the TRiP98 computations, Dr. G. Iancu for providing the proton base data, Dr. A. Wittig for her help during the patient data retrieval phase and free (as in freedom) software for making the study just possible.
The authors declare that they have no competing interests. The authors alone are responsible for the content and writing of the paper.
FA designed the study, performed initial treatment planning, designed the data evaluation process and managed the study data, authored and managed the manuscript. UJ co-designed the study, clinically revised the treatment plans, performed numerical and statistical evaluation and extensively co-authored the manuscript. REC provided clinical input during study design, performed patient selection, reviewed treatment plans and reviewed the final manuscript. WS provided substantial methodological input during study design and data evaluation and reviewed the final manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Ammazzalorso, F., Jelen, U., Engenhart-Cabillic, R. et al. Dosimetric robustness against setup errors in charged particle radiotherapy of skull base tumors. Radiat Oncol 9, 279 (2014). https://doi.org/10.1186/s13014-014-0279-2