Importance of deformable image registration and biological dose summation in planning of radiotherapy retreatments

Reirradiation - still navigating uncharted waters?

With the emergence of high-precision radiotherapy technologies such as stereotactic ablative radiotherapy (SABR), MR guided brachytherapy, image guided intensity modulated photon and proton radiotherapy and most recently daily adaptive radiotherapy, reirradiation is increasingly recognized as a viable treatment option for many patients. This includes those with recurrent, metastatic or new malignancies post initial radiotherapy. The primary challenge in reirradiation lies in balancing tumor control against the risk of severe toxicity from cumulative radiation doses to previously irradiated normal tissue. Although technology for precise delivery has advanced at a fast pace, clinical practice of reirradiation still mostly relies on individual expertise, as prospective evidence is scarce, the level of reporting in clinical studies is not standardized and of low quality - especially with respect to cumulative doses received by organs at risk. A recent ESTRO/EORTC initiative proposed a standardized definition of reirradiation and formulated general requirements for minimal reporting in clinical studies [1]. As a consequence we found it timely to convene for an international and interdisciplinary meeting with experts in the field to summarize the current evidence, identify knowledge gaps and explore which best practices can be derived for safe reirradiation. The meeting was held on 15.06.2023 in Zurich and was endorsed by the scientific societies SASRO, DEGRO and ESTRO. Here, we report on available evidence and research priorities in the field of reirradiation, as discussed during the meeting.

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.

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.

Brain Re-Irradiation Robustly Accounting for Previously Delivered Dose

(1) Background: The STRIDeR (Support Tool for Re-Irradiation Decisions guided by Radiobiology) planning pathway aims to facilitate anatomically appropriate and radiobiologically meaningful re-irradiation (reRT). This work evaluated the STRIDeR pathway for robustness compared to a more conservative manual pathway. (2) Methods: For ten high-grade glioma reRT patient cases, uncertainties were applied and cumulative doses re-summed. Geometric uncertainties of 3, 6 and 9 mm were applied to the background dose, and LQ model robustness was tested using α/β variations (values 1, 2 and 5 Gy) and the linear quadratic linear (LQL) model δ variations (values 0.1 and 0.2). STRIDeR robust optimised plans, incorporating the geometric and α/β uncertainties during optimisation, were also generated. (3) Results: The STRIDeR and manual pathways both achieved clinically acceptable plans in 8/10 cases but with statistically significant improvements in the PTV D98% (p < 0.01) for STRIDeR. Geometric and LQ robustness tests showed comparable robustness within both pathways. STRIDeR plans generated to incorporate uncertainties during optimisation resulted in a superior plan robustness with a minimal impact on PTV dose benefits. (4) Conclusions: Our results indicate that STRIDeR pathway plans achieved a similar robustness to manual pathways with improved PTV doses. Geometric and LQ model uncertainties can be incorporated into the STRIDeR pathway to facilitate robust optimisation.

Treatment plan optimisation for reirradiation

BACKGROUND

The STRIDeR (Support Tool for Re-Irradiation Decisions guided by Radiobiology) project aims to create a clinically viable re-irradiation planning pathway within a commercial treatment planning system (TPS). Such a pathway should account for previously delivered dose, voxel-by-voxel, taking fractionation effects, tissue recovery and anatomical changes into account. This work presents the workflow and technical solutions in the STRIDeR pathway.

METHODS

The pathway was implemented in RayStation (version 9B DTK) to allow an original dose distribution to be used as background dose to guide optimisation of re-irradiation plans. Organ at risk (OAR) planning objectives in equivalent dose in 2 Gy fractions (EQD2) were applied cumulatively across the original and re-irradiation treatments, with optimisation of the re-irradiation plan performed voxel-by-voxel in EQD2. Different approaches to image registration were employed to account for anatomical change. Data from 21 patients who received pelvic Stereotactic Ablative Radiotherapy (SABR) re-irradiation were used to illustrate the use of the STRIDeR workflow. STRIDeR plans were compared to those produced using a standard manual method.

RESULTS

The STRIDeR pathway resulted in clinically acceptable plans in 20/21 cases. Compared to plans produced using the laborious manual method, less constraint relaxation was required or higher re-irradiation doses could be prescribed in 3/21.

CONCLUSION

The STRIDeR pathway used background dose to guide radiobiologically meaningful, anatomically-appropriate re-irradiation treatment planning within a commercial TPS. This provides a standardised and transparent approach, offering more informed re-irradiation and improved cumulative OAR dose evaluation.

Evaluation of Biological Effective Dose in Gamma Knife Staged Stereotactic Radiosurgery for Large Brain Metastases

