A Comparison between Three-Dimensional Conformal Radiotherapy, Intensity-Modulated Radiotherapy, and Volumetric-Modulated Arc Therapy Techniques for Stereotactic Body Radiotherapy of Lung Tumors ()
1. Introduction
Lung cancer is the most common malignancy worldwide and account for the most cancer deaths among both men and women [1] . A widely used treatment modality for inoperable lung cancer, especially early-stage non- small cell lung cancer (NSCLC) is stereotactic body radiation therapy (SBRT) [2] . Small sized fields with a sharp penumbra and high fractional dose require the use of advanced immobilization, image-guided localization and more accurate treatment planning procedures. It is common for the majority of stereotactic body radiation therapy (SBRT) of lung tumors to be treated with Intensity Modulation Radiation Therapy (IMRT) [3] . The increased use of SBRT has created the need to evaluate various treatment modalities available including three- dimensional conformal radiotherapy (3D CRT), IMRT and the relatively new volumetric-modulated arc therapy (VMAT) [4] . SmartArc (SA) is one version of VMAT that modulates the intensity of the beam as it rotates around the patient. SA delivery method in combination with inverse planning has the ability to produce highly conformal dose distributions while using greatly reduced treatment times [5] . A few reports have demonstrated that SA can deliver highly conformal plans for SBRT of lung [6] - [8] . RapidArc (RA) is another novel implementation of VMAT that has been reported to have superior target dose conformity using less treatment time [3] [8] . The aim of this study is to dosimetrically compare 3D CRT, IMRT, and SA techniques that employ coplanar and non-coplanar beam arrangements for SBRT of lung tumors.
2. Methods and Materials
2.1. Planning
Nine (n = 9) consecutive SBRT lung patients previously treated at our institution were selected for this study including a patient with bilateral targets. Patients were simulated supine in the SBRT Pro-Lok immobilization device (CIVCO Medical Solutions, Orange City, IA) with abdominal compression. Patients underwent a 4D-CT scan using the Real-time Position Management (RPM) Respiratory System (Varian Medical Systems, Palo Alto, CA) on a GE LightSpeed 16-slice CT scanner (GE Healthcare, Waukesha, WI). The 4D-CT images were used to generate the internal target volume (ITV). Plans were created for treatment on a 23EX linear accelerator (LINAC) (Varian Medical Systems, Palo Alto, CA) using a 6MV photon beam.
For treatment planning, the ITV was expanded by 5mm to generate the planning target volume (PTV). The PTV volumes ranged from 6.46 cc to 40.45 cc with locations in the lung as shown in Table 1. The prescription dose (Rx) of 50 Gy is delivered in 5 fractions. The organs at risk (OAR) include the lungs, spinal cord, aorta, heart, trachea, liver, and stomach. The dose constraints are based on the TG101 report for a 5 fraction SBRT [4] . Highlights of the dose constraint from the report are shown in Table 2. Tissue heterogeneity correction and 2 mm dose grid spacing were used in planning. Care was taken to ensure that the higher isodose lines or hotspots occur within the PTV and not in adjacent normal tissue.
Pinnacle3 treatment planning system’s (Philips Medical System, Andover, MA) collapsed cone convolutionsuperposition (CCCS) algorithmwas used in creating the plans. Six (n = 6) treatment plans were created for each case: 1 coplanar (C) and 1 non-coplanar (NC) 3D CRT, IMRT, and SA and normalized such that the Rx covers 98% of the PTV. While the 3D CRT and IMRT plans consisted of 7 non-opposing fields, the SA-C plans had 2 coplanar arcs of 360˚ in the clockwise and counter-clockwise direction. The SA-NC plans consisted of 2 shortened arcs with 10˚ - 15˚ couch angles and gantry rotation spanning 160˚ - 220˚, thereby avoiding the contralateral lung. The IMRT and SA plans used identical dose-volume constraints for optimization. Direct machine parameter optimization (DMPO) feature was used during optimization thereby eliminating the need to convert the
Table 1. The size and location of the ten (n = 10) PTVs used in this study.
Table 2. Highlights of the TG-101 OAR dose constraints for a 5-fraction SBRT of lung tumor.
beam fluence generated to deliverable dose. Following optimization, a final dose calculation was performed using CCCS algorithm at 2 mm dose grid spacing.
2.2. Planning Criteria
Achieving a high degree of target conformity in SBRT lung tumors while sparing the OARs can be quite complex. Plan quality metric for the target includes the ICRU Conformity index (CI) [9] , R50%, and homogeneity index (HI). CI is defined as:
(1)
The closer the value of CI to 1.0, the better would be the dose conformity. R50% index is an indication of the low dose spillage and is defined as:
(2)
Lower R50% ratio indicates greater dose falloff and better dose conformity around the PTV. One way to calculate the homogeneity of the target dose distribution uses the difference in doses delivered to 2% and 98% of the volume of PTV represented as D2%T and D98%T, respectively and the mean dose to the PTV (DMean).
HI = (D2%T − D98%T)/DMean (3)
Values of HI closer to 0 indicate greater dose homogeneity within the volume of PTV while large values indicate more heterogeneous dose distribution. Target coverage depends on the fractional volume of the PTV receiving the Rx dose and is defined as:
(4)
With respect to OARs, the maximum dose to the spinal cord, aorta, heart, trachea, liver, and stomach were tracked. The percent volume of total lung that received 5 Gy (V5), 20 Gy (V20) as well as the normal tissue irradiated by 5 Gy and 25 Gy isodose line were recorded. Beam angles and arc lengths were chosen for the six treatment plans in order to meet 98% coverage of the PTV and the least value of total lung V5 and V20 and R50%. The data were extracted from the plan dose volume histogram (DVH) and were averaged for each patient and for each planning technique. The isodose distribution and the corresponding DVHs were shown for one representative patient in Figure 1 and Figure 2 respectively. The indices were compared across the six treatment techniques for the same beam configuration using a 2-tailed Student’s paired T-test.
Figure 1. Isodose line comparison for PTV #4. Top row, left to right: coplanar 3D CRT, IMRT, and SA. Bottom row, left to right: non-coplanar 3D CRT, IMRT, and SA.
Figure 2. DVH of the GTV, PTV, total lung, left lung, cord, and Aorta for PTV #4.
3. Results
Prior to analysis, the patient plans are sorted by the volume of PTV, as shown in Table 1 and in all the figures. Equation (1), Equation (2) and Equation (3) were used in tabulation of values of CI, R50% and HI from the 6 treatment plans, respectively. Dose conformity of the PTV coverage could be judged by the value of CI. Ideally, CI should be less than 1.2 and values up to 1.5 are considered as a minor deviation for a clinically deliverable lung SBRT plan, per RTOG 0813. As shown in Figure 3, the measured CI values < 1.5, with the exception of the 3D CRT coplanar plans of PTV# 1, 6 in which case CI values < 1.75.
The low dose spillage depends on the dose falloff gradient beyond the PTV extending into normal tissue and is measured using the R50% parameter. The measured R50% values ranged between 3.4 and 6.8, as shown in Figure 4. This can be compared against the acceptable R50% values tabulated under 2 categories―no deviation and with minor deviation in Table 1 in RTOG 0813. The generated plans had larger than a minor deviation in the R50% values of 3D CRT coplanar plan in PTVs 4, 5 and for both configuration of 3D CRT plans in PTV#9.
The target dose homogeneity was measured by the ratio of the difference between doses that irradiate 2% and 98% of the PTV and the mean dose to PTV. Shown in Figure 5 are the measured HI values that range between 0.25 and 0.6.
Shown in Table 3 is a summary of the PTV and OAR dosimetric comparison of the six treatment techniques averaged across patients. The measured data for PTV is plotted in Figure 6. IMRT and SA produced the most conformal plans with mean CI values close to unity (1.04). The mean CI values for coplanar and non-coplanar 3D CRT plans were on average significantly higher by 26% and 19% with p-value < 0.01, 0.01 respectively. The CI values of coplanar and non-coplanar SA plans were not significantly different from those of the corresponding IMRT plans.
As shown in Figure 6, the mean R50% values were lowest for the non-coplanar SA plans with a value of 4.19. The mean R50% values for the 3D CRT non-coplanar and IMRT non-coplanar plans were 13% and 16% higher,
Figure 3. Conformity index for the 6 plans generated for each of the 10 PTVs. Notice that with the exception the 3D CRT coplanar plans of PTV# 1 and 6, CI < 1.5.
Figure 4. R50% indices for the 6 plans generated for each of the 10 PTVs.
Figure 5. Homogeneity index for the 6 plans generated for each of the 10 PTVs.
Figure 6. Dosimetric comparison of treatment techniques based on the CI, HI, R50% and Coverage values averaged among the 10 patient plans.
Table 3. Summary of dosimetric comparison between the six treatment techniques expressed as mean values ± standard deviation. *and +indicates statistical significance based on Student’s 2-tailed T-test (p-value < 0.05) with respect to the corresponding coplanar (C) or non-coplanar (NC) configuration of 3D CRT, IMRT plans respectively.
respectively. The R50% values of non-coplanar SA plans were significantly different from that of both IMRT non-coplanar and 3D CRT non-coplanar plans with p-value < 0.05 and p-value < 0.02, respectively. The R50% values of coplanar SA plans were significantly different from that coplanar 3D CRT plans with p-value < 0.001, but not significantly different from the coplanar IMRT plans with p-value > 0.10.
The mean HI values range between 0.35 and 0.44 for the 6 treatment plans. Dose heterogeneity in the PTV is significantly higher in the IMRT coplanar plans than the corresponding 3D CRT and SA coplanar plans with p- value < 0.005 and p-value < 0.05, respectively.
In the measured OAR data, the quantities analyzed include V5, V20 for the total lung, and volumes of normal tissue irradiated by 5 Gy, 25 Gy as well as the maximum dose to spinal cord, aorta, heart, trachea, liver, and stomach. All six treatment techniques were able to meet the OAR dose constraints and the only significant reduction in OAR dose came from coplanar VMAT plans for a 35.5% reduction in the maximum cord dose than coplanar 3D CRT plans (p-value < 0.03). The maximum dose to other OARs were not significantly different (p- value > 0.05). The deviations in V5 and V20 of the total lung are negligible among the six treatment techniques, as shown in Figure 7, and the coplanar SA plans had the highest mean V5 value of 0.16 ± 0.06.
Regarding 5 Gy, 25 Gy isodose volumes of the normal tissue surrounding the target, there is no significant difference among the six treatment techniques. The 3D non-coplanar and SA non-coplanar plans had the lowest mean values of 5 Gy, 25 Gy volume of normal tissue respectively, while the highest values were from the SA coplanar plans and the 3D CRT coplanar plans, as displayed in Figure 8.
4. Discussion
Our study compared the 3D CRT, IMRT and SA planning techniques for SBRT of lung tumor using both cop-lanar and non-coplanar beam arrangement. Similar studies and article are available in the literature [3] [7] [10] [11] . While Holt et al. compared the efficacy of coplanar SA against IMRT for SBRT of early stage lung
Figure 7. The percentage of total lung irradiated by 5 Gy, 20 Gy isodose lines, shown as mean value ± 1 standard deviation. Notice that the differences among the 6 treatment techniques are insignificant.
Figure 8. The volume of normal tissue covered by the 5 Gy and 25 Gy isodose lines for the six treatment techniques, shown as mean values ± 1 standard deviation.
cancer [7] , Bree et al. studied delivery of higher dose to lung tumor using IMRT than that was possible with 3D-CRT [10] .
The dose distributions in Figure 1 nicely demonstrate the target dose conformity improvements with both coplanar and non-coplanar IMRT and SA than 3D CRT plans. The dose distribution in non-coplanar SA plan with 2 shorter arcs was slightly more conformal than that in coplanar SA plan with 2 complete arcs, but the improvements in CI are not significant. The conformity of the 100% isodose lines of the IMRT and SA plans were well within the clinical acceptable limit specified in RTOG 0813 protocol of CI ≤ 1.2 and upto 1.5 would be a minor deviation. 2 of the 3D-CRT plans violated the minor deviation specification. For target dose conformity, the IMRT and SA plans scored higher than 3D CRT plans, which concurred with the results from Wu et al. [12] . Values of CIs are comparable with Verbakel et al. [5] , although in contrast to the findings in that study, we did not observe any improvement in CIs with SA plans compared with that of non-coplanar IMRT plans. The number and arrangement of beams as well as the width of multileaf collimator might play a role.
Regarding R50%, both IMRT and SA plans gave tighter dose distribution than the 3D CRT plans. Specifically, R50% values of non-coplanar SA plans were significantly smaller than that of all other plans. With the exception of 2 cases, R50%values of non-coplanar SA plans were the lowest in comparison with other plans. This was expected, as non-coplanar SA spreads out the dose in non-coplanar directions while the dose deposited in the coplanar SA plan is within the plane of the arc. This is consistent with results from other similar studies [7] [8] .
Heterogeneity of target dose distribution is higher in IMRT coplanar plans than those of both 3D CRT and VMAT coplanar plans. A nasopharyngeal tumor study reported that the 7 field IMRT plans have significantly more heterogeneous dose distribution than VMAT plans [13] .
The maximum cord dose was significantly lower for the SA coplanar plans than the 3D CRT coplanar plans. This is in agreement with the findings of a study on VMAT of spinal cord [14] . The V5 of lung and normal tissue were higher with the SA plans, as expected from the distribution of dose through the plane of the arc, but well below the clinically acceptable levels. This is in line with other lung SBRT studies which had quoted similar increase to varying degrees depending on many factors including the beam arrangement [7] . The dosimetric data for the rest of the OARs were not significantly different.
The type of dose calculation algorithm and the TPS could greatly impact the dosimetric results, especially in the case of heterogeneous media and small field geometry. Pinnacle3 features the CCCS algorithm that includes scatter corrections in a heterogeneous media. The dose-calculation algorithm has been shown to result in lower conformity in the lung SBRT plans [15] . Some of the dose differences could be due to the limitations in the beam modeling of low density media in the Pinnacle3 CCCS algorithm [16] . Coupled with lack of adequate scatter contribution in the small radiation field sizes that are commonly used in SBRT beams, dose differences were observed near the interface [17] [18] .
The comparison between 3D CRT, IMRT and SA in this study is specific to SBRT of lung tumor. Limitations of this study include any systematic error due to the low patient number evaluated in this study. For IMRT, a seven-field configuration and for SA, a two-arc configuration was examined. Alternate planning techniques using more beams could be attempted. However, incorporation of such changes could increase the dosimetric quality and the plan complexity at the cost of treatment time. Long-term complication-free survival and patient outcomes are needed before definite conclusions can be made. It is advisable to analyze treatment plans of each patient to determine which plan will be the most suited. The suitability of the treatment technique could vary depending on the staging, tumor size, location of tumor, vicinity of the OARs as well as the dose tolerance criteria observed.
5. Conclusion
In this study, all six treatment techniques were able to deliver conformal SBRT plans while meeting the dose constraints of the OARs. While IMRT and SA plans scored better than 3D CRT plans with regard to dose conformity and dose fall-off, IMRT plans fared better than others in the target dose homogeneity index. SA coplanar plans provided significantly lower maximum dose to the cord. However, no one treatment technique was able to produce a plan with significantly superior plan quality indices over others for any given patient.