The impact of breathing motion versus heterogeneity effects in lung cancer treatment planning

Impact of Intra-Fractional Motion on Dose Distributions in Lung IMRT

AIM

To investigate the impact of intra-fractional motion on dose distribution in patients treated with intensity-modulated radiation therapy (IMRT) for lung cancer.

MATERIALS AND METHODS

Twenty patients who had undergone IMRT for non-small cell lung cancer were selected for this retrospective study. For each patient, a four-dimensional computed tomography (CT) image set was acquired and clinical treatment plans were developed using the average CT. Dose distributions were then re-calculated for each of the 10 phases of respiratory cycle and combined using deformable image registration to produce cumulative dose distributions that were compared with the clinical treatment plans.

RESULTS

Intra-fractional motion reduced planning target volume (PTV) coverage in all patients. The median reduction of PTV volume covered by the prescription isodose was 3.4%; D98 was reduced by 3.1 Gy. Changes in the mean lung dose were within ±0.7 Gy. V20 for the lung increased in most patients; the median increase was 1.6%. The dose to the spinal cord was unaffected by intra-fractional motion. The dose to the heart was slightly reduced in most patients. The median reduction in the mean heart dose was 0.22 Gy, and V30 was reduced by 2.5%.The maximum dose to the esophagus was also reduced in most patients, by 0.74 Gy, whereas V50 did not change significantly. The median number of points in which dose differences exceeded 3%/3 mm was 6.2%.

FINDINGS

Intra-fractional anatomical changes reduce PTV coverage compared to the coverage predicted by clinical treatment planning systems that use the average CT for dose calculation. Doses to organs at risk were mostly over-predicted.

Towards accurate dose accumulation for Step-&-Shoot IMRT: Impact of weighting schemes and temporal image resolution on the estimation of dosimetric motion effects

PURPOSE

Breathing-induced motion effects on dose distributions in radiotherapy can be analyzed using 4D CT image sequences and registration-based dose accumulation techniques. Often simplifying assumptions are made during accumulation. In this paper, we study the dosimetric impact of two aspects which may be especially critical for IMRT treatment: the weighting scheme for the dose contributions of IMRT segments at different breathing phases and the temporal resolution of 4D CT images applied for dose accumulation.

METHODS

Based on a continuous problem formulation a patient- and plan-specific scheme for weighting segment dose contributions at different breathing phases is derived for use in step-&-shoot IMRT dose accumulation. Using 4D CT data sets and treatment plans for 5 lung tumor patients, dosimetric motion effects as estimated by the derived scheme are compared to effects resulting from a common equal weighting approach. Effects of reducing the temporal image resolution are evaluated for the same patients and both weighting schemes.

RESULTS

The equal weighting approach underestimates dosimetric motion effects when considering single treatment fractions. Especially interplay effects (relative misplacement of segments due to respiratory tumor motion) for IMRT segments with only a few monitor units are insufficiently represented (local point differences >25% of the prescribed dose for larger tumor motion). The effects, however, tend to be averaged out over the entire treatment course. Regarding temporal image resolution, estimated motion effects in terms of measures of the CTV dose coverage are barely affected (in comparison to the full resolution) when using only half of the original resolution and equal weighting. In contrast, occurence and impact of interplay effects are poorly captured for some cases (large tumor motion, undersized PTV margin) for a resolution of 10/14 phases and the more accurate patient- and plan-specific dose accumulation scheme.

CONCLUSIONS

Radiobiological consequences of reported single fraction local point differences >25% of the prescribed dose are widely unclear and should be subject to future investigation. Meanwhile, if aiming at accurate and reliable estimation of dosimetric motion effects, precise weighting schemes such as the presented patient- and plan-specific scheme for step-&-shoot IMRT and full available temporal 4D CT image resolution should be applied for IMRT dose accumulation.

Extension of the NCAT phantom for the investigation of intra-fraction respiratory motion in IMRT using 4D Monte Carlo

The purpose of this work was to create a computational platform for studying motion in intensity modulated radiotherapy (IMRT). Specifically, the non-uniform rational B-spline (NURB) cardiac and torso (NCAT) phantom was modified for use in a four-dimensional Monte Carlo (4D-MC) simulation system to investigate the effect of respiratory-induced intra-fraction organ motion on IMRT dose distributions as a function of diaphragm motion, lesion size and lung density. Treatment plans for four clinical scenarios were designed: diaphragm peak-to-peak amplitude of 1 cm and 3 cm, and two lesion sizes--2 cm and 4 cm diameter placed in the lower lobe of the right lung. Lung density was changed for each phase using a conservation of mass calculation. Further, a new heterogeneous lung model was implemented and tested. Each lesion had an internal target volume (ITV) subsequently expanded by 15 mm isotropically to give the planning target volume (PTV). The PTV was prescribed to receive 72 Gy in 40 fractions. The MLC leaf sequence defined by the planning system for each patient was exported and used as input into the MC system. MC simulations using the dose planning method (DPM) code together with deformable image registration based on the NCAT deformation field were used to find a composite dose distribution for each phantom. These composite distributions were subsequently analyzed using information from the dose volume histograms (DVH). Lesion motion amplitude has the largest effect on the dose distribution. Tumor size was found to have a smaller effect and can be mitigated by ensuring the planning constraints are optimized for the tumor size. The use of a dynamic or heterogeneous lung density model over a respiratory cycle does not appear to be an important factor with a <or=0.6% change in the mean dose received by the ITV, PTV and right lung. The heterogeneous model increases the realism of the NCAT phantom and may provide more accurate simulations in radiation therapy investigations that use the phantom. This work further evaluates the NCAT phantom for use as a tool in radiation therapy research in addition to its extensive use in diagnostic imaging and nuclear medicine research. Our results indicate that the NCAT phantom, combined with 4D-MC simulations, is a useful tool in radiation therapy investigations and may allow the study of relative effects in many clinically relevant situations.

