Prospective evaluation of patient-reported anxiety and experiences with adaptive radiation therapy on an MR-linac

Purpose

An integrated magnetic resonance scanner and linear accelerator (MR-linac) was implemented with daily online adaptive radiation therapy (ART). This study evaluated patient-reported experiences with their overall hospital care as well as treatment in the MR-linac environment.

Methods

Patients pre-screened for MR eligibility and claustrophobia were referred to simulation on a 1.5 T MR-linac. Patient-reported experience measures were captured using two validated surveys. The 15-item MR-anxiety questionnaire (MR-AQ) was administered immediately after the first treatment to rate MR-related anxiety and relaxation. The 40-item satisfaction with cancer care questionnaire rating doctors, radiation therapists, the services and care organization and their outpatient experience was administered immediately after the last treatment using five-point Likert responses. Results were analyzed using descriptive statistics.

Results

205 patients were included in this analysis. Multiple sites were treated across the pelvis and abdomen with a median treatment time per fraction of 46 and 66 min respectively. Patients rated MR-related anxiety as "not at all" (87%), "somewhat" (11%), "moderately" (1%) and "very much so" (1%). Positive satisfaction responses ranged from 78 to 100% (median 93%) across all items. All radiation therapist-specific items were rated positively as 96-100%. The five lowest rated items (range 78-85%) were related to general provision of information, coordination, and communication. Overall hospital care was rated positively at 99%.

Conclusion

In this large, single-institution prospective cohort, all patients had low MR-related anxiety and completed treatment as planned despite lengthy ART treatments with the MR-linac. Patients overall were highly satisfied with their cancer care involving ART using an MR-linac.

Clinical Outcomes of Patients with Metastatic Breast Cancer Treated with Hypo-Fractionated Liver Radiotherapy

PURPOSE

To retrospectively review the clinical outcomes of patients with metastatic breast cancer (MBCa) following liver directed ablative intent radiotherapy (RT).

METHODS

Demographics, disease and treatment characteristics of patients with MBCa who received liver metastasis (LM) directed ablative RT between 2004-2020 were analysed. The primary outcome was local control (LC), secondary outcomes included overall survival (OS) and progression-free survival (PFS) analyzed by univariate (UVA) and multi-variable analysis (MVA).

RESULTS

Thirty MBCa patients with 50 LM treated with 5-10 fraction RT were identified. Median follow-up was 14.6 (range 0.9-156.2) months. Class of metastatic disease was described as induced (12 patients, 40%), repeat (15 patients, 50%) and de novo (three patients, 10%). Median size of treated LM was 3.1 cm (range 1-8.8 cm) and median biologically effective dose delivered was 122 (Q1-Q3; 98-174) Gy₃. One-year LC rate was 100%. One year and two-year survival was 89% and 63%, respectively, with size of treated LM predictive of OS (HR 1.35, p = 0.023) on UVA. Patients with induced OMD had a significantly higher rate of progression (HR 4.77, p = 0.01) on UVA, trending to significance on MVA (HR 3.23, p = 0.051).

CONCLUSIONS

Hypo-fractionated ablative liver RT in patients with MBCa provides safe, tolerable treatment with excellent LC.

The future of MRI in radiation therapy: Challenges and opportunities for the MR community

Radiation therapy is a major component of cancer treatment pathways worldwide. The main aim of this treatment is to achieve tumor control through the delivery of ionizing radiation while preserving healthy tissues for minimal radiation toxicity. Because radiation therapy relies on accurate localization of the target and surrounding tissues, imaging plays a crucial role throughout the treatment chain. In the treatment planning phase, radiological images are essential for defining target volumes and organs-at-risk, as well as providing elemental composition (e.g., electron density) information for radiation dose calculations. At treatment, onboard imaging informs patient setup and could be used to guide radiation dose placement for sites affected by motion. Imaging is also an important tool for treatment response assessment and treatment plan adaptation. MRI, with its excellent soft tissue contrast and capacity to probe functional tissue properties, holds great untapped potential for transforming treatment paradigms in radiation therapy. The MR in Radiation Therapy ISMRM Study Group was established to provide a forum within the MR community to discuss the unmet needs and fuel opportunities for further advancement of MRI for radiation therapy applications. During the summer of 2021, the study group organized its first virtual workshop, attended by a diverse international group of clinicians, scientists, and clinical physicists, to explore our predictions for the future of MRI in radiation therapy for the next 25 years. This article reviews the main findings from the event and considers the opportunities and challenges of reaching our vision for the future in this expanding field.

