Evaluation of residual abdominal tumour motion in carbon ion gated treatments through respiratory motion modelling

Real-time nonstandard-shaped gold fiducial marker tracking on x-ray fluoroscopic images for prostate radiotherapy

Objective.The prostate moves in accordance with the movement of surrounding organs. Tumor position can change by ≥3 mm during radiotherapy. Given the difficulties of visualizing the prostate fluoroscopically, fiducial markers are generally implanted into the prostate to monitor its motion during treatment. Recently, internally motion guidance methods of the prostate using a 99.5% gold/0.5% iron flexible notched wire fiducial marker (Gold Anchor® , Naslund Medical AB, Huddinge, Sweden), which requires a 22 gauge needle, has been used. However, because the notched wire can retain its linear shape, acquire a spiral shape, or roll into an irregular ball, detecting it on fluoroscopic images in real-time incurs higher computation costs.Approach.We developed a fiducial tracking algorithm to achieve real-time computation. The marker is detected on the first image frame using a shape filter that employs inter-class variance for the marker likelihood calculated by the filter, focusing on the large difference in densities between the marker and its surroundings. After the second frame, the marker is tracked by adding to the shape filter the similarity to the template cropped from the area around the marker position detected in the first frame. We retrospectively evaluated the algorithm's marker tracking accuracy for ten prostate cases, analyzing two fractions in each case.Main results.Tracking positional accuracy averaged over all patients was 0.13 ± 0.04 mm (mean ± standard deviation, Euclidean distance) and 0.25 ± 0.09 mm (95th percentile). Computation time was 2.82 ± 0.20 ms/frame averaged over all frames.Significance.Our algorithm successfully and stably tracked irregularly-shaped markers in real time.

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.

The role of multiple anatomical scenarios in plan optimization for carbon ion radiotherapy of pancreatic cancer

PURPOSE /OBJECTIVE

To quantify benefits of robust optimization on multiple 4DCT acquisitions combined with off-line treatment adaptation for neoadjuvant carbon ion therapy (CIRT) of pancreatic cancer.

MATERIAL/METHODS

For 10 previously treated patients, 4DCTs were acquired around -15 (CT_Plan), -5 (RE₁), -1 (RE₂) and +6 (RE₃) days from RT start. Treatment plans were newly optimized to a dose prescription of 38.4 Gy(RBE) (8 fractions) with a constraint of 38 Gy(RBE) to 1% of the gastrointestinal organs at risk volume (D_1%). Three strategies were tested: (A) robust optimization on CT_Plan maximum exhale (0Ex) with 3 mm set-up, 3% range uncertainty, including 30%-inhale; (B) addition of the RE₁-0Ex scenario; (C) plan recalculation at each RE_i and adaptation (RP_i) according to deviation thresholds from clinical goals. The cumulative variation of target coverage and GI-OARs doses was evaluated. Duodenum contours of all 4DCTs of each patient were registered on CT_Plan-0Ex. The capacity of pre-RT acquisitions to predict duodenum position was investigated by computing the intersection of contours at CT_plan, RE₁, or their union, with respect to subsequent 4DCTs and the CTV, coupled with increasing margin.

RESULTS

(A) No recalculation exceeded the D_1% constraint. (B) The inclusion of RE₁-0Ex in the optimization problem improved inter-fraction robustness on a patient-specific basis, but was non-significant on average. (C) Half of the plans would be re-optimized to recover target coverage and/or minimize duodenum dose, at least once. A significant difference was observed between pre-RT duodenum contours when intersecting subsequent contours, either with a margin expansion.

CONCLUSION

Anatomical variations highlighted at multiple RE_i proved that a fast and efficient online adaptation is essential to optimize treatment quality of CIRT for pancreatic cancer.

Time-resolved MRI for off-line treatment robustness evaluation in carbon-ion radiotherapy of pancreatic cancer

PURPOSE

