Simulation and experimental benchmarking of a proton pencil beam scanning nozzle model for development of MR-integrated proton therapy.

BACKGROUND: MR-integrated proton therapy is under development. It consists of the unique challenge of integrating a proton pencil beam scanning (PBS) beam line nozzle with an magnetic resonance imaging (MRI) scanner. The magnetic interaction between these two components is deemed high risk as the MR images can be degraded if there is cross-talk during beam delivery and image acquisition. PURPOSE: To create and benchmark a self-consistent proton PBS nozzle model for empowering the next stages of MR-integrated proton therapy development, namely exploring and de-risking complete integrated prototype system designs including magnetic shielding of the PBS nozzle. MATERIALS AND METHODS: Magnetic field (COMSOL Multiphysics ${\text{Multiphysics}}$ ) and radiation transport (Geant4) models of a proton PBS nozzle located at OncoRay (Dresden, Germany) were developed according to the manufacturers specifications. Geant4 simulations of the PBS process were performed by using magnetic field data generated by the COMSOL Multiphysics ${\text{Multiphysics}}$ simulations. In total 315 spots were simulated which consisted of a 40 × 30 cm 2 $40\times 30\,{\text{cm}}^{2}$ scan pattern with 5 cm spot spacings and for proton energies of 70, 100, 150, 200, and 220 MeV. Analysis of the simulated deflection at the beam isocenter plane was performed to determine the self-consistency of the model. The magnetic fringe field from a sub selection of 24 of the 315 spot simulations were directly compared with high precision magnetometer measurements. These focused on the maximum scanning setting of ± $\pm$  20 cm beam deflection as generated from the second scanning magnet in the PBS for a proton beam energy of 220 MeV. Locations along the beam line central axis (CAX) were measured at beam isocenter and downstream of 22, 47, 72, 97, and 122 cm. Horizontal off-axis positions were measured at 22 cm downstream of isocenter ( ± $\pm$  50, ± $\pm$  100, and ± $\pm$  150 cm from CAX). RESULTS: The proton PBS simulations had good spatial agreement to the theoretical values in all 315 spots examined at the beam line isocenter plane (0-2.9 mm differences or within 1.5 % of the local spot deflection amount). Careful analysis of the experimental measurements were able to isolate the changes in magnetic fields due solely to the scanning magnet contribution, and showed 1.9  ± $\pm$  1.2 µ T $\bf{\mu} {\text{T}}$ -9.4 ± $\pm$  1.2 µ T $\bf{\mu} {\text{T}}$ changes over the range of measurement locations. Direct comparison with the equivalent simulations matched within the measurement apparatus and setup uncertainty in all but one measurement point. CONCLUSIONS: For the first time a robust, accurate and self-consistent model of a proton PBS nozzle assembly has been created and successfully benchmarked for the purposes of advancing MR-integrated proton therapy research. The model will enable confidence in further simulation based work on fully integrated designs including MRI scanners and PBS nozzle magnetic shielding in order to de-risk and realize the full potential of MR-integrated proton therapy.

Patterns of practice of image guided particle therapy for cranio-spinal irradiation: A site specific multi-institutional survey of European Particle Therapy Network.

PURPOSE: To investigate the current practice patterns in image-guided particle therapy (IGPT) for cranio-spinal irradiation (CSI). METHODS: A multi-institutional survey was distributed to European particle therapy centres to analyse all aspects of IGPT. Based on the survey results, a Delphi consensus analysis was developed to define minimum requirements and optimal workflow for clinical practice. The centres participating in the institutional survey were invited to join the Delphi process. RESULTS: Eleven centres participated in the survey. Imaging for treatment planning was rather similar among the centres with Computed Tomography (CT) being the main modality. For positioning verification, 2D IGPT was more commonly used than 3D IGPT. Two centres performed routinely imaging for plan adaptation, by the rest ad hoc. Eight centres participated in the Delphi consensus analysis. The full consensus was reached on the use of CT imaging without contrast for treatment planning and the role of magnetic resonance imaging (MRI) in target and organs-at-risk delineation. There was an agreement on the necessity to perform patient position verification and correction before each isocentre. The most important outcome was the clear need for standardization and harmonization of the workflow. CONCLUSION: There were differences in CSI IGPT clinical practice among the European particle therapy centres. Moreover, the optimal workflow as identified by experts was not yet reached. There is a strong need for consensus guidelines. The state-of-the-art imaging technology and protocols need to be implemented into clinical practice to improve the quality of IGPT for CSI.

Prediction of radiation pneumonitis using the effective α/ß of lungs and heart in NSCLC patients treated with proton beam therapy.

