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

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.

Ultra-hypofractionated prostate cancer radiotherapy: Dosimetric impact of real-time intrafraction prostate motion and daily anatomical changes

PURPOSE

To retrospectively assess the differences between planned and delivered dose during ultra-hypofractionated (UHF) prostate cancer treatments, by evaluating the dosimetric impact of daily anatomical variations alone, and in combination with prostate intrafraction motion.

METHODS

Prostate intrafraction motion was recorded with a transperineal ultrasound probe in 15 patients treated by UHF radiotherapy (36.25 Gy/5 fractions). The dosimetric objective was to cover 99 % of the clinical target volume with the 100 % prescription isodose line. After treatment, planning CT (pCT) images were deformably registered onto daily Cone Beam CT to generate pseudo-CT for dose accumulation (accumulated CT, aCT). The interplay effect was accounted by synchronizing prostatic shifts and beam geometry. Finally, the shifted dose maps were accumulated (moved-accumulated CT, maCT).

RESULTS

No significant change in daily CTV volumes was observed. Conversely, CTV V100% was 98.2 ± 0.8 % and 94.7 ± 2.6 % on aCT and maCT, respectively, compared with 99.5 ± 0.2 % on pCT (p < 0.0001). Bladder volume was smaller than planned in 76 % of fractions and D5cc was 33.8 ± 3.2 Gy and 34.4 ± 3.4 Gy on aCT (p = 0.02) and maCT (p = 0.01) compared with the pCT (36.0 ± 1.1 Gy). The rectum was smaller than planned in 50.3 % of fractions, but the dosimetric differences were not statistically significant, except for D1cc, found smaller on the maCT (33.2 ± 3.2 Gy, p = 0.02) compared with the pCT (35.3 ± 0.7 Gy).

CONCLUSIONS

Anatomical variations and prostate movements had more important dosimetric impact than anatomical variations alone, although, in some cases, the two phenomena compensated. Therefore, an efficient IGRT protocol is required for treatment implementation to reduce setup errors and control intrafraction motion.

Impact of daily plan adaptation on accumulated doses in ultra-hypofractionated magnetic resonance-guided radiation therapy of prostate cancer

Background and purpose

Ultra-hypofractionated online adaptive magnetic resonance-guided radiotherapy (MRgRT) is promising for prostate cancer. However, the impact of online adaptation on target coverage and organ-at-risk (OAR) sparing at the level of accumulated dose has not yet been reported. Using deformable image registration (DIR)-based accumulation, we compared the delivered adapted dose with the simulated non-adapted dose.

Materials and methods

Twenty-three prostate cancer patients treated at two clinics with 0.35 T magnetic resonance-guided linear accelerator (MR-linac) following the same treatment protocol (5 ×  7.5 Gy with urethral sparing and daily adaptation) were included. The fraction MR images were deformably registered to the planning MR image. Both non-adapted and adapted fraction doses were accumulated with the corresponding vector fields. Two DIR approaches were implemented. PTV* (planning target volume minus urethra+2mm) D95%, CTV* (clinical target volume minus urethra) D98%, and OARs (urethra+2mm, bladder, and rectum) D0.2cc, were evaluated. Statistical significance was inferred from a two-tailed Wilcoxon signed-rank test (p < 0.05).

Results

Normalized to the baseline, the accumulated PTV* D95% increased significantly by 2.7 % ([1.5, 4.3]%) through adaptation, and the CTV* D98% by 1.2 % ([0.1, 1.7]%). For the OARs after adaptation, accumulated bladder D0.2cc decreased by 0.4 % ([-1.2, 0.4]%), urethra+2mmD0.2cc by 0.8 % ([-1.6, -0.1]%), while rectum D0.2cc increased by 2.6 % ([1.2, 4.9]%). For all patients, rectum D0.2cc was still below the clinical constraint. Results of both DIR approaches differed on average by less than 0.2 %.

Conclusions

Online adaptation in MRgRT improved target coverage and OARs sparing at the level of accumulated dose.

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 Use of MR-Guided Radiation Therapy for Head and Neck Cancer and Recommended Reporting Guidance

