Effect of intra-fraction motion on the accumulated dose for free-breathing MR-guided stereotactic body radiation therapy of renal-cell carcinoma

A systematic review of 4D magnetic resonance imaging techniques for abdominal radiotherapy treatment planning

Background and purpose

Four-dimensional magnetic resonance imaging (4DMRI) has gained interest as an alternative to the current standard for motion management four-dimensional tomography (4DCT) in abdominal radiotherapy treatment planning (RTP). This review aims to assess the 4DMRI literature in abdomen, focusing on technical considerations and the validity of using 4DMRI for patients within radiotherapy protocols.

Materials and methods

The review followed the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. A comprehensive search was performed across the Medline, Embase, Scopus, and Web of Science databases, covering all years up to December 31, 2023. The studies were grouped into two categories: 4DMRI reconstructed from 3DMRI acquisition; and 4DMRI reconstructed from multi-slice 2DMRI acquisition.

Results

A total of 39 studies met the inclusion criteria and were analysed to provide key findings. Key findings were 4DMRI had the potential to improve abdominal RTP for patients by providing accurate tumour definition and motion assessment compared to 4DCT. 4DMRI reconstructed from 3DMRI acquisition showed promise as a feasible approach for motion management in abdominal RTP regarding spatial resolution. Currently,the slice thickness achieved on 4DMRI reconstructed from multi-slice 2DMRI acquisitions was unsuitable for clinical purposes. Lastly, the current barriers for clinical implementation of 4DMRI were the limited availability of validated commercial solutions and the lack of larger cohort comparative studies to 4DCT for target delineation and plan optimisation.

Conclusion

4DMRI showed potential improvements in abdominal RTP, but standards and guidelines for the use of 4DMRI in radiotherapy were required to demonstrate clinical benefits.

Optimising the MR-Linac as a standard treatment modality

The magnetic resonance linear accelerator (MR-Linac) offers a new treatment paradigm, providing improved visualisation of targets and organs at risk while allowing for daily adaptation of treatment plans in real time. Online MR-guided adaptive treatment has reduced treatment uncertainties; however, the additional treatment time and resource requirements may be a concern. We present our experience of integrating an MR-Linac into a busy department and provide recommendations for improved clinical and resource efficiency. Furthermore, we discuss potential future technological innovations that can further optimise clinical productivity in a busy department.

Local Therapies for Hepatocellular Carcinoma and Role of MRI-Guided Adaptive Radiation Therapy

Hepatocellular carcinoma (HCC) is the most common liver tumor, with a continually rising incidence. The curative treatment for HCC is surgical resection or liver transplantation; however, only a small portion of patients are eligible due to local tumor burden or underlying liver dysfunction. Most HCC patients receive nonsurgical liver-directed therapies (LDTs), including thermal ablation, transarterial chemoembolization (TACE), transarterial radioembolization (TARE), and external beam radiation therapy (EBRT). Stereotactic ablative body radiation (SABR) is a specific type of EBRT that can precisely deliver a high dose of radiation to ablate tumor cells using a small number of treatments (or fractions, typically 5 or less). With onboard MRI imaging, MRI-guided SABR can improve therapeutic dose while minimizing normal tissue exposure. In the current review, we discuss different LDTs and compare them with EBRT, specifically SABR. The emerging MRI-guided adaptive radiation therapy has been reviewed, highlighting its advantages and potential role in HCC management.

Development of an abdominal dose accumulation tool and assessments of accumulated dose in gastrointestinal organs

Objective.Online adaptive radiotherapy has demonstrated improved dose conformality in response to inter-fraction geometric variations in the abdomen. The dosimetric impact of intra-fractional variations in anatomic configuration resulting from breathing, gastric contraction and slow configuration motion, however, have been largely ignored, leading to differences between delivered and planned. To investigate the impact of intra-fractional abdominal motions on delivered dose, anatomical deformations due to these three motion modes were extracted from dynamic MRI data using a previously developed hierarchical motion modeling methodology.Approach. Motion magnitudes were extracted from deformation fields between a reference state and all other motion states of the patient. Delivered dose estimates to various gastrointestinal organs (stomach, duodenum, small bowel and colon) were calculated on each motion state of the patient and accumulated to estimate the delivered dose to each organ for the entire treatment fraction.Main results. Across a sample of 10 patients, maximal motions of 33.6, 33.4, 47.6 and 49.2 mm were observed over 20 min for the stomach, duodenum, small bowel and colon respectively. Dose accumulation results showed that motions could lead to average increases of 2.0, 2.1, 1.1, 0.7 Gy to the maximum dose to 0.5cc (D0.5cc) and 3.0, 2.5, 1.3, 0.9 Gy to the maximum dose to 0.1cc (D0.1cc) for these organs at risk. From the 40 dose accumulations performed (10 for each organ at risk), 27 showed increases of modeled delivered dose compared to planned doses, 4 of which exceeded planned dose constraints.Significance. The use of intra-fraction motion measurements to accumulate delivered doses is feasible, and supports retrospective estimation of dose delivery to improve estimates of delivered doses, and further guide strategies for both plan adaptation as well as advances in intra-fraction motion management.

A Simulation Study of Tolerance of Breathing Amplitude Variations in Radiotherapy of Lung Cancer Using 4DCT and Time-Resolved 4DMRI

