End-to-end empirical validation of dose accumulation in MRI-guided adaptive radiotherapy for prostate cancer using an anthropomorphic deformable pelvis phantom

End-to-end testing for stereotactic radiotherapy including the development of a Multi-Modality phantom

PURPOSE

A new insert for a commercially available end-to-end test phantom was designed and in-house manufactured by 3D printing. Subsequently, the insert was tested for different stereotactic radiation therapy workflows (SRS, SBRT, FSRT, and Multimet) also in comparison to the original insert.

MATERIAL AND METHODS

Workflows contained imaging (MR, CT), treatment planning, positioning, and irradiation. Positioning accuracy was evaluated for non-coplanar x-ray, kV- and MV-CBCT systems, as well as surface guided radiation therapy. Dosimetric accuracy of the irradiation was measured with an ionization chamber at four different linear accelerators including dynamic tumor tracking for SBRT.

RESULTS

CT parameters of the insert were within the specification. For MR images, the new insert allowed quantitative analysis of the MR distortion. Positioning accuracy of the phantom with the new insert using the imaging systems of the different linacs was < 1 mm/degree also for MV-CBCT and a non-coplanar imaging system which caused > 3 mm deviation with the original insert. Deviation of point dose values was <3% for SRS, FSRT, and SBRT for both inserts. For the Multimet plans deviations exceeded 10% because the ionization chamber was not positioned in each metastasis, but in the center of phantom and treatment plan.

CONCLUSION

The in-house manufactured insert performed well in all steps of four stereotactic treatment end-to-end tests. Advantages over the commercially available alternative were seen for quantitative analysis of deformation correction in MR images, applicability for non-coplanar x-ray imaging, and dynamic tumor tracking.

Deformable anthropomorphic pelvis phantom for dose accumulation verification

Objective.The validation of deformable image registration (DIR) for contour propagation is often done using contour-based metrics. Meanwhile, dose accumulation requires evaluation of voxel mapping accuracy, which might not be accurately represented by contour-based metrics. By fabricating a deformable anthropomorphic pelvis phantom, we aim to (1) quantify the voxel mapping accuracy for various deformation scenarios, in high- and low-contrast regions, and (2) identify any correlation between dice similarity coefficient (DSC), a commonly used contour-based metric, and the voxel mapping accuracy for each organ.Approach. Four organs, i.e. pelvic bone, prostate, bladder and rectum (PBR), were 3D printed using PLA and a Polyjet digital material, and assembled. The latter three were implanted with glass bead and CT markers within or on their surfaces. Four deformation scenarios were simulated by varying the bladder and rectum volumes. For each scenario, nine DIRs with different parameters were performed on RayStation v10B. The voxel mapping accuracy was quantified by finding the discrepancy between true and mapped marker positions, termed the target registration error (TRE). Pearson correlation test was done between the DSC and mean TRE for each organ.Main results. For the first time, we fabricated a deformable phantom purely from 3D printing, which successfully reproduced realistic anatomical deformations. Overall, the voxel mapping accuracy dropped with increasing deformation magnitude, but improved when more organs were used to guide the DIR or limit the registration region. DSC was found to be a good indicator of voxel mapping accuracy for prostate and rectum, but a comparatively poorer one for bladder. DSC > 0.85/0.90 was established as the threshold of mean TRE ⩽ 0.3 cm for rectum/prostate. For bladder, extra metrics in addition to DSC should be considered.Significance. This work presented a 3D printed phantom, which enabled quantification of voxel mapping accuracy and evaluation of correlation between DSC and voxel mapping accuracy.

Real-time motion-including dose estimation of simulated multi-leaf collimator-tracked magnetic resonance-guided radiotherapy

BACKGROUND

Real-time dose estimation is a key-prerequisite to enable online intra-fraction treatment adaptation in magnetic resonance (MR)-guided radiotherapy (MRgRT). It is an essential component for the assessment of the dosimetric benefits and risks of online adaptive treatments, such as multi-leaf collimator (MLC)-tracking.

PURPOSE

We present a proof-of-concept for a software workflow for real-time dose estimation of MR-guided adaptive radiotherapy based on real-time data-streams of the linac delivery parameters and target positions.

METHODS

A software workflow, combining our in-house motion management software DynaTrack, a real-time dose calculation engine that connects to a research version of the treatment planning software (TPS) Monaco (v.6.09.00, Elekta AB, Stockholm, Sweden) was developed and evaluated. MR-guided treatment delivery on the Elekta Unity MR-linac was simulated with and without MLC-tracking for three prostate patients, previously treated on the Elekta Unity MR-linac (36.25 Gy/five fractions). Three motion scenarios were used: no motion, regular motion, and erratic prostate motion. Accumulated monitor units (MUs), centre of mass target position and MLC-leaf positions, were forwarded from DynaTrack at a rate of 25 Hz to a Monte Carlo (MC) based dose calculation engine which utilises the research GPUMCD-library (Elekta AB, Stockholm, Sweden). A rigid isocentre shift derived from the selected motion scenarios was applied to a bulk density-assigned session MR-image. The respective electron density used for treatment planning was accessed through the research Monaco TPS. The software workflow including the online dose reconstruction was validated against offline dose reconstructions. Our investigation showed that MC-based real-time dose calculations that account for all linac states (including MUs, MLC positions and target position) were infeasible, hence states were randomly sampled and used for calculation as follows; Once a new linac state was received, a dose calculation with 106 photons was started. Linac states that arrived during the time of the ongoing calculation were put into a queue. After completion of the ongoing calculation, one new linac state was randomly picked from the queue and assigned the MU accumulated from the previous state until the last sample in the queue. The queue was emptied, and the process repeated throughout treatment simulation.