Objective

Gamma knife (GK) staged stereotactic radiosurgery (Staged-SRS) has become an effective treatment option for large brain metastases (BMs); however, it has been challenging to evaluate the total dose because of tumor shrinkage between two staged sessions. This study aims to evaluate total biological effective dose (BED) in Staged-SRS, and to compare the BED with those in single-fraction SRS (SF-SRS) and hypo-fractionated SRS (HF-SRS).

Methods

Patients treated with GK Staged-SRS at a single institution were retrospectively included. The dose delivered in two sessions of Staged-SRS was summed using the deformable image registration. Each patient was replanned for SF-SRS and HF-SRS. The total BEDs were computed using the linear-quadratic model. Tumor BED_98% and brain V_84Gy2, equivalent to V_12Gy in SF-SRS, were compared between SF-SRS, HF-SRS, and Staged-SRS plans with the Wilcoxon test.

Results

Twelve patients with 24 BMs treated with GK Staged-SRS were identified. We observed significant differences (p < 0.05) in tumor BED_98% but comparable brain V_84Gy2 (p = 0.677) between the Staged-SRS and SF-SRS plans. No dosimetric advantages of Staged-SRS over HF-SRS were observed. Tumor BED_98% in the HF-SRS plans were significantly higher than those in the Staged-SRS plans (p < 0.05). Despite the larger PTVs, brain V_84Gy2 in the HF-SRS plans remained lower (p < 0.05).

Conclusion

We presented an approach to calculate the composite BEDs delivered to both tumor and normal brain tissue in Staged-SRS. Compared to SF-SRS, Staged-SRS delivers a higher dose to tumor but maintains a comparable dose to normal brain tissue. Our results did not show any dosimetric advantages of Staged-SRS over HF-SRS.

Adaptive proton therapy

Radiation therapy treatments are typically planned based on a single image set, assuming that the patient's anatomy and its position relative to the delivery system remains constant during the course of treatment. Similarly, the prescription dose assumes constant biological dose-response over the treatment course. However, variations can and do occur on multiple time scales. For treatment sites with significant intra-fractional motion, geometric changes happen over seconds or minutes, while biological considerations change over days or weeks. At an intermediate timescale, geometric changes occur between daily treatment fractions. Adaptive radiation therapy is applied to consider changes in patient anatomy during the course of fractionated treatment delivery. While traditionally adaptation has been done off-line with replanning based on new CT images, online treatment adaptation based on on-board imaging has gained momentum in recent years due to advanced imaging techniques combined with treatment delivery systems. Adaptation is particularly important in proton therapy where small changes in patient anatomy can lead to significant dose perturbations due to the dose conformality and finite range of proton beams. This review summarizes the current state-of-the-art of on-line adaptive proton therapy and identifies areas requiring further research.

Biologically consistent dose accumulation using daily patient imaging

BACKGROUND

This work addresses a basic inconsistency in the way dose is accumulated in radiotherapy when predicting the biological effect based on the linear quadratic model (LQM). To overcome this inconsistency, we introduce and evaluate the concept of the total biological dose, bEQD_d.

METHODS

Daily computed tomography imaging of nine patients treated for prostate carcinoma with intensity-modulated radiotherapy was used to compute the delivered deformed dose on the basis of deformable image registration (DIR). We compared conventional dose accumulation (DA) with the newly introduced bEQD_d, a new method of accumulating biological dose that considers each fraction dose and tissue radiobiology. We investigated the impact of the applied fractionation scheme (conventional/hypofractionated), uncertainties induced by the DIR and by the assigned α/β-value.

RESULTS

bEQD_d was systematically higher than the conventionally accumulated dose with difference hot spots of 3.3-4.9 Gy detected in six out of nine patients in regions of high dose gradient in the bladder and rectum. For hypofractionation, differences are up to 8.4 Gy. The difference amplitude was found to be in a similar range to worst-case uncertainties induced by DIR and was higher than that induced by α/β.

CONCLUSION

Using bEQD_d for dose accumulation overcomes a potential systematic inaccuracy in biological effect prediction based on accumulated dose. Highest impact is found for serial-type late responding organs at risk in dose gradient regions and for hypofractionation. Although hot spot differences are in the order of several Gray, in dose-volume parameters there is little difference compared with using conventional or biological DA. However, when local dose information is used, e.g. dose surface maps, difference hot spots can potentially change outcomes of dose-response modelling and adaptive treatment strategies.

Photons or protons for reirradiation in (non-)small cell lung cancer: Results of the multicentric ROCOCO in silico study