Although magnetic resonance imaging (MRI) has become standard diagnostic workup for head and neck malignancies and is currently recommended by most radiological societies for pharyngeal and oral carcinomas, its utilization in radiotherapy has been heterogeneous during the last decades. However, few would argue that implementing MRI for annotation of target volumes and organs at risk provides several advantages, so that implementation of the modality for this purpose is widely accepted. Today, the term MR-guidance has received a much broader meaning, including MRI for adaptive treatments, MR-gating and tracking during radiotherapy application, MR-features as biomarkers and finally MR-only workflows. First studies on treatment of head and neck cancer on commercially available dedicated hybrid-platforms (MR-linacs), with distinct common features but also differences amongst them, have also been recently reported, as well as "biological adaptation" based on evaluation of early treatment response via functional MRI-sequences such as diffusion weighted ones. Yet, all of these approaches towards head and neck treatment remain at their infancy, especially when compared to other radiotherapy indications. Moreover, the lack of standardization for reporting MR-guided radiotherapy is a major obstacle both to further progress in the field and to conduct and compare clinical trials. Goals of this article is to present and explain all different aspects of MR-guidance for radiotherapy of head and neck cancer, summarize evidence, as well as possible advantages and challenges of the method and finally provide a comprehensive reporting guidance for use in clinical routine and trials.

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.

Gating and intrafraction drift correction on a 1.5 T MR-Linac: Clinical dosimetric benefits for upper abdominal tumors

This work reports on the first seven patients treated with gating and baseline drift correction on the high-field MR-Linac system. Dosimetric analysis showed that the active motion management system improved congruence to the planned dose, efficiently mitigating detrimental effects of intrafraction motion in the upper abdomen.

Evaluation of Unity 1.5 T MR-linac plan quality in patients with prostate cancer

The Unity magnetic resonance (MR) linear accelerator (MRL) with MR-guided adaptive radiotherapy (MRgART) is capable of online MRgART where images are acquired on the treatment day and the radiation treatment plan is immediately replanned and performed. We evaluated the MRgART plan quality and plan reproducibility of the Unity MRL in patients with prostate cancer. There were five low- or moderate-risk and five high-risk patients who received 36.25 Gy or 40 Gy, respectively in five fractions. All patients underwent simulation magnetic resonance imaging (MRI) and five online adaptive MRI. We created plans for 5, 7, 9, 16, and 20 beams and for 60, 100, and 150 segments. We evaluated the target and organ doses for different number of beams and segments, respectively. Variation in dose constraint between the simulation plan and online adaptive plan was measured for each patient to assess plan reproducibility. The plan quality improved with the increasing number of beams. However, the proportion of significantly improved dose constraints decreased as the number of beams increased. For some dose parameters, there were statistically significant differences between 60 and 100 segments, and 100 and 150 segments. The plan of five beams exhibited limited reproducibility. The number of segments had minimal impact on plan reproducibility, but 60 segments sometimes failed to meet dose constraints for online adaptive plan. The optimization and delivery time increased with the number of beams and segments. We do not recommend using five or fewer beams for a reproducible and high-quality plan in the Unity MRL. In addition, many number of segments and beams may help meet dose constraints during online adaptive plan. Treatment with the Unity MRL should be performed with the appropriate number of beams and segments to achieve a good balance among plan quality, delivery time, and optimization time.

Experimental validation of multi-fraction online adaptations in magnetic resonance guided radiotherapy

Background and purpose

Radiotherapy plan verification is generally performed on the reference plan based on the pre-treatment anatomy. However, the introduction of online adaptive treatments demands a new approach, as plans are created daily on different anatomies. The aim of this study was to experimentally validate the accuracy of total doses of multi-fraction plan adaptations in magnetic resonance imaging guided radiotherapy in a phantom study, isolated from the uncertainty of deformable image registration.

Materials and methods

We experimentally verified the total dose, measured on external beam therapy 3 (EBT3) film, using a treatment with five online adapted fractions. Three series of experiments were performed, each focusing on a category of inter-fractional variation; translations, rotations and body modifications. Variations were introduced during each fraction and adapted plans were generated and irradiated. Single fraction doses and total doses over five online adapted fractions were investigated.

Results

The online adapted measurements and calculations showed a good agreement for single fractions and multi-fraction treatments for the dose profiles, gamma passing rates, dose deviations and distances to agreement. The gamma passing rate using a 2%/2 mm criterion ranged from 99.2% to 99.5% for a threshold dose of 10% of the maximum dose (Dmax) and from 96.2% to 100% for a threshold dose of 90% of Dmax, for the total translations, rotations and body modifications.

Conclusions

The total doses of multi-fraction treatments showed similar accuracies compared to single fraction treatments, indicating an accurate dosimetric outcome of a multi-fraction treatment in adaptive magnetic resonance imaging guided radiotherapy.

Gaussian Processes for real-time 3D motion and uncertainty estimation during MR-guided radiotherapy