RESULTS

On average 27% (23%-30%) of received samples were used in the real-time calculation, corresponding to a calculation time for one linac state of 148 ms. Median gamma pass rate (2%/3 mm local) was 100.0% (99.9%-100%) within the PTV volume and 99.1% (90.1%-99.4.0%) with a 15% dose cut off. Differences in PTVDmean , CTVDmean , RectumD2% , and BladderD2% (offline-online, % of prescribed dose) were below 0.64%. Beam-by-beam comparisons showed deviations below 0.07 Gy. Repeated simulations resulted in standard deviations below 0.31% and 0.12 Gy for the investigated volume and dose criteria respectively.

CONCLUSIONS

Real-time dose estimation was successfully performed using the developed software workflow for different prostate motion traces with and without MLC-tracking. Negligible dosimetric differences were seen when comparing online and offline reconstructed dose, enabling online intra-fraction treatment decisions based on estimates of the delivered dose.

Impact of daily plan adaptation on accumulated doses in ultra-hypofractionated magnetic resonance-guided radiation therapy of prostate cancer

Background and purpose

Ultra-hypofractionated online adaptive magnetic resonance-guided radiotherapy (MRgRT) is promising for prostate cancer. However, the impact of online adaptation on target coverage and organ-at-risk (OAR) sparing at the level of accumulated dose has not yet been reported. Using deformable image registration (DIR)-based accumulation, we compared the delivered adapted dose with the simulated non-adapted dose.

Materials and methods

Twenty-three prostate cancer patients treated at two clinics with 0.35 T magnetic resonance-guided linear accelerator (MR-linac) following the same treatment protocol (5 ×  7.5 Gy with urethral sparing and daily adaptation) were included. The fraction MR images were deformably registered to the planning MR image. Both non-adapted and adapted fraction doses were accumulated with the corresponding vector fields. Two DIR approaches were implemented. PTV* (planning target volume minus urethra+2mm) D95%, CTV* (clinical target volume minus urethra) D98%, and OARs (urethra+2mm, bladder, and rectum) D0.2cc, were evaluated. Statistical significance was inferred from a two-tailed Wilcoxon signed-rank test (p < 0.05).

Results

Normalized to the baseline, the accumulated PTV* D95% increased significantly by 2.7 % ([1.5, 4.3]%) through adaptation, and the CTV* D98% by 1.2 % ([0.1, 1.7]%). For the OARs after adaptation, accumulated bladder D0.2cc decreased by 0.4 % ([-1.2, 0.4]%), urethra+2mmD0.2cc by 0.8 % ([-1.6, -0.1]%), while rectum D0.2cc increased by 2.6 % ([1.2, 4.9]%). For all patients, rectum D0.2cc was still below the clinical constraint. Results of both DIR approaches differed on average by less than 0.2 %.

Conclusions

Online adaptation in MRgRT improved target coverage and OARs sparing at the level of accumulated dose.

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.

Experimental validation of a novel method of dose accumulation for the rectum

BACKGROUND

Dose-surface maps (DSMs) are an increasingly popular tool to evaluate spatial dose-outcome relationships for the rectum. Recently, DSM addition has been proposed as an alternative method of dose accumulation from deformable registration-based techniques. In this study, we performed the first experimental investigation of the accuracy at which DSM accumulation can capture the total dose delivered to a rectum's surface in the presence of inter-fraction motion.

MATERIAL AND METHODS

A custom PVC rectum phantom capable of representing typical rectum inter-fraction motion and filling variations was constructed for this project. The phantom allowed for the placement of EBT3 film sheets on the representative rectum surface to measure rectum surface dose. A multi-fraction prostate VMAT treatment was designed and delivered to the phantom in a water tank for a variety of inter-fraction motion scenarios. DSMs for each fraction were calculated in two ways using CBCT images acquired during delivery and summed to produce accumulated DSMs. Accumulated DSMs were then compared to film measurements using gamma analysis (3%/2 mm criteria). Similarity of isodose clusters between films and DSMs was also investigated.

RESULTS

Baseline agreement between film measurements and accumulated DSMs for a stationary rectum was 95.6%. Agreement between film and accumulated DSMs in the presence of different types of inter.-fraction motion was ≥92%, and isodose cluster mean distance to agreement was within 1.5 mm for most scenarios. Overall, DSM accumulation performed the best when using DSMs that accounted for changes in rectum path orientation.

