Respiratory motion management in particle therapy

A review of the clinical introduction of 4D particle therapy research concepts

Background and purpose

Many 4D particle therapy research concepts have been recently translated into clinics, however, remaining substantial differences depend on the indication and institute-related aspects. This work aims to summarise current state-of-the-art 4D particle therapy technology and outline a roadmap for future research and developments.

Material and methods

This review focused on the clinical implementation of 4D approaches for imaging, treatment planning, delivery and evaluation based on the 2021 and 2022 4D Treatment Workshops for Particle Therapy as well as a review of the most recent surveys, guidelines and scientific papers dedicated to this topic.

Results

Available technological capabilities for motion surveillance and compensation determined the course of each 4D particle treatment. 4D motion management, delivery techniques and strategies including imaging were diverse and depended on many factors. These included aspects of motion amplitude, tumour location, as well as accelerator technology driving the necessity of centre-specific dosimetric validation. Novel methodologies for X-ray based image processing and MRI for real-time tumour tracking and motion management were shown to have a large potential for online and offline adaptation schemes compensating for potential anatomical changes over the treatment course. The latest research developments were dominated by particle imaging, artificial intelligence methods and FLASH adding another level of complexity but also opportunities in the context of 4D treatments.

Conclusion

This review showed that the rapid technological advances in radiation oncology together with the available intrafractional motion management and adaptive strategies paved the way towards clinical implementation.

Characteristics of breathing-adapted gating using surface guidance for use in particle therapy: A phantom-based end-to-end test from CT simulation to dose delivery

To account for intra-fractional tumor motion during dose delivery in radiotherapy, various treatment strategies are clinically implemented such as breathing-adapted gating and irradiating the tumor during specific breathing phases. In this work, we present a comprehensive phantom-based end-to-end test of breathing-adapted gating utilizing surface guidance for use in particle therapy. A commercial dynamic thorax phantom was used to reproduce regular and irregular breathing patterns recorded by the GateRT respiratory monitoring system. The amplitudes and periods of recorded breathing patterns were analysed and compared to planned patterns (ground-truth). In addition, the mean absolute deviations (MAD) and Pearson correlation coefficients (PCC) between the measurements and ground-truth were assessed. Measurements of gated and non-gated irradiations were also analysed with respect to dosimetry and geometry, and compared to treatment planning system (TPS). Further, the latency time of beam on/off was evaluated. Compared to the ground-truth, measurements performed with GateRT showed amplitude differences between 0.03 ± 0.02 mm and 0.26 ± 0.03 mm for regular and irregular breathing patterns, whilst periods of both breathing patterns ranged with a standard deviation between 10 and 190 ms. Furthermore, the GateRT software precisely acquired breathing patterns with a maximum MAD of 0.30 ± 0.23 mm. The PCC constantly ranged between 0.998 and 1.000. Comparisons between TPS and measured dose profiles indicated absolute mean dose deviations within institutional tolerances of ±5%. Geometrical beam characteristics also varied within our institutional tolerances of 1.5 mm. The overall time delays were <60 ms and thus within both recommended tolerances published by ESTRO and AAPM of 200 and 100 ms, respectively. In this study, a non-invasive optical surface-guided workflow including image acquisition, treatment planning, patient positioning and gated irradiation at an ion-beam gantry was investigated, and shown to be clinically viable. Based on phantom measurements, our results show a clinically-appropriate spatial, temporal, and dosimetric accuracy when using surface guidance in the clinical setting, and the results comply with international and institutional guidelines and tolerances.

Selection of breathing phase number in 4D scanned proton treatment planning optimization for lung tumors

The standard four-dimensional (4D) treatment planning includes all breathing states in the optimization process, which is time-consuming. This work was aimed to optimize the number of intermediate phases needed for 4D proton treatment planning optimization to reduce the computational cost. Five 4D optimization strategies adopting different numbers of intermediate states and one three-dimensional (3D) optimization plan were studied for fifteen lung cancer patients treated with scanned protons, optimizing on all ten phases (4D_10), two extreme phases (4D_2), six phases during the exhalation stage (4D_6EX), six phases during the inhalation stage (4D_6IN), two extreme phases plus an intermediate state (4D_3) and average computed tomography image (3D), respectively. The 4D dose evaluation was conducted on all the ten phases, considering the interplay effect. The resulting doses accumulated on the reference phase were computed and compared. Compared to the 4D optimization plans, the 3D optimization plan performed inferiorly in target coverage, but superiorly in organ at risks (OARs) sparing. For the 4D optimization, all the five 4D plans showed similar performance in OARs protection. However, the 4D_6EX and 4D_6IN strategies out-performed the 4D_2 and 4D_3 plans in dose homogeneity. The computing times of the 4D_2, 4D_3, 4D_6EX and 4D_6IN approaches decreased to 32%, 41%, 66% and 67% of the 4D_10 method, respectively. Thus, our study suggested that the use of all phases during inhalation or exhalation stage might be a feasible approach substituting for the full phase strategy to reduce the calculation load while guaranteeing the plan quality for scanned proton therapy.

Patient-specific three-dimensional image reconstruction from a single X-ray projection using a convolutional neural network for on-line radiotherapy applications

Background and purpose: Radiotherapy is commonly chosen to treat thoracic and abdominal cancers. However, irradiating mobile tumors accurately is extremely complex due to the organs' breathing-related movements. Different methods have been studied and developed to treat mobile tumors properly. The combination of X-ray projection acquisition and implanted markers is used to locate the tumor in two dimensions (2D) but does not provide three-dimensional (3D) information. The aim of this work is to reconstruct a high-quality 3D computed tomography (3D-CT) image based on a single X-ray projection to locate the tumor in 3D without the need for implanted markers. Materials and Methods: Nine patients treated for a lung or liver cancer in radiotherapy were studied. For each patient, a data augmentation tool was used to create 500 new 3D-CT images from the planning four-dimensional computed tomography (4D-CT). For each 3D-CT, the corresponding digitally reconstructed radiograph was generated, and the 500 2D images were input into a convolutional neural network that then learned to reconstruct the 3D-CT. The dice score coefficient, normalized root mean squared error and difference between the ground-truth and the predicted 3D-CT images were computed and used as metrics. Results: Metrics' averages across all patients were 85.5% and 96.2% for the gross target volume, 0.04 and 0.45 Hounsfield unit (HU), respectively. Conclusions: The proposed method allows reconstruction of a 3D-CT image from a single digitally reconstructed radiograph that could be used in real-time for better tumor localization and improved treatment of mobile tumors without the need for implanted markers.

Intrafraction tumor motion monitoring and dose reconstruction for liver pencil beam scanning proton therapy

Background

Pencil beam scanning (PBS) proton therapy can provide highly conformal target dose distributions and healthy tissue sparing. However, proton therapy of hepatocellular carcinoma (HCC) is prone to dosimetrical uncertainties induced by respiratory motion. This study aims to develop intra-treatment tumor motion monitoring during respiratory gated proton therapy and combine it with motion-including dose reconstruction to estimate the delivered tumor doses for individual HCC treatment fractions.

Methods