PURPOSE: Radiation pneumonitis (RP) remains a major complication in non-small cell lung cancer (NSCLC) patients undergoing radiochemotherapy (RCHT). Traditionally, the mean lung dose (MLD) and the volume of the total lung receiving at least 20 Gy (V20Gy) are used to predict RP in patients treated with normo-fractionated photon therapy. However, other models, including the actual dose-distribution in the lungs using the effective α/ß model or a combination of radiation doses to the lungs and heart, have been proposed for predicting RP. Moreover, the models established for photons may not hold for patients treated with passively-scattered proton therapy (PSPT). Therefore, we here tested and validated novel predictive parameters for RP in NSCLC patient treated with PSPT. METHODS: Data on the occurrence of RP, structure files and dose-volume histogram parameters for lungs and heart of 96 NSCLC patients, treated with PSPT and concurrent chemotherapy, was retrospectively retrieved from prospective clinical studies of two international centers. Data was randomly split into a training set (64 patients) and a validation set (32 patients). Statistical analyses were performed using binomial logistic regression. RESULTS: The biologically effective dose (BED) of the'lungs - GTV' significantly predicted RP ≥ grade 2 in the training-set using both a univariate model (p = 0.019, AUCtrain = 0.72) and a multivariate model in combination with the effective α/ß parameter of the heart (pBED = 0.006, [Formula: see text] = 0.043, AUCtrain = 0.74). However, these results did not hold in the validation-set (AUCval = 0.52 andAUCval = 0.50, respectively). Moreover, these models were found to neither outperform a model built with the MLD (p = 0.015, AUCtrain = 0.73, AUCval = 0.51), nor a multivariate model additionally including the V20Gy of the heart (pMLD = 0.039, pV20Gy,heart = 0.58, AUCtrain = 0.74, AUCval = 0.53). CONCLUSION: Using the effective α/ß parameter of the lungs and heart we achieved similar performance to commonly used models built for photon therapy, such as MLD, in predicting RP ≥ grade 2. Therefore, prediction models developed for photon RCHT still hold for patients treated with PSPT.

MRI magnitude signal-based proton beam visualisation in water phantoms reflects composite effects of beam-induced buoyant convection and radiation chemistry.

Objective. Local magnetic resonance (MR) signal loss was previously observed during proton beam irradiation of free-floating water phantoms at ambient temperature using a research prototype in-beam magnetic resonance imaging (MRI) scanner. The emergence of this MR signal loss was hypothesised to be dependent on beam-induced convection. The aim of this study was therefore to unravel whether physical conditions allowing the development of convection must prevail for the beam-induced MRI signatures to emerge.Approach. The convection dependence of MRI magnitude signal-based proton beam visualisation was investigated in combined irradiation and imaging experiments using a gradient echo (GE)-based time-of-flight (ToF) angiography pulse sequence, which was first tested for its suitability for proton beam visualisation in free-floating water phantoms at ambient temperature. Subsequently, buoyant convection was selectively suppressed in water phantoms using either mechanical barriers or temperature control of water expansivity. The underlying contrast mechanism was further assessed using sagittal imaging and variation of T1 relaxation time-weighting.Main results. In the absence of convection-driven water flow, weak beam-induced MR signal changes occurred, whereas strong changes did occur when convection was not mechanically or thermally inhibited. Moreover, the degree of signal loss was found to change with the variation of T1-weighting. Consequently, beam-induced MR signal loss in free-floating water phantoms at ambient temperature does not exclusively originate from buoyant convection, but is caused by local composite effects of beam-induced motion and radiation chemistry resulting in a local change in the water T1 relaxation time.Significance. The identification of ToF angiography sequence-based proton beam visualisation in water phantoms to result from composite effects of beam-induced motion and radiation chemistry represents the starting point for the future elucidation of the currently unexplained motion-based MRI contrast mechanism and the identification of the proton beam-induced material change causing T1 relaxation time lengthening.

Investigation of contrast mechanisms for MRI phase signal-based proton beam visualization in water phantoms.

PURPOSE: The low sensitivity and limitation to water phantoms of convection-dependent MRI magnitude signal-based proton beam visualization hinder its in vivo applicability in MR-integrated proton beam therapy. The purpose of the present study was, therefore, to assess possible contrast mechanisms for MRI phase signal-based proton beam visualization that can potentially be exploited to enhance the sensitivity of the method and extend its applicability to tissue materials. METHODS: To assess whether proton beam-induced magnetic field perturbations, changes in material susceptibility or convection result in detectable changes in the MRI phase signal, water phantom characteristics, experiment timing, and imaging parameters were varied in combined irradiation and imaging experiments using a time-of-flight angiography pulse sequence on a prototype in-beam MRI scanner. Velocity encoding was used to further probe and quantify beam-induced convection. RESULTS: MRI phase signal-based proton beam visualization proved feasible. The observed phase difference contrast was evoked by beam-induced buoyant convection with flow velocities in the mm/s range. Proton beam-induced magnetic field perturbations or changes in magnetic susceptibility did not influence the MRI phase signal. Velocity encoding was identified as a means to enhance the detection sensitivity. CONCLUSION: Because the MRI phase difference contrast observed during proton beam irradiation of water phantoms is caused by beam-induced convection, this method will unlikely be transferable to tightly compartmentalized tissue wherein flow effects are restricted. However, strong velocity encoded pulse sequences were identified as promising candidates for the future development of MRI-based methods for water phantom-based geometric quality assurance in MR-integrated proton beam therapy.

