Precision assessment of bowel motion quantification using 3D cine-MRI for radiotherapy.

Objective. The bowel is an important organ at risk for toxicity during pelvic and abdominal radiotherapy. Identifying regions of high and low bowel motion with MRI during radiotherapy may help to understand the development of bowel toxicity, but the acquisition time of MRI is rather long. The aim of this study is to retrospectively evaluate the precision of bowel motion quantification and to estimate the minimum MRI acquisition time.Approach. We included 22 gynaecologic cancer patients receiving definitive radiotherapy with curative intent. The 10 min pre-treatment 3D cine-MRI scan consisted of 160 dynamics with an acquisition time of 3.7 s per volume. Deformable registration of consecutive images generated 159 deformation vector fields (DVFs). We defined two motion metrics, the 50th percentile vector lengths (VL50) of the complete set of DVFs was used to measure median bowel motion. The 95th percentile vector lengths (VL95) was used to quantify high motion of the bowel. The precision of these metrics was assessed by calculating their variation (interquartile range) in three different time frames, defined as subsets of 40, 80, and 120 consecutive images, corresponding to acquisition times of 2.5, 5.0, and 7.5 min, respectively.Main results. For the full 10 min scan, the minimum motion per frame of 50% of the bowel volume (M50%) ranged from 0.6-3.5 mm for the VL50 motion metric and 2.3-9.0 mm for the VL95 motion metric, across all patients. At 7.5 min scan time, the variation in M50% was less than 0.5 mm in 100% (VL50) and 95% (VL95) of the subsets. A scan time of 5.0 and 2.5 min achieved a variation within 0.5 mm in 95.2%/81% and 85.7%/57.1% of the subsets, respectively.Significance. Our 3D cine-MRI technique quantifies bowel loop motion with 95%-100% confidence with a precision of 0.5 mm variation or less, using a 7.5 min scan time.

Quantifying the reduction of respiratory motion by mechanical ventilation with MRI for radiotherapy.

BACKGROUND: Due to respiratory motion, accurate radiotherapy delivery to thoracic and abdominal tumors is challenging. We aimed to quantify the ability of mechanical ventilation to reduce respiratory motion, by measuring diaphragm motion magnitudes in the same volunteers during free breathing (FB), mechanically regularized breathing (RB) at 22 breaths per minute (brpm), variation in mean diaphragm position across multiple deep inspiration breath-holds (DIBH) and diaphragm drift during single prolonged breath-holds (PBH) in two MRI sessions. METHODS: In two sessions, MRIs were acquired from fifteen healthy volunteers who were trained to be mechanically ventilated non-invasively We measured diaphragm motion amplitudes during FB and RB, the inter-quartile range (IQR) of the variation in average diaphragm position from one measurement over five consecutive DIBHs, and diaphragm cranial drift velocities during single PBHs from inhalation (PIBH) and exhalation (PEBH) breath-holds. RESULTS: RB significantly reduced the respiratory motion amplitude by 39%, from median (range) 20.9 (10.6-41.9) mm during FB to 12.8 (6.2-23.8) mm. The median IQR for variation in average diaphragm position over multiple DIBHs was 4.2 (1.0-23.6) mm. During single PIBHs with a median duration of 7.1 (2.0-11.1) minutes, the median diaphragm cranial drift velocity was 3.0 (0.4-6.5) mm/minute. For PEBH, the median duration was 5.8 (1.8-10.2) minutes with 4.4 (1.8-15.1) mm/minute diaphragm drift velocity. CONCLUSIONS: Regularized breathing at a frequency of 22 brpm resulted in significantly smaller diaphragm motion amplitudes compared to free breathing. This would enable smaller treatment volumes in radiotherapy. Furthermore, prolonged breath-holding from inhalation and exhalation with median durations of six to seven minutes are feasible. TRIAL REGISTRATION: Medical Ethics Committee protocol NL.64693.018.18.

Accuracy of an MR-only workflow for prostate radiotherapy using semi-automatically burned-in fiducial markers.