MRI-Guided Online Adaptive Stereotactic Body Radiation Therapy of Liver and Pancreas Tumors on an MR-Linac System

Simple Summary

The hybrid magnetic resonance imaging and medical linear accelerator (MR-Linac) systems are expected to revolutionize radiation therapy, uniquely offering high quality soft-tissue contrast and fast imaging to facilitate the online re-planning and guidance of the treatment delivery. While the clinical procedures for stereotactic body radiotherapy are well-established for conventional linacs (with their strengths and weaknesses), they still require significant development and refinement for the MR-Linac systems. Adjustment of fractionation schemes including clinical goals, patient selection, organ motion management, treatment session length, staff logistics, and overall complexity of the online re-planning sessions are key factors that drive the adoption of MR-Linac technologies. In this report, we present the clinical implementation of an MRI-guided radiation therapy workflow, which was used to treat 16 upper gastro-intestinal cancer patients on a 1.5 T MR-Linac platform. The workflow was proven to be feasible for a wide range of clinical scenarios, and the overall treatment session time was significantly reduced as tasks were optimized and the clinical team gradually gained expertise.

Abstract

Purpose: To describe a comprehensive workflow for MRI-guided online adaptive stereotactic body radiation therapy (SBRT) specific to upper gastrointestinal cancer patients with abdominal compression on a 1.5T MR-Linac system. Additionally, we discuss the workflow’s clinical feasibility and early experience in the case of 16 liver and pancreas patients. Methods: Eleven patients with liver cancer and five patients with pancreas cancer were treated with online adaptive MRI-guidance under abdominal compression. Two liver patients received single-fraction treatments; the remainder plus all pancreas cancer patients received five fractions. A total of 65 treatment sessions were investigated to provide analytics relevant to the online adaptive processes. The quantification of target and organ motion as well as definition and validation of internal target volume (ITV) margins were performed via multi-contrast imaging provided by three different 2D cine sequences. The plan generation was driven by full re-optimization strategies and using T2-weighted 3D image series acquired by means of a respiratory-triggered exhale phase or a time-averaged imaging protocol. As a pre-requisite for the clinical development of the procedure, the image quality was thoroughly investigated via phantom measurements and numerical simulations specific to upper abdominal sites. The delivery of the online adaptive treatments was facilitated by real-time monitoring with 2D cine imaging. Results: Liver 1-fraction and 5-fraction online adaptive session time were on average 80 and 67.5 min, respectively. The total session length varied between 70–90 min for a single fraction and 55–90 min for five fractions. The pancreas sessions were 54–85 min long with an average session time of 68.2 min. Target visualization on the 2D cine image data varied per patient, with at least one of the 2D cine sequences providing sufficient contrast to confidently identify its location and confirm reproducibility of ITV margins. The mean/range of absolute and relative dose values for all treatment sessions evaluated with ArcCheck were 90.6/80.9–96.1% and 99/95.4–100%, respectively. Conclusion: MR-guidance is feasible for liver and pancreas tumors when abdominal compression is used to reduce organ motion, improve imaging quality, and achieve a robust intra- and inter-fraction patient setup. However, the treatment length is significantly longer than for the conventional linac, and patient compliance is paramount for the successful completion of the treatment. Opportunities for reducing the online adaptive session time should be explored. As the next steps, dose-of-the-day and dose accumulation analysis and tools are needed to enhance the workflow and to help further refine the online re-planning processes.

Quantification of magnetic susceptibility fingerprint of a 3D linearity medical device

PURPOSE

The study investigates the numerical modelling as well as experimental validation of magnetic susceptibility effects with respect to a 3D linearity phantom used for the quantification of MR image distortions.

METHODS

Magnetic field numerical simulations based on finite difference methods were conducted to generate the susceptibility (χ) model of the MRID^3D phantom. Experimental data was acquired and analyzed for eight different MR scanners to include a wide range of scanning parameters. Distortion vector fields were generated by applying a harmonic analysis based on finite elements methods. Phantom scans for the same setup but with opposite polarities of the frequency encoding gradient were processed in conjunction with the susceptibility modelling to separately quantify three field components due to gradient non-linearities (GNL), B₀ inhomogeneities and χ perturbations.