CONCLUSION

Dose accumulation performed with DSMs was found to accurately replicate total delivered dose to a rectum phantom in the presence of inter-fraction motion.

Integration of operator-validated contours in deformable image registration for dose accumulation in radiotherapy

Background and Purpose

Deformable image registration (DIR) is a core element of adaptive radiotherapy workflows, integrating daily contour propagation and/or dose accumulation in their design. Propagated contours are usually manually validated and may be edited, thereby locally invalidating the registration result. This means the registration cannot be used for dose accumulation. In this study we proposed and evaluated a novel multi-modal DIR algorithm that incorporated contour information to guide the registration. This integrates operator-validated contours with the estimated deformation vector field and warped dose.

Materials and Methods

The proposed algorithm consisted of both a normalized gradient field-based data-fidelity term on the images and an optical flow data-fidelity term on the contours. The Helmholtz-Hodge decomposition was incorporated to ensure anatomically plausible deformations. The algorithm was validated for same- and cross-contrast Magnetic Resonance (MR) image registrations, Computed Tomography (CT) registrations, and CT-to-MR registrations for different anatomies, all based on challenging clinical situations. The contour-correspondence, anatomical fidelity, registration error, and dose warping error were evaluated.

Results

The proposed contour-guided algorithm considerably and significantly increased contour overlap, decreasing the mean distance to agreement by a factor of 1.3 to 13.7, compared to the best algorithm without contour-guidance. Importantly, the registration error and dose warping error decreased significantly, by a factor of 1.2 to 2.0.

Conclusions

Our contour-guided algorithm ensured that the deformation vector field and warped quantitative information were consistent with the operator-validated contours. This provides a feasible semi-automatic strategy for spatially correct warping of quantitative information even in difficult and artefacted cases.

Comparison of different dose accumulation strategies to estimate organ doses after stereotactic magnetic resonance-guided adaptive radiotherapy

INTRODUCTION

Re-irradiation is frequently performed in the era of precision oncology, but previous doses to organs-at-risk (OAR) must be assessed to avoid cumulative overdoses. Stereotactic magnetic resonance-guided online adaptive radiotherapy (SMART) enables highly precise ablation of tumors close to OAR. However, OAR doses may change considerably during adaptive treatment, which complicates potential re-irradiation. We aimed to compare the baseline plan with different dose accumulation techniques to inform re-irradiation.

PATIENTS & METHODS

We analyzed 18 patients who received SMART to lung or liver tumors inside prospective databases. Cumulative doses were calculated inside the planning target volumes (PTV) and OAR for the adapted plans and theoretical non-adapted plans via (1) cumulative dose volume histograms (DVH sum plan) and (2) deformable image registration (DIR)-based dose accumulation to planning images (DIR sum plan). We compared cumulative dose parameters between the baseline plan, DVH sum plan and DIR sum plan using equivalent doses in 2 Gy fractions (EQD2).

RESULTS

Individual patients presented relevant increases of near-maximum doses inside the proximal bronchial tree, spinal cord, heart and gastrointestinal OAR when comparing adaptive treatment to the baseline plans. The spinal cord near-maximum doses were significantly increased in the liver patients (D2% median: baseline 6.1 Gy, DIR sum 8.1 Gy, DVH sum 8.4 Gy, p = 0.04; D0.1 cm³ median: baseline 6.1 Gy, DIR sum 8.1 Gy, DVH sum 8.5 Gy, p = 0.04). Three OAR overdoses occurred during adaptive treatment (DIR sum: 1, DVH sum: 2), and four more intense OAR overdoses would have occurred during non-adaptive treatment (DIR sum: 4, DVH sum: 3). Adaptive treatment maintained similar PTV coverages to the baseline plans, while non-adaptive treatment yielded significantly worse PTV coverages in the lung (D95% median: baseline 86.4 Gy, DIR sum 82.4 Gy, DVH sum 82.2 Gy, p = 0.006) and liver patients (D95% median: baseline 87.4 Gy, DIR sum 82.1 Gy, DVH sum 81.1 Gy, p = 0.04).

CONCLUSION

OAR doses can increase during SMART, so that re-irradiation should be planned based on dose accumulations of the adapted plans instead of the baseline plan. Cumulative dose volume histograms represent a simple and conservative dose accumulation strategy.

Feasibility of delivered dose reconstruction for MR-guided SBRT of pancreatic tumors with fast, real-time 3D cine MRI

BACKGROUND AND PURPOSE

In MR-guided SBRT of pancreatic cancer, intrafraction motion is typically monitored with (interleaved) 2D cine MRI. However, tumor surroundings are often not fully captured in these images, and motion might be distorted by through-plane movement. In this study, the feasibility of highly accelerated 3D cine MRI to reconstruct the delivered dose during MR-guided SBRT was assessed.

MATERIALS AND METHODS