Three HCC-patients were planned to receive 58 GyRBE (n=2) or 67.5 GyRBE (n=1) of exhale respiratory gated PBS proton therapy in 15 fractions. The treatment planning was based on the exhale phase of a 4-dimensional CT scan. Daily setup was based on cone-beam CT (CBCT) imaging of three implanted fiducial markers. An external marker block (RPM) on the patient's abdomen was used for exhale gating in free breathing. This study was based on 5 fractions (patient 1), 1 fraction (patient 2) and 6 fractions (patient 3) where a post-treatment control CBCT was available. After treatment, segmented 2D marker positions in the post-treatment CBCT projections provided the estimated 3D motion trajectory during the CBCT by a probability-based method. An external-internal correlation model (ECM) that estimated the tumor motion from the RPM motion was built from the synchronized RPM signal and marker motion in the CBCT. The ECM was then used to estimate intra-treatment tumor motion. Finally, the motion-including CTV dose was estimated using a dose reconstruction method that emulates tumor motion in beam's eye view as lateral spot shifts and in-depth motion as changes in the proton beam energy. The CTV homogeneity index (HI) The CTV homogeneity index (HI) was calculated as D 2 %  -  D 98 % D 50 %   × 100 % .

Results

The tumor position during spot delivery had a root-mean-square error of 1.3 mm in left-right, 2.8 mm in cranio-caudal and 1.7 mm in anterior-posterior directions compared to the planned position. On average, the CTV HI was larger than planned by 3.7%-points (range: 1.0-6.6%-points) for individual fractions and by 0.7%-points (range: 0.3-1.1%-points) for the average dose of 5 or 6 fractions.

Conclusions

A method to estimate internal tumor motion and reconstruct the motion-including fraction dose for PBS proton therapy of HCC was developed and demonstrated successfully clinically.

Management of Motion and Anatomical Variations in Charged Particle Therapy: Past, Present, and Into the Future

The major aim of radiation therapy is to provide curative or palliative treatment to cancerous malignancies while minimizing damage to healthy tissues. Charged particle radiotherapy utilizing carbon ions or protons is uniquely suited for this task due to its ability to achieve highly conformal dose distributions around the tumor volume. For these treatment modalities, uncertainties in the localization of patient anatomy due to inter- and intra-fractional motion present a heightened risk of undesired dose delivery. A diverse range of mitigation strategies have been developed and clinically implemented in various disease sites to monitor and correct for patient motion, but much work remains. This review provides an overview of current clinical practices for inter and intra-fractional motion management in charged particle therapy, including motion control, current imaging and motion tracking modalities, as well as treatment planning and delivery techniques. We also cover progress to date on emerging technologies including particle-based radiography imaging, novel treatment delivery methods such as tumor tracking and FLASH, and artificial intelligence and discuss their potential impact towards improving or increasing the challenge of motion mitigation in charged particle therapy.

Report of AAPM Task Group 290: Respiratory motion management for particle therapy

Dose uncertainty induced by respiratory motion remains a major concern for treating thoracic and abdominal lesions using particle beams. This Task Group report reviews the impact of tumor motion and dosimetric considerations in particle radiotherapy, current motion‐management techniques, and limitations for different particle‐beam delivery modes (i.e., passive scattering, uniform scanning, and pencil‐beam scanning). Furthermore, the report provides guidance and risk analysis for quality assurance of the motion‐management procedures to ensure consistency and accuracy, and discusses future development and emerging motion‐management strategies. This report supplements previously published AAPM report TG76, and considers aspects of motion management that are crucial to the accurate and safe delivery of particle‐beam therapy. To that end, this report produces general recommendations for commissioning and facility‐specific dosimetric characterization, motion assessment, treatment planning, active and passive motion‐management techniques, image guidance and related decision‐making, monitoring throughout therapy, and recommendations for vendors. Key among these recommendations are that: (1) facilities should perform thorough planning studies (using retrospective data) and develop standard operating procedures that address all aspects of therapy for any treatment site involving respiratory motion; (2) a risk‐based methodology should be adopted for quality management and ongoing process improvement.

Future Developments in Charged Particle Therapy: Improving Beam Delivery for Efficiency and Efficacy

The physical and clinical benefits of charged particle therapy (CPT) are well recognized. However, the availability of CPT and complete exploitation of dosimetric advantages are still limited by high facility costs and technological challenges. There are extensive ongoing efforts to improve upon these, which will lead to greater accessibility, superior delivery, and therefore better treatment outcomes. Yet, the issue of cost remains a primary hurdle as utility of CPT is largely driven by the affordability, complexity and performance of current technology. Modern delivery techniques are necessary but limited by extended treatment times. Several of these aspects can be addressed by developments in the beam delivery system (BDS) which determines the overall shaping and timing capabilities enabling high quality treatments. The energy layer switching time (ELST) is a limiting constraint of the BDS and a determinant of the beam delivery time (BDT), along with the accelerator and other factors. This review evaluates the delivery process in detail, presenting the limitations and developments for the BDS and related accelerator technology, toward decreasing the BDT. As extended BDT impacts motion and has dosimetric implications for treatment, we discuss avenues to minimize the ELST and overview the clinical benefits and feasibility of a large energy acceptance BDS. These developments support the possibility of advanced modalities and faster delivery for a greater range of treatment indications which could also further reduce costs. Further work to realize methodologies such as volumetric rescanning, FLASH, arc, multi-ion and online image guided therapies are discussed. In this review we examine how increased treatment efficiency and efficacy could be achieved with improvements in beam delivery and how this could lead to faster and higher quality treatments for the future of CPT.

Target motion management in breast cancer radiation therapy

Background

Over the last two decades, breast cancer remains the main cause of cancer deaths in women. To treat this type of cancer, radiation therapy (RT) has proved to be efficient. RT for breast cancer is, however, challenged by intrafractional motion caused by respiration. The problem is more severe for the left-sided breast cancer due to the proximity to the heart as an organ-at-risk. While particle therapy results in superior dose characteristics than conventional RT, due to the physics of particle interactions in the body, particle therapy is more sensitive to target motion.

Conclusions

This review highlights current and emerging strategies for the management of intrafractional target motion in breast cancer treatment with an emphasis on particle therapy, as a modern RT technique. There are major challenges associated with transferring real-time motion monitoring technologies from photon to particles beams. Surface imaging would be the dominant imaging modality for real-time intrafractional motion monitoring for breast cancer. The magnetic resonance imaging (MRI) guidance and ultra high dose rate (FLASH)-RT seem to be state-of-the-art approaches to deal with 4D RT for breast cancer.

Carbon ion radiotherapy using fiducial markers for prostate cancer in Osaka HIMAK: Treatment planning

PURPOSE

Carbon ion radiotherapy for prostate cancer was performed using two fine needle Gold Anchor (GA) markers for patient position verification in Osaka Heavy Ion Medical Accelerator in Kansai (Osaka HIMAK). The present study examined treatment plans for prostate cases using beam-specific planning target volume (bsPTV) based on the effect of the markers on dose distribution and analysis of target movements.

MATERIALS AND METHODS

Gafchromic EBT3 film was used to measure dose perturbations caused by markers. First, the relationships between the irradiated film density and absolute dose with different linear energy transfer distributions within a spread-out Bragg peak (SOBP) were confirmed. Then, to derive the effect of markers, two types of markers, including GA, were placed at the proximal, center, and distal depths within the same SOBP, and dose distributions behind the markers were measured using the films. The amount of internal motion of prostate was derived from irradiation results and analyzed to determine the margins of the bsPTV.

RESULTS

The linearity of the film densities against absolute doses was constant within the SOBP and the amount of dose perturbations caused by the markers was quantitatively estimated from the film densities. The dose perturbation close behind the markers was smallest (<10% among depths within the SOBP regardless of types of markers) and increased with depth. The effect of two types of GAs on dose distributions was small and could be ignored in the treatment planning. Based on the analysis results of internal motions of prostate, required margins of the bsPTV were found to be 8, 7, and 7 mm in left-right (LR), anterior-posterior (AP), and superior-inferior (SI) directions, respectively.

CONCLUSION

