Online daily adaptive proton therapy

An in-vivo treatment monitoring system for ion-beam radiotherapy based on 28 Timepix3 detectors

Ion-beam radiotherapy is an advanced cancer treatment modality offering steep dose gradients and a high biological effectiveness. These gradients make the therapy vulnerable to patient-setup and anatomical changes between treatment fractions, which may go unnoticed. Charged fragments from nuclear interactions of the ion beam with the patient tissue may carry information about the treatment quality. Currently, the fragments escape the patient undetected. Inter-fractional in-vivo treatment monitoring based on these charged nuclear fragments could make ion-beam therapy safer and more efficient. We developed an ion-beam monitoring system based on 28 hybrid silicon pixel detectors (Timepix3) to measure the distribution of fragment origins in three dimensions. The system design choices as well as the ion-beam monitoring performance measurements are presented in this manuscript. A spatial resolution of 4 mm along the beam axis was achieved for the measurement of individual fragment origins. Beam-range shifts of1.5 mm were identified in a clinically realistic treatment scenario with an anthropomorphic head phantom. The monitoring system is currently being used in a prospective clinical trial at the Heidelberg Ion Beam Therapy Centre for head-and-neck as well as central nervous system cancer patients.

The partial adaptation strategy for online-adaptive proton therapy: A proof of concept study in head and neck cancer patients

BACKGROUND

The accuracy of intensity-modulated proton therapy (IMPT) is greatly affected by anatomy variations that might occur during the treatment course. Online plan adaptations have been proposed as a solution to intervene promptly during a treatment session once the anatomy changes are detected. The implementation of online-adaptive proton therapy (OAPT) is still hindered by time-consuming tasks in the workflow.

PURPOSE

The study introduces the novel concept of partial adaptation and aims at investigating its feasibility as a potential solution to parallelize tasks during an OAPT workflow for saving valuable in-room time.

METHODS

The proof-of-principle simulation study includes datasets from six head and neck cancer (HNC) patients, each consisting of one planning CT (pCT) and three contoured control CTs (cCTs). Robust 3-field normo-fractionated initial IMPT plans were generated on the pCTs with a standardized field configuration, delivering 66 Gy and 54 Gy to the high-risk and low-risk clinical target volume (CTVHigh and CTVLow), respectively. For each cCT, a dose-mimicking-based partial adaptation was applied: two fields were adapted on the current anatomy taking into account the background dose of the first non-adapted field supposedly delivered in the meantime. Fraction doses on the cCTs resulting from partially adapted plans with different first (non-adapted) field assignments were compared against those from non-adapted and fully adapted plans regarding target coverage and organs at risk (OARs) sparing. The robustness of partially adapted plans was also evaluated.

RESULTS

Partially adapted plans showed comparable results to fully adapted plans and were superior to non-adapted plans for both target coverage and OAR sparing. Target coverage degradation in the non-adapted plans (median D98%: 95.9% and 97.5% for CTVLow and CTVHigh, respectively) was recovered by both partial (98.0% and 98.5%) and full adaptation (98.2% and 98.7%) in comparison to the initial plans (98.7% and 98.8%). The initial hotspot dose for the CTVHigh (median D2%: 101.8%) increased in the non-adapted plans (102.9%) and was recovered by the adaptive strategies (partial: 102.5%, full: 101.9%). The near-maximum dose (D0.01cc) to brainstem and spinal cord was within clinical constraints for all investigated dose distributions, but clearly increased for no adaptation and improved in the (both partially and fully) adapted plans with respect to the non-adapted ones. The parotids' median doses (D50) were mainly patient-specific depending on the proximity to the target region, but anyway lower for the partially and fully adapted plans compared to the non-adapted ones. OAR sparing was furthermore improved for the partially adapted plans in comparison to full adaptation. Robustness of the target dose metrics was preserved in all evaluated scenarios.

CONCLUSIONS

For OAPT of HNC patients, partial adaptation is able to generate plans of superior conformity to non-adapted plans and of comparable conformity as fully adapted plans, while having the potential to speed up the online-adaptive workflows. Thus, partial adaptation represents an intermediate approach until fast online adaptation workflows become available. Furthermore, it can be applied in workflows where online treatment verification stops the delivery and triggers an online adaptation for the remaining fraction.

A reference-point-method-based online proton treatment plan re-optimization strategy and a novel solution to planning constraint infeasibility problem

Objective. Propose a highly automated treatment plan re-optimization strategy suitable for online adaptive proton therapy. The strategy includes a rapid re-optimization method that generates quality replans and a novel solution that efficiently addresses the planning constraint infeasibility issue that can significantly prolong the re-optimization process.Approach. We propose a systematic reference point method (RPM) model that minimizes the l-infinity norm from the initial treatment plan in the daily objective space for online re-optimization. This model minimizes the largest objective value deviation among the objectives of the daily replan from their reference values, leading to a daily replan similar to the initial plan. Whether a set of planning constraints is feasible with respect to the daily anatomy cannot be known before solving the corresponding optimization problem. The conventional trial-and-error-based relaxation process can cost a significant amount of time. To that end, we propose an optimization problem that first estimates the magnitude of daily violation of each planning constraint. Guided by the violation magnitude and clinical importance of the constraints, the constraints are then iteratively converted into objectives based on their priority until the infeasibility issue is solved.Main results.The proposed RPM-based strategy generated replans similar to the offline manual replans within the online time requirement for six head and neck and four breast patients. The average targetD95and relevant organ at risk sparing parameter differences between the RPM replans and clinical offline replans were -0.23, -1.62 Gy for head and neck cases and 0.29, -0.39 Gy for breast cases. The proposed constraint relaxation solution made the RPM problem feasible after one round of relaxation for all four patients who encountered the infeasibility issue.Significance. We proposed a novel RPM-based re-optimization strategy and demonstrated its effectiveness on complex cases, regardless of whether constraint infeasibility is encountered.

Uncertainty estimation and evaluation of deformation image registration based convolutional neural networks

Objective.Fast and accurate deformable image registration (DIR), including DIR uncertainty estimation, is essential for safe and reliable clinical deployment. While recent deep learning models have shown promise in predicting DIR with its uncertainty, challenges persist in proper uncertainty evaluation and hyperparameter optimization for these methods. This work aims to develop and evaluate a model that can perform fast DIR and predict its uncertainty in seconds.Approach.This study introduces a novel probabilistic multi-resolution image registration model utilizing convolutional neural networks to estimate a multivariate normal distributed dense displacement field (DDF) in a multimodal image registration problem. To assess the quality of the DDF distribution predicted by the model, we propose a new metric based on the Kullback-Leibler divergence. The performance of our approach was evaluated against three other DIR algorithms (VoxelMorph, Monte Carlo dropout, and Monte Carlo B-spline) capable of predicting uncertainty. The evaluation of the models included not only the quality of the deformation but also the reliability of the estimated uncertainty. Our application investigated the registration of a treatment planning computed tomography (CT) to follow-up cone beam CT for daily adaptive radiotherapy.Main results.The hyperparameter tuning of the models showed a trade-off between the estimated uncertainty's reliability and the deformation's accuracy. In the optimal trade-off, our model excelled in contour propagation and uncertainty estimation (p <0.05) compared to existing uncertainty estimation models. We obtained an average dice similarity coefficient of 0.89 and a KL-divergence of 0.15.Significance.By addressing challenges in DIR uncertainty estimation and evaluation, our work showed that both the DIR and its uncertainty can be reliably predicted, paving the way for safe deployment in a clinical environment.

Extension of the cone-beam CT field-of-view using two complementary short scans

BACKGROUND

Robotic C-arm cone-beam computed tomography (CBCT) scanners provide fast in-room imaging in radiotherapy. Their mobility extends beyond performing a gantry rotation, but they might encounter obstructions to their motion which limit the gantry angle range. The axial field-of-view (FOV) of a reconstructed CBCT image depends on the acquisition geometry. When imaging a large anatomical location, such as the thorax, abdomen, or pelvis, a centered cone beam might be insufficient to acquire untruncated projection images. Some CBCT scanners can laterally displace their detector and collimate the beam to increase the FOV, but the gantry must then perform a360 ∘ $360^{\circ}$ rotation to provide complete data for reconstruction.

PURPOSE

To extend the FOV of a CBCT image with a single short scan (gantry angle range of180 ∘ + $180^{\circ}+$ fan angle) using two complementary short scans.

METHODS

We defined an acquisition protocol using two short scans during which the source follows the same trajectory and where the detector has equal and opposite tilt and/or offset between the two scans, which we refer to as complementary scans. We created virtual acquisitions using a Monte Carlo simulator on a digital anthropomorphic phantom and on a computed tomography (CT) scan of a patient abdomen. For our proposed method, each simulation produced two complementary sets of projections, which were weighted for redundancies and used to reconstruct one CBCT image. We compared the resulting images to the ground truth phantoms and simulations of conventional scans.

RESULTS

Reconstruction artifacts were slightly more prominent in the complementary scans w.r.t. a complete scan with untruncated projections but matched those in a single short scan without truncation. When analyzing reconstructed scans from simulated projections with scatter and corrected with prior CT information, we found a global agreement between complementary and conventional scan approaches.

CONCLUSIONS

When dealing with a limited range of motion of the gantry of a CBCT scanner, two complementary short scans are a technically valid alternative to a full 360∘ $^{\circ}$ scan with equal FOV. This approach enables FOV extension without collisions or hardware upgrades.

A retrospective study on the investigation of potential dosimetric benefits of online adaptive proton therapy for head and neck cancer

PURPOSE

Proton therapy is sensitive to anatomical changes, often occurring in head-and-neck (HN) cancer patients. Although multiple studies have proposed online adaptive proton therapy (APT), there is still a concern in the radiotherapy community about the necessity of online APT. We have performed a retrospective study to investigate the potential dosimetric benefits of online APT for HN patients relative to the current offline APT.

METHODS

Our retrospective study has a patient cohort of 10 cases. To mimic online APT, we re-evaluated the dose of the in-use treatment plan on patients' actual treatment anatomy captured by cone-beam CT (CBCT) for each fraction and performed a templated-based automatic replanning if needed, assuming that these were performed online before treatment delivery. Cumulative dose of the simulated online APT course was calculated and compared with that of the actual offline APT course and the designed plan dose of the initial treatment plan (referred to as nominal plan). The ProKnow scoring system was employed and adapted for our study to quantify the actual quality of both courses against our planning goals.