In this study, we investigate the use of magnetic resonance imaging (MRI) for the clinical evaluation of gating treatment robustness in carbon-ion radiotherapy (CIRT) of pancreatic cancer. Indeed, MRI allows radiation-free repeated scans and fast dynamic sequences for time-resolved (TR) imaging (cine-MRI), providing information on inter- and intra-fraction cycle-to-cycle variations of respiratory motion. MRI can therefore support treatment planning and verification, overcoming the limitations of the current clinical standard, that is, four-dimensional computed tomography (4DCT), which describes an "average" breathing cycle neglecting breathing motion variability.

METHODS

We integrated a technique to generate a virtual CT (vCT) from 3D MRI with a method for 3D reconstruction from 2D cine-MRI, to produce TR vCTs for dose recalculations. For eight patients, the method allowed evaluating inter-fraction variations at end-exhale and intra-fraction cycle-to-cycle variability within the gating window in terms of tumor displacement and dose to the target and organs at risk.

RESULTS

The median inter-fraction tumor motion was in the range 3.33-12.16 mm, but the target coverage was robust (-0.4% median D95% variation). Concerning cycle-to-cycle variations, the gating technique was effective in limiting tumor displacement (1.35 mm median gating motion) and corresponding dose variations (-3.9% median D95% variation). The larger exposure of organs at risk (duodenum and stomach) was caused by inter-fraction motion, whereas intra-fraction cycle-to-cycle dose variations were limited.

CONCLUSIONS

This study proposed a method for the generation of TR vCTs from MRI, which enabled an off-line evaluation of gating treatment robustness and suggested its feasibility to support treatment planning of pancreatic tumors in CIRT.

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.

Hadron Therapy Achievements and Challenges: The CNAO Experience

Protons and carbon ions (hadrons) have useful properties for the treatments of patients affected by oncological pathologies. They are more precise than conventional X-rays and possess radiobiological characteristics suited for treating radio-resistant or inoperable tumours. This paper gives an overview of the status of hadron therapy around the world. It focusses on the Italian National Centre for Oncological Hadron therapy (CNAO), introducing operation procedures, system performance, expansion projects, methodologies and modelling to build individualized treatments. There is growing evidence that supports safety and effectiveness of hadron therapy for a variety of clinical situations. However, there is still a lack of high-level evidence directly comparing hadron therapy with modern conventional radiotherapy techniques. The results give an overview of pre-clinical and clinical research studies and of the treatments of 3700 patients performed at CNAO. The success and development of hadron therapy is strongly associated with the creation of networks among hadron therapy facilities, clinics, universities and research institutions. These networks guarantee the growth of cultural knowledge on hadron therapy, favour the efficient recruitment of patients and present available competences for R&D (Research and Development) programmes.

Commissioning and quality assurance of a novel solution for respiratory-gated PBS proton therapy based on optical tracking of surface markers

We present the commissioning and quality assurance of our clinical protocol for respiratory gating in pencil beam scanning proton therapy for cancer patients with moving targets. In a novel approach, optical tracking has been integrated in the therapy workflow and used to monitor respiratory motion from multiple surrogates, applied on the patients' chest. The gating system was tested under a variety of experimental conditions, specific to proton therapy, to evaluate reaction time and reproducibility of dose delivery control. The system proved to be precise in the application of beam gating and allowed the mitigation of dose distortions even for large (1.4cm) motion amplitudes, provided that adequate treatment windows were selected. The total delivered dose was not affected by the use of gating, with measured integral error within 0.15cGy. Analysing high-resolution images of proton transmission, we observed negligible discrepancies in the geometric location of the dose as a function of the treatment window, with gamma pass rate greater than 95% (2%/2mm) compared to stationary conditions. Similarly, pass rate for the latter metric at the 3%/3mm level was observed above 97% for clinical treatment fields, limiting residual movement to 3mm at end-exhale. These results were confirmed in realistic clinical conditions using an anthropomorphic breathing phantom, reporting a similarly high 3%/3mm pass rate, above 98% and 94%, for regular and irregular breathing, respectively. Finally, early results from periodic QA tests of the optical tracker have shown a reliable system, with small variance observed in static and dynamic measurements.

