Retrospective self-sorted 4D-MRI for the liver

The Use of Magnetic Resonance Imaging in Radiation Therapy Treatment Simulation and Planning

Ever since its introduction as a diagnostic imaging tool the potential of magnetic resonance imaging (MRI) in radiation therapy (RT) treatment simulation and planning has been recognized. Recent technical advances have addressed many of the impediments to use of this technology and as a result have resulted in rapid and growing adoption of MRI in RT. The purpose of this article is to provide a broad review of the multiple uses of MR in the RT treatment simulation and planning process, identify several of the most used clinical scenarios in which MR is integral to the simulation and planning process, highlight existing limitations and provide multiple unmet needs thereby highlighting opportunities for the diagnostic MR imaging community to contribute and collaborate with our oncology colleagues. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 5.

MRI-Guided Radiation Therapy Systems

MR-Guided Radiation Therapy (MRIgRT) has been made possible only due to the ingenuity and commitment of commercial radiation therapy system vendors. Unlike conventional linear accelerator systems, MRIgRT systems have had to overcome significant and previously untested techniques to integrate the MRI systems with the radiation therapy delivery systems. Each of these three commercial systems has developed different approaches to integrating their MR and Linac functions. Each has also decided on a different main magnetic field strength, from 0.35T to 1.5T, as well as different design philosophies for other systems, such as the patient support assembly and treatment planning workflow. This paper is intended to provide the reader with a detailed understanding of each system's configuration so that the reader can better interpret the scientific literature concerning these commercial MRIgRT systems.

Respiratory motion management using a single rapid MRI scan for a 0.35 T MRI-Linac system

BACKGROUND

MRI has a rapidly growing role in radiation therapy (RT) for treatment planning, real-time image guidance, and beam gating (e.g., MRI-Linac). Free-breathing 4D-MRI is desirable in respiratory motion management for therapy. Moreover, high-quality 3D-MRIs without motion artifacts are needed to delineate lesions. Existing MRI methods require multiple scans with lengthy acquisition times or are limited by low spatial resolution, contrast, and signal-to-noise ratio.

PURPOSE

