Quantifying dose perturbations in high-risk prostate radiotherapy due to translational and rotational motion of prostate and pelvic lymph nodes

Background

Radiotherapy of the prostate and the pelvic lymph nodes (LN) is a part of the standard of care treatment for high-risk prostate cancer. The independent translational and rotational (i.e., six-degrees-of-freedom, [6DoF]) motion of the prostate and LN target during and between fractions can perturb the dose distribution. However, no standard dose reconstruction method accounting for differential 6DoF target motion is available.

Purpose

We present a framework for monitoring motion-induced dose perturbations for two independently moving target volumes in 6DoF. The framework was used to determine the dose perturbation for the prostate and the LN target caused by differential 6DoF motion for a cohort of high-risk prostate cancer patients. As a potential first step toward real-time dose-guided high-risk prostate radiotherapy, we furthermore investigated if the dose reconstruction was fast enough for real-time application for both targets.

Methods

Twenty high-risk prostate cancer patients were treated with 3-arc volumetric modulated arc therapy (VMAT). Kilovoltage intrafraction monitoring (KIM) with triggered kilovoltage (kV) images acquired every 3 throughout 7–10 fractions per patient was used for retrospective 6DoF intrafraction prostate motion estimation. The 6DoF interfraction LN motion was determined from a pelvic bone match between the planning CT and a post-treatment cone beam CT (CBCT). Using the retrospectively extracted motion, real-time 6DoF motion-including dose reconstruction was simulated using the in-house developed software DoseTracker. A data stream with the 6DoF target positions and linac parameters was broadcasted at a 3-Hz frequency to DoseTracker. In a continuous loop, DoseTracker calculated the target dose increments including the specified motion and, for comparison, without motion. The motion-induced change in D_99.5% for the prostate CTV (ΔD_99.5%) and in D_98% for the LN CTV (ΔD_98%) was calculated using the final cumulative dose of each fraction and averaged over all imaged fractions. The real-time reconstructed dose distribution of DoseTracker was benchmarked against a clinical treatment planning system (TPS) and it was investigated whether the calculation speed was fast enough to keep up with the incoming data stream.

Results

Translational motion was largest in cranio-caudal (CC) direction (prostate: [-5.9, +8.4] mm; LN: [-9.9; +11.0] mm) and anterior–posterior (AP) direction (prostate:[-5.6; +6.9] mm; LN: [-9.6; +11.0] mm). The pitch was the largest rotation (prostate: [-22.5; +25.2] deg; LN: [-3.9; +5.5] deg). The prostate CTV ΔD_99.5% was [-16.2; +2.5]% for single fractions and [-3.0; +1.7]% when averaged over all imaged fractions. The LN CTV ΔD_98% was [-19.8; +1.2]% for single fractions and [-3.1; +0.9]% after averaging. Mean (Standard deviation) absolute dose errors in DoseTracker of 107.8% (Std: 1.9%) for the prostate and 105.5% (Std:1.4%) for the LN were corrected during dose reconstruction by automatically calculated normalization factors. It resulted in accurate calculation of the motion-induced dose errors with relative differences between DoseTracker and TPS dose calculations of -0.1% (Std: 0.5%) (prostate CTV ΔD_99.5%) and -0.2% (Std: 0.5%) (LN CTV ΔD_98%). The DoseTracker calculation was fast enough to keep up with the incoming inputs for all but two out of 107 184 dose calculations.

Conclusion

Using the developed framework for dose perturbation monitoring, we found that the differential 6DoF target motion caused substantial dose perturbation for individual fractions, which largely averaged out after several fractions. The framework was shown to provide reliable dose calculations and a sufficiently high-dose reconstruction speed to be applicable in real-time.