Dosimetric analysis of local failures in skull-base chordoma and chondrosarcoma following pencil beam scanning proton therapy

Background

Despite combined modality treatment involving surgery and radiotherapy, a relevant proportion of skull-base chordoma and chondrosarcoma patients develop a local recurrence (LR). This study aims to analyze patterns of recurrence and correlate LR with a detailed dosimetric analysis.

Methods

222 patients were treated with proton radiotherapy for chordoma (n = 151) and chondrosarcoma (n = 71) at the PSI between 1998 and 2012. All patients underwent surgery, followed by pencil-beam scanning proton therapy to a mean dose of 72.5 ± 2.2Gy_RBE. A retrospective patterns of recurrence analysis was performed: LR were contoured on follow-up MRI, registered with planning-imaging and the overlap with initial target structures (GTV, PTV_high-dose, PTV_low-dose) was calculated. DVH parameters of planning structures and recurrences were calculated and correlated with LR using univariate and multivariate cox regression.

Results

After a median follow-up of 50 months, 35 (16%) LR were observed. Follow-up MRI imaging was available for 27 (77%) of these recurring patients. Only one (3.7%) recurrence was located completely outside the initial PTV (surgical pathway recurrence). The mean proportions of LR covered by the initial target structures were 48% (range 0–86%) for the GTV, 70% (range 0–100%) for PTV_high and 83% (range 0–100%) for PTV_low. In the univariate analysis, the following DVH parameters were significantly associated with LR: GTV(V < 66Gy_RBE, p = 0.01), GTV(volume, p = 0.02), PTV_high(max, p = 0.02), PTV_high(V < 66Gy_RBE, p = 0.03), PTV_high(V < 59Gy_RBE, p = 0.02), PTV_high(volume, p = 0.01) and GTV(D95, p = 0.05). In the multivariate analysis, only histology (chordoma vs. chondrosarcoma, p = 0.01), PTV_high(volume, p = 0.05) and GTV(V < 66Gy_RBE, p = 0.02) were independent prognostic factors for LR.

Conclusion

This study identified DVH parameters, which are associated with the risk of local recurrence after proton therapy using pencil-beam scanning for patients with skull-base chordoma and chondrosarcoma.

Anthropomorphic phantom for deformable lung and liver CT and MR imaging for radiotherapy

In this study, a functioning and ventilated anthropomorphic phantom was further enhanced for the purpose of CT and MR imaging of the lung and liver. A deformable lung, including respiratory tract was 3D printed. Within the lung's inner structures is a solid region shaped from a patient's lung tumour and six nitro-glycerine capsules as reference landmarks. The full internal mesh was coated, and the tumour filled, with polyorganosiloxane based gel. A moulded liver was created with an external casing of silicon filled with polyorganosiloxane gel and flexible plastic internal structures. The liver, fitted to the inferior portion of the right lung, moves along with the lung's ventilation. In the contralateral side, a cavity is designed to host a dosimeter, whose motion is correlated to the lung pressure. A 4DCT of the phantom was performed along with static and 4D T1 weighted MR images. The CT Hounsfield units (HU) for the flexible 3D printed material were -600-100 HU (lung and liver structures), for the polyorganosiloxane gel 30-120 HU (lung coating and liver filling) and for the silicon 650-800 HU (liver casing). The MR image intensity units were 0-40, 210-280 and 80-130, respectively. The maximum range of motion in the 4D imaging for the superior lung was 1-3.5 mm and 3.5-8 mm in the inferior portion. The liver motion was 5.5-8.0 mm at the tip and 5.7-10.0 mm at the dome. No measurable drift in motion was observed over a 2 h session and motion was reproducible over three different sessions for sin²(t), sin⁴(t) and a patient-like breathing curve with the interquartile range of amplitudes for all breathing cycles within 0.5 mm. The addition of features within the lung and of a deformable liver will allow the phantom to be used for imaging studies such as validation of 4DMRI and pseudo CT methods.

Technical Note: Benchmarking automated eye tracking and human detection for motion monitoring in ocular proton therapy

PURPOSE

Ocular proton therapy is an effective therapeutic option for patients affected with uveal melanomas. An optical eye-tracking system (ETS) aiming at noninvasive motion monitoring was developed and tested in a clinical scenario.

MATERIALS AND METHODS

The ETS estimates eye position and orientation at 25 frames per second using the three-dimensional position of pupil and cornea curvature centers identified, in the treatment room, through stereoscopic optical imaging and infrared eye illumination. Its capabilities for automatic detection of eye motion were retrospectively evaluated on 60 treatment fractions. Then, the ETS performance was benchmarked against the clinical standard based on visual control and manual beam interruption.

RESULTS

Eye-tracking system detected eye position successfully in 97% of all available frames. Eye-tracking system-based eye monitoring during therapy guarantees quicker response to involuntary eye motions than manual beam interruptions and avoids unnecessary beam interruptions.

CONCLUSIONS

