Recommendations for improved reproducibility of ADC derivation on behalf of the Elekta MRI-linac consortium image analysis working group

Background and purpose: The apparent diffusion coefficient (ADC), a potential imaging biomarker for radiotherapy response, needs to be reproducible before translation into clinical use. The aim of this study was to evaluate the multi-centre delineation- and calculation-related ADC variation and give recommendations to minimize it. Materials and methods: Nine centres received identical diffusion-weighted and anatomical magnetic resonance images of different cancerous tumours (adrenal gland, pelvic oligo metastasis, pancreas, and prostate). All centres delineated the gross tumour volume (GTV), clinical target volume (CTV), and viable tumour volume (VTV), and calculated ADCs using both their local calculation methods and each of the following calculation conditions: b-values 0–500 vs. 150–500 s/mm 2, region-of-interest (ROI)-based vs. voxel-based calculation, and mean vs. median. ADC variation was assessed using the mean coefficient of variation across delineations (CV D) and calculation methods (CV C). Absolute ADC differences between calculation conditions were evaluated using Friedman's test. Recommendations for ADC calculation were formulated based on observations and discussions within the Elekta MRI-linac consortium image analysis working group. Results: The median (range) CV D and CV C were 0.06 (0.02–0.32) and 0.17 (0.08–0.26), respectively. The ADC estimates differed 18% between b-value sets and 4% between ROI/voxel-based calculation (p-values < 0.01). No significant difference was observed between mean and median (p = 0.64). Aligning calculation conditions between centres reduced CV C to 0.04 (0.01–0.16). CV D was comparable between ROI types. Conclusion: Overall, calculation methods had a larger impact on ADC reproducibility compared to delineation. Based on the results, significant sources of variation were identified, which should be considered when initiating new studies, in particular multi-centre investigations

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

Accelerating multi-modal image registration using a supervoxel-based variational framework

For the successful completion of medical interventional procedures, several concepts, such as daily positioning compensation, dose accumulation or delineation propagation, rely on establishing a spatial coherence between planning images and images acquired at different time instants over the course of the therapy. To meet this need, image-based motion estimation and compensation relies on fast, automatic, accurate and precise registration algorithms. However, image registration quickly becomes a challenging and computationally intensive task, especially when multiple imaging modalities are involved. In the current study, a novel framework is introduced to reduce the computational overhead of variational registration methods. The proposed framework selects representative voxels of the registration process, based on a supervoxel algorithm. Costly calculations are hereby restrained to a subset of voxels, leading to a less expensive spatial regularized interpolation process. The novel framework is tested in conjunction with the recently proposed EVolution multi-modal registration method. This results in an algorithm requiring a low number of input parameters, is easily parallelizable and provides an elastic voxel-wise deformation with a subvoxel accuracy. The performance of the proposed accelerated registration method is evaluated on cross-contrast abdominal T1/T2 MR-scans undergoing a known deformation and annotated CT-images of the lung. We also analyze the ability of the method to capture slow physiological drifts during MR-guided high intensity focused ultrasound therapies and to perform multi-modal CT/MR registration in the abdomen. Results have shown that computation time can be reduced by 75% on the same hardware with no negative impact on the accuracy

Development and clinical introduction of automated radiotherapy treatment planning for prostate cancer

To develop an automated radiotherapy treatment planning and optimization workflow to efficiently create patient specifically optimized clinical grade treatment plans for prostate cancer and to implement it in clinical practice. A two-phased planning and optimization workflow was developed to automatically generate 77Gy 5-field simultaneously integrated boost intensity modulated radiation therapy (SIB-IMRT) plans for prostate cancer treatment. A retrospective planning study (n  =  100) was performed in which automatically and manually generated treatment plans were compared. A clinical pilot (n  =  21) was performed to investigate the usability of our method. Operator time for the planning process was reduced to  <5 min. The retrospective planning study showed that 98 plans met all clinical constraints. Significant improvements were made in the volume receiving 72Gy (V72Gy) for the bladder and rectum and the mean dose of the bladder and the body. A reduced plan variance was observed. During the clinical pilot 20 automatically generated plans met all constraints and 17 plans were selected for treatment. The automated radiotherapy treatment planning and optimization workflow is capable of efficiently generating patient specifically optimized and improved clinical grade plans. It has now been adopted as the current standard workflow in our clinic to generate treatment plans for prostate cancer

