Comparison of Prospectively Generated Glioma Treatment Plans Clinically Delivered on Magnetic Resonance Imaging (MRI)-Linear Accelerator (MR-Linac) Versus Conventional Linac: Predicted and Measured Skin Dose.

Introduction: Magnetic resonance imaging-linear accelerator radiotherapy is an innovative technology that requires special consideration for secondary electron interactions within the magnetic field, which can alter dose deposition at air-tissue interfaces. As part of ongoing quality assurance and quality improvement of new radiotherapy technologies, the purpose of this study was to evaluate skin dose modelled from the treatment planning systems of a magnetic resonance imaging-linear accelerator and a conventional linear accelerator, and then correlate with in vivo measurements of delivered skin dose from each linear accelerator.
Methods: In this prospective cohort study, 37 consecutive glioma patients had treatment planning completed and approved prior to radiotherapy initiation using commercial treatment planning systems: a Monte Carlo-based algorithm for magnetic resonance imaging-linear accelerator or a convolution-based algorithm for conventional linear accelerator. In vivo skin dose was measured using an optically stimulated luminescent dosimeter.
Results: Monte Carlo-based magnetic resonance imaging-linear accelerator plans and convolution-based conventional linear accelerator plans had similar dosimetric parameters for target volumes and organs-at-risk. However, magnetic resonance imaging-linear accelerator plans had 1.52 Gy higher mean dose to air cavities (P < .0001) and 1.10 Gy higher mean dose to skin (P < .0001). In vivo skin dose was 14.5% greater for magnetic resonance imaging-linear accelerator treatments (P = .0027), and was more accurately predicted by Monte Carlo-based calculation (ρ = 0.95, P < .0001) versus convolution-based (ρ = 0.80, P = .0096).
Conclusion: This is the first prospective dosimetric comparison of glioma patients clinically treated on both magnetic resonance imaging-linear accelerator and conventional linear accelerator. Our findings suggest that skin doses were significantly greater with magnetic resonance imaging-linear accelerator plans but correlated better with in vivo measurements of actual skin dose from delivered treatments. Future magnetic resonance imaging-linear accelerator planning processes are being designed to account for skin dosimetry and treatment delivery.

Introduction

Modern radiotherapy with selective use of magnetic resonance imaging (MRI) during the treatment course has shown promise for evaluating glioma tumor dynamics, predicting glioma treatment response and identifying regions for dose escalation based on functional imaging. Typically, no more than 2 to 3 MR images are obtained during a 6-week radiation treatment course due to practicalities and costs of repeat MRIs. The advent of integrated MRI-linear accelerator (MR-Linac) systems enables daily acquisition of high field strength (1.5 Tesla) diagnostic quality MRI to facilitate adaptive radiotherapy, to individualize radiation treatment plans and to incorporate functional imaging. At our institution, we have been developing the application of the Unity MR-Linac (Elekta AB) to brain tumors, and have treated over 150 glioma patients using either 3 or 6 weeks of daily fractionated standard radiotherapy with concurrent temozolomide.

The challenge of delivering radiation within a high MRI field strength environment is the influence of the Lorentz force on dose deposition. The electron return effect (ERE) refers to secondary electrons exiting tissue into air being curved back to deposit dose at the tissue surface due to the Lorentz effect, which results from the presence of the strong magnetic field placed perpendicularly to the beam direction.[3-11] Until now, ERE has been investigated in phantom studies, and in simulated dosimetric studies[13-17] with research versions of MR-Linac planning software prior to the release of the current clinical treatment planning version. Assessment of skin dose is important in radiotherapy for glioma patients because of the potential for dose-related clinical toxicities, and as a component of ongoing quality assurance and quality improvement of new radiotherapy technologies.

