Three-row MRI receive array with remote circuitry to preserve radiation transparency

Objective.Up to this point, 1.5 T linac-compatible coil array layouts have been restricted to one or two rows of coils because of the desire to place radiation-opaque circuitry adjacent to the coils and outside the window through which the linac beam travels. Such layouts can limit parallel imaging performance. The purpose of this work was to design and build a three-row array in which remotely located circuits permitted a central row of coils while preserving the radiolucent window.Approach.The remote circuits consisted of a phase shifter to cancel the phase introduced by the coaxial link between the circuit and coil, followed by standard components for tuning, matching, detuning, and preamplifier decoupling. Tests were performed to compare prototype single-channel coils with remote or local circuits, which were followed by tests comparing two and three-row arrays .Main results.The single-channel coil with the remote circuit maintained 85% SNR at depths of 30 mm or more as compared to a coil with local circuit. The three-row array provided similar SNR as the two-row array, along with geometry factor advantages for parallel imaging acceleration in the head-foot direction.Significance.The remote circuit strategy could potentially support future MR-linac arrays by allowing greater flexibility in array layout compared to those confined by local circuits, which can be leveraged for parallel imaging acceleration.

MR compatible detectors assessment for a 0.35 T MR-linac commissioning

PURPOSE

To assess a large panel of MR compatible detectors on the full range of measurements required for a 0.35 T MR-linac commissioning by using a specific statistical method represented as a continuum of comparison with the Monte Carlo (MC) TPS calculations. This study also describes the commissioning tests and the secondary MC dose calculation validation.

MATERIAL AND METHODS

Plans were created on the Viewray TPS to generate MC reference data. Absolute dose points, PDD, profiles and output factors were extracted and compared to measurements performed with ten different detectors: PTW 31010, 31021, 31022, Markus 34045 and Exradin A28 MR ionization chambers, SN Edge shielded diode, PTW 60019 microdiamond, PTW 60023 unshielded diode, EBT3 radiochromic films and LiF µcubes. Three commissioning steps consisted in comparison between calculated and measured dose: the beam model validation, the output calibration verification in four different phantoms and the commissioning tests recommended by the IAEA-TECDOC-1583.

MAIN RESULTS

The symmetry for the high resolution detectors was higher than the TPS data of about 1%. The angular responses of the PTW 60023 and the SN Edge were - 6.6 and - 11.9% compared to the PTW 31010 at 60°. The X/Y-left and the Y-right penumbras measured by the high resolution detectors were in good agreement with the TPS values except for the PTW 60023 for large field sizes. For the 0.84 × 0.83 cm2 field size, the mean deviation to the TPS of the uncorrected OF was - 1.7 ± 1.6% against - 4.0 ± 0.6% for the corrected OF whereas we found - 4.8 ± 0.8% for passive dosimeters. The mean absolute dose deviations to the TPS in different phantoms were 0 ± 0.4%, - 1.2 ± 0.6% and 0.5 ± 1.1% for the PTW 31010, PTW 31021 and Exradin A28 MR respectively.

CONCLUSIONS

The magnetic field effects on the measurements are considerably reduced at low magnetic field. The PTW 31010 ionization chamber can be used with confidence in different phantoms for commissioning and QA tests requiring absolute dose verifications. For relative measurements, the PTW 60019 presented the best agreement for the full range of field size. For the profile assessment, shielded diodes had a behaviour similar to the PTW 60019 and 60023 while the ionization chambers were the most suitable detectors for the symmetry. The output correction factors published by the IAEA TRS 483 seem to be applicable at low magnetic field pending the publication of new MR specific values.

Multileaf collimator characterization and modeling for a 1.5 T MR-linac using static synchronous and asynchronous sweeping gaps

Objective.The Elekta unity MR-linac delivers step-and-shoot intensity modulated radiotherapy plans using a multileaf collimator (MLC) based on the Agility MLC used on conventional Elekta linacs. Currently, details of the physical Unity MLC and the computational model within its treatment planning system (TPS)Monacoare lacking in published literature. Recently, a novel approach to characterize the physical properties of MLCs was introduced using dynamic synchronous and asynchronous sweeping gap (aSG) tests. Our objective was to develop a step-and-shoot version of the dynamic aSG test to characterize the Unity MLC and the computational MLC models in theMonacoandRayStationTPSs.Approach.Dynamic aSG were discretized into a step-and-shoot aSG by investigating the number of segments/sweep and the minimal number of monitor units (MU) per segment. The step-and-shoot aSG tests were compared to the dynamic aSG tests on a conventional linac at a source-to-detector distance of 143.5 cm, mimicking the Unity configuration. the step-and-shoot aSG tests were used to characterize the Unity MLC through measurements and dose calculations in both TPSs.Main results.The step-and-shoot aSGs tests with 100 segments and 5 MU/segment gave results very similar to the dynamic aSG experiments. The effective tongue-and-groove width of the Unity gradually increased up to 1.4 cm from the leaf tip end. The MLC models inRayStationandMonacoagreed with experimental data within 2.0% and 10%, respectively. The largest discrepancies inMonacowere found for aSG tests with >10 mm leaf interdigitation, which are non-typical for clinical plans.Significance.The step-and-shoot aSG tests accurately characterize the MLC in step-and-shoot delivery mode. The MLC model inRayStation2023B accurately describes the tongue-and-groove and leaf tip effects whereasMonacooverestimates the tongue-and-groove shadowing further away from the leaf tip end.

MR-linac MLC positioning QA by digitally stitching dual double-exposed films

PURPOSE

The picket fence (PF) test is highly recommended for multi-leaf collimator (MLC) quality assurance. However, since the electronic portal imaging device (EPID) on the Elekta Unity only covers a small area, it is not feasible to perform the PF test for the entire MLC. Here, we propose a technique for the PF test by stitching two double-exposed films.

METHODS

Two EBT3 films were used to encompass the entire MLC, with each one covering one half of the area. Two fields were employed to apply double exposure: a PF pattern consisting of 11 2 mm wide pickets and a 2.84 cm x 22 cm open field. The edges of the open field defined by the diaphragms were used to correct film rotation as well as align them horizontally. The PF pattern was also measured with the EPID where the pickets were used to align the films vertically. Individual leaf positions were detected on the merged film for quantitative analysis. Various MLC positioning errors were introduced to evaluate the technique's sensitivity.

RESULTS

The merged films covered 72 leaf pairs properly (four leaf pairs on both sides were outside the treatment couch). With the EPID, the leaf positioning accuracy was -0.02 ± 0.07 mm (maximum: 0.29 mm) and the picket width variation was 0.00 ± 0.03 mm (maximum: 0.11 mm); with the films, the position accuracy and width variation were -0.03 ± 0.13 mm (maximum: 0.80 mm) and 0.00 ± 0.13 mm (maximum: 0.74 mm), respectively. The EPID was able to detect errors of 0.5 mm or above with submillimeter accuracy; the films were only able to detect errors > 1.0 mm.

CONCLUSION

We developed a quantitative technique for the PF test on the Elekta Unity. The merged films covered nearly the entire MLC leaf banks. The technique exhibited clinically acceptable accuracy and sensitivity to MLC positioning errors.

Technical note: Cryostat transmission characterization for MR linac - temporal stability, clinical impact and change implementation

BACKGROUND

In the Unity MR linac (Elekta AB, Stockholm, Sweden), the radiation beam traverses the cryostat and the coil support structure. The resulting beam attenuation must be considered for output calibration and its variation with gantry angle must be characterized in the treatment planning system (TPS).

PURPOSE

The aim of this work was to investigate the impact of a change of the cryostat transmission characterization (CTC) curve, due to the helium level modification, on clinical treatment plan dosimetry and to report on the experience with the CTC curve update.

METHODS

Twenty stereotactic body radiotherapy (SBRT) treatment plans: 10 prostate and 10 oligo-metastatic cancer plans, prepared with a beam model incorporating the CTC curve acquired at installation time, were re-calculated using the model implementing CTC curve post helium top-up. To account for the CTC change as well as to align our system to the recent reference conditions recommendations, the new model was commissioned with the emphasis on the specifics associated with the treatment plan adaptation and the existence of the offline and online TPS components.

RESULTS

Average CTV mean dose reduction by 0.45% in prostate cases and average GTV mean dose reduction by 0.22% in oligo-metastatic cases was observed. Updated model validation showcased good agreement between measurements and TPS calculations.

CONCLUSIONS

The agreement between CTC measurements demonstrates its temporal constancy and robustness of the measurement method employed. A helium fill level change was shown to affect the CTC and led to a small but systematic dose calculation inaccuracy. Finally, model validation and end-to-end testing results presented, underscore the minimal impact of transitioning to the new beam model and new reference conditions.

Repeatability quantification of brain diffusion-weighted imaging for future clinical implementation at a low-field MR-linac

BACKGROUND

Longitudinal assessments of apparent diffusion coefficients (ADCs) derived from diffusion-weighted imaging (DWI) during intracranial radiotherapy at magnetic resonance imaging-guided linear accelerators (MR-linacs) could enable early response assessment by tracking tumor diffusivity changes. However, DWI pulse sequences are currently unavailable in clinical practice at low-field MR-linacs. Quantifying the in vivo repeatability of ADC measurements is a crucial step towards clinical implementation of DWI sequences but has not yet been reported on for low-field MR-linacs. This study assessed ADC measurement repeatability in a phantom and in vivo at a 0.35 T MR-linac.

METHODS

Eleven volunteers and a diffusion phantom were imaged on a 0.35 T MR-linac. Two echo-planar imaging DWI sequence variants, emphasizing high spatial resolution ("highRes") and signal-to-noise ratio ("highSNR"), were investigated. A test-retest study with an intermediate outside-scanner-break was performed to assess repeatability in the phantom and volunteers' brains. Mean ADCs within phantom vials, cerebrospinal fluid (CSF), and four brain tissue regions were compared to literature values. Absolute relative differences of mean ADCs in pre- and post-break scans were calculated for the diffusion phantom, and repeatability coefficients (RC) and relative RC (relRC) with 95% confidence intervals were determined for each region-of-interest (ROI) in volunteers.

RESULTS

Both DWI sequence variants demonstrated high repeatability, with absolute relative deviations below 1% for water, dimethyl sulfoxide, and polyethylene glycol in the diffusion phantom. RelRCs were 7% [5%, 12%] (CSF; highRes), 12% [9%, 22%] (CSF; highSNR), 9% [8%, 12%] (brain tissue ROIs; highRes), and 6% [5%, 7%] (brain tissue ROIs; highSNR), respectively. ADCs measured with the highSNR variant were consistent with literature values for volunteers, while smaller mean values were measured for the diffusion phantom. Conversely, the highRes variant underestimated ADCs compared to literature values, indicating systematic deviations.

CONCLUSIONS

High repeatability of ADC measurements in a diffusion phantom and volunteers' brains were measured at a low-field MR-linac. The highSNR variant outperformed the highRes variant in accuracy and repeatability, at the expense of an approximately doubled voxel volume. The observed high in vivo repeatability confirms the potential utility of DWI at low-field MR-linacs for early treatment response assessment.

Potential dosimetric error in the adaptive workflow of a 1.5 T MR-Linac from patient movement relative to immobilisation systems

