Automatic reconstruction of the delivered dose of the day using MR-linac treatment log files and online MR imaging

Development and validation of an MR-driven dose-of-the-day procedure for online adaptive radiotherapy in upper gastrointestinal cancer patients

Objective.To develop and validate a dose-of-the-day (DOTD) treatment plan verification procedure for liver and pancreas cancer patients treated with an magnetic resonance (MR)-Linac system.Approach.DOTD was implemented as an automated process that uses 3D datasets collected during treatment delivery. Particularly, the DOTD pipeline's input included the adapt-to-shape (ATS) plan-i.e. 3D-MR dataset acquired at beginning of online session, anatomical contours, dose distribution-and 3D-MR dataset acquired during beam-on (BON). The DOTD automated analysis included (a) ATS-to-BON image intensity-based deformable image registration (DIR), (b) ATS-to-BON contours mapping via DIR, (c) BON-to-ATS contours copying through rigid registration, (d) determining ATS-to-BON dosimetric differences, and (e) PDF report generation. The DIR process was validated by two expert reviewers. ATS-plans were recomputed on BON datasets to assess dose differences. DOTD analysis was performed retrospectively for 75 treatment fractions (12-liver and 5-pancreas patients).Main results.The accuracy of DOTD process relied on DIR and mapped contours quality. Most DIR-generated contours (99.6%) were clinically acceptable. DICE correlated with depreciation of DIR-based region of interest mapping process. The ATS-BON plan difference was found negligible (<1%). The duodenum and large bowel exhibited highest variations, 24% and 39% from fractional values, for 5-fraction liver and pancreas. For liver 1-fraction, a 62% variation was observed for duodenum.Significance.The DOTD methodology provides an automated approach to quantify 3D dosimetric differences between online plans and their delivery. This analysis offers promise as a valuable tool for plan quality assessment and decision-making in the verification stage of the online workflow.

Magnetic resonance image-guided adaptive radiotherapy enables safe CTV-to-PTV margin reduction in prostate cancer: a cine MRI motion study

Introduction

We aimed to establish if stereotactic body radiotherapy to the prostate can be delivered safely using reduced clinical target volume (CTV) to planning target volume (PTV) margins on the 1.5T MR-Linac (MRL) (Elekta, Stockholm, Sweden), in the absence of gating.

Methods

Cine images taken in 3 orthogonal planes during the delivery of prostate SBRT with 36.25 Gray (Gy) in 5 fractions on the MRL were analysed. Using the data from 20 patients, the percentage of radiotherapy (RT) delivery time where the prostate position moved beyond 1, 2, 3, 4 and 5 mm in the left-right (LR), superior-inferior (SI), anterior-posterior (AP) and any direction was calculated.

Results

The prostate moved less than 3 mm in any direction for 90% of the monitoring period in 95% of patients. On a per-fraction basis, 93% of fractions displayed motion in all directions within 3 mm for 90% of the fraction delivery time. Recurring motion patterns were observed showing that the prostate moved with shallow drift (most common), transient excursions and persistent excursions during treatment.

Conclusion

A 3 mm CTV-PTV margin is safe to use for the treatment of 5 fraction prostate SBRT on the MRL, without gating. In the context of gating this work suggests that treatment time will not be extensively lengthened when an appropriate gating window is applied.

Treatment planning evaluation and experimental validation of the magnetic resonance-based intrafraction drift correction

Background and purpose

MRI-guided online adaptive treatments can account for interfractional variations, however intrafraction motion reduces treatment accuracy. Intrafraction plan adaptation methods, such as the Intrafraction Drift Correction (IDC) or sub-fractionation, are needed. IDC uses real-time automatic monitoring of the tumor position to initiate plan adaptations by repositioning segments. IDC is a fast adaptation method that occurs only when necessary and this method could enable margin reduction. This research provides a treatment planning evaluation and experimental validation of the IDC.

Materials and methods

An in silico treatment planning evaluation was performed for 13 prostate patients mid-treatment without and with intrafraction plan adaptation (IDC and sub-fractionation). The adaptation methods were evaluated using dose volume histogram (DVH) metrics. To experimentally verify IDC a treatment was mimicked whereby a motion phantom containing an EBT3 film moved mid-treatment, followed by repositioning of segments. In addition, the delivered treatment was irradiated on a diode array phantom for plan quality assurance purposes.

Results

The planning study showed benefits for using intrafraction adaptation methods relative to no adaptation, where the IDC and sub-fractionation showed consistently improved target coverage with median target coverages of 100.0%. The experimental results verified the IDC with high minimum gamma passing rates of 99.1% and small mean dose deviations of maximum 0.3%.

Conclusion

The straightforward and fast IDC technique showed DVH metrics consistent with the sub-fractionation method using segment weight re-optimization for prostate patients. The dosimetric and geometric accuracy was shown for a full IDC workflow using film and diode array dosimetry.

Real-time motion-including dose estimation of simulated multi-leaf collimator-tracked magnetic resonance-guided radiotherapy

BACKGROUND

Real-time dose estimation is a key-prerequisite to enable online intra-fraction treatment adaptation in magnetic resonance (MR)-guided radiotherapy (MRgRT). It is an essential component for the assessment of the dosimetric benefits and risks of online adaptive treatments, such as multi-leaf collimator (MLC)-tracking.

PURPOSE

We present a proof-of-concept for a software workflow for real-time dose estimation of MR-guided adaptive radiotherapy based on real-time data-streams of the linac delivery parameters and target positions.

METHODS

A software workflow, combining our in-house motion management software DynaTrack, a real-time dose calculation engine that connects to a research version of the treatment planning software (TPS) Monaco (v.6.09.00, Elekta AB, Stockholm, Sweden) was developed and evaluated. MR-guided treatment delivery on the Elekta Unity MR-linac was simulated with and without MLC-tracking for three prostate patients, previously treated on the Elekta Unity MR-linac (36.25 Gy/five fractions). Three motion scenarios were used: no motion, regular motion, and erratic prostate motion. Accumulated monitor units (MUs), centre of mass target position and MLC-leaf positions, were forwarded from DynaTrack at a rate of 25 Hz to a Monte Carlo (MC) based dose calculation engine which utilises the research GPUMCD-library (Elekta AB, Stockholm, Sweden). A rigid isocentre shift derived from the selected motion scenarios was applied to a bulk density-assigned session MR-image. The respective electron density used for treatment planning was accessed through the research Monaco TPS. The software workflow including the online dose reconstruction was validated against offline dose reconstructions. Our investigation showed that MC-based real-time dose calculations that account for all linac states (including MUs, MLC positions and target position) were infeasible, hence states were randomly sampled and used for calculation as follows; Once a new linac state was received, a dose calculation with 106 photons was started. Linac states that arrived during the time of the ongoing calculation were put into a queue. After completion of the ongoing calculation, one new linac state was randomly picked from the queue and assigned the MU accumulated from the previous state until the last sample in the queue. The queue was emptied, and the process repeated throughout treatment simulation.

RESULTS

On average 27% (23%-30%) of received samples were used in the real-time calculation, corresponding to a calculation time for one linac state of 148 ms. Median gamma pass rate (2%/3 mm local) was 100.0% (99.9%-100%) within the PTV volume and 99.1% (90.1%-99.4.0%) with a 15% dose cut off. Differences in PTVDmean , CTVDmean , RectumD2% , and BladderD2% (offline-online, % of prescribed dose) were below 0.64%. Beam-by-beam comparisons showed deviations below 0.07 Gy. Repeated simulations resulted in standard deviations below 0.31% and 0.12 Gy for the investigated volume and dose criteria respectively.

CONCLUSIONS

Real-time dose estimation was successfully performed using the developed software workflow for different prostate motion traces with and without MLC-tracking. Negligible dosimetric differences were seen when comparing online and offline reconstructed dose, enabling online intra-fraction treatment decisions based on estimates of the delivered dose.

Adaptive ultra-hypofractionated whole-pelvic radiotherapy in high-risk and very high-risk prostate cancer on 1.5-Tesla MR-Linac: Estimated delivered dose and early toxicity results

Background

Magnetic resonance (MR)-guided ultra-hypofractionated radiotherapy with whole-pelvic irradiation (UHF-WPRT) is a novel approach to radiotherapy for patients with high-risk (HR) and very high-risk (VHR) prostate cancer (PCa). However, the inherent complexity of adaptive UHF-WPRT might inevitably result in longer on-couch time. We aimed to estimate the delivered dose, study the feasibility and safety of adaptive UHF-WPRT on a 1.5-Tesla MR-Linac.

Methods

Ten patients with clinical stage T3a-4N0-1M0-1c PCa, who consecutively received UHF-WPRT, were enrolled prospectively. The contours of the target and organ-at-risks on the position verification-MR (PV-MR), beam-on 3D-MR(Bn-MR), and post-MR (after radiotherapy delivery) were derived from the pre-MR data by deformable image registration. The physician then manually adjusted them, and dose recalculation was performed accordingly. GraphPad Prism 9 (GraphPad Prism Software Inc.) was utilized for conducting statistical analyses.

Results

