Dose tracking assessment for magnetic resonance guided adaptive radiotherapy of rectal cancers.

BACKGROUND: Magnetic resonance-guided adaptive radiotherapy (MRgART) at MR-Linac allows for plan optimisation on the MR-based synthetic CT (sCT) images, adjusting the target and organs at risk according to the patient's daily anatomy. Conversely, conventional linac image-guided radiotherapy (IGRT) involves rigid realignment of regions of interest to the daily anatomy, followed by the delivery of the reference computed tomography (CT) plan. This study aims to evaluate the effectiveness of MRgART versus IGRT for rectal cancer patients undergoing short-course radiotherapy, while also assessing the dose accumulation process to support the findings and determine its usefulness in enhancing treatment accuracy. METHODS: Nineteen rectal cancer patients treated with a 1.5 Tesla MR-Linac with a prescription dose of 25 Gy (5 Gy x 5) and undergoing daily adapted radiotherapy by plan optimization based on online MR-based sCT images, were included in this retrospective study. For each adapted plan ([Formula: see text]), a second plan ([Formula: see text]) was generated by recalculating the reference CT plan on the daily MR-based sCT images after rigid registration with the reference CT images to simulate the IGRT workflow. Dosimetry of [Formula: see text] and[Formula: see text]was compared for each fraction. Cumulative doses on the first and last fractions were evaluated for both workflows. The dosimetry per single fraction and the cumulative doses were compared using dose-volume histogram parameters. RESULTS: Ninety-five fractions delivered with MRgART were compared to corresponding simulated IGRT fractions. All MRgART fractions fulfilled the target clinical requirements. IGRT treatments did not meet the expected target coverage for 63 out of 94 fractions (67.0%), with 13 fractions showing a V95 median point percentage decrease of 2.78% (range, 1.65-4.16%), and 55 fractions exceeding the V107% threshold with a median value of 15.4 cc (range, 6.0-43.8 cc). For the bladder, the median [Formula: see text] values were 18.18 Gy for the adaptive fractions and 19.60 Gy for the IGRT fractions. Similarly the median [Formula: see text] values for the small bowel were 23.40 Gy and 25.69 Gy, respectively. No statistically significant differences were observed in the doses accumulated on the first or last fraction for the adaptive workflow, with results consistent with the single adaptive fractions. In contrast, accumulated doses in the IGRT workflow showed significant variations mitigating the high dose constraint, nevertheless, more than half of the patients still did not meet clinical requirements. CONCLUSIONS: MRgART for short-course rectal cancer treatments ensures that the dose delivered matches each fraction of the planned dose and the results are confirmed by the dose accumulation process, which therefore seems redundant. In contrast, IGRT may lead to target dose discrepancies and non-compliance with organs at risk constraints and dose accumulation can still highlight notable dosimetric differences.

Adaptive intensity modulated proton therapy using 4D robust planning: a proof-of-concept for the application of dose mimicking approach.

Objective.A four-dimensional robust optimisation (4DRO) is usually employed when the tumour respiratory motion needs to be addressed. However, it is computationally demanding, and an automated method is preferable for adaptive planning to avoid manual trial-and-error. This study proposes a 4DRO technique based on dose mimicking for automated adaptive planning.Approach.Initial plans for 4DRO intensity modulated proton therapy were created on an average CT for four patients with clinical target volume (CTV) in the lung, oesophagus, or pancreas, respectively. These plans were robustly optimised using three phases of four-dimensional computed tomography (4DCT) and accounting for setup and density uncertainties. Weekly 4DCTs were used for adaptive replanning, using a constant relative biological effectiveness (cRBE) of 1.1. Two methods were used: (1) template-based adaptive (TA) planning and (2) dose-mimicking-based adaptive (MA) planning. The plans were evaluated using variable RBE (vRBE) weighted doses and biologically consistent dose accumulation (BCDA).Main results.MA and TA plans had comparable CTV coverage except for one patient where the MA plan had a higher D98 and lower D2 but with an increased D2 in few organs at risk (OARs). CTV D98 deviations in non-adaptive plans from the initial plans were up to -7.2 percentage points (p.p.) in individual cases and -1.8 p.p. when using BCDA. For the OARs, MA plans showed a reduced mean dose and D2 compared to the TA plans, with few exceptions. The vRBE-weighted accumulated doses had a mean dose and D2 difference of up to 0.3 Gy and 0.5 Gy, respectively, in the OARs with respect to cRBE-weighted doses.Significance.MA plans indicate better performance in target coverage and OAR dose sparing compared to the TA plans in 4DRO adaptive planning. Moreover, MA method is capable of handling both forms of anatomical variation, namely, changes in density and relative shifts in the position of OARs.