We evaluated the dose reductions caused by markers and determined the margins of the bsPTV, which was applied to the treatment using fiducial markers, using the analysis results of prostate movements.

Dose Calculation Algorithms for External Radiation Therapy: An Overview for Practitioners

Radiation therapy (RT) is a constantly evolving therapeutic technique; improvements are continuously being introduced for both methodological and practical aspects. Among the features that have undergone a huge evolution in recent decades, dose calculation algorithms are still rapidly changing. This process is propelled by the awareness that the agreement between the delivered and calculated doses is of paramount relevance in RT, since it could largely affect clinical outcomes. The aim of this work is to provide an overall picture of the main dose calculation algorithms currently used in RT, summarizing their underlying physical models and mathematical bases, and highlighting their strengths and weaknesses, referring to the most recent studies on algorithm comparisons. This handy guide is meant to provide a clear and concise overview of the topic, which will prove useful in helping clinical medical physicists to perform their responsibilities more effectively and efficiently, increasing patient benefits and improving the overall quality of the management of radiation treatment.

Validation of dose distribution for liver tumors treated with real-time-image gated spot-scanning proton therapy by log data based dose reconstruction

In spot scanning proton therapy (SSPT), the spot position relative to the target may fluctuate through tumor motion even when gating the radiation by utilizing a fiducial marker. We have established a procedure that evaluates the delivered dose distribution by utilizing log data on tumor motion and spot information. The purpose of this study is to show the reliability of the dose distributions for liver tumors treated with real-time-image gated SSPT (RGPT). In the evaluation procedure, the delivered spot information and the marker position are synchronized on the basis of log data on the timing of the spot irradiation and fluoroscopic X-ray irradiation. Then a treatment planning system reconstructs the delivered dose distribution. Dose distributions accumulated for all fractions were reconstructed for eight liver cases. The log data were acquired in all 168 fractions for all eight cases. The evaluation was performed for the values of maximum dose, minimum dose, D99, and D5-D95 for the clinical target volumes (CTVs) and mean liver dose (MLD) scaled by the prescribed dose. These dosimetric parameters were statistically compared between the planned dose distribution and the reconstructed dose distribution. The mean difference of the maximum dose was 1.3% (95% confidence interval [CI]: 0.6%-2.1%). Regarding the minimum dose, the mean difference was 0.1% (95% CI: -0.5%-0.7%). The mean differences of D99, D5-D95 and MLD were below 1%. The reliability of dose distributions for liver tumors treated with RGPT-SSPT was shown by the evaluation of the accumulated dose distributions.

Combining rescanning and gating for a time-efficient treatment of mobile tumors using pencil beam scanning proton therapy

BACKGROUND AND PURPOSE

Respiratory motion during proton therapy can severely degrade dose distributions, particularly due to interplay effects when using pencil beam scanning. Combined rescanning and gating treatments for moving tumors mitigates dose degradation, but at the cost of increased treatment delivery time. The objective of this study was to identify the time efficiency of these dose degradation-motion mitigation strategies for different range of motions.

MATERIALS AND METHODS

Seventeen patients with thoracic or abdominal tumors were studied. Tumor motion amplitudes ranged from 2-30 mm. Deliveries using different combinations of rescanning and gating were simulated with a dense dose spot grid (4 × 4 × 2.5 mm³) for all patients and a sparse dose spot grid (8 × 8 × 5 mm³) for six patients with larger tumor movements (>8 mm). The resulting plans were evaluated in terms of CTV coverage and time efficiency.

RESULTS

Based on the studied patient cohort, it has been shown that for amplitudes up to 5 mm, no motion mitigation is required with a dense spot grid. For amplitudes between 5 and 10 mm, volumetric rescanning should be applied while maintaining a 100% duty cycle when using a dense spot grid. Although gating could be envisaged to reduce the target volume for intermediate motion, it has been shown that the dose to normal tissues would only be reduced marginally. Moreover, the treatment time would increase. Finally, for larger motion amplitudes, both volumetric rescanning and respiratory gating should be applied with both spot grids. In addition, it has been shown that a dense spot grid delivers better CTV dose coverage than a sparse dose grid.

CONCLUSION

Volumetric rescanning and/or respiratory gating can be used in order to effectively and efficiently mitigate dose degradation due to tumor movement.

Physics of Particle Beam and Hypofractionated Beam Delivery in NSCLC

The dosimetric advantages of particle therapy lead to significantly reduced integral dose to normal tissues, making it an attractive treatment option for body sites such as the thorax. With reduced normal tissue dose comes the potential for dose escalation, toxicity reduction, or hypofractionation. While proton and heavy ion therapy have been used extensively for NSCLC, there are challenges in planning and delivery compared with X-ray-based radiation therapy. Particularly, range uncertainties compounded by breathing motion have to be considered. This article summarizes the current state of particle therapy for NSCLC with a specific focus on the impact of dosimetric uncertainties in planning and delivery.

Liver-ultrasound-guided lung tumour tracking for scanned proton therapy: a feasibility study

Pencil beam scanned (PBS) proton therapy of lung tumours is hampered by respiratory motion and the motion-induced density changes along the beam path. In this simulation study, we aim to investigate the effectiveness of proton beam tracking for lung tumours both under ideal conditions and in conjunction with a respiratory motion model guided by real-time ultrasound imaging of the liver. Multiple-breathing-cycle 4DMRIs of the thorax and abdominal 2D ultrasound images were acquired simultaneously for five volunteers. Deformation vector fields extracted from the 4DMRI, referred to as ground truth motion, were used to generate 4DCT(MRI) data sets of two lung cancer patients, resulting in 10 data sets with variable motion patterns. Given the 4DCT(MRI) and the corresponding ultrasound images as surrogate data, a patient-specific motion model was built. The model consists of an autoregressive model and Gaussian process regression for the temporal and spatial prediction, respectively. Two-field PBS plans were optimised on the reference CTs, and 4D dose calculations (4DDC) were used to simulate dose delivery for (a) unmitigated motion, (b) ideal 2D and 3D tracking (both beam adaption and 4DDC based on ground truth motion), and (c) realistic 2D and 3D tracking (beam adaption based on motion predictions, 4DDC on ground truth motion). Model-guided tracking retrieved clinically acceptable target dose homogeneity, as seen in a substantial reduction of the D5%-D95% compared to the non-mitigated simulation. Tracking in 2D and 3D resulted in a similar improvement of the dose homogeneity, as did ideal and realistic tracking simulations. In some cases, however, the tracked deliveries resulted in a shift towards higher or lower dose levels, leading to unacceptable target over- or under-coverage. The presented motion modelling framework was shown to be an accurate motion prediction tool for the use in proton beam tracking. Tracking alone, however, may not always effectively mitigate motion effects, making it necessary to combine it with other techniques such as rescanning.

Continuous real time 3D motion reproduction using dynamic MRI and precomputed 4DCT deformation fields

Radiotherapy of mobile tumors requires specific imaging tools and models to reduce the impact of motion on the treatment. Online continuous nonionizing imaging has become possible with the recent development of magnetic resonance imaging devices combined with linear accelerators. This opens the way to new guided treatment methods based on the real-time tracking of anatomical motion. In such devices, 2D fast MR-images are well-suited to capture and predict the real-time motion of the tumor. To be used effectively in an adaptive radiotherapy, these MR images have to be combined with X-ray images such as CT, which are necessary to compute the irradiation dose deposition. We therefore developed a method combining both image modalities to track the motion on MR images and reproduce the tracked motion on a sequence of 3DCT images in real-time. It uses manually placed navigators to track organ interfaces in the image, making it possible to select anatomical object borders that are visible on both MRI and CT modalities and giving the operator precise control of the motion tracking quality. Precomputed deformation fields extracted from the 4DCT acquired in the planning phase are then used to deform existing 3DCT images to match the tracked object position, creating a new set of 3DCT images encompassing irregularities in the breathing pattern for the complete duration of the MRI acquisition. The final continuous reconstructed 4DCT image sequence reproduces the motion captured by the MRI sequence with high precision (difference below 2 mm).

