Precision IORT – Image guided intraoperative radiation therapy (igIORT) using online treatment planning including tissue heterogeneity correction

Development of a TLD-100 based set up for in vivo dosimetry in Intraoperative Electron Beam Radiation Therapy (IOERT): an experimental and clinical evaluation

PURPOSE

This study presents an in-house developed set-up enabling the placement of multiple TLDs over the target region, for in vivo dosimetry in intraoperative electron beam radiotherapy (IOERT).

METHODS AND MATERIALS

TLD-100 (LiF:MgTi) rods were calibrated at 10 Gy and their response was determined for the nominal electron energies of the Mobetron LINAC within the dose range from 4 Gy to 20 Gy. Irradiation of various set-ups was performed using the ionisation chambers (IC) ROOS parallel plate 34001 (PTW, Freiburg) and 3D Semiflex 31021 (PTW, Freiburg), a microdiamond detector 60019 (PTW, Freiburg) and EBT3 films (Ashland™) to investigate beam perturbations that may result due to the structure of the set up. EGSnrc Monte Carlo (MC) simulations evaluated the response of the TLDs in clinical beams of the available electron energies, the influence of the catheter to the TLD dose scoring and the depth dependence of the TLD dose assessment. TLD measurements in-phantoms and in patient in vivo were realised and compared to the expected doses estimated using data of water phantom measurements and 3D MC electron dose calculations of a dedicated IOERT treatment planning system (Radiance TPS- GmV, Tres Cantos, Madrid).

RESULTS

MC and measurements verified that no energy correction is needed for the used electron beams. Correction factors for the dose non-linear response were evaluated. High resolution dose measurements showed local hot spots beneath the flap. However, there are no significant perturbations of the electron beam or on the dose delivery to the targeted volume. MC simulations demonstrated no signal attenuation due to the catheter and 1 % effect of the depth of TLD measurement relative to the depth of calibration was noted. TLD measurements in phantom set-ups agreed with expected doses with less than 2.6 % in phantoms and by 1 % in patient in vivo.

CONCLUSIONS

Our results demonstrate the suitability of using the implemented TLD-based workflow for in vivo dosimetry purposes in the operation room (OR) environment.

IOeRT conventional and FLASH treatment planning system implementation exploiting fast GPU Monte Carlo: The case of breast cancer

Partial breast irradiation for the treatment of early-stage breast cancer patients can be performed by means of Intra Operative electron Radiation Therapy (IOeRT). One of the main limitations of this technique is the absence of a treatment planning system (TPS) that could greatly help in ensuring a proper coverage of the target volume during irradiation. An IOeRT TPS has been developed using a fast Monte Carlo (MC) and an ultrasound imaging system to provide the best irradiation strategy (electron beam energy, applicator position and bevel angle) and to facilitate the optimisation of dose prescription and delivery to the target volume while maximising the organs at risk sparing. The study has been performed in silico, exploiting MC simulations of a breast cancer treatment. Ultrasound-based input has been used to compute the absorbed dose maps in different irradiation strategies and a quantitative comparison between the different options was carried out using Dose Volume Histograms. The system was capable of exploring different beam energies and applicator positions in few minutes, identifying the best strategy with an overall computation time that was found to be completely compatible with clinical implementation. The systematic uncertainty related to tissue deformation during treatment delivery with respect to imaging acquisition was taken into account. The potential and feasibility of a GPU based full MC TPS implementation of IOeRT breast cancer treatments has been demonstrated in-silico. This long awaited tool will greatly improve the treatment safety and efficacy, overcoming the limits identified within the clinical trials carried out so far.

First clinical application of image-guided intraoperative electron radiation therapy with real time intraoperative dose calculation in recurrent rectal cancer: technical procedure