In total, we collected 188 MR scans (50 pre-MR, 50 PV-MR, 44 Bn-MR, and 44 post-MR scans). With median 59 min, the mean prostate clinical target volume (CTV)-V100% was 98.59% ± 2.74%, and the mean pelvic CTVp-V100% relative percentages of all scans was 99.60% ± 1.18%. The median V29 Gy change in the rectal wall was -2% (-18% to 20%). With a median follow-up of 9 months, no patient had acute Common Terminology Criteria for Adverse Events (CTCAE) grade 2 or more severe genitourinary (GU) or gastrointestinal (GI) toxicities (0%).

Conclusion

UHF-RT to the prostate and the whole pelvis with concomitant boost to positive nodes using an Adapt-To-Shape (ATS) workflow was technically feasible for patients with HR and VHR PCa, presenting only mild GU and GI toxicities. The estimated target dose during the beam-on phase was clinically acceptable based on the 3D-MR-based dosimetry analysis.

Clinical trial registration

Chinese Clinical Trial Registry ChiCTR2000033382.

Respiratory motion modelling for MR-guided lung cancer radiotherapy: model development and geometric accuracy evaluation

Objective.Respiratory motion of lung tumours and adjacent structures is challenging for radiotherapy. Online MR-imaging cannot currently provide real-time volumetric information of the moving patient anatomy, therefore limiting precise dose delivery, delivered dose reconstruction, and downstream adaptation methods.Approach.We tailor a respiratory motion modelling framework towards an MR-Linac workflow to estimate the time-resolved 4D motion from real-time data. We develop a multi-slice acquisition scheme which acquires thick, overlapping 2D motion-slices in different locations and orientations, interleaved with 2D surrogate-slices from a fixed location. The framework fits a motion model directly to the input data without the need for sorting or binning to account for inter- and intra-cycle variation of the breathing motion. The framework alternates between model fitting and motion-compensated super-resolution image reconstruction to recover a high-quality motion-free image and a motion model. The fitted model can then estimate the 4D motion from 2D surrogate-slices. The framework is applied to four simulated anthropomorphic datasets and evaluated against known ground truth anatomy and motion. Clinical applicability is demonstrated by applying our framework to eight datasets acquired on an MR-Linac from four lung cancer patients.Main results.The framework accurately reconstructs high-quality motion-compensated 3D images with 2 mm3isotropic voxels. For the simulated case with the largest target motion, the motion model achieved a mean deformation field error of 1.13 mm. For the patient cases residual error registrations estimate the model error to be 1.07 mm (1.64 mm), 0.91 mm (1.32 mm), and 0.88 mm (1.33 mm) in superior-inferior, anterior-posterior, and left-right directions respectively for the building (application) data.Significance.The motion modelling framework estimates the patient motion with high accuracy and accurately reconstructs the anatomy. The image acquisition scheme can be flexibly integrated into an MR-Linac workflow whilst maintaining the capability of online motion-management strategies based on cine imaging such as target tracking and/or gating.

Impact of technological advances in treatment planning, image guidance, and treatment delivery on target margin design for prostate cancer radiotherapy: an updated review

Recent innovations in image guidance, treatment delivery, and adaptive radiotherapy (RT) have created a new paradigm for planning target volume (PTV) margin design for patients with prostate cancer. We performed a review of the recent literature on PTV margin selection and design for intact prostate RT, excluding post-operative RT, brachytherapy, and proton therapy. Our review describes the increased focus on prostate and seminal vesicles as heterogenous deforming structures with further emergence of intra-prostatic GTV boost and concurrent pelvic lymph node treatment. To capture recent innovations, we highlight the evolution in cone beam CT guidance, and increasing use of MRI for improved target delineation and image registration and supporting online adaptive RT. Moreover, we summarize new and evolving image-guidance treatment platforms as well as recent reports of novel immobilization strategies and motion tracking. Our report also captures recent implementations of artificial intelligence to support image guidance and adaptive RT. To characterize the clinical impact of PTV margin changes via model-based risk estimates and clinical trials, we highlight recent high impact reports. Our report focusses on topics in the context of PTV margins but also showcase studies attempting to move beyond the PTV margin recipes with robust optimization and probabilistic planning approaches. Although guidelines exist for target margins conventional using CT-based image guidance, further validation is required to understand the optimal margins for online adaptation either alone or combined with real-time motion compensation to minimize systematic and random uncertainties in the treatment of patients with prostate cancer.

Simulating an intra-fraction adaptive workflow to enable PTV margin reduction in MRIgART volumetric modulated arc therapy for prostate SBRT

Purpose

This study simulates a novel prostate SBRT intra-fraction re-optimization workflow in MRIgART to account for prostate intra-fraction motion and evaluates the dosimetric benefit of reducing PTV margins.

Materials and methods

VMAT prostate SBRT treatment plans were created for 10 patients using two different PTV margins, one with a 5 mm margin except 3 mm posteriorly (standard) and another using uniform 2 mm margins (reduced). All plans were prescribed to 36.25 Gy in 5 fractions and adapted onto each daily MRI dataset. An intra-fraction adaptive workflow was simulated for the reduced margin group by synchronizing the radiation delivery with target position from cine MRI imaging. Intra-fraction delivered dose was reconstructed and prostate DVH metrics were evaluated under three conditions for the reduced margin plans: Without motion compensation (no-adapt), with a single adapt prior to treatment (ATP), and lastly for intra-fraction re-optimization during delivery (intra). Bladder and rectum DVH metrics were compared between the standard and reduced margin plans.

Results

As expected, rectum V18 Gy was reduced by 4.4 ± 3.9%, D1cc was reduced by 12.2 ± 6.8% (3.4 ± 2.3 Gy), while bladder reductions were 7.8 ± 5.6% for V18 Gy, and 9.6 ± 7.3% (3.4 ± 2.5 Gy) for D1cc for the reduced margin reference plans compared to the standard PTV margin. For the intrafraction replanning approach, average intra-fraction optimization times were 40.0 ± 2.9 seconds, less than the time to deliver one of the four VMAT arcs (104.4 ± 9.3 seconds) used for treatment delivery. When accounting for intra-fraction motion, prostate V36.25 Gy was on average 96.5 ± 4.0%, 99.1 ± 1.3%, and 99.6 ± 0.4 for the non-adapt, ATP, and intra-adapt groups, respectively. The minimum dose received by the prostate was less than 95% of the prescription dose in 84%, 36%, and 10% of fractions, for the non-adapt, ATP, and intra-adapt groups, respectively.

Conclusions

Intra-fraction re-optimization improves prostate coverage, specifically the minimum dose to the prostate, and enables PTV margin reduction and subsequent OAR sparing. Fast re-optimizations enable uninterrupted treatment delivery.

MRI-Guided Radiotherapy for Prostate Cancer: Seeing is Believing

The advent of MRI guided radiotherapy (MRIgRT) offers enormous promise in the treatment of prostate cancer. The MR-linac offers men the opportunity to receive daily MR imaging to guide and influence their radiotherapy treatment. This review focuses on the advantages that MRIgRT potentially offers as well as any potential disadvantages to MRIgRT that may have been recognized thus far. Ongoing clinical trials evaluating this novel treatment platform for the treatment of prostate cancer are also discussed.

Investigation of autosegmentation techniques on T2-weighted MRI for off-line dose reconstruction in MR-linac workflow for head and neck cancers

BACKGROUND

In order to accurately accumulate delivered dose for head and neck cancer patients treated with the Adapt to Position workflow on the 1.5T magnetic resonance imaging (MRI)-linear accelerator (MR-linac), the low-resolution T2-weighted MRIs used for daily setup must be segmented to enable reconstruction of the delivered dose at each fraction.

PURPOSE

In this pilot study, we evaluate various autosegmentation methods for head and neck organs at risk (OARs) on on-board setup MRIs from the MR-linac for off-line reconstruction of delivered dose.

METHODS

Seven OARs (parotid glands, submandibular glands, mandible, spinal cord, and brainstem) were contoured on 43 images by seven observers each. Ground truth contours were generated using a simultaneous truth and performance level estimation (STAPLE) algorithm. Twenty total autosegmentation methods were evaluated in ADMIRE: 1-9) atlas-based autosegmentation using a population atlas library (PAL) of 5/10/15 patients with STAPLE, patch fusion (PF), random forest (RF) for label fusion; 10-19) autosegmentation using images from a patient's 1-4 prior fractions (individualized patient prior [IPP]) using STAPLE/PF/RF; 20) deep learning (DL) (3D ResUNet trained on 43 ground truth structure sets plus 45 contoured by one observer). Execution time was measured for each method. Autosegmented structures were compared to ground truth structures using the Dice similarity coefficient, mean surface distance (MSD), Hausdorff distance (HD), and Jaccard index (JI). For each metric and OAR, performance was compared to the inter-observer variability using Dunn's test with control. Methods were compared pairwise using the Steel-Dwass test for each metric pooled across all OARs. Further dosimetric analysis was performed on three high-performing autosegmentation methods (DL, IPP with RF and 4 fractions [IPP_RF_4], IPP with 1 fraction [IPP_1]), and one low-performing (PAL with STAPLE and 5 atlases [PAL_ST_5]). For five patients, delivered doses from clinical plans were recalculated on setup images with ground truth and autosegmented structure sets. Differences in maximum and mean dose to each structure between the ground truth and autosegmented structures were calculated and correlated with geometric metrics.