BACKGROUND: The benefit of MR-only workflow compared to current CT-based workflow for prostate radiotherapy is reduction of systematic errors in the radiotherapy chain by 2-3 mm. Nowadays, MRI is used for target delineation while CT is needed for position verification. In MR-only workflows, MRI based synthetic CT (sCT) replaces CT. Intraprostatic fiducial markers (FMs) are used as a surrogate for the position of the prostate improving targeting. However, FMs are not visible on sCT. Therefore, a semi-automatic method for burning-in FMs on sCT was developed. Accuracy of MR-only workflow using semi-automatically burned-in FMs was assessed and compared to CT/MR workflow. METHODS: Thirty-one prostate cancer patients receiving radiotherapy, underwent an additional MR sequence (mDIXON) to create an sCT for MR-only workflow simulation. Three sources of accuracy in the CT/MR- and MR-only workflow were investigated. To compare image registrations for target delineation, the inter-observer error (IOE) of FM-based CT-to-MR image registrations and soft-tissue-based MR-to-MR image registrations were determined on twenty patients. Secondly, the inter-observer variation of the resulting FM positions was determined on twenty patients. Thirdly, on 26 patients CBCTs were retrospectively registered on sCT with burned-in FMs and compared to CT-CBCT registrations. RESULTS: Image registration for target delineation shows a three times smaller IOE for MR-only workflow compared to CT/MR workflow. All observers agreed in correctly identifying all FMs for 18 out of 20 patients (90%). The IOE in CC direction of the center of mass (COM) position of the markers was within the CT slice thickness (2.5 mm), the IOE in AP and RL direction were below 1.0 mm and 1.5 mm, respectively. Registrations for IGRT position verification in MR-only workflow compared to CT/MR workflow were equivalent in RL-, CC- and AP-direction, except for a significant difference for random error in rotation. CONCLUSIONS: MR-only workflow using sCT with burned-in FMs is an improvement compared to the current CT/MR workflow, with a three times smaller inter observer error in CT-MR registration and comparable CBCT registration results between CT and sCT reference scans. Trial registry Medical Research Involving Human Subjects Act (WMO) does apply to this study and was approved by the Medical Ethics review Committee of the Academic Medical Center. Registration number: NL65414.018.18. Date of registration: 21-08-2018.

Dosimetric comparison of library of plans and online MRI-guided radiotherapy of cervical cancer in the presence of intrafraction anatomical changes.

BACKGROUND: Online magnetic resonance imaging (MRI)-guided radiotherapy of cervical cancer has the potential to further reduce dose to organs at risk (OAR) as compared to a library of plans (LOP) approach. This study presents a dosimetric comparison of an MRI-guided strategy with a LOP strategy taking intrafraction anatomical changes into account. METHODS: The 14 patients included in this study were treated with chemo radiation at our institute and received weekly MRIs after informed consent. The MRI-guided strategy consisted of treatment plans created on the weekly sagittal MRI with 3 mm and 5 mm planning target volume (PTV) margin for clinical target volume (CTV) cervix-uterus (MRI_3mm and MRI_5mm). The plans for the LOP strategy were based on interpolations of CTV cervix-uterus on pretreatment full and empty bladder scans. Dose volume histogram (DVH) parameters were compared for targets and OARs as delineated on the weekly transversal MRI, which was acquired on average 10 min after the sagittal MRI. RESULTS: For the MRI_5mm strategy D98% of the high-risk CTV was at least 95% for all weekly MRIs of all patients, while for the LOP and MRI_3mm strategy this requirement was not satisfied for at least one weekly MRI for 1 and 3 patients, respectively. The average reduction of the volume of the reference dose (95% of the prescribed dose) as compared to the LOP strategy was 464 cm3 for the MRI_3mm strategy, and 422 cm3 for the MRI_5mm strategy. The bowel bag constraint V40Gy < 350 cm3 was violated for 13 patients for the LOP strategy and for 5 patients for both MRI_3mm and MRI_5mm strategy. CONCLUSIONS: With online MRI-guided radiotherapy of cervical cancer considerable sparing of OARs can be achieved. If a new treatment plan can be generated and delivered within 10 min, an online MRI-guided strategy with a 5 mm PTV margin for CTV cervix-uterus is sufficient to account for intrafraction anatomical changes. TRIAL REGISTRATION: NL44492.018.13.

