TrueBeam Low Dose Rate Investigation for Pulsed Reduced Dose Rate IMRT ()
1. Introduction
Pulsed reduced dose-rate radiotherapy (PRDR) has become a powerful treatment for recurrent diseases, which comprise approximately 15% - 20% of treatments [1] [2] [3] . Introduced as an extension of continuous low dose rate therapies, PRDR combines the advantages of accelerator-based delivery and low dose hyper-radio-sensitivity by separating conventional dose fractions into 0.2 Gy pulses separated by three-minute intervals to achieve a time averaged dose rate of 0.0667 Gy/min [4] [5] . This hyper-radio-sensitivity is caused by accumulation of tumor cells in the G2 and M phases of the cell cycle during treatment.
Due to the complex nature of these treatments, PRDR using 3D-conformal techniques can be very difficult. Often patients have been previously irradiated; as such, nearby critical structures may have already received at or close to their toxicity limit. Researchers are actively investigating ways to apply PRDR to more advanced techniques, such as static field Intensity Modulated Radiation Therapy (IMRT) [6] , Volumetric Modulated Arc Therapy (VMAT) [7] [8] , and tomotherapy [9] [10] . The primary challenge with IMRT-based PRDR, however, is that IMRT fields are difficult to separate into 0.2 Gy pulses. Currently available treatment planning systems do not include optimization constraints that specifically limit dose rate. Instead, treatments must be split into multiple (often ten or more) beams or arcs each delivering less than 0.2 Gy, or the same beams/arcs must be delivered repeated times with each sub-fraction less than 0.2 Gy [7] .
When delivering PRDR treatments on the Varian® TrueBeam™ system (Varian Medical Systems, Palo Alto, CA), an alternative technique can be used to reduce the time averaged dose rate. Instead of splitting the total dose into 0.2 Gy pulses, a lower repetition rate can be used to achieve a time averaged dose rate of 0.0667 Gy/min. The TrueBeam system can run as low as 5 MU/min. This greatly simplifies treatment planning of PRDR cases, as conventional (normal fractionation) beam arrangements and optimization techniques can be used. After optimization, the dose from each beam to the target can be computed and the dose rate adjusted during treatment delivery to achieve the desired target dose rate.
For VMAT and sliding window IMRT techniques, reducing the repetition rate will significantly slow down the gantry and MLC leaf motion, respectively. Even with static field step and shoot-type plans, the low repetition rate may affect the operation of the servo mechanisms that maintain beam stability parameters such as flatness, symmetry, and dose per MU. Typical IMRT commissioning practices also may not specifically include testing delivery at these low repetition rates [11] .
Furthermore, performing patient specific plan delivery Quality Assurance (QA) at such a low repetition rate can also be problematic. Jursinic et al. [12] demonstrated that some diodes display a reduced sensitivity relative to ion chamber measurements as the dose rate is reduced. Ahmed et al. [13] conducted a similar experiment on a TrueBeam system and Sun Nuclear® ArcCHECK® measurement phantom (Sun Nuclear Corporation, Melbourne FL), finding that the ArcCHECK under-responded by −2.2% relative to ion chamber measurements at 5 MU/min. This adds additional challenges when commissioning a PRDR program; prior to performing end-to-end testing on an IMRT QA phantom, users must first re-validate their QA system at these repetition rates and factor their devices’ limitations into their results.
After commissioning is completed, there is an additional logistical concern as well. Conducting PRDR IMRT delivery QA measurements at the true low dose rate requires at least 30 minutes per run. Plans often require more than one measurement as well (either for troubleshooting or to increase the detector spatial frequency/area), making it difficult to conduct these measurements after hours in a busy clinic.
The purpose of this study was to demonstrate a procedure by which other clinicians could follow to develop PRDR at their clinics. The procedure starts with a series of tests to validate the delivery system at low repetition rates, then to characterize the performance of the IMRT QA system, and finally a series of end-to-end tests of patient plans that establish equivalence between IMRT delivery at PRDR and conventional repetition rates. To the author’s knowledge, this is also the first reported characterization of the ScandiDos Delta4® phantom (ScandiDos AB, Uppsala Sweden) at low repetition rates.
2. Methods
All measurements were performed on the TrueBeam 2.5 treatment system using the conventional flattened 6 MV beam. IMRT QA phantom measurements were acquired on the Delta4 phantom and October 2016 software. Where stated that measurements were performed at each dose rate, acquisitions were repeated at 5, 10, 15, 20, 40, 60, 80, 100, and 600 MU/min. Unless otherwise noted, static beams were delivered using a 20 × 20 cm2 open field, IEC 0˚ gantry/collimator positions, with the device set up to 100 cm Source to Surface Distance (SSD) using the optical distance indicator and aligned to the treatment system central axis using the light field, verified by the lasers.
2.1. Beam Characteristic Equivalence
Monitor unit dose, linearity, and beam flatness/symmetry reproducibility with dose rate were evaluated using the Sun Nuclear IC Profiler™ device and v3.4.2 software. The IC Profiler was chosen as it has been previously demonstrated to be equivalent to water tank measurements detecting changes in these metrics [14] . The treatment system was then configured to deliver 50 MU at each dose rate. For each measurement, the central axis integral counts, IEC X/Y area symmetry, and diagonal flatness, statistics were computed and recorded.
The area symmetry is defined by Equation (1) [15] , where and are the trapezoidal integrated area under the left and right central 80% field, respectively:
(1)
The diagonal flatness is defined as the average dose at 80% of the field width along each diagonal sub-axis normalized by the central axis [15] , defined by Equation (2) where is the response at distance d from the central axis along diagonal sub-axis i:
(2)
The symmetry and flatness statistics from 5 to 100 MU/min tested were tested for equivalence to their respective 600 MU/min values using a paired Two One-Sided Test (TOST) by Schuirmann [16] and Westlake [17] . The robust t-test of Yuen and Dixon [18] [19] was used as it makes no assumption of normality. An equivalence margin of 0.5% was used. In addition, the results were tested for linear dependence with dose rate using an Iterative re-Weighted Least Squares (IWLS) robust regression was performed using an MM-type estimator and Tukey bi-square [20] re-descending score function described in Yohai [21] and Koller and Stahel [22] .
Next, an IBA FC65-P Farmer chamber (Ion Beam Applications SA, Louvain-la-Neuve Belgium) was placed at 1.5 cm depth in solid water. The chamber was then connected to a Standard Imaging® SuperMax electrometer (Standard Imaging Inc., Middleton, WI). The collected charge was recorded for 5, 10, 15, 25, and 50 MU delivery using a 20 × 20 cm2 open field, 0˚ gantry, and 5 MU/min dose rate and corrected for temperature and pressure. The MU linearity residual errors were computed as the observed relative corrected charge minus the predicted value given a linear fit through 50 MU. The standard error was computed as the root mean square of all residuals.
2.2. Delta 4 Static Beam Dose Rate Equivalence
Delta 4 response reproducibility with dose rate was evaluated by comparing the Delta 4 response relative to an ionization chamber. Using the same ionization chamber setup as before, the collected charge was recorded for 50 MU at each dose rate. The measurements were then repeated with the Delta 4 aligned to isocenter. For each measurement, mean dose across the central 6 × 6 cm2 detectors was recorded and divided by the ionization chamber response. The ratios were then normalized to 600 MU/min and tested for equivalence and linear dependence as before using a 0.2% equivalence margin. Verification of MU linearity was also repeated with the Delta 4.
2.3. Delta 4 Step and Shoot IMRT Equivalence
To validate the equivalence of Delta 4 IMRT QA results with dose rate, three PRDR static step and shoot IMRT QA plans were delivered at low dose rates and 100 MU/min. Each plan was optimized in the Pinnacle3 v9.8 Treatment Planning System (Philips Healthcare, Madison, WI). A summary of the beams in each plan are provided in Table 1. All patients were brain treatments, as this is the most common indication for PRDR at the University of Wisconsin-Madison. Plans were chosen with varying target sizes, creating a range of field sizes from 27.5 to 196 cm3. The average number of MU per segment ranged from 3.1 to 12.3, with a mean of 4.7 MU. Dose rate ranged from 5 to 20 MU/min.
Next, each plan was copied to the Delta 4 and re-calculated using a 2 × 2 × 2 mm3 resolution. Each plan was then delivered to the Delta 4 at the low dose rate and at 100 MU/min and analyzed in the Delta 4 software. The right parietal lobe plan was measured both 5 MU/min and 10 MU/min to provide additional statistical power. The median absolute dose in the target region defined by diodes above 50% of the maximum dose and mean Gamma index [23] using 2% (of maximum dose), 2 mm distance to agreement, and 20% threshold criteria were then recorded. The dose difference and mean Gamma at low and high dose rates were tested for equivalence using a paired TOST, with equivalence margins of 0.2% and 0.05, respectively.
Table 1. Summary of step and shoot IMRT plan parameters.
Finally, to evaluate if dose rate affects the variance in IMRT QA results the right parietal lobe field 5 MU/min, 10 MU/min, and 100 MU/min beam measurements were each repeated three times. The standard deviation of each measurement set was then computed and tested for equivalence using the same paired TOST.
3. Results and Discussion
The following sections detail the results of the tests described above. Statistical analysis was conducted in R v.3.3.2 (R Foundation, https://www.r-project.org). A significance level of 0.05 (95% confidence) was established a priori for hypothesis testing.
3.1. Beam Characteristic Equivalence
Figure 1 illustrates the central axis dose, diagonal flatness, and area symmetry (along X and Y axes) with dose rate. At 600 MU/min, these values are consistent with peers’ results from other TrueBeam clinical systems [24] . The central axis dose was found to be equivalent with doserate, ε = 0.5%, 95% CI [−0.09% 0.11%], p < 0.001, with no linear dependence observed, p = 0.111. The diagonal flatness was not equivalent, ε = 0.5%, 95% CI [−0.91% −0.35%], p = 0.795, and instead did display a linear dependence, p = 0.023, adj. R2 = 0.73. This is likely due to the increased noise in the measured profiles as the dose rate decreased; the Y axis profile at 5 MU/min is plotted (in magenta) over the other dose rate profiles in Figure 2.
Area symmetry along the X axis was equivalent, ε = 0.5%, 95% CI [−0.37% 0.08%], p = 0.009, although with an inverse linear relationship with dose rate, p = 0.020, adj. R2 = 0.58. Area symmetry along the Y axis was equivalent, ε = 0.5%, 95% CI [0.22% 0.45%], p = 0.014, with no statistical linear dependence, p = 0.972.
Figure 1. Beam characteristics (dose, flatness, and symmetry) with dose rate.
Figure 2. Beam Y axis profile at different dose rates. The magenta (noisy) profile is 5 MU/min.
Changes in symmetry with dose rate suggest that the dose servos have a slight effect on beam steering. That said, AAPM Task Group 142 [25] recommends period QA flatness and symmetry changes from baseline to not exceed 1% from baseline; at 1% equivalence margin all beam characteristics were statistically significant including flatness, p < 0.019. Finally, for MU linearity, the residual standard and maximum errors were 0.25% and 0.55% at 5 MU/min between 5 - 50 MU, respectively. The residual errors are plotted in Figure 3.
Figure 3. MU linearity residual errors for ion chamber and Delta 4.
Figure 4. Delta 4 response relative to ion chamber with dose rate for a 20 × 20 field, normalized to 600 MU/min.
3.2. Delta 4 Static Beam Dose Equivalence
Figure 4 illustrates the Delta4 mean response (within the central 6 × 6 cm2 de-
tectors) relative to ion chamber response under the same conditions as a function of dose rate, normalized to 1.000 at 600 MU/min. The low dose rates (5 - 100 MU/min) were found to be equivalent to 600 MU/min, ε = 0.2%, 95% CI [−0.20% −0.03%], p = 0.049, as well as no linear dependence with dose rate, p = 0.763.
These findings suggest that the Delta4 diode design does not suffer from the low dose rate under-response reported with the ArcCHECK [13] and other diodes [12] . According to the vendor, the proprietary design of the diodes was optimized to maximize response constancy with delivery parameter variability, including PRF rate (personal communication, January 2017). For MU linearity, the residual standard and maximum errors were 0.40% and 0.76%, respectively. Residual errors are plotted in Figure 3.
3.3. Delta4 Step and Shoot IMRT Equivalence
Paired measurements (at 5 - 20 MU/min and 100 MU/min) of each beam in Table 1 were collected for equivalence testing (n = 41), and are plotted in Figure 5. The median absolute dose was found to be equivalent between the two dose rate groups, ε = 0.2%, 95% CI [−0.015% 0.13%], p = 0.001, while the mean 2%/2 mm Gamma index was also found to be equivalent, ε = 0.05, 95% CI [−0.0034 0.029], p < 0.001. The maximum difference in median absolute dose between groups was 0.8%. IWLS robust regression yielded a significant correlation, p < 0.001, adj. R2 = 0.99.
Figure 5. Scatter plot of per-field Delta 4 median absolute dose difference (relative to TPS) at 5 - 20 MU/min versus 100 MU/min. The dashed line indicates identical response.
Furthermore, the variance in measurement was also found to be equivalent between dose rate groups (n = 5) for the median absolute dose, ε = 0.2%, 95% CI [−0.0056% 0.080%], p < 0.001, and Gamma index metrics, ε = 0.05, 95% CI [−0.0084 0.025], p = 0.003.
Finally, using the planned low dose rate, the composite median dose difference (between measured and calculated by the treatment planning system) in the target region and 3% (of maximum dose), 3 mm, 20% threshold Gamma index met the clinic’s IMRT QA acceptance criteria for all plans. The mean, minimum, and maximum median dose differences were 1.4%, 0.9%, and 2.2%, while the Gamma index pass rates were 99.3%, 98.2%, and 100%, respectively. It is notable that all PRDR plans yielded a median difference greater than one; unfortunately, not enough plans were measured to determine if the distribution of PRDR plan results statistically differ from other types of treatments. As this technique is employed at this center and more plans are created, this is an opportunity for future study.
4. Conclusion
This report demonstrates the capability of the TrueBeam system to deliver step and shoot PRDR IMRT plans and establishes equivalence between IMRT QA at 100 MU/min and lower dose rates. Static beam measurements found the TrueBeam system to achieve MU constancy, linearity, and symmetry within 0.5% margin. The Delta 4 phantom also responded equivalently within a 0.2% margin between ion chamber measurements for static beams and between 100 MU/min and lower dose rates for step and shoot IMRT.