OBJECTIVE

Locally recurrent disease is of increasing concern in (non-)small cell lung cancer [(N)SCLC] patients. Local reirradiation with photons or particles may be of benefit to these patients. In this multicentre in silico trial performed within the Radiation Oncology Collaborative Comparison (ROCOCO) consortium, the doses to the target volumes and organs at risk (OARs) were compared when using several photon and proton techniques in patients with recurrent localised lung cancer scheduled to undergo reirradiation.

METHODS

24 consecutive patients with a second primary (N)SCLC or recurrent disease after curative-intent, standard fractionated radio(chemo)therapy were included in this study. The target volumes and OARs were centrally contoured and distributed to the participating ROCOCO sites. Remaining doses to the OARs were calculated on an individual patient's basis. Treatment planning was performed by the participating site using the clinical treatment planning system and associated beam characteristics.

RESULTS

Treatment plans for all modalities (five photon and two proton plans per patient) were available for 22 patients (N = 154 plans). 3D-conformal photon therapy and double-scattered proton therapy delivered significantly lower doses to the target volumes. The highly conformal techniques, i.e., intensity modulated radiation therapy (IMRT), volumetric modulated arc therapy (VMAT), CyberKnife, TomoTherapy and intensity-modulated proton therapy (IMPT), reached the highest doses in the target volumes. Of these, IMPT was able to statistically significantly decrease the radiation doses to the OARs.

CONCLUSION

Highly conformal photon and proton beam techniques enable high-dose reirradiation of the target volume. They, however, significantly differ in the dose deposited in the OARs. The therapeutic options, i.e., reirradiation or systemic therapy, need to be carefully weighed and discussed with the patients.

ADVANCES IN KNOWLEDGE

Highly conformal photon and proton beam techniques enable high-dose reirradiation of the target volume. In light of the abilities of the various highly conformal techniques to spare specific OARs, the therapeutic options need to be carefully weighed and patients included in the decision-making process.

Normal tissue complication probability (NTCP) models for modern radiation therapy

Mathematical models of normal tissue complication probability (NTCP) able to robustly predict radiation-induced morbidities (RIM) play an essential role in the identification of a personalized optimal plan, and represent the key to maximizing the benefits of technological advances in radiation therapy (RT). Most modern RT techniques pose, however, new challenges in estimating the risk of RIM. The aim of this report is to schematically review NTCP models in the framework of advanced radiation therapy techniques. Issues relevant to hypofractionated stereotactic body RT and ion beam therapy are critically reviewed. Reirradiation scenarios for new or recurrent malignances and NTCP are also illustrated. A new phenomenological approach to predict RIM is suggested.

Deformable image registration for radiation therapy: principle, methods, applications and evaluation

Background: Deformable image registration (DIR) is increasingly used in the field of radiation therapy (RT) to account for anatomical deformations. The aims of this paper are to describe the main applications of DIR in RT and discuss current DIR evaluation methods. Methods: Articles on DIR published from January 2000 to October 2018 were extracted from PubMed and Science Direct. Our search was restricted to articles that report data obtained from humans, were written in English, and address DIR methods for RT. A total of 207 articles were selected from among 2506 identified in the search process. Results: At planning, DIR is used for organ delineation using atlas-based segmentation, deformation-based planning target volume definition, functional planning and magnetic resonance imaging-based dose calculation. In image-guided RT, DIR is used for contour propagation and dose calculation on per-treatment imaging. DIR is also used to determine the accumulated dose from fraction to fraction in external beam RT and brachytherapy, both for dose reporting and adaptive RT. In the case of re-irradiation, DIR can be used to estimate the cumulated dose of the two irradiations. Finally, DIR can be used to predict toxicity in voxel-wise population analysis. However, the evaluation of DIR remains an open issue, especially when dealing with complex cases such as the disappearance of matter. To quantify DIR uncertainties, most evaluation methods are limited to geometry-based metrics. Software companies have now integrated DIR tools into treatment planning systems for clinical use, such as contour propagation and fraction dose accumulation. Conclusions: DIR is increasingly important in RT applications, from planning to toxicity prediction. DIR is routinely used to reduce the workload of contour propagation. However, its use for complex dosimetric applications must be carefully evaluated by combining quantitative and qualitative analyses.

Deformable Registration for Dose Accumulation