RESULTS

The average score of the nominal plans over the 10 cases is 41.0, while those of the actual offline APT course and our simulated online course is 25.8 and 37.5, respectively. Compared to the offline APT course, our online course improved dose quality for all cases, with the score improvement ranging from 0.4 to 26.9 and an average improvement of 11.7.

CONCLUSION

The results of our retrospective study have demonstrated that online APT can better address anatomical changes for HN cancer patients than the current offline replanning practice. The advanced artificial intelligence based automatic replanning technology presents a promising avenue for extending potential benefits of online APT.

Feasibility of Monte Carlo dropout-based uncertainty maps to evaluate deep learning-based synthetic CTs for adaptive proton therapy

BACKGROUND

Deep learning has shown promising results to generate MRI-based synthetic CTs and to enable accurate proton dose calculations on MRIs. For clinical implementation of synthetic CTs, quality assurance tools that verify their quality and reliability are required but still lacking.

PURPOSE

This study aims to evaluate the predictive value of uncertainty maps generated with Monte Carlo dropout (MCD) for verifying proton dose calculations on deep-learning-based synthetic CTs (sCTs) derived from MRIs in online adaptive proton therapy.

METHODS

Two deep-learning models (DCNN and cycleGAN) were trained for CT image synthesis using 101 paired CT-MR images. sCT images were generated using MCD for each model by performing 10 inferences with activated dropout layers. The final sCT was obtained by averaging the inferred sCTs, while the uncertainty map was obtained from the HU variance corresponding to each voxel of 10 sCTs. The resulting uncertainty maps were compared to the observed HU-, range-, WET-, and dose-error maps between the sCT and planning CT. For range and WET errors, the generated uncertainty maps were projected along the 90-degree angle. To evaluate the dose distribution, a mask based on the 5%-isodose curve was applied to only include voxels along the beam paths. Pearson's correlation coefficients were calculated to determine the correlation between the uncertainty maps and HUs, range, WET, and dose errors. To evaluate the dosimetric accuracy of synthetic CTs, clinical proton treatment plans were recalculated and compared to the pCTs RESULTS: Evaluation of the correlation showed an average of r = 0.92 ± 0.03 and r = 0.92 ± 0.03 for errors between uncertainty-HU, r = 0.66 ± 0.09 and r = 0.62 ± 0.06 between uncertainty-range, r = 0.64 ± 0.06 and r = 0.58 ± 0.07 between uncertainty-WET, and r = 0.65 ± 0.09 and r = 0.67 ± 0.07 between uncertainty and dose difference for DCNN and cycleGAN model, respectively. Dosimetric comparison for target volumes showed an average 3%/3 mm gamma pass rate of 99.76 ± 0.43 (DCNN) and 99.10 ± 1.27 (cycleGAN).

CONCLUSION

The observed correlations between uncertainty maps and the various metrics (HU, range, WET, and dose errors) demonstrated the potential of MCD-based uncertainty maps as a reliable QA tool to evaluate the accuracy of deep learning-based sCTs.

Adaptive radiotherapy for head and neck cancer: Pitfalls and possibilities from the radiation oncologist's point of view

BACKGROUND

Patients with head and neck cancer (HNC) may experience substantial anatomical changes during the course of radiotherapy treatment. The implementation of adaptive radiotherapy (ART) proves effective in managing the consequent impact on the planned dose distribution.

METHODS

This narrative literature review comprehensively discusses the diverse strategies of ART in HNC and the documented dosimetric and clinical advantages associated with these approaches, while also addressing the current challenges for integration of ART into clinical practice.

RESULTS AND CONCLUSION

Although based on mainly non-randomized and retrospective trials, there is accumulating evidence that ART has the potential to reduce toxicity and improve quality of life and tumor control in HNC patients treated with RT. However, several questions remain regarding accurate patient selection, the ideal frequency and timing of replanning, and the appropriate way for image registration and dose calculation. Well-designed randomized prospective trials, with a predetermined protocol for both image registration and dose summation, are urgently needed to further investigate the dosimetric and clinical benefits of ART.

Using the gamma-index analysis for inter-fractional comparison of in-beam PET images for head-and-neck treatment monitoring in proton therapy: A Monte Carlo simulation study

GOAL

In-beam Positron Emission Tomography (PET) is a technique for in-vivo non-invasive treatment monitoring for proton therapy. To detect anatomical changes in patients with PET, various analysis methods exist, but their clinical interpretation is problematic. The goal of this work is to investigate whether the gamma-index analysis, widely used for dose comparisons, is an appropriate tool for comparing in-beam PET distributions. Focusing on a head-and-neck patient, we investigate whether the gamma-index map and the passing rate are sensitive to progressive anatomical changes.

METHODS/MATERIALS

We simulated a treatment course of a proton therapy patient using FLUKA Monte Carlo simulations. Gradual emptying of the sinonasal cavity was modeled through a series of artificially modified CT scans. The in-beam PET activity distributions from three fields were evaluated, simulating a planar dual head geometry. We applied the 3D-gamma evaluation method to compare the PET images with a reference image without changes. Various tolerance criteria and parameters were tested, and results were compared to the CT-scans.

RESULTS

Based on 210 MC simulations we identified appropriate parameters for the gamma-index analysis. Tolerance values of 3 mm/3% and 2 mm/2% were suited for comparison of simulated in-beam PET distributions. The gamma passing rate decreased with increasing volume change for all fields.

CONCLUSION

The gamma-index analysis was found to be a useful tool for comparing simulated in-beam PET images, sensitive to sinonasal cavity emptying. Monitoring the gamma passing rate behavior over the treatment course is useful to detect anatomical changes occurring during the treatment course.

Synthetic CT imaging for PET monitoring in proton therapy: a simulation study

Objective.This study addresses a fundamental limitation of in-beam positron emission tomography (IB-PET) in proton therapy: the lack of direct anatomical representation in the images it produces. We aim to overcome this shortcoming by pioneering the application of deep learning techniques to create synthetic control CT images (sCT) from combining IB-PET and planning CT scan data.Approach.We conducted simulations involving six patients who underwent irradiation with proton beams. Leveraging the architecture of a visual transformer (ViT) neural network, we developed a model to generate sCT images of these patients using the planning CT scans and the inter-fractional simulated PET activity maps during irradiation. To evaluate the model's performance, a comparison was conducted between the sCT images produced by the ViT model and the authentic control CT images-serving as the benchmark.Main results.The structural similarity index was computed at a mean value across all patients of 0.91, while the mean absolute error measured 22 Hounsfield Units (HU). Root mean squared error and peak signal-to-noise ratio values were 56 HU and 30 dB, respectively. The Dice similarity coefficient exhibited a value of 0.98. These values are comparable to or exceed those found in the literature. More than 70% of the synthetic morphological changes were found to be geometrically compatible with the ones reported in the real control CT scan.Significance.Our study presents an innovative approach to surface the hidden anatomical information of IB-PET in proton therapy. Our ViT-based model successfully generates sCT images from inter-fractional PET data and planning CT scans. Our model's performance stands on par with existing models relying on input from cone beam CT or magnetic resonance imaging, which contain more anatomical information than activity maps.

Generation and evaluation of anatomy-preserving virtual CT for online adaptive proton therapy

BACKGROUND

Daily CTs generated by CBCT correction are required for daily replanning in online-adaptive proton therapy (APT) to effectively deal with inter-fractional changes. Out of the currently available methods, the suitability of a daily CT generation method for proton dose calculation also depends on the anatomical site.

PURPOSE

We propose an anatomy-preserving virtual CT (APvCT) method as a hybrid method of CBCT correction, which is especially suitable for large anatomy deformations. The accuracy of the hybrid method was assessed by comparison with the corrected CBCT (cCBCT) and virtual CT (vCT) methods in the context of online APT.

METHODS

Seventy-one daily CBCTs of four prostate cancer patients treated with intensity modulated proton therapy (IMPT) were converted to daily CTs using cCBCT, vCT, and the newly proposed APvCT method. In APvCT, planning CT (pCT) were mapped to CBCT geometry using deformable image registration with boundary conditions on controlling regions of interest (ROIs) created with deep learning segmentation on cCBCT. The relative frequency distribution (RFD) of HU, mass density and stopping power ratio (SPR) values were assessed and compared with the pCT. The ROIs in the APvCT and vCT were compared with cCBCT in terms of Dice similarity coefficient (DSC) and mean distance-to-agreement (mDTA). For each patient, a robustly optimized IMPT plan was created on the pCT and subsequent daily adaptive plans on daily CTs. For dose distribution comparison on the same anatomy, the daily adaptive plans on cCBCT and vCT were recalculated on the corresponding APvCT. The dose distributions were compared in terms of isodose volumes and 3D global gamma-index passing rate (GPR) at γ(2%, 2 mm) criterion.

RESULTS

For all patients, no noticeable difference in RFDs was observed amongst APvCT, vCT, and pCT except in cCBCT, which showed a noticeable difference. The minimum DSC value was 0.96 and 0.39 for contours in APvCT and vCT respectively. The average value of mDTA for APvCT was 0.01 cm for clinical target volume and ≤0.01 cm for organs at risk, which increased to 0.18 cm and ≤0.52 cm for vCT. The mean GPR value was 90.9%, 64.5%, and 67.0% for APvCT versus cCBCT, vCT versus cCBCT, and APvCT versus vCT, respectively. When recalculated on APvCT, the adaptive cCBCT and vCT plans resulted in mean GPRs of 89.5 ± 5.1% and 65.9 ± 19.1%, respectively. The mean DSC values for 80.0%, 90.0%, 95.0%, 98.0%, and 100.0% isodose volumes were 0.97, 0.97, 0.97, 0.95, and 0.91 for recalculated cCBCT plans, and 0.89, 0.88, 0.87, 0.85, and 0.81 for recalculated vCT plans. Hausdorff distance for the 100.0% isodose volume in some cases of recalculated cCBCT plans on APvCT exceeded 1.00 cm.