Dynamic gating window technique for the reduction of dosimetric error in respiratory-gated spot-scanning particle therapy: An initial phantom study using patient tumor trajectory data

Spot-scanning particle therapy possesses advantages, such as high conformity to the target and efficient energy utilization compared with those of the passive scattering irradiation technique. However, this irradiation technique is sensitive to target motion. In the current clinical situation, some motion management techniques, such as respiratory-gated irradiation, which uses an external or internal surrogate, have been clinically applied. In surrogate-based gating, the size of the gating window is fixed during the treatment in the current treatment system. In this study, we propose a dynamic gating window technique, which optimizes the size of gating window for each spot by considering a possible dosimetric error. The effectiveness of the dynamic gating window technique was evaluated by simulating irradiation using a moving target in a water phantom. In dosimetric characteristics comparison, the dynamic gating window technique exhibited better performance in all evaluation volumes with different effective depths compared with that of the fixed gate approach. The variation of dosimetric characteristics according to the target depth was small in dynamic gate compared to fixed gate. These results suggest that the dynamic gating window technique can maintain an acceptable dose distribution regardless of the target depth. The overall gating efficiency of the dynamic gate was approximately equal or greater than that of the fixed gating window. In dynamic gate, as the target depth becomes shallower, the gating efficiency will be reduced, although dosimetric characteristics will be maintained regardless of the target depth. The results of this study suggest that the proposed gating technique may potentially improve the dose distribution. However, additional evaluations should be undertaken in the future to determine clinical applicability by assuming the specifications of the treatment system and clinical situation.

Modeling RBE-weighted dose variations in irregularly moving abdominal targets treated with carbon ion beams

PURPOSE

To model four-dimensional (4D) relative biological effectiveness (RBE)-weighted dose variations in abdominal lesions treated with scanned carbon ion beam in case of irregular breathing motion.

METHODS

The proposed method, referred to as bioWED method, combines the simulation of tumor motion in a patient- and beam-specific water equivalent depth (WED)-space with RBE modeling, aiming at the estimation of RBE-weighted dose changes due to respiratory motion. The method was validated on a phantom, simulating gated and free breathing dose delivery, and on a patient case, for which free breathing irradiation was assumed and both amplitude and baseline breathing irregularities were simulated through a respiratory motion model. We quantified (a) the effect of motion on the equivalent uniform dose (EUD) and the RBE-weighted dose-volume histograms (DVH), by comparing the planned dose distribution with "ground truth" 4D RBE-weighted doses computed using 4D computed tomography data, and (ii) the estimation error, by comparing the doses estimated with the bioWED method to "ground truth" 4D RBE-weighted doses.

RESULTS

In the phantom validation, the estimation error on the EUD was limited with respect to the motion effect and the median estimation error on relevant RBE-weighted DVH metrics remained within 5%. In the patient study, the estimation error as computed on the EUD was smaller than the corresponding motion effect, exhibiting the largest values in the baseline irregularity simulation. However, the median estimation error over all simulations was below 3.2% considering relevant DVH metrics.

CONCLUSIONS

In the evaluated cases, the bioWED method showed proper accuracy when compared to deformable image registration-based 4D dose calculation. Therefore, it can be seen as a tool to test treatment plan robustness against irregular breathing motion, although its accuracy decreases as a function of increasing soft tissue deformation and should be evaluated on a larger patient dataset.

Managing treatment-related uncertainties in proton beam radiotherapy for gastrointestinal cancers

In recent years, there has been rapid adaption of proton beam radiotherapy (RT) for treatment of various malignancies in the gastrointestinal (GI) tract, with increasing number of institutions implementing intensity modulated proton therapy (IMPT). We review the progress and existing literature regarding the technical aspects of RT planning for IMPT, and the existing tools that can help with the management of uncertainties which may impact the daily delivery of proton therapy. We provide an in-depth discussion regarding range uncertainties, dose calculations, image guidance requirements, organ and body cavity filling consideration, implanted devices and hardware, use of fiducials, breathing motion evaluations and both active and passive motion management methods, interplay effect, general IMPT treatment planning considerations including robustness plan evaluation and optimization, and finally plan monitoring and adaptation. These advances have improved confidence in delivery of IMPT for patients with GI malignancies under various scenarios.

Current status of proton therapy techniques for lung cancer

Proton beams have been used for cancer treatment for more than 28 years, and several technological advancements have been made to achieve improved clinical outcomes by delivering more accurate and conformal doses to the target cancer cells while minimizing the dose to normal tissues. The state-of-the-art intensity modulated proton therapy is now prevailing as a major treatment technique in proton facilities worldwide, but still faces many challenges in being applied to the lung. Thus, in this article, the current status of proton therapy technique is reviewed and issues regarding the relevant uncertainty in proton therapy in the lung are summarized.

RIDOS: A new system for online computation of the delivered dose distributions in scanning ion beam therapy

PURPOSE

To describe a new system for scanned ion beam therapy, named RIDOS (Real-time Ion DOse planning and delivery System), which performs real time delivered dose verification integrating the information from a clinical beam monitoring system with a Graphic Processing Unit (GPU) based dose calculation in patient Computed Tomography.

METHODS

A benchmarked dose computation algorithm for scanned ion beams has been parallelized and adapted to run on a GPU architecture. A workstation equipped with a NVIDIA GPU has been interfaced through a National Instruments PXI-crate with the dose delivery system of the Italian National Center of Oncological Hadrontherapy (CNAO) to receive in real-time the measured beam parameters. Data from a patient monitoring system are also collected to associate the respiratory phases with each spot during the delivery of the dose. Using both measured and planned spot properties, RIDOS evaluates during the few seconds of inter-spill time the cumulative delivered and prescribed dose distributions and compares them through a fast γ-index algorithm.

RESULTS

The accuracy of the GPU-based algorithms was assessed against the CPU-based ones and the differences were found below 1‰. The cumulative planned and delivered doses are computed at the end of each spill in about 300 ms, while the dose comparison takes approximatively 400 ms. The whole operation provides the results before the next spill starts.

CONCLUSIONS

RIDOS system is able to provide a fast computation of the delivered dose in the inter-spill time of the CNAO facility and allows to monitor online the dose deposition accuracy all along the treatment.

Clinical implementation of respiratory-gated spot-scanning proton therapy: An efficiency analysis of active motion management

Purpose

The aim of this work is to describe the clinical implementation of respiratory‐gated spot‐scanning proton therapy (SSPT) for the treatment of thoracic and abdominal moving targets. The experience of our institution is summarized, from initial acceptance and commissioning tests to the development of standard clinical operating procedures for simulation, motion assessment, motion mitigation, treatment planning, and gated SSPT treatment delivery.

Materials and methods

A custom respiratory gating interface incorporating the Real‐Time Position Management System (RPM, Varian Medical Systems, Inc., Palo Alto, CA, USA) was developed in‐house for our synchrotron‐based delivery system. To assess gating performance, a motion phantom and radiochromic films were used to compare gated vs nongated delivery. Site‐specific treatment planning protocols and conservative motion cutoffs were developed, allowing for free‐breathing (FB), breath‐holding (BH), or phase‐gating (Ph‐G). Room usage efficiency of BH and Ph‐G treatments was retrospectively evaluated using beam delivery data retrieved from our record and verify system and DICOM files from patient‐specific quality assurance (QA) procedures.