A multi-institutional comparison of retrospective deformable dose accumulation for online adaptive magnetic resonance-guided radiotherapy.

Background and Purpose: Application of different deformable dose accumulation (DDA) solutions makes institutional comparisons after online-adaptive magnetic resonance-guided radiotherapy (OA-MRgRT) challenging. The aim of this multi-institutional study was to analyze accuracy and agreement of DDA-implementations in OA-MRgRT. Material and Methods: One gold standard (GS) case deformed with a biomechanical-model and five clinical cases consisting of prostate (2x), cervix, liver, and lymph node cancer, treated with OA-MRgRT, were analyzed. Six centers conducted DDA using institutional implementations. Deformable image registration (DIR) and DDA results were compared using the contour metrics Dice Similarity Coefficient (DSC), surface-DSC, Hausdorff-distance (HD95%), and accumulated dose-volume histograms (DVHs) analyzed via intraclass correlation coefficient (ICC) and clinical dosimetric criteria (CDC). Results: For the GS, median DDA errors ranged from 0.0 to 2.8 Gy across contours and implementations. DIR of clinical cases resulted in DSC > 0.8 for up to 81.3% of contours and a variability of surface-DSC values depending on the implementation. Maximum HD95%=73.3 mm was found for duodenum in the liver case. Although DVH ICC > 0.90 was found after DDA for all but two contours, relevant absolute CDC differences were observed in clinical cases: Prostate I/II showed maximum differences in bladder V28Gy (10.2/7.6%), while for cervix, liver, and lymph node the highest differences were found for rectum D2cm3 (2.8 Gy), duodenum Dmax (7.1 Gy), and rectum D0.5cm3 (4.6 Gy). Conclusion: Overall, high agreement was found between the different DIR and DDA implementations. Case- and algorithm-dependent differences were observed, leading to potentially clinically relevant results. Larger studies are needed to define future DDA-guidelines.

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.

Quantification of deformable image registration uncertainties for dose accumulation on head and neck cancer proton treatments.

PURPOSE: Head and neck cancer (HNC) patients in radiotherapy require adaptive treatment plans due to anatomical changes. Deformable image registration (DIR) is used in adaptive radiotherapy, e.g. for deformable dose accumulation (DDA). However, DIR's ill-posedness necessitates addressing uncertainties, often overlooked in clinical implementations. DIR's further clinical implementation is hindered by missing quantitative commissioning and quality assurance tools. This study evaluates one pathway for more quantitative DDA uncertainties. METHODS: For five HNC patients, each with multiple repeated CTs acquired during treatment, a simultaneous-integrated boost (SIB) plan was optimized. Recalculated doses were warped individually using multiple DIRs from repeated to reference CTs, and voxel-by-voxel dose ranges determined an error-bar for DDA. Followed by evaluating, a previously proposed early-stage DDA uncertainty estimation method tested for lung cancer, which combines geometric DIR uncertainties, dose gradients and their directional dependence, in the context of HNC. RESULTS: Applying multiple DIRs show dose differences, pronounced in high dose gradient regions. The patient with largest anatomical changes (-13.1 % in ROI body volume), exhibited 33 % maximum uncertainty in contralateral parotid, with 54 % of voxels presenting an uncertainty >5 %. Accumulation over multiple CTs partially mitigated uncertainties. The estimation approach predicted 92.6 % of voxels within ±5 % to the reference dose uncertainty across all patients. CONCLUSIONS: DIR variations impact accumulated doses, emphasizing DDA uncertainty quantification's importance for HNC patients. Multiple DIR dose warping aids in quantifying DDA uncertainties. An estimation approach previously described for lung cancer was successfully validated for HNC, for SIB plans, presenting different dose gradients, and for accumulated treatments.

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.