As patient breathing irregularities can introduce a large uncertainty in targeting the internal tumor volume (ITV) of lung cancer patients, and thereby affect treatment quality, this study evaluates dose tolerance of tumor motion amplitude variations in ITV-based volumetric modulated arc therapy (VMAT). A motion-incorporated planning technique was employed to simulate treatment delivery of 10 lung cancer patients' clinical VMAT plans using original and three scaling-up (by 0.5, 1.0, and 2.0 cm) motion waveforms from single-breath four-dimensional computed tomography (4DCT) and multi-breath time-resolved 4D magnetic resonance imaging (TR-4DMRI). The planning tumor volume (PTV = ITV + 5 mm margin) dose coverage (PTV D95%) was evaluated. The repeated waveforms were used to move the isocenter in sync with the clinical leaf motion and gantry rotation. The continuous VMAT arcs were broken down into many static beam fields at the control points (2°-interval) and the composite plan represented the motion-incorporated VMAT plan. Eight motion-incorporated plans per patient were simulated and the plan with the native 4DCT waveform was used as a control. The first (D95% ≤ 95%) and second (D95% ≤ 90%) plan breaching points due to motion amplitude increase were identified and analyzed. The PTV D95% in the motion-incorporated plans was 99.4 ± 1.0% using 4DCT, closely agreeing with the corresponding ITV-based VMAT plan (PTV D95% = 100%). Tumor motion irregularities were observed in TR-4DMRI and triggered D95% ≤ 95% in one case. For small tumors, 4 mm extra motion triggered D95% ≤ 95%, and 6-8 mm triggered D95% ≤ 90%. For large tumors, 14 mm and 21 mm extra motions triggered the first and second breaching points, respectively. This study has demonstrated that PTV D95% breaching points may occur for small tumors during treatment delivery. Clinically, it is important to monitor and avoid systematic motion increase, including baseline drift, and large random motion spikes through threshold-based beam gating.

Stereotactic Body Radiotherapy for Localized Kidney Cancer

PURPOSE OF REVIEW

Stereotactic body radiation therapy (SBRT) is increasingly utilized in the management of localized kidney cancers, particularly for patients who are not surgical candidates. Herein, we provide a narrative review of SBRT in the management of localized kidney cancers.

RECENT FINDINGS

Recent prospective studies and multi-institutional retrospective studies highlight the safety and efficacy of SBRT in the management of renal tumors, a disease previously thought to be radioresistant. Studies have shown that local control is greater than 90% with rare grade 3 or 4 toxicity and no grade 5 toxicity. SBRT can be utilized successfully in the treatment of large kidney tumors (> 5 cm). New techniques such as MRI-guided radiation therapy may further improve the therapeutic ratio. However, randomized clinical trials are necessary to confirm the optimal dosing schedule and compare outcomes with nephrectomy, which remains the standard of care in suitable patients. Advances in SBRT have made this modality a safe and effective treatment option in the management of localized kidney cancers.

Improving the clinical workflow of a MR-Linac by dosimetric evaluation of synthetic CT

Adaptive radiotherapy performed on the daily magnetic resonance imaging (MRI) is an option to improve the treatment quality. In the adapt-to-shape workflow of 1.5-T MR-Linac, the contours of structures are adjusted on the basis of patient daily MRI, and the adapted plan is recalculated on the MRI-based synthetic computed tomography (syCT) generated by bulk density assignment. Because dosimetric accuracy of this strategy is a priority and requires evaluation, this study aims to explore the usefulness of adding an assessment of dosimetric errors associated with recalculation on syCT to the clinical workflow. Sixty-one patients, with various tumor sites, treated using a 1.5-T MR-Linac were included in this study. In Monaco V5.4, the target and organs at risk (OARs) were contoured, and a reference CT plan that contains information about the outlined contours, their average electron density (ED), and the priority of ED assignment was generated. To evaluate the dosimetric error of syCT caused by the inherent approximation within bulk density assignment, the reference CT plan was recalculated on the syCT obtained from the reference CT by forcing all contoured structures to their mean ED defined on the reference plan. The dose–volume histogram (DVH) and dose distribution of the CT and syCT plan were compared. The causes of dosimetric discrepancies were investigated, and the reference plan was reworked to minimize errors if needed. For 54 patients, gamma analysis of the dose distribution on syCT and CT show a median pass rate of 99.7% and 98.5% with the criteria of 3%/3 mm and 2%/2 mm, respectively. DVH difference of targets and OARs remained less than 1.5% or 1 Gy. For the remaining patients, factors (i.e., inappropriate ED assignments) influenced the dosimetric agreement of the syCT vs. CT reference DVH by up to 21%. The causes of the errors were promptly identified, and the DVH dosimetry was realigned except for two lung treatments for which a significant discrepancy remained. The recalculation on the syCT obtained from the planning CT is a powerful tool to assess and decrease the minimal error committed during the adaptive plan on the MRI-based syCT.

Time-course assessment of 3D-image distortion on the 1.5 T Marlin/Elekta Unity MR-LINAC

PURPOSE

The efficacy of MR-guided radiotherapy on a MR-LINAC (MR-L) is dependent on the geometric accuracy of its MR images over clinically relevant Fields-of-View (FOVs). Our objectives were to: evaluate gradient non-linearity (GNL) on the Elekta Unity MR-L across time via 76 weekly measurements of 3D-distortion over concentrically larger diameter spherical volumes (DSVs); quantify distortion measurement error; and assess the temporal stability of spatial distortion using statistical process control (SPC).

METHODS

MR-image distortion was assessed using a large-FOV 3D-phantom containing 1932 markers embedded in seven parallel plates, spaced 25 mm × 25 mm in- and 55 mm through-plane. Automatically analyzed T1 images yielded distortions in 200, 300, 400 and 500 mm concentric DSVs. Distortion measurement error was evaluated using median absolute difference analysis of imaging repeatability tests.