Results

More than 70 patients were treated using active motion management between the launch of our motion mitigation program in October 2015 and the end date of data collection of this study in January 2018. During acceptance procedures, we found that overall system latency is clinically‐suitable for Ph‐G. Regarding room usage efficiency, the average number of energy layers delivered per minute was <10 for Ph‐G, 10‐15 for BH and ≥15 for FB, making Ph‐G the slowest treatment modality. When comparing to continuous delivery measured during pretreatment QA procedures, the median values of BH treatment time were extended from 6.6 to 9.3 min (+48%). Ph‐G treatments were extended from 7.3 to 13.0 min (+82%).

Conclusions

Active motion management has been crucial to the overall success of our SSPT program. Nevertheless, our conservative approach has come with an efficiency cost that is more noticeable in Ph‐G treatments and should be considered in decision‐making.

A quality assurance for respiratory gated proton irradiation with range modulation wheel

The purpose of this study was to provide periodic quality assurance (QA) methods for respiratory-gated proton beam with a range modulation wheel (RMW) and to clarify the characteristics and long-term stability of the respiratory-gated proton beam. A two-dimensional detector array and a solid water phantom were used to measure absolute dose, spread-out Bragg peak (SOBP) width and proton range for monthly QA. SOBP width and proton range were measured using an oblique incidence beam to the lateral side of a solid water phantom and compared between with and without a gating proton beam. To measure the delay time of beam-on/off for annual QA, we collected the beam-on/off signals and the dose monitor-detected pulse. We analyzed the results of monthly QA over a 15-month period and investigated the delay time by machine signal analysis. The dose deviations at proximal, SOBP center and distal points were -0.083 ± 0.25%, 0.026 ± 0.20%, and -0.083 ± 0.35%, respectively. The maximum dose deviation between with and without respiratory gating was -0.95% at the distal point and other deviations were within ±0.5%. Proximal and SOBP center doses showed the same trend over a 15-month period. Delay times of beam-on/off for 200 MeV/SOBP 16 cm were 140.5 ± 0.8 ms and 22.3 ± 13.0 ms, respectively. Delay times for 160 MeV/SOBP 10 cm were 167.5 ± 15.1 ms and 19.1 ± 9.8 ms. Our beam delivery system with the RMW showed sufficient stability for respiratory-gated proton therapy and the system did not show dependency on the energy and the respiratory wave form. The delay times of beam-on/off were within expectations. The proposed QA methods will be useful for managing the quality of respiratory-gated proton beams and other beam delivery systems.

Commissioning of a respiratory gating system involving a pressure sensor in carbon-ion scanning radiotherapy

This study reports the commissioning methodology and results of a respiratory gating system [AZ - 733 V/733 VI (Anzai Medical Co., Japan)] using a pressure sensor in carbon-ion scanning radiotherapy. Commissioning includes choosing a location and method for pressure sensor installation, delay time measurement of the system, and the final flow test. Additionally, we proposed a methodology for the determination of a threshold level of generating an on/off gate for the beam to the respiratory waveform, which is important for clinical application. Regarding the location and method for installation of the pressure sensor, the actual person's abdomen, back of the body position, and supine/prone positioning were checked. By comparing the motion between the pressure sensor output and the reference LED sensor motion, the chest rear surface was shown to be unsuitable for the sensor installation, due to noise in the signal caused by the cardiac beat. Regarding delay time measurement of the system, measurements were performed for the following four steps: (a). Actual motion to wave signal generation; (b). Wave signal to gate signal generation; (c). Gate signal to beam on/off signal generation; (d). Beam on/off signal to the beam irradiation. The total delay time measured was 46 ms (beam on)/33 ms (beam off); these were within the prescribed tolerance time (<100 ms). Regarding the final flow test, an end-to-end test was performed with a patient verification system using an actual carbon-ion beam; the respiratory gating irradiation was successfully performed, in accordance with the intended timing. Finally, regarding the method for determining the threshold level of the gate generation of the respiration waveform, the target motion obtained from 4D-CT was assumed to be correlated with the waveform obtained from the pressure sensor; it was used to determine the threshold value in amplitude direction.

A novel and individualized robust optimization method using normalized dose interval volume constraints (NDIVC) for intensity-modulated proton radiotherapy

PURPOSE

Intensity-modulated proton therapy (IMPT) is known to be sensitive to patient setup and range uncertainty issues. Multiple robust optimization methods have been developed to mitigate the impact of these uncertainties. Here, we propose a new robust optimization method, which provides an alternative way of robust optimization in IMPT, and is clinically practical, which will enable users to control the balance between nominal plan quality and plan robustness in a user-defined fashion.

METHOD

We calculated nine individual dose distributions which corresponded to one nominal and eight extreme scenarios caused by patient setup and proton beam's range uncertainties. For each voxel, the normalized dose interval (NDI) is defined as the full dose range variation divided by the maximum dose in all uncertainty scenarios (NDI = [max - min dose]/max dose), which was then used to calculate the normalized dose interval volume histogram (NDIVH) curves. The areas under the NDIVH curves were used to quantify plan robustness. A normalized dose interval volume constraint (NDIVC) applied to the target was incorporated to specify the desired robustness which was user-defined. Users could then explore the trade-off between nominal plan quality and plan robustness by adjusting the position of the NDIVCs on the NDIVH curves freely. We benchmarked our method using one lung, five head and neck (H&N), and three prostate cases by comparing our results to those derived using the voxel-wise worst-case robust optimization.

RESULTS

Using the benchmark cases, our new method achieved quality IMPT plans comparable to those derived from the voxel-wise worst-case robust optimization for both nominal plan quality and plan robustness in general; even more conformal and more homogeneous target dose distributions in some cases, if proper NDIVCs were applied. The AUC under NDIVH, as a precise quantitative index of plan robustness, was consistent with DVH bandwidths. Additionally, we demonstrated the feasibility of adjusting the position of NDIVCs in the NDIVH curves which allowed users to explore the trade-off between nominal plan quality and plan robustness.

CONCLUSIONS

The NDIVH-based robust optimization method provided a novel and individualized way of robust optimization in IMPT, and enables users to adjust the balance between nominal plan quality and plan robustness in a user-defined fashion. This method is applicable for continued improvement and developing the next generation of IMPT planning algorithms in the future.

A Monte-Carlo-based and GPU-accelerated 4D-dose calculator for a pencil beam scanning proton therapy system

PURPOSE

The presence of respiratory motion during radiation treatment leads to degradation of the expected dose distribution, both for target coverage and healthy tissue sparing, particularly for techniques like pencil beam scanning proton therapy which have dynamic delivery systems. While tools exist to estimate this degraded four-dimensional (4D) dose, they typically have one or more deficiencies such as not including the particular effects from a dynamic delivery, using analytical dose calculations, and/or using nonphysical dose-accumulation methods. This work presents a clinically useful 4D-dose calculator that addresses each of these shortcomings.

METHODS

To quickly compute the 4D dose, the three main tasks of the calculator were run on graphics processing units (GPUs). These tasks were (a) simulating the delivery of the plan using measured delivery parameters to distribute the plan amongst 4DCT phases characterizing the patient breathing, (b) using an in-house Monte Carlo simulation (MC) dose calculator to determine the dose delivered to each breathing phase, and (c) accumulating the doses from the various breathing phases onto a single phase for evaluation. The accumulation was performed by individually transferring the energy and mass of dose-grid subvoxels, a technique that models the transfer of dose in a more physically realistic manner. The calculator was run on three test cases, with lung, esophagus, and liver targets, respectively, to assess the various uncertainties in the beam delivery simulation as well as to characterize the dose-accumulation technique.