In magnetic resonance- (MR-) based adaptive workflows for an MR-linac, the treatment plan is optimized and recalculated online using the daily MR images. The Unity MR-linac is supplied with a patient positioning device (ppd) using pelvic and abdomen thermoplastic masks attached to a board with high-density components. This study highlights the dosimetric effect of using this in such workflows when there are relative patient-ppd displacements, as these are not visualized on MR imaging and the treatment planning system assumes the patient is fixed relative to the ppd. The online adapted plans of two example rectum cancer patients treated at a Unity MR-linac were perturbed by introducing relative patient-ppd displacements, and the effect was evaluated on plan dosimetry. Forty-eight perturbed clinical adapted plans were recalculated, based on online MR-based synthetic computed tomography, and compared with the original plans, using dose-volume histogram parameters and gamma analysis. The target volume covered by the prescribed dose ( D pre ) and by at least 107% of D pre varied up to - 1.87% and + 3.67%, respectively for 0.5 cm displacements, and to - 3.18% and + 4.96% for 2 cm displacements; whilst 2%-2 mm gamma analysis showed a median value of 92.9%. The use of a patient positioning system with high-density components in a Unity MR-based online adaptive treatment workflow can introduce unrecognized errors in plan dosimetry and it is recommended not to use such a device for such treatments, without modifying the device and the workflow, followed by careful clinical evaluation, or alternatively to use other immobilization methods.

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.

MR-guided stereotactic radiotherapy of infra-diaphragmatic oligometastases: Evaluation of toxicity and dosimetric parameters

BACKGROUND AND PURPOSE

The SOFT trial is a prospective, multicenter, phase 2 trial investigating magnetic resonance (MR)-guided stereotactic ablative radiotherapy (SABR) for abdominal, soft tissue metastases in patients with oligometastatic disease (OMD) (clinicaltrials.gov ID NCT04407897). We present the primary endpoint analysis of 1-year treatment-related toxicity (TRAE).

MATERIALS AND METHODS

Patients with up to five oligometastases from non-hematological cancers were eligible for inclusion. A risk-adapted strategy prioritized fixed organs at risk (OAR) constraints over target coverage. Fractionation schemes were 45-67.5 Gy in 3-8 fractions. The primary endpoint was grade ≥ 4 TRAE within 12 months post-SABR. The association between the risk of gastrointestinal (GI) toxicity and clinical and dosimetric parameters was tested using a normal tissue complication probability model.

RESULTS

We included 121 patients with 147 oligometastatic targets, mainly located in the liver (41 %), lymph nodes (35 %), or adrenal glands (14 %). Nearly half of all targets (48 %, n = 71) were within 10 mm of a radiosensitive OAR. No grade 4 or 5 TRAEs, 3.5 % grade 3 TRAEs, and 43.7 % grade 2 TRAEs were reported within the first year of follow-up. We found a significant association between grade ≥ 2 GI toxicity and the parameters GI OAR D0.1cc, D1cc, and D20cc.

CONCLUSION

In this phase II study of MR-guided SABR of oligometastases in the infra-diaphragmatic region, we found a low incidence of toxicity despite half of the lesions being within 10 mm of a radiosensitive OAR. GI OAR D0.1cc, D1cc, and D20cc were associated with grade ≥ 2 GI toxicity.

Is it necessary to perform measurement-based patient-specific quality assurance for online adaptive radiotherapy with Elekta Unity MR-Linac?

This study aimed to investigate the necessity of measurement-based patient-specific quality assurance (PSQA) for online adaptive radiotherapy by analyzing measurement-based PSQA results and calculation-based 3D independent dose verification results with Elekta Unity MR-Linac. There are two workflows for Elekta Unity enabled in the treatment planning system: adapt to position (ATP) and adapt to shape (ATS). ATP plans are those which have relatively slighter shifts from reference plans by adjusting beam shapes or weights, whereas ATS plans are the new plans optimized from the beginning with probable re-contouring targets and organs-at-risk. PSQA gamma passing rates were measured using an MR-compatible ArcCHECK diode array for 78 reference plans and corresponding 208 adaptive plans (129 ATP plans and 79 ATS plans) of Elekta Unity. Subsequently, the relationships between ATP, or ATS plans and reference plans were evaluated separately. The Pearson's r correlation coefficients between ATP or ATS adaptive plans and corresponding reference plans were also characterized using regression analysis. Moreover, the Bland-Altman plot method was used to describe the agreement of PSQA results between ATP or ATS adaptive plans and reference plans. Additionally, Monte Carlo-based independent dose verification software ArcherQA was used to perform secondary dose check for adaptive plans. For ArcCHECK measurements, the average gamma passing rates (ArcCHECK vs. TPS) of PSQA (3%/2 mm criterion) were 99.51% ± 0.88% and 99.43% ± 0.54% for ATP and ATS plans, respectively, which were higher than the corresponding reference plans 99.34% ± 1.04% (p < 0.05) and 99.20% ± 0.71% (p < 0.05), respectively. The Pearson's r correlation coefficients were 0.720 between ATP and reference plans and 0.300 between ATS and reference plans with ArcCHECK, respectively. Furthermore, >95% of data points of differences between both ATP and ATS plans and reference plans were within ±2σ (standard deviation) of the mean difference between adaptive and reference plans with ArcCHECK measurements. With ArcherQA calculation, the average gamma passing rates (ArcherQA vs. TPS) were 98.23% ± 1.64% and 98.15% ± 1.07% for ATP and ATS adaptive plans, separately. It might be unnecessary to perform measurement-based PSQA for both ATP and ATS adaptive plans for Unity if the gamma passing rates of both measurements of corresponding reference plans and independent dose verification of adaptive plans have high gamma passing rates. Periodic machine QA and verification of adaptive plans were recommended to ensure treatment safety.

Dose Calculation Accuracy of Beam Models in RadCalc for a 1.5 T MR-Linac

The purpose of this study is to evaluate RadCalc, an independent dose verification software, for patient-specific quality assurance (PSQA) in online adaptive planning with a magnetic resonance linear accelerator (MR-linac) of a 1.5 T. Version 7.1.4 of RadCalc to introduce the capability to establish a beam model that incorporates MR field characteristics. A total of six models were established, with one using manufacturer-provided data and the others differing in percentage depth dose (PDD) data sources. Overall, two models utilized PDD data from the treatment planning system (TPS), and three used commissioned PDD data from gantry angles of 0° and 270°. Simple tests on a virtual water phantom assessed dose-calculation accuracy, revealing percentage differences ranging from -0.5% to -20.6%. Excluding models with significant differences, clinical tests on 575 adaptive plans (prostate, liver, and breast) showed percentage differences of -0.51%, 1.12%, and 4.10%, respectively. The doses calculated using RadCalc demonstrated similar trends to those of the PSQA-based measurements. The newly released version of RadCalc enables beam modeling that considers the characteristics of the 1.5 T magnetic field. The accuracy of the software in calculating doses at 1.5 T magnetic fields has been verified, thereby making it a reliable and effective tool for PSQA in adaptive plans.

A virtual phantom for patient-specific QA On A 1.5T MR-linac

Create a virtual ArcCHECK-MR phantom, customized for a 1.5T MR-linac, with consideration of the different density regions within the quality assurance (QA) phantom, aiming to streamline the utilization of this specialized QA device. A virtual phantom was constructed in the treatment planning system (TPS) to replicate the ArcCHECK-MR's composition, consisting of five distinct layers: "Outer" (representing the outer PMMA ring), "Complex" (simulating the printed circuit boards), "Detectors" (encompassing the detector area), "Inner" (signifying the inner PMMA ring) and "Insert" (representing the PMMA insert). These layers were defined based on geometric data and represented as contour points on a set of dummy CT images. Additionally, a setup platform was integrated as contoured structures. To determine the relative electron density (RED) values of the external and internal PMMA components, measurements were taken at 25 points in the insert using an ion chamber. A novel method for establishing the exit/entrance dose ratio (EEDR) for ArcCHECK-MR was introduced. The RED of higher density region was derived by evaluating the local gamma index passing rate results with criteria of 2% dose difference and 2 mm distance-to-agreement. The performance of the virtual phantom was assessed for Unity 7 FFF beams with a 1.5T magnetic field. The radii of the five ring structures within the virtual phantom measured 133.0 mm, 110.0 mm, 103.4 mm, 100.0 mm, and 75.0 mm for the "Outer," "Complex," "Detectors," "Inner" and "Insert" regions, respectively. The RED values were as follows: ArcCHECK-MR PMMA had a RED of 1.130, "Detectors" were assumed to have a RED of 1.000, "Complex" had a RED of 1.200, and the setup QA phantom justified a RED of 1.350. Early validation results demonstrate that the 5-layer virtual phantom, when compared to the commonly used bulk overridden phantom, offers improved capability in MR-linac environments. This enhancement led to an increase in passing rates for the local gamma index by approximately 5 ∼ 6%, when applying the criteria of 2%, 2 mm. We have successfully generated a virtual representation of the distinct regions within the ArcCHECK-MR using a TPS, addressing the challenges associated with its use in conjunction with a 1.5T MR-linac. We consistently observed favorable local gamma index passing rates across two 1.5T MR-linac and ArcCHECK-MR unit combinations. This approach has the potential to minimize uncertainties in the creation of the QA phantom for ArcCHECK-MR across various institutions.

Dosimetric comparison of glioblastoma radiotherapy treatment plans on a low-field MRI-guided linear accelerator compared to conventional C-arm linear accelerator

INTRODUCTION

Magnetic resonance imaging (MRI)-guided radiation therapy has proven to provide many benefits such as real-time tracking, dose escalation, and the ability to perform online adaptive therapy. The objective of this study is to compare curative treatment plans for glioblastoma tumors on a low-field MR-guided linac vs a C-arm linac and evaluate if they are comparable in terms of coverage, organ at risk sparing, delivery time, and deliverability.

METHODS

This is a retrospective study that consisted of 15 previously treated patients who received radiation therapy for glioblastoma on a C-arm linac. The CT simulation data used for the original clinical plans was imported into the MR-linac treatment planning system (TPS) and utilized for MR-linac plan generation. The plans were evaluated utilizing the dose volumetric histogram (DVH) and isodose lines, then compared in terms of plan quality consisting of PTV coverage, dose distributions, and OAR constraints. Statistical analysis was performed to compare differences between the two planning techniques. QA was performed on a subset of the plans to verify deliverability.

RESULTS

Plans generated on the MR-linac were more heterogenous compared to C-arm linac plans. A statistically significant difference was found in the homogeneity index (HI) and the PTV V105% volume (cc) values. The volume of the normal brain receiving 30 Gy also showed a statistically significant (p = 0.0479) difference, where on average an additional 41.5 cc's of the normal brain tissue received 30 Gy in the MR-linac plans. The maximum dose to the normal brain structure also increased in the MR-linac plans on average by 2.6 Gy (p = 0.0002). Similarly, the average maximum dose to the scalp 4 mm structure was 6.5 Gy higher in the MR-linac plans compared to C-arm linac plans (p = 0.0103). The total MU's were higher in the MR-linac plans compared to the C-arm linac plans (p = 0.0015).

CONCLUSIONS

Both MR-linac and C-arm linac plans met constraints for PTV coverage and OAR sparing, were deliverable, and resulted to be clinically acceptable. However, our study showed that MR-linac plans were not as conformal or as homogenous as C-arm linac plans utilizing noncoplanar beams.

Simulating an intra-fraction adaptive workflow to enable PTV margin reduction in MRIgART volumetric modulated arc therapy for prostate SBRT

Purpose

This study simulates a novel prostate SBRT intra-fraction re-optimization workflow in MRIgART to account for prostate intra-fraction motion and evaluates the dosimetric benefit of reducing PTV margins.

Materials and methods

VMAT prostate SBRT treatment plans were created for 10 patients using two different PTV margins, one with a 5 mm margin except 3 mm posteriorly (standard) and another using uniform 2 mm margins (reduced). All plans were prescribed to 36.25 Gy in 5 fractions and adapted onto each daily MRI dataset. An intra-fraction adaptive workflow was simulated for the reduced margin group by synchronizing the radiation delivery with target position from cine MRI imaging. Intra-fraction delivered dose was reconstructed and prostate DVH metrics were evaluated under three conditions for the reduced margin plans: Without motion compensation (no-adapt), with a single adapt prior to treatment (ATP), and lastly for intra-fraction re-optimization during delivery (intra). Bladder and rectum DVH metrics were compared between the standard and reduced margin plans.