An MRI framework for respiratory motion modelling validation

Introduction

Respiratory motion models establish a correspondence between respiratory‐correlated (RC) 4‐dimensional (4D) imaging and respiratory surrogates, to estimate time‐resolved (TR) 3D breathing motion. To evaluate the performance of motion models on real patient data, a validation framework based on magnetic resonance imaging (MRI) is proposed, entailing the use of RC 4DMRI to build the model, and on both (i) TR 2D cine‐MRI and (ii) additional 4DMRI data for testing intra‐/inter‐fraction breathing motion variability.

Methods

Repeated MRI data were acquired in 7 patients with abdominal lesions. The considered model relied on deformable image registration (DIR) for building the model and compensating for inter‐fraction baseline variations. Both 2D and 3D validation were performed, by comparing model estimations with the ground truth 2D cine‐MRI and 4DMRI respiratory phases, respectively.

Results

The median DIR error was comparable to the voxel size (1.33 × 1.33 × 5 mm³), with higher values in the presence of large inter‐fraction motion (median value: 2.97 mm). In the 2D validation, the median estimation error on anatomical landmarks’ position resulted below 4 mm in every scenario, whereas in the 3D validation it was 1.33 mm and 4.21 mm when testing intra‐ and inter‐fraction motion, respectively. The range of motion described in the cine‐MRI was comparable to the motion of the building 4DMRI, being always above the estimation error. Overall, the model performance was dependent on DIR error, presenting reduced accuracy when inter‐fraction baseline variations occurred.

Conclusions

Results suggest the potential of the proposed framework in evaluating global motion models for organ motion management in MRI‐guided radiotherapy.

A Modular System for Treating Moving Anatomical Targets With Scanned Ion Beams at Multiple Facilities: Pre-Clinical Testing for Quality and Safety of Beam Delivery

Background

Quality management and safety are integral to modern radiotherapy. New radiotherapy technologies require new consensus guidelines on quality and safety. Established analysis strategies, such as the failure modes and effects analysis (FMEA) and incident learning systems have been developed as tools to assess the safety of several types of radiation therapies. An extensive literature documents the widespread application of risk analysis methods to photon radiation therapy. Relatively little attention has been paid to performing risk analyses of nascent radiation therapy systems to treat moving tumors with scanned heavy ion beams. The purpose of this study was to apply a comprehensive safety analysis strategy to a motion-synchronized dose delivery system (M-DDS) for ion therapy.

Methods

We applied a risk analysis method to new treatment planning and treatment delivery processes with scanned heavy ion beams. The processes utilize a prototype, modular dose delivery system, currently undergoing preclinical testing, that provides new capabilities for treating moving anatomy. Each step in the treatment process was listed in a process map, potential errors for each step were identified and scored using the risk probability number in an FMEA, and the possible causes of each error were described in a fault tree analysis. Solutions were identified to mitigate the risk of these errors, including permanent corrective actions, periodic quality assurance (QA) tests, and patient specific QA (PSQA) tests. Each solution was tested experimentally.

Results

The analysis revealed 58 potential errors that could compromise beam delivery quality or safety. Each of the 14 binary (pass-or-fail) tests passed. Each of the nine QA and four PSQA tests were within anticipated clinical specifications. The modular M-DDS was modified accordingly, and was found to function at two centers.

Conclusion

We have applied a comprehensive risk analysis strategy to the M-DDS and shown that it is a clinically viable motion mitigation strategy. The described strategy can be utilized at any ion therapy center that operates with the modular M-DDS. The approach can also be adapted for use at other facilities and can be combined with existing safety analysis systems.

Investigating the use of virtual 4DCT from 4DMRI in gated carbon ion radiation therapy of abdominal tumors