RESULTS

DL and IPP methods performed best overall, all significantly outperforming inter-observer variability and with no significant difference between methods in pairwise comparison. PAL methods performed worst overall; most were not significantly different from the inter-observer variability or from each other. DL was the fastest method (33 s per case) and PAL methods the slowest (3.7-13.8 min per case). Execution time increased with a number of prior fractions/atlases for IPP and PAL. For DL, IPP_1, and IPP_RF_4, the majority (95%) of dose differences were within ± 250 cGy from ground truth, but outlier differences up to 785 cGy occurred. Dose differences were much higher for PAL_ST_5, with outlier differences up to 1920 cGy. Dose differences showed weak but significant correlations with all geometric metrics (R2 between 0.030 and 0.314).

CONCLUSIONS

The autosegmentation methods offering the best combination of performance and execution time are DL and IPP_1. Dose reconstruction on on-board T2-weighted MRIs is feasible with autosegmented structures with minimal dosimetric variation from ground truth, but contours should be visually inspected prior to dose reconstruction in an end-to-end dose accumulation workflow.

The Quality Assurance of a 1.5 T MR-Linac

The recent introduction of a commercial 1.5 T MR-linac system has considerably improved the image quality of the patient acquired in the treatment unit as well as enabling online adaptive radiation therapy (oART) treatment strategies. Quality Assurance (QA) of this new technology requires new methodology that allows for the high field MR in a linac environment. The presence of the magnetic field requires special attention to the phantoms, detectors, and tools to perform QA. Due to the design of the system, the integrated megavoltage imager (MVI) is essential for radiation beam calibrations and QA. Additionally, the alignment between the MR image system and the radiation isocenter must be checked. The MR-linac system has vendor-supplied phantoms for calibration and QA tests. However, users have developed their own routine QA systems to independently check that the machine is performing as required, as to ensure we are able to deliver the intended dose with sufficient certainty. The aim of this work is therefore to review the MR-linac specific QA procedures reported in the literature.

Continuous time-resolved estimated synthetic 4D-CTs for dose reconstruction of lung tumor treatments at a 0.35 T MR-linac

Objective.To experimentally validate a method to create continuous time-resolved estimated synthetic 4D-computed tomography datasets (tresCTs) based on orthogonal cine MRI data for lung cancer treatments at a magnetic resonance imaging (MRI) guided linear accelerator (MR-linac).Approach.A breathing porcine lung phantom was scanned at a CT scanner and 0.35 T MR-linac. Orthogonal cine MRI series (sagittal/coronal orientation) at 7.3 Hz, intersecting tumor-mimicking gelatin nodules, were deformably registered to mid-exhale 3D-CT and 3D-MRI datasets. The time-resolved deformation vector fields were extrapolated to 3D and applied to a reference synthetic 3D-CT image (sCTref), while accounting for breathing phase-dependent lung density variations, to create 82 s long tresCTs at 3.65 Hz. Ten tresCTs were created for ten tracked nodules with different motion patterns in two lungs. For each dataset, a treatment plan was created on the mid-exhale phase of a measured ground truth (GT) respiratory-correlated 4D-CT dataset with the tracked nodule as gross tumor volume (GTV). Each plan was recalculated on the GT 4D-CT, randomly sampled tresCT, and static sCTrefimages. Dose distributions for corresponding breathing phases were compared in gamma (2%/2 mm) and dose-volume histogram (DVH) parameter analyses.Main results.The mean gamma pass rate between all tresCT and GT 4D-CT dose distributions was 98.6%. The mean absolute relative deviations of the tresCT with respect to GT DVH parameters were 1.9%, 1.0%, and 1.4% for the GTVD98%,D50%, andD2%, respectively, 1.0% for the remaining nodulesD50%, and 1.5% for the lungV20Gy. The gamma pass rate for the tresCTs was significantly larger (p< 0.01), and the GTVD50%deviations with respect to the GT were significantly smaller (p< 0.01) than for the sCTref.Significance.The results suggest that tresCTs could be valuable for time-resolved reconstruction and intrafractional accumulation of the dose to the GTV for lung cancer patients treated at MR-linacs in the future.

Impact of abdominal compression on intra-fractional motion and delivered dose in magnetic resonance image-guided adaptive radiation ablation of adrenal gland metastases

PURPOSE

The current study investigated the impact of abdominal compression on motion and the delivered dose during non-gated, magnetic resonance image (MRI)-guided radiation ablation of adrenal gland metastases.

METHODS

Thirty-one patients with adrenal gland metastases treated to 45-60 Gy in 3-8 fractions on a 1.5 T MRI-linac were included in the study. The patients were breathing freely (n = 14) or with motion restricted by using an abdominal compression belt (n = 17). The time-resolved position of the target in online 2D cine MR images acquired during treatment was assessed and used to estimate the dose delivered to the GTV and abutting luminal organs at risk (OAR).

RESULTS

The median (range) 3D root-mean-square target position error was significantly higher in patients treated without a compression belt [2.9 (1.9-5.6) mm] compared to patients using the belt [2.1 (1.2-3.5) mm] (P < 0.01). The median (range) GTV V95% was significantly reduced from planned 98.6 (65.9-100) % to delivered 96.5 (64.5-99.9) % due to motion (P < 0.01). Most prominent dose reductions were found in patients showing either large target drift or respiration motion and were mainly treated without abdominal compression. Motion did not lead to an increased number of constraint violations for luminal OAR.

CONCLUSIONS

Acceptable target coverage and dose to OAR was observed in the vast majority of patients despite intra-fractional motion during adaptive MRI-guided radiation ablation. The use of abdominal compression significantly reduced the target position error and prevented the most prominent target coverage degradations and is, therefore, recommended as motion management at MRI-linacs.

IPEM topical report: guidance for the clinical implementation of online treatment monitoring solutions for IMRT/VMAT

This report provides guidance for the implementation of online treatment monitoring (OTM) solutions in radiotherapy (RT), with a focus on modulated treatments. Support is provided covering the implementation process, from identification of an OTM solution to local implementation strategy. Guidance has been developed by a RT special interest group (RTSIG) working party (WP) on behalf of the Institute of Physics and Engineering in Medicine (IPEM). Recommendations within the report are derived from the experience of the WP members (in consultation with manufacturers, vendors and user groups), existing guidance or legislation and a UK survey conducted in 2020 (Stevenset al2021). OTM is an inclusive term representing any system capable of providing a direct or inferred measurement of the delivered dose to a RT patient. Information on each type of OTM is provided but, commensurate with UK demand, guidance is largely influenced byin vivodosimetry methods utilising the electronic portal imager device (EPID). Sections are included on the choice of OTM solutions, acceptance and commissioning methods with recommendations on routine quality control, analytical methods and tolerance setting, clinical introduction and staffing/resource requirements. The guidance aims to give a practical solution to sensitivity and specificity testing. Functionality is provided for the user to introduce known errors into treatment plans for local testing. Receiver operating characteristic analysis is discussed as a tool to performance assess OTM systems. OTM solutions can help verify the correct delivery of radiotherapy treatment. Furthermore, modern systems are increasingly capable of providing clinical decision-making information which can impact the course of a patient's treatment. However, technical limitations persist. It is not within the scope of this guidance to critique each available solution, but the user is encouraged to carefully consider workflow and engage with manufacturers in resolving compatibility issues.

Accelerating 4D image reconstruction for magnetic resonance-guided radiotherapy

Background and purpose

Physiological motion impacts the dose delivered to tumours and vital organs in external beam radiotherapy and particularly in particle therapy. The excellent soft-tissue demarcation of 4D magnetic resonance imaging (4D-MRI) could inform on intra-fractional motion, but long image reconstruction times hinder its use in online treatment adaptation. Here we employ techniques from high-performance computing to reduce 4D-MRI reconstruction times below two minutes to facilitate their use in MR-guided radiotherapy.

Material and methods

Four patients with pancreatic adenocarcinoma were scanned with a radial stack-of-stars gradient echo sequence on a 1.5T MR-Linac. Fast parallelised open-source implementations of the extra-dimensional golden-angle radial sparse parallel algorithm were developed for central processing unit (CPU) and graphics processing unit (GPU) architectures. We assessed the impact of architecture, oversampling and respiratory binning strategy on 4D-MRI reconstruction time and compared images using the structural similarity (SSIM) index against a MATLAB reference implementation. Scaling and bottlenecks for the different architectures were studied using multi-GPU systems.

Results

All reconstructed 4D-MRI were identical to the reference implementation (SSIM > 0.99). Images reconstructed with overlapping respiratory bins were sharper at the cost of longer reconstruction times. The CPU  + GPU implementation was over 17 times faster than the reference implementation, reconstructing images in 60 ± 1 s and hyper-scaled using multiple GPUs.

Conclusion

Respiratory-resolved 4D-MRI reconstruction times can be reduced using high-performance computing methods for online workflows in MR-guided radiotherapy with potential applications in particle therapy.

Impact of log file source and data frequency on accuracy of log file-based patient specific quality assurance