Direct visualization of proton beam irradiation effects in liquids by MRI.

The main advantage proton beams offer over photon beams in radiation therapy of cancer patients is the dose maximum at their finite range, yielding a reduction in the dose deposited in healthy tissues surrounding the tumor. Since no direct method exists to measure the beam's range during dose delivery, safety margins around the tumor are applied, compromising the dose conformality and reducing the targeting accuracy. Here, we demonstrate that online MRI can visualize the proton beam and reveal its range during irradiation of liquid-filled phantoms. A clear dependence on beam energy and current was found. These results stimulate research into novel MRI-detectable beam signatures and already find application in the geometric quality assurance for magnetic resonance-integrated proton therapy systems currently under development.

Technical note: Experimental dosimetric characterization of proton pencil beam distortion in a perpendicular magnetic field of an in-beam MR scanner.

BACKGROUND: As it promises more precise and conformal radiation treatments, magnetic resonance imaging-integrated proton therapy (MRiPT) is seen as a next step in image guidance for proton therapy. The Lorentz force, which affects the course of the proton pencil beams, presents a problem for beam delivery in the presence of a magnetic field. PURPOSE: To investigate the influence of the 0.32-T perpendicular magnetic field of an MR scanner on the delivery of proton pencil beams inside an MRiPT prototype system. METHODS: An MRiPT prototype comprising of a horizontal pencil beam scanning beam line and an open 0.32-T MR scanner was used to evaluate the impact of the vertical magnetic field on proton beam deflection and dose spot pattern deformation. Three different proton energies (100, 150, and 220 MeV) and two spot map sizes (15 × 15 and 30 × 20 cm2 ) at four locations along the beam path without and with magnetic field were measured. Pencil-beam dose spots were measured using EBT3 films and a 2D scintillation detector. To study the magnetic field effects, a 2D Gaussian fit was applied to each individual dose spot to determine the central position ( X , Y ) $(X,Y)$ , minimum and maximum lateral standard deviation ( σ m i n $\sigma _{min}$ and σ m a x $\sigma _{max}$ ), orientation (θ), and the eccentricity (ε). RESULTS: The dose spots were subjected to three simultaneous effects: (a) lateral horizontal beam deflection, (b) asymmetric trapezoidal deformation of the dose spot pattern, and (c) deformation and rotation of individual dose spots. The strongest effects were observed at a proton energy of 100 MeV with a horizontal beam deflection of 14-186 mm along the beam path. Within the central imaging field of the MR scanner, the maximum relative dose spot size σ m a x $\sigma _{max}$ decreased by up to 3.66%, while σ m i n $\sigma _{min}$ increased by a maximum of 2.15%. The largest decrease and increase in the eccentricity of the dose spots were 0.08 and 0.02, respectively. The spot orientation θ was rotated by a maximum of 5.39°. At the higher proton energies, the same effects were still seen, although to a lesser degree. CONCLUSIONS: The effect of an MRiPT prototype's magnetic field on the proton beam path, dose spot pattern, and dose spot form has been measured for the first time. The findings show that the impact of the MF must be appropriately recognized in a future MRiPT treatment planning system. The results emphasize the need for additional research (e.g., effect of magnetic field on proton beams with range shifters and impact of MR imaging sequences) before MRiPT applications can be employed to treat patients.

Reduction of intrafraction pancreas motion using an abdominal corset compatible with proton therapy and MRI.

Background and purpose: Motion mitigation is of crucial importance in particle therapy (PT) of patients with abdominal tumors to ensure high-precision irradiation. Magnetic resonance imaging (MRI) is an excellent modality for target volume delineation and motion estimation of mobile soft-tissue tumors. Thus, the aims of this study were to develop an MRI- and PT-compatible abdominal compression device, to investigate its effect on pancreas motion reduction, and to evaluate patient tolerability and acceptance. Materials and methods: In a prospective clinical study, 16 patients with abdominal tumors received an individualized polyethylene-based abdominal corset. Pancreas motion was analyzed using time- and phase resolved MRI scans (orthogonal 2D-cine and 4D MRI) with and without compression by the corset. The pancreas was manually segmented in each MRI data set and the population-averaged center-of-mass motion in inferior-superior (IS), anterior-posterior (AP) and left-right (LR) directions was determined. A questionnaire was developed to investigate the level of patient acceptance of the corset, which the patients completed after acquisition of the planning computed tomography (CT) and MRI scans. Results: The corset was found to reduce pancreas motion predominantly in IS direction by on average 47 % - 51 % as found in the 2D-cine and 4D MRI data, respectively, while motion in the AP and LR direction was not significantly reduced. Most patients reported no discomfort when wearing the corset. Conclusion: An MRI- and PT-compatible individualized abdominal corset was presented, which substantially reduced breathing-induced pancreas motion and can be safely applied with no additional discomfort for the patients. The corset has been successfully integrated into our in-house clinical workflow for PT of tumors of the upper abdomen.