In this prospective dosimetric study, our objective was to investigate glioma patients treated on the prospective MOMENTUM trial who received at least one fraction of radiotherapy on both the MR-Linac and a conventional Linac (Elekta AB): (1) the dosimetric impact of the MR-Linac's magnetic field on the target, organs-at-risk (OARs), skin and tissues surrounding air cavities, and; (2) to compare skin dose modelled from the treatment planning software with in vivo skin dose measurements.

Patients and Methods

This prospective cohort study has been approved by the research ethics board of Sunnybrook Health Sciences Centre (Toronto, Canada) until March 3, 2023 (Study Number: 3039). The design and reporting of this study conforms to established STROBE guidelines, and all patient details were de-identified for analysis and presentation in this manuscript. Consecutive glioma patients starting radiotherapy between July 2019 and February 2021 provided informed written consent prior to enrollment on the prospective MOMENTUM trial. All patients completed maximal safe surgical resection followed by either 3 or 6 weeks of adjuvant radiotherapy according to standard practices.[20, 21] Patients were treated on the Unity MR-Linac (Elekta AB) and our standard practice was to create a clinical backup plan on a conventional cone beam CT (CBCT)-guided Linac equipped with the Agility head (Elekta AB) during treatment planning. On days when the MR-Linac was scheduled for maintenance or had a service issue such that treatment would otherwise be delayed, patients were irradiated on the conventional Linac for that particular radiotherapy fraction. All MR-Linac treatment plans and conventional Linac clinical backup plans were generated as per standard of care and deemed clinically acceptable prior to initiation.

All patients were simulated with a treatment planning CT (Brilliance, Philips Healthcare) with a slice thickness of 1.0 mm, and immobilized with a thermoplastic head immobilization mask (CIVCO Medical Solutions). MRI simulation at 1.5 Tesla (Ingenia, Philips Healthcare) was also performed in the treatment position with the immobilization mask applied. Standard volumetric axial T1 gadolinium-enhanced and T2 FLAIR sequences with a slice thickness of 1.0 mm were acquired. CT images and MRI sequences were registered and fused based on rigid mutual information registration within a region of interest box defined around the tumor.

Treatment planning for the MR-Linac was completed on the Monaco treatment planning system (TPS) (Monaco v5.40, Elekta AB). Delineation of OARs, gross tumor volume (GTV), and clinical target volume (CTV) were defined based on glioma consensus contouring guidelines. Planning target volume (PTV) was generated as a 4-mm isotropic expansion of the CTV. The Monte Carlo algorithm accounts for magnetic field effects for the MR-Linac. A TPS based on a convolution-superposition dose calculation algorithm (Pinnacle v9.8, Philips Healthcare) was used for planning on the conventional Linac (Figure 1). All efforts to accurately model skin surface dose using Monaco and Pinnacle were independently benchmarked, and beam models were all commissioned as per TG-53 and TRS-430 guidelines, which recommended dose and spatial tolerances for percentage depth dose and penumbra for a range of field sizes within the beam, outside of the beam, in the penumbra, and at the surface. In particular, Monaco's Monte Carlo dose calculation algorithm has superior agreement when compared to Pinnacle's convolution-superposition dose calculation algorithm in quantifying measured dose compared to TPS dose in the buildup region.

Figure 1.

Clinically delivered treatment plans for a representative study patient diagnosed with glioblastoma, and prescribed 60 Gy in 30 fractions. Top row depicts MR-Linac axial (A), coronal (B), and sagittal (C) dose distributions from Monaco. Bottom row depicts conventional Linac axial (D), coronal (E), and sagittal (F) dose distributions generated from Pinnacle. PTV is depicted in blue colorwash. Multi-color isodose lines are quantified in cGy.

MR-Linac plans used ≥9 coplanar non-opposing step-and-shoot Intensity-Modulated Radiation Therapy (IMRT) 7 MV photon beams, with a maximum dose rate of 425 MU/min, and beam collimation using 80 MLC leaf pairs with a leaf width of 7.175 mm at isocenter. Plans were optimized to achieve at least 90% of the prescription dose covering 99% of PTV (D99%>90%), less than 5% of PTV getting 105% of the prescription dose (V105%<5%), and less than 110% of the prescription dose to 0.1cc volume of PTV (D0.1cc<110%), as outlined by the MR-Linac consortium clinical treatment planning document.