Respiratory motion during radiotherapy causes uncertainty in the tumor's location, which is typically addressed by an increased radiation area and a decreased dose. As a result, the treatments' efficacy is reduced. The recently proposed hybrid MR-linac scanner holds the promise to efficiently deal with such respiratory motion through real-time adaptive MR-guided radiotherapy (MRgRT). For MRgRT, motion-fields should be estimated from MR-data and the radiotherapy plan should be adapted in real-time according to the estimated motion-fields. All of this should be performed with a total latency of maximally 200 ms, including data acquisition and reconstruction. A measure of confidence in such estimated motion-fields is highly desirable, for instance to ensure the patient's safety in case of unexpected and undesirable motion. In this work, we propose a framework based on Gaussian Processes to infer 3D motion-fields and uncertainty maps in real-time from only three readouts of MR-data. We demonstrated an inference frame rate up to 69 Hz including data acquisition and reconstruction, thereby exploiting the limited amount of required MR-data. Additionally, we designed a rejection criterion based on the motion-field uncertainty maps to demonstrate the framework's potential for quality assurance. The framework was validated in silico and in vivo on healthy volunteer data (n=5) acquired using an MR-linac, thereby taking into account different breathing patterns and controlled bulk motion. Results indicate end-point-errors with a 75th percentile below 1 mm in silico, and a correct detection of erroneous motion estimates with the rejection criterion. Altogether, the results show the potential of the framework for application in real-time MR-guided radiotherapy with an MR-linac.

Applicability and usage of dose mapping/accumulation in radiotherapy

Dose mapping/accumulation (DMA) is a topic in radiotherapy (RT) for years, but has not yet found its widespread way into clinical RT routine. During the ESTRO Physics workshop 2021 on "commissioning and quality assurance of deformable image registration (DIR) for current and future RT applications", we built a working group on DMA from which we present the results of our discussions in this article. Our aim in this manuscript is to shed light on the current situation of DMA in RT and to highlight the issues that hinder consciously integrating it into clinical RT routine. As a first outcome of our discussions, we present a scheme where representative RT use cases are positioned, considering expected anatomical variations and the impact of dose mapping uncertainties on patient safety, which we have named the DMA landscape (DMAL). This tool is useful for future reference when DMA applications get closer to clinical day-to-day use. Secondly, we discussed current challenges, lightly touching on first-order effects (related to the impact of DIR uncertainties in dose mapping), and focusing in detail on second-order effects often dismissed in the current literature (as resampling and interpolation, quality assurance considerations, and radiobiological issues). Finally, we developed recommendations, and guidelines for vendors and users. Our main point include: Strive for context-driven DIR (by considering their impact on clinical decisions/judgements) rather than perfect DIR; be conscious of the limitations of the implemented DIR algorithm; and consider when dose mapping (with properly quantified uncertainties) is a better alternative than no mapping.

A new target localization method for image-guided radiation therapy of prostate cancer

In imaged-guided radiation therapy (IGRT), target localization is usually done with rigid-body registration based on anatomy matching. Problems arise when the target volume can only be matched partially due to inter-fractional organ motion and deformation, resulting in deteriorated target coverage and critical structure sparing. A new target localization method is investigated in which the treatment target volume is aligned with the prescription isodose surface. Our study included 15 prostate patients previously treated with intensity-modulated radiation therapy (IMRT). Patient setup and target localization were performed using a CT-on-rails system before and after the IMRT treatment. IMRT plans were generated on the original simulation CTs (15) and the same MUs and leaf sequences were used to compute the dose distributions on post-treatment CTs (98) with the isocenter adjustments based on either anatomical structure matching or prescription isodose surface alignment. When patients were aligned with the traditional anatomy matching method, the dose to 95% of the CTV, D₉₅, received 74.0 - 77.6 Gy and the minimum CTV dose, D_min, was 61.9 - 71.6 Gy, respectively, in the cumulative dose distributions. The rectal dose-volume constraints were violated in 35.7% of the treatment fractions. When patients were aligned using the new localization method, the dose to 95% of the CTV, D₉₅, received 74.0 - 78.2 Gy and the minimum CTV dose, D_min, was 68.4 - 71.6 Gy, respectively, in the cumulative dose distributions. The rectal dose-volume constraints were violated in 17.3% of the treatment fractions. Traditional IGRT target localization based on anatomy matching is effective for population-based PTV margins but not ideal for those patients with large inter-fractional prostate rotation/deformation due to large rectal and bladder volume variation. The new method using the prescription isodose surface to align the target volume could improve the target coverage and rectal sparing for these patients, which can be implemented clinically to improve target dose delivery accuracy.

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.

Magnetic Resonance-Guided Adaptive Radiation therapy for Prostate Cancer: The First Results from the MOMENTUM study-An International Registry for the Evidence-Based Introduction of Magnetic Resonance-Guided Adaptive Radiation Therapy

PURPOSE