A 3D cine MRI sequence was developed for fast, time-resolved 4D imaging, featuring a low spatial resolution that allows for rapid volumetric imaging at 430 ms. The 3D cines were acquired during the entire beam-on time of 23 fractions of online adaptive MR-guided SBRT for pancreatic tumors on a 1.5 T MR-Linac. A 3D deformation vector field (DVF) was extracted for every cine dynamic using deformable image registration. Next, these DVFs were used to warp the partial dose delivered in the time interval between consecutive cine acquisitions. The warped dose plans were summed to obtain a total delivered dose. The delivered dose was also calculated under various motion correction strategies. Key DVH parameters of the GTV, duodenum, small bowel and stomach were extracted from the delivered dose and compared to the planned dose. The uncertainty of the calculated DVFs was determined with the inverse consistency error (ICE) in the high-dose regions.

RESULTS

The mean (SD) relative (ratio delivered/planned) D_99% of the GTV was 0.94 (0.06), and the mean (SD) relative D_0.5cc of the duodenum, small bowel, and stomach were respectively 0.98 (0.04), 1.00 (0.07), and 0.98 (0.06). In the fractions with the lowest delivered tumor coverage, it was found that significant lateral drifts had occurred. The DVFs used for dose warping had a low uncertainty with a mean (SD) ICE of 0.65 (0.07) mm.

CONCLUSION

We employed a fast, real-time 3D cine MRI sequence for dose reconstruction in the upper abdomen, and demonstrated that accurate DVFs, acquired directly from these images, can be used for dose warping. The reconstructed delivered dose showed only a modest degradation of tumor coverage, mostly attainable to baseline drifts. This emphasizes the need for motion monitoring and development of intrafraction treatment adaptation solutions, such as baseline drift corrections.

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.

Dosimetric impact of deformable image registration using radiophotoluminescent glass dosimeters with a physical geometric phantom

PURPOSE

To study the dosimetry impact of deformable image registration (DIR) using radiophotoluminescent glass dosimeter (RPLD) and custom developed phantom with various inserts.

METHODS

The phantom was developed to facilitate simultaneous evaluation of geometric and dosimetric accuracy of DIR. Four computed tomography (CT) images of the phantom were acquired with four different configurations. Four volumetric modulated arc therapy (VMAT) plans were computed for different phantom. Two different patterns were applied to combination of four phantom configurations. RPLD dose measurement was combined between corresponding two phantom configurations. DIR-based dose accumulation was calculated between corresponding two CT images with two commercial DIR software and various DIR parameter settings, and an open source software. Accumulated dose calculated using DIR was then compared with measured dose using RPLD.

RESULTS

The mean ± standard deviation (SD) of dose difference was 2.71 ± 0.23% (range, 2.22%-3.01%) for tumor-proxy and 3.74 ± 0.79% (range, 1.56%-4.83%) for rectum-proxy. The mean ± SD of target registration error (TRE) was 1.66 ± 1.36 mm (range, 0.03-4.43 mm) for tumor-proxy and 6.87 ± 5.49 mm (range, 0.54-17.47 mm) for rectum-proxy. These results suggested that DIR accuracy had wide range among DIR parameter setting.

CONCLUSIONS

The dose difference observed in our study was 3% for tumor-proxy and within 5% for rectum-proxy. The custom developed physical phantom with inserts showed potential for accurate evaluation of DIR-based dose accumulation. The prospect of simultaneous evaluation of geometric and dosimetric DIR accuracy in a single phantom may be useful for validation of DIR for clinical use.

Interfractional dose accumulation for MR-guided liver SBRT: Variation among algorithms is highly patient- and fraction-dependent

BACKGROUND AND PURPOSE

Daily plan adaptations could take the dose delivered in previous fractions into account. Due to high dose delivered per fraction, low number of fractions, steep dose gradients, and large interfractional organ deformations, this might be particularly important for liver SBRT. This study investigates inter-algorithm variation of interfractional dose accumulation for MR-guided liver SBRT.

MATERIALS AND METHODS

We assessed 27 consecutive MR-guided liver SBRT treatments of 67.5 Gy in three (n = 15) or 50 Gy in five fractions (n = 12), both prescribed to the GTV. We calculated fraction doses on daily patient anatomy, warped these doses to the simulation MRI using seven different algorithms, and accumulated the warped doses. Thus, we obtained differences in planned doses and warped or accumulated doses for each algorithm. This enabled us to calculate the inter-algorithm variations in warped doses per fraction and in accumulated doses per treatment course.

RESULTS

The four intensity-based algorithms were more consistent with planned PTV dose than affine or contour-based algorithms. The mean (range) variation of the dose difference for PTV D_95% due to dose warping by these intensity-based algorithms was 10.4 percentage points (0.3 to 43.7) between fractions and 8.6 (0.3 to 24.9) between accumulated treatment doses. As seen by these ranges, the variation was very dependent on the patient and the fraction being analyzed. Nevertheless, no correlations between patient or plan characteristics on the one hand and inter-algorithm dose warping variation on the other hand was found.

CONCLUSION

Inter-algorithm dose accumulation variation is highly patient- and fraction-dependent for MR-guided liver SBRT. We advise against trusting a single algorithm for dose accumulation in liver SBRT.