RESULTS

Four-dimensional doses were successfully computed for the three test cases with computation times ranging from 4-6 min on a server with eight NVIDIA Titan X graphics cards; the most time-consuming component was the MC dose engine. The subvoxel-based dose-accumulation technique produced stable 4D-dose distributions at subvoxel scales of 0.5-1.0 mm without impairing the total computation time. The uncertainties in the beam delivery simulation led to moderate variations of the dose-volume histograms for these cases; the variations were reduced by implementing repainting or phase-gating motion mitigation techniques in the calculator.

CONCLUSIONS

A MC-based and GPU-accelerated 4D-dose calculator was developed to estimate the effects of respiratory motion on pencil beam scanning proton therapy treatments. After future validation, the calculator could be used to assess treatment plans and its quick runtime would make it easily usable in a future 4D-robust optimization system.

Effectiveness of different rescanning techniques for scanned proton radiotherapy in lung cancer patients

Non-small cell lung cancer (NSCLC) is a tumour type thought to be well-suited for proton radiotherapy. However, the lung region poses many problems related to organ motion and can for actively scanned beams induce severe interplay effects. In this study we investigate four mitigating rescanning techniques: (1) volumetric rescanning, (2) layered rescanning, (3) breath-sampled (BS) layered rescanning, and (4) continuous breath-sampled (CBS) layered rescanning. The breath-sampled methods will spread the layer rescans over a full breathing cycle, resulting in an improved averaging effect at the expense of longer treatment times. In CBS, we aim at further improving the averaging by delivering as many rescans as possible within one breathing cycle. The interplay effect was evaluated for 4D robustly optimized treatment plans (with and without rescanning) for seven NSCLC patients in the treatment planning system RayStation. The optimization and final dose calculation used a Monte Carlo dose engine to account for the density heterogeneities in the lung region. A realistic treatment delivery time structure given from the IBA ScanAlgo simulation tool served as basis for the interplay evaluation. Both slow (2.0 s) and fast (0.1 s) energy switching times were simulated. For all seven studied patients, rescanning improves the dose conformity to the target. The general trend is that the breath-sampled techniques are superior to layered and volumetric rescanning with respect to both target coverage and variability in dose to OARs. The spacing between rescans in our breath-sampled techniques is set at planning, based on the average breathing cycle length obtained in conjunction with CT acquisition. For moderately varied breathing cycle lengths between planning and delivery (up to 15%), the breath-sampled techniques still mitigate the interplay effect well. This shows the potential for smooth implementation at the clinic without additional motion monitoring equipment.

Impact of Spot Size and Spacing on the Quality of Robustly Optimized Intensity Modulated Proton Therapy Plans for Lung Cancer

PURPOSE

To investigate how spot size and spacing affect plan quality, robustness, and interplay effects of robustly optimized intensity modulated proton therapy (IMPT) for lung cancer.

METHODS AND MATERIALS

Two robustly optimized IMPT plans were created for 10 lung cancer patients: first by a large-spot machine with in-air energy-dependent large spot size at isocenter (σ: 6-15 mm) and spacing (1.3 σ), and second by a small-spot machine with in-air energy-dependent small spot size (σ: 2-6 mm) and spacing (5 mm). Both plans were generated by optimizing radiation dose to internal target volume on averaged 4-dimensional computed tomography scans using an in-house-developed IMPT planning system. The dose-volume histograms band method was used to evaluate plan robustness. Dose evaluation software was developed to model time-dependent spot delivery to incorporate interplay effects with randomized starting phases for each field per fraction. Patient anatomy voxels were mapped phase-to-phase via deformable image registration, and doses were scored using in-house-developed software. Dose-volume histogram indices, including internal target volume dose coverage, homogeneity, and organs at risk (OARs) sparing, were compared using the Wilcoxon signed-rank test.

RESULTS

Compared with the large-spot machine, the small-spot machine resulted in significantly lower heart and esophagus mean doses, with comparable target dose coverage, homogeneity, and protection of other OARs. Plan robustness was comparable for targets and most OARs. With interplay effects considered, significantly lower heart and esophagus mean doses with comparable target dose coverage and homogeneity were observed using smaller spots.

CONCLUSIONS

Robust optimization with a small spot-machine significantly improves heart and esophagus sparing, with comparable plan robustness and interplay effects compared with robust optimization with a large-spot machine. A small-spot machine uses a larger number of spots to cover the same tumors compared with a large-spot machine, which gives the planning system more freedom to compensate for the higher sensitivity to uncertainties and interplay effects for lung cancer treatments.

Robust optimization in IMPT using quadratic objective functions to account for the minimum MU constraint

PURPOSE

Currently, in clinical practice of intensity-modulated proton therapy (IMPT), the influence of the minimum monitor unit (MU) constraint is taken into account through postprocessing after the optimization is completed. This may degrade the plan quality and plan robustness. This study aims to mitigate the impact of the minimum MU constraint directly during the plan robust optimization.

METHODS AND MATERIALS

Cao et al. have demonstrated a two-stage method to account for the minimum MU constraint using linear programming without the impact of uncertainties considered. In this study, we took the minimum MU constraint into consideration using quadratic optimization and simultaneously had the impact of uncertainties considered using robust optimization. We evaluated our method using seven cancer patients with different machine settings.

RESULT

The new method achieved better plan quality than the conventional method. The D_95% of the clinical target volume (CTV) normalized to the prescription dose was (mean [min-max]): (99.4% [99.2%-99.6%]) vs. (99.2% [98.6%-99.6%]). Plan robustness derived from these two methods was comparable. For all seven patients, the CTV dose-volume histogram band gap (narrower band gap means more robust plans) at D_95% normalized to the prescription dose was (mean [min-max]): (1.5% [0.5%-4.3%]) vs. (1.2% [0.6%-3.8%]).

CONCLUSION

Our new method of incorporating the minimum MU constraint directly into the plan robust optimization can produce machine-deliverable plans with better tumor coverage while maintaining high-plan robustness.

Advanced Radiation DOSimetry phantom (ARDOS): a versatile breathing phantom for 4D radiation therapy and medical imaging

A novel breathing phantom was designed for being used in conventional and ion-beam radiotherapy as well as for medical imaging. Accurate dose delivery and patient safety are aimed to be verified for four-dimensional (4D) treatment techniques compensating for breathing-induced tumor motion. The phantom includes anthropomorphic components representing an average human thorax. It consists of real tissue equivalent materials to fulfill the requirements for dosimetric experiments and imaging purposes. The different parts of the torso (lungs, chest wall, and ribs) and the tumor can move independently. Simple regular movements, as well as more advanced patient-specific breathing cycles are feasible while a reproducible setup can be guaranteed. The phantom provides the flexibility to use different types of dosimetric devices and was designed in a way that it is robust, transportable and easy to handle. Tolerance levels and the reliability of the phantom setup were determined in combination with tests on motion accuracy and reproducibility by using infrared optical tracking technology. Different imaging was performed including positron emission tomography imaging, 4D computed tomography as well as real-time in-room imaging. The initial dosimetric benchmarking studies were performed in a photon beam where dose parameters are predictable and the dosimetric procedures well established.

Empowering Intensity Modulated Proton Therapy Through Physics and Technology: An Overview