Magnetic resonance (MR)-guided radiation therapy (MRgRT) is a new technique for treatment of localized prostate cancer (PCa). We report the 12-month outcomes for the first PCa patients treated within an international consortium (the MOMENTUM study) on a 1.5T MR-Linac system with ultrahypofractionated radiation therapy.

METHODS AND MATERIALS

Patients treated with 5 × 7.25 Gy were identified. Prostate specific antigen-level, physician-reported toxicity (Common Terminology Criteria for Adverse Events [CTCAE]), and patient-reported outcomes (Quality of Life Questionnaire PR25 and Quality of Life Questionnaire C30 questionnaires) were recorded at baseline and at 3, 6, and 12 months of follow-up (FU). Pairwise comparative statistics were conducted to compare outcomes between baseline and FU.

RESULTS

The study included 425 patients with localized PCa (11.4% low, 82.0% intermediate, and 6.6% high-risk), and 365, 313, and 186 patients reached 3-, 6-, and 12-months FU, respectively. Median prostate specific antigen level declined significantly to 1.2 ng/mL and 0.1 ng/mL at 12 months FU for the nonandrogen deprivation therapy (ADT) and ADT group, respectively. The peak of genitourinary and gastrointestinal CTCAE toxicity was reported at 3 months FU, with 18.7% and 1.7% grade ≥2, respectively. The QLQ-PR25 questionnaire outcomes showed significant deterioration in urinary domain score at all FU moments, from 8.3 (interquartile range [IQR], 4.1-16.6) at baseline to 12.4 (IQR, 8.3-24.8; P = .005) at 3 months, 12.4 (IQR, 8.3-20.8; P = .018;) at 6 months, and 12.4 (IQR, 8.3-20.8; P = .001) at 12 months. For the non-ADT group, physician- and patient-reported erectile function worsened significantly between baseline and 12 months FU.

CONCLUSIONS

Ultrahypofractionated MR-guided radiation therapy for localized PCa using a 1.5T MR-Linac is effective and safe. The peak of CTCAE genitourinary and gastrointestinal toxicity was reported at 3 months FU. Furthermore, for patients without ADT, a significant increase in CTCAE erectile dysfunction was reported at 12 months FU. These data are useful for educating patients on expected outcomes and informing study design of future comparative-effectiveness studies.

Daily online contouring and re-planning versus translation-only correction in neurovascular-sparing magnetic resonance-guided radiotherapy for localized prostate cancer

Neurovascular bundle (NVB) and internal pudendal artery (IPA) sparing during magnetic resonance-guided radiotherapy (MRgRT) for prostate cancer aims for preservation of erectile function. Our present workflow involves daily online contouring and re-planning on a 1.5 T MR-linac, as alternative to conventional (rigid) translation-only corrections of the prostate. We compared planned dose for the NVB and IPA between strategies. Total planned dose was significantly lower with daily online contouring and re-planning for the NVB, but not for the IPA. For the NVB and IPA, the intrapatient difference between highest and lowest fraction dose was significantly smaller for the contouring and re-planning plans.

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.

Dosimetric impact of intrafraction motion under abdominal compression during MR-guided SBRT for (Peri-) pancreatic tumors

Objective. Intrafraction motion is a major concern for the safety and effectiveness of high dose stereotactic body radiotherapy (SBRT) in the upper abdomen. In this study, the impact of the intrafraction motion on the delivered dose was assessed in a patient group that underwent MR-guided radiotherapy for upper abdominal malignancies with an abdominal corset. Approach. Fast online 2D cine MRI was used to extract tumor motion during beam-on time. These tumor motion profiles were combined with linac log files to reconstruct the delivered dose in 89 fractions of MR-guided SBRT in twenty patients. Aside the measured tumor motion, motion profiles were also simulated for a wide range of respiratory amplitudes and drifts, and their subsequent dosimetric impact was calculated in every fraction. Main results. The average (SD) D 99% of the gross tumor volume (GTV), relative to the planned D 99%, was 0.98 (0.03). The average (SD) relative D 0.5cc of the duodenum, small bowel and stomach was 0.99 (0.03), 1.00 (0.03), and 0.97 (0.05), respectively. No correlation of respiratory amplitude with dosimetric impact was observed. Fractions with larger baseline drifts generally led to a larger uncertainty of dosimetric impact on the GTV and organs at risk (OAR). The simulations yielded that the delivered dose is highly dependent on the direction of on baseline drift. Especially in anatomies where the OARs are closely abutting the GTV, even modest LR or AP drifts can lead to substantial deviations from the planned dose. Significance. The vast majority of the fractions was only modestly impacted by intrafraction motion, increasing our confidence that MR-guided SBRT with abdominal compression can be safely executed for patients with abdominal tumors, without the use of gating or tracking strategies.

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.

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.

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 Adaptive Radiotherapy for OAR Sparing in Head and Neck Cancers