Evaluation of DIR schemes on tumor/organ with progressive shrinkage: accuracy of tumor/organ internal tissue tracking during the radiation treatment

PURPOSE

Accuracy of intratumoral treatment dose accumulation and response assessment highly depends on the accuracy of a DIR method. However, achievable accuracy of the existing DIR methods for tumor/organ with large and progressive shrinkage during the radiotherapy course have not been explored. This study aimed to use a bio-tissue phantom to quantify the achievable accuracy of different DIR schemes.

MATERIALS /METHODS

A fresh porcine liver was used for phantom material. Sixty gold markers were implanted on the surface and inside of the liver. To simulate the progressive radiation-induced tumor/organ shrinkage, the phantom was heated using a microwave oven incrementally from 30s to 200s in 8 phases. For each phase, the phantom was scanned by CT. Two extra image sets were generated from the original images: 1) the image set with overriding the high-density gold markers (feature image); 2) the image set with overriding the entire phantom to the mean soft tissue intensity (featureless image). Ten DIR schemes were evaluated to mimic clinical treatment situations of tumor/critical organ with respect to their surface and internal condition of image features, availability of intermediate feedback images and DIR methods. The internal marker's positions were utilized to evaluate DIR accuracy quantified by target registration error (TRE).

RESULTS

Volume reduction was about 20% ∼ 40% of the initial volume after 90s ∼ 200s of the heating. Without image features on the surface and inside of the phantom, the hybrid-DIR (image-based DIR followed by biomechanical model-based refinement) with the surface constraint achieved the registration TRE from 2.6 ± 1.2mm to 5.3 ± 2.6mm proportional to the %volume shrinkage. Meanwhile, the hybrid-DIR with the surface-marker constraint achieved the TRE from 2.4 ± 1.2mm to 2.6 ± 1.0mm. If both the surface and internal image features would be viable on the feedback images, the achievable accuracy could be minimal with the TRE from 1.6±0.9mm to 1.9 ± 1.2mm.

CONCLUSIONS

Standard DIR methods cannot guarantee intratumoral tissue registration accuracy for tumor/organ with large progressive shrinkage. Achievable accuracy with using the hybrid DIR method is highly dependent on the surface registration accuracy. If the surface registration mean TRE can be controlled within 2mm, the mean TRE of internal tissue can be controlled within 3mm.

MIRSIG position paper: the use of image registration and fusion algorithms in radiotherapy

The report of the American Association of Physicists in Medicine (AAPM) Task Group No. 132 published in 2017 reviewed rigid image registration and deformable image registration (DIR) approaches and solutions to provide recommendations for quality assurance and quality control of clinical image registration and fusion techniques in radiotherapy. However, that report did not include the use of DIR for advanced applications such as dose warping or warping of other matrices of interest. Considering that DIR warping tools are now readily available, discussions were hosted by the Medical Image Registration Special Interest Group (MIRSIG) of the Australasian College of Physical Scientists & Engineers in Medicine in 2018 to form a consensus on best practice guidelines. This position statement authored by MIRSIG endorses the recommendations of the report of AAPM task group 132 and expands on the best practice advice from the ‘Deforming to Best Practice’ MIRSIG publication to provide guidelines on the use of DIR for advanced applications.

Dosimetric considerations for moldable silicone composites used in radiotherapy applications

Due to their many favorable characteristics, moldable silicone (MS) composites have gained popularity in medicine and recently, in radiotherapy applications. We investigate the dosimetric properties of silicones in radiotherapy beams and determine their suitability as water substitutes for constructing boluses and phantoms. Two types of silicones were assessed (ρ= 1.04 g/cm³ and ρ= 1.07 g/cm³). Various dosimetric properties were characterized, including the relative electron density, the relative mean mass energy‐absorption coefficient, and the relative mean mass restricted stopping power. Silicone slabs with thickness of 1.5 cm and 5.0 cm were molded to mimic a bolus setup and a phantom setup, respectively. Measurements were conducted for Co‐60 and 6 MV photon beams, and 6 MeV electron beams. The doses at 1.5 cm and 5.0 cm depths in MS were measured with solid water (SW) backscatter material (D_MS–SW), and with a full MS setup (D_MS–MS), then compared with doses at the same depths in a full SW setup (D_SW–SW). Relative doses were reported as D_MS–SW/D_MS–SW and D_MS–MS/D_SW–SW. Experimental results were verified using Monaco treatment planning system dose calculations and Monte Carlo EGSnrc simulations. Film measurements showed varying dose ratios according to MS and beam types. For photon beams, the bolus setup D_MS–SW/D_SW–SW exhibited a 5% relative dose reduction. The dose for 6 MV beams was reduced by nearly 2% in a full MS setup. Up to 2% dose increase in both scenarios was observed for electron beams. Compared with dose in SW, an interface of MS–SW can cause relatively high differences. We conclude that it is important to characterize a particular silicone's properties in a given beam quality prior to clinical use. Because silicone compositions vary between manufacturers and differ from water/SW, accurate dosimetry using these materials requires consideration of the reported differences.