Results

As expected, rectum V18 Gy was reduced by 4.4 ± 3.9%, D1cc was reduced by 12.2 ± 6.8% (3.4 ± 2.3 Gy), while bladder reductions were 7.8 ± 5.6% for V18 Gy, and 9.6 ± 7.3% (3.4 ± 2.5 Gy) for D1cc for the reduced margin reference plans compared to the standard PTV margin. For the intrafraction replanning approach, average intra-fraction optimization times were 40.0 ± 2.9 seconds, less than the time to deliver one of the four VMAT arcs (104.4 ± 9.3 seconds) used for treatment delivery. When accounting for intra-fraction motion, prostate V36.25 Gy was on average 96.5 ± 4.0%, 99.1 ± 1.3%, and 99.6 ± 0.4 for the non-adapt, ATP, and intra-adapt groups, respectively. The minimum dose received by the prostate was less than 95% of the prescription dose in 84%, 36%, and 10% of fractions, for the non-adapt, ATP, and intra-adapt groups, respectively.

Conclusions

Intra-fraction re-optimization improves prostate coverage, specifically the minimum dose to the prostate, and enables PTV margin reduction and subsequent OAR sparing. Fast re-optimizations enable uninterrupted treatment delivery.

Prostate volume variation during 1.5T MR-guided adaptive stereotactic body radiotherapy (SBRT) and correlation with treatment toxicity

INTRODUCTION

To evaluate prostate volume change during daily-adaptive prostate SBRT on 1.5 T MR-linac and to correlate it with treatment toxicity.

METHODS

a series of patients affected by low-to-intermediate risk prostate cancer was treated by 5-fraction SBRT within a prospective study (Prot. n° 23748). Total dose was 35 Gy and 36.25 Gy delivered every day or on alternate days. Treatment toxicity was recorded with the following patient reported outcomes (PROMs): IPSS, ICIQ-SF, and EPIC-26.

RESULTS

254 patients were included in the analysis. Baseline median CTV volume was 55 cc (range 15.3-163.3). Mean prostate volume were 58.9 cc, and 62.7 cc at first and last fraction respectively (mean volume increase 6.4 %; p = <0.0001). We observed prostate swelling (mean 15.4 % increase) in 50 % of cases, stable volume (≤5% volume change) in 39 % of patients, and prostate shrinkage in 11 % of cases (mean 12.2 % reduction). Baseline CTV > 55 cc showed a trend towards higher CTV shrinkage (-10.5 % versus -14.5 %; p = 0.052). We found no correlation between CTV change and PROMs. Prostate swelling was generally compensated by the planned PTV expansion, even though the mean setup volume dropped from 47.4 cc to 38.9 cc at last fraction, with few cases not covered by initial setup margins.

CONCLUSION

The present study reported a significant prostate volume change during prostate SBRT on 1.5T MR-linac. We observed both prostate swelling in half of cases and few cases of prostate shrinkage. No correlations were found with PROMs in this population treatment with daily-adaptive strategy.

The impact of margin reduction on radiation dose distribution of ultra-hypofractionated prostate radiotherapy utilizing a 1.5-T MR-Linac

BACKGROUND

We examined the effects of reducing the planning target volume (PTV) margin in MR-guided radiotherapy (MRgRT) on the distribution of radiation dose to target volumes and organs-at-risk (OARs). Thus, we compared MR-Linac (MRL) plans with and without reduced margin and intensity-modulated radiotherapy (IMRT) plan with conventional linac for low-risk prostate cancer patients receiving 36.25 Gy in five fractions of ultra-hypofractionated radiation therapy.

MATERIALS AND METHODS

Twenty low-risk prostate cancer patients treated with 1.5 T MR-Linac were evaluated. The same planning CT images were used for four plans: the MRL-R plan with reduced margin planning target volume (PTV-R) and the MRL-N plan with normal margin PTV (PTV-N), which is also used for IMRT plan. In four plans, PTV doses, organs-at-risk (OARs) doses, the homogeneity index (HI), and monitor units were compared.

RESULTS

All plans met the criteria for PTV coverage and OARs dose constraints. The maximum and mean PTV doses were significantly higher in the MRL-R and MRL-N plans compared to the IMRT plan. The HI was lowest in the IMRT plan (0.040 ± 0.013) and highest in the MRL-N plan (0.055 ± 0.012; p < 0.001). There was no significant difference in the PTV dosimetric parameters between the MRL-R and the MRL-N plans. The high doses in the rectum was significantly lower in the MRL-R compared to other plans. The bladder V36.25 Gy was significantly lower in the MRL-R plan (2.43 ± 1.87 Gy) compared to MRL-N (4.50 ± 2.42 Gy; p < 0.001), and IMRT plans (4.76 ± 2.77 Gy; p < 0.001). There was no significant difference in the low-dose volumes of the body, maximum femur doses, or monitor units across each plan.

CONCLUSIONS

Ultra-hypofractionated MR-guided RT with 1.5 T MRL is dosimetrically feasible for patients with prostate cancer. The improved soft tissue contrast and the online adaptive plan for 1.5 T MR-Linac allows for PTV margin reduction resulted in a significant dose reduction in OARs.

Effectiveness of visual biofeedback-guided respiratory-correlated 4D-MRI for radiotherapy guidance on the MR-linac

PURPOSE

Respiratory-correlated 4D-MRI may provide motion characteristics for radiotherapy but is susceptible to irregular breathing. This study investigated the effectiveness of visual biofeedback (VBF) guidance for breathing regularization during 4D-MRI acquisitions on an MR-linac.

METHODS

A simultaneous multislice-accelerated 4D-MRI sequence was interleaved with a one-dimensional respiratory navigator (1D-RNAV) in 10 healthy volunteers on a 1.5T Unity MR-linac (Elekta AB, Stockholm, Sweden). Volunteer-specific breathing amplitudes and periods were derived from the 1D-RNAV signal obtained during unguided 4D-MRI acquisitions. These were used for the guidance waveform, while the 1D-RNAV positions were overlayed as VBF. VBF effectiveness was quantified by calculating the change in coefficient of variation (CV diff $$ {\mathrm{CV}}^{\mathrm{diff}} $$ ) for the breathing amplitude and period, the position SD of end-exhale, end-inhale and midposition locations, and the agreement between the 1D-RNAV signals and guidance waveforms. The 4D-MRI quality was assessed by quantifying amounts of missing data.

RESULTS

VBF had an average latency of 520 ± 2 ms. VBF reduced median breathing variations by 18% to 35% (amplitude) and 29% to 57% (period). Median position SD reductions ranged from -3% to 35% (end-exhale), 29% to 38% (end-inhale), and 25% to 37% (midposition). Average differences between guidance waveforms and 1D-RNAV signals were 0.0 s (period) and +1.7 mm (amplitude). VBF also decreased the median amount of missing data by 11% and 29%.

CONCLUSION

A VBF system was successfully implemented, and all volunteers were able to adapt to the guidance waveform. VBF during 4D-MRI acquisitions drastically reduced breathing variability but had limited effect on missing data in respiratory-correlated 4D-MRI.

Investigation of autosegmentation techniques on T2-weighted MRI for off-line dose reconstruction in MR-linac workflow for head and neck cancers

BACKGROUND

In order to accurately accumulate delivered dose for head and neck cancer patients treated with the Adapt to Position workflow on the 1.5T magnetic resonance imaging (MRI)-linear accelerator (MR-linac), the low-resolution T2-weighted MRIs used for daily setup must be segmented to enable reconstruction of the delivered dose at each fraction.

PURPOSE

In this pilot study, we evaluate various autosegmentation methods for head and neck organs at risk (OARs) on on-board setup MRIs from the MR-linac for off-line reconstruction of delivered dose.

METHODS

Seven OARs (parotid glands, submandibular glands, mandible, spinal cord, and brainstem) were contoured on 43 images by seven observers each. Ground truth contours were generated using a simultaneous truth and performance level estimation (STAPLE) algorithm. Twenty total autosegmentation methods were evaluated in ADMIRE: 1-9) atlas-based autosegmentation using a population atlas library (PAL) of 5/10/15 patients with STAPLE, patch fusion (PF), random forest (RF) for label fusion; 10-19) autosegmentation using images from a patient's 1-4 prior fractions (individualized patient prior [IPP]) using STAPLE/PF/RF; 20) deep learning (DL) (3D ResUNet trained on 43 ground truth structure sets plus 45 contoured by one observer). Execution time was measured for each method. Autosegmented structures were compared to ground truth structures using the Dice similarity coefficient, mean surface distance (MSD), Hausdorff distance (HD), and Jaccard index (JI). For each metric and OAR, performance was compared to the inter-observer variability using Dunn's test with control. Methods were compared pairwise using the Steel-Dwass test for each metric pooled across all OARs. Further dosimetric analysis was performed on three high-performing autosegmentation methods (DL, IPP with RF and 4 fractions [IPP_RF_4], IPP with 1 fraction [IPP_1]), and one low-performing (PAL with STAPLE and 5 atlases [PAL_ST_5]). For five patients, delivered doses from clinical plans were recalculated on setup images with ground truth and autosegmented structure sets. Differences in maximum and mean dose to each structure between the ground truth and autosegmented structures were calculated and correlated with geometric metrics.

RESULTS

DL and IPP methods performed best overall, all significantly outperforming inter-observer variability and with no significant difference between methods in pairwise comparison. PAL methods performed worst overall; most were not significantly different from the inter-observer variability or from each other. DL was the fastest method (33 s per case) and PAL methods the slowest (3.7-13.8 min per case). Execution time increased with a number of prior fractions/atlases for IPP and PAL. For DL, IPP_1, and IPP_RF_4, the majority (95%) of dose differences were within ± 250 cGy from ground truth, but outlier differences up to 785 cGy occurred. Dose differences were much higher for PAL_ST_5, with outlier differences up to 1920 cGy. Dose differences showed weak but significant correlations with all geometric metrics (R2 between 0.030 and 0.314).

CONCLUSIONS

The autosegmentation methods offering the best combination of performance and execution time are DL and IPP_1. Dose reconstruction on on-board T2-weighted MRIs is feasible with autosegmented structures with minimal dosimetric variation from ground truth, but contours should be visually inspected prior to dose reconstruction in an end-to-end dose accumulation workflow.

Impact of intrafraction motion in pancreatic cancer treatments with MR-guided adaptive radiation therapy

Purpose

The total time of radiation treatment delivery for pancreatic cancer patients with daily online adaptive radiation therapy (ART) on an MR-Linac can range from 50 to 90 min. During this period, the target and normal tissues undergo changes due to respiration and physiologic organ motion. We evaluated the dosimetric impact of the intrafraction physiological organ changes.

Methods

Ten locally advanced pancreatic cancer patients were treated with 50 Gy in five fractions with intensity-modulated respiratory-gated radiation therapy on a 0.35-T MR-Linac. Patients received both pre- and post-treatment volumetric MRIs for each fraction. Gastrointestinal organs at risk (GI-OARs) were delineated on the pre-treatment MRI during the online ART process and retrospectively on the post-treatment MRI. The treated dose distribution for each adaptive plan was assessed on the post-treatment anatomy. Prescribed dose volume histogram metrics for the scheduled plan on the pre-treatment anatomy, the adapted plan on the pre-treatment anatomy, and the adapted plan on post-treatment anatomy were compared to the OAR-defined criteria for adaptation: the volume of the GI-OAR receiving greater than 33 Gy (V33Gy) should be ≤1 cubic centimeter.