A novel amplitude binning strategy to handle irregular breathing during 4DMRI acquisition: improved imaging for radiotherapy purposes.

BACKGROUND: For radiotherapy of abdominal cancer, four-dimensional magnetic resonance imaging (4DMRI) is desirable for tumor definition and the assessment of tumor and organ motion. However, irregular breathing gives rise to image artifacts. We developed a outlier rejection strategy resulting in a 4DMRI with reduced image artifacts in the presence of irregular breathing. METHODS: We obtained 2D T2-weighted single-shot turbo spin echo images, with an interleaved 1D navigator acquisition to obtain the respiratory signal during free breathing imaging in 2 patients and 12 healthy volunteers. Prior to binning, upper and lower inclusion thresholds were chosen such that 95% of the acquired images were included, while minimizing the distance between the thresholds (inclusion range (IR)). We compared our strategy (Min95) with three commonly applied strategies: phase binning with all images included (Phase), amplitude binning with all images included (MaxIE), and amplitude binning with the thresholds set as the mean end-inhale and mean end-exhale diaphragm positions (MeanIE). We compared 4DMRI quality based on: Data included (DI); percentage of images remaining after outlier rejection. Reconstruction completeness (RC); percentage of bin-slice combinations containing at least one image after binning. Intra-bin variation (IBV); interquartile range of the diaphragm position within the bin-slice combination, averaged over three central slices and ten respiratory bins. IR. Image smoothness (S); quantified by fitting a parabola to the diaphragm profile in a sagittal plane of the reconstructed 4DMRI. A two-sided Wilcoxon's signed-rank test was used to test for significance in differences between the Min95 strategy and the Phase, MaxIE, and MeanIE strategies. RESULTS: Based on the fourteen subjects, the Min95 binning strategy outperformed the other strategies with a mean RC of 95.5%, mean IBV of 1.6 mm, mean IR of 15.1 mm and a mean S of 0.90. The Phase strategy showed a poor mean IBV of 6.2 mm and the MaxIE strategy showed a poor mean RC of 85.6%, resulting in image artifacts (mean S of 0.76). The MeanIE strategy demonstrated a mean DI of 85.6%. CONCLUSIONS: Our Min95 reconstruction strategy resulted in a 4DMRI with less artifacts and more precise diaphragm position reconstruction compared to the other strategies. TRIAL REGISTRATION: Volunteers: protocol W15_373#16.007; patients: protocol NL47713.018.14.

Beam configuration selection for robust intensity-modulated proton therapy in cervical cancer using Pareto front comparison.

The Pareto front reflects the optimal trade-offs between conflicting objectives and can be used to quantify the effect of different beam configurations on plan robustness and dose-volume histogram parameters. Therefore, our aim was to develop and implement a method to automatically approach the Pareto front in robust intensity-modulated proton therapy (IMPT) planning. Additionally, clinically relevant Pareto fronts based on different beam configurations will be derived and compared to enable beam configuration selection in cervical cancer proton therapy. A method to iteratively approach the Pareto front by automatically generating robustly optimized IMPT plans was developed. To verify plan quality, IMPT plans were evaluated on robustness by simulating range and position errors and recalculating the dose. For five retrospectively selected cervical cancer patients, this method was applied for IMPT plans with three different beam configurations using two, three and four beams. 3D Pareto fronts were optimized on target coverage (CTV D(99%)) and OAR doses (rectum V30Gy; bladder V40Gy). Per patient, proportions of non-approved IMPT plans were determined and differences between patient-specific Pareto fronts were quantified in terms of CTV D(99%), rectum V(30Gy) and bladder V(40Gy) to perform beam configuration selection. Per patient and beam configuration, Pareto fronts were successfully sampled based on 200 IMPT plans of which on average 29% were non-approved plans. In all patients, IMPT plans based on the 2-beam set-up were completely dominated by plans with the 3-beam and 4-beam configuration. Compared to the 3-beam set-up, the 4-beam set-up increased the median CTV D(99%) on average by 0.2 Gy and decreased the median rectum V(30Gy) and median bladder V(40Gy) on average by 3.6% and 1.3%, respectively. This study demonstrates a method to automatically derive Pareto fronts in robust IMPT planning. For all patients, the defined four-beam configuration was found optimal in terms of plan robustness, target coverage and OAR sparing.