Performing phantom measurements for patient-specific quality assurance (PSQA) adds a significant amount of time to the adaptive radiotherapy procedure. Log file based PSQA can be used to increase the efficiency of this process. This study compared the dosimetric accuracy of high-frequency linear accelerator (Linac) log files and low-frequency log data stored in the oncology information system (OIS). Thirty patients were included, that were recently treated in the head and neck (HN), brain, and prostate region with volumetric modulated arc therapy (VMAT) and an additional ten patients treated using stereotactic body radiation therapy (SBRT) with 3D-conformal radiotherapy (3D-CRT) technique. Log data containing a single fraction were used to calculate the dose distributions. The dosimetric differences between Linac log files and OIS logs were evaluated with a gamma analysis with 2%/2 mm criterion and dose threshold of 30%. The original treatment plan was used as a reference. Moreover, DVH parameters of D_98%, D_50%, and D_2% of the planning-target volume (PTV) and dose to several organs at risk (OARs) were reported. Significant differences in dose distributions between the two log types and the original dose were observed for PTV D_98% and D_2% (r < 0.001) for HN cases, PTV D_98% (r = 0.005) for brain cases, and PTV D_50% (r = 0.015) for prostate cases. No significant differences were found between the two log types with respect to D_50%. The root mean square (RMS) error of the leaf positions of the OIS log was approximately twice the RMS error of the Linac log file for VMAT plans, but identical for 3D-CRT plans. The relationship between the gamma pass rate and the RMS error showed a moderate correlation for the Linac log files (r = -0.58, p < 0.001) and strong correlation for OIS logs (r = -0.71, p < 0.001). Furthermore, all doses calculated using Linac log files and OIS log data had a GPR >90% for an RMS error < 3.3 mm. Based on these findings, a tolerance limit of RMS error of 3.3 mm for considering OIS log based PSQA was established. Nevertheless, the OIS log data quality should be improved to achieve adequate PSQA.

Online adaptive radiotherapy and dose delivery accuracy: A retrospective analysis

PURPOSE

With online adaptive radiotherapy (ART), patient-specific quality assurance (PSQA) testing cannot be performed prior to delivery of the adapted treatment plan. Consequently, the dose delivery accuracy of adapted plans (i.e., the ability of the system to interpret and deliver the treatment as planned) are not initially verified. We investigated the variation in dose delivery accuracy of ART on the MRIdian 0.35 T MR-linac (Viewray Inc., Oakwood, USA) between initial plans and their respective adapted plans, by analyzing PSQA results.

METHODS

We considered the two main digestive localizations treated with ART (liver and pancreas). A total of 124 PSQA results acquired with the ArcCHECK (Sun Nuclear Corporation, Melbourne, USA) multidetector system were analyzed. PSQA result variations between the initial plans and their respective adapted plans were statistically investigated and compared with the variation in MU number.

RESULTS

For the liver, limited deterioration in PSQA results was observed, and was within the limits of clinical tolerance (Initial = 98.2%, Adapted = 98.2%, p = 0.4503). For pancreas plans, only a few significant deteriorations extending beyond the limits of clinical tolerance were observed and were due to specific, complex anatomical configurations (Initial = 97.3%, Adapted = 96.5%, p = 0.0721). In parallel, we observed an influence of the increase in MU number on the PSQA results.

CONCLUSION

We show that the dose delivery accuracy of adapted plans, in terms of PSQA results, is preserved in ART processes on the 0.35 T MR-linac. Respecting good practices, and minimizing the increase in MU number can help to preserve the accuracy of delivery of adapted plans as compared to their respective initial plans.

Dose accumulation for MR-guided adaptive radiotherapy: From practical considerations to state-of-the-art clinical implementation

MRI-linear accelerator (MR-linac) devices have been introduced into clinical practice in recent years and have enabled MR-guided adaptive radiation therapy (MRgART). However, by accounting for anatomical changes throughout radiation therapy (RT) and delivering different treatment plans at each fraction, adaptive radiation therapy (ART) highlights several challenges in terms of calculating the total delivered dose. Dose accumulation strategies-which typically involve deformable image registration between planning images, deformable dose mapping, and voxel-wise dose summation-can be employed for ART to estimate the delivered dose. In MRgART, plan adaptation on MRI instead of CT necessitates additional considerations in the dose accumulation process because MRI pixel values do not contain the quantitative information used for dose calculation. In this review, we discuss considerations for dose accumulation specific to MRgART and in relation to current MR-linac clinical workflows. We present a general dose accumulation framework for MRgART and discuss relevant quality assurance criteria. Finally, we highlight the clinical importance of dose accumulation in the ART era as well as the possible ways in which dose accumulation can transform clinical practice and improve our ability to deliver personalized RT.

Assessment of delivered dose in prostate cancer patients treated with ultra-hypofractionated radiotherapy on 1.5-Tesla MR-Linac

Objective

To quantitatively characterize the dosimetric effects of long on-couch time in prostate cancer patients treated with adaptive ultra-hypofractionated radiotherapy (UHF-RT) on 1.5-Tesla magnetic resonance (MR)-linac.

Materials and methods

Seventeen patients consecutively treated with UHF-RT on a 1.5-T MR-linac were recruited. A 36.25 Gy dose in five fractions was delivered every other day with a boost of 40 Gy to the whole prostate. We collected data for the following stages: pre-MR, position verification-MR (PV-MR) in the Adapt-To-Shape (ATS) workflow, and 3D-MR during the beam-on phase (Bn-MR) and at the end of RT (post-MR). The target and organ-at-risk contours in the PV-MR, Bn-MR, and post-MR stages were projected from the pre-MR data by deformable image registration and manually adapted by the physician, followed by dose recalculation for the ATS plan.

Results

Overall, 290 MR scans were collected (85 pre-MR, 85 PV-MR, 49 Bn-MR and 71 post-MR scans). With a median on-couch time of 49 minutes, the mean planning target volume (PTV)-V_95% of all scans was 97.83 ± 0.13%. The corresponding mean clinical target volume (CTV)-V_100% was 99.93 ± 0.30%, 99.32 ± 1.20%, 98.59 ± 1.84%, and 98.69 ± 1.85%. With excellent prostate-V_100% dose coverage, the main reason for lower CTV-V_100% was slight underdosing of seminal vesicles (SVs). The median V_{29 Gy} change in the rectal wall was -1% (-20%-17%). The V_{29 Gy} of the rectal wall increased by >15% was observed in one scan. A slight increase in the high dose of bladder wall was noted due to gradual bladder growth during the workflow.

Conclusions

This 3D-MR-based dosimetry analysis demonstrated clinically acceptable estimated dose coverage of target volumes during the beam-on period with adaptive ATS workflow on 1.5-T MR-linac, albeit with a relatively long on-couch time. The 3-mm CTV-PTV margin was adequate for prostate irradiation but occasionally insufficient for SVs. More attention should be paid to restricting high-dose RT to the rectal wall when optimizing the ATS plan.

Intra-fraction motion of pelvic oligometastases and feasibility of PTV margin reduction using MRI guided adaptive radiotherapy

Purpose

This study assesses the impact of intra-fraction motion and PTV margin size on target coverage for patients undergoing radiation treatment of pelvic oligometastases. Dosimetric sparing of the bowel as a function of the PTV margin is also evaluated.

Materials and methods

Seven patients with pelvic oligometastases previously treated on our MR-linac (35 Gy in 5 fractions) were included in this study. Retrospective adaptive plans were created for each fraction on the daily MRI datasets using PTV margins of 5 mm, 3 mm, and 2 mm. Dosimetric constraint violations and GTV coverage were measured as a function of PTV margin size. The impact of intra-fraction motion on GTV coverage was assessed by tracking the GTV position on the cine MR images acquired during treatment delivery and creating an intra-fraction dose distribution for each IMRT beam. The intra-fraction dose was accumulated for each fraction to determine the total dose delivered to the target for each PTV size.

Results

All OAR constraints were achieved in 85.7%, 94.3%, and 100.0% of fractions when using 5 mm, 3 mm, and 2 mm PTV margins while scaling to 95% PTV coverage. Compared to plans with a 5 mm PTV margin, there was a 27.4 ± 12.3% (4.0 ± 2.2 Gy) and an 18.5 ± 7.3% (2.7 ± 1.4 Gy) reduction in the bowel D_0.5cc dose for 2 mm and 3 mm PTV margins, respectively. The target dose (GTV V_{35 Gy}) was on average 100.0 ± 0.1% (99.6 - 100%), 99.6 ± 1.0% (97.2 - 100%), and 99.0 ± 1.4% (95.0 - 100%), among all fractions for the 5 mm, 3 mm, and 2 mm PTV margins on the adaptive plans when accounting for intra-fraction motion, respectively.

Conclusion

A 2 mm PTV margin achieved a minimum of 95% GTV coverage while reducing the dose to the bowel for all patients.

MRI-guided adaptive radiotherapy for prostate cancer: When do we need to account for intra-fraction motion?

A shift of the daily plan can mitigate target position changes that occur between daily MR acquisition and treatment for MR-linac radiotherapy, but increases the session time. We demonstrated that our workflow strategy and decision-making process, to determine whether a subsequent shift is necessary, is appropriate.