Conventional Linac plans were generated using ≥7 coplanar non-opposing step-and-shoot IMRT 6 MV photon beams, with an additional non-coplanar beam if additional dose optimization was needed, as per standard of care. Maximum dose rate was 500 MU/min, and beam collimation used 80 MLC leaf pairs with a leaf width of 5 mm at isocenter. Plans were optimized to achieve at least 95% of PTV covered by 95% of the prescription dose (V95%≥95%), and at most 1% of PTV getting 105% of the prescription dose (V105%≤1%). All Monte Carlo dose calculations were as per institutional protocol to ensure that the total uncertainty at the dose reference point is less than 0.5%. Dose to OARs were constrained according to our institutional protocols (Supplemental Material, Table 1).

MR-Linac Monaco plans and conventional Linac Pinnacle plans for each patient were evaluated using cumulative dose-volume histogram (DVH) (Figure 2) parameters chosen a priori. For targets, dosimetric parameters of interest were extracted and compared from reference plans for the PTV, CTV, and GTV, including: (1) Volumes of tissue receiving at least 100%, 95%, 50% of the prescription dose (V100%, V95%, V50%); and (2) Minimum doses to 2%, 5%, 50%, 95%, 98% of the target volume (D2%, D5%, D50%, D95%, D98%). Dose homogeneity was assessed using homogeneity index (HI; Supplemental Material, Equation 1), with lower HI values (closer to 0) indicating more homogeneous dose distribution. Dose falloff outside of the target was assessed using gradient index (GI; Supplemental Material, Equation 2), with lower GI values characterizing a sharper dose falloff outside of the target. Dose conformity was assessed using conformity index (CI; Supplemental Material, Equation 3), with CI values closer to 1.0 signifying better conformity. OARs for each patients’ clinical treatment plans were extracted and compared including: Maximum dose to 0.03cc volume of Lens, Globe, Optic chiasm, Optic nerve (Lens D0.03cc, Globe D0.03cc, Optic chiasm D0.03cc, Optic nerve D0.03cc), Maximum dose to 0.01cc volume of Brainstem (Brainstem D0.01cc), and Mean dose to each cochlea (Cochlea Dmean). The maximum dose to a small but finite volume of an OAR was defined in order to partially manage the statistical noise inherent in a Monte Carlo TPS.

Figure 2.

Dose-volume histograms of planning target volume and organs-at-risk (A), as well as tissues surrounding air cavities and at skin surface (B) for a representative study patient. Delivered MR-Linac treatment (solid lines) is compared with conventional Linac treatment (dashed lines).

For our study, skin volumes were generated as a 5 mm rim of tissue contracted from the patient body surface (Supplemental Material, Figure 1). Tissues around air cavities were generated as a 5 mm rim of tissue expanded from air cavity volumes (eg, nasopharynx, paranasal sinuses). Dosimetric parameters analyzed included: Mean dose and Maximum dose to 2cc volumes of tissues surrounding air cavities and skin contours (Air cavity Dmean, Air cavity D2cc, Skin Dmean, Skin D2cc), as well as the volume of skin receiving 20 Gy (Skin V20Gy).