Intraoperative radiation therapy (IORT) is a radiation technique applying a single fraction with a high dose during surgery. We report the first abdomino-pelvic application of an image-guided intraoperative electron radiation therapy with intraoperative real time dose calculation based on the individual intraoperative patient anatomy. A patient suffering from locoregionally recurrent rectal cancer after treatment with neoadjuvant re-chemoradiation was chosen for this approach. After surgical removal of the recurrence, an adequate IORT applicator was placed as usual. A novel mobile imaging device (ImagingRing, MedPhoton) was positioned around the patient covering the region to be treated with the IORT-applicator in place. It allowed the acquisition of three-dimensional intraoperative cone-beam computed tomography images suitable for dose calculation using an automated scaling (heuristic object and head scatter as well as hardening corrections) of Hounsfield units. After image acquisition confirmed the correct applicator position, the images were transferred to our treatment planning system for intraoperative dose calculation. Treatment could be accomplished using the calculated dose distribution. We herein describe the details of the procedure including necessary adjustments in the typically used IORT equipment and work flow. We further discuss the pros and cons of this new approach generally overcoming a decade long limitation of IORT procedures as well as future perspectives regarding IORT treatments.

In vivodosimetry in cancer patients undergoing intraoperative radiation therapy

In vivodosimetry (IVD) is an important tool in external beam radiotherapy (EBRT) to detect major errors by assessing differences between expected and delivered dose and to record the received dose by individual patients. Also, in intraoperative radiation therapy (IORT), IVD is highly relevant to register the delivered dose. This is especially relevant in low-risk breast cancer patients since a high dose of IORT is delivered in a single fraction. In contrast to EBRT, online treatment planning based on intraoperative imaging is only under development for IORT. Up to date, two commercial treatment planning systems proposed intraoperative ultrasound or in-room cone-beam CT for real-time IORT planning. This makes IVD even more important because of the possibility for real-time treatment adaptation. Here, we summarize recent developments and applications of IVD methods for IORT in clinical practice, highlighting important contributions and identifying specific challenges such as a treatment planning system for IORT. HDR brachytherapy as a delivery technique was not considered. We add IVD for ultrahigh dose rate (FLASH) radiotherapy that promises to improve the treatment efficacy, when compared to conventional radiotherapy by limiting the rate of toxicity while maintaining similar tumour control probabilities. To date, FLASH IORT is not yet in clinical use.

Organ absorbed doses in the IORT treatment of breast cancer with the INTRABEAM device: a Monte-Carlo study

Purpose: The current prescription and the assessment of the delivered absorbed dose in intraoperative radiation therapy (IORT) with the INTRABEAM system rely mainly on depth-dose measurements in water. The accuracy of this approach is limited because tissue heterogeneity is ignored. It is also difficult to accurately determine the dose delivered to the patient experimentally as the steep dose gradient is highly sensitive to geometric errors. Our goal is to determine the dose to the target volume and the organs at risk of a clinical breast cancer patient from treatment with the system.Methods: A homogeneous water-equivalent CT dataset was derived from the preoperative CT scan of a patient by setting all materials in the patient volume as water-equivalent. This homogeneous CT data represents the current assumption of a homogenous patient, while the original CT data is considered the ground truth. An in-house Monte Carlo algorithm was used to simulate the delivered dose in both setups for a prescribed treatment dose of 20 Gy to the surface of the 3.5 cm diameter spherical applicator.Results: The doses received by 2% (D2%) of the target volume for the homogeneous and heterogeneous geometries are 16.26 Gy and 9.33 Gy, respectively. The D2% for the heart are 0.035 Gy and 0.119 Gy for the homogeneous and heterogeneous geometries, respectively. This trend is also observed for the other organs at risk.Conclusions: The assumption of a homogeneous patient overestimates the dose to the target volume and underestimates the doses to the organs at risk.

Implementation of pencil beam redefinition algorithm (PBRA) for intraoperative electron radiation therapy (IOERT) treatment planning

PURPOSE