PURPOSE

To generate virtual 4DCT from 4DMRI with field of view (FOV) extended to the entire involved patient anatomy, in order to evaluate its use in carbon ion radiation therapy (CIRT) of the abdominal site in a clinical scenario.

MATERIALS AND METHODS

The virtual 4DCT was generated by deforming a reference CT in order to (1) match the anatomy depicted in the 4DMRI within its FOV, by calculating deformation fields with deformable image registration to describe inter-fractional and breathing motion, and (2) obtain physically plausible deformation outside of the 4DMRI FOV, by propagating and modulating the previously obtained deformation fields. The implemented method was validated on a digital anthropomorphic phantom, for which a ground truth (GT) 4DCT was available. A CIRT treatment plan was optimized at the end-exhale reference CT and the RBE-weighted dose distribution was recalculated on both the virtual and GT 4DCTs. The method estimation error was quantified by comparing the virtual and GT 4DCTs and the corresponding recomputed doses. The method was then evaluated on 8 patients with pancreas or liver tumors treated with CIRT using respiratory gating at end-exhale. The clinical treatment plans adopted at the National Center for Oncological Hadrontherapy (CNAO, Pavia, Italy) were considered and the dose distribution was recomputed on all respiratory phases of the planning and virtual 4DCTs. By comparing the two datasets and the corresponding dose distributions, the geometrical and dosimetric impact of organ motion was assessed.

RESULTS

For the phantom, the error outside of the 4DMRI FOV was up to 4.5mm, but it remained sub-millimetric in correspondence to the target within the 4DMRI FOV. Although the impact of motion on the target D_95% resulted in variations ranging from 22% to 90% between the planned dose and the doses recomputed on the GT 4DCT phases, the corresponding estimation error was ≤2.2%. In the patient cases, the variation of the baseline tumor position between the planning and the virtual end-exhale CTs presented a median (interquartile range) value of 6.0 (4.9) mm. For baseline variations larger than 5mm, the tumor D_95% variation between the plan and the dose recomputed on the end-exhale virtual CT resulted larger than 10%. Median variations higher than 10% in the target D_95% and gastro-intestinal OARs D_2% were quantified at the end-inhale, whereas close to the end-exhale phase, limited variations of relevant dose metrics were found for both tumor and OARs.

CONCLUSIONS

The negligible impact of the geometrical inaccuracy in the estimated anatomy outside of the 4DMRI FOV on the overall dosimetric accuracy suggests the feasibility of virtual 4DCT with extended FOV in CIRT of the abdominal site. In the analyzed patient group, inter-fractional variations such as baseline variation and breathing variability were quantified, demonstrating the method capability to support treatment planning in gated CIRT of the abdominal site.

Accuracy of real-time respiratory motion tracking and time delay of gating radiotherapy based on optical surface imaging technique

Background

Surface-guided radiation therapy (SGRT) employs a non-invasive real-time optical surface imaging (OSI) technique for patient surface motion monitoring during radiotherapy. The main purpose of this study is to verify the real-time tracking accuracy of SGRT for respiratory motion and provide a fitting method to detect the time delay of gating.

Methods

A respiratory motion phantom was utilized to simulate respiratory motion using 17 cosine breathing pattern curves with various periods and amplitudes. The motion tracking of the phantom was performed by the Catalyst™ system. The tracking accuracy of the system (with period and amplitude variations) was evaluated by analyzing the adjusted coefficient of determination (A_R²) and root mean square error (RMSE). Furthermore, 13 actual respiratory curves, which were categorized into regular and irregular patterns, were selected and then simulated by the phantom. The Fourier transform was applied to the respiratory curves, and tracking accuracy was compared through the quantitative analyses of curve similarity using the Pearson correlation coefficient (PCC). In addition, the time delay of amplitude-based respiratory-gating radiotherapy based on the OSI system with various beam hold times was tested using film dosimetry for the Elekta Versa-HD and Varian Edge linacs. A dose convolution-fitting method was provided to accurately measure the beam-on and beam-off time delays.