Technical Note: ADAM PETer - An anthropomorphic, deformable and multimodality pelvis phantom with positron emission tomography extension for radiotherapy.

OBJECTIVE: To develop an anthropomorphic, deformable and multimodal pelvis phantom with positron emission tomography extension for radiotherapy (ADAM PETer). METHODS: The design of ADAM PETer was based on our previous pelvis phantom (ADAM) and extended for compatibility with PET and use in 3T magnetic resonance imaging (MRI). The formerly manually manufactured silicon organ surrogates were replaced by three-dimensional (3D) printed organ shells. Two intraprostatic lesions, four iliac lymph node metastases and two pelvic bone metastases were added to simulate prostate cancer as multifocal and metastatic disease. Radiological properties [computed tomography (CT) and 3T MRI] of cortical bone, bone marrow and adipose tissue were simulated by heavy gypsum, a mixture of Vaseline and K2 HPO4 and peanut oil, respectively. For soft tissues, agarose gels with varying concentrations of agarose, gadolinium (Gd) and sodium fluoride (NaF) were developed. The agarose gels were doped with patient-specific activity concentrations of a Fluorine-18 labelled compound and then filled into the 3D printed organ shells of prostate lesions, lymph node and bone metastases. The phantom was imaged at a dual energy CT and a 3T PET/MRI scanner. RESULTS: The compositions of the soft tissue surrogates are the following (given as mass fractions of agarose[w%]/NaF[w%]/Gd[w%]): Muscle (4/1/0.027), prostate (1.35/4.2/0.011), prostate lesions (2.25/4.2/0.0085), lymph node and bone metastases (1.4/4.2/0.025). In all imaging modalities, the phantom simulates human contrast. Intraprostatic lesions appear hypointense as compared to the surrounding normal prostate tissue in T2-weighted MRI. The PET signal of all tumors can be localized as focal spots at their respective site. Activity concentrations of 12.0 kBq/mL (prostate lesion), 12.4 kBq/mL (lymph nodes) and 39.5 kBq/mL (bone metastases) were measured. CONCLUSION: The ADAM PETer pelvis phantom can be used as multimodal, anthropomorphic model for CT, 3T-MRI and PET measurements. It will be central to simulate and optimize the technical workflow for the integration of PET/MRI-based radiation treatment planning of prostate cancer patients.

Characterization of magnetic interference and image artefacts during simultaneous in-beam MR imaging and proton pencil beam scanning.

For the first time, a low-field open magnetic resonance (MR) scanner was combined with a proton pencil beam scanning (PBS) research beamline. The aim of this study was to characterize the magnetic fringe fields produced by the PBS system and measure their effects on MR image quality during simultaneous PBS irradiation and image acquisition. A magnetic field camera measured the change in central resonance frequency (Δf res) and magnetic field homogeneity (ΔMFH) of the B0 field of the MR scanner during operation of the beam transport and scanning magnets. The beam energy was varied between 70 - 220 MeV and beam scanning was performed along the central horizontal and vertical axis of a 48 × 24 cm2 radiation field. The time structure of the scanning magnets' fringe fields was simultaneously recorded by a tri-axial Hall probe. MR imaging experiments were conducted using the ACR (American College of Radiology) Small MRI Phantom and a spoiled gradient echo pulse sequence during simultaneous volumetric irradiation. Computer simulations were performed to predict the effects of B 0 field perturbations due to PBS irradiation on MR image formation in k-space. Setting the beam transport magnets, horizontal and vertical scanning magnets resulted in a maximum Δf res of 50, 235 and 4 Hz, respectively. The ΔMFH was less than 3 parts per million for all measurements. MR images acquired during beam energy variation and vertical beam scanning showed no visual loss in image quality. However, MR images acquired during horizontal beam scanning showed severe coherent ghosting artefacts in phase encoding direction. Both simulated and measured k-space phase maps prove that these artefacts are caused by phase-offsets. This study shows first experimental evidence that simultaneous in-beam MR imaging during proton PBS irradiation is subject to severe loss of image quality in the absence of magnetic decoupling between the PBS and MR system.

CT-based attenuation correction of whole-body radiotherapy treatment positioning devices in PET/MRI hybrid imaging.