Multi-centre evaluation of variation in cumulative dose assessment in reirradiation scenarios.

BACKGROUND AND PURPOSE: Safe reirradiation relies on assessment of cumulative doses to organs at risk (OARs) across multiple treatments. Different clinical pathways can result in inconsistent estimates. Here, we quantified the consistency of cumulative dose to OARs across multi-centre clinical pathways. MATERIAL AND METHODS: We provided DICOM planning CT, structures and doses for two reirradiation cases: head & neck (HN) and lung. Participants followed their standard pathway to assess the cumulative physical and EQD2 doses (with provided α/ß values), and submitted DVH metrics and a description of their pathways. Participants could also submit physical dose distributions from Course 1 mapped onto the CT of Course 2 using their best available tools. To assess isolated impact of image registrations, a single observer accumulated each submitted spatially mapped physical dose for every participating centre. RESULTS: Cumulative dose assessment was performed by 24 participants. Pathways included rigid (n = 15), or deformable (n = 5) image registration-based 3D dose summation, visual inspection of isodose line contours (n = 1), or summation of dose metrics extracted from each course (n = 3). Largest variations were observed in near-maximum cumulative doses (25.4 - 41.8 Gy for HN, 2.4 - 33.8 Gy for lung OARs), with lower variations in volume/dose metrics to large organs. A standardised process involving spatial mapping of the first course dose to the second course CT followed by summation improved consistency for most near-maximum dose metrics in both cases. CONCLUSION: Large variations highlight the uncertainty in reporting cumulative doses in reirradiation scenarios, with implications for outcome analysis and understanding of published doses. Using a standardised workflow potentially including spatially mapped doses improves consistency in determination of accumulated dose in reirradiation scenarios.

Radiation therapy margin reduction for patients with localized prostate cancer: A prospective study of the dosimetric impact and quality of life.

OBJECTIVES: To investigate the impact of reducing Clinical Target Volume (CTV) to Planning Target Volume (PTV) margins on delivered radiation therapy (RT) dose and patient reported quality-of-life (QOL) for patients with localized prostate cancer. METHODS: Twenty patients were included in a single institution IRB-approved prospective study. Nine were planned with reduced margins (4 mm at prostate/rectum interface, 5 mm elsewhere), and 11 with standard margins (6/10 mm). Cumulative delivered dose was calculated using deformable dose accumulation. Each daily CBCT dataset was deformed to the planning CT (pCT), dose was computed, and accumulated on the resampled pCT using a parameter-optimized, B-spline algorithm (Elastix, ITK/VTK). EPIC-26 patient reported QOL was prospectively collected pre-treatment, post-treatment, and at 2-, 6-, 12-, 18-, 24-, 36-, 48-, and 60-month follow-ups. Post -RT QOL scores were baseline corrected and standardized to a [0-100] scale using EPIC-26 methodology. Correlations between QOL scores and dosimetric parameters were investigated, and the overall QOL differences between the two groups (QOLMargin-reduced -QOLcontrol ) were calculated. RESULTS: The median QOL follow-up length for the 20 patients was 48 months. Difference between delivered dose and planned dose did not reach statistical significance (p > 0.1) for both targets and organs at risk between the two groups. At 4 years post-RT, standardized mean QOLMargin-reduced -QOLcontrol were improved for Urinary Incontinence, Urinary Irritative/Obstructive, Bowel, and Sexual EPIC domains by 3.5, 14.8, 10.2, and 16.1, respectively (higher values better). The control group showed larger PTV/rectum and PTV/bladder intersection volumes (7.2 ± 5.8, 18.2 ± 8.1 cc) than the margin-reduced group (2.6 ± 1.8, 12.5 ± 8.3 cc), though the dose to these intersection volumes did not reach statistical significance (p > 0.1) between the groups. PTV/rectum intersection volume showed a moderate correlation (r = -0.56, p < 0.05) to Bowel EPIC domain. CONCLUSIONS: Results of this prospective study showed that margin-reduced group exhibited clinically meaningful improvement of QOL without compromising the target dose coverage.