We developed a novel method to obtain motion-resolved 4D-MRIs and motion-integrated 3D-MRI reconstruction using a single rapid (35-45 s scan on a 0.35 T MRI-Linac.

METHODS

Golden-angle radial stack-of-stars MRI scans were acquired from a respiratory motion phantom and 12 healthy volunteers (n = 12) on a 0.35 T MRI-Linac. A self-navigated method was employed to detect respiratory motion using 2000 (acquisition time = 5-7 min) and the first 200 spokes (acquisition time = 35-45 s). Multi-coil non-uniform fast Fourier transform (MCNUFFT), compressed sensing (CS), and deep-learning Phase2Phase (P2P) methods were employed to reconstruct motion-resolved 4D-MRI using 2000 spokes (MCNUFFT2000) and 200 spokes (CS200 and P2P200). Deformable motion vector fields (MVFs) were computed from the 4D-MRIs and used to reconstruct motion-corrected 3D-MRIs with the MOtion Transformation Integrated forward-Fourier (MOTIF) method. Image quality was evaluated quantitatively using the structural similarity index measure (SSIM) and the root mean square error (RMSE), and qualitatively in a blinded radiological review.

RESULTS

Evaluation using the respiratory motion phantom experiment showed that the proposed method reversed the effects of motion blurring and restored edge sharpness. In the human study, P2P200 had smaller inaccuracy in MVFs estimation than CS200. P2P200 had significantly greater SSIMs (p < 0.0001) and smaller RMSEs (p < 0.001) than CS200 in motion-resolved 4D-MRI and motion-corrected 3D-MRI. The radiological review found that MOTIF 3D-MRIs using MCNUFFT2000 exhibited the highest image quality (scoring > 8 out of 10), followed by P2P200 (scoring > 5 out of 10), and then motion-uncorrected (scoring < 3 out of 10) in sharpness, contrast, and artifact-freeness.

CONCLUSIONS

We have successfully demonstrated a method for respiratory motion management for MRI-guided RT. The method integrated self-navigated respiratory motion detection, deep-learning P2P 4D-MRI reconstruction, and a motion integrated reconstruction (MOTIF) for 3D-MRI using a single rapid MRI scan (35-45 s) on a 0.35 T MRI-Linac system.

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.

Contrast-enhanced 4D-MRI for internal target volume generation in treatment planning for liver tumors

BACKGROUND

Liver tumors are often invisible on four-dimensional commuted tomography (4D-CT). Imperfect imaging surrogates are used to estimate the tumor motion. Here, we assessed multiple 4D magnetic resonance (MR) binning algorithms for directly visualizing liver tumor motion for radiotherapy planning.

METHODS

Patients were simulated using a 3 Tesla MR and CT scanner. Three prototype binning algorithms (phase, amplitude, and two-directional) were applied to the 4D-MRIs, and the image quality was assessed using a qualitative clarity score and quantitative sharpness score. Radiation plans were generated for internal target volumes (ITVs) derived using 4D-MRI and 4D-CT, and the dosimetry of targets were compared. Paired t-tests were used to compare sharpness scores and dosimetric data.

RESULTS

Twelve patients with 17 liver tumors were scanned between May and November 2021. Compared to phase binning, two-directional demonstrated equal or better clarity and sharpness scores (end-expiration: 0.33 vs. 0.38, p=0.018, end-inspiration: 0.28 vs. 0.31, p=0.010). Compared to amplitude binning, two-directional binning captured hysteresis of ≥3 mm in 35% of patients. Evaluation of dosimetry CT-optimized plans revealed that PTV coverage of MR-derived targets were significantly lower than CT-derived targets (PTV receiving 90% of prescription: 75.56% vs. 89.38%, p=0.002).

CONCLUSION

Using contrast-enhanced 4D-MRI is feasible for directly delineating liver tumors throughout the respiratory cycle. The current standard of using radiation plans optimized for 4D-CT-derived targets achieved lower coverage of directly visualized MRI targets, suggesting that adopting MRI for motion management may improve radiation treatment of liver lesions and reduce the risk of marginal misses.

Technical feasibility and clinical evaluation of 4D-MRI guided liver SBRT on the MR-linac

PURPOSE

Image-guided stereotactic body radiation therapy (SBRT) is an important local treatment for liver metastases. MRI-guidance enables direct tumor visualization, eliminating fiducial marker implantation. The purpose of this study was to test technical feasibility of our 4D-MRI guided liver SBRT workflow. Additionally, intra-fraction target motion and consequent target-coverage were studied.

MATERIALS&METHODS

Patients with liver metastases were included in this sub-study of the prospective UMBRELLA clinical trial. Patients received mid-position (midP) SBRT. The daily adapt-to-position workflow included localization, verification and intra-fraction tumor midP monitoring using 4D-MRI. Technical feasibility was established based on persistence of the treatment protocol, treatment time ≤1 hour, no geographical miss and no unexpected acute toxicity grade >3. All 4D-MRIs were registered to the planning midP-CT and tumor midP and amplitude were calculated. Additionally, delivered target dose was accumulated incorporating the 4D-MRI intra-fraction tumor motion and evaluated with Monte-Carlo error simulations.

RESULTS

20 patients with liver metastases were included and treated with 4D-MRI guided SBRT. Feasibility criteria were met in all-but-one patient. No grade ≥3 acute toxicity was observed. Group mean (M), systematic and random midP-drifts were 2.4mm, 2.6mm and 3.1mm in CC-direction. 4D-MRI tumor CC-amplitudes were reduced compared to the simulation 4D-CT (M=-1.9mm) and decreased during treatment (M=-1.4mm). Dose accumulation showed adquate target-coverage on a population level.

CONCLUSION

We successfully demonstrated technical feasibility of 4D-MRI guided SBRT in a cohort of 20 patients with liver metastases. However, substantial midposition drifts occurred which stress the need for intra-fraction motion management strategies to further increase the precision of treatment delivery.

Master protocol trial design for technical feasibility of MR-guided radiotherapy

The master protocol trial design aims to increase efficiency in terms of trial infrastructure and protocol administration which may accelerate development of (technical) innovations in radiation oncology. A master protocol to study feasibility of techniques/software for MR-guided adaptive radiotherapy with the MR-Linac is described and discussed.

Evaluating differences in respiratory motion estimates during radiotherapy: a single planning 4DMRI versus daily 4DMRI

Background

In radiotherapy, respiratory-induced tumor motion is typically measured using a single four-dimensional computed tomography acquisition (4DCT). Irregular breathing leads to inaccurate motion estimates, potentially resulting in undertreatment of the tumor and unnecessary dose to healthy tissue. The aim of the research was to determine if a daily pre-treatment 4DMRI-strategy led to a significantly improved motion estimate compared to single planning 4DMRI (with or without outlier rejection).

Methods

4DMRI data sets from 10 healthy volunteers were acquired. The first acquisition simulated a planning MRI, the respiratory motion estimate (constructed from the respiratory signal, i.e. the 1D navigator) was compared to the respiratory signal in the subsequent scans (simulating 5–29 treatment fractions). The same procedure was performed using the first acquisition of each day as an estimate for the subsequent acquisitions that day (2 per day, 4–20 per volunteer), simulating a daily MRI strategy. This was done for three outlier strategies: no outlier rejection (NoOR); excluding 5% of the respiratory signal whilst minimizing the range (Min95) and excluding the datapoints outside the mean end-inhalation and end-exhalation positions (MeanIE).

Results

The planning MRI median motion estimates were 27 mm for NoOR, 18 mm for Min95, and 13 mm for MeanIE. The daily MRI median motion estimates were 29 mm for NoOR, 19 mm for Min95 and 15 mm for MeanIE. The percentage of time outside the motion estimate were for the planning MRI: 2%, 10% and 32% for NoOR, Min95 and MeanIE respectively. These values were reduced with the daily MRI strategy: 0%, 6% and 17%. Applying Min95 accounted for a 30% decrease in motion estimate compared to NoOR.

Conclusion

A daily MRI improved the estimation of respiratory motion as compared to a single 4D (planning) MRI significantly. Combining the Min95 technique with a daily 4DMRI resulted in a decrease of inclusion time of 6% with a 30% decrease of motion. Outlier rejection alone on a planning MRI often led to underestimation of the movement and could potentially lead to an underdosage.

Trial registration: protocol W15_373#16.007

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-021-01915-1.

Validation of a 4D-MRI guided liver stereotactic body radiation therapy strategy for implementation on the MR-linac

Purpose. Accurate tumor localization for image-guided liver stereotactic body radiation therapy (SBRT) is challenging due to respiratory motion and poor tumor visibility on conventional x-ray based images. Novel integrated MRI and radiotherapy systems enable direct in-room tumor visualization, potentially increasing treatment accuracy. As these systems currently do not provide a 4D image-guided radiotherapy strategy, we developed a 4D-MRI guided liver SBRT workflow and validated all steps for implementation on the Unity MR-linac.Materials and Methods. The proposed workflow consists of five steps: (1) acquisition of a daily 4D-MRI scan, (2) 4D-MRI to mid-position planning-CT rigid tumor registration, (3) calculation of daily tumor midP misalignment, (4) plan adaptation using adapt-to-position (ATP) with segment-weights optimization and (5) adapted plan delivery. The workflow was first validated in a motion phantom, performing regular motion at different baselines (±5 to ±10 mm) and patient-derived respiratory signals with varying degrees of irregularity. 4D-MRI derived respiratory signals and 4D-MRI to planning CT registrations were compared to the phantom input, and gamma and dose-area-histogram analyses were performed on the delivered dose distributions on film. Additionally, 4D-MRI to CT registration performance was evaluated in patient images using the full-circle method (transitivity analysis). Plan adaption was further analyzedin-silicoby creating adapted treatment plans for 15 patients with oligometastatic liver disease.Results. Phantom trajectories could be reliably extracted from 4D-MRI scans and 4D-MRI to CT registration showed submillimeter accuracy. The DAH-analysis demonstrated excellent coverage of the dose evaluation structures GTV and GTV_TD. The median daily rigid 4D-MRI to midP-CT registration precision in patient images was <2 mm. The ATP strategy restored the target dose without increased exposure to the OARs and plan quality was independent from 3D shift distance in the range of 1-26 mm.Conclusions. The proposed 4D-MRI guided strategy showed excellent performance in all workflow tests in preparation of the clinical introduction on the Unity MR-linac.

Simultaneous multi-slice accelerated 4D-MRI for radiotherapy guidance

4D-MRI is becoming increasingly important for daily guidance of thoracic and abdominal radiotherapy. This study exploits the simultaneous multi-slice (SMS) technique to accelerate the acquisition of a balanced turbo field echo (bTFE) and a turbo spin echo (TSE) coronal 4D-MRI sequence performed on 1.5 T MRI scanners. SMS single-shot bTFE and TSE sequences were developed to acquire a stack of 52 coronal 2D images over 30 dynamics. Simultaneously excited slices were separated by half the field of view. Slices intersecting with the liver-lung interface were used as navigator slices. For each navigator slice location, an end-exhale dynamic was automatically identified, and used to derive the self-sorting signal by rigidly registering the remaining dynamics. Navigator slices were sorted into 10 amplitude bins, and the temporal relationship of simultaneously excited slices was used to generate sorted 4D-MRIs for 12 healthy volunteers. The self-sorting signal was validated using anin vivopeak-to-peak motion analysis. The smoothness of the liver-lung interface was quantified by comparing to sagittal cine images acquired directly after the SMS-4D-MRI sequence. To ensure compatibility with the MR-linac radiotherapy workflow, the 4D-MRIs were transformed into 3D mid-position (MidP) images using deformable image registration. Consistency of the deformable vector fields was quantified in terms of the distance discordance metric (DDM) in the body. The SMS-4D-TSE sequence was additionally acquired for 3 lung cancer patients to investigate tumor visibility. SMS-4D-MRI acquisition and processing took approximately 7 min. 4D-MRI reconstruction was possible for 26 out of 27 acquired datasets. Missing data in the sorted 4D-MRIs varied from 4%-26% for the volunteers and varied from 8%-24% for the patients. Peak-to-peak (SD) amplitudes analysis agreed within 1.8 (1.1) mm and 0.9 (0.4) mm between the sorted 4D-MRIs and the self-sorting signals of the volunteers and patients, respectively. Liver-lung interface smoothness was found to be in the range of 0.6-3.1 mm for volunteers. The percentage of DDM values smaller than 2 mm was in the range of 85%-89% and 86%-92% for the volunteers and patients, respectively. Lung tumors were clearly visibility in the SMS-4D-TSE images and MidP images. Two fast SMS-accelerated 4D-MRI sequences were developed resulting in T₂/T₁or T₂weighted contrast. The SMS-4D-MRIs and derived 3D MidP-MRIs yielded anatomically plausible images and good tumor visibility. SMS-4D-MRI is therefore a strong candidate to be used for treatment simulation and daily guidance of thoracic and abdominal MR-guided radiotherapy.

Slice-stacking T2-weighted MRI for fast determination of internal target volume for liver tumor

Background

To investigate the feasibility of generating maximum intensity projection (MIP) images to determine internal target volume (ITV) using slice-stacking MRI (SS-MRI) technique.

Methods

Slice-stacking is a technique which applies a multi-slice MRI acquisition to generate a 3D MIP for ITV contouring, without reconstructing 4D-MRI. 4D digital extended cardiac-torso (XCAT) phantom was used to generate MIP images with sequential 2D HASTE sequence, with different tumor diameters (10, 30 and 50 mm) and with simulated regular and irregular (patient) breathing motions. A reference MIP was generated using all acquisition images. Consecutive repetitions were then used to generate MIP to analyze the relationship between Dice's similarity coefficient (DSC) and the number of repetitions, and the relationship between the relative ITV volume difference and the number of repetitions. Images from XCAT phantom and from three hepatic carcinoma patients were collected in this study to demonstrate the feasibility of this technique.

Results

For both regular and irregular breathing motion, the average DSC of ITV is >0.94 and the average relative ITV volume difference is <10% (approximately 0.15 cm³) when using 5 repeated scanning images to reconstruct MIP for tumor diameter of 10 mm. As tumor diameter increases, the DSC of ITV is >0.97 and the relative ITV volume difference is <5% for regular breathing motion, and the DSC of ITV is >0.97 and the relative ITV volume difference is <5.5% for irregular breathing motion when using 5 repeated scanning images to reconstruct MIP. In patient image study, the mean relative ITV volume difference is <3% and the mean DSC is 0.99 when using 5 repeated scanning images to reconstruct MIP.

Conclusions

The number of scans required to generate tumor ITV for slice-stacking method (5-7 repetition) is 3-4 times less than that of 4D-MRI (15-20 repetitions). It is feasible to generate a fast clinically acceptable ITV using slice-stacking method with sequential 2D MR images.

Effects of Gadolinium MRI Contrast Agents on DNA Damage and Cell Survival when Used in Combination with Radiation

Gadolinium is a commonly used contrast agent for magnetic resonance imaging (MRI). The goal of this work was to determine how MRI contrast agents affect radiosensitivity for tumour cells. Using a 225kVp X-ray cabinet source, immunofluorescence and clonogenic assays were performed on six cancer cell lines: lung (H460), pancreas (MiaPaCa2), prostate (DU145), breast (MCF7), brain (U87) and liver (HEPG2). Dotarem® contrast agent, at concentrations of 0.2, 2 and 20 mM, was used to determine its effect on DNA damage and cell survival. Measurements were performed using inductively coupled plasma mass spectrometry (ICP-MS) to determine the amount of gadolinium taken up by each cell line for each concentration. A statistically significant increase in DNA damage was seen for all cell lines at a dose of 1 Gy for concentrations of 2 and 20 mM, at 1 h postirradiation. At 24 h postirradiation, most of the DNA damage had been repaired, with approximately 90% repair for almost all doses of radiation and concentrations of Dotarem. Clonogenic results showed no statistically significant decrease in cell survival for any cell line or concentration. Uptake measurements showed cell line-specific variations in uptake, with MCF7 and HEPG2 cells having a high percentage uptake compared to other cell lines, with 151.4 ± 0.3 × 10-15 g and 194.8 ± 0.4 × 10-15 g per cell, respectively, at 2 mM Dotarem concentration. In this work, a variability in gadolinium uptake was observed between cell lines. A significant increase was seen in initial levels of DNA damage after 1 Gy irradiation for all six cancer cell lines; however, no significant decrease in cell survival was seen with the clonogenic assay. The observation of high levels of repair suggest that while initial levels of DNA damage are increased, this damage is almost entirely repaired within 24 h, and does not affect the ability of cells to survive and produce colonies.

Retrospective four-dimensional magnetic resonance imaging of liver: Method development

Purpose of our research was to develop a four‐dimensional (4D) magnetic resonance imaging (MRI) method of liver. Requirements of the method were to create a clinical procedure with acceptable imaging time and sufficient temporal and spatial accuracy. The method should produce useful planning image sets for stereotactic body radiation therapy delivery both during breath‐hold and in free breathing. The purpose of the method was to improve the localization of liver metastasis. The method was validated with phantom tests. Imaging parameters were optimized to create a 4D dataset compressed to one respiratory cycle of the whole liver with clinically reasonable level of image contrast and artifacts. Five healthy volunteers were imaged with T2‐weighted SSFSE research sequence. The respiratory surrogate signal was observed by the linear navigator interleaved with the anatomical liver images. The navigator was set on head‐feet — direction on the superior surface of the liver to detect the edge of diaphragm. The navigator signal and 2D liver image data were retrospectively processed with a self‐developed MATLAB algorithm. A deformable phantom for 4D imaging tests was constructed by combining deformable tissue‐equivalent material and a commercial programmable motor unit of the 4D phantom with a clinically relevant range of deformation patterns. 4D Computed Tomography images were used as reference to validate the MRI protocol. The best compromise of reasonable accuracy and imaging time was found with 2D T2‐weighted SSFSE imaging sequence using parameters: TR = 500–550 ms, images/slices = 20, slice thickness = 3 mm. Then, image processing with number of respiratory phases = 8 constructed accurate 4D images of liver. We have developed the 4D‐MRI method visualizing liver motions three‐dimensionally in one representative respiratory cycle. From phantom tests it was found that the spatial agreement to 4D‐CT is within 2 mm that is considered sufficient for clinical applications.

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.

MRI basics for radiation oncologists

•

Issues of MRI that are relevant for radiation oncologists are addressed.

•

Radiation oncology requires dedicated scan protocols.

•

Use of diagnostic protocols is not recommended for radiotherapy.

•

MR images must be made in treatment position with the standard positioning devices.

•

Safety screening prior to entering the MRI room is crucial.

MRI is increasingly used in radiation oncology to facilitate tumor and organ-at-risk delineation and image guidance. In this review, we address issues of MRI that are relevant for radiation oncologists when interpreting MR images offered for radiotherapy. Whether MRI is used in combination with CT or in an MRI-only workflow, it is generally necessary to ensure that MR images are acquired in treatment position, using the positioning and fixation devices that are commonly applied in radiotherapy. For target delineation, often a series of separate image sets are used with distinct image contrasts, acquired within a single exam. MR images can suffer from image distortions. While this can be avoided with dedicated scan protocols, in a diagnostic setting geometrical fidelity is less relevant and is therefore less accounted for. Since geometrical fidelity is of utmost importance in radiation oncology, it requires dedicated scan protocols. The strong magnetic field of an MRI scanner and the use of radiofrequency radiation can cause safety hazards if not properly addressed. Safety screening is crucial for every patient and every operator prior to entering the MRI room.

A fast volumetric 4D-MRI with sub-second frame rate for abdominal motion monitoring and characterization in MRI-guided radiotherapy

Background

To propose a fast volumetric 4D-MRI based on 3D pulse sequence acquisition for abdominal motion monitoring and characterization in MRI-guided radiotherapy (MRgRT).

Methods

A 3D spoiled gradient echo sequence volumetric interpolated breath-hold examination (VIBE) [repetition time/echo time (TR/TE) =0.53/1.57 ms, flip-angle =5°, receiver bandwidth (RBW) =1,400 Hz/voxel] based 4D-MRI acquisition, accelerated by 4-fold controlled aliasing in parallel imaging results in higher acceleration (CAIPIRINHA), named CAIPIRINHA-VIBE 4D-MRI, was implemented on a 1.5T MRI simulator (MR-sim) and applied for abdominal imaging of nine healthy volunteers under free breathing. One hundred and forty-four dynamics of the entire abdomen volume (56 slices), in total 8,064 (144×56) images with a voxel size of 2.7×2.7×4.0 mm³, were acquired in 89 s for 4D-MRI. This CAIPIRINHA-VIBE 4D-MRI was qualitatively compared with a 2D half-Fourier acquisition single-shot turbo spin-echo (2D-HASTE) based 4D-MRI. The motions of liver dome, kidney and spleen were analyzed using the CAIPIRINHA-VIBE 4D-MRI data. The kidney motion was quantitatively characterized in terms of motion range and the correlations between left and right kidneys.

Results

CAIPIRINHA-VIBE 4D-MRI was successfully conducted in all subjects. CAIPIRINHA-VIBE 4D-MRI exhibited much higher effective volumetric temporal resolution (0.615 vs. ~5 s/volume) and better reconstructed volume consistency than 2D-HASTE 4D-MRI. CAIPIRINHA-VIBE 4D-MRI was able to characterize the respiratory motion of abdominal organs simultaneously in three orthogonal directions, and could potentially be used for whole abdomen deformable motion tracking. Renal motion range was most pronounced in superior-inferior (SI) direction (L: 10.03±2.65 mm; R: 10.38±2.80 mm), significantly larger (P<0.001) than that in anterior-posterior (AP) and the least in left-right (LR) directions. Right kidney had significantly larger mobility (4.18±2.19 vs. 2.32±1.34 mm, P=0.045) than left kidney in AP, but not in LR and SI directions. The Pearson correlation coefficients r between left and right kidney motion were 0.5063 (P=0.164), 0.6624 (P=0.052) and 0.5752 (P=0.105) in LR, AP and SI correspondingly. The correlation of renal motion in SI and AP was found significant in right kidney (r=0.843, P=0.004) but not in left kidney (r=0.467, P=0.205).

Conclusions

A fast volumetric 4D-MRI was implemented for abdominal motion monitoring in MRgRT. A sub-second volumetric temporal resolution of 0.615 s, covering the entire abdomen, was demonstrated for respiratory motion monitoring and characterization. This technique holds potentials for MRgRT applications.

A novel amplitude binning strategy to handle irregular breathing during 4DMRI acquisition: improved imaging for radiotherapy purposes

Background

For radiotherapy of abdominal cancer, four-dimensional magnetic resonance imaging (4DMRI) is desirable for tumor definition and the assessment of tumor and organ motion. However, irregular breathing gives rise to image artifacts. We developed a outlier rejection strategy resulting in a 4DMRI with reduced image artifacts in the presence of irregular breathing.

Methods

We obtained 2D T2-weighted single-shot turbo spin echo images, with an interleaved 1D navigator acquisition to obtain the respiratory signal during free breathing imaging in 2 patients and 12 healthy volunteers. Prior to binning, upper and lower inclusion thresholds were chosen such that 95% of the acquired images were included, while minimizing the distance between the thresholds (inclusion range (IR)). We compared our strategy (Min95) with three commonly applied strategies: phase binning with all images included (Phase), amplitude binning with all images included (MaxIE), and amplitude binning with the thresholds set as the mean end-inhale and mean end-exhale diaphragm positions (MeanIE). We compared 4DMRI quality based on:

Data included (DI); percentage of images remaining after outlier rejection.

Reconstruction completeness (RC); percentage of bin-slice combinations containing at least one image after binning.

Intra-bin variation (IBV); interquartile range of the diaphragm position within the bin-slice combination, averaged over three central slices and ten respiratory bins.

IR.

Image smoothness (S); quantified by fitting a parabola to the diaphragm profile in a sagittal plane of the reconstructed 4DMRI.

A two-sided Wilcoxon’s signed-rank test was used to test for significance in differences between the Min95 strategy and the Phase, MaxIE, and MeanIE strategies.

Results

Based on the fourteen subjects, the Min95 binning strategy outperformed the other strategies with a mean RC of 95.5%, mean IBV of 1.6 mm, mean IR of 15.1 mm and a mean S of 0.90. The Phase strategy showed a poor mean IBV of 6.2 mm and the MaxIE strategy showed a poor mean RC of 85.6%, resulting in image artifacts (mean S of 0.76). The MeanIE strategy demonstrated a mean DI of 85.6%.

Conclusions

Our Min95 reconstruction strategy resulted in a 4DMRI with less artifacts and more precise diaphragm position reconstruction compared to the other strategies.

Trial registration

Volunteers: protocol W15_373#16.007; patients: protocol NL47713.018.14

Electronic supplementary material

The online version of this article (10.1186/s13014-019-1279-z) contains supplementary material, which is available to authorized users.

Delineation of a Cardiac Planning Organ-At-Risk Volume Using Real-Time Magnetic Resonance Imaging for Cardiac Protection in Thoracic and Breast Radiation Therapy

PURPOSE

Cardiac radiation is associated with cardiotoxicity in patients with thoracic and breast malignancies. We conducted a prospective study using cine magnetic resonance imaging (MRI) scans to evaluate heart motion. We hypothesized that cine MRI could be used to define population-based cardiac planning organ-at-risk volumes (PRV).

METHODS AND MATERIALS

A total of 16 real-time acquisitions were obtained per subject on a 1.5 Tesla MRI (Philips Ingenia). Planar cine MRI was performed in 4 sequential sagittal and coronal planes at free-breathing (FB) and deep-inspiratory breath hold (DIBH). In-plane cardiac motion was assessed using a scale-invariant feature transformation-based algorithm. Subject-specific pixel motion ranges were defined in anteroposterior (AP), left-right (LR), and superoinferior (SI) planes. Averages of the 98% and 67% of the maximum ranges of pixel displacement were defined by subject, then averaged across the cohort to calculate PRV expansions at FB and DIBH.

RESULTS

Data from 20 subjects with a total of 3120 image frames collected per subject in coronal and sagittal planes at DIBH and FB, and 62,400 total frames were analyzed. Cohort averages of 98% of the maximum cardiac motion ranges comprised margin expansions of 12.5 ± 1.1 mm SI, 5.8 ± 1.2 mm AP, and 6.6 ± 1.0 mm LR at FB and 6.7 ± 1.5 mm SI, 4.7 ± 1.3 mm AP, and 5.3 ± 1.3 mm LR at DIBH. Margins for 67% of the maximum range comprised 7.7 ± 0.7 mm SI, 3.2 ± 0.6 mm AP, and 3.7 ± 0.6 mm LR at FB and 4.1 ± 0.9 mm SI, 2.7 ± 0.8 mm AP, and 3.2 ± 0.8 mm LR at DIBH. Subsequently, these margins were simplified to form PRVs for treatment planning.

CONCLUSIONS

We implemented scale-invariant feature transformation-based motion tracking for analysis of the cardiac cine MRI scans to quantify motion and create cohort-based cardiac PRVs to improve cardioprotection in breast and thoracic radiation.

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.

A Self-Sorting Coronal 4D-MRI Method for Daily Image Guidance of Liver Lesions on an MR-LINAC

PURPOSE

Novel hybrid MR-LINAC devices provide MRI's superior soft-tissue contrast in the treatment room and thus have the potential to increase accuracy of liver stereotactic body radiation therapy (SBRT). Requirements for daily position verification using 4-dimensional MRI include tumor visibility and short acquisition-reconstruction time (preferably <5 min). The proposed method provides fast acquisition-reconstruction time and the flexibility to vary T1- and T2-weighting, using standard imaging sequences for straightforward implementation on an MR-LINAC.

METHODS AND MATERIALS

Images were acquired using a coronal 2-dimensional, multislice, single-shot turbo spin-echo (TSE) and turbo field-echo (TFE) sequence, which were repeated 30 times. An image-based self-sorting signal (ImS) was extracted from the data, and rigid registration of the diaphragm per slice position was performed and corrected for amplitude variation in the anteroposterior direction. Data were sorted into 10 bins according to amplitude and phase. ImS was validated in 4 healthy volunteers against a navigator signal. Positional variations within bins, missing data, and smoothness of the liver dome were compared between amplitude and phase binning in 10 volunteers. Tumor contrast and registration were investigated in 3 patients.

RESULTS

Each ImS was found to be in excellent agreement with the navigator signal with a correlation coefficient of >0.95 and binning differences of <1 bin. Better liver dome smoothness per bin in case of amplitude binning compared with that in phase binning (2.0-2.6 mm vs 2.4-3.7 mm, respectively) is a tradeoff for more missing data (3.5%-17.5% vs 3.5%-4.7%, respectively). Liver lesions were visible in almost all coronal TSE and TFE images, but the lesion boundary was better defined in the TSE images. Rigid registrations could be performed on the tumor area.

CONCLUSIONS

An efficient self-sorted 4-dimensional MRI method was developed and validated using standard sequences and fast reconstruction on a LINAC-integrated MRI scanner providing good tumor visibility for daily image-guided liver stereotactic body radiation therapy.