Results

A_R² and RMSE for the cosine curves were 0.9990–0.9996 and 0.110–0.241 mm for periods ranging from 1 s to 10 s and 0.9990–0.9994 and 0.059–0.175 mm for amplitudes ranging from 3 mm to 15 mm. The PCC for the actual respiratory curves ranged from 0.9955 to 0.9994, which was not significantly affected by breathing patterns. For gating radiotherapy, the average beam-on and beam-off time delays were 1664 ± 72 and 25 ± 30 ms for Versa-HD and 303 ± 45 and 34 ± 25 ms for Edge, respectively. The time delay was relatively stable as the beam hold time increased.

Conclusions

The OSI technique provides high accuracy for respiratory motion tracking. The proposed dose convolution-fitting method can accurately measure the time delay of respiratory-gating radiotherapy. When the OSI technique is used for respiratory-gating radiotherapy, the time delay for the beam-on is considerably longer than the beam-off.

Organ motion quantification and margins evaluation in carbon ion therapy of abdominal lesions

PURPOSE

In image-guided particle radiotherapy of abdominal lesions, respiratory motion hinders treatment accuracy. In this study, 2D cineMRI data were used to quantify the tumor (GTV) motion and to evaluate the clinical approach based on deriving an internal target volume (ITV) from a planning 4DCT for gating treatments.

METHODS

Seven patients with abdominal lesions were treated with carbon-ion therapy at the National Centre of Oncological Hadron-therapy (Italy). The MR scan was performed on the same day of the 4DCT acquisition. For four patients, an additional MR was acquired approximately after 1 week. The cineMRI combined with deformable image registration algorithm was used to quantify tumor motion. Afterwards, two ITVs were defined considering (1) all phases (ITVFB) and (2) only phases within the gating window (ITVG), and then compared with the clinical (4DCT-derived) ITVs (ITVCG and ITVCFB).

RESULTS

Tumor residual motion estimated by cineMRI data in the two MRI sessions resulted not significantly different from 4DCT, although cineMRI accounted for cycle-to-cycle variations. The ITV normalized for the GTV median values were higher for ITVFB with respect to ITVG, ITVCFB and ITVCG. The Hausdorff distances with respect to the GTV were up to 10.55 mm, 3.13 mm, 5.56 mm and 2.51 mm, for ITVFB, ITVG, ITVCFB and ITVCG, respectively. According to both metrics, ITVCG and ITVG were not found significantly different.

CONCLUSIONS

CineMRI acquisitions allowed to quantify organ motion without delivering additional dose to the patient and to verify treatment margins in gated carbon-ion therapy of abdominal lesions.

Perturbation analysis of 4D dose distribution for scanned carbon-ion beam radiotherapy

PURPOSE

To evaluate the patients' set-up error-induced perturbation effects on 4D dose distributions (4DDD) of range-adapted internal target volume-based (raITV) treatment plan using lung and liver 4DCT data sets.

METHODS

We enrolled 20 patients with lung and liver cancer treated with respiratory-gated carbon-ion beam scanning therapy. PTVs were generated by adding a 2 mm range-adapted set-up margin on the raITVs. Set-up errors were simulated by shifting the beam isocenter in three translational directions of ±2 mm, ±4 mm, and ±6 mm. 4DDDs were calculated for both nominal and isocenter-shifted situations. Dose metrics of CTV dose coverage (D95) and normal tissue sparing were evaluated. Statistical significance with p < 0.01 was considered by Wilcoxon signed rank test.

RESULTS

The CTV dose coverage was more sensitive to set-up errors for lung cases than for liver cases, and more serious in superior-inferior direction. The sufficient CTV-D95 > 98% could be achieved with set-up errors less than ±2 mm in all shift directions both for lung and liver cases. With the increase of set-up error, the CTV dose coverage decreased gradually. The clinical criterial of CTV-D95 > 95% could not be fulfilled with set-up error reached to ±4 mm for lung cases, and ±6 mm for liver cases. OAR doses did not have a significant difference with each set-up error for both lung and liver cases.