Quantification of image distortions on the Utrecht interstitial CT/MR brachytherapy applicator at 3T MRI.

PURPOSE: To quantify distortions on MR images of the Utrecht interstitial CT/MR applicator at a field strength of 3T using an MRI-only method. MATERIALS AND METHODS: An MR-compatible phantom suspending the applicator in water was built and imaged on a Philips Ingenia 3T MRI scanner. A map of the magnetic field (B0) was calculated from multiecho images and used to quantify the field inhomogeneity. The expected displacements of the applicator could be quantified using the measured field inhomogeneity and sequence bandwidth. Additionally, two scans were acquired using opposing readout gradients. These scans were rigidly matched and their displacement was compared with the expected displacements from the B0 map. These same methods were applied in 4 patients. By rigid matching of the scans acquired with opposing readout direction the applicator displacement due to image distortion from B0 inhomogeneity as well as patient movement and organ deformation was determined. RESULTS: According to the B0 map, the displacement on the intrauterine device of the plastic brachytherapy applicator was <0.4 mm for both the phantom and patients. Displacements obtained by the opposing readout method were ≤0.8 mm for each patient with a mean ± SD over the patients of 0.3 ± 0.1 mm. CONCLUSION: The results of our study indicate that the B0 method agrees with the opposing readout method. Displacements caused by magnetic field inhomogeneity on 3T MRI were small compared with displacements due to patient movement and organ deformation.

The dosimetric impact of leaf interdigitation and leaf width on VMAT treatment planning in Pinnacle: comparing Pareto fronts.

To evaluate in an objective way the effect of leaf interdigitation and leaf width on volumetric modulated arc therapy plans in Pinnacle. Three multileaf collimators (MLCs) were modeled: two 10 mm leaf width MLCs, with and without interdigitating leafs, and a 5 mm leaf width MLC with interdigitating leafs. Three rectum patients and three prostate patients were used for the planning study. In order to compare treatment techniques in an objective way, a Pareto front comparison was carried out. 200 plans were generated in an automated way, per patient per MLC model, resulting in a total of 3600 plans. From these plans, Pareto-optimal plans were selected which were evaluated for various dosimetric variables. The capability of leaf interdigitation showed little dosimetric impact on the treatment plans, when comparing the 10 mm leaf width MLC with and without leaf interdigitation. When comparing the 10 mm leaf width MLC with the 5 mm leaf width MLC, both with interdigitating leafs, improvement in plan quality was observed. For both patient groups, the integral dose was reduced by 0.6 J for the thin MLC. For the prostate patients, the mean dose to the anal sphincter was reduced by 1.8 Gy and the conformity of the V(95%) was reduced by 0.02 using the thin MLC. The V(65%) of the rectum was reduced by 0.1% and the dose homogeneity with 1.5%. For rectum patients, the mean dose to the bowel was reduced by 1.4 Gy and the mean dose to the bladder with 0.8 Gy for the thin MLC. The conformity of the V(95%) was equivalent for the 10 and 5 mm leaf width MLCs for the rectum patients. We have objectively compared three types of MLCs in a planning study for prostate and rectum patients by analyzing Pareto-optimal plans which were generated in an automated way. Interdigitation of MLC leafs does not generate better plans using the SmartArc algorithm in Pinnacle. Changing the MLC leaf width from 10 to 5 mm generates better treatment plans although the clinical relevance remains to be proven.