Percentage depth dose evaluation in heterogeneous media using thermoluminescent dosimetry

The purpose of this study is to investigate the influence of lung heterogeneity inside a soft tissue phantom on percentage depth dose (PDD). PDD curves were obtained experimentally using LiF:Mg,Ti (TLD‐100) thermoluminescent detectors and applying Eclipse treatment planning system algorithms Batho, modified Batho (M‐Batho or BMod), equivalent TAR (E‐TAR or EQTAR), and anisotropic analytical algorithm (AAA) for a 15 MV photon beam and field sizes of 1×1,2×2,5×5, and 10×10cm2. Monte Carlo simulations were performed using the DOSRZnrc user code of EGSnrc. The experimental results agree with Monte Carlo simulations for all irradiation field sizes. Comparisons with Monte Carlo calculations show that the AAA algorithm provides the best simulations of PDD curves for all field sizes investigated. However, even this algorithm cannot accurately predict PDD values in the lung for field sizes of 1×1 and 2×2cm2. An overdosage in the lung of about 40% and 20% is calculated by the AAA algorithm close to the interface soft tissue/lung for 1×1 and 2×2cm2 field sizes, respectively. It was demonstrated that differences of 100% between Monte Carlo results and the algorithms Batho, modified Batho, and equivalent TAR responses may exist inside the lung region for the 1×1cm2 field.

PACS number: 87.55.kd

Design and testing of a simulation framework for dosimetric motion studies integrating an anthropomorphic computational phantom into four-dimensional Monte Carlo

We have designed a simulation framework for motion studies in radiation therapy by integrating the anthropomorphic NCAT phantom into a 4D Monte Carlo dose calculation engine based on DPM. Representing an artifact-free environment, the system can be used to identify class solutions as a function of geometric and dosimetric parameters. A pilot dynamic conformal study for three lesions ( approximately 2.0 cm) in the right lung was performed (70 Gy prescription dose). Tumor motion changed as a function of tumor location, according to the anthropomorphic deformable motion model. Conformal plans were simulated with 0 to 2 cm margin for the aperture, with additional 0.5 cm for beam penumbra. The dosimetric effects of intensity modulated radiotherapy (IMRT) vs. conformal treatments were compared in a static case. Results show that the Monte Carlo simulation framework can model tumor tracking in deformable anatomy with high accuracy, providing absolute doses for IMRT and conformal radiation therapy. A target underdosage of up to 3.67 Gy (lower lung) was highlighted in the composite dose distribution mapped at exhale. Such effects depend on tumor location and treatment margin and are affected by lung deformation and ribcage motion. In summary, the complexity in the irradiation of moving targets has been reduced to a controlled simulation environment, where several treatment options can be accurately modeled and quantified The implemented tools will be utilized for extensive motion study in lung/liver irradiation.

Dosimetric impact of motion in free-breathing and gated lung radiotherapy: a 4D Monte Carlo study of intrafraction and interfraction effects

The purpose of this study was to investigate if interfraction and intrafraction motion in free-breathing and gated lung IMRT can lead to systematic dose differences between 3DCT and 4DCT. Dosimetric effects were studied considering the breathing pattern of three patients monitored during the course of their treatment and an in-house developed 4D Monte Carlo framework. Imaging data were taken in free-breathing and in cine mode for both 3D and 4D acquisition. Treatment planning for IMRT delivery was done based on the free-breathing data with the CORVUS (North American Scientific, Chatsworth, CA) planning system. The dose distributions as a function of phase in the breathing cycle were combined using deformable image registration. The study focused on (a) assessing the accuracy of the CORVUS pencil beam algorithm with Monte Carlo dose calculation in the lung, (b) evaluating the dosimetric effect of motion on the individual breathing phases of the respiratory cycle, and (c) assessing intrafraction and interfraction motion effects during free-breathing or gated radiotherapy. The comparison between (a) the planning system and the Monte Carlo system shows that the pencil beam algorithm underestimates the dose in low-density regions, such as lung tissue, and overestimates the dose in high-density regions, such as bone, by 5% or more of the prescribed dose (corresponding to approximately 3-5 Gy for the cases considered). For the patients studied this could have a significant impact on the dose volume histograms for the target structures depending on the margin added to the clinical target volume (CTV) to produce either the planning target (PTV) or internal target volume (ITV). The dose differences between (b) phases in the breathing cycle and the free-breathing case were shown to be negligible for all phases except for the inhale phase, where an underdosage of the tumor by as much as 9.3 Gy relative to the free-breathing was observed. The large difference was due to breathing-induced motion/deformation affecting the soft/lung tissue density and motion of the bone structures (such as the rib cage) in and out of the beam. Intrafraction and interfraction dosimetric differences between (c) free-breathing and gated delivery were found to be small. However, more significant dosimetric differences, of the order of 3%-5%, were observed between the dose calculations based on static CT (3DCT) and the ones based on time-resolved CT (4DCT). These differences are a consequence of the larger contribution of the inhale phase in the 3DCT data than in the 4DCT.