As deformable image registration makes its way into the clinical routine, the summation of doses from fractionated treatment regimens to evaluate cumulative doses to targets and healthy tissues is also becoming a frequently utilized tool in the context of image-guided adaptive radiotherapy. Accounting for daily geometric changes using deformable image registration and dose accumulation potentially enables a better understanding of dose-volume-effect relationships, with the goal of translation of this knowledge to personalization of treatment, to further enhance treatment outcomes. Treatment adaptation involving image deformation requires patient-specific quality assurance of the image registration and dose accumulation processes, to ensure that uncertainties in the 3D dose distributions are identified and appreciated from a clinical relevance perspective. While much research has been devoted to identifying and managing the uncertainties associated with deformable image registration and dose accumulation approaches, there are still many unanswered questions. Here, we provide a review of current deformable image registration and dose accumulation techniques, and related clinical application. We also discuss salient issues that need to be deliberated when applying deformable algorithms for dose mapping and accumulation in the context of adaptive radiotherapy and response assessment.

A Systematic Post-QUANTEC Review of Tolerance Doses for Late Toxicity After Prostate Cancer Radiation Therapy

PURPOSE

The aims of this study were to systematically review tolerance doses for late distinct gastrointestinal (GI), genitourinary (GU), and sexual dysfunction (SD) symptoms after external beam radiation therapy (EBRT) alone and treatments involving brachytherapy (BT) for prostate cancer after Quantitative Analysis of Normal Tissue Effects in the Clinic (QUANTEC) and ultimately to perform quantitative syntheses of identified dose/volume tolerances represented by dose-volume histogram (DVH) thresholds, that is, statistically significant (P ≤ .05) cutoff points between symptomatic and asymptomatic patients in a certain study.

METHODS AND MATERIALS

PubMed was scrutinized for full-text articles in English after QUANTEC (January 1, 2010). The inclusion criteria were randomized controlled trials, case-control studies, or cohort studies with tolerance doses for late distinct symptoms ≥3 months after primary radiation therapy for prostate cancer (N > 30). All DVH thresholds were converted into equivalent doses in 2-Gy fractions (EQD2α/β) and were fitted with a linear or linear-quadratic function (goodness of fit, R2). The review was registered on PROSPERO (CRD42016042464).

RESULTS

From 33 identified studies, which included 36 to 746 patients per symptom domain, the majority of dose/volume tolerances were derived for GI toxicity after EBRT alone (GI, 97 thresholds; GU, 8 thresholds; SD, 1 threshold). For 5 symptoms (defecation urgency, diarrhea, fecal incontinence, proctitis, and rectal bleeding), relationships between dose/volume tolerances across studies (R2 = 0.93 [0.82-1.00]), and across symptoms, leading to a curve for overall GI toxicity (R2 = 0.98), could be determined. For these symptoms, mainly rectal thresholds were found throughout low and high doses (10 Gy ≤ equivalent dose in 2-Gy fractions using α/β = 3Gy (EQD23) ≤ 50 Gy and 55 Gy ≤ EQD23 ≤ 78 Gy, respectively). For BT with or without EBRT, dose/volume tolerances were also mainly identified for GI toxicity (GI, 14 thresholds; GU, 4 thresholds; SD, 2 thresholds) with the largest number of DVH thresholds concerning rectal bleeding (5 thresholds).

CONCLUSIONS

Updated dose/volume tolerances after QUANTEC were found for 17 GI, GU, or SD symptoms. A DVH curve described the relationship between dose/volume tolerances across 5 GI symptoms after EBRT alone. Restricting treatments for EBRT alone using the lower boundaries of this curve is likely to limit overall GI toxicity, but this should be explored prospectively. Dose/volume tolerances for GU and SD toxicity after EBRT alone and after BT with or without EBRT were scarce and support further research including data-sharing initiatives to untangle the dose/volume relationships for these symptoms.

Re-irradiation volumetric modulated arc therapy optimization based on cumulative biologically effective dose objectives

The objective of this note is to introduce a clinical tool that generates ideal base plan dose distributions to enable re‐irradiation volumetric modulated arc therapy (VMAT) optimization based on cumulative biological effective dose objectives for specific organs at risk (OARs). The tool is demonstrated with a lung cancer case that required re‐irradiation at our clinic. First, previous treatment dose is deformed onto the retreatment computed tomography (CT) using commercial software. Then, the in‐house Matlab tool alters the deformed previous dose using radiobiological concepts on a voxel‐by‐voxel manner to generate an ideal base plan dose distribution. Ideal base plans that were generated using the in‐house Matlab tool were compatible with the Varian Eclipse™ treatment planning system. The tool enabled optimization of VMAT re‐irradiation plans using cumulative dose limits for OARs and all OAR cumulative dose objectives were met on the first optimization for the recurrent lung cancer case tested.