Similar to other radiation therapy techniques, intraoperative electron radiation therapy (IOERT) can benefit from an online treatment planning system (TPS). Among all the analytical electron dose calculation algorithms, pencil beam redefinition algorithm (PBRA) has shown an acceptable accuracy in inhomogeneities. The input dataset for PBRA includes electron planar fluence, mean direction and root mean square (RMS) spread about the mean direction which had been introduced based on the conventional linear accelerator geometry in former studies. Herein, three methods for implementing PBRA for IOERT system are presented.

METHODS

The initialization parameters were identified using Monte Carlo (MC) simulation of a dedicated IOERT system equipped with a cylindrical 10 cm applicator, irradiating a water phantom. Phase space distribution of electrons was recorded on a plane below the applicator. The input dataset was extracted for 2 × 2 mm² pixels and energy bin width of 1 MeV.

RESULTS

PBRA was implemented with three initialization methods and compared to MC. The 3D gamma analysis of the algorithm with the Formula method, which was in best agreement with MC in a simple water phantom, showed passing rates of more than 99 % for all nominal energies and it was 97.1 % for 8 MeV in the presence of protecting disk and irregular surface. Implementing PBRA on CUDA C++ resulted in 5 s run time for 8 MeV nominal energy in a water phantom.

CONCLUSIONS

The agreement between PBRA dose calculation and MC is promising for the development of an intraoperative TPS for IOERT.

Treatment Planning in Intraoperative Radiation Therapy (IORT): Where Should We Go?

As opposed to external beam radiation therapy (EBRT), treatment planning systems (TPS) dedicated to intraoperative radiation therapy (IORT) were not subject to radical modifications in the last two decades. However, new treatment regimens such as ultrahigh dose rates and combination with multiple treatment modalities, as well as the prospected availability of dedicated in-room imaging, call for important new features in the next generation of treatment planning systems in IORT. Dosimetric accuracy should be guaranteed by means of advanced dose calculation algorithms, capable of modelling complex scattering phenomena and accounting for the non-tissue equivalent materials used to shape and compensate electron beams. Kilovoltage X-ray based IORT also presents special needs, including the correct description of extremely steep dose gradients and the accurate simulation of applicators. TPSs dedicated to IORT should also allow real-time imaging to be used for treatment adaptation at the time of irradiation. Other features implemented in TPSs should include deformable registration and capability of radiobiological planning, especially if unconventional irradiation schemes are used. Finally, patient safety requires that the multiple features be integrated in a comprehensive system in order to facilitate control of the whole process.

A Comparison between Electron Gamma Shower, National Research Council/Easy Particle Propagation (EGSnrc/Epp) and Monte Carlo N-Particle Transport Code (MCNP) in Simulation of the INTRABEAM ® System with Spherical Applicators

BACKGROUND

Online Monte Carlo (MC) treatment planning is very crucial to increase the precision of intraoperative radiotherapy (IORT). However, the performance of MC methods depends on the geometries and energies used for the problem under study.

OBJECTIVE

This study aimed to compare the performance of MC N-Particle Transport Code version 4c (MCNP4c) and Electron Gamma Shower, National Research Council/easy particle propagation (EGSnrc/Epp) MC codes using similar geometry of an INTRABEAM^® system.

MATERIAL AND METHODS

This simulation study was done by increasing the number of particles and compared the performance of MCNP4c and EGSnrc/Epp simulations using an INTRABEAM^® system with 1.5 and 5 cm diameter spherical applicators. A comparison of these two codes was done using simulation time, statistical uncertainty, and relative depth-dose values obtained after doing the simulation by each MC code.

RESULTS

The statistical uncertainties for the MCNP4c and EGSnrc/Epp MC codes were below 2% and 0.5%, respectively. 1e9 particles were simulated in 117.89 hours using MCNP4c but a much greater number of particles (5e10 particles) were simulated in a shorter time of 90.26 hours using EGSnrc/Epp MC code. No significant deviations were found in the calculated relative depth-dose values for both in the presence and absence of an air gap between MCNP4c and EGSnrc/Epp MC codes. Nevertheless, the EGSnrc/Epp MC code was found to be speedier and more efficient to achieve accurate statistical precision than MCNP4c.