Clinical application of a sub-fractionation workflow for intrafraction re-planning during prostate radiotherapy treatment on a 1.5 Tesla MR-Linac: A practical method to mitigate intrafraction motion

BACKGROUND

Intrafraction motion during radiotherapy limits margin reduction and dose escalation. Magnetic resonance (MR)-guided linear accelerators (MR-Linac) have emphasised this issue by enabling intrafraction imaging. We present and clinically apply a new workflow to counteract systematic intrafraction motion during MR-guided stereotactic body radiotherapy (SBRT).

MATERIALS AND METHODS

With the sub-fractionation workflow, the daily dose is delivered in multiple sequential parts (sub-fractions), each adapted to the latest anatomy. As each sub-fractionation treatment plan complies with the dose constraints, no online dose accumulation is required. Imaging and treatment planning are executed in parallel with dose delivery to minimise dead time, enabling an efficient workflow. The workflow was implemented on a 1.5 T MR-Linac and applied in 15 prostate cancer (PCa) patients treated with 5 × 7.25 Gy in two sub-fractions of 3.625 Gy (10 × 3.625 Gy in total). Intrafraction clinical target volume (CTV) motion was determined and compared to a workflow with single-plan delivery. Furthermore, required planning target volume (PTV) margins were determined.

RESULTS

Average on-table time was 42.7 min. Except for two fractions, all fractions were delivered within 60 min. Average intrafraction 3D CTV displacement (±standard deviation) was 1.1 mm (± 0.7) with the sub-fractionation workflow, whereas this was up to 3.5 mm (± 2.4) without sub-fractionation. Calculated PTV margins required with sub-fractionation were 1.0 mm (left-right), 2.4 mm (cranial-caudal), and 2.6 mm (anterior-posterior).

CONCLUSION

Feasibility of the sub-fractionation workflow was demonstrated in 15 PCa patients treated with two sub-fractions on a 1.5 T MR-Linac. The workflow allows for significant PTV margin reduction in these patients by reducing systematic intrafraction motion during SBRT.

Dosimetric robustness of lung tumor photon radiotherapy evaluated from multiple event CT imaging

PURPOSE

Intrafractional respiratory motion is a concern for lung tumor radiotherapy but full evaluation of its impact is hampered by the lack of images representing the true motion. This study presents a novel evaluation using free-breathing images acquired over realistic treatment times to study the dosimetric impact of respiratory motion in photon radiotherapy.

METHODS

Cine-CT images of 14 patients with lung cancer acquired during eight minutes of free-breathing at three occasions were used to simulate dose tracking of four different planning methods. These methods aimed to deliver 54 Gy in three fractions to D_50% of the target and were denoted as robust 4D (RB4), homogeneous fluence to the ITV (FLU) and an isodose prescription to the ITV with a high central dose (ISD), concurrently renormalized (IRN). Differences in dose coverage probability and homogeneity between the methods were quantified. Correlations between underdosage and attributes regarding the tumor and its motion were investigated.

RESULTS

Despite tumor motion amplitudes being larger than in the 4DCT all but FLU achieved the intended CTV D_50% for the cohort average. For all methods but IRN at least 93% of the patients would have received 95% of the intended dose. No differences in D_50% were found between RB4 and ISD nor IRN. However, RB4 led to better homogeneity.

CONCLUSIONS

Tumor motion in free-breathing not covered by the 4DCT had a small impact on dose. The RB4 is recommended for planning of free-breathing treatments. No factor was found that consistently correlated dose degradation with patient or motion attributes.

Assessment of intrafractional prostate motion and its dosimetric impact in MRI-guided online adaptive radiotherapy with gating

PURPOSE

This study aimed to evaluate the intrafractional prostate motion captured during gated magnetic resonance imaging (MRI)-guided online adaptive radiotherapy for prostate cancer and analyze its impact on the delivered dose as well as the effect of gating.

METHODS

Sagittal 2D cine-MRI scans were acquired at 4 Hz during treatment at a ViewRay MRIdian (ViewRay Inc., Oakwood Village, OH, USA) MR linac. Prostate shifts in anterior-posterior (AP) and superior-inferior (SI) directions were extracted separately. Using the static dose cloud approximation, the planned fractional dose was shifted according to the 2D gated motion (residual motion in gating window) to estimate the delivered dose by superimposing and averaging the shifted dose volumes. The dose of a hypothetical non-gated delivery was reconstructed similarly using the non-gated motion. For the clinical target volume (CTV), rectum, and bladder, dose-volume histogram parameters of the planned and reconstructed doses were compared.

RESULTS

In total, 174 fractions (15.7 h of cine-MRI) from 10 patients were evaluated. The average (±1 σ) non-gated prostate motion was 0.6 ± 1.0 mm in the AP and 0.0 ± 0.6 mm in the SI direction with respect to the centroid position of the gating boundary. 95% of the shifts were within [-3.5, 2.7] mm in the AP and [-2.9, 3.2] mm in the SI direction. For the gated treatment and averaged over all fractions, CTV D_98% decreased by less than 2% for all patients. The rectum and the bladder D_2% increased by less than 3% and 0.5%, respectively. Doses reconstructed for gated and non-gated delivery were similar for most fractions.

CONCLUSION

A pipeline for extraction of prostate motion during gated MRI-guided radiotherapy based on 2D cine-MRI was implemented. The 2D motion data enabled an approximate estimation of the delivered dose. For the majority of fractions, the benefit of gating was negligible, and clinical dosimetric constraints were met, indicating safety of the currently adopted gated MRI-guided treatment workflow.

Impact of intrafraction changes in delivered dose of the day for prostate cancer patients treated with stereotactic body radiotherapy via MR-Linac

•

Beam on MR acquisition on the MR-Linac can be used to compute DDOTD.

•

Intrafraction motion via volumetric variability of OARs can impact dosimetry.

•

Computation of the DDOTD may help inform prospective fractions for SBRT prostate.

Purpose

The purpose of this study is to evaluate the impact of intrafraction pelvic motion by comparing the adapted plan dose (APD) and the computed delivered dose of the day (DDOTD) for patients with prostate cancer (PCa) treated with SBRT on the MR-Linac.

Methods

Twenty patients with PCa treated with MR-guided adaptive SBRT were included. A 9-field IMRT distribution was adapted based on the anatomy of the day to deliver a total prescription dose of 3000 cGy in 5 fractions to the prostate plus a 5 mm isotropic margin. Prostate, bladder, and rectum were re-contoured on the MR-image acquired during treatment delivery (MR_BO). DDOTD was computed by propagating the dose from the daily adapted plan generated during treatment onto the MR_BO.

Results

Target coverage was met for all fractions, however, computed DDOTD was significantly less than the APD (p < 0.05). During an average treatment of 53 min, mean bladder volume increased by 116%, which led to a significant decrease in the DDOTD bladder D40% (p < 0.001). However, DDOTD to bladder 5 cc was significantly higher (p < 0.001) than APD. Rectum intrafraction changes were observed based on a volume change of −20% to 83% and presence of significant dose changes from APD to DDOTD for rectum D20% (p < 0.05) and D1cc (p < 0.0001).

Conclusions

Intrafraction motion observed during prostate SBRT treatment on the MR-Linac have dosimetric impacts on both the target and organs at risk. Post-treatment computation using DDOTD may inform adaptation beyond anatomic changes in subsequent treatment fractions to best capitalize on MR-Linac technology and widen the therapeutic index of SBRT for PCa.

MR-guided beam gating: Residual motion, gating efficiency and dose reconstruction for stereotactic treatments of the liver and lung

PURPOSE

This study aims to investigate the efficiency and the geometric as well as the dosimetric benefit of magnetic-resonance guided beam gating for stereotactic treatments in moving organs.

METHOD

Patients treated with MR-guided (MRIdian system) SBRT for lung (n = 10) and liver (n = 10) targets were analyzed. Breath-hold gating was performed based on lesion tracking in sagittal cine MRI images. The target offset from the geometric center of the gating window with and without gating was evaluated. A dose reconstruction workflow based on convolution of these 2D position-probability maps and the daily 3D dose distribution was used to estimate the daily delivered dose including motion. The dose to the clinical target volume (CTV) and to a 2-cm ring structure around the planning target volume were evaluated.

RESULTS

The applied gating protocol resulted in a mean (±standard deviation) gating efficiency of 55%±16%. Over all patients, the mean target offset (2D-root-mean-square error) was 8.3 ± 4.3 mm, which reduced to 2.4 ± 0.6 mm during gating. The dose reconstruction showed a mean deviation in CTV coverage (D95) from the static plans of -1.7%±1.8% with gating and -12.0%±8.4% if no gating would have been used. The mean dose (Dmean) in the ring structure, with respect to the static plans, showed mean deviations of -0.1%±0.3% with gating and -1.6%±1.8% without gating.

CONCLUSION

The MRIdian system enables gating based on the inner anatomy and the implemented dose reconstruction workflow demonstrated geometric robust delivery of the planned radiation doses.

A hybrid 2D/4D-MRI methodology using simultaneous multislice imaging for radiotherapy guidance

Purpose