CONCLUSIONS

APvCT contours showed good agreement with reference contours of cCBCT which indicates anatomy preservation in APvCT. A vCT with erroneous anatomy can result in an incorrect adaptive plan. Further, slightly lower values of GPR between the APvCT and cCBCT-based adaptive plans can be explained by the difference in the cCBCT's SPR RFD from the pCT.

Uncertainty-aware MR-based CT synthesis for robust proton therapy planning of brain tumour

BACKGROUND AND PURPOSE

Deep learning techniques excel in MR-based CT synthesis, but missing uncertainty prediction limits its clinical use in proton therapy. We developed an uncertainty-aware framework and evaluated its efficiency in robust proton planning.

MATERIALS AND METHODS

A conditional generative-adversarial network was trained on 64 brain tumour patients with paired MR-CT images to generate synthetic CTs (sCT) from combined T1-T2 MRs of three orthogonal planes. A Bayesian neural network predicts Laplacian distributions for all voxels with parameters (μ, b). A robust proton plan was optimized using three sCTs of μ and μ±b. The dosimetric differences between the plan from sCT (sPlan) and the recalculated plan (rPlan) on planning CT (pCT) were quantified for each patient. The uncertainty-aware robust plan was compared to conventional robust (global ± 3 %) and non-robust plans.

RESULTS

In 8-fold cross-validation, sCT-pCT image differences (Mean-Absolute-Error) were 80.84 ± 9.84HU (body), 35.78 ± 6.07HU (soft tissues) and 221.88 ± 31.69HU (bones), with Dice scores of 90.33 ± 2.43 %, 95.13 ± 0.80 %, and 85.53 ± 4.16 %, respectively. The uncertainty distribution positively correlated with absolute prediction error (Correlation Coefficient: 0.62 ± 0.01). The uncertainty-conditioned robust optimisation improved the rPlan-sPlan agreement, e.g., D95 absolute difference (CTV) was 1.10 ± 1.24 % compared to conventional (1.64 ± 2.71 %) and non-robust (2.08 ± 2.96 %) optimisation. This trend was consistent across all target and organs-at-risk indexes.

CONCLUSION

The enhanced framework incorporates 3D uncertainty prediction and generates high-quality sCTs from MR images. The framework also facilitates conditioned robust optimisation, bolstering proton plan robustness against network prediction errors. The innovative feature of uncertainty visualisation and robust analyses contribute to evaluating sCT clinical utility for individual patients.

Proton Pencil-Beam Scanning Stereotactic Body Radiation Therapy and Hypofractionated Radiation Therapy for Thoracic Malignancies: Patterns of Practice Survey and Recommendations for Future Development from NRG Oncology and PTCOG

Stereotactic body radiation therapy (SBRT) and hypofractionation using pencil-beam scanning (PBS) proton therapy (PBSPT) is an attractive option for thoracic malignancies. Combining the advantages of target coverage conformity and critical organ sparing from both PBSPT and SBRT, this new delivery technique has great potential to improve the therapeutic ratio, particularly for tumors near critical organs. Safe and effective implementation of PBSPT SBRT/hypofractionation to treat thoracic malignancies is more challenging than the conventionally-fractionated PBSPT due to concerns of amplified uncertainties at the larger dose per fraction. NRG Oncology and Particle Therapy Cooperative Group (PTCOG) Thoracic Subcommittee surveyed US proton centers to identify practice patterns of thoracic PBSPT SBRT/hypofractionation. From these patterns, we present recommendations for future technical development of proton SBRT/hypofractionation for thoracic treatment. Amongst other points, the recommendations highlight the need for volumetric image guidance and multiple CT-based robust optimization and robustness tools to minimize further the impact of uncertainties associated with respiratory motion. Advances in direct motion analysis techniques are urgently needed to supplement current motion management techniques.

The Role of Proton Therapy for Prostate Cancer in the Setting of Hip Prosthesis

PURPOSE

Given that the current standard of proton therapy (PT) for prostate cancer is through bilateral beams, this modality is typically avoided when it comes to treatment of patients with hip prosthesis. The purpose of this study was to evaluate whether novel PT methods, i.e., anterior proton beams and proton arc therapy (PArc), could be feasible options to treat this patient subpopulation. We evaluate PT methods in the context of dosimetry and robustness and compare with standard of practice volumetric modulated arc therapy (VMAT) to explore any potential benefits.

METHODS

Two PT and one VMAT treatment plans were retrospectively created for 10 patients who participated in a clinical trial with a weekly repeat CT (rCT) imaging component. All plans were robustly optimized and featured: (1) combination anterior oblique and lateral proton beams (AoL), (2) PArc, and (3) VMAT. All patients had hydrogel spacers in place, which enabled safe application of anterior proton beams. The planned dose was 70 Gy (RBE) to the entire prostate gland and 50 Gy (RBE) to the proximal seminal vesicles in 28 fractions. Along with plan dose-volume metrics, robustness to setup and interfractional variations were evaluated using the weekly rCT images. The linear energy transfer (LET)-weighted dose was evaluated for PArc plans to ensure urethra sparing given the typical high-LET region at the end of range.

RESULTS

Both PT methods were dosimetrically feasible and provided reduction of some key OAR metrics compared to VMAT except for penile bulb, while providing equally good target coverage. Significant differences in median rectum V35 (22-25%), penile bulb Dmean (5 Gy), rectum V61 (2%), right femoral head Dmean (5 Gy), and bladder V39 (4%) were found between PT and VMAT. All plans were equally robust to variations. LET-weighted dose in urethra was equivalent to the physical dose for PArc plans and hence no added urethral toxicity was expected.

CONCLUSIONS

PT for treatment of prostate cancer patients with hip prosthesis is feasible and equivalent or potentially superior to VMAT in quality in some cases. The choice of radiotherapy regimen can be personalized based on patient characteristics to achieve the best treatment outcome.

Optically stimulated luminescence dosimeters for simultaneous measurement of point dose and dose-weighted LET in an adaptive proton therapy workflow

Purpose

To demonstrate the suitability of optically stimulated luminescence detectors (OSLDs) for accurate simultaneous measurement of the absolute point dose and dose-weighted linear energy transfer (LETD) in an anthropomorphic phantom for experimental validation of daily adaptive proton therapy.

Methods

A clinically realistic intensity-modulated proton therapy (IMPT) treatment plan was created based on a CT of an anthropomorphic head-and-neck phantom made of tissue-equivalent material. The IMPT plan was optimized with three fields to deliver a uniform dose to the target volume covering the OSLDs. Different scenarios representing inter-fractional anatomical changes were created by modifying the phantom. An online adaptive proton therapy workflow was used to recover the daily dose distribution and account for the applied geometry changes. To validate the adaptive workflow, measurements were performed by irradiating Al2O3:C OSLDs inside the phantom. In addition to the measurements, retrospective Monte Carlo simulations were performed to compare the absolute dose and dose-averaged LET (LETD) delivered to the OSLDs.

Results

The online adaptive proton therapy workflow was shown to recover significant degradation in dose conformity resulting from large anatomical and positioning deviations from the reference plan. The Monte Carlo simulations were in close agreement with the OSLD measurements, with an average relative error of 1.4% for doses and 3.2% for LETD. The use of OSLDs for LET determination allowed for a correction for the ionization quenched response.

Conclusion

The OSLDs appear to be an excellent detector for simultaneously assessing dose and LET distributions in proton irradiation of an anthropomorphic phantom. The OSLDs can be cut to almost any size and shape, making them ideal for in-phantom measurements to probe the radiation quality and dose in a predefined region of interest. Although we have presented the results obtained in the experimental validation of an adaptive proton therapy workflow, the same approach can be generalized and used for a variety of clinical innovations and workflow developments that require accurate assessment of point dose and/or average LET.

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.

Near-field coded-mask technique and its potential for proton therapy monitoring

Objective.Prompt-gamma imaging encompasses several approaches to the online monitoring of the beam range or deposited dose distribution in proton therapy. We test one of the imaging techniques - a coded mask approach - both experimentally and via simulations.Approach.Two imaging setups have been investigated experimentally. Each of them comprised a structured tungsten collimator in the form of a modified uniformly redundant array mask and a LYSO:Ce scintillation detector of fine granularity. The setups differed in detector dimensions and operation mode (1D or 2D imaging). A series of measurements with radioactive sources have been conducted, testing the performance of the setups for near-field gamma imaging. Additionally, Monte Carlo simulations of a larger setup of the same type were conducted, investigating its performance with a realistic gamma source distribution occurring during proton therapy.Main results.The images of point-like sources reconstructed from two small-scale prototypes' data using the maximum-likelihood expectation maximisation algorithm constitute the experimental proof of principle for the near-field coded-mask imaging modality, both in the 1D and the 2D mode. Their precision allowed us to calibrate out certain systematic offsets appearing due to the limited alignment accuracy of setup elements. The simulation of the full-scale setup yielded a mean distal falloff retrieval precision of 0.72 mm in the studies for beam energy range 89.5-107.9 MeV and with 1 × 108protons (a typical number for distal spots). The implemented algorithm of image reconstruction is relatively fast-a typical procedure needs several seconds.Significance.Coded-mask imaging appears a valid option for proton therapy monitoring. The results of simulations let us conclude that the proposed full-scale setup is competitive with the knife-edge-shaped and the multi-parallel slit cameras investigated by other groups.

Review and recommendations on deformable image registration uncertainties for radiotherapy applications

Deformable image registration (DIR) is a versatile tool used in many applications in radiotherapy (RT). DIR algorithms have been implemented in many commercial treatment planning systems providing accessible and easy-to-use solutions. However, the geometric uncertainty of DIR can be large and difficult to quantify, resulting in barriers to clinical practice. Currently, there is no agreement in the RT community on how to quantify these uncertainties and determine thresholds that distinguish a good DIR result from a poor one. This review summarises the current literature on sources of DIR uncertainties and their impact on RT applications. Recommendations are provided on how to handle these uncertainties for patient-specific use, commissioning, and research. Recommendations are also provided for developers and vendors to help users to understand DIR uncertainties and make the application of DIR in RT safer and more reliable.