CONCLUSIONS

The range-adapted set-up margin successfully prevented dose degradation of 4DDDs in the presence of the same magnitude of set-up error for raITV-based carbon-ion beam scanning therapy.

A machine learning-based real-time tumor tracking system for fluoroscopic gating of lung radiotherapy

To improve respiratory-gated radiotherapy accuracy, we developed a machine learning approach for markerless tumor tracking and evaluated it using lung cancer patient data. Digitally reconstructed radiography (DRR) datasets were generated using planning 4DCT data. Tumor positions were selected on respective DRR images to place the GTV center of gravity in the center of each DRR. DRR subimages around the tumor regions were cropped so that the subimage size was defined by tumor size. Training data were then classified into two groups: positive (including tumor) and negative (not including tumor) samples. Machine learning parameters were optimized by the extremely randomized tree method. For the tracking stage, a machine learning algorithm was generated to provide a tumor likelihood map using fluoroscopic images. Prior probability tumor positions were also calculated using the previous two frames. Tumor position was then estimated by calculating maximum probability on the tumor likelihood map and prior probability tumor positions. We acquired treatment planning 4DCT images in eight patients. Digital fluoroscopic imaging systems on either side of the vertical irradiation port allowed fluoroscopic image acquisition during treatment delivery. Each fluoroscopic dataset was acquired at 15 frames per second. We evaluated the tracking accuracy and computation times. Tracking positional accuracy averaged over all patients was 1.03 ± 0.34 mm (mean ± standard deviation, Euclidean distance) and 1.76 ± 0.71 mm ([Formula: see text] percentile). Computation time was 28.66 ± 1.89 ms/frame averaged over all frames. Our markerless algorithm successfully estimated tumor position in real time.

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.

Virtual 4DCT from 4DMRI for the management of respiratory motion in carbon ion therapy of abdominal tumors

PURPOSE

To evaluate a method for generating virtual four-dimensional computed tomography (4DCT) from four-dimensional magnetic resonance imaging (4DMRI) data in carbon ion radiotherapy with pencil beam scanning for abdominal tumors.

METHODS

Deformable image registration is used to: (a) register each respiratory phase of the 4DMRI to the end-exhale MRI; (b) register the reference end-exhale CT to the end-exhale MRI volume; (c) generate the virtual 4DCT by warping the registered CT according to the obtained deformation fields. A respiratory-gated carbon ion treatment plan is optimized on the planning 4DCT and the corresponding dose distribution is recalculated on the virtual 4DCT. The method was validated on a digital anthropomorphic phantom and tested on eight patients (18 acquisitions). For the phantom, a ground truth dataset was available to assess the method performances from the geometrical and dosimetric standpoints. For the patients, the virtual 4DCT was compared with the planning 4DCT.

RESULTS

In the phantom, the method exhibits a geometrical accuracy within the voxel size and Dose Volume Histograms deviations up to 3.3% for target V_95% (mean dose difference ≤ 0.2% of the prescription dose, gamma pass rate > 98%). For patients, the virtual and the planning 4DCTs show good agreement at end-exhale (3% median D_95% difference), whereas other respiratory phases exhibit moderate motion variability with consequent dose discrepancies, confirming the need for motion mitigation strategies during treatment.

CONCLUSIONS

The virtual 4DCT approach is feasible to evaluate treatment plan robustness against intra- and interfraction motion in carbon ion therapy delivered at the abdominal site.

Real-time control of respiratory motion: Beyond radiation therapy