Review and recommendations on deformable image registration uncertainties for radiotherapy applications.

Deformable image registration (DIR) is a versatile tool used in many applications in radiotherapy (RT). DIR algorithms have been implemented in many commercial treatment planning systems providing accessible and easy-to-use solutions. However, the geometric uncertainty of DIR can be large and difficult to quantify, resulting in barriers to clinical practice. Currently, there is no agreement in the RT community on how to quantify these uncertainties and determine thresholds that distinguish a good DIR result from a poor one. This review summarises the current literature on sources of DIR uncertainties and their impact on RT applications. Recommendations are provided on how to handle these uncertainties for patient-specific use, commissioning, and research. Recommendations are also provided for developers and vendors to help users to understand DIR uncertainties and make the application of DIR in RT safer and more reliable.

Fractionation versus Adaptation for Compensation of Target Volume Changes during Online Adaptive Radiotherapy for Bladder Cancer: Answers from a Prospective Registry.

Online adaptive radiotherapy (ART) allows adaptation of the dose distribution to the anatomy captured by with pre-adaptation imaging. ART is time-consuming, and thus intra-fractional deformations can occur. This prospective registry study analyzed the effects of intra-fraction deformations of clinical target volume (CTV) on the equivalent uniform dose (EUDCTV) of focal bladder cancer radiotherapy. Using margins of 5-10 mm around CTV on pre-adaptation imaging, intra-fraction CTV-deformations found in a second imaging study reduced the 10th percentile of EUDCTV values per fraction from 101.1% to 63.2% of the prescribed dose. Dose accumulation across fractions of a series was determined with deformable-image registration and worst-case dose accumulation that maximizes the correlation of cold spots. A strong fractionation effect was demonstrated-the EUDCTV was above 95% and 92.5% as determined by the two abovementioned accumulation methods, respectively, for all series of dose fractions. A comparison of both methods showed that the fractionation effect caused the EUDCTV of a series to be insensitive to EUDCTV-declines per dose fraction, and this could be explained by the small size and spatial variations of cold spots. Therefore, ART for each dose fraction is unnecessary, and selective ART for fractions with large inter-fractional deformations alone is sufficient for maintaining a high EUDCTV for a radiotherapy series.

Deep learning based uncertainty prediction of deformable image registration for contour propagation and dose accumulation in online adaptive radiotherapy.

Objective.Online adaptive radiotherapy aims to fully leverage the advantages of highly conformal therapy by reducing anatomical and set-up uncertainty, thereby alleviating the need for robust treatments. This requires extensive automation, among which is the use of deformable image registration (DIR) for contour propagation and dose accumulation. However, inconsistencies in DIR solutions between different algorithms have caused distrust, hampering its direct clinical use. This work aims to enable the clinical use of DIR by developing deep learning methods to predict DIR uncertainty and propagating it into clinically usable metrics.Approach.Supervised and unsupervised neural networks were trained to predict the Gaussian uncertainty of a given deformable vector field (DVF). Since both methods rely on different assumptions, their predictions differ and were further merged into a combined model. The resulting normally distributed DVFs can be directly sampled to propagate the uncertainty into contour and accumulated dose uncertainty.Main results.The unsupervised and combined models can accurately predict the uncertainty in the manually annotated landmarks on the DIRLAB dataset. Furthermore, for 5 patients with lung cancer, the propagation of the predicted DVF uncertainty into contour uncertainty yielded for both methods anexpected calibration errorof less than 3%. Additionally, theprobabilisticly accumulated dose volume histograms(DVH) encompass well the accumulated proton therapy doses using 5 different DIR algorithms. It was additionally shown that the unsupervised model can be used for different DIR algorithms without the need for retraining.Significance.Our work presents first-of-a-kind deep learning methods to predict the uncertainty of the DIR process. The methods are fast, yield high-quality uncertainty estimates and are useable for different algorithms and applications. This allows clinics to use DIR uncertainty in their workflows without the need to change their DIR implementation.

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.

CBCT-based deformable dose accumulation of external beam radiotherapy in cervical cancer.