Technical note: Evaluation of deep learning based synthetic CTs clinical readiness for dose and NTCP driven head and neck adaptive proton therapy

BACKGROUND

Adaptive proton therapy workflows rely on accurate imaging throughout the treatment course. Our centre currently utilizes weekly repeat CTs (rCTs) for treatment monitoring and plan adaptations. However, deep learning-based methods have recently shown to successfully correct CBCT images, which suffer from severe imaging artifacts, and generate high quality synthetic CT (sCT) images which enable CBCT-based proton dose calculations.

PURPOSE

To compare daily CBCT-based sCT images to planning CTs (pCT) and rCTs of head and neck (HN) cancer patients to investigate the dosimetric accuracy of CBCT-based sCTs in a scenario mimicking actual clinical practice.

METHODS

Data of 56 HN cancer patients, previously treated with proton therapy was used to generate 1.962 sCT images, using a previously developed and trained deep convolutional neural network. Clinical IMPT treatment plans were recalculated on the pCT, weekly rCTs and daily sCTs. The dosimetric accuracy of sCTs was compared to same day rCTs and the initial planning CT. As a reference, rCTs were also compared to pCTs. The dose difference between sCTs and rCTs/pCT was quantified by calculating the D98 difference for target volumes and Dmean difference for organs-at-risk. To investigate the clinical relevancy of possible dose differences, NTCP values were calculated for dysphagia and xerostomia.

RESULTS

For target volumes, only minor dose differences were found for sCT versus rCT and sCT versus pCT, with dose differences mostly within ±1.5%. Larger dose differences were observed in OARs, where a general shift towards positive differences was found, with the largest difference in the left parotid gland. Delta NTCP values for grade 2 dysphagia and xerostomia were within ±2.5% for 90% of the sCTs.

CONCLUSIONS

Target doses showed high similarity between rCTs and sCTs. Further investigations are required to identify the origin of the dose differences at OAR levels and its relevance in clinical decision making.

Deep learning-based fast denoising of Monte Carlo dose calculation in carbon ion radiotherapy

BACKGROUND

Plan verification is one of the important steps of quality assurance (QA) in carbon ion radiotherapy. Conventional methods of plan verification are based on phantom measurement, which is labor-intensive and time-consuming. Although the plan verification method based on Monte Carlo (MC) simulation provides a more accurate modeling of the physics, it is also time-consuming when simulating with a large number of particles. Therefore, how to ensure the accuracy of simulation results while reducing simulation time is the current difficulty and focus.

PURPOSE

The purpose of this work was to evaluate the feasibility of using deep learning-based MC denoising method to accelerate carbon-ion radiotherapy plan verification.

METHODS

Three models, including CycleGAN, 3DUNet and GhostUNet with Ghost module, were used to denoise the 1 × 106 carbon ions-based MC dose distribution to the accuracy of 1 × 108 carbon ions-based dose distribution. The CycleGAN's generator, 3DUNet and GhostUNet were all derived from the 3DUNet network. A total of 59 cases including 29 patients with head-and-neck cancers and 30 patients with lung cancers were collected, and 48 cases were randomly selected as the training set of the CycleGAN network and six cases as the test set. For the 3DUNet and GhostUNet models, the numbers of training set, validation set, and test set were 47, 6, and 6, respectively. Finally, the three models were evaluated qualitatively and quantitatively using RMSE and three-dimensional gamma analysis (3 mm, 3%).

RESULTS

The three end-to-end trained models could be used for denoising the 1 × 106 carbon ions-based dose distribution, and their generalization was proved. The GhostUNet obtained the lowest RMSE value of 0.075, indicating the smallest difference between its denoised and 1 × 108 carbon ions-based dose distributions. The average gamma passing rate (GPR) between the GhostUNet denoising-based versus 1 × 108 carbon ions-based dose distributions was 99.1%, higher than that of the CycleGAN at 94.3% and the 3DUNet at 96.2%. Among the three models, the GhostUNet model had the fewest parameters (4.27 million) and the shortest training time (99 s per epoch) but achieved the best denoising results.

CONCLUSION

The end-to-end deep network GhostUNet outperforms the CycleGAN, 3DUNet models in denoising MC dose distributions for carbon ion radiotherapy. The network requires less than 5 s to denoise a sample of MC simulation with few particles to obtain a qualitative and quantitative result comparable to the dose distribution simulated by MC with relatively large number particles, offering a significant reduction in computation time.

Deep learning framework to improve the quality of cone-beam computed tomography for radiotherapy scenarios

BACKGROUND

The application of cone-beam computed tomography (CBCT) in image-guided radiotherapy and adaptive radiotherapy remains limited due to its poor image quality.

PURPOSE

In this study, we aim to develop a deep learning framework to generate high-quality CBCT images for therapeutic applications.

METHODS

The synthetic CT (sCT) generation from the CBCT was proposed using a transformer-based network with a hybrid loss function. The network was trained and validated using the data from 176 patients to produce a general model that can be extensively applied to enhance CBCT images. After the first therapy, each patient can receive paired CBCT/planning CT (pCT) scans, and the obtained data were used to fine-tune the general model for further improvement. For subsequent treatment, a patient-specific, personalized model was made available. In total, 34 patients were examined for general model testing, and another six patients who underwent rescanned pCT scan were used for personalized model training and testing.

RESULTS

The general model decreased the mean absolute error (MAE) from 135 HU to 59 HU as compared to the CBCT. The hybrid loss function demonstrated superior performance in CT number correction and noise/artifacts reduction. The proposed transformer-based network also showed superior power in CT number correction compared to the classical convolutional neural network. The personalized model showed improvement based on the general model in some details, and the MAE was reduced from 59 HU (for the general model) to 57 HU (p < 0.05 Wilcoxon signed-rank test).

CONCLUSION

We established a deep learning framework based on transformer for clinical needs. The deep learning model demonstrated potential for continuous improvement with the help of a suggested personalized training strategy compatible with the clinical workflow.

Detectability of Anatomical Changes With Prompt-Gamma Imaging: First Systematic Evaluation of Clinical Application During Prostate-Cancer Proton Therapy

PURPOSE

The development of online-adaptive proton therapy (PT) is essential to overcome limitations encountered by day-to-day variations of the patient's anatomy. Range verification could play an essential role in an online feedback loop for the detection of treatment deviations such as anatomical changes. Here, we present the results of the first systematic patient study regarding the detectability of anatomical changes by a prompt-gamma imaging (PGI) slit-camera system.

METHODS AND MATERIALS

For 15 patients with prostate cancer, PGI measurements were performed during 105 fractions (201 fields) with in-room control computed tomography (CT)acquisitions. Field-wise doses on control CT scans were manually classified as whether showing relevant or non-relevant anatomical changes. This manual classification of the treatment fields was then used to establish an automatic field-wise ground truth based on spot-wise dosimetric range shifts, which were retrieved from integrated depth-dose (IDD) profiles. To determine the detection capability of anatomical changes with PGI, spot-wise PGI-based range shifts were initially compared with corresponding dosimetric IDD range shifts. As final endpoint, the agreement of a developed field-wise PGI classification model with the field-wise ground truth was determined. Therefore, the PGI model was optimized and tested for a subcohort of 131 and 70 treatment fields, respectively.

RESULTS

The correlation between PGI and IDD range shifts was high, ρpearson = 0.67 (p < 0.01). Field-wise binary PGI classification resulted in an area under the curve of 0.72 and 0.80 for training and test cohorts, respectively. The model detected relevant anatomical changes in the independent test cohort, with a sensitivity and specificity of 74% and 79%, respectively.

CONCLUSIONS

For the first time, evidence of the detection capability of anatomical changes in prostate-cancer PT from clinically acquired PGI data is shown. This emphasizes the benefit of PGI-based range verification and demonstrates its potential for online-adaptive PT.

Comparison and evaluation of different deep learning models of synthetic CT generation from CBCT for nasopharynx cancer adaptive proton therapy

BACKGROUND

Cone-beam computed tomography (CBCT) scanning is used for patient setup in image-guided radiotherapy. However, its inaccurate CT numbers limit its applicability in dose calculation and treatment planning.

PURPOSE

This study compares four deep learning methods for generating synthetic CT (sCT) to determine which method is more appropriate and offers potential for further clinical exploration in adaptive proton therapy for nasopharynx cancer.

METHODS

CBCTs and deformed planning CT (dCT) from 75 patients (60/5/10 for training, validation and testing) were used to compare cycle-consistent Generative Adversarial Network (cycleGAN), Unet, Unet+cycleGAN and conditionalGenerative Adversarial Network (cGAN) for sCT generation. The sCT images generated by each method were evaluated against dCT images using mean absolute error (MAE), structural similarity (SSIM), peak signal-to-noise ratio (PSNR), spatial non-uniformity (SNU) and radial averaging in the frequency domain. In addition, dosimetric accuracy was assessed through gamma analysis, differences in water equivalent thickness (WET), and dose-volume histogram metrics.

RESULTS

The cGAN model has demonstrated optimal performance in the four models across various indicators. In terms of image quality under global condition, the average MAE has been reduced to 16.39HU, SSIM has increased to 95.24%, and PSNR has increased to 28.98. Regarding dosimetric accuracy, the gamma passing rate (2%/2 mm) has reached 99.02%, and the WET difference is only 1.28 mm. The D95 value of CTVs coverage and Dmax value of spinal cord, brainstem show no significant differences between dCT and sCT generated by cGAN model.

CONCLUSIONS

The cGAN model has been shown to be a more suitable approach for generating sCT using CBCT, considering its characteristics and concepts. The resulting sCT has the potential for application in adaptive proton therapy.

Comparing Predicted Toxicities between Hypofractionated Proton and Photon Radiotherapy of Liver Cancer Patients with Different Adaptive Schemes