Motion management in radiation oncology is an important aspect of modern treatment planning and delivery. Special attention has been paid to control respiratory motion in recent years. However, other medical procedures related to both diagnosis and treatment are likely to benefit from the explicit control of breathing motion. Quantitative imaging - including increasingly important tools in radiology and nuclear medicine - is among the fields where a rapid development of motion control is most likely, due to the need for quantification accuracy. Emerging treatment modalities like focussed-ultrasound tumor ablation are also likely to benefit from a significant evolution of motion control in the near future. In the present article an overview of available respiratory motion systems along with ongoing research in this area is provided. Furthermore, an attempt is made to envision some of the most expected developments in this field in the near future.

Preoperative chemotherapy and carbon ions therapy for treatment of resectable and borderline resectable pancreatic adenocarcinoma: a prospective, phase II, multicentre, single-arm study

Background

Pancreatic adenocarcinoma is a high-mortality neoplasm with a documented 5-years-overall survival around 5%. In the last decades, a real breakthrough in the treatment of the disease has not been achieved. Here we propose a prospective, phase II, multicentre, single-arm study aiming to assess the efficacy and the feasibility of a therapeutic protocol combining chemotherapy, carbon ion therapy and surgery for resectable and borderline resectable pancreatic adenocarcinoma.

Method

The purpose of this trial (PIOPPO Protocol) is to assess the efficacy and the feasibility of 3 cycles of FOLFIRINOX neoadjuvant chemotherapy followed by a short-course of carbon ion radiotherapy (CIRT) for resectable or borderline resectable pancreatic adenocarcinoma patients. Primary outcome of this study is the assessment of local progression free survival (L-PFS). The calculation of sample size is based on the analysis of the primary endpoint “progression free survival” according to Fleming’s Procedure.

Discussion

Very preliminary results provide initial evidence of the feasibility of the combined chemotherapy and CIRT in the neoadjuvant setting for resectable or borderline resectable pancreatic cancer. Completion of the accrual and long term results are awaited to see if this combination of treatment is advisable and will provide the expected benefits.

Trial registration

ClinicalTrials.gov Identifier: NCT03822936 registered on January 2019.

Real-time intrafraction motion monitoring in external beam radiotherapy

Radiotherapy (RT) aims to deliver a spatially conformal dose of radiation to tumours while maximizing the dose sparing to healthy tissues. However, the internal patient anatomy is constantly moving due to respiratory, cardiac, gastrointestinal and urinary activity. The long term goal of the RT community to 'see what we treat, as we treat' and to act on this information instantaneously has resulted in rapid technological innovation. Specialized treatment machines, such as robotic or gimbal-steered linear accelerators (linac) with in-room imaging suites, have been developed specifically for real-time treatment adaptation. Additional equipment, such as stereoscopic kilovoltage (kV) imaging, ultrasound transducers and electromagnetic transponders, has been developed for intrafraction motion monitoring on conventional linacs. Magnetic resonance imaging (MRI) has been integrated with cobalt treatment units and more recently with linacs. In addition to hardware innovation, software development has played a substantial role in the development of motion monitoring methods based on respiratory motion surrogates and planar kV or Megavoltage (MV) imaging that is available on standard equipped linacs. In this paper, we review and compare the different intrafraction motion monitoring methods proposed in the literature and demonstrated in real-time on clinical data as well as their possible future developments. We then discuss general considerations on validation and quality assurance for clinical implementation. Besides photon RT, particle therapy is increasingly used to treat moving targets. However, transferring motion monitoring technologies from linacs to particle beam lines presents substantial challenges. Lessons learned from the implementation of real-time intrafraction monitoring for photon RT will be used as a basis to discuss the implementation of these methods for particle RT.

Validation of a model for physical dose variations in irregularly moving targets treated with carbon ion beams

PURPOSE

In particle therapy, conventional treatment planning systems rely on an imaging representation of the irradiated region to compute the dose. For irregular breathing, when an imaging dataset describing the actual motion is not available, a different approach for dose estimation is needed. To this aim, we validate a method for the estimation of physical dose variations in gated carbon ion treatments, providing also a demonstration of the feasibility of physical dose metrics to assess the method performance. Finally, we describe a sample use case, in which this method is used to assess plan robustness with respect to undetected irregular tumor motion.