Considering the clinical potential of protons attributable to their physical characteristics, interest in proton therapy has increased greatly in this century, as has the number of proton therapy installations. Until recently, passively scattered proton therapy was used almost entirely. Notably, the overall clinical results to date have not shown a convincing benefit of protons over photons. A rapid transition is now occurring with the implementation of the most advanced form of proton therapy, intensity modulated proton therapy (IMPT). IMPT is superior to passively scattered proton therapy and intensity modulated radiation therapy (IMRT) dosimetrically. However, numerous limitations exist in the present IMPT methods. In particular, compared with IMRT, IMPT is highly vulnerable to various uncertainties. In this overview we identify three major areas of current limitations of IMPT: treatment planning, treatment delivery, and motion management, and discuss current and future efforts for improvement. For treatment planning, we need to reduce uncertainties in proton range and in computed dose distributions, improve robust planning and optimization, enhance adaptive treatment planning and delivery, and consider how to exploit the variability in the relative biological effectiveness of protons for clinical benefit. The quality of proton therapy also depends on the characteristics of the IMPT delivery systems and image guidance. Efforts are needed to optimize the beamlet spot size for both improved dose conformality and faster delivery. For the latter, faster energy switching time and increased dose rate are also needed. Real-time in-room volumetric imaging for guiding IMPT is in its early stages with cone beam computed tomography (CT) and CT-on-rails, and continued improvements are anticipated. In addition, imaging of the proton beams themselves, using, for instance, prompt γ emissions, is being developed to determine the proton range and to reduce range uncertainty. With the realization of the advances described above, we posit that IMPT, thus empowered, will lead to substantially improved clinical results.

Comparison of photon volumetric modulated arc therapy, intensity-modulated proton therapy, and intensity-modulated carbon ion therapy for delivery of hypo-fractionated thoracic radiotherapy

Purpose

The aim of the present study was to compare the dose distribution generated from photon volumetric modulated arc therapy (VMAT), intensity modulated proton therapy (IMPT), and intensity modulated carbon ion therapy (IMCIT) in the delivery of hypo-fractionated thoracic radiotherapy.

Methods and materials

Ten selected patients who underwent thoracic particle therapy between 2015 and 2016 were re-planned to receive a relative biological effectiveness (RBE) weighted dose of 60 Gy (i.e., GyE) in 15 fractions delivered with VMAT, IMPT, or IMCIT with the same optimization criteria. Treatment plans were then compared.

Results

There were no significant differences in target volume dose coverage or dose conformity, except improved D₉₅ was found with IMCIT compared with VMAT (p = 0.01), and IMCIT was significantly better than IMPT in all target volume dose parameters. Particle therapy led to more prominent lung sparing at low doses, and this result was most prominent with IMCIT (p < 0.05). Improved sparing of other thoracic organs at risk (OARs) was observed with particle therapy, and IMCIT further lowered the D_1cc and D_5cc for major blood vessels, as compared with IMPT (p = 0.01).

Conclusion

Although it was comparable to VMAT, IMCIT led to significantly better tumor target dose coverage and conformity than did IMPT. Particle therapy, compared with VMAT, improved thoracic OAR sparing. IMCIT, compared with IMPT, may further improve normal lung and major blood vessel sparing under limited respiratory motion.

Electronic supplementary material

The online version of this article (doi:10.1186/s13014-017-0866-0) contains supplementary material, which is available to authorized users.

A study of the beam-specific interplay effect in proton pencil beam scanning delivery in lung cancer

BACKGROUND

For lung tumors with large motion amplitudes, the use of proton pencil beam scanning (PBS) can produce large dose errors. In this study, we assess under what circumstances PBS can be used to treat lung cancer patients who exhibit large tumor motion, based on the quantification of tumor motion and the dose interplay.

MATERIAL AND METHODS

PBS plans were optimized on average 4DCT datasets using a beam-specific PTV method for 10 consecutive patients with locally advanced non-small-cell-lung-cancer (NSCLC) treated with proton therapy to 6660/180 cGy. End inhalation (CT0) and end exhalation (CT50) were selected as the two extreme scenarios to acquire the relative stopping power ratio difference (Δrsp) for a respiration cycle. The water equivalent difference (ΔWET) per radiological path was calculated from the surface of patient to the iCTV by integrating the Δrsp of each voxel. The magnitude of motion of voxels within the target follows a quasi-Gaussian distribution. A motion index (MI (>5mm WET)), defined as the percentage of target voxels with an absolute integral ΔWET larger than 5 mm, was adopted as a metric to characterize interplay. To simulate the treatment process, 4D dose was calculated by accumulating the spot dose on the corresponding respiration phase to the reference phase CT50 by deformable image registration based on spot timing and patient breathing phase.

RESULTS

The study indicated that the magnitude of target underdose in a single fraction plan is proportional to the MI (p < .001), with larger motion equating to greater dose degradation and standard deviations. The target homogeneity, minimum, maximum and mean dose in the 4D dose accumulations of 37 fractions varied as a function of MI.

CONCLUSIONS

This study demonstrated that MI can predict the level of dose degradation, which potentially serves as a clinical decision tool to assess whether lung cancer patients are potentially suitable to receive PBS treatment.

Evaluation of Dosimetric Robustness of Carbon Ion Boost Therapy for Anal Carcinoma

Purpose

The radiation therapy treatment outcome of human papillomavirus-negative anal carcinoma may be improved by the biological effectiveness of carbon ions. However, abdominal tissue motion can compromise the precision of carbon ion therapy. This work aims to evaluate the dosimetric feasibility of carbon ion boost (CIB) therapy for anal carcinoma.

Materials and Methods

An algorithm to generate computed tomographies based on daily magnetic resonance imaging data and deformable image registration was developed. By means of this algorithm, fractional computed tomography data for 54 treatment fractions for 3 different patients with anal carcinoma were derived. The dose for a sequential CIB (CIB_seq) treatment plan was recalculated on the fractional computed tomography data and accumulated over the number of fractions. The resulting dose distributions were compared to standard intensity-modulated radiation therapy treatment with an integrated photon boost.

Results

For the investigated patient cases, similar dosimetric results for CIB_seq treatment and for intensity-modulated radiation therapy with an integrated photon boost were found. For CIB_seq treatment, bladder-filling variation had the strongest influence on the dose distribution. However, the detrimental effects on the mean target dose remained below 1 Gy (RBE) as compared to photon therapy.

Conclusion

This study shows the dosimetric feasibility of CIB therapy for anal carcinoma for the first time and gives reason for clinical exploitation of the enhanced biological effect of carbon ions for patients with human papillomavirus-negative anal cancer.

Robustness of 4D-optimized scanned carbon ion beam therapy against interfractional changes in lung cancer

BACKGROUND AND PURPOSE

Moving targets could be conformally treated with actively scanned carbon ion beams using 4D-optimization. As this heavily exploits 4D-CTs, an important question is whether the conformity also upholds in the context of interfractional changes, i.e. variable positioning, anatomy and breathing patterns.

MATERIALS AND METHODS

In 4 lung cancer patients, 6 weekly 4D-CTs were available. 4D-CTs and their phases were non-rigidly registered to propagate contours and 4D-doses. On the first 4D-CT, a 4D-optimized plan delivering a uniform dose to each motion phase (total dose 9.4Gy(RBE)) was simulated, as well as an ITV plan for comparison. On the five following 4D-CTs, 4D-dose was forward calculated and evaluated for target coverage and conformity. Variable uniform (3-7mm) and range margins (2mm/%) were investigated.

RESULTS

For all patients, target coverage (V95>95% accumulated over 5 fractions) could be achieved, but with variable margin size weakly depending on motion amplitude and range changes. The same margins were also necessary for ITV plans, which lead to lower conformity and higher integral doses.