With the availability of MRI linacs, online adaptive intensity modulated radiotherapy (IMRT) has become a treatment option for liver cancer patients, often combined with hypofractionation. Intensity modulated proton therapy (IMPT) has the potential to reduce the dose to healthy tissue, but it is particularly sensitive to changes in the beam path and might therefore benefit from online adaptation. This study compares the normal tissue complication probabilities (NTCPs) for liver and duodenal toxicity for adaptive and non-adaptive IMRT and IMPT treatments of liver cancer patients. Adaptive and non-adaptive IMRT and IMPT plans were optimized to 50 Gy (RBE = 1.1 for IMPT) in five fractions for 10 liver cancer patients, using the original MRI linac images and physician-drawn structures. Three liver NTCP models were used to predict radiation-induced liver disease, an increase in albumin-bilirubin level, and a Child-Pugh score increase of more than 2. Additionally, three duodenal NTCP models were used to predict gastric bleeding, gastrointestinal (GI) toxicity with grades >3, and duodenal toxicity grades 2-4. NTCPs were calculated for adaptive and non-adaptive IMRT and IMPT treatments. In general, IMRT showed higher NTCP values than IMPT and the differences were often significant. However, the differences between adaptive and non-adaptive treatment schemes were not significant, indicating that the NTCP benefit of adaptive treatment regimens is expected to be smaller than the expected difference between IMRT and IMPT.

Dosimetric comparison of autocontouring techniques for online adaptive proton therapy

Objective.Anatomical and daily set-up uncertainties impede high precision delivery of proton therapy. With online adaptation, the daily plan is reoptimized on an image taken shortly before the treatment, reducing these uncertainties and, hence, allowing a more accurate delivery. This reoptimization requires target and organs-at-risk (OAR) contours on the daily image, which need to be delineated automatically since manual contouring is too slow. Whereas multiple methods for autocontouring exist, none of them are fully accurate, which affects the daily dose. This work aims to quantify the magnitude of this dosimetric effect for four contouring techniques.Approach.Plans reoptimized on automatic contours are compared with plans reoptimized on manual contours. The methods include rigid and deformable registration (DIR), deep-learning based segmentation and patient-specific segmentation.Main results.It was found that independently of the contouring method, the dosimetric influence of usingautomaticOARcontoursis small (<5% prescribed dose in most cases), with DIR yielding the best results. Contrarily, the dosimetric effect of using theautomatic target contourwas larger (>5% prescribed dose in most cases), indicating that manual verification of that contour remains necessary. However, when compared to non-adaptive therapy, the dose differences caused by automatically contouring the target were small and target coverage was improved, especially for DIR.Significance.The results show that manual adjustment of OARs is rarely necessary and that several autocontouring techniques are directly usable. Contrarily, manual adjustment of the target is important. This allows prioritizing tasks during time-critical online adaptive proton therapy and therefore supports its further clinical implementation.

Clinical evaluation of synthetic computed tomography methods in adaptive proton therapy of lung cancer patients

Background and purpose

Efficient workflows for adaptive proton therapy are of high importance. This study evaluated the possibility to replace repeat-CTs (reCTs) with synthetic CTs (sCTs), created based on cone-beam CTs (CBCTs), for flagging the need of plan adaptations in intensity-modulated proton therapy (IMPT) treatment of lung cancer patients.

Materials and methods

Forty-two IMPT patients were retrospectively included. For each patient, one CBCT and a same-day reCT were included. Two commercial sCT methods were applied; one based on CBCT number correction (Cor-sCT), and one based on deformable image registration (DIR-sCT). The clinical reCT workflow (deformable contour propagation and robust dose re-computation) was performed on the reCT as well as the two sCTs. The deformed target contours on the reCT/sCTs were checked by radiation oncologists and edited if needed. A dose-volume-histogram triggered plan adaptation method was compared between the reCT and the sCTs; patients needing a plan adaptation on the reCT but not on the sCT were denoted false negatives. As secondary evaluation, dose-volume-histogram comparison and gamma analysis (2%/2mm) were performed between the reCT and sCTs.

Results

There were five false negatives, two for Cor-sCT and three for DIR-sCT. However, three of these were only minor, and one was caused by tumour position differences between the reCT and CBCT and not by sCT quality issues. An average gamma pass rate of 93% was obtained for both sCT methods.

Conclusion

Both sCT methods were judged to be of clinical quality and valuable for reducing the amount of reCT acquisitions.

Image guidance and interfractional anatomical variation in paediatric abdominal radiotherapy

OBJECTIVES

To identify variables predicting interfractional anatomical variations measured with cone-beam CT (CBCT) throughout abdominal paediatric radiotherapy, and to assess the potential of surface-guided radiotherapy (SGRT) to monitor these changes.

METHODS

Metrics of variation in gastrointestinal (GI) gas volume and separation of the body contour and abdominal wall were calculated from 21 planning CTs and 77 weekly CBCTs for 21 abdominal neuroblastoma patients (median 4 years, range: 2 - 19 years). Age, sex, feeding tubes, and general anaesthesia (GA) were explored as predictive variables for anatomical variation. Furthermore, GI gas variation was correlated with changes in body and abdominal wall separation, as well as simulated SGRT metrics of translational and rotational corrections between CT/CBCT.

RESULTS

GI gas volumes varied 74 ± 54 ml across all scans, while body and abdominal wall separation varied 2.0 ± 0.7 mm and 4.1 ± 1.5 mm from planning, respectively. Patients < 3.5 years (p = 0.04) and treated under GA (p < 0.01) experienced greater GI gas variation; GA was the strongest predictor in multivariate analysis (p < 0.01). Absence of feeding tubes was linked to greater body contour variation (p = 0.03). GI gas variation correlated with body (R = 0.53) and abdominal wall (R = 0.63) changes. The strongest correlations with SGRT metrics were found for anterior-posterior translation (R = 0.65) and rotation of the left-right axis (R = -0.36).

CONCLUSIONS

Young age, GA, and absence of feeding tubes were linked to stronger interfractional anatomical variation and are likely indicative of patients benefiting from adaptive/robust planning pathways. Our data suggest a role for SGRT to inform the need for CBCT at each treatment fraction in this patient group.

ADVANCES IN KNOWLEDGE

This is the first study to suggest the potential role of SGRT for the management of internal interfractional anatomical variation in paediatric abdominal radiotherapy.

Online adaptive planning methods for intensity-modulated radiotherapy

Online adaptive radiation therapy aims at adapting a patient's treatment plan to their current anatomy to account for inter-fraction variations before daily treatment delivery. As this process needs to be accomplished while the patient is immobilized on the treatment couch, it requires time-efficient adaptive planning methods to generate a quality daily treatment plan rapidly. The conventional planning methods do not meet the time requirement of online adaptive radiation therapy because they often involve excessive human intervention, significantly prolonging the planning phase. This article reviews the planning strategies employed by current commercial online adaptive radiation therapy systems, research on online adaptive planning, and artificial intelligence's potential application to online adaptive planning.

A hybrid multi-particle approach to range assessment-based treatment verification in particle therapy

Particle therapy (PT) used for cancer treatment can spare healthy tissue and reduce treatment toxicity. However, full exploitation of the dosimetric advantages of PT is not yet possible due to range uncertainties, warranting development of range-monitoring techniques. This study proposes a novel range-monitoring technique introducing the yet unexplored concept of simultaneous detection and imaging of fast neutrons and prompt-gamma rays produced in beam-tissue interactions. A quasi-monolithic organic detector array is proposed, and its feasibility for detecting range shifts in the context of proton therapy is explored through Monte Carlo simulations of realistic patient models and detector resolution effects. The results indicate that range shifts of [Formula: see text] can be detected at relatively low proton intensities ([Formula: see text] protons/spot) when spatial information obtained through imaging of both particle species are used simultaneously. This study lays the foundation for multi-particle detection and imaging systems in the context of range verification in PT.

A survey of practice patterns for adaptive particle therapy for interfractional changes

Background and purpose

Anatomical changes may compromise the planned target coverage and organs-at-risk dose in particle therapy. This study reports on the practice patterns for adaptive particle therapy (APT) to evaluate current clinical practice and wishes and barriers to further implementation.

Materials and methods

An institutional questionnaire was distributed to PT centres worldwide (7/2020-6/2021) asking which type of APT was used, details of the workflow, and what the wishes and barriers to implementation were. Seventy centres from 17 countries participated. A three-round Delphi consensus analysis (10/2022) among the authors followed to define recommendations on required actions and future vision.

Results

Out of the 68 clinically operational centres, 84% were users of APT for at least one treatment site with head and neck being most common. APT was mostly performed offline with only two online APT users (plan-library). No centre used online daily re-planning. Daily 3D imaging was used for APT by 19% of users. Sixty-eight percent of users had plans to increase their use or change their technique for APT. The main barrier was "lack of integrated and efficient workflows". Automation and speed, reliable dose deformation for dose accumulation and higher quality of in-room volumetric imaging were identified as the most urgent task for clinical implementation of online daily APT.

Conclusion

Offline APT was implemented by the majority of PT centres. Joint efforts between industry research and clinics are needed to translate innovations into efficient and clinically feasible workflows for broad-scale implementation of online APT.

Large anatomical changes in head-and-neck cancers – a dosimetric comparison of online and offline adaptive proton therapy

Purpose

This work evaluates an online adaptive (OA) workflow for head-and-neck (H&N) intensity-modulated proton therapy (IMPT) and compares it with full offline replanning (FOR) in patients with large anatomical changes.

Methods

IMPT treatment plans are created retrospectively for a cohort of eight H&N cancer patients that previously required replanning during the course of treatment due to large anatomical changes. Daily cone-beam CTs (CBCT) are acquired and corrected for scatter, resulting in 253 analyzed fractions. To simulate the FOR workflow, nominal plans are created on the planning-CT and delivered until a repeated-CT is acquired; at this point, a new plan is created on the repeated-CT. To simulate the OA workflow, nominal plans are created on the planning-CT and adapted at each fraction using a simple beamlet weight-tuning technique. Dose distributions are calculated on the CBCTs with Monte Carlo for both delivery methods. The total treatment dose is accumulated on the planning-CT.

Results

Daily OA improved target coverage compared to FOR despite using smaller target margins. In the high-risk CTV, the median D₉₈ degradation was 1.1 % and 2.1 % for OA and FOR, respectively. In the low-risk CTV, the same metrics yield 1.3 % and 5.2 % for OA and FOR, respectively. Smaller setup margins of OA reduced the dose to all OARs, which was most relevant for the parotid glands.

Conclusion