CONCLUSION

Therefore, in all comparisons criteria used, EGSnrc/Epp MC code is much better than MCNP4c MC code for simulating an INTRABEAM^® system.

Optimized method for in vivo dosimetry with small films in pelvic IOERT for rectal cancer

PURPOSE

Intra-Operative Electron Radiation Therapy (IOERT) is used to treat rectal cancer at our institution, and in vivo measurements with Gafchromic EBT3® films were introduced as quality assurance. The purpose of this work was to quantify the uncertainties associated with digitization of very small EBT3 films irradiated simultaneously, in order to optimize in vivo dosimetry for IOERT.

METHODS

Film samples of different sizes - M1 (5×5cm²), M2 (1.5×1.5 cm²), M3 (1.0×1.5 cm²) and M4 (0.75×1.5 cm²) - were used to quantify typical variations (uncertainties) due to scanner fluctuations, misalignment, film inhomogeneity, long-term effect of film cutting, small rotations, film curling, edge effects and the influence of opaque templates. Fitting functions and temporal validity of sensitometric curves were also assessed.

RESULTS

Film curling, intra-film variability and scanner fluctuations are important effects that need to be minimized or considered in the uncertainty budget. Small rotations, misalignments and film cutting have little or no influence on the readings. Most fitting functions perform well, but the quantity used for dose quantification determines over- or under-valuation of dose in the long term. Edge effects and the influence of opaque templates need to be well understood, to allow optimization of methodology to the intended purpose.

CONCLUSION

The proposed method allows practical and simultaneous digitization of up to ten small irradiated film samples, with an experimental uncertainty of 1%.

Dosimetric characteristics of the INTRABEAM ® system with spherical applicators in the presence of air gaps and tissue heterogeneities

The main aim of this study was to investigate the dosimetric characteristics of the INTRABEAM ® system in the presence of air gaps between the surface of applicators (APs) and tumor bed. Additionally, the effect of tissue heterogeneities was another focus. Investigating the dosimetric characteristics of the INTRABEAM® system is essential to deliver the required dose to the tumor bed correctly and reduce the delivered dose to the ribs and lung. Choosing the correct AP size and fitting it to the lumpectomy cavity is essential to remove the effect of air gaps and avoid inaccurate dose delivery. Consequently, the Geant4 toolkit was used to simulate the INTRABEAM ® system with spherical APs of various sizes. The wall effect of the ion chamber (IC) PTW 34013 used in the present study was checked. The simulations were validated in comparison with measurements, and then used to calculate any inaccuracies in dose delivery in the presence of 4- and 10-mm air gaps between the surface of the APs and the tumor bed. Also, the doses received due to tissue heterogeneities were characterized. It turned out that measurements and simulations were approximately in agreement (± 2%) for all sizes of APs. The perturbation factor introduced by the IC due to differences in graphite-coated polyethylene and air as compared to the phantom material was approximately equal to one for all AP. The greatest relative dose delivery difference was observed for an AP with a diameter of 1.5 cm, i.e., 44% and 70% in the presence of 4- and 10-mm air gaps, respectively. In contrast, the lowest relative dose delivery difference was observed for an AP with a diameter of 5 cm, i.e., 24% and 42% in the presence of 4- and 10-mm air gaps, respectively. Increasing APs size showed a decrease in relative dose delivery difference due to the presence of air gaps. In addition, the undesired dose received by the ribs turned out to be higher when a treatment site closer to the ribs was assumed. The undesired dose received by the ribs increased as the AP size increased. The lung dose turned out to be decreased due to the shielding effect of the ribs, small lung density, and long separation distance from the AP surface.