Respiratory motion management is important in abdominothoracic radiotherapy. Fast imaging of the tumor can facilitate multileaf collimator (MLC) tracking that allows for smaller treatment margins, while repeatedly imaging the full field‐of‐view is necessary for 4D dose accumulation. This study introduces a hybrid 2D/4D‐MRI methodology that can be used for simultaneous MLC tracking and dose accumulation on a 1.5 T Unity MR‐linac (Elekta AB, Stockholm, Sweden).

Methods

We developed a hybrid 2D/4D‐MRI methodology that uses a simultaneous multislice (SMS) accelerated MRI sequence, which acquires two coronal slices simultaneously and repeatedly cycles through slice positions over the image volume. As a result, the fast 2D imaging can be used prospectively for MLC tracking and the SMS slices can be sorted retrospectively into respiratory‐correlated 4D‐MRIs for dose accumulation. Data were acquired in five healthy volunteers with an SMS‐bTFE and SMS‐TSE MRI sequence. For each sequence, a prebeam dataset and a beam‐on dataset were acquired simulating the two phases of MR‐linac treatments. Prebeam data were used to generate a 4D‐based motion model and a reference mid‐position volume, while beam‐on data were used for real‐time motion extraction and reconstruction of beam‐on 4D‐MRIs. In addition, an in‐silico computational phantom was used for validation of the hybrid 2D/4D‐MRI methodology. MLC tracking experiments were performed with the developed methodology, for which real‐time SMS data reconstruction was enabled on the scanner. A 15‐beam 8× 7.5 Gy intensity‐modulated radiotherapy plan for lung stereotactic body radiotherapy with isotropic 3 mm GTV‐to‐PTV margins was created. Dosimetry experiments were performed using a 4D motion phantom. The latency between target motion and updating the radiation beam was determined and compensated. Local gamma analyses were performed to quantify dose differences compared to a static reference delivery, and dose area histograms (DAHs) were used to quantify the GTV and PTV coverage.

Results

In‐vivo data acquisition and MLC tracking experiments were successfully performed with the developed hybrid 2D/4D‐MRI methodology. Real‐time liver–lung interface motion estimation had a Pearson's correlation of 0.996 (in‐vivo) and 0.998 (in‐silico). A median (5th–95th percentile) error of 0.0 (−0.9 to 0.7) mm and 0.0 (−0.2 to 0.2) mm was found for real‐time motion estimation for in‐vivo and in‐silico, respectively. Target motion prediction beyond the liver–lung interface had a median root mean square error of 1.6 mm (in‐vivo) and 0.5 mm (in‐silico). Beam‐on 4D MRI reconstruction required a median amount of data equal to an acquisition time of 2:21–3:17 min, which was 20% less data compared to the prebeam‐derived 4D‐MRI. System latency was reduced from 501 ± 12 ms to −1 ± 3 ms (SMS‐TSE) and from 398 ± 10 ms to −10 ± 4 ms (SMS‐bTFE) by a linear regression prediction filter. The local gamma analysis agreed within −3.8% to 3.3% (SMS‐bTFE) and −5.3% to 10% (SMS‐TSE) with a reference MRI sequence. The DAHs revealed a relative D98% GTV coverage between 97% and 100% (SMS‐bTFE) and 100% and 101% (SMS‐TSE) compared to the static reference.

Conclusions

The presented 2D/4D‐MRI methodology demonstrated the potential for accurately extracting real‐time motion for MLC tracking in abdominothoracic radiotherapy, while simultaneously reconstructing contiguous respiratory‐correlated 4D‐MRIs for dose accumulation.

Predicting gamma evaluation results of patient-specific head and neck volumetric-modulated arc therapy quality assurance based on multileaf collimator patterns and fluence map features: A feasibility study

The purpose of this study was to develop a predictive model for patient-specific VMAT QA results using multileaf collimator (MLC) effect and texture analysis. The MLC speed, acceleration and texture analysis features were extracted from 106 VMAT plans as predictors. Gamma passing rate (GPR) was collected as a response class with gamma criteria of 2%/2 mm and 3%/2 mm. The model was trained using two machine learning methods: AdaBoost classification and bagged regression trees model. GPR was classified into the "PASS" and "FAIL" for the classification model using the institutional warning level. The accuracy of the model was assessed using sensitivity and specificity. In addition, the accuracy of the regression model was determined using the difference between predicted and measured GPR. For the AdaBoost classification model, the sensitivity/specificity was 94.12%/100% and 63.63%/53.13% at gamma criteria of 2%/2 mm and 3%/2 mm, respectively. For the bagged regression trees model, the sensitivity/specificity was 94.12%/91.89% and 61.18%/68.75% at gamma criteria of 2%/2 mm and 3%/2 mm, respectively. The root mean square error (RMSE) of difference between predicted and measured GPR was found at 2.44 and 1.22 for gamma criteria of 2%/2 mm and 3%/2 mm, respectively. The promising result was found at tighter gamma criteria 2%/2 mm with 94.12% sensitivity (both bagged regression trees and AdaBoost classification model) and 100% specificity (AdaBoost classification model).

Offline and online LSTM networks for respiratory motion prediction in MR-guided radiotherapy

Objective. Gated beam delivery is the current clinical practice for respiratory motion compensation in MR-guided radiotherapy, and further research is ongoing to implement tracking. To manage intra-fractional motion using multileaf collimator tracking the total system latency needs to be accounted for in real-time. In this study, long short-term memory (LSTM) networks were optimized for the prediction of superior–inferior tumor centroid positions extracted from clinically acquired 2D cine MRIs. Approach. We used 88 patients treated at the University Hospital of the LMU Munich for training and validation (70 patients, 13.1 h), and for testing (18 patients, 3.0 h). Three patients treated at Fondazione Policlinico Universitario Agostino Gemelli were used as a second testing set (1.5 h). The performance of the LSTMs in terms of root mean square error (RMSE) was compared to baseline linear regression (LR) models for forecasted time spans of 250 ms, 500 ms and 750 ms. Both the LSTM and the LR were trained with offline (offline LSTM and offline LR) and online schemes (offline+online LSTM and online LR), the latter to allow for continuous adaptation to recent respiratory patterns. Main results. We found the offline+online LSTM to perform best for all investigated forecasts. Specifically, when predicting 500 ms ahead it achieved a mean RMSE of 1.20 mm and 1.00 mm, while the best performing LR model achieved a mean RMSE of 1.42 mm and 1.22 mm for the LMU and Gemelli testing set, respectively. Significance. This indicates that LSTM networks have potential as respiratory motion predictors and that continuous online re-optimization can enhance their performance.

Accumulated bladder wall dose is correlated with patient-reported acute urinary toxicity in prostate cancer patients treated with stereotactic, daily adaptive MR-guided radiotherapy

BACKGROUND AND PURPOSE

Magnetic resonance (MR)-guided linear accelerators (MR-Linac) enable accurate estimation of delivered doses through dose accumulation using daily MR images and treatment plans. We aimed to assess the association between the accumulated bladder (wall) dose and patient-reported acute urinary toxicity in prostate cancer (PCa) patients treated with stereotactic body radiation therapy (SBRT).

MATERIALS AND METHODS

One-hundred-and-thirty PCa patients treated on a 1.5T MR-Linac were included. Patients filled out International Prostate Symptom Scores (IPSS) questionnaires at baseline, 1 month, and 3 months post-treatment. Deformable image registration-based dose accumulation was performed to reconstruct the delivered dose. Dose parameters for both bladder and bladder wall were correlated with a clinically relevant increase in IPSS (≥10 points) and/or start of alpha-blockers within 3 months using logistic regression.

RESULTS

Thirty-nine patients (30%) experienced a clinically relevant IPSS increase and/or started with alpha-blockers. Bladder D5cm³, V_10-35Gy (in %), and D_mean and Bladder wall V_10-35Gy (cm³ and %) and D_mean were correlated with the outcome (odds ratios 1.04-1.33, p-values 0.001-0.044). Corrected for baseline characteristics, bladder V_10-35Gy (in %) and D_mean and bladder wall V_10-35Gy (cm³ and %) and D_mean were still correlated with the outcome (odds ratios 1.04-1.30, p-values 0.001-0.028). Bladder wall parameters generally showed larger AUC values.

CONCLUSION

This is the first study to assess the correlation between accumulated bladder wall dose and patient-reported urinary toxicity in PCa patients treated with MR-guided SBRT. The dose to the bladder wall is a promising parameter for prediction of patient-reported urinary toxicity and therefore warrants prospective validation and consideration in treatment planning.

MR-guided radiotherapy for prostate cancer: state of the art and future perspectives

Advances in radiotherapy technology have increased precision of treatment delivery and in some tumour types, improved cure rates and decreased side effects. A new generation of radiotherapy machines, hybrids of an MRI scanner and a linear accelerator, has the potential to further transform the practice of radiation therapy in some cancers. Facilitating superior image quality and the ability to change the dose distribution online on a daily basis (termed "daily adaptive replanning"), MRI-guided radiotherapy machines allow for new possibilities including increasing dose, for hard to treat cancers, and more selective sparing of healthy tissues, where toxicity reduction is the key priority.These machines have already been used to treat most types of cancer, although experience is still in its infancy. This review summarises the potential and current evidence for MRI-guided radiotherapy, with a predominant focus on prostate cancer. Current advantages and disadvantages are discussed including a realistic appraisal of the likely potential to improve patient outcomes. In addition, horizon scanning for near-term possibilities for research and development will hopefully delineate the potential role for this technology over the next decade.