Results

Across the 50 adapted plans for the 10 patients studied, 70% were adapted to meet the duodenum constraint, 74% for the stomach, 12% for the colon, and 48% for the small bowel. Owing to intrafraction organ motion, at the time of post-treatment imaging, the adaptive criteria were exceeded for the duodenum in 62% of fractions, the stomach in 36%, the colon in 10%, and the small bowel in 48%. Compared to the scheduled plan, the post-treatment plans showed a decrease in the V33Gy, demonstrating the benefit of plan adaptation for 66% of the fractions for the duodenum, 95% for the stomach, 100% for the colon, and 79% for the small bowel.

Conclusion

Post-treatment images demonstrated that over the course of the adaptive plan generation and delivery, the GI-OARs moved from their isotoxic low-dose region and nearer to the dose-escalated high-dose region, exceeding dose-volume constraints. Intrafraction motion can have a significant dosimetric impact; therefore, measures to mitigate this motion are needed. Despite consistent intrafraction motion, plan adaptation still provides a dosimetric benefit.

Evaluation of Unity 1.5 T MR-linac plan quality in patients with prostate cancer

The Unity magnetic resonance (MR) linear accelerator (MRL) with MR-guided adaptive radiotherapy (MRgART) is capable of online MRgART where images are acquired on the treatment day and the radiation treatment plan is immediately replanned and performed. We evaluated the MRgART plan quality and plan reproducibility of the Unity MRL in patients with prostate cancer. There were five low- or moderate-risk and five high-risk patients who received 36.25 Gy or 40 Gy, respectively in five fractions. All patients underwent simulation magnetic resonance imaging (MRI) and five online adaptive MRI. We created plans for 5, 7, 9, 16, and 20 beams and for 60, 100, and 150 segments. We evaluated the target and organ doses for different number of beams and segments, respectively. Variation in dose constraint between the simulation plan and online adaptive plan was measured for each patient to assess plan reproducibility. The plan quality improved with the increasing number of beams. However, the proportion of significantly improved dose constraints decreased as the number of beams increased. For some dose parameters, there were statistically significant differences between 60 and 100 segments, and 100 and 150 segments. The plan of five beams exhibited limited reproducibility. The number of segments had minimal impact on plan reproducibility, but 60 segments sometimes failed to meet dose constraints for online adaptive plan. The optimization and delivery time increased with the number of beams and segments. We do not recommend using five or fewer beams for a reproducible and high-quality plan in the Unity MRL. In addition, many number of segments and beams may help meet dose constraints during online adaptive plan. Treatment with the Unity MRL should be performed with the appropriate number of beams and segments to achieve a good balance among plan quality, delivery time, and optimization time.

Continuous time-resolved estimated synthetic 4D-CTs for dose reconstruction of lung tumor treatments at a 0.35 T MR-linac

Objective.To experimentally validate a method to create continuous time-resolved estimated synthetic 4D-computed tomography datasets (tresCTs) based on orthogonal cine MRI data for lung cancer treatments at a magnetic resonance imaging (MRI) guided linear accelerator (MR-linac).Approach.A breathing porcine lung phantom was scanned at a CT scanner and 0.35 T MR-linac. Orthogonal cine MRI series (sagittal/coronal orientation) at 7.3 Hz, intersecting tumor-mimicking gelatin nodules, were deformably registered to mid-exhale 3D-CT and 3D-MRI datasets. The time-resolved deformation vector fields were extrapolated to 3D and applied to a reference synthetic 3D-CT image (sCTref), while accounting for breathing phase-dependent lung density variations, to create 82 s long tresCTs at 3.65 Hz. Ten tresCTs were created for ten tracked nodules with different motion patterns in two lungs. For each dataset, a treatment plan was created on the mid-exhale phase of a measured ground truth (GT) respiratory-correlated 4D-CT dataset with the tracked nodule as gross tumor volume (GTV). Each plan was recalculated on the GT 4D-CT, randomly sampled tresCT, and static sCTrefimages. Dose distributions for corresponding breathing phases were compared in gamma (2%/2 mm) and dose-volume histogram (DVH) parameter analyses.Main results.The mean gamma pass rate between all tresCT and GT 4D-CT dose distributions was 98.6%. The mean absolute relative deviations of the tresCT with respect to GT DVH parameters were 1.9%, 1.0%, and 1.4% for the GTVD98%,D50%, andD2%, respectively, 1.0% for the remaining nodulesD50%, and 1.5% for the lungV20Gy. The gamma pass rate for the tresCTs was significantly larger (p< 0.01), and the GTVD50%deviations with respect to the GT were significantly smaller (p< 0.01) than for the sCTref.Significance.The results suggest that tresCTs could be valuable for time-resolved reconstruction and intrafractional accumulation of the dose to the GTV for lung cancer patients treated at MR-linacs in the future.

Dose-response dependencies of OSL dosimeters in conventional linacs and 1.5T MR-linacs: an experimental and Monte Carlo study

Objective. This work focuses on the optically stimulated luminescence dosimetry (OSLD) dose-response characterization, with emphasis on 1.5T MR-Linacs.Approach. Throughout this study, the nanoDots OSLDs (Landauer, USA) were considered. In groups of three, the mean OSLD response was measured in a conventional linac and an MR-Linac under various irradiation conditions to investigate (i) dose-response linearity with and without the 1.5T magnetic field, (ii) signal fading rate and its dependencies, (iii) beam quality, detector orientation and dose rate dependencies in a conventional linac, (iii) potential MR imaging related effects on OSLD response and (iv) detector orientation dependence in an MR-Linac. Monte Carlo calculations were performed to further quantify angular dependence after rotating the detector around its central axis parallel to the magnetic field, and determine the magnetic field correction factors,kB,Q,for all cardinal detector orientations.Main results. OSLD dose-response supralinearity in an MR-Linac setting was found to agree within uncertainties with the corresponding one in a conventional linac, for the axial detector orientation investigated. Signal fading rate does not depend on irradiation conditions for the range of 3-30 d considered. OSLD angular (orientation) dependence is more pronounced under the presence of a magnetic field. OSLDs irradiated with and without real-time T2w MR imaging enabled during irradiation yielded the same response within uncertainties.kB,Qvalues were determined for all three cardinal orientations. Corrections needed reached up to 6.4%. However, if OSLDs are calibrated in the axial orientation and then irradiated in an MR-Linac placed again in the axial orientation (perpendicular to the magnetic field), then simulations suggest thatkB,Qcan be considered unity within uncertainties, irrespective of the incident beam angle.Significance. This work contributes towards OSLD dose-response characterization and relevant correction factors availability. OSLDs are suitable for QA checks in MR-based beam gating applications andin vivodosimetry in MR-Linacs.

Experimental validation of multi-fraction online adaptations in magnetic resonance guided radiotherapy

Background and purpose

Radiotherapy plan verification is generally performed on the reference plan based on the pre-treatment anatomy. However, the introduction of online adaptive treatments demands a new approach, as plans are created daily on different anatomies. The aim of this study was to experimentally validate the accuracy of total doses of multi-fraction plan adaptations in magnetic resonance imaging guided radiotherapy in a phantom study, isolated from the uncertainty of deformable image registration.

Materials and methods

We experimentally verified the total dose, measured on external beam therapy 3 (EBT3) film, using a treatment with five online adapted fractions. Three series of experiments were performed, each focusing on a category of inter-fractional variation; translations, rotations and body modifications. Variations were introduced during each fraction and adapted plans were generated and irradiated. Single fraction doses and total doses over five online adapted fractions were investigated.

Results

The online adapted measurements and calculations showed a good agreement for single fractions and multi-fraction treatments for the dose profiles, gamma passing rates, dose deviations and distances to agreement. The gamma passing rate using a 2%/2 mm criterion ranged from 99.2% to 99.5% for a threshold dose of 10% of the maximum dose (Dmax) and from 96.2% to 100% for a threshold dose of 90% of Dmax, for the total translations, rotations and body modifications.

Conclusions

The total doses of multi-fraction treatments showed similar accuracies compared to single fraction treatments, indicating an accurate dosimetric outcome of a multi-fraction treatment in adaptive magnetic resonance imaging guided radiotherapy.

Accumulated dose implications from systematic dose-rate transients in gated treatments with Viewray MRIdian accelerators

The combination of magnetic resonance (MR) imaging and linear accelerators (linacs) into MR-Linacs enables continuous MR imaging and advanced gated treatments of patients. Previously, a dose-rate transient (∼8% reduced dose rate during the initial 0.5 s of each beam) was identified for a Viewray MRIdian MR-Linac (Klavsenet al2022Radiation Measurement106759). Here, the dose-rate transient is studied in more detail at four linacs of the same type at different hospitals. The implications of dose-rate transients were examined for gated treatments. The dose-rate transients were investigated using dose-per pulse measurements with organic plastic scintillators in three experiments: (i) A gated treatment with the scintillator placed in a moving target in a dynamic phantom, (ii) a gated treatment with the same dynamic conditions but with the scintillator placed in a stationary target, and (iii) measurements in a water-equivalent material to examine beam quality deviations at a dose-per-pulse basis. Gated treatments (i) compared with non-gated treatments with a static target in the same setup showed a broadening of accumulated dose profiles due to motion (dose smearing). The linac with the largest dose-rate transient had a reduced accumulated dose of up to (3.1 ± 0.65) % in the center of the PTV due to the combined dose smearing and dose-rate transient effect. Dose-rate transients were found to vary between different machines. Two MR-Linacs showed initial dose-rate transients that could not be identified from conventional linearity tests. The source of the transients includes an initial change in photon fluence rate and an initial change in x-ray beam quality. For gated treatments, this caused a reduction of more than 1% dose delivered at the central part of the beam for the studied, cyclic-motion treatment plan. Quality assurance of this effect should be considered when gated treatment with the Viewray MRIdian is implemented clinically.

Performance assessment of variant UNet-based deep-learning dose engines for MR-Linac-based prostate IMRT plans

Objective. UNet-based deep-learning (DL) architectures are promising dose engines for traditional linear accelerator (Linac) models. Current UNet-based engines, however, were designed differently with various strategies, making it challenging to fairly compare the results from different studies. The objective of this study is to thoroughly evaluate the performance of UNet-based models on magnetic-resonance (MR)-Linac-based intensity-modulated radiation therapy (IMRT) dose calculations.Approach. The UNet-based models, including the standard-UNet, cascaded-UNet, dense-dilated-UNet, residual-UNet, HD-UNet, and attention-aware-UNet, were implemented. The model input is patient CT and IMRT field dose in water, and the output is patient dose calculated by DL model. The reference dose was calculated by the Monaco Monte Carlo module. Twenty training and ten test cases of prostate patients were included. The accuracy of the DL-calculated doses was measured using gamma analysis, and the calculation efficiency was evaluated by inference time.Results. All the studied models effectively corrected low-accuracy doses in water to high-accuracy patient doses in a magnetic field. The gamma passing rates between reference and DL-calculated doses were over 86% (1%/1 mm), 98% (2%/2 mm), and 99% (3%/3 mm) for all the models. The inference times ranged from 0.03 (graphics processing unit) to 7.5 (central processing unit) seconds. Each model demonstrated different strengths in calculation accuracy and efficiency; Res-UNet achieved the highest accuracy, HD-UNet offered high accuracy with the fewest parameters but the longest inference, dense-dilated-UNet was consistently accurate regardless of model levels, standard-UNet had the shortest inference but relatively lower accuracy, and the others showed average performance. Therefore, the best-performing model would depend on the specific clinical needs and available computational resources.Significance. The feasibility of using common UNet-based models for MR-Linac-based dose calculations has been explored in this study. By using the same model input type, patient training data, and computing environment, a fair assessment of the models' performance was present.