Simple Summary

Normal tissue toxicities in head and neck cancer persist as a cause of decreased quality of life and are associated with poorer treatment outcomes. The aim of this article is to review organ at risk (OAR) sparing approaches available in MR-guided adaptive radiotherapy and present future developments which hope to improve treatment outcomes. Increasing the spatial conformity of dose distributions in radiotherapy is an important first step in reducing normal tissue toxicities, and MR-guided treatment devices presents a new opportunity to use biological information to drive treatment decisions on a personalized basis.

Abstract

MR-linac devices offer the potential for advancements in radiotherapy (RT) treatment of head and neck cancer (HNC) by using daily MR imaging performed at the time and setup of treatment delivery. This article aims to present a review of current adaptive RT (ART) methods on MR-Linac devices directed towards the sparing of organs at risk (OAR) and a view of future adaptive techniques seeking to improve the therapeutic ratio. This ratio expresses the relationship between the probability of tumor control and the probability of normal tissue damage and is thus an important conceptual metric of success in the sparing of OARs. Increasing spatial conformity of dose distributions to target volume and OARs is an initial step in achieving therapeutic improvements, followed by the use of imaging and clinical biomarkers to inform the clinical decision-making process in an ART paradigm. Pre-clinical and clinical findings support the incorporation of biomarkers into ART protocols and investment into further research to explore imaging biomarkers by taking advantage of the daily MR imaging workflow. A coherent understanding of this road map for RT in HNC is critical for directing future research efforts related to sparing OARs using image-guided radiotherapy (IGRT).

Experimental investigation of dynamic real-time rotation-including dose reconstruction during prostate tracking radiotherapy

BACKGROUND

Hypofractionation in prostate radiotherapy is of increasing interest. Steep dose gradients and a large weight on each individual fraction emphasize the need for motion management. Real-time motion management techniques such as multi-leaf collimator (MLC) tracking or couch tracking typically adjust for translational motion while rotations remain uncompensated with unknown dosimetric impact.

PURPOSE

The purpose of this study is to demonstrate and validate dynamic real-time rotation-including dose reconstruction during radiotherapy experiments with and without MLC and couch tracking.

METHODS

Real-time dose reconstruction was performed using the in-house developed software DoseTracker. DoseTracker receives streamed target positions and accelerator parameters during treatment delivery and uses a pencil beam algorithm with water density assumption to reconstruct the dose in a moving target. DoseTracker's ability to reconstruct motion-induced dose errors in a dynamically rotating and translating target was investigated during three different scenarios: (1) no motion compensation and translational motion correction with (2) MLC tracking and (3) couch tracking. In each scenario, dose reconstruction was performed online and in real-time during delivery of two dual-arc volumetric modulated arc therapy (VMAT) prostate plans with a prescribed fraction dose of 7 Gy to the prostate and simultaneous intraprostatic lesion boosts with doses of at least 8 Gy, but up to 10 Gy as long as the organs-at-risk dose constraints were fulfilled. The plans were delivered to a pelvis phantom that replicated three patient-measured motion traces using a rotational insert with 21 layers of EBT3 film spaced 2.5 mm apart. DoseTracker repeatedly calculated the actual motion-including dose increment and the planned static dose increment since the last calculation in 84500 points in the film stack. The experiments were performed with a TrueBeam accelerator with MLC and couch tracking based on electromagnetic transponders embedded in the film stack. The motion-induced dose error was quantified as the difference between the final cumulative dose with motion and without motion using the 2D 2%/2mm γ-failure rate and the difference in dose to 95% of the clinical target volume (CTV ΔD_95% ) and the gross target volume (GTV ΔD_95% ) as well as the difference in dose to 0.1 cm³ of the urethra, bladder, and rectum (ΔD_0.1CC ). The motion-induced errors were compared between dose reconstructions and film measurements.

RESULTS

The dose was reconstructed in all calculation points at a mean frequency of 4.7 Hz. The root-mean-square difference between real-time reconstructed and film measured motion-induced errors was 3.1%-points (γ-failure rate), 0.13 Gy (CTV ΔD_95% ), 0.23 Gy (GTV ΔD_95% ), 0.19 Gy (urethra ΔD_0.1CC ), 0.09 Gy (bladder ΔD_0.1CC ), and 0.07 Gy (rectum ΔD_0.1CC ).

CONCLUSIONS

In a series of phantom experiments, online real-time rotation-including dose reconstruction was performed for the first time. The calculated motion-induced errors agreed well with film measurements. The dose reconstruction provides a valuable tool for monitoring dose delivery and investigating the efficacy of advanced motion-compensation techniques in the presence of translational and rotational motion. This article is protected by copyright. All rights reserved.