CONCLUSION

4D-optimization appears feasible also under interfractional changes and maintains a dosimetric advantage over less conformal ITV irradiations. Further studies are needed to identify patients benefiting most from the technically more complex 4D-optimization.

Feasibility Study on Cardiac Arrhythmia Ablation Using High-Energy Heavy Ion Beams

High-energy ion beams are successfully used in cancer therapy and precisely deliver high doses of ionizing radiation to small deep-seated target volumes. A similar noninvasive treatment modality for cardiac arrhythmias was tested here. This study used high-energy carbon ions for ablation of cardiac tissue in pigs. Doses of 25, 40, and 55 Gy were applied in forced-breath-hold to the atrioventricular junction, left atrial pulmonary vein junction, and freewall left ventricle of intact animals. Procedural success was tracked by (1.) in-beam positron-emission tomography (PET) imaging; (2.) intracardiac voltage mapping with visible lesion on ultrasound; (3.) lesion outcomes in pathohistolgy. High doses (40–55 Gy) caused slowing and interruption of cardiac impulse propagation. Target fibrosis was the main mediator of the ablation effect. In irradiated tissue, apoptosis was present after 3, but not 6 months. Our study shows feasibility to use high-energy ion beams for creation of cardiac lesions that chronically interrupt cardiac conduction.

Respiration-gated fast-rescanning carbon-ion radiotherapy

Phase-controlled rescanning of the carbon-ion beam offers fast and precise dose application with decreased irradiation of normal tissue. However, organ movement with respiration remains a unique challenge. Technological development has enabled the simultaneous application of beam-energy-modulated markerless phase-controlled rescanning with respiration gating, allowing scanning treatment of respiration-mobile tumors with carbon. A total of 10 patients with tumors in the liver or lung were treated in a feasibility study at our facility using this combination. At a median of 10.5 months, follow-up examination including computed tomography/magnetic resonance imaging revealed no grade 2+ acute adverse effects with this new therapy. Two patients with complex disease experienced local recurrence, which may be improved with increased dose delivery. One patient died of unrelated causes. All other patients are alive with good control at the time of writing. Though long-term observation is pending, these are promising initial results for use of the carbon-beam phase-controlled rescanning method in respiration-mobile disease.

Optimization of an on-board imaging system for extremely rapid radiation therapy

PURPOSE

Next-generation extremely rapid radiation therapy systems could mitigate the need for motion management, improve patient comfort during the treatment, and increase patient throughput for cost effectiveness. Such systems require an on-board imaging system that is competitively priced, fast, and of sufficiently high quality to allow good registration between the image taken on the day of treatment and the image taken the day of treatment planning. In this study, three different detectors for a custom on-board CT system were investigated to select the best design for integration with an extremely rapid radiation therapy system.

METHODS

Three different CT detectors are proposed: low-resolution (all 4×4 mm pixels), medium-resolution (a combination of 4×4 mm pixels and 2×2 mm pixels), and high-resolution (all 1×1 mm pixels). An in-house program was used to generate projection images of a numerical anthropomorphic phantom and to reconstruct the projections into CT datasets, henceforth called "realistic" images. Scatter was calculated using a separate Monte Carlo simulation, and the model included an antiscatter grid and bowtie filter. Diagnostic-quality images of the phantom were generated to represent the patient scan at the time of treatment planning. Commercial deformable registration software was used to register the diagnostic-quality scan to images produced by the various on-board detector configurations. The deformation fields were compared against a "gold standard" deformation field generated by registering initial and deformed images of the numerical phantoms that were used to make the diagnostic and treatment-day images. Registrations of on-board imaging system data were judged by the amount their deformation fields differed from the corresponding gold standard deformation fields--the smaller the difference, the better the system. To evaluate the registrations, the pointwise distance between gold standard and realistic registration deformation fields was computed.

RESULTS

By most global metrics (e.g., mean, median, and maximum pointwise distance), the high-resolution detector had the best performance but the medium-resolution detector was comparable. For all medium- and high-resolution detector registrations, mean error between the realistic and gold standard deformation fields was less than 4 mm. By pointwise metrics (e.g., tracking a small lesion), the high- and medium-resolution detectors performed similarly. For these detectors, the smallest error between the realistic and gold standard registrations was 0.6 mm and the largest error was 3.6 mm.

CONCLUSIONS

The medium-resolution CT detector was selected as the best for an extremely rapid radiation therapy system. In essentially all test cases, data from this detector produced a significantly better registration than data from the low-resolution detector and a comparable registration to data from the high-resolution detector. The medium-resolution detector provides an appropriate compromise between registration accuracy and system cost.

Motion mitigation for lung cancer patients treated with active scanning proton therapy

PURPOSE

Motion interplay can affect the tumor dose in scanned proton beam therapy. This study assesses the ability of rescanning and gating to mitigate interplay effects during lung treatments.

METHODS

The treatments of five lung cancer patients [48 Gy(RBE)/4fx] with varying tumor size (21.1-82.3 cm(3)) and motion amplitude (2.9-30.6 mm) were simulated employing 4D Monte Carlo. The authors investigated two spot sizes (σ ∼ 12 and ∼ 3 mm), three rescanning techniques (layered, volumetric, breath-sampled volumetric) and respiratory gating with a 30% duty cycle.

RESULTS

For 4/5 patients, layered rescanning 6/2 times (for the small/large spot size) maintains equivalent uniform dose within the target >98% for a single fraction. Breath sampling the timing of rescanning is ∼ 2 times more effective than the same number of continuous rescans. Volumetric rescanning is sensitive to synchronization effects, which was observed in 3/5 patients, though not for layered rescanning. For the large spot size, rescanning compared favorably with gating in terms of time requirements, i.e., 2x-rescanning is on average a factor ∼ 2.6 faster than gating for this scenario. For the small spot size however, 6x-rescanning takes on average 65% longer compared to gating. Rescanning has no effect on normal lung V20 and mean lung dose (MLD), though it reduces the maximum lung dose by on average 6.9 ± 2.4/16.7 ± 12.2 Gy(RBE) for the large and small spot sizes, respectively. Gating leads to a similar reduction in maximum dose and additionally reduces V20 and MLD. Breath-sampled rescanning is most successful in reducing the maximum dose to the normal lung.

CONCLUSIONS

Both rescanning (2-6 times, depending on the beam size) as well as gating was able to mitigate interplay effects in the target for 4/5 patients studied. Layered rescanning is superior to volumetric rescanning, as the latter suffers from synchronization effects in 3/5 patients studied. Gating minimizes the irradiated volume of normal lung more efficiently, while breath-sampled rescanning is superior in reducing maximum doses to organs at risk.

Estimation of organ doses and effective doses in image-guided respiration-gated radiotherapy

Dose conformity in thoracic and abdominal ion-beam radiotherapy is degraded by respiratory motion. To improve conformity, an image-guided respiration-gated system can be used in the treatment room. The purpose of this study was to estimate the organ doses and effective doses to patients from an image-guided respiration-gated system. Glass dosemeters were inserted into an adult anthropomorphic phantom and were attached to the surface on the phantom. The phantom was placed on the treatment couch, and the imaging dose from fluoroscopy was evaluated. In addition to the organ doses, the effective doses were also estimated according to the ICRP Publication 103. The irradiation time is over 3-5 min per beam angle. When image acquisition conditions were assumed for thoracic treatment, the effective doses and maximal skin doses were 0.48-0.79 mSv and 5.9-9.9 mGy, respectively. The estimated doses can be the base data for considering radiological protection in the radiotherapy.