Quantifications of Intensity-Modulated Radiation Therapy Plan Complexities in Magnetic Resonance Image Guided Radiotherapy Systems
Article information
Abstract
Background
In this study, the complexities of step-and-shoot intensity-modulated radiation therapy (IMRT) plans in magnetic resonance-guided radiation therapy systems were evaluated.
Materials and Methods
Overall, 194 verification plans from the abdomen, prostate, and breast sites were collected using a 60Co-based ViewRay radiotherapy system (ViewRay Inc., Cleveland, OH, USA). Various plan complexity metrics (PCMs) were calculated for each verification plan, including the modulation complexity score (MCS), plan-averaged beam area (PA), plan-averaged beam irregularity, plan-averaged edge (PE), plan-averaged beam modulation, number of segments, average area among all segments (AA/Seg), and total beam-on time (TBT). The plan deliverability was quantified in terms of gamma passing rates (GPRs) with a 1 mm/2% criterion, and the Pearson correlation coefficients between GPRs and various PCMs were analyzed.
Results and Discussion
For the abdomen, prostate, and breast groups, the average GPRs with the 1 mm/2% criterion were 77.8±6.0%, 79.8±4.9%, and 84.7±7.3%; PCMs were 0.263, 0.271, and 0.386; PAs were 15.001, 18.779, and 35.683; PEs were 1.575, 1.444, and 1.028; AA/Segs were 15.37, 19.89, and 36.64; and TBTs were 18.86, 19.33, and 5.91 minutes, respectively. The various PCMs, i.e., MCS, PA, PE, AA/Seg, and TBT, showed statistically significant Pearson correlation coefficients of 0.416, 0.627, −0.541, 0.635, and −0.397, respectively, with GPRs.
Conclusion
The area-related metrics exhibited strong correlations with GPRs. Moreover, the AA/Seg metric can be used to estimate the IMRT plan accuracy without beam delivery in the 60Co-based ViewRay radiotherapy system.
Introduction
Although intensity-modulated radiation therapy (IMRT) delivers prescription doses to target volumes while sparing the normal tissues [1–3], it is often associated with high degrees of modulations including severe multi-leaf collimator (MLC) movements, numerous monitor units (MUs), and the frequent use of small and irregular beam apertures [4–6]. This complicated treatment techniques may cause both delivery errors and dose calculations [7–9]. Therefore, the accuracy of the plan delivery should be verified in advance of the treatment, and a gamma evaluation has been widely adopted as a pre-treatment quality assurance (QA) [10, 11]. While gamma analysis is practical by comparing the calculated and measured dose distributions, there exists several limitations such as clinical irrelevances and consumptions of human and clinical resources [12–14]. To substitute or support the gamma analysis platform, several studies have attempted to characterize degrees of plan modulation for IMRT and volumetric modulated arc therapy (VMAT) plans [4, 6, 7, 15–17].
IMRT and VMAT are known to be capable of allowing dose escalation and reduced toxicity to the organs at risk when compared with three-dimensional (3D) conformal radiation therapy [18, 19]. Although step-and-shoot (SS) IMRT is less complex [19], many studies have attempted to quantify the degree of modulation and complexity for the SS IMRT plans. McNiven et al. [4] suggested the modulation complexity score (MCS) to quantify the complexity and deliverability of SS IMRT plans, and reported improved sensitivity using the MCS compared with the number of MUs. Du et al. [15] devised the plan-averaged area (PA), the plan-averaged irregularity (PI), the plan-averaged modulation (PM), and plan-normalized monitor unit metrics to evaluate the complexity of the SS IMRT plans quantitatively. They demonstrated a negative correlation between the minimum segment area and PM, indicating that the small beam apertures in spine stereotactic plans could induce large discrepancies between the planned and measured doses. The complexity of the created static MLC openings (defined as regions with 5% dose difference between the measurement and calculation), devised by Gotstedt et al. [20], was correlated to the aperture-based complexity metrics and dependent on the penumbra region relative to the area of the opening.
Recently, the introduction of magnetic resonance-guided radiation therapy (MRgRT) has enabled real-time respiratory gating without an external surrogate [21–24]. The ViewRay system (ViewRay Inc., Cleveland, OH, USA) is the world’s first commercial MRgRT system that integrates an MR imaging system, a radiation therapy (RT) delivery system, and an integrated treatment planning system [22–25]. The RT system consists of three 60Co heads (120° apart), each of which is equipped with 60 double-focused MLCs without collimator jaws, with the beams delivered using the SS technique [22–26]. While this unique beam delivery method emphasizes the quantification of the plan complexity, this has rarely been studied for the 60Co-based ViewRay system to date. Owing to the distinctive process of beam delivery and unique collimator geometry, i.e., 60 double-focused MLCs without a primary collimator [22–26], those dedicated to the ViewRay system should be individually evaluated. Recently, Liu et al. [26] attempted to quantify the treatment plan complexity for this system, demonstrating a linear relationship between the total treatment time and PM.
Due to current pre-treatment QA mostly relies on the measurement-based gamma analysis, estimation of gamma passing rate (GPR) with complexity metrics could provide foresights in terms of accurate delivery. When plans with low GPR estimation, they can be modified without actual QA beam delivery. The increased use of adaptive RT, especially in MRgRT systems, also necessities the prediction of the GPR in advance of treatment [27, 28]. The online adaptive plan has different plan parameters compared with the initial plan such as gantry angles, segment shapes, and beam-on times, whose deliverability should be re-verified by comparing the recalculated dose distribution and delivery for the re-optimized plan. However, the fraction-specific QA for each adaptive plan is impractical due to potential anatomical changes when the patient is transferred from the table [29–32].
In this study, we evaluated various plan complexity metrics (PCMs) subject to the SS IMRT in the MRgRT system with a relatively large number of dataset in multiple disease groups. The impacts of PCMs on delivery accuracies were demonstrated with correlation analysis across PCMs, GPRs, and parameter differences between planned and actual parameters.
Materials and Methods
1. Data Selection and Treatment Planning
Following approval by the Institutional Review Board of Seoul National University Hospital (No. 1906-061-1039), 194 verification plans from 113 patients were retrospectively collected in this study. Treatment sites included abdomen, prostate, and breast cancer patients. Abdomen and prostate cancer patients were treated with stereotactic ablative body radiotherapy (SABR) schemes, while breast cancer patients were consisted of accelerated partial breast irradiation (APBI) plans with hypofractionated IMRT. Patients with pancreatic and liver cancers are categorized into abdomen groups.
All patients were scanned with the Brilliance CT Big Bore (Phillips, Amsterdam, The Netherlands) and ViewRay 0.35 T magnetic resonance imaging system on the same day. The computed tomography (CT) images were rigidly deformed to MR images for dose calculations with an electron density map. The treatment plans were established using the MRIdian treatment planning system (ViewRay Inc.) with the SS IMRT techniques. The value of IMRT efficiency was set to 0.5, and prescriptions and optimization objectives were described as follows. For abdomen groups, the initial prescription dose was 50 Gy in five fractions. If any of doses received by 1 cc volume (D1cc) of a duodenum, a bowel, and a stomach exceeded 33 Gy, the prescription dose was reduced until their tolerance levels were met and optimization was repeated [33]. Also, dose-volumetric constraints of D15cc less than 32 Gy for a heart and D700cc less than 21 Gy for a normal liver were applied [34, 35]. For prostate cancer patients, the prescription dose was 36.25 Gy in five fractions, and volume receiving more than 100% of the prescription dose (V100%), V90%, V80%, and V50% less than 5%, 10%, 20%, and 50%, respectively, were used for bowel constraints [33]. For breast cancer patients, they were prescribed by 38.5 Gy in 10 fractions, and constraints of V20Gy less than 20% for an ipsilateral lung, and mean dose less than 10 Gy for a heart were applied [35].
The verification plans were then generated by replacing patients’ images with the MR-compatible ArcCHECK (Sun Nuclear Corp., Melbourne, FL, USA) phantom images. The reference dose distributions were calculated by using a Monte Carlo simulation-based dose calculation algorithm with a 3 mm grid resolution supported by the MRIdian treatment planning system. Calculations with magnetic field option were applied.
2. Measurements and Gamma Analysis
Each dose distribution was measured with the MR-compatible ArcCHECK consisting of the SunPoint diode detector (Sun Nuclear Corp.). The spacing and active volume of the detector are 10 mm and 0.019 mm3, respectively. Prior to the measurements, array calibration and absolute dose calibration of the ArcCHECK were performed according to the manufacturer’s calibration procedure. Two-dimensional gamma analysis was performed using the SNC patient software (version 6.6; Sun Nuclear Corp.) with a 10% threshold. The dose distribution was automatically interpolated with 1 mm dose grid in the software. The global gamma analysis was performed with gamma criterion of 1 mm/2%.
3. Plan Complexity Metrics
As the radiation isotope (60Co) is used as a source of the ViewRay system, the beam delivery is governed by the beam-on time. In this study, to account for this property in the PCMs, the MU weighting factor in the initial definition was replaced by the time weighting factor to emphasize the segment and the field with a longer beam-on time. Each plan consists of multiple gantry angles called as “field,” which possess their own beam aperture shapes (segment). In Table 1, the modified definitions of the PCMs are summarized with respect to the segment phases, field phases, and plan phases.
MCS suggested by McNiven et al. [4], and PA, PI, plan-averaged edge (PE), and PM suggested by Du et al. [15], and Younge et al. [36] were utilized. Same as characteristics of the initial definitions, plans with high complexity possessed low MCS and PA values, and high PI, PE, and PM. In additions, an average aperture area per each segment (AA/Seg), and total beam-on time (TBT) along with total beam-on time per daily dose (Time/DD) were analyzed. To investigate the impact of segments with small area, and short beam-on time, the percent proportion of segment with beam-on times less than 3, 5, and 10 seconds (%Time/Seg<3, %Time/Seg<5, and %Time/Seg<10), and average areas among all segments less than 2.5, 5.0, 7.5, and 10.0 cm2 (%SA/Seg<2.5, %SA/Seg<5.0, %SA/Seg<7.5, and %SA/Seg<10.0) were used for analysis [37].
4. Statistical Analysis
To find the significant differences on PCMs, F-test was firstly performed to determine whether two values possess an equality of variances, and then t-test was followed to find the statistical significances. The Pearson correlation coefficient between GPRs and various PCMs were calculated. Only coefficients with statistical significances (p<0.001) are presented.
Results and Discussion
The GPRs with a 1 mm/2% criterion were 77.8±6.0%, 79.8± 4.9%, and 84.7±7.3%, for abdomen, prostate, and breast groups, respectively. Breast IMRT plans exhibited significantly higher GPRs than those of all the other SABR plans, on average. It seems appropriate for SNC software interpolates the calculated dose distribution with 1 mm grid as gamma analysis between the dose distribution calculated with 2.5 mm grid (1.25 mm interpolation) and 1.25 mm grid showed 97.9% GPRs by using 1% and 0.1 mm criterion, and 0% threshold [38].
The representative aperture shapes along with the aperture area (AA), the aperture irregularity (AI), the aperture edge (AE), and the time weighting factor for each segment in the abdomen SABR plans are presented in Fig. 1. The abdomen SABR plan exhibited higher complexity than the APBI IMRT plans with smaller and much irregular beam apertures, resulting in a lower beam area metric (BA) and higher beam irregularity (BI), beam edge (BE), and beam modulation (BM) in the field phase.
Statistics for PCMs according to the treatment sites, along with statistical significances in PCMs between breast IMRT and SABR groups, and those between abdomen SABR and prostate SABR groups were also provided. Breast IMRT plans showed the statistically higher MCS, PA, and AA/Seg, and lower PI, PE, PM, TBT, and Time/DD compared to SABR plans (Table 2).
The Pearson correlation coefficients between various PCMs and GPRs with a 1 mm/2% criterion are presented in Table 3, and scatter plots were provided in Fig. 2. The majority of PCMs showed significant correlations with GPRs except for PI, PM, and #Seg. The majority of area-related metrics such as PA, PE, AA/Seg, %SA/Seg<2.5, %SA/Seg<5.0, %SA/Seg<7.5, and %SA/Seg<10.0 showed significant correlations with GPRs (p<0.001). Furthermore, segments containing short beam-on time also showed significant correlations with GPRs (p<0.001). A strongest correlation with the GPR was occurred with AA/Seg metric showing Pearson correlation coefficient of 0.635.
Study results revealed area-related metrics tended to possess strong correlations with GPRs. In small fields, the accuracies of calculations and measurements became could become poor due to the lack of lateral charged particle equilibrium and overlapping penumbras [39–41]. As several researches reported the presence of the MR could impact the dosimetric performances [42–45], the inaccuracies in small fields might also be occurred same as the absence of the MR fields [46]. When using highly modulated fields with large number of small segments, inaccurate output factors may cause uncertainties in dose calculation [47]. Our study results showed concordant tendencies with previous studies showing plans with small PA or AA/Seg possessed low GPRs. Specifically, those impacts are more powerful than those with time parameters.
Despite the lack of contradictory correlation tendencies compared to another literature, the correlation strength between the PCMs and GPRs was relatively low [5]. This might be because the number of fields was greater than that reported in a previous study [5], which in turn resulted in the complexity for each field being potentially smeared out by utilizing higher number of fields. Future studies will focus on a more precise and direct evaluation of the plans only consisting of complex beam apertures. Even so, the study outcomes present the large proportion of segments with small apertures and short beam-on times might reduce the deliverability of the plans.
Most previous studies on patient-specific QA for online adaptive RT plans did not perform a conventional measurement-based QA. In the study by Acharya et al. [29], their in-house QA tools provided comparisons of important plan parameters between initial and re-optimized plans along with independent Monte Carlo-based dose calculations and 3D gamma analysis [48]. Lamb et al. [49] discussed the challenges of measurement-based QA for online adaptive plans due to potential alternations in the anatomy when the patient is transferred from the table. They retrospectively delivered 20 plans of adaptive fractions to a QA phantom, all of which passed the tolerance level (i.e., >95% GPR with a 3 mm/3% criterion). An automated process to check the consistency of contour volumes, segment beam-on time, and dose-volumetric parameters is also a meaningful approach to mitigate unexpected errors [49]. Several alternative approaches in fraction-specific measurement-based QA are routinely used clinically or are under investigation such as secondary dose calculations, automated parameter comparisons between the initial and re-optimized plans, and MR-compatible electronic portal imaging detector-based QA workflows [29, 30, 48, 50].
Although 60Co ViewRay system has been replaced to magnetic resonance linear accelerator (MR-Linac) system recently, there are still several institutions using this system. Furthermore, the study outcomes highlights that the large amount of small segment aperture could aggravate the plan delivery accuracy, which can also be applied in all SS IMRT with MR-Linac system.
Conclusion
The area-related metrics, such as PA, AA/Seg, and average areas less than specific area, exhibited the strong correlations with GPRs. As all metrics can be obtained at the planning stage, plans with small AA/Seg might be modified without actual plan delivery in the 60Co-based MRgRT system.
Notes
Conflict of Interest
No potential conflict of interest relevant to this article was reported.
Author Contribution
Conceptualization: Chun M. Data curation: Kwon O. Formal analysis: Kwon O. Funding acquisition: Chun M. Methodology: Chun M. Visualization: Kwon O. Writing - original draft: Chun M. Writing - review & editing: Park JM, Kim J. Investigation: Kim J. Resources: Kwon O. Software: Chun M. Supervision: Park JM, Kim J. Validation: Chun M, Kwon O.
Acknowledgements
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (No. NRF-2019R1F1A1041944).