On-line daily plan optimization combined with a virtual couch shift procedure to address intrafraction motion in prostate magnetic resonance guided radiotherapy

Background and purpose

In daily adaptive magnetic resonance (MR)-guided radiotherapy, plans are adapted based on the patient's daily anatomy. During this adaptation phase, prostate intrafraction motion (IM) can occur. The aim of this study was to investigate the efficacy of always applying a subsequent virtual couch shift (VCS) to counter IM that occurred during the daily contour and plan adaption (CPa) procedure.

Material and Methods

One hundred fifty patients with low and intermediate risk prostate cancer were treated with 5x7.25 Gy fractions on a 1.5 T MR-Linac. In each fraction, contour adaptation and dose re-optimization was performed using the session’s first MR-scan. IM that occurred here was countered using two methods. One patient group had selective VCS (sVCS) applied if the CTV reached outside the PTV on a second MR acquired during plan optimization. The other group had always VCS (aVCS) applied for any prostate shift greater than 1 mm. Remaining IM during beam delivery was determined using 3D cine-MR.

Results

Percentage of fractions where a VCS was applied was 28% (sVCS) vs 78% (aVCS). Always applying VCS significantly reduced influences of systematic prostate IM. Population random and systematic median values in all translations directions were lower for the aVCS than sVCS group, but not for the population random cranial-caudal direction.

Conclusion

Applying VCS after daily CPa reduced impact of systematic prostate drift in especially the posterior and caudal translation direction. However, due to the continuous and stochastical nature of prostate IM, margin reduction below 4 mm requires fast intrafraction plan adaption methods.

Seminal vesicle intrafraction motion during the delivery of radiotherapy sessions on a 1.5 T MR-Linac

PURPOSE

To evaluate seminal vesicle (SV) intrafraction motion using cinematic magnetic resonance imaging (cine-MR) during the delivery of online adaptive MR-Linac radiotherapy fractions, in preparation of MR-guided extremely hypofractionated radiotherapy for intermediate to high-risk prostate cancer patients.

MATERIAL AND METHODS

Fifty prostate cancer patients were treated with 5 × 7.25 Gy on a 1.5 Tesla MR-Linac. 3D Cine-MR imaging was started simultaneously and acquired over the full beam-on period. Intrafraction motion in this cine-MR was determined for each SV separately with a previously validated soft-tissue contrast-based tracking algorithm. Motion statistics and coverage probability for the SVs and prostate were determined based on the obtained results.

RESULTS

SV motion was automatically determined during the beam-on period (approx. 10 min) for 247 fractions. SV intrafraction motion shows larger spread than prostate intrafraction motion and increases over time. This difference is especially evident in the anterior and cranial translation directions. Significant difference in rotation about the left-right axis was found, with larger rotation for the SVs than the prostate. Intra-fraction coverage probability of 99% can be achieved when using 5 mm isometric expansion for the left and right SV and 3 mm for the prostate.

CONCLUSION

This is the first study to investigate SV intrafraction motion during MR-guided RT sessions on an MR-Linac. We have shown that high quality 3D cine-MR imaging and SV tracking during RT is feasible with beam-on. The tracking method as described may be used as input for a fast replanning algorithm, which allows for intrafraction plan adaptation.

An investigation of using log-file analysis for automated patient-specific quality assurance in MRgRT

Objective

Adaptive radiation therapy (ART) is an integral part of MR‐guided RT (MRgRT), requiring a new RT plan for each treatment fraction and resulting in a significant increase in patient‐specific quality assurance (PSQA). This study investigates the possibility of using treatment log‐file for automated PSQA.

Method

All treatment plans were delivered in 1.5T Unity MR‐Linac (Elekta). A Unity compatible version of LinacView (Standard Imaging) was commissioned to automatically monitor and analyze the log‐files. A total of 220 fields were delivered and measured by ArcCheck^®‐MR (Sun Nuclear) and LinacView. Thirty incorrectly matched fields were also delivered to check for error detection sensitivity. The gamma analysis, γ, with 3%, 3 mm criteria was used in both ArcCheck^®‐MR and LinacView. Additionally, the gantry angle, jaws, and multileaf collimators (MLC) positions reported in the log‐file were compared with plan positions using TG‐142 criteria.

Result

The γ (3%, 3 mm) for the 190 plans were found to be between the range of 72.5%–100.0% and 95.4%–100.0% for ArcCheck^®‐MR and LinacVeiw, respectively. All the delivered gantry angle and jaws were found to be within 0.2° and 2 mm. MLCs that were outside the guard leaves or under the diaphragms were found to have more than 1.0 mm discrepancy. This was attributed to the linac internal override for these MLCs and had no dosimetric impact. Excluding these discrepancies, all MLC positions were found to be within 1.0 mm. The γ (3%, 3 mm) for the 30 incorrectly matched fields were found to be 3.9%–84.8% and 0.1%–64.4% for ArcCheck^®‐MR and LinacVeiw, respectively.

Conclusion

Significant ranked correlation demonstrates the automated log‐file analysis can be used for PSQA and expedite the ART workflow. Ongoing PSQA will be compared with log‐file analysis to investigate the longer term reproducibility and correlation.

Automatic 3D Monte-Carlo-based secondary dose calculation for online verification of 1.5 T magnetic resonance imaging guided radiotherapy

•

First implementation of an independent 3D-secondary dose calculation (3D-SDC).

•

Validation of the 3D-SDC solution using patient plans and experimental plan QA.

•

Online SDC of central targets is feasible with a median calculation time of 1:23 min.

•

Peripheral targets with small beam numbers need alternative validation strategies.

Background and purpose

Hybrid magnetic resonance linear accelerator (MR-Linac) systems represent a novel technology for online adaptive radiotherapy. 3D secondary dose calculation (SDC) of online adapted plans is required to assure patient safety. Currently, no 3D-SDC solution is available for 1.5T MR-Linac systems. Therefore, the aim of this project was to develop and validate a method for online automatic 3D-SDC for adaptive MR-Linac treatments.

Materials and methods

An accelerator head model was designed for an 1.5T MR-Linac system, neglecting the magnetic field. The use of this model for online 3D-SDC of MR-Linac plans was validated in a three-step process: (1) comparison to measured beam data, (2) investigation of performance and limitations in a planning phantom and (3) clinical validation using n = 100 patient plans from different tumor entities, comparing the developed 3D-SDC with experimental plan QA.

Results

The developed model showed median gamma passing rates compared to MR-Linac base data of 84.7%, 100% and 99.1% for crossplane, inplane and depth-dose-profiles, respectively. Comparison of 3D-SDC and full dose calculation in a planning phantom revealed that with ⩾5 beams gamma passing rates >95% can be achieved for central target locations. With a median calculation time of 1:23 min, 3D-SDC of online adapted clinical MR-Linac plans demonstrated a median gamma passing rate of 98.9% compared to full dose calculation, whereas experimental plan QA reached 99.5%.

Conclusion

Here, we describe the first technical 3D-SDC solution for online adaptive MR-guided radiotherapy. For clinical situations with peripheral targets and a small number of beams additional verification appears necessary. Further improvement may include 3D-SDC with consideration of the magnetic field.

ESTRO-ACROP recommendations on the clinical implementation of hybrid MR-linac systems in radiation oncology

Online magnetic resonance-guided radiotherapy (oMRgRT) represents one of the most innovative applications of current image-guided radiation therapy (IGRT). The revolutionary concept of oMRgRT systems is the ability to acquire MR images for adaptive treatment planning and also online imaging during treatment delivery. The daily adaptive planning strategies allow to improve targeting accuracy while avoiding critical structures. This ESTRO-ACROP recommendation aims to provide an overview of available systems and guidance for best practice in the implementation phase of hybrid MR-linac systems. Unlike the implementation of other radiotherapy techniques, oMRgRT adds the MR environment to the daily practice of radiotherapy, which might be a new experience for many centers. New issues and challenges that need to be thoroughly explored before starting clinical treatments will be highlighted.

MR-Guided Radiotherapy for Liver Malignancies

MR guided radiotherapy represents one of the most promising recent technological innovations in the field. The possibility to better visualize therapy volumes, coupled with the innovative online adaptive radiotherapy and motion management approaches, paves the way to more efficient treatment delivery and may be translated in better clinical outcomes both in terms of response and reduced toxicity. The aim of this review is to present the existing evidence about MRgRT applications for liver malignancies, discussing the potential clinical advantages and the current pitfalls of this new technology.

Automatic dosimetric verification of online adapted plans on the Unity MR-Linac using 3D EPID dosimetry

BACKGROUND AND PURPOSE

The Unity MR-Linac is equipped with an EPID, the images from which contain information about the dose delivered to the patient. The purpose of this study was to introduce a framework for the automatic dosimetric verification of online adapted plans using 3D EPID dosimetry and to present the obtained dosimetric results.

MATERIALS AND METHODS