Minimum non-isotropic and asymmetric margins for taking into account intrafraction prostate motion during moderately hypofractionated radiotherapy

PURPOSE

To investigate the impact on dose distribution of intrafraction motion during moderate hypofractionated prostate cancer treatments and to estimate minimum non-isotropic and asymmetric (NI-AS) treatment margins taking motion into account.

METHODS

Prostate intrafraction 3D displacements were recorded with a transperineal ultrasound probe and were evaluated in 46 prostate cancer patients (876 fractions) treated by moderate hypofractionated radiation therapy (60 Gy in 20 fractions). For 18 patients (346 fractions), treatment plans were recomputed increasing CTV-to-PTV margins from 0 to 6 mm with an auto-planning optimization algorithm. Dose distribution was estimated using the voxel shifting method by displacing CTV structure according to the retrieved movements. Time-dependent margins were finally calculated using both van Herk's formula and the voxel shifting method.

RESULTS

Mean intrafraction prostate displacements observed were -0.02 ± 0.52 mm, 0.27 ± 0.78 mm and -0.43 ± 1.06 mm in left-right, supero-inferior and antero-posterior directions, respectively. The CTV dosimetric coverage increased with increased CTV-to-PTV margins but it decreased with time. Hence using van Herk's formula, after 7 min of treatment, a margin of 0.4 and 0.5 mm was needed in left and right, 1.5 and 0.7 mm in inferior and superior and 1.1 and 3.2 mm in anterior and posterior directions, respectively. Conversely, using the voxel shifting method, a margin of 0 mm was needed in left-right, 2 mm in superior, 3 mm in inferior and anterior and 5 mm in posterior directions, respectively. With this latter NI-AS margin strategy, the dosimetric target coverage was equivalent to the one obtained with a 5 mm homogeneous margin.

CONCLUSIONS

NI-AS margins would be required to optimally take into account intrafraction motion.

Injection of hydrogel spacer increased maximal intrafractional prostate motion in anterior and superior directions during volumetric modulated arc therapy-stereotactic body radiation therapy for prostate cancer

Background

The aim of this study was to clarify the association between intrafractional prostate shift and hydrogel spacer.

Methods

Thirty-eight patients who received definitive volumetric modulated arc therapy (VMAT)-stereotactic body radiation therapy (SBRT) for prostate cancer with prostate motion monitoring in our institution in 2018–2019 were retrospectively evaluated. In order to move the rectum away from the prostate, hydrogel spacer (SpaceOAR system, Boston Scientific, Marlborough, the United States) injection was proposed to the patients as an option in case of meeting the indication of use. We monitored intrafractional prostate motion by using a 4-dimensional (4D) transperineal ultrasound device: the Clarity 4D ultrasound system (Elekta AB). The deviation of the prostate was monitored in each direction: superior-inferior, left–right, and anterior–posterior. We also calculated the vector length. The maximum intrafractional displacement (MID) per fraction for each direction was detected and mean of MIDs was calculated per patient. The MIDs in the non-spacer group and the spacer group were compared using the unpaired t-test.

Results

We reviewed 33 fractions in eight patients as the spacer group and 148 fractions in 30 patients as the non-spacer group. The superior MID was 0.47 ± 0.07 (mean ± SE) mm versus 0.97 ± 0.24 mm (P = 0.014), the inferior MID was 1.07 ± 0.11 mm versus 1.03 ± 0.25 mm (P = 0.88), the left MID was 0.74 ± 0.08 mm versus 0.87 ± 0.27 mm (P = 0.55), the right MID was 0.67 ± 0.08 mm versus 0.92 ± 0.21 mm (P = 0.17), the anterior MID was 0.45 ± 0.06 mm versus 1.16 ± 0.35 mm (P = 0.0023), and the posterior MID was 1.57 ± 0.17 mm versus 1.37 ± 0.22 mm (P = 0.56) in the non-spacer group and the spacer group, respectively. The max of VL was 2.24 ± 0.19 mm versus 2.89 ± 0.62 mm (P = 0.19), respectively.

Conclusions

Our findings suggest that maximum intrafractional prostate motion during VMAT-SBRT was larger in patients with hydrogel spacer injection in the superior and anterior directions. Since this difference seemed not to disturb the dosimetric advantage of the hydrogel spacer, we do not recommend routine avoidance of the hydrogel spacer use.

Real-Time Non-Rigid 3D Respiratory Motion Estimation for MR-Guided Radiotherapy Using MR-MOTUS