Experimental validation of daily adaptive proton therapy

Anatomical changes during proton therapy require rapid treatment plan adaption to mitigate the associated dosimetric impact. This in turn requires a highly efficient workflow that minimizes the time between imaging and delivery. At the Paul Scherrer Institute, we have developed an online adaptive workflow, which is specifically designed for treatments in the skull-base/cranium, with the focus set on simplicity and minimizing changes to the conventional workflow. The dosimetric and timing performance of this daily adaptive proton therapy (DAPT) workflow has been experimentally investigated using an in-house developed DAPT software and specifically developed anthropomorphic phantom. After a standard treatment preparation, which includes the generation of a template plan, the treatment can then be adapted each day, based on daily imaging acquired on an in-room CT. The template structures are then rigidly propagated to this CT and the daily plan is fully re-optimized using the same field arrangement, DVH constraints and optimization settings of the template plan. After a dedicated plan QA, the daily plan is delivered. To minimize the time between imaging and delivery, clinically integrated software for efficient execution of all online adaption steps, as well as tools for comprehensive and automated QA checks, have been developed. Film measurements of an end-to-end validation of a multi-fraction DAPT treatment showed high agreement to the calculated doses. Gamma pass rates with a 3%/3 mm criteria were >92% when comparing the measured dose to the template plan. Additionally, a gamma pass rate >99% was found comparing measurements to the Monte Carlo dose of the daily plans reconstructed from the logfile, accumulated over the delivered fractions. With this, we experimentally demonstrate that the described adaptive workflow can be delivered accurately in a timescale similar to a standard delivery.

Technical Note: On the feasibility of performing dosimetry in target and organ at risk using polymer dosimetry gel and thermoluminescence detectors in an anthropomorphic, deformable, and multimodal pelvis phantom

OBJECTIVE

To assess the feasibility of performing dose measurements in the target (prostate) and an adjacent organ at risk (rectum) using polymer dosimetry gel and thermoluminescence detectors (TLDs) in an anthropomorphic, deformable, and multimodal pelvis phantom (ADAM PETer).

METHODS

The 3D printed prostate organ surrogate of the ADAM PETer phantom was filled with polymer dosimetry gel. Nine TLD600 (LiF:Mg,Ti) were installed in 3 × 3 rows on a specifically designed 3D-printed TLD holder. The TLD holder was inserted into the rectum at the level of the prostate and fixed by a partially inflated endorectal balloon. Computed tomography (CT) images were taken and treatment planning was performed. A prescribed dose of 4.5 Gy was delivered to the planning target volume (PTV). The doses measured by the dosimetry gel in the prostate and the TLDs in the rectum ("measured dose") were compared to the doses calculated by the treatment planning system ("planned dose") on a voxel-by-voxel basis.

RESULTS

In the prostate organ surrogate, the 3D-γ-index was 97.7% for the 3% dose difference and 3 mm distance to agreement criterium. In the center of the prostate organ surrogate, measured and planned doses showed only minor deviations (<0.1 Gy, corresponding to a percentage error of 2.22%). On the edges of the prostate, slight differences between planned and measured doses were detected with a maximum deviation of 0.24 Gy, corresponding to 5.3% of the prescribed dose. The difference between planned and measured doses in the TLDs was on average 0.08 Gy (range: 0.02-0.21 Gy), corresponding to 1.78% of the prescribed dose (range: 0.44%-4.67%).

CONCLUSIONS

The present study demonstrates the feasibility of using polymer dosimetry gel and TLDs for 3D and 1D dose measurements in the prostate and the rectum organ surrogates in an anthropomorphic, deformable and multimodal phantom. The described methodology might offer new perspectives for end-to-end tests in image-guided adaptive radiotherapy workflows.

A strategy to prevent a temperature-induced MRI artifact in warm liquid phantoms due to convection currents

MRI phantom studies often fail to mimic the temperature of the human body, which can negatively impact accuracy. An artifact induced by increasing temperature in liquid phantoms was observed, presenting a significant challenge to temperature-controlled experiments. In this study we characterize and provide a solution to eliminate this temperature-induced MRI artifact. Low concentration (0.5-2.5 mM) agar phantoms were prepared. Utilizing a temperature-controlled phantom holder, T₁ - and T₂ -weighted structural images were acquired at 7 T along with quantitative B₀ , B₁ , T₁ , T₂ and ADC maps at both 25 and 37°C. Additionally, computer simulations were conducted to demonstrate the fluid flow and thermal flux patterns in water to provide an insight into the origins of the artifact. Evidence from computer simulation and quantitative MRI strongly suggest the artifact was caused by heat transfer in the form of natural convection leading to structured patterns of signal loss in MR images. The artifact was present up to agar concentrations of 1.5 mM (T₁ = 3068 ± 16 ms, T₂ = 1052 ± 20 ms, ADC = 2.29 ± 0.36 × 10^-3 mm² /s at 25°C; T₁ = 3928 ± 44 ms, T₂ = 1122 ± 24 ms, ADC = 2.64 ± 0.49 × 10^-3 mm² /s at 37°C), above which point increased sample viscosity no longer allows for convection currents, thereby eliminating the artifact. The methodology described in this work simplifies quantitative MR acquisition of liquid phantoms at physiological temperature by suppressing convection currents with relatively small changes to intrinsic MR parameters (T₁ increased by 1.4% and T₂ decreased by 17% for 1.5 mM agar at 25°C).