RESULTS

Over the measurement period absolute time-averaged distortion varied between: dr = 0.30 - 0.49 mm, 0.53 - 0.80 mm, 1.0 - 1.4 mm and 2.28 - 2.37 mm, for DSVs 200, 300, 400 and 500 mm at the 98^th percentile level. Repeatability tests showed that imaging/repositioning introduces negligible error: mean ≤ 0.02 mm (max ≤ 0.3 mm). SPC analysis showed image distortion was stable across all DSVs; however, noticeable changes in GNL were observed following servicing at the one-year mark.

CONCLUSIONS

Image distortion on the MR-L is in the sub-millimeter range for DSVs ≤ 300 mm and stable across time, with SPC analysis indicating all measurements remain within control for each DSV.

Assessing organ at risk position variation and its impact on delivered dose in kidney SABR

Background

Delivered organs at risk (OARs) dose may vary from planned dose due to interfraction and intrafraction motion during kidney SABR treatment. Cases of bowel stricture requiring surgery post SABR treatment were reported in our institution. This study aims to provide strategies to reduce dose deposited to OARs during SABR treatment and mitigate risk of gastrointestinal toxicity.

Methods

Small bowel (SB), large bowel (LB) and stomach (STO) were delineated on the last cone beam CT (CBCT) acquired before any dose had been delivered (PRE CBCT) and on the first CBCT acquired after any dose had been delivered (MID CBCT). OAR interfraction and intrafraction motion were estimated from the shortest distance between OAR and the internal target volume (ITV). Adaptive radiation therapy (ART) was used if dose limits were exceeded by projecting the planned dose on the anatomy of the day.

Results

In 36 patients, OARs were segmented on 76 PRE CBCTs and 30 MID CBCTs. Interfraction motion was larger than intrafraction motion in STO (p-value = 0.04) but was similar in SB (p-value = 0.8) and LB (p-value = 0.2). LB was inside the planned 100% isodose in all PRE CBCTs and MID CBCTs in the three patients that suffered from bowel stricture. SB D0.03cc was exceeded in 8 fractions (4 patients). LB D1.5cc was exceeded in 4 fractions (2 patients). Doses to OARs were lowered and limits were all met with ART on the anatomy of the day.

Conclusions

Interfraction motion was responsible for OARs overdosage. Dose limits were respected by using ART with the anatomy of the day.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-022-02041-2.

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.

Technical Note: End-to-end verification of an MR-Linac using a dynamic motion phantom

PURPOSE

MR-Linac integrates an MRI scanner and a linear accelerator to provide adaptive radiation treatment. Superior tissue contrast and real-time imaging can give the clinicians confidence to reduce the margins of the planning target volume (PTV). The purpose of this study was to verify the dosimetric accuracy of an MR-Linac system in treating a moving target and assess the error with different motion patterns and adaptation methods.

METHODS

We performed an end-to-end test for Elekta Unity (Elekta) using the 4D Dynamic Thorax Phantom (CIRS MRgRT 008Z), comparing the measured and planned dose. The moving phantom had four measurement locations in the tumor, liver, kidney, and spinal cord regions with a PTW30013 ion chamber. For seven different motion patterns, we first acquired simulation CT using a slow-scanning protocol, based on which we generated reference plans. The treatment technique was the standard intensity-modulated radiation therapy (IMRT). We tested both adaptation workflows: the Adapt-to-Position (ATP) and the Adapt-to-Shape (ATS). The three-dimensional (3D) distribution was measured using a diode array phantom (Sun Nuclear Inc.) to check the dose distribution accuracy as part of the routine QA process. We also performed end-to-end tests on a conventional Linac. Finally, we used SPSS Statistics 22.0 (Inc., Chicago, IL, USA) for data analysis.

RESULTS

All pretreatment reference plans and delivered plans had excellent QA results with a better than 95% passing rate of relative gamma analysis (2%/2 mm criteria). The adaptive planning for MR-Linac produced quality plans. The measured dose in the target agreed with the calculated dose.

CONCLUSIONS

The adaptive treatment on the MR-Linac system investigated met the expected performance with tumor motions. The outline of the target could be visualized and accurately contoured on the 3D MR for online planning. Under different motion patterns, the difference between the measured and calculated dose was acceptable clinically.

Adaptive Magnetic Resonance-Guided Stereotactic Body Radiotherapy: The Next Step in the Treatment of Renal Cell Carcinoma

Adaptive MR-guided radiotherapy (MRgRT) is a new treatment paradigm and its role as a non-invasive treatment option for renal cell carcinoma is evolving. The early clinical experience to date shows that real-time plan adaptation based on the daily MRI anatomy can lead to improved target coverage and normal tissue sparing. Continued technological innovations will further mitigate the challenges of organ motion and enable more advanced treatment adaptation, and potentially lead to enhanced oncologic outcomes and preservation of renal function. Future applications look promising to make a positive clinical impact and further the personalization of radiotherapy in the management of renal cell carcinoma.

Technical Note: Cumulative dose modeling for organ motion management in MRI-guided radiation therapy

PURPOSE

To develop a method for continuous online dose accumulation during irradiation in MRI-guided radiation therapy (MRgRT) and to demonstrate its application in evaluating the impact of internal organ motion on cumulative dose.

METHODS