OBJECTIVE: To implement computed tomography (CT)-based attenuation maps of radiotherapy (RT) positioning hardware and radiofrequency (RF) coils to enable hybrid positron emission tomography/magnetic resonance imaging (PET/MRI)-based RT treatment planning. MATERIALS AND METHODS: The RT positioning hardware consisted of a flat RT table overlay, coil holders for abdominal scans, coil holders for head and neck scans and an MRI compatible hip and leg immobilization device. CT images of each hardware element were acquired on a CT scanner. Based on the CT images, attenuation maps of the devices were created. Validation measurements were performed on a PET/MR scanner using a 68Ge phantom (48 MBq, 10 min scan time). Scans with each device in treatment position were performed. Then, reference scans containing only the phantom were taken. The scans were reconstructed online (at the PET/MRI scanner) and offline (via e7tools on a PC) using identical reconstruction parameters. Average reconstructed activity concentrations of the device and reference scans were compared. RESULTS: The device attenuation maps were successfully implemented. The RT positioning devices caused an average decrease of reconstructed PET activity concentration in the range between -8.3 ± 2.1% (mean ± SD) (head and neck coil holder with coils) to -1.0 ± 0.5% (abdominal coil holder). With attenuation correction taking into account RT hardware, these values were reduced to -2.0 ± 1.2% and -0.6 ± 0.5%, respectively. The results of the offline and online reconstructions were nearly identical, with a difference of up to 0.2%. CONCLUSION: The decrease in reconstructed activity concentration caused by the RT positioning devices is clinically relevant and can successfully be corrected using CT-based attenuation maps. Both the offline and online reconstruction methods are viable options.

Advanced Sandwich Composite Cores for Patient Support in Advanced Clinical Imaging and Oncology Treatment.

Ongoing advances in both imaging and treatment for oncology purposes have seen a significant rise in the use of not only the individual imaging modalities, but also their combination in single systems such as Positron Emission Tomography combined with Computed Tomography (PET-CT) and PET-MRI (Magnetic Resonance Imaging) when planning for advanced oncology treatment, the most demanding of which is proton therapy. This has identified issues in the availability of suitable materials upon which to support the patient undergoing imaging and treatment owing to the differing requirements for each of the techniques. Sandwich composites are often selected to solve this issue but there is little information regarding optimum materials for their cores. In this paper, we presented a range of materials which are suitable for such purposes and evaluated the performance for use in terms of PET signal attenuation, proton beam stopping, MRI signal shading and X-Ray CT visibility. We found that Extruded Polystyrene offers the best compromise for patient support and positioning structures across all modalities tested, allowing for significant savings in treatment planning time and delivering more efficient treatment with lower margins.

MR-guided proton therapy: a review and a preview.

BACKGROUND: The targeting accuracy of proton therapy (PT) for moving soft-tissue tumours is expected to greatly improve by real-time magnetic resonance imaging (MRI) guidance. The integration of MRI and PT at the treatment isocenter would offer the opportunity of combining the unparalleled soft-tissue contrast and real-time imaging capabilities of MRI with the most conformal dose distribution and best dose steering capability provided by modern PT. However, hybrid systems for MR-integrated PT (MRiPT) have not been realized so far due to a number of hitherto open technological challenges. In recent years, various research groups have started addressing these challenges and exploring the technical feasibility and clinical potential of MRiPT. The aim of this contribution is to review the different aspects of MRiPT, to report on the status quo and to identify important future research topics. METHODS: Four aspects currently under study and their future directions are discussed: modelling and experimental investigations of electromagnetic interactions between the MRI and PT systems, integration of MRiPT workflows in clinical facilities, proton dose calculation algorithms in magnetic fields, and MRI-only based proton treatment planning approaches. CONCLUSIONS: Although MRiPT is still in its infancy, significant progress on all four aspects has been made, showing promising results that justify further efforts for research and development to be undertaken. First non-clinical research solutions have recently been realized and are being thoroughly characterized. The prospect that first prototype MRiPT systems for clinical use will likely exist within the next 5 to 10 years seems realistic, but requires significant work to be performed by collaborative efforts of research groups and industrial partners.

First application of a high-resolution silicon detector for proton beam Bragg peak detection in a 0.95 T magnetic field.