Background: Delivered radiotherapy doses do not exactly match those planned for a course of treatment, largely due to inter-fraction changes in anatomy. In this study, accumulated delivered dose was calculated for a sample of cervical cancer patients, by deformably registering daily cone beam computed tomography (CBCT) images to the planning computed tomography (CT) scan. Planned and accumulated doses were compared for the clinical target volume (CTV), bladder, and rectum.Material and Methods: For 10 patients receiving 45 Gy in 25 fractions of external beam radiotherapy, daily dose distributions were calculated on CBCT. These images were deformed onto the planning CT and the dose was accumulated using Velocity 4.1 (Varian Medical Systems, Palo Alto, USA). The quality of deformable image registration was evaluated visually and by calculating Dice similarity coefficients and mean distance to agreement.Results: V95%>99% was achieved for the primary CTV in 9/10 patients for the planned dose distribution and 7/10 patients for the accumulated dose distribution. Primary CTV coverage by 95% of the prescription dose was reduced in one patient, due to an increase in anterior-posterior separation. Comparison of planned and accumulated dose volume histograms (DVHs) for the bladder and rectum found agreement within 5% at low and intermediate doses, but differences exceeded 20% at higher doses. Direct addition of CBCT DVHs was seen to be a poor estimate for the accumulated DVH at higher doses.Conclusion: Computation of delivered radiotherapy dose that accounts for inter-fraction anatomical changes is important for establishing dose-effect relationships. Updating delivered dose distributions after each fraction would support informed clinical decision making on any potential treatment interventions.

Patient-specific quality assurance for deformable IMRT/IMPT dose accumulation: Proposition and validation of energy conservation based validation criterion.

BACKGROUND: Deformable image registration (DIR)-based dose accumulation (DDA) is regularly used in adaptive radiotherapy research. However, the applicability and reliability of DDA for direct clinical usage are still being debated. One primary concern is the validity of DDA, particularly for scenarios with substantial anatomical changes, for which energy-conservation problems were observed in conceptual studies. PURPOSE: We present and validate an energy-conservation (EC)-based DDA validation workflow and further investigate its usefulness for actual patient data, specifically for lung cancer cases. METHODS: For five non-small cell lung cancer (NSCLC) patients, DDA based on five selective DIR methods were calculated for five different treatment plans, which include one intensity-modulated photon therapy (IMRT), two intensity-modulated proton therapy (IMPT), and two combined proton-photon therapy (CPPT) plans. All plans were optimized on the planning CT (planCT) acquired in deep inspiration breath-hold (DIBH) and were re-optimized on the repeated DIBH CTs of three later fractions. The resulting fractional doses were warped back to the planCT using each DIR. An EC-based validation of the accumulation process was implemented and applied to all DDA results. Correlations between relative organ mass/volume variations and the extent of EC violation were then studied using Bayesian linear regression (BLR). RESULTS: For most OARs, EC violation within 10% is observed. However, for the PTVs and GTVs with substantial regression, severe overestimation of the fractional energy was found regardless of treatment type and applied DIR method. BLR results show that EC violation is linearly correlated to the relative mass variation (R^2 > 0.95) and volume variation (R^2 > 0.60). CONCLUSION: DDA results should be used with caution in regions with high mass/volume variation for intensity-based DIRs. EC-based validation is a useful approach to provide patient-specific quality assurance of the validity of DDA in radiotherapy.

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