Development and clinical introduction of automated radiotherapy treatment planning for prostate cancer

To develop an automated radiotherapy treatment planning and optimization workflow to efficiently create patient specifically optimized clinical grade treatment plans for prostate cancer and to implement it in clinical practice. A two-phased planning and optimization workflow was developed to automatically generate 77Gy 5-field simultaneously integrated boost intensity modulated radiation therapy (SIB-IMRT) plans for prostate cancer treatment. A retrospective planning study (n  =  100) was performed in which automatically and manually generated treatment plans were compared. A clinical pilot (n  =  21) was performed to investigate the usability of our method. Operator time for the planning process was reduced to  <5 min. The retrospective planning study showed that 98 plans met all clinical constraints. Significant improvements were made in the volume receiving 72Gy (V72Gy) for the bladder and rectum and the mean dose of the bladder and the body. A reduced plan variance was observed. During the clinical pilot 20 automatically generated plans met all constraints and 17 plans were selected for treatment. The automated radiotherapy treatment planning and optimization workflow is capable of efficiently generating patient specifically optimized and improved clinical grade plans. It has now been adopted as the current standard workflow in our clinic to generate treatment plans for prostate cancer

Towards fast online intrafraction replanning for free-breathing stereotactic body radiation therapy with the MR-linac

The hybrid MRI-radiotherapy machines, like the MR-linac (Elekta AB, Stockholm, Sweden) installed at the UMC Utrecht (Utrecht, The Netherlands), will be able to provide real-time patient imaging during treatment. In order to take advantage of the system's capabilities and enable online adaptive treatments, a new generation of software should be developed, ranging from motion estimation to treatment plan adaptation. In this work we present a proof of principle adaptive pipeline designed for high precision stereotactic body radiation therapy (SBRT) suitable for sites affected by respiratory motion, like renal cell carcinoma (RCC). We utilized our research MRL treatment planning system (MRLTP) to simulate a single fraction 25 Gy free-breathing SBRT treatment for RCC by performing inter-beam replanning for two patients and one volunteer. The simulated pipeline included a combination of (pre-beam) 4D-MRI and (online) 2D cine-MR acquisitions. The 4DMRI was used to generate the mid-position reference volume, while the cine-MRI, via an in-house motion model, provided three-dimensional (3D) deformable vector fields (DVFs) describing the anatomical changes during treatment. During the treatment fraction, at an inter-beam interval, the mid-position volume of the patient was updated and the delivered dose was accurately reconstructed on the underlying motion calculated by the model. Fast online replanning, targeting the latest anatomy and incorporating the previously delivered dose was then simulated with MRLTP. The adaptive treatment was compared to a conventional mid-position SBRT plan with a 3 mm planning target volume margin reconstructed on the same motion trace. We demonstrate that our system produced tighter dose distributions and thus spared the healthy tissue, while delivering more dose to the target. The pipeline was able to account for baseline variations/drifts that occurred during treatment ensuring target coverage at the end of the treatment fraction

Intrafraction Motion Management of Renal Cell Carcinoma With Magnetic Resonance Imaging-Guided Stereotactic Body Radiation Therapy