PURPOSE: To report on experimental results of a high spatial resolution silicon-based detector exposed to therapeutic quality proton beams in a 0.95 T transverse magnetic field. These experimental results are important for the development of accurate and novel dosimetry methods in future potential real-time MRI-guided proton therapy systems. METHODS: A permanent magnet device was utilized to generate a 0.95 T magnetic field over a 4 × 20 × 15 cm3 volume. Within this volume, a high-resolution silicon diode array detector was positioned inside a PMMA phantom of 4 × 15 × 12 cm3 . This detector contains two orthogonal strips containing 505 sensitive volumes spaced at 0.2 mm apart. Proton beams collimated to a circle of 10 mm diameter with nominal energies of 90 MeV, 110 MeV, and 125 MeV were incident on the detector from an edge-on orientation. This allows for a measurement of the Bragg peak at 0.2 mm spatial resolution in both the depth and lateral profile directions. The impact of the magnetic field on the proton beams, that is, a small deflection was also investigated. A Geant4 Monte Carlo simulation was performed of the experimental setup to aid in interpretation of the results. RESULTS: The nominal Bragg peak for each proton energy was successfully observed with a 0.2 mm spatial resolution in the 0.95 T transverse magnetic field in both a depth and lateral profiles. The proton beam deflection (at 0.95 T) was a consistent 2 ±0.5 mm at the center of the magnetic volume for each beam energy. However, a pristine Bragg peak was not observed for each energy. This was caused by the detector packaging having small air gaps between layers of the phantom material surrounding the diode array. These air gaps act to degrade the shape of the Bragg peak, and further to this, the nonwater equivalent silicon chip acts to separate the Bragg peak into multiple peaks depending on the proton path taken. Overall, a promising performance of the silicon detector array was observed, however, with a qualitative assessment rather than a robust quantitative dosimetric evaluation at this stage of development. CONCLUSIONS: For the first time, a high-resolution silicon-based radiation detector has been used to measure proton beam Bragg peak deflections in a phantom due to a strong magnetic field. Future efforts are required to optimize the detector packaging to strengthen the robustness of the dosimetric quantities obtained from the detector. Such high-resolution silicon diode arrays may be useful in future efforts in MRI-guided proton therapy research.

Performing clinical 18F-FDG-PET/MRI of the mediastinum optimising a dedicated, patient-friendly protocol.

OBJECTIVE: To construct a mediastinal-specific fluorine-18-fluorodeoxyglucose (F-FDG)-PET/MR protocol with high-quality MRI of minimal acquisition-time and comparable diagnostic value to F-FDG-PET/computed tomography (CT). MATERIALS AND METHODS: Fifteen healthy participants received PET/MRI and 10 patients with mediastinal tumours (eight non-small-cell lung, two oesophageal cancer) received F-FDG-PET/MRI immediately after F-FDG-PET/CT. Sequences volume interpolated breath-hold examination (T1-VIBE) and Half-Fourier acquisition single-shot turbo spin echo (T2-HASTE) were optimised by varying the parameters: breath-hold (BH, end-expiration), fat suppression (spectral adiabatic inversion recovery), and ECG-triggering (ECG, end-diastole). Image quality (IQ) of each sequence-variation was qualitatively scored by medical experts and quantitatively assessed by calculating signal-to-noise ratios, contrast relative to muscle, standardized-uptake-value, and tumour-to-blood ratios. Patient comfort was evaluated on patients' experience. Diagnostic accuracy of F-FDG-PET/MRI was compared to F-FDG-PET/CT, in reference to histopathology/cytopathology. RESULTS: ECG-triggered T1-VIBE images showed the highest signal-to-noise ratio (P < 0.01) and the largest contrast between mediastinal soft-tissues, regardless of BH or free-breathing acquisition. IQ of ECG-triggered T1-VIBE scans in BH were scored qualitatively highest with good reader agreement (κ = 0.62). IQ of T2-HASTE was not significantly affected by BH acquisition (P > 0.9). Qualitative IQ of T1-VIBE and T2-HASTE declined after spectral adiabatic inversion recovery fat-suppression. All patients could maintain BH at end-expiration and reported no discomfort. Diagnostic performance of F-FDG-PET/MR was not significantly different from F-FDG-PET/CT with comparable staging, standardized-uptake-values, and tumour-to-blood ratios. However, T-status was more often over-staged on F-FDG-PET/CT, while N-status was more frequently under-staged on F-FDG-PET/MR. CONCLUSION: ECG-triggered T1-VIBE sequences acquired during short, multiple BHs are recommended for mediastinal imaging using F-FDG-PET/MR. With dedicated protocols, F-FDG-PET/MRI will be useful in thoracic oncology and aid in diagnostic evaluation and tailored treatment decision-making.

Comparison of pancreatic respiratory motion management with three abdominal corsets for particle radiation therapy: Case study.