An intensity-modulated radiation therapy (IMRT) treatment plan is partitioned into its unique apertures. Dose for each planned aperture is calculated using Monte Carlo dose simulation on each phase of a four-dimensional computed tomography (4D-CT) dataset. Deformable image registration is then performed both (a) between each frame of a cine-MRI acquisition obtained during treatment and a reference frame, and (b) between each volume of the 4D-CT phases and a reference phase. These registrations are used to associate each cine image with a 4D-CT phase. Additionally, for each 4D-CT phase, the deformation vector field (DVF) is used to warp the pre-calculated dose volumes per aperture onto the reference CT dataset. To estimate the dose volume delivered during each frame of the cine-MRI acquisition, we retrieve the pre-calculated warped dose volume for the delivered aperture on the associated 4D-CT phase and adjust it by a rigid translation to account for baseline drift and instances where motion on the cine image exceeds the amplitude observed between 4D-CT phases.

RESULTS

The proposed dose accumulation method is retrospectively applied to a liver cancer case previously treated on an MRgRT platform. Cumulative dose estimated for free-breathing and respiration-gated delivery is compared against dose calculated on static anatomy. In this sample case, the target minimum dose and D 98 varied by as much as 5% and 7%, respectively.

CONCLUSION

We demonstrate a technique suitable for continuous online dose accumulation during MRgRT. In contrast to other approaches, dose is pre-calculated per aperture and phase and then retrieved based on a mapping scheme between cine MRI and 4D-CT datasets, aiming at reducing the computational burden for potential real-time applications.

In-Air Electron Streaming Effect for Esophageal Cancer Radiotherapy With a 1.5 T Perpendicular Magnetic Field: A Treatment Planning Study

Purpose

To investigate the in-air out-of-field electron streaming effect (ESE) for esophageal cancer radiotherapy in the presence of 1.5 T perpendicular magnetic field.

Methods

Ten esophageal cancer patients treated with conventional Linac were retrospectively enrolled into a cohort of this study, with the prescription of 4,400 cGy/20 fx. All cases received IMRT replanning using Elekta Unity MR-Linac specified Monaco system, denoted as primary plan. To visualize the in-air dose outside the body in Monaco system, an auxiliary structure was created by extending the external structure. For each case, another comparable plan with no magnetic field was created using the same planning parameters. The plan was also recalculated by placing a bolus upon the neck and chin area to investigate its shielding effect for ESE. Dosimetric evaluations of the out-of-field neck and chin skin area and statistical analysis for these plans were then performed.

Results

Out-of-field ESE was also observed in esophageal cancer treatment planning under 1.5 T magnetic field, while totally absent for plans with no magnetic field. On average, the maximum dose to the neck and chin skin area of the primary plan (657.92 ± 69.07 cGy) was higher than that of plan with no magnetic field (281.78 ± 36.59 cGy, p = 0.005) and plan with bolus (398.43 ± 69.19 cGy, p = 0.007). DVH metrics D_1cc (the minimum dose to 1 cc volume) of the neck and chin skin for primary plan was 382.06 ± 44.14 cGy, which can be reduced to 212.42 ± 23.65 cGy by using the 1 cm bolus (with p = 0.005), even lower than the plan without magnetic field (214.45 ± 23.82, p = 0.005). No statistically significant difference of the neck and chin skin dose between the plan with bolus and plan with no magnetic field was observed (all with p > 0.05).

Conclusion

For MRI guided esophageal cancer radiotherapy, a relatively high out-of-field neck and chin skin doses will be introduced by ESE in the presence of magnetic field. It is therefore recommended to take this into account during the planning phase. Adding bolus could effectively reduce the ESE dose contributions, achieve the shielding effect almost equivalent to the scenario with no magnetic field. Further explorations of measurement verifications for the ESE dose distributions are required.

Dosimetric feasibility of hypofractionation for SBRT treatment of lymph node oligometastases on the 1.5T MR-linac

PURPOSE

At our department, MR-guided stereotactic body radiation therapy (SBRT) using the 1.5T MR-linac system (Unity, Elekta AB, Stockholm, Sweden) has been initiated for patients with lymph node oligometastases. Superior soft tissue contrast and the possibility for online plan adaptation on the Unity may allow for hypofractionated treatment. The purpose of this study was to investigate the dosimetric feasibility and compare the plan quality of different hypofractionated schemes.

METHODS AND MATERIALS

Data was used from 12 patients with single lymph node oligometastases (10 pelvic, 2 para-aortic), which were all treated on the Unity with a prescribed dose of 5x7 Gy to 95% of the PTV. Hypofractionation was investigated for 3x10 Gy and 1x20 Gy schemes (all 60 Gy BED α/β = 10). The pre-treatment plans were evaluated based on dose criteria and plan quality. If all criteria were met, the number of online adapted plans which also met all dose criteria was investigated. For pre-treatment plans meeting the criteria for all three fractionation schemes, the plan quality after online adaptation was compared using the four parameters described in the NRG-BR001 phase 1 trial.

RESULTS

Pre-treatment plans met all clinical criteria for the three different fractionation schemes in 10, 9 and 6 cases. 50/50, 45/45 17/30 of the corresponding online adapted plans met all criteria, respectively. Violations were primarily caused by surrounding organs at risk overlapping or adjacent to the PTV. The 1x20 Gy treatment plans were, in general, of lesser quality than the 5x7 Gy and 3x10 Gy plans.

CONCLUSION

Hypofractionated radiotherapy for lymph node oligometastases on the 1.5T MR-linac is feasible based on dose criteria and plan quality metrics. The location of the target relative to critical structures should be considered in choosing the most suitable fractionation scheme. Especially for single fraction treatment, meeting all dose criteria in the pre-treatment situation does not guarantee that this also applies during online treatment.