Real-time myocardial landmark tracking for MRI-guided cardiac radio-ablation using Gaussian Processes

The high speed of cardiorespiratory motion introduces a unique challenge for cardiac stereotactic radio-ablation (STAR) treatments with the MR-linac. Such treatments require tracking myocardial landmarks with a maximum latency of 100 ms, which includes the acquisition of the required data. The aim of this study is to present a new method that enables tracking myocardial landmarks from few readouts of MRI data, thereby achieving a latency sufficient for STAR treatments. We present a tracking framework that requires few readouts of k-space data as input, which can be acquired at least an order of magnitude faster than MR-images. Combined with the real-time tracking speed of a probabilistic machine learning framework called Gaussian Processes, this allows to track myocardial landmarks with a sufficiently low latency for cardiac STAR guidance. This includes both the acquisition of required data, and the tracking inference. The framework is demonstrated in 2D on a motion phantom, and in vivo on volunteers and a ventricular tachycardia (arrhythmia) patient. Moreover, the feasibility of an extension to 3D was demonstrated by in silico 3D experiments with a digital motion phantom. The framework was compared with template matching - a reference, image-based, method - and linear regression methods. Results indicate an order of magnitude lower total latency (<10 ms) for the proposed framework in comparison with alternative methods. The root-mean-square-distances and mean end-point-distance with the reference tracking method was less than 0.8 mm for all experiments, showing excellent (sub-voxel) agreement. The high accuracy in combination with a total latency of less than 10 ms - including data acquisition and processing - make the proposed method a suitable candidate for tracking during STAR treatments. Additionally, the probabilistic nature of the Gaussian Processes also gives access to real-time prediction uncertainties, which could prove useful for real-time quality assurance during treatments.

Feasibility study of adaptive radiotherapy for esophageal cancer using artificial intelligence autosegmentation based on MR-Linac

Objective

We proposed a scheme for automatic patient-specific segmentation in Magnetic Resonance (MR)-guided online adaptive radiotherapy based on daily updated, small-sample deep learning models to address the time-consuming delineation of the region of interest (ROI) in the adapt-to-shape (ATS) workflow. Additionally, we verified its feasibility in adaptive radiation therapy for esophageal cancer (EC).

Methods

Nine patients with EC who were treated with an MR-Linac were prospectively enrolled. The actual adapt-to-position (ATP) workflow and simulated ATS workflow were performed, the latter of which was embedded with a deep learning autosegmentation (AS) model. The first three treatment fractions of the manual delineations were used as input data to predict the next fraction segmentation, which was modified and then used as training data to update the model daily, forming a cyclic training process. Then, the system was validated in terms of delineation accuracy, time, and dosimetric benefit. Additionally, the air cavity in the esophagus and sternum were added to the ATS workflow (producing ATS+), and the dosimetric variations were assessed.

Results

The mean AS time was 1.40 [1.10-1.78 min]. The Dice similarity coefficient (DSC) of the AS model gradually approached 1; after four training sessions, the DSCs of all ROIs reached a mean value of 0.9 or more. Furthermore, the planning target volume (PTV) of the ATS plan showed a smaller heterogeneity index than that of the ATP plan. Additionally, V5 and V10 in the lungs and heart were greater in the ATS+ group than in the ATS group.

Conclusion

The accuracy and speed of artificial intelligence-based AS in the ATS workflow met the clinical radiation therapy needs of EC. This allowed the ATS workflow to achieve a similar speed to the ATP workflow while maintaining its dosimetric advantage. Fast and precise online ATS treatment ensured an adequate dose to the PTV while reducing the dose to the heart and lungs.

Precision radiotherapy using MR-linac for pancreatic neuroendocrine tumors in MEN1 patients (PRIME): a protocol for a phase I-II trial, and systematic review on available evidence for radiotherapy of pNETs

Background

Surgical resection is the standard of care for the treatment of pancreatic neuro-endocrine tumors (pNETs) in patients with Multiple Endocrine Neoplasia Type 1 (MEN1). However, surgery can cause significant short- and long-term morbidity. Magnetic resonance-guided radiotherapy (MRgRT) is a potential effective treatment with little side effects. With traditional radiotherapy techniques, irradiation of pancreatic tumors to high dose levels was hampered by poor visibility of the tumor during treatment. MRgRT uses onboard MRI to guide the treatment, thereby enabling delivery of ablative irradiation doses to the tumor, while sparing surrounding tissues. In this study, we describe results from a systematic review assessing efficacy of radiotherapy in pNET and present the protocol of the PRIME study.

Methods

PubMed, Embase and Cochrane Library were searched for articles assessing efficacy and side effects of radiotherapy for the treatment of pNETs. Risk of bias was assessed using the ROBINS-I Risk of Bias Tool for observational studies. Descriptive statistics were used to describe results of included trials.

Results

Four studies comprising of 33 patients treated by conventional radiotherapy were included. Despite the heterogeneity of studies, radiotherapy appeared to be effective for the treatment of pNETs with most patients responding (45.5%) or stabilizing (42.4%) in tumor size.

Conclusion and trial design

Due to the limited literature available and concerns about damage to surrounding tissue, conventional radiotherapy is currently little used for pNETs. The PRIME study is a phase I-II trial with a single arm prospective cohort study design, investigating the efficacy of MRgRT in MEN1 patients with pNET. MEN1 patients with growing pNETs with a size between 1.0 and 3.0 cm without malignant features are eligible for inclusion. Patients are treated with 40 Gy in 5 fractions on the pNET, using online adaptive MRgRT on a 1.5T MR-linac. The primary endpoint is the change in tumor size at MRI 12 months follow-up. Secondary endpoints include radiotoxicity, quality of life, endocrine and exocrine pancreas function, resection rate, metastatic free and overall survival. When MRgRT is found effective with low radiotoxicity, it could reduce the need for surgery for pNET and preserve quality of life.

Systematic Review Registration

PROSPERO https://clinicaltrials.gov/, (CRD42022325542).

Gaussian Processes for real-time 3D motion and uncertainty estimation during MR-guided radiotherapy

Respiratory motion during radiotherapy causes uncertainty in the tumor's location, which is typically addressed by an increased radiation area and a decreased dose. As a result, the treatments' efficacy is reduced. The recently proposed hybrid MR-linac scanner holds the promise to efficiently deal with such respiratory motion through real-time adaptive MR-guided radiotherapy (MRgRT). For MRgRT, motion-fields should be estimated from MR-data and the radiotherapy plan should be adapted in real-time according to the estimated motion-fields. All of this should be performed with a total latency of maximally 200 ms, including data acquisition and reconstruction. A measure of confidence in such estimated motion-fields is highly desirable, for instance to ensure the patient's safety in case of unexpected and undesirable motion. In this work, we propose a framework based on Gaussian Processes to infer 3D motion-fields and uncertainty maps in real-time from only three readouts of MR-data. We demonstrated an inference frame rate up to 69 Hz including data acquisition and reconstruction, thereby exploiting the limited amount of required MR-data. Additionally, we designed a rejection criterion based on the motion-field uncertainty maps to demonstrate the framework's potential for quality assurance. The framework was validated in silico and in vivo on healthy volunteer data (n=5) acquired using an MR-linac, thereby taking into account different breathing patterns and controlled bulk motion. Results indicate end-point-errors with a 75th percentile below 1 mm in silico, and a correct detection of erroneous motion estimates with the rejection criterion. Altogether, the results show the potential of the framework for application in real-time MR-guided radiotherapy with an MR-linac.

Influence of different factors on registration error in a 1.5 T MR-guided linac

Purpose. Accurate image registration is an important step in online image-guided adaptive radiotherapy. The aim of this study was to investigate the effects of different factors on registration accuracy in a magnetic resonance (MR)-guided adaptive radiotherapy workflow.Materials and Methods. A thorax motion phantom was used to obtain computed tomography (CT) simulations in 8 different motion modes and to generate 8 reference plans. Daily pretreatment online MR images were obtained at 5 different positions in each reference plan. Online MR and CT simulations were separately registered using bone structures and the gross tumor volume (GTV) as ROIs, and the image shift distance was recorded by the online treatment planning system. The difference between the shift distance and the real isocentric distance was the registration error. The registration error was analyzed, and the effects of the setup position, motion mode and ROI selection on the registration error were investigated by multivariate analysis of variance.Result. The minimum values of registration error (ΔX, ΔY, ΔZ) were -1.90 mm, -2.70 mm and -2.40 mm, respectively, and the maximum values were 1.70 mm, 4.30 mm and -0.90 mm. ΔY showed the maximum mean standard deviation of 1.25 mm, and ΔZshowed the minimum mean standard deviation of 0.27 mm. The standard deviation of the registration error is largest in the inferior/superior direction. The motion mode of the phantom and ROI selection were significantly correlated with ΔX, ΔY, and ΔZ(p< 0.05).Conclusion. The registration result with the spine as the selected ROI was better than that with the GTV as the ROI. In 1.5 T MR-linac clinical treatment, more attention should be given to patient movement repeatability and to controlling the intrafractional motion as much as possible. It is not recommended to make the GTV-PTV margin expansion less than 2 mm for MR-linac.

Online adaptive radiotherapy and dose delivery accuracy: A retrospective analysis

PURPOSE

With online adaptive radiotherapy (ART), patient-specific quality assurance (PSQA) testing cannot be performed prior to delivery of the adapted treatment plan. Consequently, the dose delivery accuracy of adapted plans (i.e., the ability of the system to interpret and deliver the treatment as planned) are not initially verified. We investigated the variation in dose delivery accuracy of ART on the MRIdian 0.35 T MR-linac (Viewray Inc., Oakwood, USA) between initial plans and their respective adapted plans, by analyzing PSQA results.

METHODS

We considered the two main digestive localizations treated with ART (liver and pancreas). A total of 124 PSQA results acquired with the ArcCHECK (Sun Nuclear Corporation, Melbourne, USA) multidetector system were analyzed. PSQA result variations between the initial plans and their respective adapted plans were statistically investigated and compared with the variation in MU number.

RESULTS

For the liver, limited deterioration in PSQA results was observed, and was within the limits of clinical tolerance (Initial = 98.2%, Adapted = 98.2%, p = 0.4503). For pancreas plans, only a few significant deteriorations extending beyond the limits of clinical tolerance were observed and were due to specific, complex anatomical configurations (Initial = 97.3%, Adapted = 96.5%, p = 0.0721). In parallel, we observed an influence of the increase in MU number on the PSQA results.

CONCLUSION

We show that the dose delivery accuracy of adapted plans, in terms of PSQA results, is preserved in ART processes on the 0.35 T MR-linac. Respecting good practices, and minimizing the increase in MU number can help to preserve the accuracy of delivery of adapted plans as compared to their respective initial plans.

Linac-based versus MR-guided SBRT for localized prostate cancer: a comparative evaluation of acute tolerability

AIM

This study aims to compare acute toxicity of prostate cancer (PCa) stereotactic body radiotherapy (SBRT) delivered by MR-guided radiotherapy (MRgRT) with 1.5-T MR-linac or by volumetric modulated arc (VMAT) with conventional linac.

METHODS

Patients with low-to-favorable intermediate risk class PCa were treated with exclusive SBRT (35 Gy in five fractions). Patients treated with MRgRT were enrolled in an Ethical Committee (EC) approved trial (Prot. n° 23,748), while patients treated with conventional linac were enrolled in an EC approved phase II trial (n° SBRT PROG112CESC). The primary end-point was the acute toxicity. Patients were included in the analysis if they had at least 6 months of follow-up for the primary end-point evaluation. Toxicity assessment was performed according to CTCAE v5.0 scale. International Prostatic Symptoms Score (IPSS) was also performed.

RESULTS