For 10 randomly selected patients, in vivo skin dose measurements during one fraction of MR-Linac treatment and one fraction of conventional Linac treatment were obtained for the same patient using Optically Stimulated Luminescent Dosimeters (OSLDs; nanoDots, Landauer) placed in a defined location on the patient's skin near the PTV (Supplemental Material, Figure 2). Our institution has done extensive OSLD surface dose measurements previously to characterize and benchmark Monaco's Monte Carlo dose calculation engine and has cross-validated OSLD measurements with other methods.[9, 12, 25] Ideal OSLD positions were selected in the TPS at locations near the PTV, where the thermoplastic mask touches the skin. The TPS coordinates of the OSLD relative to the treatment isocenter was noted. The thermoplastic mask was setup on the treatment unit, and with the aid of a 3-dimensional laser positioning system, the OSLD locations were identified on the mask and the OSLD taped to the inside. OSLD dose point in the TPS was defined as the placement location on the treatment unit, at 0.6 mm depth from the patient surface. In the TPS, the patient surface was defined as the 0.6 g/cc boundary between tissue and air. A point dose was selected at a depth of 0.6 mm, with the mean dose within a radius of 0.1 mm about the point reported as the calculated dose. The 0.6 mm depth is the water equivalent depth of an OSLD taking into account its plastic casing, and the accuracy was rated to be ± 3%. Dose delivered to skin, as measured by the OSLD, were correlated to both corresponding clinical treatment plans to determine which TPS best predicted the in vivo dose measurement.

Descriptive statistics were used to summarize the dosimetric parameters. Assumption of normality was assessed using the Smirnov–Kolmogorov test. Student's t-paired test or Wilcoxon signed-rank test, as appropriate based on normality assumption, were used to compare dosimetric parameters between MR-Linac plans and conventional Linac plans. Spearman's correlation was used to assess the relationship between in vivo OSLD measurements and TPS skin dose. All statistical tests were 2-sided, and threshold used for statistical significance was P < .05. Statistical analyses were performed using version 9.4 of the SAS system for Windows (2002-2012 SAS Institute, Inc.).

Results

Patient Population

Thirty-seven glioma patients treated with at least one fraction on both the MR-Linac and conventional Linac met study criteria and were analyzed. Baseline patient characteristics of patients included in the analyses are listed in Table 1.

Table 1.

Patient Characteristics.

Characteristic	Total (N = 37)
Age
<60 years	22 (59.5%)
60-69 years	7 (18.9%)
≥70 years	8 (21.6%)
Gender
Male	25 (67.6%)
Female	12 (32.4%)
WHO grade
II	9 (24.3%)
III	4 (10.8%)
IV	24 (64.9%)
Prescription
40 Gy in 15 fractions	8 (21.6%)
54 Gy in 30 fractions	15 (40.5%)
60 Gy in 30 fractions	13 (35.1%)
59.4 Gy in 33 fractions	1 (2.7%)

Target Volumes and Organs-at-Risk

Comparisons of PTV dosimetric parameters between MR-Linac and conventional Linac plans are detailed in Table 2. There is no statistical difference in PTV V100%, V95%, V50%, D98%, and D95% (P > .05). Small but statistically significant differences were observed in mean PTV D50% (P < .0001), D5% (P < .0001), and D2% (P < .0001) for MR-Linac plans compared with conventional Linac. MR-Linac dose distributions had lower HI (P = .0007), better CI (P = .0064), and equivalent GI (P > .05) compared to conventional Linac.

Table 2.

Dosimetric Impact of Magnetic Field on Target Volume Parameters in Treated Glioma Patients.

PTV	MR-Linac (Mean ± ME)	Conventional Linac (Mean ± ME)	Paired difference (Mean ± ME)	P-value
V100% (cc)	265.94 ± 41.15	265.28 ± 41.43	0.66 ± 3.87	.9705
V95% (cc)	278.48 ± 42.19	269.44 ± 43.00	9.04 ± 14.54	.2920
V50% (cc)	280.07 ± 42.19	279.96 ± 42.17	0.11 ± 0.18	.5000
D98% (Gy)	52.19 ± 2.46	51.57 ± 2.43	0.62 ± 0.66	.5786
D95% (Gy)	53.11 ± 2.50	52.67 ± 2.47	0.44 ± 0.46	.1090
D50% (Gy)	54.67 ± 2.57	54.39 ± 2.57	0.27 ± 0.10	<.0001
D5% (Gy)	55.92 ± 2.63	55.17 ± 2.60	0.75 ± 0.12	<.0001
D2% (Gy)	56.28 ± 2.65	55.42 ± 2.61	0.86 ± 0.12	<.0001
Homogeneity index	1.05 ± 0.01	1.05 ± 0.02	0.01 ± 0.01	.0007
Gradient index	1.07 ± 0.04	1.22 ± 0.33	−0.14 ± 0.29	.8303
Conformity index	0.87 ± 0.02	0.83 ± 0.03	0.04 ± 0.03	.0064