The MR-Linac is a combination of an MR-scanner and radiotherapy linear accelerator (Linac) which holds the promise to increase the precision of radiotherapy treatments with MR-guided radiotherapy by monitoring motion during radiotherapy with MRI, and adjusting the radiotherapy plan accordingly. Optimal MR-guidance for respiratory motion during radiotherapy requires MR-based 3D motion estimation with a latency of 200-500 ms. Currently this is still challenging since typical methods rely on MR-images, and are therefore limited by the 3D MR-imaging latency. In this work, we present a method to perform non-rigid 3D respiratory motion estimation with 170 ms latency, including both acquisition and reconstruction. The proposed method called real-time low-rank MR-MOTUS reconstructs motion-fields directly from k -space data, and leverages an explicit low-rank decomposition of motion-fields to split the large scale 3D+t motion-field reconstruction problem posed in our previous work into two parts: (I) a medium-scale offline preparation phase and (II) a small-scale online inference phase which exploits the results of the offline phase for real-time computations. The method was validated on free-breathing data of five volunteers, acquired with a 1.5T Elekta Unity MR-Linac. Results show that the reconstructed 3D motion-field are anatomically plausible, highly correlated with a self-navigation motion surrogate ( R=0.975 ±0.0110 ), and can be reconstructed with a total latency of 170 ms that is sufficient for real-time MR-guided abdominal radiotherapy.

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.

[Primary treatment of prostate cancer using 1.5 T MR-linear accelerator]

Hintergrund

Der potenzielle Nutzen des verbesserten Weichteilkontrastes von MR-Sequenzen gegenüber der Computertomographie (CT) für die Radiotherapie des Prostatakarzinoms ist bekannt und führt zu konsistenteren und kleineren Zielvolumina sowie verbesserter Risikoorganschonung. Hybridgeräte aus Magnetresonanztomographie (MRT) und Linearbeschleuniger (MR-Linac) stellen eine neue vielversprechende Erweiterung der radioonkologischen Therapieoptionen dar.

Material und Methoden

Dieser Artikel gibt eine Übersicht über bisherige Erfahrungen, Indikationen, Vorteile und Herausforderungen für die Radiotherapie des primären Prostatakarzinoms mit dem 1,5-T-MR-Linac.

Ergebnisse

Alle strahlentherapeutischen Therapieindikationen für das primäre Prostatakarzinom können mit dem 1,5-T-MR-Linac abgedeckt werden. Die potenziellen Vorteile umfassen die tägliche MR-basierte Lagekontrolle in Bestrahlungsposition und die Möglichkeit der täglichen Echtzeitanpassung des Bestrahlungsplans an die aktuelle Anatomie der Beckenorgane (adaptive Strahlentherapie). Zusätzlich werden am 1,5-T-MR-Linac funktionelle MRT-Sequenzen für individuelles Response-Assessment für die Therapieanpassung untersucht. Dadurch soll das therapeutische Fenster weiter optimiert werden. Herausforderungen stellen u. a. die technische Komplexität und die Dauer der Behandlungssitzung dar.

Schlussfolgerung

Der 1,5-T-MR-Linac erweitert das radioonkologische Spektrum in der Therapie des Prostatakarzinoms und bietet Vorteile durch tagesaktuelle MRT-basierte Zielvolumendefinition und Planadaptation. Weitere klinische Untersuchungen sind notwendig, um die Patienten zu identifizieren, die von der Behandlung am MR-Linac gegenüber anderen strahlentherapeutischen Methoden besonders profitieren.

Technical Challenges of Real-Time Adaptive MR-Guided Radiotherapy

In the past few years, radiotherapy (RT) has experienced a major technological innovation with the development of hybrid machines combining magnetic resonance (MR) imaging and linear accelerators. This new technology for MR-guided cancer treatment has the potential to revolutionize the field of adaptive RT due to the opportunity to provide high-resolution, real-time MR imaging before and during treatment application. However, from a technical point of view, several challenges remain which need to be tackled to ensure safe and robust real-time adaptive MR-guided RT delivery. In this manuscript, several technical challenges to MR-guided RT are discussed. Starting with magnetic field strength tradeoffs, the potential and limitations for purely MR-based RT workflows are discussed. Furthermore, the current status of real-time 3D MR imaging and its potential for real-time RT are summarized. Finally, the potential of quantitative MR imaging for future biological RT adaptation is highlighted.

Quantitative investigation of dose accumulation errors from intra-fraction motion in MRgRT for prostate cancer