Purpose: This study assesses the impact of intra-fraction motion and PTV margin size on target coverage for patients undergoing radiation treatment of pelvic oligometastases. Dosimetric sparing of the bowel as a function of the PTV margin is also evaluated. Materials and methods: Seven patients with pelvic oligometastases previously treated on our MR-linac (35 Gy in 5 fractions) were included in this study. Retrospective adaptive plans were created for each fraction on the daily MRI datasets using PTV margins of 5 mm, 3 mm, and 2 mm. Dosimetric constraint violations and GTV coverage were measured as a function of PTV margin size. The impact of intra-fraction motion on GTV coverage was assessed by tracking the GTV position on the cine MR images acquired during treatment delivery and creating an intra-fraction dose distribution for each IMRT beam. The intra-fraction dose was accumulated for each fraction to determine the total dose delivered to the target for each PTV size. Results: All OAR constraints were achieved in 85.7%, 94.3%, and 100.0% of fractions when using 5 mm, 3 mm, and 2 mm PTV margins while scaling to 95% PTV coverage. Compared to plans with a 5 mm PTV margin, there was a 27.4 ± 12.3% (4.0 ± 2.2 Gy) and an 18.5 ± 7.3% (2.7 ± 1.4 Gy) reduction in the bowel D0.5cc dose for 2 mm and 3 mm PTV margins, respectively. The target dose (GTV V35 Gy) was on average 100.0 ± 0.1% (99.6 - 100%), 99.6 ± 1.0% (97.2 - 100%), and 99.0 ± 1.4% (95.0 - 100%), among all fractions for the 5 mm, 3 mm, and 2 mm PTV margins on the adaptive plans when accounting for intra-fraction motion, respectively. Conclusion: A 2 mm PTV margin achieved a minimum of 95% GTV coverage while reducing the dose to the bowel for all patients.

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.

Practical and technical key challenges in head and neck adaptive radiotherapy: The GORTEC point of view.

Anatomical variations occur during head and neck (H&N) radiotherapy (RT) treatment. These variations may result in underdosage to the target volume or overdosage to the organ at risk. Replanning during the treatment course can be triggered to overcome this issue. Due to technological, methodological and clinical evolutions, tools for adaptive RT (ART) are becoming increasingly sophisticated. The aim of this paper is to give an overview of the key steps of an H&N ART workflow and tools from the point of view of a group of French-speaking medical physicists and physicians (from GORTEC). Focuses are made on image registration, segmentation, estimation of the delivered dose of the day, workflow and quality assurance for an implementation of H&N offline and online ART. Practical recommendations are given to assist physicians and medical physicists in a clinical workflow.

Liver SBRT dose accumulation to assess the impact of anatomic variations on normal tissue doses and toxicity in patients treated with concurrent sorafenib.

BACKGROUND AND PURPOSE: Unexpected liver volume reductions occurred during trials of liver SBRT and concurrent sorafenib. The aims were to accumulate liver SBRT doses to assess the impact of these anatomic variations on normal tissue dose parameters and toxicity. MATERIALS AND METHODS: Thirty-two patients with hepatocellular carcinoma (HCC) or metastases treated on trials of liver SBRT (30-57 Gy, 6 fractions) and concurrent sorafenib were analyzed. SBRT doses were accumulated using biomechanical deformable registration of daily cone-beam CT. Dose deviations (accumulated-planned) for normal tissues were compared for patients with liver volume reductions > 100 cc versus stable volumes, and accumulated doses were reported for three patients with grade 3-5 luminal gastrointestinal toxicities. RESULTS: Patients with reduced (N = 12) liver volumes had larger mean deviations of 0.4-1.3 Gy in normal tissues, versus -0.2-0.4 Gy for stable cases (N = 20), P > 0.05. Deviations > 5% of the prescribed dose occurred in both groups. Two HCC patients with toxicities to small and large bowel had liver volume reductions and deviations to the maximum dose of 4% (accumulated 36.9 Gy) and 3% (accumulated 33.4 Gy) to these organs respectively. Another HCC patient with a toxicity of unknown location plus tumor rupture, had stable liver volumes and deviations to luminal organs of -6% to 4.5% (accumulated < 30.5 Gy). CONCLUSION: Liver volume reductions during SBRT and concurrent sorafenib were associated with larger increases in accumulated dose to normal tissues versus stable liver volumes. These dosimetric changes may have further contributed to toxicities in HCC patients who have higher baseline risks.

Sigmoid dose accumulation and reporting for multifractionated brachytherapy for cervical cancer: Methodological development of sigmoid points through virtual endoscopic method.