ME = margin of error for 95% CI; Paired difference = MR-Linac − conventional Linac.

Comparison of OAR dosimetric parameters between MR-Linac and conventional Linac plans is shown in Figure 2a, and summarized in Table 3. There was no statistical difference in the D0.03cc for the optic chiasm and optic nerves (P > .05), and no statistical difference in the mean dose to each cochlea (P > .05). Small but statistically significant differences were observed in clinically delivered D0.03cc for brainstem (P = .0307), D0.03cc for each globe (P < .0010), and D0.03cc for each lens (P < .0001) for the MR-Linac.

Table 3.

Dosimetric Impact of Magnetic Field on Organs-at-Risk (OARs) in Treated Glioma Patients.

OARs	MR-Linac (Mean ± ME)	Conventional Linac (Mean ± ME)	Paired difference (Mean ± ME)	P-value
Brainstem D0.1cc (Gy)	47.27 ± 4.48	46.72 ± 4.58	0.55 ± 0.79	.0307
Optic Chiasm D0.03cc (Gy)	43.88 ± 4.78	42.78 ± 5.24	1.11 ± 1.17	.5635
Optic Nerve Right D0.03cc (Gy)	35.47 ± 5.64	35.36 ± 6.09	0.11 ± 1.38	.1804
Optic Nerve Left D0.03cc (Gy)	34.29 ± 5.87	34.64 ± 6.14	−0.35 ± 1.26	.4858
Globe Right D0.03cc (Gy)	24.90 ± 4.51	21.87 ± 4.61	3.03 ± 1.69	.0009
Globe Left D0.03cc (Gy)	22.26 ± 3.54	18.56 ± 3.48	3.70 ± 1.98	.0005
Cochlea Right Dmean (Gy)	18.25 ± 6.70	18.56 ± 6.53	−0.30 ± 1.78	.5597
Cochlea Left Dmean (Gy)	15.52 ± 6.94	16.39 ± 7.05	−0.87 ± 1.80	.5046
Lens Right D0.03cc (Gy)	7.15 ± 0.79	5.25 ± 0.79	1.89 ± 0.67	<.0001
Lens Left D0.03cc (Gy)	7.43 ± 0.75	5.29 ± 0.85	2.14 ± 0.71	<.0001

ME = margin of error for 95% CI; Paired difference = MR-Linac − conventional Linac.

Tissues Around Air Cavities and at Skin Surface

The difference in dosimetric parameters for tissues surrounding air cavities and skin between MR-Linac and conventional Linac plans are shown in Figure 2b, and summarized in Table 4. For tissues surrounding air cavities, MR-Linac plans had statistically significant higher Dmean (P < .0001) and D2cc (P = .0007) compared to conventional Linac. For skin tissues, MR-Linac plans had statistically significant higher Dmean (P < .0001) and V20Gy (P = .0001) compared to conventional Linac, but no difference in D2cc (P = .7975).

Table 4.

Dosimetric Impact of Magnetic Field on Tissues Surrounding Air Cavities and at Skin Surface in Treated Glioma Patients.