Dose accumulation for personalized stereotactic MR-guided adaptive radiation therapy in prostate cancer

BACKGROUND AND PURPOSE

Adaptive MR-guided radiotherapy (MRgRT) is an innovative approach for delivering stereotactic body radiotherapy (SBRT) in prostate cancer (PC). Despite the increased clinical use of SBRT for PC, there is limited data on the relation between the actual delivered dose and toxicity. We aimed to identify dose parameters based on the total accumulated delivered bladder dose (DOSE_ACC^TX). Furthermore, for future personalization, we studied whether prospective accumulation of the first 3 of 5 fractions (DOSE_ACC^3FR) could be used as a representative of DOSE_ACC^TX.

MATERIALS AND METHODS

We deployed a recently validated deformable image registration-based dose accumulation strategy to reconstruct DOSE_ACC^TX and DOSE_ACC^3FR in 101 PC patients treated with stereotactic MRgRT. IPSS scores at baseline, end of MRgRT, at 6 and 12 weeks after treatment were analyzed to identify a clinically relevant increase of acute urinary symptoms. A receiver operator characteristic curve analysis was used to investigate the correlation of an increase in IPSS and bladder DOSE_ACC^TX (range V₅-V_{36.25 Gy,} D_1cc, D_5cc) and DOSE_ACC^3FR (range V₆-V_{21.8 Gy,} D_1cc, D_5cc) parameters.

RESULTS

A clinically relevant increase in IPSS in the three months following MRgRT was observed in 25 patients. The V_20Gy-32Gy from DOSE_ACC^TX and V_15Gy-18Gy from DOSE_ACC^3FR showed good correlation with IPSS increase with area under the curve (AUC) values ranging from 0.71 to 0.75. In contrast, baseline dosimetry showed a poor correlation with AUC values between 0.53 and 0.62.

CONCLUSION

DOSE_ACC^TX was superior to baseline dosimetry in predicting acute urinary symptoms. Because DOSE_ACC^3FR also showed good correlation, this can potentially be used to optimize MRgRT for the remaining fractions.

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

OBJECTIVE

To develop an anthropomorphic, deformable and multimodal pelvis phantom with positron emission tomography extension for radiotherapy (ADAM PETer).

METHODS

The design of ADAM PETer was based on our previous pelvis phantom (ADAM) and extended for compatibility with PET and use in 3T magnetic resonance imaging (MRI). The formerly manually manufactured silicon organ surrogates were replaced by three-dimensional (3D) printed organ shells. Two intraprostatic lesions, four iliac lymph node metastases and two pelvic bone metastases were added to simulate prostate cancer as multifocal and metastatic disease. Radiological properties [computed tomography (CT) and 3T MRI] of cortical bone, bone marrow and adipose tissue were simulated by heavy gypsum, a mixture of Vaseline and K₂ HPO₄ and peanut oil, respectively. For soft tissues, agarose gels with varying concentrations of agarose, gadolinium (Gd) and sodium fluoride (NaF) were developed. The agarose gels were doped with patient-specific activity concentrations of a Fluorine-18 labelled compound and then filled into the 3D printed organ shells of prostate lesions, lymph node and bone metastases. The phantom was imaged at a dual energy CT and a 3T PET/MRI scanner.

RESULTS

The compositions of the soft tissue surrogates are the following (given as mass fractions of agarose[w%]/NaF[w%]/Gd[w%]): Muscle (4/1/0.027), prostate (1.35/4.2/0.011), prostate lesions (2.25/4.2/0.0085), lymph node and bone metastases (1.4/4.2/0.025). In all imaging modalities, the phantom simulates human contrast. Intraprostatic lesions appear hypointense as compared to the surrounding normal prostate tissue in T2-weighted MRI. The PET signal of all tumors can be localized as focal spots at their respective site. Activity concentrations of 12.0 kBq/mL (prostate lesion), 12.4 kBq/mL (lymph nodes) and 39.5 kBq/mL (bone metastases) were measured.

CONCLUSION

The ADAM PETer pelvis phantom can be used as multimodal, anthropomorphic model for CT, 3T-MRI and PET measurements. It will be central to simulate and optimize the technical workflow for the integration of PET/MRI-based radiation treatment planning of prostate cancer patients.

Dosimetric impact of intrafraction prostate rotation and accuracy of gating, multi-leaf collimator tracking and couch tracking to manage rotation: An end-to-end validation using volumetric film measurements

BACKGROUND AND PURPOSE