Daily OA can maintain prescription doses and constraints over the course of fractionated treatment, even in cases of large anatomical changes, reducing the necessity for manual replanning in H&N IMPT.

Dosimetric impact of adaptive proton therapy in head and neck cancer - A review

Background

Intensity Modulated Proton Therapy (IMPT) in head and neck cancer (HNC) is susceptible to anatomical changes and patient set-up inaccuracies during the radiotherapy course, which can cause discrepancies between planned and delivered dose. The discrepancies can be counteracted by adaptive replanning strategies. This article reviews the observed dosimetric impact of adaptive proton therapy (APT) and the timing to perform a plan adaptation in IMPT in HNC.

Methods

A literature search of articles published in PubMed/MEDLINE, EMBASE and Web of Science from January 2010 to March 2022 was performed. Among a total of 59 records assessed for possible eligibility, ten articles were included in this review.

Results

Included studies reported on target coverage deterioration in IMPT plans during the RT course, which was recovered with the application of an APT approach. All APT plans showed an average improved target coverage for the high- and low-dose targets as compared to the accumulated dose on the planned plans. Dose improvements up to 2.5 Gy (3.5 %) and up to 4.0 Gy (7.1 %) in the D98 of the high- and low dose targets were observed with APT. Doses to the organs at risk (OARs) remained equal or decreased slightly after APT was applied. In the included studies, APT was largely performed once, which resulted in the largest target coverage improvement, but eventual additional APT improved the target coverage further. There is no data showing what is the most appropriate timing for APT.

Conclusion

APT during IMPT for HNC patients improves target coverage. The largest improvement in target coverage was found with a single adaptive intervention, and an eventual second or more frequent APT application improved the target coverage further. Doses to the OARs remained equal or decreased slightly after applying APT. The most optimal timing for APT is yet to be determined.

Proton Beam Therapy and Photon-Based Magnetic Resonance Image-Guided Radiation Therapy: The Next Frontiers of Radiation Therapy for Hepatocellular Carcinoma

External beam radiation therapy (EBRT) has increasingly been utilized in the treatment of hepatocellular carcinoma (HCC) due to technological advances with positive clinical outcomes. Innovations in EBRT include improved image guidance, motion management, treatment planning, and highly conformal techniques such as intensity-modulated radiation therapy (IMRT) and stereotactic body radiation therapy (SBRT). Moreover, proton beam therapy (PBT) and magnetic resonance image-guided radiation therapy (MRgRT) have expanded the capabilities of EBRT. PBT offers the advantage of minimizing low- and moderate-dose radiation to the surrounding normal tissue, thereby preserving uninvolved liver and allowing for dose escalation. MRgRT provides the advantage of improved soft tissue delineation compared to computerized tomography (CT) guidance. Additionally, MRgRT with online adaptive therapy is particularly useful for addressing motion not otherwise managed and reducing high-dose radiation to the normal tissue such as the stomach and bowel. PBT and online adaptive MRgRT are emerging technological advancements in EBRT that may provide a significant clinical benefit for patients with HCC.

GPU accelerated Monte Carlo scoring of positron emitting isotopes produced during proton therapy for PET verification

Objective.Verification of delivered proton therapy treatments is essential for reaping the many benefits of the modality, with the most widely proposedin vivoverification technique being the imaging of positron emitting isotopes generated in the patient during treatment using positron emission tomography (PET). The purpose of this work is to reduce the computational resources and time required for simulation of patient activation during proton therapy using the GPU accelerated Monte Carlo code FRED, and to validate the predicted activity against the widely used Monte Carlo code GATE.Approach.We implement a continuous scoring approach for the production of positron emitting isotopes within FRED version 5.59.9. We simulate treatment plans delivered to 95 head and neck patients at Centrum Cyklotronowe Bronowice using this GPU implementation, and verify the accuracy using the Monte Carlo toolkit GATE version 9.0.Main results.We report an average reduction in computational time by a factor of 50 when using a local system with 2 GPUs as opposed to a large compute cluster utilising between 200 to 700 CPU threads, enabling simulation of patient activity within an average of 2.9 min as opposed to 146 min. All simulated plans are in good agreement across the two Monte Carlo codes. The two codes agree within a maximum of 0.95σon a voxel-by-voxel basis for the prediction of 7 different isotopes across 472 simulated fields delivered to 95 patients, with the average deviation over all fields being 6.4 × 10^-3σ.Significance.The implementation of activation calculations in the GPU accelerated Monte Carlo code FRED provides fast and reliable simulation of patient activation following proton therapy, allowing for research and development of clinical applications of range verification for this treatment modality using PET to proceed at a rapid pace.

SWFT-Net: a deep learning framework for efficient fine-tuning spot weights towards adaptive proton therapy

Objective. One critical task for adaptive proton therapy is how to perform spot weight re-tuning and reoptimize plan, both of which are time-consuming and labor intensive. We proposed a deep learning framework (SWFT-Net) to speed up such a task, a starting point for us to move towards online adaptive proton therapy.Approach. For a H&N patient case, a reference intensity modulated proton therapy plan was generated. For data augmentation, spot weights were modified to generate three datasets (DS10, DS30, DS50), corresponding to different levels of weight adjustment. For each dataset, the samples were split into the training and testing groups at a ratio of 8:2 (6400 for training, 1706 for testing). To ease the difficulty of machine learning, the residuals of dose maps and spot weights (i.e. difference relative to a reference) were used as inputs and outputs, respectively. Quantitative analyses were performed in terms of normalized root mean square error (NRMSE) of spot weights, Gamma passing rate and dose difference within the PTV.Main results. The SWFT-Net is able to generate an adapted plan in less than a second with a NVIDIA GeForce RTX 3090 GPU. For the 1706 samples in the testing dataset, the NRMSE is 0.41% (DS10), 1.05% (DS30) and 2.04% (DS50), respectively. Cold/hot spots in the dose maps after adaptation are observed. The mean relative dose difference is 0.64% (DS10), 0.92% (DS30) and 0.88% (DS50), respectively. For all three datasets, the mean Gamma passing rate is consistently over 95% for both 1 mm/1% and 3 mm/3% settings.Significance. The proposed SWFT-Net is a promising tool to help realize adaptive proton therapy. It can be used as an alternative tool to other spot fine-tuning optimization algorithms, likely demonstrating superior performance in terms of speed, accuracy, robustness and minimum human interaction. This study lays down a foundation for us to move further incorporating other factors such as daily anatomical changes and propagated PTVs, and develop a truly online adaptive workflow in proton therapy.

An online adaptive plan library approach for intensity modulated proton therapy for head and neck cancer

BACKGROUND AND PURPOSE

In intensity modulated proton therapy (IMPT), the impact of setup errors and anatomical changes is commonly mitigated by robust optimization with population-based setup robustness (SR) settings and offline replanning. In this study we propose and evaluate an alternative approach based on daily plan selection from patient-specific pre-treatment established plan libraries (PLs). Clinical implementation of the PL strategy would be rather straightforward compared to daily online re-planning.

MATERIALS AND METHODS

For 15 head-and-neck cancer patients, the planning CT was used to generate a PL with 5 plans, robustly optimized for increasing SR: 0, 1, 2, 3, 5 mm, and 3% range robustness. Repeat CTs (rCTs) and realistic setup and range uncertainty distributions were used for simulation of treatment courses for the PL approach, treatments with fixed SR (fSR₃) and a trigger-based offline adaptive schedule for 3 mm SR (fSR₃OfA). Daily plan selection in the PL approach was based only on recomputed dose to the CTV on the rCT.

RESULTS

Compared to using fSR₃ and fSR₃OfA, the risk of xerostomia grade ≥ II & III and dysphagia ≥ grade III were significantly reduced with the PL. For 6/15 patients the risk of xerostomia and/or dysphagia ≥ grade II could be reduced by > 2% by using PL. For the other patients, adherence to target coverage constraints was often improved. fSR₃OfA resulted in significantly improved coverage compared to PL for selected patients.

CONCLUSION

The proposed PL approach resulted in overall reduced NTCPs compared to fSR₃ and fSR₃OfA at limited cost in target coverage.

Use of CBCT plus plan robustness for reducing QACT frequency in intensity-modulated proton therapy: Head-and-neck cases

PURPOSE

Anatomic variation has a significant dosimetric impact in intensity-modulated proton therapy. Weekly or biweekly computed tomography (CT) scans, called quality assurance CTs (QACTs), are used to monitor anatomic and resultant dose changes to determine whether adaptive plans are needed. Frequent CT scans result in unwanted QACT dose and increased clinical workloads. This study proposed utilizing patient setup cone-beam CTs (CBCTs) and treatment plan robustness to reduce the frequency of QACTs.

METHODS

We retrospectively analyzed data from 27 patients with head-and-neck cancer, including 594 CBCTs, 136 QACTs, and 19 adaptive plans. For each CBCT, water-equivalent thickness (WET) along the pencil-beam path was calculated. For each treatment plan, the WET of the first-day CBCT was used as the reference, and the mean WET changes (ΔWET) in each following CBCT was used as the surrogate of proton range change. Using CBCTs acquired prior to a QACT, we predicted the ΔWET on the QACT day by a linear regression model. The impact of range change on target dose was calculated as the predicted ΔWET weighted by the monitor units of each field. In addition, plan robustness was estimated from the robust dose-volume histograms (DVHs) and combined with ΔWET to reduce QACT frequency. Robustness was estimated from the distance between the DVH curves of the nominal and worst scenarios.

RESULTS

When the estimated mean ΔWET was <6.5 mm (or <7.5 mm if the robustness was >95%), the QACT could be skipped without missing any adaptive planning; otherwise a QACT was required. Overall, 41% of QACTs could be eliminated when ΔWET was <6.5 mm and 56% when ΔWET was <7.5 mm, and robustness was >95%. At least one QACT could have been omitted in 25 of the 27 cases under skipping thresholds at ΔWETs <7.5 mm and R > 95%.

CONCLUSION

This study suggests that the number of QACTs can be greatly reduced by calculating range change in patient setup CBCTs and can be further reduced by combining this information with analyses of plan robustness.

Deep learning-based 4D-synthetic CTs from sparse-view CBCTs for dose calculations in adaptive proton therapy