A modern review of the uncertainties in volumetric imaging of respiratory-induced target motion in lung radiotherapy

Radiotherapy has become a critical component for the treatment of all stages and types of lung cancer, often times being the primary gateway to a cure. However, given that radiation can cause harmful side effects depending on how much surrounding healthy tissue is exposed, treatment of the lung can be particularly challenging due to the presence of moving targets. Careful implementation of every step in the radiotherapy process is absolutely integral for attaining optimal clinical outcomes. With the advent and now widespread use of stereotactic body radiation therapy (SBRT), where extremely large doses are delivered, accurate, and precise dose targeting is especially vital to achieve an optimal risk to benefit ratio. This has largely become possible due to the rapid development of image-guided technology. Although imaging is critical to the success of radiotherapy, it can often be plagued with uncertainties due to respiratory-induced target motion. There has and continues to be an immense research effort aimed at acknowledging and addressing these uncertainties to further our abilities to more precisely target radiation treatment. Thus, the goal of this article is to provide a detailed review of the prevailing uncertainties that remain to be investigated across the different imaging modalities, as well as to highlight the more modern solutions to imaging motion and their role in addressing the current challenges.

Reliability of ITV approach to varying treatment fraction time: a retrospective analysis based on 2D cine MR images

BACKGROUND

Internal Target Volume (ITV) is one of the most common strategies to passively manage tumour motion in Radiotherapy (RT). The reliability of this approach is based on the assumption that the tumour motion estimated during pre-treatment 4D Computed Tomography (CT) acquisition is representative of the motion during the whole RT treatment. With the introduction of Magnetic Resonance-guided RT (MRgRT), it has become possible to monitor tumour motion during the treatment and verify this assumption. Aim of this study was to investigate the reliability of the ITV approach with respect to the treatment fraction time (TFT) in abdominal and thoracic lesions.

METHODS

A total of 12 thoracic and 15 abdominal lesions was analysed. Before treatment, a 10-phase 4DCT was acquired and ITV margins were estimated considering the envelope of the lesion contoured on the different 4DCT phases. All patients underwent MRgRT treatment in free-breathing, monitoring the tumour position on a sagittal plane with 4 frames per second (sec). ITV margins were projected on the tumour trajectory and the percentage of treatment time in which the tumour was inside the ITV (%TT) was measured to varying of TFT. The ITV approach was considered moderately reliable when %TT ≥ 90% and strongly reliable when %TT ≥ 95%. Additional ITV margins required to achieve %TT ≥ 95% were also calculated.

RESULTS

In the analysed cohort of patients, ITV strategy can be considered strongly reliable only for lung lesions with TFT ≤ 7 min (min). The ITV strategy can be considered only moderately reliable for abdominal lesions, and additional margins are required to obtain %TT ≥ 95%. Considering a TFT ≤ 4 min, additional margins of 2 mm in cranio-caudal (CC) and 1 mm in antero-posterior (AP) are suggested for pancreatic lesions, 3 mm in CC and 2 mm in AP for renal and liver ones.

CONCLUSIONS

On the basis of the analysed cases, the ITV approach appears to be reliable in the thorax, while it results more challenging in the abdomen, due to the higher uncertainty in ITV definition and to the observed larger intra and inter-fraction motion variability. The addition of extra margins based on the TFT may represent a valid tool to compensate such limitations.

Interest of positron-emission tomography and magnetic resonance imaging for radiotherapy planning and control

Computed tomography (CT) in the treatment position is currently indispensable for planning radiation therapy. Other imaging modalities, such as magnetic resonance imaging (MRI) and positron emission-tomography (PET), can be used to improve the definition of the tumour and/or healthy tissue but also to provide functional data of the target volume. Accurate image registration is essential for treatment planning, so MRI and PET scans should be registered at the planning CT scan. Hybrid PET/MRI scans with a hard plane can be used but pose the problem of the absence of CT scans. Finally, techniques for moving the patient on a rigid air-cushioned table allow PET/CT/MRI scans to be performed in the treatment position while limiting the patient's movements exist. At the same time, the advent of MRI-linear accelerator systems allows to redefine image-guided radiotherapy and to propose treatments with daily recalculation of the dose. The place of PET during treatment remains more confidential and currently only in research and prototype status. The same development of imaging during radiotherapy is underway in proton therapy.

Time-resolved volumetric MRI in MRI-guided radiotherapy: an in silico comparative analysis

MRI-treatment units enable 2D cine-MRI centred in the tumour for motion detection in radiotherapy, but they lack 3D information due to spatio-temporal limits. To derive time-resolved 3D information, different approaches have been proposed in the literature, but a rigorous comparison among these strategies has not yet been performed. The goal of this study is to quantitatively investigate five published strategies that derive time-resolved volumetric MRI in MRI-guided radiotherapy: Propagation, out-of-plane motion compensation, Fayad model, ROI-based model and Stemkens model. Comparisons were performed using an MRI digital phantom generated with six different patient-derived motion signals and tumour-shapes. An average 4D cycle was generated as well as 2D cine-MRI data with corresponding 3D in-room ground truth. Quantitative analysis was performed by comparing the estimated 3D volume to the ground truth available for each 2D cine-MRI sample. A grouped patient statistical analysis was performed to evaluate the performance of the selected methods, in case of tumour tracking or motion estimation of the whole anatomy. Analyses were also performed based on patient characteristics. Quantitative ranking of the investigated methods highlighted that Propagation and ROI-based model strategies achieved an overall median tumour centre of mass 3D distance from the ground truth of 1.1 mm and 1.3 mm, respectively, and a diaphragm distance below 1.6 mm. Higher errors and variabilities were instead obtained for other methods, which lack the ability to compensate for in-room variations and to account for regional changes. These results were especially evident when further analysing patient characteristics, where errors above 2 mm/5 mm in tumour/diaphragm were found for more irregular breathing patterns in case of out-of-plane motion compensation, Fayad and Stemkens models. These findings suggest the potential of the proposed in silico framework to develop and compare strategies to estimate time-resolved 3DMRI in MRI-guided radiotherapy.