The framework was active during the delivery of 1207 online adapted plans corresponding to 127 clinical IMRT treatments (74 prostate, 19 rectum, 19 liver and 15 lymph node oligometastases). EPID reconstructed dose distributions in the patient geometry were calculated automatically and then compared to the dose distributions calculated online by the treatment planning system (TPS). The comparison was performed by γ-analysis (3% global/2mm/10% threshold) and by the difference in median dose to the high-dose volume (ΔHDV_D50). 85% for γ-pass rate and 5% for ΔHDV_D50 were used as tolerance limit values.

RESULTS

93% of the online plans were verified automatically by the framework. Missing EPID data was the reason for automation failure. 91% of the verified plans were within tolerance.

CONCLUSION

Automatic dosimetric verification of online adapted plans on the Unity MR-Linac is feasible using in vivo 3D EPID dosimetry. Almost all online adapted plans were approved automatically by the framework. This newly developed framework is a major step forward towards the clinical implementation of a permanent safety net for the entire online adaptive workflow.

Feasibility of MR-guided ultrahypofractionated radiotherapy in 5, 2 or 1 fractions for prostate cancer

The drive towards hypofractionated prostate radiotherapy is motivated by a low alpha/beta ratio for prostate cancer (1 to 3 Gy) compared to surrounding organs at risk, implying an improved therapeutic ratio with increasing dose per fraction. Early evidence from studies of ultrahypofractionated (UHF) prostate HDR brachytherapy has shown good tolerability in terms of normal tissue toxicities and clinical outcomes similar to conventional fractionation schedules. MR-guided stereotactic body radiotherapy (SBRT) with online plan adaptation and real-time tumour imaging may enable UHF doses to be delivered to the prostate safely, without the invasiveness of brachytherapy. The feasibility of UHF prostate treatment planning for the Unity MR-Linac (MRL, Elekta AB, Stockholm) was investigated for target prescriptions and planning constraints derived from the HDR brachytherapy and SBRT literature. Monaco 5.40 (Elekta) was used to generate MRL step-and-shoot IMRT plans for three dose fractionation protocols (5, 2 and 1 fractions), for ten randomly selected previously treated prostate cancer patients. Of the ten plans per UHF scheme, all clinical goals were met in all cases for 5 fractions, and in six cases for both 2 and 1 fraction schemes. PTV D95% was compromised by up to 6.4% and 3.9% of the associated target dose for 2 and 1 fraction plans respectively. There were two cases of PTV D95% compromise greater than a 5% dose decrease for the 2 fraction plans. The study suggests feasibility of the UHF treatment planning approaches if combined with real-time motion mitigation strategies.

Consistent and invertible deformation vector fields for a breathing anthropomorphic phantom: a post-processing framework for the XCAT phantom

Breathing motion is challenging for radiotherapy planning and delivery. This requires advanced four-dimensional (4D) imaging and motion mitigation strategies and associated validation tools with known deformations. Numerical phantoms such as the XCAT provide reproducible and realistic data for simulation-based validation. However, the XCAT generates partially inconsistent and non-invertible deformations where tumours remain rigid and structures can move through each other. We address these limitations by post-processing the XCAT deformation vector fields (DVF) to generate a breathing phantom with realistic motion and quantifiable deformation. An open-source post-processing framework was developed that corrects and inverts the XCAT-DVFs while preserving sliding motion between organs. Those post-processed DVFs are used to warp the first XCAT-generated image to consecutive time points providing a 4D phantom with a tumour that moves consistently with the anatomy, the ability to scale lung density as well as consistent and invertible DVFs. For a regularly breathing case, the inverse consistency of the DVFs was verified and the tumour motion was compared to the original XCAT. The generated phantom and DVFs were used to validate a motion-including dose reconstruction (MIDR) method using isocenter shifts to emulate rigid motion. Differences between the reconstructed doses with and without lung density scaling were evaluated. The post-processing framework produced DVFs with a maximum [Formula: see text]-percentile inverse-consistency error of 0.02 mm. The generated phantom preserved the dominant sliding motion between the chest wall and inner organs. The tumour of the original XCAT phantom preserved its trajectory while deforming consistently with the underlying tissue. The MIDR was compared to the ground truth dose reconstruction illustrating its limitations. MIDR with and without lung density scaling resulted in small dose differences up to 1 Gy (prescription 54 Gy). The proposed open-source post-processing framework overcomes important limitations of the original XCAT phantom and makes it applicable to a wider range of validation applications within radiotherapy.

Prostate intrafraction motion during the preparation and delivery of MR-guided radiotherapy sessions on a 1.5T MR-Linac

PURPOSE

To evaluate prostate intrafraction motion using MRI during the full course of online adaptive MR-Linac radiotherapy (RT) fractions, in preparation of MR-guided extremely hypofractionated RT.

MATERIAL AND METHODS

Five low and intermediate risk prostate cancer patients were treated with 20 × 3.1 Gy fractions on a 1.5T MR-Linac. Each fraction, initial MRI (Pre) scans were obtained at the start of every treatment session. Pre-treatment planning MRI contours were propagated and adapted to this Pre scan after which plan re-optimization was started in the treatment planning system followed by dose delivery. 3D Cine-MR imaging was started simultaneously with beam-on and acquired over the full beam-on period. Prostate intrafraction motion in this cine-MR was determined with a previously validated soft-tissue contrast based tracking algorithm. In addition, absolute accuracy of the method was determined using a 4D phantom.

RESULTS

Prostate motion was completely automatically determined over the full on-couch period (approx. 45 min) with no identified mis-registrations. The translation 95% confidence intervals are within clinically applied margins of 5 mm, and plan adaption for intrafraction motion was required in only 4 out of 100 fractions.

CONCLUSION

This is the first study to investigate prostate intrafraction motions during entire MR-guided RT sessions on an MR-Linac. We have shown that high quality 3D cine-MR imaging and prostate tracking during RT is feasible with beam-on. The clinically applied margins of 5 mm have proven to be sufficient for these treatments and may potentially be further reduced using intrafraction plan adaptation guided by cine-MR imaging.

Delivered dose quantification in prostate radiotherapy using online 3D cine imaging and treatment log files on a combined 1.5T magnetic resonance imaging and linear accelerator system

BACKGROUND AND PURPOSE

Monitoring the intrafraction motion and its impact on the planned dose distribution is of crucial importance in radiotherapy. In this work we quantify the delivered dose for the first prostate patients treated on a combined 1.5T Magnetic Resonance Imaging (MRI) and linear accelerator system in our clinic based on online 3D cine-MR and treatment log files.

MATERIALS AND METHODS

A prostate intrafraction motion trace was obtained with a soft-tissue based rigid registration method with six degrees of freedom from 3D cine-MR dynamics with a temporal resolution of 8.5-16.9 s. For each fraction, all dynamics were also registered to the daily MR image used during the online treatment planning, enabling the mapping to this reference point. Moreover, each fraction's treatment log file was used to extract the timestamped machine parameters during delivery and assign it to the appropriate dynamic volume. These partial plans to dynamic volume combinations were calculated and summed to yield the delivered fraction dose. The planned and delivered dose distributions were compared among all patients for a total of 100 fractions.

RESULTS

The clinical target volume underwent on average a decrease of 2.2% ± 2.9% in terms of D99% coverage while bladder V62Gy was increased by 1.6% ± 2.3% and rectum V62Gy decreased by 0.2% ± 2.2%.

CONCLUSIONS

The first MR-linac dose reconstruction results based on prostate tracking from intrafraction 3D cine-MR and treatment log files are presented. Such a pipeline is essential for online adaptation especially as we progress to MRI-guided extremely hypofractionated treatments.

Dosimetric effect of intrafraction motion and different localization strategies in prostate SBRT

The aim of this study was to evaluate the dosimetric effect of continuous motion monitoring based localization (Calypso, Varian Medical Systems), gating and intrafraction motion correction in prostate SBRT. Delivered doses were modelled by reconstructing motion inclusive dose distributions for different localization strategies. Actually delivered dose (strategy A) utilized initial Calypso localization, CBCT and additional pre-treatment motion correction by kV-imaging and Calypso, and gating during the irradiation. The effect of gating was investigated by simulating non-gated treatments (strategy B). Additionally, non-gated and single image-guided (CBCT) localization was simulated (strategy C). A total of 308 fractions from 22 patients were reconstructed. The dosimetric effect was evaluated by comparing motion inclusive target and risk organ dose-volume parameters to planned values. Motion induced dose deficits were seen mainly in PTV and CTV to PTV margin regions, whereas CTV dose deficits were small in all strategies: mean ± SD difference in CTVD99% was -0.3 ± 0.4%, -0.4 ± 0.6% and -0.7 ± 1.2% in strategies A, B and C, respectively. Largest dose deficits were seen in individual fractions for strategy C (maximum dose reductions were -29.0% and -7.1% for PTVD95% and CTVD99%, respectively). The benefit of gating was minor, if additional motion correction was applied immediately prior to irradiation. Continuous motion monitoring based localization and motion correction ensured the target coverage and minimized the OAR exposure for every fraction and is recommended to use in prostate SBRT. The study is part of clinical trial NCT02319239.