Tissues surrounding air cavities and at skin surface	MR-Linac (Mean ± ME)	Conventional Linac (Mean ± ME)	Paired difference (Mean ± ME)	P-value
Air Cavity Dmean (Gy)	12.05 ± 2.27	10.53 ± 2.08	1.52 ± 0.60	<.0001
Air Cavity D2cc (Gy)	41.78 ± 6.10	40.55 ± 6.44	1.23 ± 0.98	.0007
Skin Dmean (Gy)	8.64 ± 1.99	7.54 ± 2.00	1.10 ± 0.54	<.0001
Skin D2cc (Gy)	45.55 ± 2.80	45.45 ± 3.67	0.11 ± 0.85	.7975
Skin V20Gy (cm3)	85.59 ± 13.70	66.55 ± 12.16	19.04 ± 11.45	.0001

ME = margin of error for 95% CI; Paired difference = MR-Linac − conventional Linac.

Differences between predicted skin dose from the TPS and delivered skin dose measured from OSLD are quantified in Table 5. Mean skin dose as determined by the TPS was 7.8 Gy greater on MR-Linac plans (95% CI 3.7-11.8 Gy) compared with conventional Linac (P = .0020). Mean in vivo OSLD skin measurements were 4.6 Gy greater (14.5% higher) on MR-Linac treatments (95% CI 2.1-7.1 Gy) compared with conventional Linac (P = .0027). There was a 1.4 Gy mean difference between MR-Linac Monaco modelled skin dose and in vivo OSLD measured skin dose. In contrast, there was a 4.0 Gy mean difference between conventional Linac Pinnacle modelled skin dose and in vivo OSLD measured skin dose. Furthermore, MR-Linac Monaco predicted in vivo OSLD measured skin dose more accurately (Spearman's correlation ρ = 0.9500, P < .0001). By comparison, there was a weaker association between conventional Linac Pinnacle modelled skin dose and in vivo OSLD measured skin dose (Spearman's correlation ρ = 0.8000, P < .0096).

Table 5.

Comparison of Predicted Skin Dose From TPS Software and Delivered Skin Dose Measured From OSLD.

Skin dose	MR-Linac (Mean ± ME)	Conventional Linac (Mean ± ME)	Mean difference (Mean ± ME)	P-value
TPS (Gy)	35.0 ± 7.7	27.2 ± 6.8	7.8 ± 4.1	.0020
OSLD (Gy)	36.3 ± 9.6	31.7 ± 8.8	4.6 ± 2.5	.0027
Mean difference	1.4 ± 2.1	4.0 ± 5.2
P-value	Spearman ρ = 0.9500, P < .0001	Spearman ρ = 0.8000, P < .0096

ME = margin of error for 95% CI.

Discussion

This is the first prospective dosimetric comparison of glioma patients clinically treated on both MR-Linac and conventional Linac and has several novel components. Our study's patient population uniquely represents a large cohort of glioma patients who received at least one fraction on both MR-Linac and conventional Linac based on clinically approved radiotherapy plans. Furthermore, radiotherapy plans were prospectively generated prior to treatment and delivered on both MR-Linac and conventional Linac, in contrast to previous studies of simulated plans that were retrospectively generated for dosimetric comparison. Lastly and importantly, we performed in vivo skin dose measurements and determined that MR-Linac treatment plans had significantly greater skin doses, but were also better correlated with in vivo skin dose measurements.

We demonstrate that Monaco is able to accurately generate safe MR-Linac radiotherapy treatment plans for glioma patients that achieve planning objectives. We observed that with coplanar beam arrangements, MR-Linac treatments have lower homogeneity, but higher dose conformity and equivalent dose falloff outside of the target, when compared with conventional Linac. Conventional Linac with more coplanar beams, or additional non-coplanar beams could potentially achieve better dose falloff, but at the expense of greater monitor units and longer treatment duration. MR-Linac treatment plans had more heterogeneous dose distributions, which is consistent with the observed small but statistically significant increase in PTV D50%, D5%, and D2%. This is also consistent with previous reports showing higher heterogeneity and higher median V100% for MR-Linac plans compared with conventional Linac. Similarly, a very small but statistically significant increase in Brainstem D0.1cc, each Globe D0.03cc, and each Lens D0.03cc was observed in MR-Linac plans. Over the course of a patient's entire radiotherapy regimen, the absolute summed dose difference was <1 Gy for PTV parameters, <1 Gy for brainstem, approximately 3 Gy for each globe, and approximately 2 Gy for each lens. Since MR-Linac treatment plans are adapted to position relative to intracranial soft tissues every fraction, the exact location of these minimally higher dose regions varies geospatially every fraction, which may negate their effects when accumulated over the treatment course. All MR-Linac treatments achieved standard planning objectives and dose constraints, and it is unlikely that these small differences translate into clinically relevant outcomes.