ReconSocket: a low-latency raw data streaming interface for real-time MRI-guided radiotherapy

With the recent advent of hybrid MRI-guided radiotherapy systems, continuous intra-fraction MR imaging for motion monitoring has become feasible. The ability to perform real-time custom image reconstructions is however often lacking. In this work we present a low-latency streaming solution, ReconSocket, which provides a real-time stream of k-space data from the magnetic resonance imaging (MRI) to custom reconstruction servers. We determined the performance of the data streaming by measuring the streaming latency (i.e. non-zero time delay due to data transfer and processing) and jitter (i.e. deviations from periodicity) using an ultra-fast 1D MRI acquisition of a moving phantom. Simultaneously, its position was recorded with near-zero time delay. The feasibility of low-latency custom reconstructions was tested by measuring the imaging latency (i.e. time delay between physical change and appearance of that change on the image) for several non-Cartesian 2D and 3D acquisitions using an in-house implemented reconstruction server. The measured streaming latency of the ReconSocket interface was [Formula: see text] ms. 98% of the incoming data packets arrived within a jitter range of 367 [Formula: see text]s. This shows that the ReconSocket interface can provide reliable real-time access to MRI data, acquired during the course of a MRI-guided radiotherapy fraction. The total imaging latency was measured to be 221 ms (2D) and 3889 ms (3D) for exemplary acquisitions, using the custom image reconstruction server. These imaging latencies are approximately equal to half of the temporal footprint (T _acq /2) of the respective 2D and 3D golden-angle radial sequences. For radial sequences, it was previously showed that T _acq /2 is the expected contribution of only the data acquisition to the total imaging latency. Indeed, the contribution of the non-Cartesian reconstruction to the total imaging latency was minor (<10%): 21 ms for 2D, 300 ms for 3D, indicating that the acquisition, i.e. the physical encoding of the image itself is the major contributor to the total imaging latency.

Adaptive radiotherapy: The Elekta Unity MR-linac concept

BACKGROUND AND PURPOSE

The promise of the MR-linac is that one can visualize all anatomical changes during the course of radiotherapy and hence adapt the treatment plan in order to always have the optimal treatment. Yet, there is a trade-off to be made between the time spent for adapting the treatment plan against the dosimetric gain. In this work, the various daily plan adaptation methods will be presented and applied on a variety of tumour sites. The aim is to provide an insight in the behavior of the state-of-the-art 1.5 T MRI guided on-line adaptive radiotherapy methods.

MATERIALS AND METHODS

To explore the different available plan adaptation workflows and methods, we have simulated online plan adaptation for five cases with varying levels of inter-fraction motion, regions of interest and target sizes: prostate, rectum, esophagus and lymph node oligometastases (single and multiple target). The plans were evaluated based on the clinical dose constraints and the optimization time was measured.

RESULTS

The time needed for plan adaptation ranged between 17 and 485 s. More advanced plan adaptation methods generally resulted in more plans that met the clinical dose criteria. Violations were often caused by insufficient PTV coverage or, for the multiple lymph node case, a too high dose to OAR in the vicinity of the PTV. With full online replanning it was possible to create plans that met all clinical dose constraints for all cases.

CONCLUSION

Daily full online replanning is the most robust adaptive planning method for Unity. It is feasible for specific sites in clinically acceptable times. Faster methods are available, but before applying these, the specific use cases should be explored dosimetrically.

MR-guidance in clinical reality: current treatment challenges and future perspectives

Magnetic Resonance-guided radiotherapy (MRgRT) marks the beginning of a new era. MR is a versatile and suitable imaging modality for radiotherapy, as it enables direct visualization of the tumor and the surrounding organs at risk. Moreover, MRgRT provides real-time imaging to characterize and eventually track anatomical motion. Nevertheless, the successful translation of new technologies into clinical practice remains challenging. To date, the initial availability of next-generation hybrid MR-linac (MRL) systems is still limited and therefore, the focus of the present preview was on the initial applicability in current clinical practice and on future perspectives of this new technology for different treatment sites.MRgRT can be considered a groundbreaking new technology that is capable of creating new perspectives towards an individualized, patient-oriented planning and treatment approach, especially due to the ability to use daily online adaptation strategies. Furthermore, MRL systems overcome the limitations of conventional image-guided radiotherapy, especially in soft tissue, where target and organs at risk need accurate definition. Nevertheless, some concerns remain regarding the additional time needed to re-optimize dose distributions online, the reliability of the gating and tracking procedures and the interpretation of functional MR imaging markers and their potential changes during the course of treatment. Due to its continuous technological improvement and rapid clinical large-scale application in several anatomical settings, further studies may confirm the potential disruptive role of MRgRT in the evolving oncological environment.

Assessment of 3D motion modeling performance for dose accumulation mapping on the MR-linac by simultaneous multislice MRI