Both gating and tracking can mitigate the deteriorating dosimetric impact of intrafraction translation during prostate stereotactic body radiotherapy (SBRT). However, their ability to manage intrafraction rotation has not yet been thoroughly investigated. The dosimetric accuracy of gating, MLC tracking and couch tracking to manage intrafraction prostate rotation was investigated.

MATERIALS AND METHODS

Treatment plans for end-to-end tests of prostate SBRT with focal boosting were generated for a dynamic anthropomorphic pelvis phantom. The phantom applied internal lateral rotation (up to 25°) and coupled vertical and longitudinal translation of a radiochromic film stack that was used for dose measurements. Dose was delivered for each plan while the phantom applied motion according to three typical prostate motion traces without compensation (i), with gating (ii), with MLC tracking (iii) or with couch tracking (iv). Measured doses for the four motion compensation strategies were compared with the planned dose in terms of γ-index analysis, target coverage and organs at risk (OAR) sparing.

RESULTS

Intrafraction rotation reduced the 3%(global)/2mm γ-index passing rate (γPR) for the prostate target volume by median (range) -33.2% (-68.6%, -4.1%) when no motion compensation was applied. The use of motion compensation improved the γPR by 13.2% (-0.4%, 32.9%) for gating, by 6.0% (-0.8%, 27.7%) for MLC tracking and by 11.1% (1.2%, 22.9%) for couch tracking. The three compensation techniques improved the target coverage in most cases. Gating showed better OAR sparing than MLC tracking or couch tracking.

CONCLUSIONS

Compensation of intrafraction prostate rotation with gating, MLC tracking and couch tracking was investigated experimentally for the first time. All three techniques improved the dosimetric accuracy, but residual motion-related dose errors remained due to the lack of rotation correction.

Adaptive Radiation Therapy (ART) Strategies and Technical Considerations: A State of the ART Review From NRG Oncology

The integration of adaptive radiation therapy (ART), or modifying the treatment plan during the treatment course, is becoming more widely available in clinical practice. ART offers strong potential for minimizing treatment-related toxicity while escalating or de-escalating target doses based on the dose to organs at risk. Yet, ART workflows add complexity into the radiation therapy planning and delivery process that may introduce additional uncertainties. This work sought to review presently available ART workflows and technological considerations such as image quality, deformable image registration, and dose accumulation. Quality assurance considerations for ART components and minimum recommendations are described. Personnel and workflow efficiency recommendations are provided, as is a summary of currently available clinical evidence supporting the implementation of ART. Finally, to guide future clinical trial protocols, an example ART physician directive and a physics template following standard NRG Oncology protocol is provided.

PET and MRI guided adaptive radiotherapy: Rational, feasibility and benefit

Adaptive radiotherapy (ART) corresponds to various replanning strategies aiming to correct for anatomical variations occurring during the course of radiotherapy. The goal of the article was to report the rational, feasibility and benefit of using PET and/or MRI to guide this ART strategy in various tumor localizations. The anatomical modifications defined by scanner taking into account tumour mobility and volume variation are not always sufficient to optimise treatment. The contribution of functional imaging by PET or the precision of soft tissue by MRI makes it possible to consider optimized ART. Today, the most important data for both PET and MRI are for lung, head and neck, cervical and prostate cancers. PET and MRI guided ART appears feasible and safe, however in a very limited clinical experience. Phase I/II studies should be therefore performed, before proposing cost-effectiveness comparisons in randomized trials and before using the approach in routine practice.

Quantifying the dose accumulation uncertainty after deformable image registration in head-and-neck radiotherapy

BACKGROUND AND PURPOSE

Deformable image registration (DIR) facilitated dose reconstruction and accumulation can be applied to assess delivered dose and verify the validity of the treatment plan during treatment. This retrospective study used in silico deformations based on clinically observed anatomical changes as ground truth to investigate the uncertainty of reconstructed and accumulated dose in head-and-neck radiotherapy (HNRT).

MATERIALS AND METHODS

A planning CT (pCT), cone beam CT (CBCT) from week one of treatment and three later CBCTs were selected for 12 HNRT patients. These images were used to generate in silico reference CBCTs and deformation vector fields (DVFs) as ground truth with B-spline DIR. Inverse consistency (IC) of voxels was assessed by determining their net displacement after successive application of the forward and backward DVF. The reconstructed dose based on demons DIR was compared to the ground truth to assess the structure-specific uncertainties of this DIR algorithm for inverse consistent and inverse inconsistent voxels.

RESULTS

Overall, 98.5% of voxels were inverse consistent with the 95% level of confidence range for dose reconstruction of a single fraction equal to [-2.3%; +2.1%], [-10.2%; +15.2%] and [-9.5%; +12.5%] relative to their planned dose for target structures, critical organs at risk (OARs) and non-critical OARs, respectively. Inverse inconsistent voxels generally showed a higher level of uncertainty.

CONCLUSION

The uncertainty in accumulated dose using DIR can be accurately quantified and incorporated in dose-volume histograms (DVHs). This method can be used to prospectively assess the adequacy of target coverage during treatment in an objective manner.