We also quantitatively characterized the impact of the MR-Linac's magnetic field on delivered dose to skin and tissue surrounding air cavities. Compared to conventional Linac, we observed that MR-Linac treatments had 1.52 Gy higher Dmean (P < .0001), and 1.23 Gy higher D2cc (P = .0007) for tissues surrounding air cavities. Skin D2cc was not statistically different (P > .05), but skin Dmean was 1.10 Gy higher (P < .0001), and skin V20Gy was 19.04 cm3 higher (P = .0001) with MR-Linac treatment. This is consistent with recent preliminary studies investigating the effect of the MR-Linac's magnetic field on radiotherapy treatment.[12-17] Tseng et al used Monaco to retrospectively generate MR-Linac plans with 9 coplanar non-opposing IMRT beams on 24 patients with intact single brain metastases, and found MR-Linac had 0.08 Gy higher Dmean and 0.6 Gy higher D2cc for skin, and 0.07 Gy higher Dmean and 0.3 Gy higher D2cc for tissues around air cavities. Schrenk et al used an open-source Monte Carlo-based TPS to retrospectively generate plans in the presence of a magnetic field with ≥7 coplanar non-opposing 3D-CRT and IMRT beams on 4 patients with non-small cell lung cancer, and found that the presence of the perpendicular magnetic field increased mean dose to tissues surrounding the lung air cavity by 0.5 Gy (18.5%). Nachbar et al used Monaco to prospectively generate a MR-Linac plan with 7 coplanar non-opposing IMRT beams on 1 breast cancer patient, and found MR-Linac had 0.3 Gy higher D2cc and 15.3% higher V35Gy for skin tissue. Xia et al used Monaco to retrospectively generate MR-Linac plans with 9 coplanar non-opposing equally spaced IMRT beams on 10 patients with hypopharyngeal cancer, and found MR-Linac had 1.30-1.81 Gy higher Dmean and 1.68-5.43 Gy higher Dmax for skin, and no difference in dose to tissues surrounding air cavities except for higher Dmax to larynx and trachea. Boldrini et al used MRIdian TPS to retrospectively generate MR-Linac plans with 7 coplanar non-opposing equally-spaced IMRT beams on 10 patients with locally advanced rectal cancer, and found MR-Linac had higher skin dose and higher PTV V105% (14.8%) compared with conventional Linac VMAT with 2 full coplanar arcs (5.0%) and conventional Linac IMRT using 5 coplanar beams (7.3%). Taken together, these studies demonstrated that MR-Linac plans have small increases in dose to skin and tissues surrounding air cavities, and are consistent with our findings.

A potential limitation of our study is the variability in dose fractionations used. However, plan evaluation was performed with pairwise comparisons between MR-Linac and conventional Linac treatments for each patient and are independent of absolute values. Second, there may be potential uncertainty in OSLD measurements caused by placement, air gaps, and surface effects. To mitigate this, a single OSLD measurement was obtained from the same patient's MR-Linac and conventional Linac treatments using standardized technique, although we acknowledge that reproducibility could be assessed by performing additional measurements. Third, we acknowledge that the MR-Linac and conventional Linac utilize different beam energies, different beam parameters, and different TPS to generate standard of care radiotherapy treatment plans that are clinically acceptable. There are differences in how each TPS models patient surface and calculates surface dosimetry. The true effect of the magnetic field on modelled skin dose remains unclear, and highlights the importance of OSLD measurements to properly compare skin surface doses between conventional Linac and MR-Linac treatments. Thus, our objective was to compare the skin dose from standard of care plans on each Linac individually and not an investigation to understand the underlying causes of differences in plan quality. Finally, each patient's clinically delivered treatment plan was analyzed on the latest version of clinical Monaco TPS, and we acknowledge that future iterations of clinical Monaco may have even better dose optimization, dose engines, and dose modelling.