RESULTS

The numerical modelling showed a significant range of χ value of up to 8.23 ppm, with a mean value of 2.9 ppm. The χ perturbations were found to be mostly present at the end plates of the cylindrical phantom design. The simulations also showed that setup rotations of up to 10° introduced only negligible variations in the χ model of less than 0.1 ppm. This allows for a straightforward practical implementation of the modelling as a single lookup table. After correcting for the χ perturbations, the B₀ inhomogeneities were derived and found to be in good agreement with either the MR system manufacturer specifications or experimental data available in the literature.

CONCLUSIONS

It is possible to accurately model the magnetic susceptibility signature of a 3D linearity device and remove it as a post-processing correction step. This is important as the procedure unlocks the ability of determining both the GNL field and B₀ map of the scanner without the need of extra acquisitions or phantoms.

MRI evaluation of normal tissue deformation and breathing motion under an abdominal compression device

Purpose

Abdominal compression can minimize breathing motion in stereotactic radiotherapy, though it may impact the positioning of dose‐limiting normal tissues. This study quantified the reproducibility of abdominal normal tissues and respiratory motion with the use of an abdominal compression device using MR imaging.

Methods

Twenty healthy volunteers had repeat MR over 3 days under an abdominal compression plate device. Normal tissues were delineated on daily axial T2‐weighted MR and compared on days 2 and 3 relative to day 1, after adjusting for baseline shifts relative to bony anatomy. Inter‐fraction organ deformation was computed using deformable registration of axial T2 images. Deformation > 5 mm was assumed to be clinically relevant. Inter‐fraction respiratory amplitude changes and intra‐fraction baseline drifts during imaging were quantified on daily orthogonal cine‐MR (70 s each), and changes > 3 mm were assumed to be relevant.

Results

On axial MR, the mean inter‐fraction normal tissue deformation was > 5 mm for all organs (range 5.1–13.4 mm). Inter‐fraction compression device misplacements > 5 mm and changes in stomach volume > 50% occurred at a rate of 93% and 38%, respectively, in one or more directions and were associated with larger adjacent organ deformation, in particular for the duodenum. On cine‐MR, inter‐fraction amplitude changes > 3 mm on day 2 and 3 relative to day 1 occurred at a rate of < 12.5% (mean superior–inferior change was 1.6 mm). Intra‐fraction baseline drifts > 3 mm during any cine‐MR acquisition occurred at a rate of 23% (mean superior–inferior changes was 2.4 mm).

Conclusions

Respiratory motion under abdominal compression is reproducible in most subjects within 3 mm. However, inter‐fraction deformations greater than 5 mm in normal tissues were common and larger than inter‐ and intra‐fraction respiratory changes. Deformations were driven mostly by variable stomach contents and device positioning. The magnitude of this motion may impact normal tissue dosimetry during stereotactic radiotherapy.

Extensive Unpredictable Pancreas Cancer Inter-fraction Motion: A Case Report

We present a case of locally advanced pancreatic cancer with duodenal invasion treated with consolidative chemoradiation, where extensive unpredictable interfraction motion was observed. Initially, two attempts were made to treat with volumetric modulated arc therapy technique. However, due to substantial interfractional motion of the pancreatic head mass relative to the regional nodal areas, the patient was eventually replanned and treated with a four-field box technique. This case highlights the difficulty in delivering conformal radiation to the pancreas and quantifies the movement of the target, the adjacent biliary stent, and regional nodes.

Technical Note: Harmonic analysis applied to MR image distortion fields specific to arbitrarily shaped volumes

PURPOSE

Magnetic resonance imaging is expected to play a more important role in radiation therapy given the recent developments in MR-guided technologies. MR images need to consistently show high spatial accuracy to facilitate RT-specific tasks such as treatment planning and in-room guidance. The present study investigates a new harmonic analysis method for the characterization of complex three-dimensional (3D) fields derived from MR images affected by system-related distortions.

METHODS

An interior Dirichlet problem based on solving the Laplace equation with boundary conditions (BCs) was formulated for the case of a 3D distortion field. The second-order boundary value problem (BVP) was solved using a finite elements method (FEM) for several quadratic geometries - that is, sphere, cylinder, cuboid, D-shaped, and ellipsoid. To stress-test the method and generalize it, the BVP was also solved for more complex surfaces such as a Reuleaux 9-gon and the MR imaging volume of a scanner featuring a high degree of surface irregularities. The BCs were formatted from reference experimental data collected with a linearity phantom featuring a volumetric grid structure. The method was validated by comparing the harmonic analysis results with the corresponding experimental reference fields.