Hybrid MR-linac systems enable intrafraction motion monitoring during radiation therapy. Since time-resolved 3D MRI is still challenging, various motion models have been developed that rely on time-resolved 2D imaging. Continuous validation of these models is important for accurate dose accumulation mapping. In this study we used 2D simultaneous multislice (SMS) imaging to improve the PCA-based motion modeling method developed previously (Stemkens et al 2016 Phys. Med. Biol. 61 5335-55). From the additional simultaneously acquired slices, several independent motion models could be generated, which allowed for an assessment of the sensitivity of the motion model to the location of the time-resolved 2D slices. Additionally, the best model could be chosen at every time-point, increasing the method's robustness. Imaging experiments were performed in six healthy volunteers using three simultaneous slices, which generated three independent models per volunteer. For each model the motion traces of the liver tip and both kidneys were estimated. We found that the location of the 2D slices influenced the model's error in five volunteers significantly with a p -value  <0.05, and that selecting the best model at every time-point can improve the method. This allows for more accurate and robust motion characterization in MR-guided radiotherapy.

A dual-purpose MRI acquisition to combine 4D-MRI and dynamic contrast-enhanced imaging for abdominal radiotherapy planning

For successful abdominal radiotherapy it is crucial to have a clear tumor definition and an accurate characterization of the motion. While dynamic contrast-enhanced (DCE) MRI aids tumor visualization, it is often hampered by motion artifacts. 4D-MRI characterizes this motion, but often lacks the contrast to clearly visualize the tumor. This dual requirement is challenging due to time constraints. Here, we propose combining both into a single acquisition by reconstructing the data in various ways in order to achieve both high spatio-temporal resolution DCE-MRI and accurate 4D-MRI motion estimates. A 5 min T₁-weigthed DCE acquisition was collected in five renal-cell carcinoma patients and simulated in a digital phantom. Data were acquired continuously using a 3D golden angle radial stack-of-stars acquisition. This enabled three types of reconstruction; (1) a high spatio-temporal resolution DCE time series, (2) a 5D reconstruction and (3) a contrast-enhanced 4D-MRI for motion characterization. Motion extracted from the 4D- and 5D-MRI was compared with a separately acquired 4D-MRI and additional 2D cine MR imaging. Simulations on XCAT showed that 5D reconstructions severely underestimated motion ([Formula: see text]), whereas contrast-enhanced 4D-MRI only underestimated motion by [Formula: see text]. This was confirmed in the in vivo data where motion of the contrast-enhanced 4D-MRI was approximately [Formula: see text] smaller than the motion in the 2D cine MRI (5.8 mm versus 6.5 mm), but equal to a separately acquired 4D-MRI (5.8 mm versus 5.9 mm). 5D reconstructions underestimated the motion by more than [Formula: see text], but minimized respiratory-induced blurring in the contrast enhanced images. DCE time-series demonstrated clear tumor visualization and contrast enhancement throughout the entire field-of-view. DCE- and 4D-MRI can be integrated into a single acquisition that enables different reconstructions with complementary information for abdominal radiotherapy planning and, in an MRI-guided treatment, precise motion information, input for motion models, and rapid feedback on the contrast enhancement.

Multiresolution radial MRI to reduce IDLE time in pre-beam imaging on an MR-Linac (MR-RIDDLE)

Online adaptive MR-guided radiation therapy improves treatment quality at the expense of considerable longer treatment time. The treatment lengthening partially originates from the preparatory (pre-beam) MR imaging required to encode all the information needed for contour propagation, contour adaptation and replanning. MRI requires several minutes of scan time before the encoded information is converted to usable images, which results in long idle times before the first clinical tasks are performed. In this study we propose a novel imaging sequence, called MR-RIDDLE, that reduces the idle time and therefore speeds-up the workflow in online MR-guided radiation therapy. MR-RIDDLE enables multiresolution image reconstruction to commence during data acquisition where low resolution images are available within one minute, after which the data collection continuous for subsequent high-resolution image updates. We demonstrate that the low resolution images can be used to accurately propagate contours from the pre-treatment scan. For abdominothoracic tumours MR-RIDDLE inherently captures a motion-blurred representation of the mid-position, which we were able to deblur using a combination of an internal motion surrogate and auto-adaptive soft-gating filters. Our results demonstrate that MR-RIDDLE provides a robust, flexible and time-efficient strategy for pre-beam imaging, even for cases with large respiratory movements or baseline shifts within the acquisition. We anticipate that this novel concept of parallelising the MR imaging and the clinical tasks has the potential to considerably speed-up and streamline the online MR-guided radiation therapy workflow.

A ROI-based global motion model established on 4DCT and 2D cine-MRI data for MRI-guidance in radiation therapy

In-room magnetic resonance imaging (MRI) allows the acquisition of fast 2D cine-MRI centered in the tumor for advanced motion management in radiotherapy. To achieve 3D information during treatment, patient-specific motion models can be considered the most viable solution. However, conventional global motion models are built using a single motion surrogate, independently from the anatomical location. In this work, we present a novel motion model based on regions of interest (ROIs) established on 4D computed tomography (4DCT) and 2D cine-MRI, aiming at accurately compensating for changes during treatment. In the planning phase, a motion model is built on a 4DCT dataset, through 3D deformable image registration (DIR). ROIs are then defined and correlated with motion fields derived by 2D DIR between CT slices centered in the tumor. In the treatment phase, the model is applied to in-room cine-MRI data to compensate for organ motion in a multi-modal framework, aiming at estimating a time-resolved 3DCT. The method is validated on a digital phantom and tested on two lung patients. Analysis is performed by considering different anatomical planes (coronal, sagittal and a combination of the two) and evaluating the performance of the method on tumor and diaphragm. For the phantom study, the ROI-based model results in a uniform median error on both diaphragm and tumor below 1.5 mm. For what concerns patients, median errors on both diaphragm and tumor are around 2 mm (maximum patient resolution), confirming the capability of the method to regionally compensate for motion. A novel ROI-based motion model is proposed as an integral part of an envisioned clinical MRI-guided workflow aiming at enhanced image guidance compared to conventional strategies.