Potential strategies to mitigate skin dose in MR-Linac treatment include increasing the number of beam angles,[25, 29, 30] using opposing beam configurations, using VMAT, using partial arcs, using smaller margins on target volumes, and specifically using a skin structure objective during planning optimization to minimize skin dose. Kim et al investigated the effects of different beam arrangements on skin dosimetry for partial breast radiotherapy using MR-Linac, and demonstrated an 11% to 18% increase in skin D1cc and an 146% to 149% increase in skin V30 with the addition of the magnetic field. However, increasing the number of beam angles, and going from IMRT to VMAT reduced the skin dose in the presence of the magnetic field. Bainbridge et al investigated the effect of a PTV margin reduction from 7 to 3 mm on skin dosimetry for lung cancer radiotherapy using MR-Linac. They demonstrated that MR-Linac plans with 7 mm margins had 0.4 Gy higher Dmean skin dose and 1.7 Gy higher D2cc skin dose compared to conventional Linac plans with 7 mm margins. However, MR-Linac plans with a reduced margin (3 mm) alleviated this issue, while also decreasing other OAR metrics, and allowed isotoxic dose escalation. In glioblastoma, a phase I/II trial used a reduced 5 mm CTV margin with hypofractionated regimens and demonstrated survival outcomes similar to that of historical control. Therefore, MRI-guided treatment of gliomas with a reduced CTV margin and an adaptive framework not only potentially mitigates skin dose effects on the MR-Linac, but may allow for opportunities to decrease toxicities and incorporate boost strategies based on functional imaging. Finally, since this study demonstrates that Monaco's Monte Carlo dose algorithm can accurately model the near-surface dose, using the TPS IMRT optimizer with skin as an avoidance structure can reliably decrease skin dose. Additional work to incorporate skin dose mitigation strategies into planning processes is being developed. Future work to prospectively assess dose-related clinical toxicity outcomes, such as wound dehiscence, permanent alopecia, and radiation dermatitis is important. Further studies investigating clinical scenarios where higher dose is required to skin and superficial targets are warranted as well.

Conclusion

This is the first prospective dosimetric comparison of glioma patients clinically treated with at least one fraction on both MR-Linac and conventional Linac, combined with in vivo correlation to delivered dose on each platform. We observed minimal dosimetric impact of the MR-Linac's magnetic field for target volumes and standard OARs. However, higher doses to tissues at skin surface and surrounding air cavities were observed for clinically delivered MR-Linac treatments. In vivo correlation of delivered skin dose suggests more accurate prediction in MR-Linac plans. Future MR-Linac planning processes are being designed to account for skin dosimetry and treatment delivery. Further work investigating MR-Linac dose-related clinical toxicities, and also clinical scenarios requiring higher dose to skin and superficial targets are important.

Click here for additional data file.

Supplemental material, sj-docx-1-tct-10.1177_15330338221124695 for Comparison of Prospectively Generated Glioma Treatment Plans Clinically Delivered on Magnetic Resonance Imaging (MRI)-Linear Accelerator (MR-Linac) Versus Conventional Linac: Predicted and Measured Skin Dose by Michael H. Wang, Anthony Kim, Mark Ruschin, Hendrick Tan, Hany Soliman, Sten Myrehaug, Jay Detsky, Zain Husain, Eshetu G. Atenafu, Brian Keller, Arjun Sahgal and Chia-Lin Tseng in Technology in Cancer Research & Treatment