RESULTS

The harmonic fields were found to be in good agreement with the baseline experimental data for all geometries investigated. In the case of quadratic domains, the percentage of sampling points with residual values larger than 1 mm was 0.5% and 0.2% for the axial components and vector magnitude, respectively. For the general case of a domain defined by the available MR imaging field of view, the reference data showed a peak distortion of about 1 mm and 79% of the sampling points carried a distortion magnitude larger than 1 mm (tolerance intrinsic to the experimental data). The upper limits of the residual values after comparison with the harmonic fields showed max and mean of 1.4 and 0.25 mm, respectively, with only 1.5% of sampling points exceeding 1 mm.

CONCLUSIONS

A novel harmonic analysis approach relying on finite element methods was introduced and validated for multiple volumes with surface shape functions ranging from simple to highly complex. Since a boundary value problem is solved the method requires input data from only the surface of the desired domain of interest. It is believed that the harmonic method will facilitate (a) the design of new phantoms dedicated for the quantitation of MR image distortions in large volumes and (b) an integrative approach of combining multiple imaging tests specific to radiotherapy into a single test object for routine imaging quality control.

Investigation of the 4D composite MR image distortion field associated with tumor motion for MR-guided radiotherapy

PURPOSE

Magnetic resonance (MR) images are affected by geometric distortions due to the specifics of the MR scanner and patient anatomy. Quantifying the distortions associated with mobile tumors is particularly challenging due to real anatomical changes in the tumor's volume, shape, and relative location within the MR imaging volume. In this study, the authors investigate the 4D composite distortion field, which combines the effects of the susceptibility-induced and system-related distortion fields, experienced by mobile lung tumors.

METHODS

The susceptibility (χ) effects were numerically simulated for two specific scenarios: (a) a full motion cycle of a lung tumor due to breathing as depicted on ten phases of a 4D CBCT data set and (b) varying the tumor size and location in lung tissue via a synthetically generated sphere with variable diameter (4-80 mm). The χ simulation procedure relied on the segmentation and generation of 3D susceptibility (χ) masks and computation of the magnetic field by means of finite difference methods. A system-related distortion field, determined with a phantom and image processing algorithm, was used as a reference. The 4D composite distortion field was generated as the vector summation of the χ-induced and system-related fields. The analysis was performed for two orientations of the main magnetic field (B0), which correspond to several MRIgRT system configurations. Specifically, B0 was set along the z-axis as in the case of a cylindrical-bore scanner and in the (x,y)-plane as for a biplanar MR. Computations were also performed for a full revolution at 15° increments in the case of a rotating biplanar magnet. Histograms and metrics such as maximum, mean, and range were used to evaluate the characteristics of the 4D distortion field.

RESULTS

The χ-induced field depends on the change in volume and shape of the moving tumor as well as the local surrounding anatomy. In the case of system-related distortions, the tumor experiences increased field perturbations as it moves further away from the MR isocenter. For a mobile lung tumor, the 4D composite field, corresponding to a 1.5 T field and a readout gradient of 5 mT/m, amounts to 3.0 and 2.8 mm for the MRIgRT system designs featuring B0 oriented along the z-axis (cylindrical-bore scanner) and in the (x,y)-plane (biplanar scanner), respectively. For a rotating biplanar scanner, the composite distortion field varied nonlinearly with the rotation angle. Overall, the dominant contribution to the composite field was from the system-related distortion field. The tumor centroid experienced a systematic shift of 2 mm and showed a negligible perturbation for different B0 values. The dependency on the tumor size was also investigated, namely the max values varied from 1.2 to 2.5 mm for spherical volumes with a diameter between 4 and 80 mm.

CONCLUSIONS

The composite distortion field requires adequate quantification for lung radiation therapy applications such as treatment planning, pretreatment patient setup verification, and real-time treatment delivery guidance. For certain scenarios such as small tumor volumes, the spatial distortions may be corrected by applying systematic shifts derived from a single tumor motion phase. In the case of high readout gradients common to fast imaging applications, the χ distortions were found to be less than 1 mm irrespective of scanner configuration.

