Treatment planning
A total of 17 patients (5 for brain, 3 for esophagus, 3 for lungs, 3 for vertebra and 3 for rectum) previously treated with a conventional linac (Elekta Versa HDTM) at our institute were randomly selected from the clinical archive. Among the patients, 11 were men and 6 women, and the median age was 54 years (range from 31 to 74). In a routine clinical procedure, a Computed Tomography (CT) simulation (Siemens SomatomTM, Munich, Germany) was acquired using 140-kVp X rays, 80-cm field of view, and 0.3 cm uniform slice thickness. The target volume and organs at risk (OAR) were delineated by a radiation oncologist with specialty in that treatment site. The clinical plans were designed using Monaco (v5.11). The same planning CT and structure set were used for re-planning on Monaco (v5.40.01). The re-planning for simulated MR-Linac treatment was done by the same dosimetrists following the same institutional protocols.
The impact of magnetic field on dose calculation was taken into account in Monaco (v5.40.01) with the GPUMCD algorithm [14]. The GPUMCD calculation would regress to the results of x-ray voxel Monte Carlo (XVMC) [15] without magnetic field present for the clinical machine for the patients studies. The Unity beam model also accounted for transmission through the cryostat, the couch and receiver coils with a complicated set of attenuation parameters determined at the commissioning. The beam filtration and the nominal beam energy of the MR-Linac also differ from the 6MV photon beams with a Versa HD. Table 1 presents the main parameters of the two Elekta machines (Unity MR-Linac and Versa HD) for comparison.
Table 1 Differences in irradiation geometry between a conventional Versa HD linac and the MR-Linac
Device specifications
|
Versa HD
|
MR-Linac
|
Static magnetic field
|
-
|
1.5 T
|
Nominal beam energy
|
6 MV FFF
|
7 MV FFF
|
Additional beam filtration
|
-
|
Cryostat
|
MLC leaf width at isocenter
|
5.0 mm
|
7.15 mm
|
MLC leaf travel direction
|
Arbitrary
|
Cranio-caudal
|
Source-to-axis distance
|
100 cm
|
143.5 cm
|
Isocenter position relative to patient
|
Variable
|
Fixed at bore center
|
The prescription dose were 60 Gy in 30 fractions for brain plans, 44 Gy in 20 fractions for esophagus plans, 45 Gy in 20 fractions for lung plans, 36 Gy in 10 fractions for vertebra plans and 50 Gy in 25 fractions for rectum plans. All the plans were designed with 5 coplanar fields except the brain cases, which were planned with 9 fields evenly spaced gantry angles. Since currently Elekta Unity only supports the step and shoot (ss) technique for IMRT, while the clinical plans used dynamic MLC (dMLC), for studying the delivery efficiency and TPS usability, reference ss-IMRT plans with Versa HD was generated for all cases. Based on our experience, the difference in the plan quality between dMLC-IMRT and ss-IMRT is negligible if the same beam setup and optimization parameters are used. The step-and-shoot plans for both MR-Linac and Versa HD were all limited to a maximum of 250 segments, while the plans for dMLC were limited to a maximum of 20 control points per beam. Other key parameters for the three groups plans for example nominal energy, dose calculation algorithm, grid spacing, statistical uncertainty, minimum segment area, minimum segment width, minimum MU per segment were presented in table 2.
Table 2 Calculation and segmentation parameters for the three plan groups
Plan parameters
|
MRL-IMRT
|
ss-IMRT
|
dMLC-IMRT
|
Energy
|
7 MV FFF
|
6 MV FFF
|
6 MV FFF
|
Algorithm
|
GPUMCD
|
XVMC
|
XVMC
|
IMRT technique
|
Step-and-shoot
|
Step-and-shoot
|
Dynamic MLC
|
Grid spacing (cm)
|
0.3
|
0.3
|
0.3
|
Statistical uncertainty (%) Per Control Point
|
3
|
3
|
3
|
Minimum segment area (cm2)
|
2
|
2
|
-
|
Minimum segment width (cm)
|
0.5
|
0.5
|
0.5
|
Minimum MU/segment
|
4
|
4
|
-
|
Maximum # segments per plan
|
250
|
250
|
-
|
Maximum # of Control Points per beam
|
-
|
-
|
20
|
Plan evaluation and comparison
Plan evaluations were conducted by the attending physicians and clinical physicists originally assigned to the cases using the Dose-volume histograms (DVH) metrics based on the same institutional protocols. Parameters such as homogeneity index (HI) and conformity index (CI) were used to evaluate the targets dose homogeneity and conformity, while the mean dose (Dmean) and maximum point dose (Dmax) were used for OARs evaluation. The HI and CI are defined as[16,17]:
HI = (D2%-D98%)/D50% (1)
CI = (TVPV×TVPV )/(VPTV×VTV) (2)
where D2%, D98% and D50% represent the minimum dose covering 2%, 98% and 50% of the target volume, respectively. VTV is the treatment volume of the body received the prescribed dose, VPTV is the volume of PTV, and TVPV is the target volume covered by the prescribed dose. The lower the HI value means the better the homogeneity[17]. CI is normally used to quantitatively measure the conformality of the dose distribution relative to the target volume, which denotes the ratio of reference dose received by targets and normal tissue. The CI value closer to 1 means a better conformality[18].
For OARs, the Dmax and Dmean were used to evaluate for a serial organ, while the Dmean and/or Vx (% OARs volume receiving x Gy) were used to evaluate a parallel organ. For the skin dose evaluation, the interested volume is defined as the shell volume with the interior surface as a 4 mm contraction from the body contour. The V10 of the unspecified normal tissue was calculated to evaluate the volume of the low dose region, which defined as body subtracted target volume. For assessing the impact of ERE to the lung, the maximum dose in a 5 mm thick layer inside of the lung surface was evaluated [11]. Table 3 presents the dosimetric parameters evaluated for the treatment sites studied.
Table 3 DVH metrics and evaluation parameters for various treatment sites plans.
Brain
|
Esophagus
|
Lungs
|
Vertebra
|
Rectum
|
PGTV
|
HI
|
PGTV
|
HI
|
PGTV
|
HI
|
PGTV
|
HI
|
PGTV
|
HI
|
CI
|
CI
|
CI
|
CI
|
CI
|
PCTV
|
HI
|
PCTV
|
HI
|
PCTV
|
HI
|
PCTV
|
HI
|
PCTV
|
HI
|
CI
|
CI
|
CI
|
CI
|
CI
|
Brain Stem
|
Dmax
|
Lungs
|
V5
|
Lungs
|
V5
|
Skin
|
Dmax
|
Bladder
|
Dmax
|
Dmean
|
V20
|
V20
|
NT
|
V10
|
Dmean
|
Spinal Cord
|
Dmax
|
V30
|
V30
|
|
|
Small
Bowel
|
Dmax
|
Dmean
|
Dmean
|
Dmean
|
|
|
colon
|
Dmax
|
Optic Chiasm
|
Dmax
|
Spinal Cord
|
Dmax
|
Esophagus
|
Dmax
|
|
|
Femoral head
|
V30
|
Dmean
|
Dmean
|
V40
|
|
|
V40
|
Optic Nerve
|
Dmax
|
Heart
|
V30
|
SpinalCord
|
Dmax
|
|
|
Dmean
|
Dmean
|
V40
|
Dmean
|
|
|
Skin
|
Dmax
|
Eyes
|
Dmean
|
Lung surface
|
Dmax
|
Heart
|
V30
|
|
|
NT
|
V10
|
Lens
|
Dmax
|
Skin
|
Dmax
|
V40
|
|
|
|
|
Skin
|
Dmax
|
NT
|
V10
|
Lung surface
|
Dmax
|
|
|
|
|
NT
|
V10
|
|
|
Skin
|
Dmax
|
|
|
|
|
|
|
|
|
NT
|
V10
|
|
|
|
|
NT: Normal Tissue; Dmax: Maximum dose; Dmean: Mean dose; Vx: % volume receiving x Gy.
The TPS usability is assessed by the time it took for generating an IMRT plan from initiating optimization stage 1 to the end of final dose calculation for each method. Additionally, the total number of MUs for each plan was also recorded to provide an estimate of the deliverability[19].
All analyses were performed using IBM SPSS (v25) statistical software (IBM Corporation, Armonk, NY, USA). Two-sided t tests were carried out and P<0.05 was considered statistically significant. Blinded reviews by the physicians and physicists were performed to assess the plan quality and delivery efficiency.