Eye-tracking system shows promise for on-line monitoring of eye motion. Its introduction in the clinical workflow will guarantee a swifter treatment course for the patient and the clinical personnel.

Online daily adaptive proton therapy

It is recognized that the use of a single plan calculated on an image acquired some time before the treatment is generally insufficient to accurately represent the daily dose to the target and to the organs at risk. This is particularly true for protons, due to the physical finite range. Although this characteristic enables the generation of steep dose gradients, which is essential for highly conformal radiotherapy, it also tightens the dependency of the delivered dose to the range accuracy. In particular, the use of an outdated patient anatomy is one of the most significant sources of range inaccuracy, thus affecting the quality of the planned dose distribution. A plan should be ideally adapted as soon as anatomical variations occur, ideally online. In this review, we describe in detail the different steps of the adaptive workflow and discuss the challenges and corresponding state-of-the art developments in particular for an online adaptive strategy.

Evaluation of the ray-casting analytical algorithm for pencil beam scanning proton therapy

For pencil beam scanned (PBS) proton therapy, analytical dose calculation engines are still typically used for the optimisation process, and often for the final evaluation of the plan. Recently however, the suitability of analytical calculations for planning PBS treatments has been questioned. Conceptually, the two main approaches for these analytical dose calculations are the ray-casting (RC) and the pencil-beam (PB) method. In this study, we compare dose distributions and dosimetric indices, calculated on both the clinical dose calculation grid and as a function of dose grid resolution, to Monte Carlo (MC) calculations. The analysis is done using a comprehensive set of clinical plans which represent a wide choice of treatment sites. When analysing dose difference histograms for relative treatment plans, pencil beam calculations with double grid resolution perform best, with on average 97.7%/91.9% (RC), 97.9%/92.7% (RC, double grid resolution), 97.6%/91.0% (PB) and 98.6%/94.0% (PB, double grid resolution) of voxels agreeing within  ±5%/±  3% between the analytical and the MC calculations. Even though these point-to-point dose comparison shows differences between analytical and MC calculations, for all algorithms, clinically relevant dosimetric indices agree within  ±4% for the PTV and within  ±5% for critical organs. While the clinical agreement depends on the treatment site, there is no substantial difference of indices between the different algorithms. The pencil-beam approach however comes at a higher computational cost than the ray-casting calculation. In conclusion, we would recommend using the ray-casting algorithm for fast dose optimization and subsequently combine it with one MC calculation to scale the absolute dose and assure the quality of the treatment plan.

Gas and seismicity within the Istanbul seismic gap

Understanding micro-seismicity is a critical question for earthquake hazard assessment. Since the devastating earthquakes of Izmit and Duzce in 1999, the seismicity along the submerged section of North Anatolian Fault within the Sea of Marmara (comprising the “Istanbul seismic gap”) has been extensively studied in order to infer its mechanical behaviour (creeping vs locked). So far, the seismicity has been interpreted only in terms of being tectonic-driven, although the Main Marmara Fault (MMF) is known to strike across multiple hydrocarbon gas sources. Here, we show that a large number of the aftershocks that followed the M 5.1 earthquake of July, 25^th 2011 in the western Sea of Marmara, occurred within a zone of gas overpressuring in the 1.5–5 km depth range, from where pressurized gas is expected to migrate along the MMF, up to the surface sediment layers. Hence, gas-related processes should also be considered for a complete interpretation of the micro-seismicity (~M < 3) within the Istanbul offshore domain.

A study of lateral fall-off (penumbra) optimisation for pencil beam scanning (PBS) proton therapy

The lateral fall-off is crucial for sparing organs at risk in proton therapy. It is therefore of high importance to minimize the penumbra for pencil beam scanning (PBS). Three optimisation approaches are investigated: edge-collimated uniformly weighted spots (collimation), pencil beam optimisation of uncollimated pencil beams (edge-enhancement) and the optimisation of edge collimated pencil beams (collimated edge-enhancement). To deliver energies below 70 MeV, these strategies are evaluated in combination with the following pre-absorber methods: field specific fixed thickness pre-absorption (fixed), range specific, fixed thickness pre-absorption (automatic) and range specific, variable thickness pre-absorption (variable). All techniques are evaluated by Monte Carlo simulated square fields in a water tank. For a typical air gap of 10 cm, without pre-absorber collimation reduces the penumbra only for water equivalent ranges between 4-11 cm by up to 2.2 mm. The sharpest lateral fall-off is achieved through collimated edge-enhancement, which lowers the penumbra down to 2.8 mm. When using a pre-absorber, the sharpest fall-offs are obtained when combining collimated edge-enhancement with a variable pre-absorber. For edge-enhancement and large air gaps, it is crucial to minimize the amount of material in the beam. For small air gaps however, the superior phase space of higher energetic beams can be employed when more material is used. In conclusion, collimated edge-enhancement combined with the variable pre-absorber is the recommended setting to minimize the lateral penumbra for PBS. Without collimator, it would be favourable to use a variable pre-absorber for large air gaps and an automatic pre-absorber for small air gaps.