SU-D-213CD-06: Workflow and Safety Systems of a Linac-MR Sim-Brachytherapy MRgRT™ Facility

PURPOSE

To develop the operational workflow and safety systems of a magnetic resonance-guided radiotherapy system (MRgRT™), which comprises an MR scanner on rails that travels between a linac vault, MR simulation room and brachytherapy suite.

METHODS

To develop a safe and streamlined clinical workflow, we conducted a comprehensive process review based on a layered approach to overall MRgRT safety that included i) facility design, (ii) workflow iii) system design and interlocks and iv) policies and procedures. We applied existing guidelines for MR and radiation safety, and employed system-level failure modes and effects analyses to design the MRgRT facility and clinical procedures.

RESULTS

In the MRgRT system configuration, the MR and treatment systems are physically decoupled and used independently requiring novel administration of existing MR and radiation guidelines. A key element for the safe operation of the moving MR unit is the concept that all three rooms represent zone 4 areas (American College of Radiology guidelines). Using this concept, we applied MR guidelines to develop safe procedures for the overall suite, including screening of all persons entering the suite in zone 2 and control of ferromagnetic materials. We generated a clinical workflow that ensures expedient and safe transition between MR imaging and treatment delivery in both the linac and brachytherapy rooms. In addition, we designed emergency protocols for MRgRT, which helped drive requirements for the facility and system design, e.g., need for an accessible MR-safe stretcher.

CONCLUSIONS

We designed the first comprehensive description of the MRgRT workflow, interlocking systems and safety procedures. With this layered approach to safety, we addressed critical aspects regarding safe operation and workflow for the system and provided multiple redundancies for key processes. Coupled with customized staff training, the proposed design ensures the safe operation of the MRgRT facility. This work has received research personnel support from IMRIS.

Implications of tissue magnetic susceptibility-related distortion on the rotating magnet in an MR-linac design

PURPOSE

One of the recently published concepts that combine the soft-tissue imaging capabilities of MRI with external beam radiotherapy involves the rigid coupling of a linac with a rotating biplanar low-field MR imaging system. While such a system would prevent possible image distortion resulting from relative motion between the magnet and the linac, the rotation of the magnet around the patient can itself introduce possibilities for image distortion that need to be addressed. While there are straightforward techniques in the literature for correcting distortions from gradient nonlinearities and nonuniform magnetic fields during image reconstruction, the correction of distortions related to tissue magnetic susceptibility is more complex. This work investigates the extent of this latter distortion type under the regime of a rotating magnetic field.

METHODS

CT images covering patient anatomy in the head, lung, and male pelvic regions were obtained and segmented into components of air, bone, and soft tissue. Each of these three components was assigned bulk magnetic susceptibility values in accordance with those found in the literature. A finite-difference algorithm was then implemented to solve for magnetic field distortion maps should the anatomies be placed in the uniform polarizing field of an MR system. The algorithm was repeated multiple times as the polarizing field was rotated axially about the virtual patient in 15 degrees increments. In this way, a map of maximum distortion, and the range of distortion as the magnetic field is rotated about each anatomical region could be determined. The consequence of these susceptibility distortions in terms of geometric signal shift was calculated for 0.2 T, as well as another low-field system (0.5 T), and a higher field 1.5 T system for comparison, using the assumption of a frequency encoding gradient strength of 5 mT/m.

RESULTS

At 0.2 T, the susceptibility-related distortion was limited to less than 0.5 mm given an encoding gradient strength of 5 mT/m or higher. To maintain this same level of geometric accuracy, the 0.5 T system would require a moderately higher minimum gradient strength of 11 mT/m, and at a typical MR field strength of 1.5 T this minimum gradient strength would increase to 33 mT/m. The influence of magnetic susceptibility on mean frequency shift as the field orientation was rotated was also investigated and found to account for less than half a millimeter at 1.5 T, and negligible for low-field systems.

CONCLUSIONS

A study of three sites (head, lung, and prostate) that are vulnerable to magnetic susceptibility-related distortions were studied, and showed that in the context of a rotating polarizing magnet, low-field systems can maintain geometric accuracy of 0.5 mm with at most moderate limitations on sequence parameters. This conclusion will likely apply only to endogenous tissues, as implanted materials such as titanium can create field distortions much in excess of what may normally be induced in the body. Items containing such materials (hip prostheses, for example) will require individual scrutiny.