Purpose: One of the major challenges in stereotactic body radiation therapy (SBRT) of renal cell carcinoma is internal motion during treatment. Previous literature has aimed to mitigate the effects of motion by expanding the treatment margins or respiratory tracking. Online magnetic resonance imaging (MRI)-guided radiation therapy has the potential to further improve the treatment of renal cell carcinoma by direct visualization of the tumor during treatment. The efficacy of 2 motion management techniques were assessed: tumor trailing and respiratory tracking. The simulation of a single-fraction, MRI-based SBRT was performed to quantify intrafraction motion and assess the efficacy of the different motion management strategies. Methods and materials: Fifteen patients were included in the study. At the beginning and end of the scanning protocol, 2 cine-MRI scans were acquired to assess cyclic respiratory motion. In addition, 3-dimensional spoiled gradient echo scans were acquired at 4 different time points to assess the slow drifts over 25 minutes. The systematic and random errors owing to intrafraction drift were calculated, as well as the random error induced by respiratory motion. The motion margins were calculated for tumor trailing and respiratory tracking and compared with the margin when no motion compensation would be performed to assess the relative efficacy of each technique. Results: The largest respiratory tumor motion was observed along the caudo-cranial direction with a median 95% maximum amplitude of approximately 12 mm. Σ DRIFT, σ DRIFT, and σ RESP were determined to be 1.0 mm 1.8 mm, and 3.8 mm, respectively. Without mechanical immobilization, intrafraction drift accounted for 75% of the total intrafraction motion margin for online midposition-based SBRT treatments. Conclusions: The contribution of intrafraction drift to the total internal motion margin is much larger than periodic respiratory motion. This makes tumor trailing a viable option to consider on the MRI linac because it allows for 3-dimensional MRI acquisitions during beam delivery, which simplifies the introduction of new techniques, such as dose accumulation and online intrafraction replanning

Dosimetric impact of soft-tissue based intrafraction motion from 3D cine-MR in prostate SBRT

To investigate the dosimetric impact of intrafraction translation and rotation motion of the prostate, as extracted from daily acquired post-treatment 3D cine-MR based on soft-tissue contrast, in extremely hypofractionated (SBRT) prostate patients. Accurate dose reconstruction is performed by using a prostate intrafraction motion trace which is obtained with a soft-tissue based rigid registration method on 3D cine-MR dynamics with a temporal resolution of 11 s. The recorded motion of each time-point was applied to the planning CT, resulting in the respective dynamic volume used for dose calculation. For each treatment fraction, the treatment delivery record was generated by proportionally splitting the plan into 11 s intervals based on the delivered monitor units. For each fraction the doses of all partial plan/dynamic volume combinations were calculated and were summed to lead to the motion-affected fraction dose. Finally, for each patient the five fraction doses were summed, yielding the total treatment dose. Both daily and total doses were compared to the original reference dose of the respective patient to assess the impact of the intrafraction motion. Depending on the underlying motion of the prostate, different types of motion-affected dose distributions were observed. The planning target volumes (PTVs) ensured CTV_30 (seminal vesicles) D99% coverage for all patients, CTV_35 (prostate corpus) coverage for 97% of the patients and GTV_50 (local boost) for 83% of the patients when compared against the strict planning target D99% value. The dosimetric impact due to prostate intrafraction motion in extremely hypofractionated treatments was determined. The presented study is an essential step towards establishing the actual delivered dose to the patient during radiotherapy fractions

Fiducial marker based intra-fraction motion assessment on cine-MR for MR-linac treatment of prostate cancer

We have developed a method to determine intrafraction motion of the prostate through automatic fiducial marker (FM) tracking on 3D cine-magnetic resonance (MR) images with high spatial and temporal resolution. Twenty-nine patients undergoing prostate stereotactic body radiotherapy (SBRT), with four implanted cylindrical gold FMs, had cine-MR imaging sessions after each of five weekly fractions. Each cine-MR examination consisted of 55 sequentially obtained 3D datasets ('dynamics'), acquired over a 11 s period, covering a total of 10 min. FM locations in the first dynamic were manually identified by a clinician, FM centers in subsequent dynamics were automatically determined. Center of mass (COM) translations and rotations were determined by calculating the rigid transformations between the FM template of the first and subsequent dynamics. The algorithm was applied to 7315 dynamics over 133 scans of 29 patients and the obtained results were validated by comparing the COM locations recorded by the clinician at the halfway-dynamic (after 5 min) and end dynamic (after 10 min). The mean COM translations at 10 min were X: 0.0 0.8 mm, Y: 1.0 1.9 mm and Z: 0.9 2.0 mm. The mean rotation results at 10 min were X: 0.1 3.9°, Y: 0.0 1.3° and Z: 0.1 1.2°. The tracking success rate was 97.7% with a mean 3D COM error of 1.1 mm. We have developed a robust, fast and accurate FM tracking algorithm for cine-MR data, which allows for continuous monitoring of prostate motion during MR-guided radiotherapy (MRgRT). These results will be used to validate automatic prostate tracking based on soft-tissue contrast