BACKGROUND

Time-resolved 4D cone beam-computed tomography (4D-CBCT) allows a daily assessment of patient anatomy and respiratory motion. However, 4D-CBCTs suffer from imaging artifacts that affect the CT number accuracy and prevent accurate proton dose calculations. Deep learning can be used to correct CT numbers and generate synthetic CTs (sCTs) that can enable CBCT-based proton dose calculations.

PURPOSE

In this work, sparse view 4D-CBCTs were converted into 4D-sCT utilizing a deep convolutional neural network (DCNN). 4D-sCTs were evaluated in terms of image quality and dosimetric accuracy to determine if accurate proton dose calculations for adaptive proton therapy workflows of lung cancer patients are feasible.

METHODS

A dataset of 45 thoracic cancer patients was utilized to train and evaluate a DCNN to generate 4D-sCTs, based on sparse view 4D-CBCTs reconstructed from projections acquired with a 3D acquisition protocol. Mean absolute error (MAE) and mean error were used as metrics to evaluate the image quality of single phases and average 4D-sCTs against 4D-CTs acquired on the same day. The dosimetric accuracy was checked globally (gamma analysis) and locally for target volumes and organs-at-risk (OARs) (lung, heart, and esophagus). Furthermore, 4D-sCTs were also compared to 3D-sCTs. To evaluate CT number accuracy, proton radiography simulations in 4D-sCT and 4D-CTs were compared in terms of range errors. The clinical suitability of 4D-sCTs was demonstrated by performing a 4D dose reconstruction using patient specific treatment delivery log files and breathing signals.

RESULTS

4D-sCTs resulted in average MAEs of 48.1 ± 6.5 HU (single phase) and 37.7 ± 6.2 HU (average). The global dosimetric evaluation showed gamma pass ratios of 92.3% ± 3.2% (single phase) and 94.4% ± 2.1% (average). The clinical target volume showed high agreement in D₉₈ between 4D-CT and 4D-sCT, with differences below 2.4% for all patients. Larger dose differences were observed in mean doses of OARs (up to 8.4%). The comparison with 3D-sCTs showed no substantial image quality and dosimetric differences for the 4D-sCT average. Individual 4D-sCT phases showed slightly lower dosimetric accuracy. The range error evaluation revealed that lung tissues cause range errors about three times higher than the other tissues.

CONCLUSION

In this study, we have investigated the accuracy of deep learning-based 4D-sCTs for daily dose calculations in adaptive proton therapy. Despite image quality differences between 4D-sCTs and 3D-sCTs, comparable dosimetric accuracy was observed globally and locally. Further improvement of 3D and 4D lung sCTs could be achieved by increasing CT number accuracy in lung tissues.

Low-dose CT allows for accurate proton therapy dose calculation and plan optimization

Objective. Protons offer a more conformal dose delivery compared to photons, yet they are sensitive to anatomical changes over the course of treatment. To minimize range uncertainties due to anatomical variations, a new CT acquisition at every treatment session would be paramount to enable daily dose calculation and subsequent plan adaptation. However, the series of CT scans results in an additional accumulated patient dose. Reducing CT radiation dose and thereby decreasing the potential risk of radiation exposure to patients is desirable, however, lowering the CT dose results in a lower signal-to-noise ratio and therefore in a reduced quality image. We hypothesized that the signal-to-noise ratio provided by conventional CT protocols is higher than needed for proton dose distribution estimation. In this study, we aim to investigate the effect of CT imaging dose reduction on proton therapy dose calculations and plan optimization. Approach. To verify our hypothesis, a CT dose reduction simulation tool has been developed and validated to simulate lower-dose CT scans from an existing standard-dose scan. The simulated lower-dose CTs were then used for proton dose calculation and plan optimization and the results were compared with those of the standard-dose scan. The same strategy was adopted to investigate the effect of CT dose reduction on water equivalent thickness (WET) calculation to quantify CT noise accumulation during integration along the beam. Main results. The similarity between the dose distributions acquired from the low-dose and standard-dose CTs was evaluated by the dose-volume histogram and the 3D Gamma analysis. The results on an anthropomorphic head phantom and three patient cases indicate that CT imaging dose reduction up to 90% does not have a significant effect on proton dose calculation and plan optimization. The relative error was employed to evaluate the similarity between WET maps and was found to be less than 1% after reducing the CT imaging dose by 90%. Significance. The results suggest the possibility of using low-dose CT for proton therapy dose estimation, since the dose distributions acquired from the standard-dose and low-dose CTs are clinically equivalent.

In-vivo range verification analysis with in-beam PET data for patients treated with proton therapy at CNAO

Morphological changes that may arise through a treatment course are probably one of the most significant sources of range uncertainty in proton therapy. Non-invasive in-vivo treatment monitoring is useful to increase treatment quality. The INSIDE in-beam Positron Emission Tomography (PET) scanner performs in-vivo range monitoring in proton and carbon therapy treatments at the National Center of Oncological Hadrontherapy (CNAO). It is currently in a clinical trial (ID: NCT03662373) and has acquired in-beam PET data during the treatment of various patients. In this work we analyze the in-beam PET (IB-PET) data of eight patients treated with proton therapy at CNAO. The goal of the analysis is twofold. First, we assess the level of experimental fluctuations in inter-fractional range differences (sensitivity) of the INSIDE PET system by studying patients without morphological changes. Second, we use the obtained results to see whether we can observe anomalously large range variations in patients where morphological changes have occurred. The sensitivity of the INSIDE IB-PET scanner was quantified as the standard deviation of the range difference distributions observed for six patients that did not show morphological changes. Inter-fractional range variations with respect to a reference distribution were estimated using the Most-Likely-Shift (MLS) method. To establish the efficacy of this method, we made a comparison with the Beam’s Eye View (BEV) method. For patients showing no morphological changes in the control CT the average range variation standard deviation was found to be 2.5 mm with the MLS method and 2.3 mm with the BEV method. On the other hand, for patients where some small anatomical changes occurred, we found larger standard deviation values. In these patients we evaluated where anomalous range differences were found and compared them with the CT. We found that the identified regions were mostly in agreement with the morphological changes seen in the CT scan.

Robust treatment planning in scanned carbon-ion radiotherapy for pancreatic cancer: Clinical verification using in-room computed tomography images

Purpose

Carbon-ion beam (C-beam) has a sharp dose distribution called the Bragg peak. Carbon-ion radiation therapy, such as stereotactic body radiotherapy in photon radiotherapy, can be completed in a short period by concentrating the radiation dose on the tumor while minimizing the dose to organs at-risk. However, the stopping position of C-beam is sensitive to density variations along the beam path and such variations can lower the tumor dose as well as cause the delivery of an unexpectedly high dose to the organs at risk. We evaluated the clinical efficacy of a robust planning technique considering gastrointestinal gas (G-gas) to deliver accurate radiation doses in carbon-ion radiotherapy for pancreatic cancer.

Materials and methods

We focused on the computed tomography (CT) value replacement method. Replacement signifies the overwriting of CT values in the CT images. The most effective replacement method for robust treatment planning was determined by verifying the effects of the three replacement patterns. We selected 10 consecutive patients. Pattern 1 replaces the CT value of the G-gas contours with the value of the region without G-gas (P1). This condition indicates a no-gas state. Pattern 2 replaces each gastrointestinal contour using the mean CT value of each contour (P2). The effect of G-gas was included in the replacement value. Pattern 3 indicates no replacement (P3). We analyzed variations in the target coverage (TC) and homogeneity index (HI) from the initial plan using in-room CT images. We then performed correlation analysis on the variations in G-gas, TC, and HI to evaluate the robustness against G-gas.

Results

Analysis of variations in TC and HI revealed a significant difference between P1 and P3 and between P2 and P3. Although no statistically significant difference was observed between P1 and P2, variations, including the median, tended to be fewer in P2. The correlation analyses for G-gas, TC, and HI showed that P2 was less likely to be affected by G-gas.

Conclusion

For a treatment plan that is robust to G-gas, P2 mean replacement method should be used. This method does not necessitate any particular software or equipment, and is convenient to implement in clinical practice.

A more effective CT synthesizer using transformers for cone-beam CT-guided adaptive radiotherapy

Purpose

The challenge of cone-beam computed tomography (CBCT) is its low image quality, which limits its application for adaptive radiotherapy (ART). Despite recent substantial improvement in CBCT imaging using the deep learning method, the image quality still needs to be improved for effective ART application. Spurred by the advantages of transformers, which employs multi-head attention mechanisms to capture long-range contextual relations between image pixels, we proposed a novel transformer-based network (called TransCBCT) to generate synthetic CT (sCT) from CBCT. This study aimed to further improve the accuracy and efficiency of ART.

Materials and methods

In this study, 91 patients diagnosed with prostate cancer were enrolled. We constructed a transformer-based hierarchical encoder–decoder structure with skip connection, called TransCBCT. The network also employed several convolutional layers to capture local context. The proposed TransCBCT was trained and validated on 6,144 paired CBCT/deformed CT images from 76 patients and tested on 1,026 paired images from 15 patients. The performance of the proposed TransCBCT was compared with a widely recognized style transferring deep learning method, the cycle-consistent adversarial network (CycleGAN). We evaluated the image quality and clinical value (application in auto-segmentation and dose calculation) for ART need.

Results

TransCBCT had superior performance in generating sCT from CBCT. The mean absolute error of TransCBCT was 28.8 ± 16.7 HU, compared to 66.5 ± 13.2 for raw CBCT, and 34.3 ± 17.3 for CycleGAN. It can preserve the structure of raw CBCT and reduce artifacts. When applied in auto-segmentation, the Dice similarity coefficients of bladder and rectum between auto-segmentation and oncologist manual contours were 0.92 and 0.84 for TransCBCT, respectively, compared to 0.90 and 0.83 for CycleGAN. When applied in dose calculation, the gamma passing rate (1%/1 mm criterion) was 97.5% ± 1.1% for TransCBCT, compared to 96.9% ± 1.8% for CycleGAN.

Conclusions

The proposed TransCBCT can effectively generate sCT for CBCT. It has the potential to improve radiotherapy accuracy.