METHODS

The method entails the definition of a patient- and beam-specific water equivalent depth (WED) space, the simulation of motion as a translation equal to tumor displacement, and the reconstruction of the altered dose. We validated the approach using four-dimensional computed tomographies (4DCTs) and clinical plans in 12 patients, treated with respiratory gated carbon ion beams at the National Centre for Oncological Hadrontherapy (Pavia, Italy). Using the end-exhale CT and dose distribution as a reference, the physical dose delivered at the end-inhale tumor position was estimated and compared to the ground-truth dose recalculation on the end-inhale CT. Biologically effective and physical dose variations between the plan and the recalculation were compared as well. As a use case, we evaluated dose changes caused by simulated irregular tumor motion, that is, linear and nonlinear baseline shifts and/or amplitude variations with hysteresis.

RESULTS

The ratio between biologically effective and physical equivalent uniform dose (EUD) variations due to end-exhale to end-inhale motion was less than one for 96% of investigated structures. In the validation study, we found a median error corresponding to a 14% EUD overestimation for the tumor and 4% EUD underestimation for a subgroup of organs at risk, together with a high EUD variation due to motion [median 352% EUD variation between end-exhale and end-inhale doses in the planning tumor volume (PTV)]. Considering relevant dose-volume histogram (DVH) metrics, the median difference between estimated and ground truth doses was ≤ 4%. Gamma analysis between estimated and recalculated dose distributions resulted in a pass rate > 80% for 83% of the target volumes. For the two patients selected for the sample use case, a patient-specific assessment of the method performance was performed on the 4DCT and it was possible to relate EUD variations of both tumor and organs at risk to the simulated target motion.

CONCLUSIONS

The physical dose distribution was found to be more sensitive to motion with respect to the biologically effective one, suggesting the suitability of the physical dose metrics for the WED-space method validation. We showed that the method can compensate for intra-fractional tumor motion with proper accuracy in the selected patient group, although its use is recommended when limited deformations are expected. In conclusion, the WED-space method can provide simulations of dose alteration due to irregular breathing when imaging data are lacking, and, once integrated with relative biological effectiveness (RBE) modeling, it would be useful in evaluating the robustness of carbon ion treatment plans.

Experimental verification of a two-dimensional respiratory motion compensation system with ultrasound tracking technique in radiation therapy

This study proposed respiratory motion compensation system (RMCS) combined with an ultrasound image tracking algorithm (UITA) to compensate for respiration-induced tumor motion during radiotherapy, and to address the problem of inaccurate radiation dose delivery caused by respiratory movement. This study used an ultrasound imaging system to monitor respiratory movements combined with the proposed UITA and RMCS for tracking and compensation of the respiratory motion. Respiratory motion compensation was performed using prerecorded human respiratory motion signals and also sinusoidal signals. A linear accelerator was used to deliver radiation doses to GAFchromic EBT3 dosimetry film, and the conformity index (CI), root-mean-square error, compensation rate (CR), and planning target volume (PTV) were used to evaluate the tracking and compensation performance of the proposed system. Human respiratory pattern signals were captured using the UITA and compensated by the RMCS, which yielded CR values of 34-78%. In addition, the maximum coronal area of the PTV ranged from 85.53 mm² to 351.11 mm² (uncompensated), which reduced to from 17.72 mm² to 66.17 mm² after compensation, with an area reduction ratio of up to 90%. In real-time monitoring of the respiration compensation state, the CI values for 85% and 90% isodose areas increased to 0.7 and 0.68, respectively. The proposed UITA and RMCS can reduce the movement of the tracked target relative to the LINAC in radiation therapy, thereby reducing the required size of the PTV margin and increasing the effect of the radiation dose received by the treatment target.