A total of 135 patients were included in the analysis. Seventy-two (53.3%) were treated with MR-linac and 63 (46.7%) with conventional linac. The median initial PSA before RT was 6.1 ng/ml (range 0.49-19). Globally, acute G1, G2, and G3 toxicity occurred in 39 (28.8%), 20 (14.5%), and 5 (3.7%) patients. At the univariate analysis, acute G1 toxicity did not differ between MR-linac and conventional linac (26.4% versus 31.8%), as well as G2 toxicity (12.5% versus 17.5%; p = 0.52). Acute G2 gastrointestinal (GI) toxicity occurred in 7% and 12.5% of cases in MR-linac and conventional linac group, respectively (p = 0.06), while acute G2 genitourinary toxicity occurred in 11% and 12.8% in MR-linac and conventional linac, respectively (p = 0.82). The median IPSS before and after SBRT was 3 (1-16) and 5 (1-18). Acute G3 toxicity occurred in two cases in the MR-linac and three cases in the conventional linac group (p = n.s.).

CONCLUSION

Prostate SBRT with 1.5-T MR-linac is feasible and safe. Compared to conventional linac, MRgRT might to potentially reduce the overall G1 acute toxicity at 6 months, and seems to show a trend toward a lower incidence of grade 2 GI toxicity. A longer follow-up is necessary to assess the late efficacy and toxicity.

High-resolution entry and exit surface dosimetry in a 1.5 T MR-linac

The magnetic field of a transverse MR-linac alters electron trajectories as the photon beam transits through materials, causing lower doses at flat entry surfaces and increased doses at flat beam-exiting surfaces. This study investigated the response of a MOSFET detector, known as the MOSkin™, for high-resolution surface and near-surface percentage depth dose measurements on an Elekta Unity. Simulations with Geant4 and the Monaco treatment planning system (TPS), and EBT-3 film measurements, were also performed for comparison. Measured MOSkin™ entry surface doses, relative to D_max, were (9.9 ± 0.2)%, (10.1 ± 0.3)%, (11.3 ± 0.6)%, (12.9 ± 1.0)%, and (13.4 ± 1.0)% for 1 × 1 cm², 3 × 3 cm², 5 × 5 cm², 10 × 10 cm², and 22 × 22 cm² fields, respectively. For the investigated fields, the maximum percent differences of Geant4, TPS, and film doses extrapolated and interpolated to a depth suitable for skin dose assessment at the beam entry, relative to MOSkin™ measurements at an equivalent depth were 1.0%, 2.8%, and 14.3%, respectively, and at a WED of 199.67 mm at the beam exit, 3.2%, 3.7% and 5.7%, respectively. The largest measured increase in exit dose, due to the electron return effect, was 15.4% for the 10 × 10 cm² field size using the MOSkin™ and 17.9% for the 22 × 22 cm² field size, using Geant4 calculations. The results presented in the study validate the suitability of the MOSkin™ detector for transverse MR-linac surface dosimetry.

Experimental and Monte Carlo-based determination of magnetic field correction factors k B , Q $k_{B,Q}$ in high-energy photon fields for two ionization chambers

BACKGROUND

The integration of magnetic resonance tomography into clinical linear accelerators provides high-contrast, real-time imaging during treatment and facilitates online-adaptive workflows in radiation therapy treatments. The associated magnetic field also bends the trajectories of charged particles via the Lorentz force, which may alter the dose distribution in a patient or a phantom and affects the dose response of dosimetry detectors.

PURPOSE

To perform an experimental and Monte Carlo-based determination of correction factors k B , Q $k_{B,Q}$ , which correct the response of ion chambers in the presence of external magnetic fields in high-energy photon fields.

METHODS

The response variation of two different types of ion chambers (Sun Nuclear SNC125c and SNC600c) in strong external magnetic fields was investigated experimentally and by Monte Carlo simulations. The experimental data were acquired at the German National Metrology Institute, PTB, using a clinical linear accelerator with a nominal photon energy of 6 MV and an external electromagnet capable of generating magnetic flux densities of up to 1.5 T in opposite directions. The Monte Carlo simulation geometries corresponded to the experimental setup and additionally to the reference conditions of IAEA TRS-398. For the latter, the Monte Carlo simulations were performed with two different photon spectra: the 6 MV spectrum of the linear accelerator used for the experimental data acquisition and a 7 MV spectrum of a commercial MRI-linear accelerator. In each simulation geometry, three different orientations of the external magnetic field, the beam direction and the chamber orientation were investigated.

RESULTS

Good agreement was achieved between Monte Carlo simulations and measurements with the SNC125c and SNC600c ionization chambers, with a mean deviation of 0.3% and 0.6%, respectively. The magnitude of the correction factor k B , Q $k_{B,Q}$ strongly depends on the chamber volume and on the orientation of the chamber axis relative to the external magnetic field and the beam directions. It is greater for the SNC600c chamber with a volume of 0.6 cm³ than for the SNC125c chamber with a volume of 0.1 cm³ . When the magnetic field direction and the chamber axis coincide, and they are perpendicular to the beam direction, the ion chambers exhibit a calculated overresponse of less than 0.7(6)% (SNC600c) and 0.3(4)% (SNC125c) at 1.5 T and less than 0.3(0)% (SNC600c) and 0.1(3)% (SNC125c) for 0.35 T for nominal beam energies of 6 MV and 7 MV. This chamber orientation should be preferred, as k B , Q $k_{B,Q}$ may increase significantly in other chamber orientations. Due to the special geometry of the guard ring, no dead-volume effects have been observed in any orientation studied. The results show an intra-type variation of 0.17% and 0.07% standard uncertainty (k=1) for the SNC125c and SNC600c, respectively.

CONCLUSION

Magnetic field correction factors k B , Q $k_{B,Q}$ for two different ion chambers and for typical clinical photon beam qualities were presented and compared with the few data existing in the literature. The correction factors may be applied in clinical reference dosimetry for existing MRI-linear accelerators.

Evaluation of MRI-only based online adaptive radiotherapy of abdominal region on MR-linac

PURPOSE

A hybrid magnetic resonance linear accelerator (MRL) can perform magnetic resonance imaging (MRI) with high soft-tissue contrast to be used for online adaptive radiotherapy (oART). To obtain electron densities needed for the oART dose calculation, a computed tomography (CT) is often deformably registered to MRI. Our aim was to evaluate an MRI-only based synthetic CT (sCT) generation as an alternative to the deformed CT (dCT)-based oART in the abdominal region.

METHODS

The study data consisted of 57 patients who were treated on a 0.35 T MRL system mainly for abdominal tumors. Simulation MRI-CT pairs of 43 patients were used for training and validation of a prototype convolutional neural network sCT-generation algorithm, based on HighRes3DNet, for the abdominal region. For remaining test patients, sCT images were produced from simulation MRIs and daily MRIs. The dCT-based plans were re-calculated on sCT with identical calculation parameters. The sCT and dCT were compared in terms of geometric agreement and calculated dose.

RESULTS

The mean and one standard deviation of the geometric agreement metrics over dCT-sCT-pairs were: mean error of 8 ± 10 HU, mean absolute error of 49 ± 10 HU, and Dice similarity coefficient of 55 ± 12%, 60 ± 5%, and 82 ± 15% for bone, fat, and lung tissues, respectively. The dose differences between the sCT and dCT-based dose for planning target volumes were 0.5 ± 0.9%, 0.6 ± 0.8%, and 0.5 ± 0.8% at D_2% , D_50% , and D_98% in physical dose and 0.8 ± 1.4%, 0.8 ± 1.2%, and 0.6 ± 1.1% in biologically effective dose (BED). For organs-at-risk, the dose differences of all evaluated dose-volume histogram points were within [-4.5%, 7.8%] and [-1.1 Gy, 3.5 Gy] in both physical dose and BED.

CONCLUSIONS

The geometric agreement metrics were within typically reported values and most average relative dose differences were within 1%. Thus, an MRI-only sCT-based approach is a promising alternative to the current clinical practice of the abdominal oART on MRL.

Adaptive radiotherapy for breast cancer

Research in the field of local and locoregional breast cancer radiotherapy aims to maintain excellent oncological outcomes while reducing treatment-related toxicity. Adaptive radiotherapy (ART) considers variations in target and organs at risk (OARs) anatomy occurring during the treatment course and integrates these in re-optimized treatment plans. Exploiting ART routinely in clinic may result in smaller target volumes and better OAR sparing, which may lead to reduction of acute as well as late toxicities. In this review MR-guided and CT-guided ART for breast cancer patients according to different clinical scenarios (neoadjuvant and adjuvant partial breast irradiation, whole breast, chest wall and regional nodal irradiation) are reviewed and their advantages as well as challenging aspects discussed.

Comparison of MR-guided radiotherapy accumulated doses for central lung tumors with non-adaptive and online adaptive proton therapy

BACKGROUND

Stereotactic body radiation therapy (SBRT) of central lung tumors with photon or proton therapy has a risk of increased toxicity. Treatment planning studies comparing accumulated doses for state-of-the-art treatment techniques, such as MR-guided radiotherapy (MRgRT) and intensity modulated proton therapy (IMPT), are currently lacking.

PURPOSE

We conducted a comparison of accumulated doses for MRgRT, robustly optimized non-adaptive IMPT, and online adaptive IMPT for central lung tumors. A special focus was set on analyzing the accumulated doses to the bronchial tree, a parameter linked to high-grade toxicities.

METHODS

Data of 18 early-stage central lung tumor patients, treated at a 0.35 T MR-linac in eight or five fractions, were analyzed. Three gated treatment scenarios were compared: (S1) online adaptive MRgRT, (S2) non-adaptive IMPT, and (S3) online adaptive IMPT. The treatment plans were recalculated or reoptimized on the daily imaging data acquired during MRgRT, and accumulated over all treatment fractions. Accumulated dose-volume histogram (DVH) parameters of the gross tumor volume (GTV), lung, heart, and organs-at-risk (OARs) within 2 cm of the planning target volume (PTV) were extracted for each scenario and compared in Wilcoxon signed-rank tests between S1 & S2, and S1 & S3.

RESULTS

The accumulated GTV D_98% was above the prescribed dose for all patients and scenarios. Significant reductions (p < 0.05) of the mean ipsilateral lung dose (S2: -8%; S3: -23%) and mean heart dose (S2: -79%; S3: -83%) were observed for both proton scenarios compared to S1. The bronchial tree D_0.1cc was significantly lower for S3 (S1: 48.1 Gy; S3: 39.2 Gy; p = 0.005), but not significantly different for S2 (S2: 45.0 Gy; p = 0.094), compared to S1. The D_0.1cc for S2 and S3 compared to S1 was significantly (p < 0.05) smaller for OARs within 1-2 cm of the PTV (S1: 30.2 Gy; S2: 24.6 Gy; S3: 23.1 Gy), but not significantly different for OARs within 1 cm of the PTV.

CONCLUSIONS

A significant dose sparing potential of non-adaptive and online adaptive proton therapy compared to MRgRT for OARs in close, but not direct proximity of central lung tumors was identified. The near-maximum dose to the bronchial tree was not significantly different for MRgRT and non-adaptive IMPT. Online adaptive IMPT achieved significantly lower doses to the bronchial tree compared to MRgRT.

Evaluation of real-time tumor contour prediction using LSTM networks for MR-guided radiotherapy

BACKGROUND AND PURPOSE

Magnetic resonance imaging guided radiotherapy (MRgRT) with deformable multileaf collimator (MLC) tracking would allow to tackle both rigid displacement and tumor deformation without prolonging treatment. However, the system latency must be accounted for by predicting future tumor contours in real-time. We compared the performance of three artificial intelligence (AI) algorithms based on long short-term memory (LSTM) modules for the prediction of 2D-contours 500ms into the future.