MRI-guidance for motion management in external beam radiotherapy: current status and future challenges

High precision conformal radiotherapy requires sophisticated imaging techniques to aid in target localisation for planning and treatment, particularly when organ motion due to respiration is involved. X-ray based imaging is a well-established standard for radiotherapy treatments. Over the last few years, the ability of magnetic resonance imaging (MRI) to provide radiation-free images with high-resolution and superb soft tissue contrast has highlighted the potential of this imaging modality for radiotherapy treatment planning and motion management. In addition, these advantageous properties motivated several recent developments towards combined MRI radiation therapy treatment units, enabling in-room MRI-guidance and treatment adaptation. The aim of this review is to provide an overview of the state-of-the-art in MRI-based image guidance for organ motion management in external beam radiotherapy. Methodological aspects of MRI for organ motion management are reviewed and their application in treatment planning, in-room guidance and adaptive radiotherapy described. Finally, a roadmap for an optimal use of MRI-guidance is highlighted and future challenges are discussed.

The radiobiological impact of motion tracking of liver, pancreas and kidney SBRT tumors in a MR-linac

The purpose of this work is to evaluate and quantify the potential radiobiological advantages of tumor tracking using the MR-linac for three disease sites: liver, pancreas and kidney. From each disease site, three patients were selected and 4DCT data sets were used. We applied two planning methods using the Monaco treatment planning system (Elekta AB,Stockholm,Sweden): (1) the conventional ITV method using a 6MV Agility beam and (2) a simulated tracking method using MLC GTV tracking with a 7MV MR-linac beam model incorporating a 1.5 T transverse magnetic field. A 5 mm isotropic PTV margin was added to the ITV or the GTV, and 95% of the PTV volume received 100% of the prescription dose. To evaluate the potential radiobiological advantages of tumor tracking, the normal tissue complication probabilities (NTCPs) were calculated for each organ at risk (OAR) using the Layman Kutcher Burman (LKB) model. The average reduction in the target volume, due to tracking, was 31.1%, 26.3% and 26.9% for liver, pancreas and kidney patients, respectively. For each OAR, the % differences in NTCP between the two methods were calculated. The mean 2 Gy equivalent OAR dose for all patients was less than 29.1 Gy, below which the NTCP for most OARs was not sensitive to equivalent uniform dose (EUD). As a result, a NTCP benefit, due to tracking, was observed in 26% of the data. For all three disease sites, the maximum NTCP improvements were for the normal kidney, the bowels, and the duodenum, with reductions in associated toxicities of 79% (radiation nephropathy), 69% (stricture/fistula) and 25% (ulceration), respectively. This study demonstrates the potential benefit of using a MR-linac tracking system to reduce NTCPs. The normal kidney, the bowels and the duodenum showed the largest NTCP improvements. This, in part, is due to the rapid changes in NTCP for small EUD changes.

Predicting tumour motion during the whole radiotherapy treatment: a systematic approach for thoracic and abdominal lesions based on real time MR

INTRODUCTION

Aim of this study was to investigate the ability of pre-treatment four dimensional computed tomography (4DCT) to capture respiratory-motion observed in thoracic and abdominal lesions during treatment. Treatment motion was acquired using full-treatment cine-MR acquisitions. Results of this analysis were compared to the ability of 30 seconds (s) cine Magnetic Resonance (MR) to estimate the same parameters.

METHODS

A 4DCT and 30 s cine-MR (ViewRay, USA) were acquired on the simulation day for 7 thoracic and 13 abdominal lesions. Mean amplitude, intra- and inter-fraction amplitude variability, and baseline drift were extracted from the full treatment data acquired by 2D cine-MR, and correlated to the motion on pre-treatment 30 s cine-MR and 4DCT. Using the full treatment data, safety margins on the ITV, necessary to account for all motion variability from 4DCT observed during treatment, were calculated. Mean treatment amplitudes were 2 ± 1 mm and 5 ± 3 mm in the anteroposterior (AP) and craniocaudal (CC) direction, respectively. Differences between mean amplitude during treatment and amplitude on 4DCT or during 30 s cine-MR were not significant, but 30 s cine-MR was more accurate than 4DCT. Intra-fraction amplitude variability was positively correlated with both 30 s cine-MR and 4DCT amplitude. Inter-fraction amplitude variability was minimal.

RESULTS

Mean baseline drift over all fractions and patients equalled 1 ± 1 mm in both CC and AP direction, but drifts per fraction up to 16 mm (CC) and 12 mm (AP) were observed. Margins necessary on the ITV ranged from 0 to 8 mm in CC and 0 to 5 mm in AP direction. Neither amplitude on 4DCT nor during 30 s cine MR is correlated to the magnitude of drift or the necessary margins in both directions.

CONCLUSION

Lesions moving with small amplitude show limited amplitude variability throughout treatment, making passive motion management strategies seem adequate. However, other variations such as baseline drifts and shifts still cause significant geometrical uncertainty, favouring real-time monitoring and an active approach for all lesions influenced by respiratory motion.