Integrating Structure Propagation Uncertainties in the Optimization of Online Adaptive Proton Therapy Plans

Currently, adaptive strategies require time- and resource-intensive manual structure corrections. This study compares different strategies: optimization without manual structure correction, adaptation with physician-drawn structures, and no adaptation. Strategies were compared for 16 patients with pancreas, liver, and head and neck (HN) cancer with 1-5 repeated images during treatment: 'reference adaptation', with structures drawn by a physician; 'single-DIR adaptation', using a single set of deformably propagated structures; 'multi-DIR adaptation', using robust planning with multiple deformed structure sets; 'conservative adaptation', using the intersection and union of all deformed structures; 'probabilistic adaptation', using the probability of a voxel belonging to the structure in the optimization weight; and 'no adaptation'. Plans were evaluated using reference structures and compared using a scoring system. The reference adaptation with physician-drawn structures performed best, and no adaptation performed the worst. For pancreas and liver patients, adaptation with a single DIR improved the plan quality over no adaptation. For HN patients, integrating structure uncertainties brought an additional benefit. If resources for manual structure corrections would prevent online adaptation, manual correction could be replaced by a fast 'plausibility check', and plans could be adapted with correction-free adaptation strategies. Including structure uncertainties in the optimization has the potential to make online adaptation more automatable.

Reflections on beam configuration optimization for intensity-modulated proton therapy

Presumably, intensity-modulated proton radiotherapy (IMPT) is the most powerful form of proton radiotherapy. In the current state of the art, IMPT beam configurations (i.e. the number of beams and their directions) are, in general, chosen subjectively based on prior experience and practicality. Beam configuration optimization (BCO) for IMPT could, in theory, significantly enhance IMPT’s therapeutic potential. However, BCO is complex and highly computer resource-intensive. Some algorithms for BCO have been developed for intensity-modulated photon therapy (IMRT). They are rarely used clinically mainly because the large number of beams typically employed in IMRT renders BCO essentially unnecessary. Moreover, in the newer form of IMRT, volumetric modulated arc therapy, there are no individual static beams. BCO is of greater importance for IMPT because it typically employs a very small number of beams (2-4) and, when the number of beams is small, BCO is critical for improving plan quality. However, the unique properties and requirements of protons, particularly in IMPT, make BCO challenging. Protons are more sensitive than photons to anatomic changes, exhibit variable relative biological effectiveness along their paths, and, as recently discovered, may spare the immune system. Such factors must be considered in IMPT BCO, though doing so would make BCO more resource intensive and make it more challenging to extend BCO algorithms developed for IMRT to IMPT. A limited amount of research in IMPT BCO has been conducted; however, considerable additional work is needed for its further development to make it truly effective and computationally practical. This article aims to provide a review of existing BCO algorithms, most of which were developed for IMRT, and addresses important requirements specific to BCO for IMPT optimization that necessitate the modification of existing approaches or the development of new effective and efficient ones.

Proton Therapy for Squamous Cell Carcinoma of the Head and Neck: Early Clinical Experience and Current Challenges

Simple Summary

Proton therapy is a promising type of radiation therapy used to destroy tumor cells. It has the potential to further improve the outcomes for patients with head and neck cancer since it allows to minimize the radiation dose to vital structures around the tumor, leading to less toxicity. This paper describes the current experience worldwide with proton therapy in head and neck cancer.

Abstract

Proton therapy (PT) is a promising development in radiation oncology, with the potential to further improve outcomes for patients with squamous cell carcinoma of the head and neck (HNSCC). By utilizing the finite range of protons, healthy tissue can be spared from beam exit doses that would otherwise be irradiated with photon-based treatments. Current evidence on PT for HNSCC is limited to comparative dosimetric analyses and retrospective single-institution series. As a consequence, the recognized indications for the reimbursement of PT remain scarce in most countries. Nevertheless, approximately 100 PT centers are in operation worldwide, and initial experiences for HNSCC are being reported. This review aims to summarize the results of the early clinical experience with PT for HNSCC and the challenges that are currently faced.

Millisecond speed deep learning based proton dose calculation with Monte Carlo accuracy

Objective.Next generation online and real-time adaptive radiotherapy workflows require precise particle transport simulations in sub-second times, which is unfeasible with current analytical pencil beam algorithms (PBA) or Monte Carlo (MC) methods. We present a deep learning based millisecond speed dose calculation algorithm (DoTA) accurately predicting the dose deposited by mono-energetic proton pencil beams for arbitrary energies and patient geometries.Approach.Given the forward-scattering nature of protons, we frame 3D particle transport as modeling a sequence of 2D geometries in the beam's eye view. DoTA combines convolutional neural networks extracting spatial features (e.g. tissue and density contrasts) with a transformer self-attention backbone that routes information between the sequence of geometry slices and a vector representing the beam's energy, and is trained to predict low noise MC simulations of proton beamlets using 80 000 different head and neck, lung, and prostate geometries.Main results.Predicting beamlet doses in 5 ± 4.9 ms with a very high gamma pass rate of 99.37 ± 1.17% (1%, 3 mm) compared to the ground truth MC calculations, DoTA significantly improves upon analytical pencil beam algorithms both in precision and speed. Offering MC accuracy 100 times faster than PBAs for pencil beams, our model calculates full treatment plan doses in 10-15 s depending on the number of beamlets (800-2200 in our plans), achieving a 99.70 ± 0.14% (2%, 2 mm) gamma pass rate across 9 test patients.Significance.Outperforming all previous analytical pencil beam and deep learning based approaches, DoTA represents a new state of the art in data-driven dose calculation and can directly compete with the speed of even commercial GPU MC approaches. Providing the sub-second speed required for adaptive treatments, straightforward implementations could offer similar benefits to other steps of the radiotherapy workflow or other modalities such as helium or carbon treatments.

GPU-accelerated Monte Carlo-based online adaptive proton therapy - a feasibility study

PURPOSE

To develop an online Graphic-Processing-Unit (GPU)-accelerated Monte-Carlo-based adaptive radiation therapy (ART) workflow for pencil beam scanning (PBS) proton therapy to address inter-fraction anatomical changes in patients treated with PBS.

METHODS AND MATERIALS

A four-step workflow was developed using our in-house developed GPU-accelerated Monte-Carlo-based treatment planning system to implement online Monte-Carlo-based ART for PBS. The first step conducts diffeomorphic demon-based deformable image registration (DIR) to propagate contours on the initial planning CT (pCT) to the verification CT (vCT) to form a new structure set. The second step performs forward dose calculation of the initial plan on the vCT with the propagated contours after manual approval (possible modifications involved). The third step triggers a re-optimization of the plan depending on whether the verification dose meets the clinical requirements or not. A robust evaluation will be done for both the verification plan in the second step and the re-opotimized plan in the third step. The fourth step involves a two-stage (before and after delivery) patient specific quality assurance (PSQA) of the re-optimized plan. The before-delivery PSQA is to compare the plan dose to the dose calculated using an independent fast open-source Monte Carlo code, MCsquare. The after-delivery PSQA is to compare the plan dose to the dose re-calculated using the log file (spot MU, spot position, and spot energy) collected during the delivery. Jaccard index (JI), Dice similarity coefficients (DSCs), and Hausdorff distance (HD) were used to assess the quality of the propagated contours in the first step. A commercial plan evaluation software, ClearCheck™, was integrated into the workflow to carry out efficient plan evaluation. 3D Gamma analysis was used during the fourth step to ensure the accuracy of the plan dose from re-optimization. Three patients with three different disease sites were chosen to evaluate the feasibility of the online ART workflow for PBS.

RESULTS

For all three patients, the propagated contours were found to have good volume conformance [JI (lowest-highest: 0.833-0.983) and DSC (0.909-0.992)] but sub-optimal boundary coincidence [HD (2.37-20.76 mm)] for organs at risk (OARs). The verification dose evaluated by ClearCheck™ showed significant degradation of the target coverage due to the inter-fractional anatomical changes. Re-optimization on the vCT resulted in great improvement of the plan quality to a clinically acceptable level. 3D Gamma analyses of PSQA confirmed the accuracy of the plan dose before delivery (mean Gamma index = 98.74% with a threshold of 2%/2 mm/10%), and after delivery based on the log files (mean Gamma index = 99.05% with a threshold of 2%/2 mm/10%). The average time cost for the complete execution of the workflow was around 858 seconds, excluding the time for manual intervention.

CONCLUSION

The proposed online ART workflow for PBS was demonstrated to be efficient and effective by generating a re-optimized plan that significantly improved the plan quality. This article is protected by copyright. All rights reserved.

A two-step method to improve image quality of CBCT with phantom-based supervised and patient-based unsupervised learning strategies

Objective.In this study, we aimed to develop deep learning framework to improve cone-beam computed tomography (CBCT) image quality for adaptive radiation therapy (ART) applications.Approach.Paired CBCT and planning CT images of 2 pelvic phantoms and 91 patients (15 patients for testing) diagnosed with prostate cancer were included in this study. First, well-matched images of rigid phantoms were used to train a U-net, which is the supervised learning strategy to reduce serious artifacts. Second, the phantom-trained U-net generated intermediate CT images from the patient CBCT images. Finally, a cycle-consistent generative adversarial network (CycleGAN) was trained with intermediate CT images and deformed planning CT images, which is the unsupervised learning strategy to learn the style of the patient images for further improvement. When testing or applying the trained model on patient CBCT images, the intermediate CT images were generated from the original CBCT image by U-net, and then the synthetic CT images were generated by the generator of CycleGAN with intermediate CT images as input. The performance was compared with conventional methods (U-net/CycleGAN alone trained with patient images) on the test set.Results.The proposed two-step method effectively improved the CBCT image quality to the level of CT scans. It outperformed conventional methods for region-of-interest contouring and HU calibration, which are important to ART applications. Compared with the U-net alone, it maintained the structure of CBCT. Compared with CycleGAN alone, our method improved the accuracy of CT number and effectively reduced the artifacts, making it more helpful for identifying the clinical target volume.Significance.This novel two-step method improves CBCT image quality by combining phantom-based supervised and patient-based unsupervised learning strategies. It has immense potential to be integrated into the ART workflow to improve radiotherapy accuracy.