MATERIALS AND METHODS

Models were trained (52 patients, 3.1h of motion), validated (18 patients, 0.6h) and tested (18 patients, 1.1h) with cine MRs from patients treated at one institution. Additionally, we used three patients (2.9h) treated at another institution as second testing set. We implemented 1) a classical LSTM network (LSTM-shift) predicting tumor centroid positions in superior-inferior and anterior-posterior direction which are used to shift the last observed tumor contour. The LSTM-shift model was optimized both in an offline and online fashion. We also implemented 2) a convolutional LSTM model (ConvLSTM) to directly predict future tumor contours and 3) a convolutional LSTM combined with spatial transformer layers (ConvLSTM-STL) to predict displacement fields used to warp the last tumor contour.

RESULTS

The online LSTM-shift model was found to perform slightly better than the offline LSTM-shift and significantly better than the ConvLSTM and ConvLSTM-STL. It achieved a 50% Hausdorff distance of 1.2mm and 1.0mm for the two testing sets, respectively. Larger motion ranges were found to lead to more substantial performance differences across the models.

CONCLUSION

LSTM networks predicting future centroids and shifting the last tumor contour are the most suitable for tumor contour prediction. The obtained accuracy would allow to reduce residual tracking errors during MRgRT with deformable MLC-tracking.

Understanding the Benefit of Magnetic Resonance-guided Adaptive Radiotherapy in Rectal Cancer Patients: a Single-centre Study

AIMS

Neoadjuvant chemoradiotherapy followed by surgery is the mainstay of treatment for patients with rectal cancer. Standard clinical target volume (CTV) to planning target volume (PTV) margins of 10 mm are used to accommodate inter- and intrafraction motion of target. Treating on magnetic resonance-integrated linear accelerators (MR-linacs) allows for online manual recontouring and adaptation (MRgART) enabling the reduction of PTV margins. The aim of this study was to investigate motion of the primary CTV (CTVA; gross tumour volume and macroscopic nodes with 10 mm expansion to cover microscopic disease) in order to develop a simultaneous integrated boost protocol for use on MR-linacs.

MATERIALS AND METHODS

Patients suitable for neoadjuvant chemoradiotherapy were recruited for treatment on MR-linac using a two-phase technique; only the five phase 1 fractions on MR-linac were used for analysis. Intrafraction motion of CTVA was measured between pre-treatment and post-treatment MRI scans. In MRgART, isotropically expanded pre-treatment PTV margins from 1 to 10 mm were rigidly propagated to post-treatment MRI to determine overlap with 95% of CTVA. The PTV margin was considered acceptable if overlap was >95% in 90% of fractions. To understand the benefit of MRgART, the same methodology was repeated using a reference computed tomography planning scan for pre-treatment imaging.

RESULTS

In total, nine patients were recruited between January 2018 and December 2020 with T3a-T4, N0-N2, M0 disease. Forty-five fractions were analysed in total. The median motion across all planes was 0 mm, demonstrating minimal intrafraction motion. A PTV margin of 3 and 5mm was found to be acceptable in 96 and 98% of fractions, respectively. When comparing to the computed tomography reference scan, the analysis found that PTV margins to 5 and 10 mm only acceptably covered 51 and 76% of fractions, respectively.

CONCLUSION

PTV margins can be reduced to 3-5 mm in MRgART for rectal cancer treatment on MR-linac within an simultaneous integrated boost protocol.

Dose accumulation for MR-guided adaptive radiotherapy: From practical considerations to state-of-the-art clinical implementation

MRI-linear accelerator (MR-linac) devices have been introduced into clinical practice in recent years and have enabled MR-guided adaptive radiation therapy (MRgART). However, by accounting for anatomical changes throughout radiation therapy (RT) and delivering different treatment plans at each fraction, adaptive radiation therapy (ART) highlights several challenges in terms of calculating the total delivered dose. Dose accumulation strategies-which typically involve deformable image registration between planning images, deformable dose mapping, and voxel-wise dose summation-can be employed for ART to estimate the delivered dose. In MRgART, plan adaptation on MRI instead of CT necessitates additional considerations in the dose accumulation process because MRI pixel values do not contain the quantitative information used for dose calculation. In this review, we discuss considerations for dose accumulation specific to MRgART and in relation to current MR-linac clinical workflows. We present a general dose accumulation framework for MRgART and discuss relevant quality assurance criteria. Finally, we highlight the clinical importance of dose accumulation in the ART era as well as the possible ways in which dose accumulation can transform clinical practice and improve our ability to deliver personalized RT.

MRI of radiation chemistry: First images and investigation of potential mechanisms

BACKGROUND

Paramagnetic species such as O2 and free radicals can enhance T1 and T2 relaxation times. If the change in relaxation time is sufficiently large, the contrast will be generated in magnetic resonance images. Since radiation is known to be capable of altering the concentration of O2 and free radicals during water radiolysis, it may be possible for radiation to induce MR signal change.

PURPOSE

We present the first reported instance of x-ray-induced MR signal changes in water phantoms and investigate potential paramagnetic relaxation enhancement mechanisms associated with radiation chemistry changes in oxygen and free radical concentrations.

METHODS

Images of water and 10 mM coumarin phantoms were acquired on a 0.35 T MR-linac before, during, and after a dose delivery of 80 Gy using an inversion-recovery dual-echo sequence with water nullified. Radiation chemistry simulations of these conditions were performed to calculate changes in oxygen and free radical concentrations. Published relaxivity values were then applied to calculate the resulting T1 change, and analytical MR signal equations were used to calculate the associated signal change.

RESULTS

Compared to pre-irradiation reference images, water phantom images taken during and after irradiation showed little to no change, while coumarin phantom images showed a small signal loss in the irradiated region with a contrast-to-noise ratio (CNR) of 1.0-2.5. Radiation chemistry simulations found oxygen depletion of -11 µM in water and -31 µM in coumarin, resulting in a T1 lengthening of 24 ms and 68 ms respectively, and a simulated CNR of 1.0 and 2.8 respectively. This change was consistent with observations in both direction and magnitude. Steady-state superoxide, hydroxyl, hydroperoxyl, and hydrogen radical concentrations were found to contribute less than 1 ms of T1 change.

CONCLUSION

Observed radiation-induced MR signal changes were dominated by an oxygen depletion mechanism. Free radicals were concluded to play a minor secondary role under steady-state conditions. Future applications may include in vivo FLASH treatment verification but would require an MR sequence with a better signal-to-noise ratio and higher temporal resolution than the one used in this study.

Intra-fraction motion of pelvic oligometastases and feasibility of PTV margin reduction using MRI guided adaptive radiotherapy

Purpose

This study assesses the impact of intra-fraction motion and PTV margin size on target coverage for patients undergoing radiation treatment of pelvic oligometastases. Dosimetric sparing of the bowel as a function of the PTV margin is also evaluated.

Materials and methods

Seven patients with pelvic oligometastases previously treated on our MR-linac (35 Gy in 5 fractions) were included in this study. Retrospective adaptive plans were created for each fraction on the daily MRI datasets using PTV margins of 5 mm, 3 mm, and 2 mm. Dosimetric constraint violations and GTV coverage were measured as a function of PTV margin size. The impact of intra-fraction motion on GTV coverage was assessed by tracking the GTV position on the cine MR images acquired during treatment delivery and creating an intra-fraction dose distribution for each IMRT beam. The intra-fraction dose was accumulated for each fraction to determine the total dose delivered to the target for each PTV size.

Results

All OAR constraints were achieved in 85.7%, 94.3%, and 100.0% of fractions when using 5 mm, 3 mm, and 2 mm PTV margins while scaling to 95% PTV coverage. Compared to plans with a 5 mm PTV margin, there was a 27.4 ± 12.3% (4.0 ± 2.2 Gy) and an 18.5 ± 7.3% (2.7 ± 1.4 Gy) reduction in the bowel D_0.5cc dose for 2 mm and 3 mm PTV margins, respectively. The target dose (GTV V_{35 Gy}) was on average 100.0 ± 0.1% (99.6 - 100%), 99.6 ± 1.0% (97.2 - 100%), and 99.0 ± 1.4% (95.0 - 100%), among all fractions for the 5 mm, 3 mm, and 2 mm PTV margins on the adaptive plans when accounting for intra-fraction motion, respectively.

Conclusion

A 2 mm PTV margin achieved a minimum of 95% GTV coverage while reducing the dose to the bowel for all patients.

On the feasibility of cardiac substructure sparing in magnetic resonance imaging guided stereotactic lung radiotherapy

BACKGROUND

Lung stereotactic body radiotherapy (SBRT) has proven an effective treatment for medically inoperable lung tumors, even for (ultra-)central tumors. Recently, there has been growing interest in radiation-induced cardiac toxicity in lung radiotherapy. More specifically, dose to cardiac (sub-)structures (CS) was found to correlate with survival after radiotherapy.

PURPOSE

Our goal is first, to investigate the percentage of patients who require CS sparing in an magnetic resonance imaging guided lung SBRT workflow, and second, to quantify how successful implementation of cardiac sparing would be.

METHODS

The patient cohort consists of 34 patients with stage II-IV lung cancer who were treated with SBRT between 2017 and 2020. A mid-position computed tomography (CT) image was used to create treatment plans for the 1.5 T Unity MR-linac (Elekta AB, Stockholm, Sweden) following clinical templates. Under guidance of a cardio-thoracic radiologist, 11 CS were contoured manually for each patient. Dose constraints for five CS were extracted from the literature. Patients were stratified according to their need for cardiac sparing depending on the CS dose in their non-CS constrained MR-linac treatment plans. Cardiac sparing treatment plans (CSPs) were then created and dosimetrically compared with their non-CS constrained treatment plan counterparts. CSPs complied with the departmental constraints and were considered successful when fulfilling all CS constraints, and partially successful if some CS constraints could be fulfilled. Predictors for the need for and feasibility of cardiac sparing were explored, specifically planning target volume (PTV) size, cranio-caudal (CC) distance, 3D distance, and in-field overlap volume histograms (iOVH).

RESULTS

47% of the patients (16 out of 34) were in need of cardiac sparing. A successful CSP could be created for 62.5% (10 out of 16) of these patients. Partially successful CSPs still complied with two to four CS constraints. No significant difference in dose to organs at risk (OARs) or targets was identified between CSPs and the corresponding non-CS constrained MR-linac plans. The need for cardiac sparing was found to correlate with distance in the CC direction between target and all of the individual CS (Mann-Whitney U-test p-values <10^-6 ). iOVHs revealed that complying with dose constraints for CS is primarily determined by in-plane distance and secondarily by PTV size.

CONCLUSION

We demonstrated that CS can be successfully spared in lung SBRT on the MR-linac for most of this patient cohort, without compromising doses to the tumor or to other OARs. CC distance between the target and CS can be used to predict the need for cardiac sparing. iOVHs, in combination with PTV size, can be used to predict if cardiac sparing will be successful for all constrained CS except the left ventricle.

Characterisation of a split gradient coil design induced systemic imaging artefact on 0.35 T MR-linac systems