Accurate spatial dose delivery in radiotherapy is frequently complicated due to changes in the patient's internal anatomy during and in-between therapy segments. The recent introduction of hybrid MRI radiotherapy systems allows unequaled soft-tissue visualization during radiation delivery and can be used for dose reconstruction to quantify the impact of motion. To this end, knowledge of anatomical deformations obtained from continuous monitoring during treatment has to be combined with information on the spatio-temporal dose delivery to perform motion-compensated dose accumulation (MCDA). Here, the influence of the choice of deformable image registration algorithm, dose warping strategy, and magnetic resonance image resolution and signal-to-noise-ratio on the resulting MCDA is investigated. For a quantitative investigation, four 4D MRI-datasets representing typical patient observed motion patterns are generated using finite element modeling and serve as a gold standard. Energy delivery is simulated intra-fractionally in the deformed image space and, subsequently, MCDA-processed. Finally, the results are substantiated by comparing MCDA strategies on clinically acquired patient data. It is shown that MCDA is needed for correct quantitative dose reconstruction. For prostate treatments, using the energy per mass transfer dose warping strategy has the largest influence on decreasing dose estimation errors.

Proof-of-concept delivery of intensity modulated arc therapy on the Elekta Unity 1.5 T MR-linac

In this work we present the first delivery of intensity modulated arc therapy on the Elekta Unity 1.5 T MR-linac. The machine's current intensity modulated radiation therapy based control system was modified suitably to enable dynamic delivery of radiation, for the purpose of exploring MRI-guided radiation therapy adaptation modes in a research setting. The proof-of-concept feasibility was demonstrated by planning and delivering two types of plans, each investigating the performance of different parts of a dynamic treatment. A series of fixed-speed arc plans was used to show the high-speed capabilities of the gantry during radiation, while several fully modulated prostate plans-optimised following the volumetric modulated arc therapy approach-were delivered in order to establish the performance of its multi-leaf collimator and diaphragms. These plans were delivered to Delta⁴ Phantom+ MR and film phantoms passing the clinical quality assurance criteria used in our clinic. In addition, we also performed some initial MR imaging experiments during dynamic therapy, demonstrating that the impact of radiation and moving gantry/collimator components on the image quality is negligible. These results show that arc therapy is feasible on the Elekta Unity system. The machine's high performance components enable dynamic delivery during fast gantry rotation and can be controlled in a stable fashion to deliver fully modulated plans.

Single-fraction prostate stereotactic body radiotherapy: Dose reconstruction with electromagnetic intrafraction motion tracking

PURPOSE

To reconstruct the dose delivered during single-fraction urethra-sparing prostate stereotactic body radiotherapy (SBRT) accounting for intrafraction motion monitored by intraprostatic electromagnetic transponders (EMT).

METHODS

We analyzed data of 15 patients included in the phase I/II "ONE SHOT" trial and treated with a single fraction of 19 Gy to the planning target volume (PTV) and 17 Gy to the urethra planning risk volume. During delivery, prostate motion was tracked with implanted EMT. SBRT was interrupted when a 3-mm threshold was trespassed and corrected unless the offset was transient. Motion-encoded reconstructed (MER) plans were obtained by splitting the original plans into multiple sub-beams with isocenter shifts based on recorded EMT positions, mimicking prostate motion during treatment. We analyzed intrafraction motion and compared MER to planned doses.

RESULTS

The median EMT motion range (±SD) during delivery was 0.26 ± 0.09, 0.22 ± 0.14 and 0.18 ± 0.10 cm in the antero-posterior, supero-inferior, and left-right axes, respectively. Treatment interruptions were needed for 8 patients because of target motion beyond limits in the antero-posterior (n = 6) and/or supero-inferior directions (n = 4). Comparing MER vs. original plan there was a median relative dose difference of -1.9% (range, -7.9 to -1.0%) and of +0.5% (-0.3-1.7%) for PTV D_98% and D_2%, respectively. The clinical target volume remained sufficiently covered with a median D_98% difference of -0.3% (-1.6-0.5%). Bladder and rectum dosimetric parameters showed significant differences between original and MER plans, but mostly remained within acceptable limits.

CONCLUSIONS

The dosimetric impact of intrafraction prostate motion was minimal for target coverage for single-fraction prostate SBRT with real-time electromagnetic tracking combined with beam gating.

Quantitative magnetic resonance imaging on hybrid magnetic resonance linear accelerators: Perspective on technical and clinical validation

Many preclinical and clinical observations support that functional magnetic resonance imaging (MRI), such as diffusion weighted (DW) and dynamic contrast enhanced (DCE) MRI, might have a predictive value for radiotherapy. The aim of this review was to assess the current status of quantitative MRI on hybrid MR-Linacs. In a literature research, four publications were identified, investigating technical feasibility, accuracy, repeatability and reproducibility of DW and DCE-MRI in phantoms and first patients. Accuracy and short term repeatability was < 5% for DW-MRI in current MR-Linac systems. Consequently, quantitative imaging providing accurate and reproducible functional information seems possible in MR-Linacs.

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.