BACKGROUND AND PURPOSE: Abdominal organ motion seriously compromises the targeting accuracy for particle therapy in patients with pancreatic adenocarcinoma. This study compares three different abdominal corsets regarding their ability to reduce pancreatic motion and their potential usability in particle therapy. MATERIALS AND METHODS: A patient-individualized polyurethane (PU), a semi-individualized polyethylene (PE), and a patient-individualized three-dimensional-scan based polyethylene (3D-PE) corset were manufactured for one healthy volunteer. Time-resolved volumetric four-dimensional-magnetic resonance imaging (4D-MRI) and single-slice two-dimensional (2D) cine-MRI scans were acquired on two consecutive days to compare free-breathing motion patterns with and without corsets. The corset material properties, such as thickness variance, material homogeneity in Hounsfield units (HU) on computed tomography (CT) scans, and manufacturing features were compared. The water equivalent ratio (WER) of corset material samples was measured using a multi-layer ionization chamber for proton energies of 150 and 200 MeV. RESULTS: All corsets reduced the pancreatic motion on average by 9.6 mm in inferior-superior and by 3.2 mm in anterior-posterior direction. With corset, the breathing frequency was approximately doubled and the day-to-day motion variations were reduced. The WER measurements showed an average value of 0.993 and 0.956 for the PE and 3DPE corset, respectively, and of 0.298 for the PU corset. The PE and 3DPE corsets showed a constant thickness of 2.8 ± 0.2 and 3.8 ± 0.2 mm, respectively and a homogeneous material composition with a standard deviation (SD) of 31 and 32 HU, respectively. The PU corset showed a variable thickness of 4.2 - 25.6 mm and a heterogeneous structure with air inclusions with an SD of 113 HU. CONCLUSION: Abdominal corsets may be effective devices to reduce pancreatic motion. For particle therapy, PE-based corsets are preferred over PU-based corset due to their material homogeneity and constant thickness.

Detectability and structural stability of a liquid fiducial marker in fresh ex vivo pancreas tumour resection specimens on CT and 3T MRI.

PURPOSE: To test the detectability of a liquid fiducial marker injected into ex vivo pancreas tumour tissue on magnetic resonance imaging (MRI) and computed tomography (CT). Furthermore, its injection performance using different needle sizes and its structural stability after fixation in formaldehyde were investigated. METHODS: Liquid fiducial markers with a volume of 20-100 µl were injected into freshly resected pancreas specimens of three patients with suspected adenocarcinoma. X­ray guided injection was performed using different needle sizes (18 G, 22 G, 25 G). The specimens were scanned on MRI and CT with clinical protocols. The markers were segmented on CT by signal thresholding. Marker detectability in MRI was assessed in the registered segmentations. Marker volume on CT was compared to the injected volume as a measure of backflow. RESULTS: Markers with a volume ≥20 µl were detected as hyperintensity on X­ray and CT. On T1- and T2-weighted 3T MRI, marker sizes ranging from 20-100 µl were visible as hypointensity. Since most markers were non-spherical, MRI detectability was poor and their differentiation from hypointensities caused by air cavities or surgical clips was only feasible with a reference CT. Marker backflow was only observed when using an 18-G needle. A volume decrease of 6.6 ± 13.0% was observed after 24 h in formaldehyde and, with the exception of one instance, no wash-out occurred. CONCLUSIONS: The liquid fiducial marker injected in ex vivo pancreatic resection specimen was visible as hyperintensity on kV X­ray and CT and as hypointensity on MRI. The marker's size was stable in formaldehyde. A marker volume of ≥50 µL is recommended in clinically used MRI sequences. In vivo injection is expected to improve the markers sphericity due to persisting metabolism and thereby enhance detectability on MRI.

Commissioning of a 4D MRI phantom for use in MR-guided radiotherapy.

PURPOSE: Systems for integrated magnetic resonance guided radiation therapy (MRgRT) provide real-time and online MRI guidance for unequaled targeting performance of moving tumors and organs at risk. The clinical introduction of such systems requires dedicated methods for commissioning and routine machine quality assurance (QA). The aim of the study was to develop a commissioning protocol and method for automatic quantification of target motion and geometric accuracy using a 4D MRI motion phantom. MATERIALS AND METHODS: The commissioning was performed on a clinically used 3 T MR scanner. The phantom was positioned on a flat tabletop overlay using an in-house constructed base plate for a quick and reproducible setup. The torso-shaped phantom body, which was filled with mineral oil as signal generating medium, included a 3D grid structure for image distortion analysis and a cylindrical thru-hole in which a software-controlled moving rod with a hypo-intense background gel and a decentralized hyper-intense target simulated 3D organ motion patterns. To allow for sequence optimization, MR relaxometry was performed to determine the longitudinal T1 and transverse T2 relaxation times of both target and background gel in the movable cylinder. The geometric image distortion was determined as the mean and maximum 3D Euclidean distance (Δmean , Δmax ) of grid points determined by nonrigid registration of a 3D spoiled gradient echo MRI scan and a CT scan. Sinusoidal 1D/2D/3D motion trajectories, varying in amplitude and frequency, as well as an exemplary 1D MR navigator diaphragm motion pattern extracted from a healthy volunteer scan, were scanned by means of 2D cine MRI and 4D MRI. Target positions were automatically extracted from 2D cine MRI using an in-house developed software tool. RESULTS: The base plate enabled a reproducible setup with a deviation of <1 mm in all directions. Relaxometry yielded T1 /T2 values for target and background gel of 208.1 ± 2.8/30.5 ± 4.7 ms and 871 ± 36/13.4  ±  1.3 ms, respectively. The 3D geometric image distortion increased with distance from the magnetic isocenter, with Δmean  = 0.58 ± 0.30 mm and Δmax  = 1.31 mm. The frequencies of the reconstructed motion patterns agreed with the preset values within 0.5%, whereas the reconstructed amplitudes showed a maximum deviation to the preset amplitudes of <0.5 mm in AP/LR direction and <0.3 mm in IS direction. CONCLUSION: A method and protocol for commissioning of a 4D MRI motion phantom on a 3 T MR scanner for MRgRT was developed. High-contrast and geometrically reliable 2D cine MR images of the phantom's moving target could be obtained. The preset motion parameters could be extracted with sufficient spatio-temporal accuracy from 2D cine MRI in all motion directions. The overall 3D geometric image distortion of <1.31 mm within the phantom grid confirms geometric accuracy of the clinically utilized 3D spoiled gradient echo sequence. The method developed can be used for routine QA tests of spatio-temporally resolved MRI data in MRgRT.