Objective. The aim of this work was to highlight and characterize a systemic 'star-like' artefact inherent to the low field 0.35 T MRIdian MR-linac system, a magnetic resonance guided radiotherapy device. This artefact is induced by the original split gradients coils design. This design causes a surjection of the intensity gradient inZ(or head-feet) direction. This artefact appears on every sequence with phase encoding in the head-feet direction.Approach. Basic gradient echo sequence and clinical mandatory bSSFP sequence were used. Three setups using manufacturer provided QA phantoms were designed: two including the linearity control grid used for the characterisation and a third including two homogeneity control spheres dedicated to the artefact management in a more clinical like situation. The presence of the artefact was checked in four different MRidian sites. The tested parameters based on the literature were: phase encoding orientation, slab selectivity, excitation bandwidth (BW_RF), acceleration factor (R) and phase/slab oversampling (PO/SO).Main results. The position of this artefact is constant and reproducible over the tested MRIdian sites. The typical singularity saturated dot or star is visible even with the 3D slab-selection enabled. A management is proposed by decreasing the BW_RF, theRin head-feet direction and increasing the PO/SO. The oversampling can be optimized using a formula to anticipate the location of artefact in the field of view.Significance. The star-like artefact has been well characterised. A manageable solution comes at the cost of acquisition time. Observed in clinical cases, the artefact may degrade the images used for the RT planning and repositioning during the treatment unless corrected.

1.5T MR-Guided Daily-Adaptive SBRT for Prostate Cancer: Preliminary Report of Toxicity and Quality of Life of the First 100 Patients

Purpose: The present study reports the preliminary outcomes in terms of adverse events and quality of life in the first 100 patients treated with 1.5T MR-guided daily-adaptive stereotactic body radiotherapy for prostate cancer. Methods: From October 2019 to December 2020, 100 patients, enrolled in a prospective study, received MR-guided SBRT for prostate cancer. Rectal spacer insertion was optional and administered in 37 patients. In total, 32 patients received androgen deprivation therapy in accordance with international guidelines. A prospective collection of data regarding toxicity and quality of life was performed. Results: The median age was 71 years (range, 52-84). The median total dose delivered was 35 Gy (35-36.25 Gy) in five sessions, either on alternate days (n = 25) or consecutive days (n = 75). For acute toxicity, we recorded: seven cases of acute G2 urinary pain and four cases of G2 gastrointestinal events. The median follow-up was 12 months (3-20), recording three late G2 urinary events and one G3 case, consisting of a patient who required a TURP 8 months after the treatment. For gastrointestinal toxicity, we observed 3 G ≥ 2 GI events, including one patient who received argon laser therapy for radiation-induced proctitis. Up to the last follow-up, all patients are alive and with no evidence of biochemical relapse, except for an M1 low-volume patient in distant progression two months after radiotherapy. QoL evaluation reported a substantial resolution of any discomfort within the second follow-up after radiotherapy, with the only exception being sexual items. Notably, after one year, global health items were improved compared to the baseline assessment. Conclusions: This study reports very promising outcomes in terms of adverse events and QoL, supporting the role of 1.5T MR-guided SBRT for prostate cancer. To date, this series is one of the first and largest available in the literature. Long-term results are warranted.

Tuning the optimal diffusion-weighted MRI parameters on a 0.35-T MR-Linac for clinical implementation: A phantom study

PURPOSE

This study aims to assess the quality of a new diffusion-weighted imaging (DWI) sequence implemented on an MR-Linac MRIdian system, evaluating and optimizing the acquisition parameters to explore the possibility of clinically implementing a DWI acquisition protocol in a 0.35-T MR-Linac.

MATERIALS AND METHODS

All the performed analyses have been carried out on two types of phantoms: a homogeneous 24-cm diameter polymethylmethacrylate (PMMA) sphere (SP) and a homemade phantom (HMP) constating in a PMMA cylinder filled with distilled water with empty sockets into which five cylindrical vials filled with five different concentrations of methylcellulose water solutions have been inserted. SP was used to evaluate the dependence of diffusion gradient inhomogeneity artifacts on gantry position. Four diffusion sequences with b-values of 500 s/mm2 and 3 averages have been acquired: three with diffusion gradients in the three main directions (phase direction, read direction, slice direction) and one with the diffusion gradients switched off. The dependence of diffusion image uniformity and SNR on the number of averages in the MR sequences was also investigated to determine the optimal number of averages. Finally, the ADC values of HMP have been computed and then compared between images acquired in the scanners at 0.35 and 1.5 T.

RESULTS

In order to acquire high-quality artifact-free DWI images, the "slice" gradient direction has been identified to be the optimal one and 0° to be the best gradient angle. Both the SNR ratio and the uniformity increase with the number of averages. A threshold value of 80 for SNR and 85% for uniformity was adopted to choose the best number of averages. By making a compromise between time and quality and limiting the number of b-values, it is possible to reduce the acquisition time to 78 s. The Passing-Bablok test showed that the two methods, with 0.35 and 1.5 T scanners, led to similar results.

CONCLUSION

The quality of the DWI has been accurately evaluated in relation to different sequence parameters, and optimal parameters have been identified to select a clinical protocol for the acquisition of ADC maps sustainable in the workflow of a hybrid radiotherapy system with a 0.35-T MRI scanner.

1.5 T MR-Guided Daily Adapted SBRT on Lymph Node Oligometastases from Prostate Cancer

Introduction: The aim of our study was to evaluate the efficacy and toxicity of a daily adaptive MR-guided SBRT on 1.5 T MR-linac in patients affected by lymph node oligometastases from PCa. Materials and Methods: The present study is a prospective observational study conducted in a single institution (protocol n°: MRI/LINAC n. 23748). Patients with oligometastatic lymph nodes from PCa treated with daily adaptive MR-guided SBRT on 1.5 T MR-linac were included in the study. There was a minimum required follow-up of 3 months after SBRT. The primary end-point was local progression-free survival (LPFS). The secondary end-points were: nodal progression-free survival (NPFS), progression-free survival (PFS), and toxicity. Results: A total of 118 lymph node oligometastases from PCa were treated with daily adaptive 1.5 T MR-guided SBRT in 63 oligometastatic patients. Of the patients, 63.5% were oligorecurrent and 36.5% were oligoprogressive. The two-year LPFS was 90.7%. The median NPFS was 22.3 months and the 2-year NPFS was 46.5%. Receiving hormone therapy before SBRT was correlated with a lower NPFS at the multivariate analysis (1 y NPFS 87.1% versus 42.8%; p = 0.002-HR 0.199, 95% CI 0.073-0.549). Furthermore, the oligorecurrent state during ADT was correlated with a lower NPFS than was the oligoprogressive state. The median PFS was 10.3 months and the 2-year PFS was 32.4%. Patients treated with hormone therapy before SBRT had a significantly lower 1-year PFS the others (28% versus 70.4%; p = 0.01-HR 0.259, 95% CI 0.117-0.574). No acute and late toxicities occurred during treatment. Conclusions: The present study is the largest prospective study of 1.5 T lymph node SBRT on MR-linac in patients with PCa. Lymph node SBRT by 1.5 T MR-linac provides high local control rates with an excellent toxicity profile.

Predicting necessity of daily online adaptive replanning based on wavelet image features for MRI guided adaptive radiation therapy

PURPOSE

Online adaptive replanning (OLAR) is generally labor-intensive and time-consuming during MRI-guided adaptive radiation therapy (MRgART). This work aims to develop a method to determine OLAR necessity during MRgART.

METHODS

A machine learning classifier was developed to predict OLAR necessity based on wavelet multiscale texture features extracted from daily MRIs and was trained and tested with data from 119 daily MRI datasets acquired during MRgART for 24 pancreatic cancer patients treated on a 1.5 T MR-Linac. Spearman correlations, interclass correlation (ICC), coefficient of variance (COV), t-test (p < 0.05), self-organized map (SOM) and maximum stable extremal region (MSER) algorithm were used to determine candidate features, which were used to build the prediction models using Bayesian classifiers. The model performance was judged using the AUC of the ROC curve.

RESULTS

Spearman correlation identified 123 features that were not redundant (r < 0.9). Of them 82 showed high ICC for repositioning > 0.6, 67 had a COV greater than 9% for OLAR. Among the 38 features passed the t-test, 25 passed the SOM and 12 passed the MSER. These final 12 features were used to build the classifier model. The combination of 2-3 features at a time was used to build the classifier models. The best performing model was a 3-feature combination, which can predict OLAR necessity with a CV-AUC of 0.98.

CONCLUSIONS

A machine learning classifier model based on the wavelet features extracted from daily MRI for pancreatic cancer was developed to automatically and objectively determine if OLAR is necessary for a treatment fraction avoiding unnecessary effort during MRgART.

MRI-GUIDED RADIOTHERAPY FOR PROSTATE CANCER: A NEW PARADIGM

Radiotherapy is one of the key treatment modalities for primary prostate cancer. During the last decade, significant advances were made in radiotherapy technology leading to increasing both physical and biological precision. Being a loco-regional treatment approach, radiotherapy requires accurate target dose deposition while sparing surrounding healthy tissue. Conventional radiotherapy is based on computerized tomography (CT) images both for radiotherapy planning and image-guidance, however, shortcomings of CT as soft tissue imaging tool are well known. Nowadays, our ability to further escalate radiotherapy dose using hypofractionation is limited by uncertainties in CT-based image guidance and verification. Magnetic resonance imaging (MRI) is a well established imaging method for pelvic organs. In prostate cancer specifically, MRI accurately depicts prostate zonal anatomy, rectum, bladder, and pelvic floor structures with previously unseen precision owing to its sharp soft tissue contrast. The advantages of including MRI in the clinical workflow of prostate cancer radiotherapy are multifold. MRI allows for true adaptive radiotherapy to unfold based on daily MRI images taken before, during and after each radiotherapy fraction. It enables accurate dose escalation to the prostate and intraprostatic tumor lesions. Technically, MRI high-strength magnetic field and linear accelerator high energy electromagnetic beams are hardly compatible, and important efforts were made to overcome these technical challenges and integrate MRI and linear accelerator into one single treatment device, called MRI-linac. Different systems are produced by two leading vendors in the field and currently, there are around 100 MRI-linacs worldwide in clinical operations. In this narrative review paper, we discuss historical perspective of image guidance in radiotherapy, basic elements of MRI, current clinical developments in MRI-guided prostate cancer radiotherapy, and challenges associated with the use of MRI-linac in clinical practice.

First experimental demonstration of VMAT combined with MLC tracking for single and multi fraction lung SBRT on an MR-linac

BACKGROUND AND PURPOSE

VMAT is not currently available on MR-linacs but could maximize plan conformality. To mitigate respiration without compromising delivery efficiency, MRI-guided MLC tumour tracking was recently developed for the 1.5 T Unity MR-linac (Elekta AB, Stockholm, Sweden) in combination with IMRT. Here, we provide a first experimental demonstration of VMAT+MLC tracking for several lung SBRT indications.

MATERIALS AND METHODS

We created central patient and phantom VMAT plans (8×7.5 Gy, 2 arcs) and we created peripheral phantom plans (3×18 & 1×34 Gy, 4 arcs). A motion phantom mimicked subject-recorded respiratory motion (A‾=11 mm, f‾=0.33 Hz, drift‾=0.3 mm/min). This was monitored using 2D-cine MRI at 4 Hz to continuously realign the beam with the target. VMAT+MLC tracking performance was evaluated using 2D film dosimetry and a novel motion-encoded and time-resolved pseudo-3D dosimetry approach.

RESULTS

We found an MLC leaf and jaw end-to-end latency of 328.05(±3.78) ms and 317.33(±4.64) ms, which was mitigated by a predictor. The VMAT plans required maximum MLC speeds of 12.1 cm/s and MLC tracking superimposes an additional 1.48 cm/s. A local 2%/1 mm gamma analysis with a static measurement as reference, revealed pass-rates of 28-46% without MLC tracking and 88-100% with MLC tracking for the 2D film analysis. Similarly the pseudo-3D gamma passing-rates increased from 22-77% to 92-100%. The dose area histograms show that MLC tracking increased the GTV D_98% by 5-20% and the PTV D_95% by 7-24%, giving similar target coverage as their respective static reference.

CONCLUSION

MRI-guided VMAT+MLC tracking is technically feasible on the MR-linac and results in highly conformal dose distribution.