PURPOSE: The sigmoid is an important organ at risk for gynecological brachytherapy (BT). However, the reliability of localization of high-dose regions during multi-fractionated treatment is limited. This work reports the methodological development of sigmoid points to summate multi-fractionated doses. METHODS AND MATERIAL: Ten paired MRI data sets of ring-based intracavitary brachytherapy were obtained. Simulating a virtual endoscope, a reference line was created along the central axis of the anorectosigmoid for each implant. A trendline was generated, and linear dose was determined. Three-dimensional (3D) coordinates of high-dose regions were identified, and overlap was determined. In the next step, 3D coordinates of high-dose sigmoid points were localized in reference to cervical os and re-verified for location in reference to sigmoid lumen and corroboration with 2cc doses. With minor modifications, sigmoid points were proposed. RESULTS: In 6 of 10 patients, high-dose regions co-localized in subsequent fractions of BT. Three high-dose regions were identified along the sigmoid length and proposed as sigmoid points in reference to cervical os. (S1'= 0.5 cm right, 1.5 cm posterior, and 2.4 cm cranial; S2' = 0.3 cm anterior and 4.5 cm cranial; S3' = 2.7 cm left, 3 cm anterior, and 3.6 cm cranial to the cervical os). S1' and S2' were located in the sigmoid in 70% and 60% of data sets. The mean difference between D2cc and S1'/S2' was 0.30 Gy and 1.06 Gy respectively. S3' had limited corroboration to sigmoid lumen or 2 cc doses. The points S1' and S2' were further modified (minor) for applicability and proposed as sigmoid points 1 and 2 (SP1 0.5 right,1.5 posterior and 2.5 cm cranial to cervical os and SP2 (0.5 cm anterior and 4.5 cm cranial to cervical os)). CONCLUSION: SP1 and SP 2 are proposed as a surrogate for 2 cc sigmoid doses and may provide a method of reliable inter-fraction dose summation. This pilot work requires further validation.

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

BACKGROUND: Stereotactic body radiation therapy (SBRT) of central lung tumors with photon or proton therapy has a risk of increased toxicity. Treatment planning studies comparing accumulated doses for state-of-the-art treatment techniques, such as MR-guided radiotherapy (MRgRT) and intensity modulated proton therapy (IMPT), are currently lacking. PURPOSE: We conducted a comparison of accumulated doses for MRgRT, robustly optimized non-adaptive IMPT, and online adaptive IMPT for central lung tumors. A special focus was set on analyzing the accumulated doses to the bronchial tree, a parameter linked to high-grade toxicities. METHODS: Data of 18 early-stage central lung tumor patients, treated at a 0.35 T MR-linac in eight or five fractions, were analyzed. Three gated treatment scenarios were compared: (S1) online adaptive MRgRT, (S2) non-adaptive IMPT, and (S3) online adaptive IMPT. The treatment plans were recalculated or reoptimized on the daily imaging data acquired during MRgRT, and accumulated over all treatment fractions. Accumulated dose-volume histogram (DVH) parameters of the gross tumor volume (GTV), lung, heart, and organs-at-risk (OARs) within 2 cm of the planning target volume (PTV) were extracted for each scenario and compared in Wilcoxon signed-rank tests between S1 & S2, and S1 & S3. RESULTS: The accumulated GTV D98% was above the prescribed dose for all patients and scenarios. Significant reductions (p < 0.05) of the mean ipsilateral lung dose (S2: -8%; S3: -23%) and mean heart dose (S2: -79%; S3: -83%) were observed for both proton scenarios compared to S1. The bronchial tree D0.1cc was significantly lower for S3 (S1: 48.1 Gy; S3: 39.2 Gy; p = 0.005), but not significantly different for S2 (S2: 45.0 Gy; p = 0.094), compared to S1. The D0.1cc for S2 and S3 compared to S1 was significantly (p < 0.05) smaller for OARs within 1-2 cm of the PTV (S1: 30.2 Gy; S2: 24.6 Gy; S3: 23.1 Gy), but not significantly different for OARs within 1 cm of the PTV. CONCLUSIONS: A significant dose sparing potential of non-adaptive and online adaptive proton therapy compared to MRgRT for OARs in close, but not direct proximity of central lung tumors was identified. The near-maximum dose to the bronchial tree was not significantly different for MRgRT and non-adaptive IMPT. Online adaptive IMPT achieved significantly lower doses to the bronchial tree compared to MRgRT.