Characterizing geometrical accuracy in clinically optimised 7T and 3T magnetic resonance images for high-precision radiation treatment of brain tumours.

BACKGROUND AND PURPOSE: In neuro-oncology, high spatial accuracy is needed for clinically acceptable high-precision radiation treatment planning (RTP). In this study, the clinical applicability of anatomically optimised 7-Tesla (7T) MR images for reliable RTP is assessed with respect to standard clinical imaging modalities. MATERIALS AND METHODS: System- and phantom-related geometrical distortion (GD) were quantified on clinically-relevant MR sequences at 7T and 3T, and on CT images using a dedicated anthropomorphic head phantom incorporating a 3D grid-structure, creating 436 points-of-interest. Global GD was assessed by mean absolute deviation (MADGlobal). Local GD relative to the magnetic isocentre was assessed by MADLocal. Using 3D displacement vectors of individual points-of-interest, GD maps were created. For clinically acceptable radiotherapy, 7T images need to meet the criteria for accurate dose delivery (GD < 1 mm) and present comparable GD as tolerated in clinically standard 3T MR/CT-based RTP. RESULTS: MADGlobal in 7T and 3T images ranged from 0.3 to 2.2 mm and 0.2-0.8 mm, respectively. MADLocal increased with increasing distance from the isocentre, showed an anisotropic distribution, and was significantly larger in 7T MR sequences (MADLocal = 0.2-1.2 mm) than in 3T (MADLocal = 0.1-0.7 mm) (p < 0.05). Significant differences in GD were detected between 7T images (p < 0.001). However, maximum MADLocal remained ≤1 mm within 68.7 mm diameter spherical volume. No significant differences in GD were found between 7T and 3T protocols near the isocentre. CONCLUSIONS: System- and phantom-related GD remained ≤1 mm in central brain regions, suggesting that 7T MR images could be implemented in radiotherapy with clinically acceptable spatial accuracy and equally tolerated GD as in 3T MR/CT-based RTP. For peripheral regions, GD should be incorporated in safety margins for treatment uncertainties. Moreover, the effects of sequence-related factors on GD needs further investigation to obtain RTP-specific MR protocols.

Ano-rectal wall dose-surface maps localize the dosimetric benefit of hydrogel rectum spacers in prostate cancer radiotherapy.

BACKGROUND AND PURPOSE: To evaluate spatial differences in dose distributions of the ano-rectal wall (ARW) using dose-surface maps (DSMs) between prostate cancer patients receiving intensity-modulated radiation therapy with and without implantable rectum spacer (IMRT+IRS; IMRT-IRS, respectively), and to correlate this with late gastro-intestinal (GI) toxicities using validated spatial and non-spatial normal-tissue complication probability (NTCP) models. MATERIALS AND METHODS: For 26 patients DSMs of the ARW were generated. From the DSMs various shape-based dose measures were calculated at different dose levels: lateral extent, longitudinal extent, and eccentricity. The contiguity of the ARW dose distribution was assessed by the contiguous-DSH (cDSH). Predicted complication rates between IMRT+IRS and IMRT-IRS plans were assessed using a spatial NTCP model and compared against a non-spatial NTCP model. RESULTS: Dose surface maps are generated for prostate radiotherapy using an IRS. Lateral extent, longitudinal extent and cDSH were significantly lower in IMRT+IRS than for IMRT-IRS at high-dose levels. Largest significant differences were observed for cDSH at dose levels >50 Gy, followed by lateral extent at doses >57 Gy, and longitudinal extent in anterior and superior-inferior directions. Significant decreases (p = 0.01) in median rectal and anal NTCPs (respectively, Gr 2 late rectal bleeding and subjective sphincter control) were predicted when using an IRS. CONCLUSIONS: Local-dose effects are predicted to be significantly reduced by an IRS. The spatial NTCP model predicts a significant decrease in Gr 2 late rectal bleeding and subjective sphincter control. Dose constraints can be improved for current clinical treatment planning.