Image registration, contour propagation and dose accumulation of external beam and brachytherapy in gynecological radiotherapy

An offline adaptive planning method based on delivered accumulated dose for brachytherapy in cervical cancer

Background and purpose

In current clinical practice, independent treatment plan optimization for each fraction of brachytherapy might not be able to fully leverage the dosimetric advantage of the cervical cancer radiotherapy combining external beam radiotherapy (EBRT) and brachytherapy (BT). This study proposed an offline adaptive planning method based on accumulated dose for BT, aiming to improve the total dose distribution of the combined radiotherapy.

Methods and materials

This study retrospectively reviewed nine cervical cancer patients treated with EBRT followed by high-dose-rate BT. For each BT fraction, we used a multi-metric deformable image registration method to accumulate the dose distributions of previously delivered EBRT and BT. The accumulated dose distribution was then imported into a customized commercial BT treatment planning system as a background in the adaptive dose optimization. Main dosimetric parameters of the target and organs at risk (OARs) were compared between the adaptive BT (ABT) and conventional BT (CBT) planning methods.

Results

For approximately 70 % of the BT fractions, the ABT plans have lower D2cc to the bladder or rectum compared with the CBT plans. In terms of total dose evaluation, the ABT planning method resulted in a decrease in mean values of D2cc, V60 and V50 for the bladder (-1.9 ± 2.0 GyEDQ2, -1.2 ± 1.2 %, and -0.9 ± 1.1 %) and rectum (-2.1 ± 1.8 GyEQD2, -1.2 ± 1.2 %, and -1.4 ± 1.3 %).

Conclusion

The offline adaptive planning method could help decrease the doses to OARs and improve the total dose distribution of combined radiotherapy, showing promising prospects for clinical use.

Physics-based data augmentation for improved training of cone-beam computed tomography auto-segmentation of the female pelvis

Background and Purpose

Labeling cone-beam computed tomography (CBCT) images is challenging due to poor image quality. Training auto-segmentation models without labelled data often involves deep-learning to generate synthetic CBCTs (sCBCT) from planning CTs (pCT), which can result in anatomical mismatches and inaccurate labels. To prevent this issue, this study assesses an auto-segmentation model for female pelvic CBCT scans exclusively trained on delineated pCTs, which were transformed into sCBCT using a physics-driven approach.

Materials and Methods

To replicate CBCT noise and artefacts, a physics-driven sCBCT (Ph-sCBCT) was synthesized from pCT images using water-phantom CBCT scans. A 3D nn-UNet model was trained for auto-segmentation of cervical cancer CBCTs using Ph-sCBCT images with pCT contours. This study included female pelvic patients: 63 for training, 16 for validation and 20 each for testing on Ph-sCBCTs and clinical CBCTs. Auto-segmentations of bladder, rectum and clinical target volume (CTV) were evaluated using Dice Similarity Coefficient (DSC) and 95th percentile Hausdorff Distance (HD95). Initial evaluation occurred on Ph-sCBCTs before testing generalizability on clinical CBCTs.

Results

The model auto-segmentation performed well on Ph-sCBCT images and generalized well to clinical CBCTs, yielding median DSC's of 0.96 and 0.94 for the bladder, 0.88 and 0.81 for the rectum, and 0.89 and 0.82 for the CTV on Ph-sCBCT and clinical CBCT, respectively. Median HD95's for the CTV were 5 mm on Ph-sCBCT and 7 mm on clinical CBCT.

Conclusions

This study demonstrates the successful training of auto-segmentation model for female pelvic CBCT images, without necessarily delineating CBCTs manually.

Associations of dose to the urethra and long-term patient-reported outcomes after radiotherapy with EBRT and HDR brachytherapy boost for prostate cancer

Purpose

Implications of radiation dose exposure to the urethra on urinary morbidity after prostate radiotherapy are poorly understood, especially by long-term patient-reported outcomes (PRO). Therefore, our primary objective was to investigate associations of urethral dose and long-term patient-reported urinary morbidity after external beam radiotherapy and high-dose rate brachytherapy boost for prostate cancer.

Materials and methods

We conducted a pre-registered (https://doi.org/10.17605/OSF.IO/A6DC3) cross-sectional study at a tertiary academic center including a consecutive sample of patients being at least two years after treatment. Primary outcome measurements included urinary domains of the EPIC-26 questionnaire. Their associations with predefined urethral dose levels were assessed by univariable analyses (Pearson's correlation) and by predefined multivariable analyses (multiple regression). Sample size calculation was based on a predefined multivariable model. A p-value < 0.05 was considered statistically significant.

Results

Among 277 screened patients, 113 patients were alive, eligible, consented, and provided PRO. The median time passed since radiotherapy was 4 years. Per univariable analysis, a higher near maximum point dose of the urethra (DU0.1cc) was associated with worse urinary incontinence (r = -0.32; CI = -0.48 - -0.13; p < 0.001) and worse overall urinary function (r = -0.21; CI = -0.38 - -0.03; p = 0.02) of the respective EPIC-26 domains. Per predefined multivariable analysis, DU0.1cc and urinary incontinence remained significantly associated (B = -0.005; CI = -0.008 - -0.002; p = 0.003). These associations were only present, when very high DU0.1 cc above 137 Gy were kept in the analysis.

Conclusions

Very high urethral near point doses appear to be associated with worse long-term patient-reported urinary morbidity after radiotherapy for prostate cancer. Urethral dose should be considered in practice and future trials to potentially minimize long-term urinary morbidity.

Trial registration

The study protocol was pre-registered prior to patient accrual on the Open Science Framework (https://doi.org/10.17605/OSF.IO/A6DC3).

Inversion-recovery ultrashort-echo-time (IR-UTE) MRI-based detection of radiation dose heterogeneity in gynecologic cancer patients treated with HDR brachytherapy

PURPOSE

To evaluate the relationship between delivered radiation (RT) and post-RT inversion-recovery ultrashort-echo-time (IR-UTE) MRI signal-intensity (SI) in gynecologic cancer patients treated with high-dose-rate (HDR) brachytherapy (BT).

METHODS

Seven patients underwent whole-pelvis RT (WPRT) followed by BT to the high-risk clinical target volume (HR-CTV). MR images were acquired at three time-points; pre-RT, post-WPRT/pre-BT, and 3-6 months post-BT. Diffuse-fibrosis (FDiffuse) was imaged with a non-contrast dual-echo IR (inversion time [TI] = 60 ms) UTE research application, with image-subtraction of the later echo, only retaining the ultrashort-echo SI. Dense-fibrosis (FDense) imaging utilized single-echo Late-Gadolinium-Enhanced IR-UTE, acquired ∼ 15 min post-Gadavist injection. Resulting FDiffuse and FDense SI were normalized to the corresponding gluteal-muscle SI. Images were deformably registered between time-points based on normal tissue anatomy. The remnant tumor at both time-points was segmented using multi-parametric MRI. Contours corresponding to the 50%, 100%, 150%, and 200% isodose lines (IDLs) of the prescription BT-dose were created. Mean FDiffuse and FDense SI within (i) each IDL contour and (ii) the remnant tumor were calculated. Post-BT FDiffuse and FDense SI were correlated with prescribed BT-dose. To determine the relationship between BT-dose and IR-UTE SI, the differences in the post-BT FDense across IDLs was determined using paired t-tests with Bonferroni correction.

RESULTS

FDense was higher in regions of higher dose for 6/7 patients, with mean ± SD values of 357 ± 103% and 331 ± 97% (p = .03) in the 100% and 50% IDL, respectively. FDense was higher in regions of higher dose in the responsive regions with mean ± SD values of 380 ± 122% and 356 ± 135% (p = .03) in the 150% and 50% IDL, respectively. Within the segmented remnant tumor, an increase in prescribed dose correlated with an increase in FDense post-BT (n = 5, r = .89, p = .04). Post-BT FDiffuse inversely correlated (n = 7, r = -.83, p = .02) with prescribed BT-dose within the 100% IDL.

CONCLUSIONS

Results suggest that FDense SI 3-6 months post-BT is a sensitive measure of tissue response to heterogeneous BT radiation-dose. Future studies will validate whether FDiffuse and FDense are accurate biomarkers of fibrotic radiation response.

Improvement of accumulated dose distribution in combined cervical cancer radiotherapy with deep learning-based dose prediction

Purpose

Difficulties remain in dose optimization and evaluation of cervical cancer radiotherapy that combines external beam radiotherapy (EBRT) and brachytherapy (BT). This study estimates and improves the accumulated dose distribution of EBRT and BT with deep learning-based dose prediction.

Materials and methods

A total of 30 patients treated with combined cervical cancer radiotherapy were enrolled in this study. The dose distributions of EBRT and BT plans were accumulated using commercial deformable image registration. A ResNet-101-based deep learning model was trained to predict pixel-wise dose distributions. To test the role of the predicted accumulated dose in clinic, each EBRT plan was designed using conventional method and then redesigned referencing the predicted accumulated dose distribution. Bladder and rectum dosimetric parameters and normal tissue complication probability (NTCP) values were calculated and compared between the conventional and redesigned accumulated doses.

Results

The redesigned accumulated doses showed a decrease in mean values of V50, V60, and D2cc for the bladder (-3.02%, -1.71%, and -1.19 Gy, respectively) and rectum (-4.82%, -1.97%, and -4.13 Gy, respectively). The mean NTCP values for the bladder and rectum were also decreased by 0.02‰ and 0.98%, respectively. All values had statistically significant differences (p < 0.01), except for the bladder D2cc (p = 0.112).

Conclusion

This study realized accumulated dose prediction for combined cervical cancer radiotherapy without knowing the BT dose. The predicted dose served as a reference for EBRT treatment planning, leading to a superior accumulated dose distribution and lower NTCP values.

Comparison of image registration methods in patients with non-melanoma skin cancer treated with superficial brachytherapy

The accumulated dose from sequential treatments of metachronous non-melanoma skin cancer can be assessed using image registration, although guidelines for selecting the appropriate algorithm are lacking. This study shows the impact of rigid (RIR), deformable (DIR) and deformable structure-based (SDIR) algorithms on the skin dose. DIR increased: the maximum dose (39.2 Gy vs 9.4 Gy), the dose to 0.1 cm3 (16.4 Gy vs 7.8 Gy) and the dose to 2 cm3 (7.6 Gy vs 5.7 Gy). RIR only affected the maximum dose, which increased to 17.0 Gy. SDIR correctly translated the dose maps, as none of the parameters changed significantly.

AAPM Task Group Report 267: A joint AAPM GEC-ESTRO report on biophysical models and tools for the planning and evaluation of brachytherapy

Brachytherapy utilizes a multitude of radioactive sources and treatment techniques that often exhibit widely different spatial and temporal dose delivery patterns. Biophysical models, capable of modeling the key interacting effects of dose delivery patterns with the underlying cellular processes of the irradiated tissues, can be a potentially useful tool for elucidating the radiobiological effects of complex brachytherapy dose delivery patterns and for comparing their relative clinical effectiveness. While the biophysical models have been used largely in research settings by experts, it has also been used increasingly by clinical medical physicists over the last two decades. A good understanding of the potentials and limitations of the biophysical models and their intended use is critically important in the widespread use of these models. To facilitate meaningful and consistent use of biophysical models in brachytherapy, Task Group 267 (TG-267) was formed jointly with the American Association of Physics in Medicine (AAPM) and The Groupe Européen de Curiethérapie and the European Society for Radiotherapy & Oncology (GEC-ESTRO) to review the existing biophysical models, model parameters, and their use in selected brachytherapy modalities and to develop practice guidelines for clinical medical physicists regarding the selection, use, and interpretation of biophysical models. The report provides an overview of the clinical background and the rationale for the development of biophysical models in radiation oncology and, particularly, in brachytherapy; a summary of the results of literature review of the existing biophysical models that have been used in brachytherapy; a focused discussion of the applications of relevant biophysical models for five selected brachytherapy modalities; and the task group recommendations on the use, reporting, and implementation of biophysical models for brachytherapy treatment planning and evaluation. The report concludes with discussions on the challenges and opportunities in using biophysical models for brachytherapy and with an outlook for future developments.

QUIZ: An arbitrary volumetric point matching method for medical image registration

Rigid pre-registration involving local-global matching or other large deformation scenarios is crucial. Current popular methods rely on unsupervised learning based on grayscale similarity, but under circumstances where different poses lead to varying tissue structures, or where image quality is poor, these methods tend to exhibit instability and inaccuracies. In this study, we propose a novel method for medical image registration based on arbitrary voxel point of interest matching, called query point quizzer (QUIZ). QUIZ focuses on the correspondence between local-global matching points, specifically employing CNN for feature extraction and utilizing the Transformer architecture for global point matching queries, followed by applying average displacement for local image rigid transformation.We have validated this approach on a large deformation dataset of cervical cancer patients, with results indicating substantially smaller deviations compared to state-of-the-art methods. Remarkably, even for cross-modality subjects, it achieves results surpassing the current state-of-the-art.

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.

A feasibility study of adaptive radiation therapy for postprostatectomy prostate cancer

Postoperative prostate radiotherapy requires large planning target volume (PTV) margins to account for motion and deformation of the prostate bed. Adaptive radiation therapy (ART) can incorporate image-guidance data to personalize PTVs that maintain coverage while reducing toxicity. We present feasibility and dosimetry results of a prospective study of postprostatectomy ART. Twenty-one patients were treated with single-adaptation ART. Conventional treatments were delivered for fractions 1 to 6 and adapted plans for the remaining 27 fractions. Clinical target volumes (CTVs) and small bowel delineated on fraction 1 to 4 CBCT were used to generate adapted PTVs and planning organ-at-risk (OAR) volumes for adapted plans. PTV volume and OAR dose were compared between ART and conventional using Wilcoxon signed-rank tests. Weekly CBCT were used to assess the fraction of CTV covered by PTV, CTV D99, and small bowel D1cc. Clinical metrics were compared using a Student's t-test (p < 0.05 significant). Offline adaptive planning required 1.9 ± 0.4 days (mean ± SD). ART decreased mean adapted PTV volume 61 ± 37 cc and bladder wall D50 compared with conventional treatment (p < 0.01). The CTV was fully covered for 96% (97%) of fractions with ART (conventional). Reconstructing dose on weekly CBCT, a nonsignificant reduction in CTV D99 was observed with ART (94%) compared to conventional (96%). Reduced CTV D99 with ART was significantly correlated with large anterior-posterior rectal diameter on simulation CT. ART reduced the number of fractions exceeding our institution's small bowel D1c limit from 14% to 7%. This study has demonstrated the feasibility of offline ART for post-prostatectomy cancer. ART facilitates PTV volume reduction while maintaining reasonable CTV coverage and can reduce the dose to adjacent normal tissues.

EviGUIDE - a tool for evidence-based decision making in image-guided adaptive brachytherapy for cervical cancer

PURPOSE

To develop a novel decision-support system for radiation oncology that incorporates clinical, treatment and outcome data, as well as outcome models from a large clinical trial on magnetic resonance image-guided adaptive brachytherapy (MR-IGABT) for locally advanced cervical cancer (LACC).

METHODS

A system, called EviGUIDE, was developed that combines dosimetric information from the treatment planning system, patient and treatment characteristics, and established tumor control probability (TCP), and normal tissue complication probability (NTCP) models, to predict clinical outcome of radiotherapy treatment of LACC. Six Cox Proportional Hazards models based on data from 1341 patients of the EMBRACE-I study have been integrated. One TCP model for local tumor control, and five NTCP models for OAR morbidities.

RESULTS

EviGUIDE incorporates TCP-NTCP graphs to help users visualize the clinical impact of different treatment plans and provides feedback on achievable doses based on a large reference population. It enables holistic assessment of the interplay between multiple clinical endpoints and tumour and treatment variables. Retrospective analysis of 45 patients treated with MR-IGABT showed that there exists a sub-cohort of patients (20%) with increased risk factors, that could greatly benefit from the quantitative and visual feedback.

CONCLUSION

A novel digital concept was developed that can enhance clinical decision- making and facilitate personalized treatment. It serves as a proof of concept for a new generation of decision support systems in radiation oncology, which incorporate outcome models and high-quality reference data, and aids the dissemination of evidence-based knowledge about optimal treatment and serve as a blueprint for other sites in radiation oncology.

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.

Artificial intelligence applications in brachytherapy: A literature review

PURPOSE

Artificial intelligence (AI) has the potential to simplify and optimize various steps of the brachytherapy workflow, and this literature review aims to provide an overview of the work done in this field.

METHODS AND MATERIALS

We conducted a literature search in June 2022 on PubMed, Embase, and Cochrane for papers that proposed AI applications in brachytherapy.

RESULTS

A total of 80 papers satisfied inclusion/exclusion criteria. These papers were categorized as follows: segmentation (24), registration and image processing (6), preplanning (13), dose prediction and treatment planning (11), applicator/catheter/needle reconstruction (16), and quality assurance (10). AI techniques ranged from classical models such as support vector machines and decision tree-based learning to newer techniques such as U-Net and deep reinforcement learning, and were applied to facilitate small steps of a process (e.g., optimizing applicator selection) or even automate the entire step of the workflow (e.g., end-to-end preplanning). Many of these algorithms demonstrated human-level performance and offer significant improvements in speed.

CONCLUSIONS

AI has potential to augment, automate, and/or accelerate many steps of the brachytherapy workflow. We recommend that future studies adhere to standard reporting guidelines. We also stress the importance of using larger sample sizes and reporting results using clinically interpretable measures.

AAPM medical physics practice guideline 13.a: HDR brachytherapy, part A

The American Association of Physicists in Medicine (AAPM) is a nonprofit professional society whose primary purposes are to advance the science, education, and professional practice of medical physics. The AAPM has more than 8000 members and is the principal organization of medical physicists in the United States. The AAPM will periodically define new practice guidelines for medical physics practice to help advance the science of medical physics and to improve the quality of service to patients throughout the United States. Existing medical physics practice guidelines (MPPGs) will be reviewed for the purpose of revision or renewal, as appropriate, on their fifth anniversary or sooner. Each medical physics practice guideline represents a policy statement by the AAPM, has undergone a thorough consensus process in which it has been subjected to extensive review, and requires the approval of the Professional Council. The medical physics practice guidelines recognize that the safe and effective use of diagnostic and therapeutic radiology requires specific training, skills, and techniques, as described in each document. Reproduction or modification of the published practice guidelines and technical standards by those entities not providing these services is not authorized. The following terms are used in the AAPM practice guidelines: (1) Must and must not: Used to indicate that adherence to the recommendation is considered necessary to conform to this practice guideline. (2) Should and should not: Used to indicate a prudent practice to which exceptions may occasionally be made in appropriate circumstances. Approved by AAPM's Executive Committee April 28, 2022.

Applying Multi-Metric Deformable Image Registration for Dose Accumulation in Combined Cervical Cancer Radiotherapy

(1) Purpose: Challenges remain in dose accumulation for cervical cancer radiotherapy combined with external beam radiotherapy (EBRT) and brachytherapy (BT) as there are many large and complex organ deformations between different treatments. This study aims to improve deformable image registration (DIR) accuracy with the introduction of multi-metric objectives for dose accumulation of EBRT and BT. (2) Materials and methods: Twenty cervical cancer patients treated with EBRT (45-50 Gy/25 fractions) and high-dose-rate BT (≥20 Gy in 4 fractions) were included for DIR. The multi-metric DIR algorithm included an intensity-based metric, three contour-based metrics, and a penalty term. Nonrigid B-spine transformation was used to transform the planning CT images from EBRT to the first BT, with a six-level resolution registration strategy. To evaluate its performance, the multi-metric DIR was compared with a hybrid DIR provided by commercial software. The DIR accuracy was measured by the Dice similarity coefficient (DSC) and Hausdorff distance (HD) between deformed and reference organ contours. The accumulated maximum dose of 2 cc (D_2cc) of the bladder and rectum was calculated and compared to simply addition of D_2cc from EBRT and BT (ΔD_2cc). (3) Results: The mean DSC of all organ contours for the multi-metric DIR were significantly higher than those for the hybrid DIR (p ≤ 0.011). In total, 70% of patients had DSC > 0.8 using the multi-metric DIR, while 15% of patients had DSC > 0.8 using the commercial hybrid DIR. The mean ΔD_2cc of the bladder and rectum for the multi-metric DIR were 3.25 ± 2.29 and 3.54 ± 2.02 Gy_EQD2, respectively, whereas those for the hybrid DIR were 2.68 ± 2.56 and 2.32 ± 3.25 Gy_EQD2, respectively. The multi-metric DIR resulted in a much lower proportion of unrealistic D_2cc than the hybrid DIR (2.5% vs. 17.5%). (4) Conclusions: Compared with the commercial hybrid DIR, the introduced multi-metric DIR significantly improved the registration accuracy and resulted in a more reasonable accumulated dose distribution.

Neural network-assisted automated image registration for MRI-guided adaptive brachytherapy in cervical cancer

PURPOSE

In image-guided adaptive brachytherapy (IGABT) a quantitative evaluation of the dosimetric changes between fractions due to anatomical variations, can be implemented via rigid registration of images from subsequent fractions based on the applicator as a reference structure. With available treatment planning systems (TPS), this is a manual and time-consuming process. The aim of this retrospective study was to automate this process. A neural network (NN) was trained to predict the applicator structure from MR images. The resulting segmentation was used to automatically register MR-volumes.

MATERIAL AND METHODS

DICOM images and plans of 56 patients treated for cervical cancer with high dose-rate (HDR) brachytherapy were used in the study. A 2D and a 3D NN were trained to segment applicator structures on clinical T2-weighted MRI datasets. Different rigid registration algorithms were investigated and compared. To evaluate a fully automatic registration workflow, the NN-predicted applicator segmentations (AS) were used for rigid image registration with the best performing algorithm. The DICE coefficient and mean distance error between dwell positions (MDE) were used to evaluate segmentation and registration performance.

RESULTS

The mean DICE coefficient for the predicted AS was 0.70 ± 0.07 and 0.58 ± 0.04 for the 3D NN and 2D NN, respectively. Registration algorithms achieved MDE errors from 8.1 ± 3.7 mm (worst) to 0.7 ± 0.5 mm (best), using ground-truth AS. Using the predicted AS from the 3D NN together with the best registration algorithm, an MDE of 2.7 ± 1.4 mm was achieved.

CONCLUSION

Using a combination of deep learning models and state of the art image registration techniques has been demonstrated to be a promising solution for automatic image registration in IGABT. In combination with auto-contouring of organs at risk, the auto-registration workflow from this study could become part of an online-dosimetric interfraction evaluation workflow in the future.

Potential utility of cone-beam CT-guided adaptive radiotherapy under end-exhalation breath-hold conditions for pancreatic cancer

PURPOSE

The purpose of this study was to demonstrate the potential utility of cone-beam computed tomography (CBCT)-guided online adaptive radiotherapy (ART) under end-exhalation breath-hold (EE-BH) conditions for pancreatic cancer (PC).

METHODS

Eleven PC patients who underwent 15-fraction volumetric-modulated arc therapy under EE-BH conditions were included. Planning CT images and daily 165 CBCT images were imported into a dedicated treatment planning system. The prescription dose was set to 48 Gy in 15 fractions. The reference plan was automatically generated along with predefined clinical goals. After segmentation was completed on CBCT images, two different plans were generated: One was an adapted (ADP) plan in which re-optimization was performed on the anatomy of the day, and the other was a scheduled (SCH) plan, which was the same as the reference plan. The dose distributions calculated using the synthetic CT created from both planning CT and CBCT were compared between the two plans. Independent calculation-based quality assurance was also performed for the ADP plans, with a gamma passing rate of 3%/3 mm.

RESULTS

All clinical goals were successfully achieved during the reference plan generation. Of the 165 sessions, gross tumor volume D_98% and clinical target volume D_98% were higher in 100 (60.1%) and 122 (74.0%) ADP fractions. In each fraction, the V_3 Gy  < 1 cm³ of the stomach and duodenum was violated in 47 (28.5%) and 48 (29.1%), respectively, of the SCH fractions, whereas no violations were observed in the ADP fractions. There were statistically significant differences in the dose-volume indices between the SCH and ADP fractions (p < 0.05). The gamma passing rates were above 95% in all ADP fractions.

CONCLUSIONS

The CBCT-guided online ART under EE-BH conditions successfully reduced the dose to the stomach and duodenum while maintaining target coverage.

Review on Treatment Planning Systems for Cervix Brachytherapy (Interventional Radiotherapy): Some Desirable and Convenient Practical Aspects to Be Implemented from Radiation Oncologist and Medical Physics Perspectives

Intracavitary brachytherapy (BT, Interventional Radiotherapy, IRT), plays an essential role in the curative intent of locally advanced cervical cancer, for which the conventional approach involves external beam radiotherapy with concurrent chemotherapy followed by BT. This work aims to review the different methodologies used by commercially available treatment planning systems (TPSs) in exclusive magnetic resonance imaging-based (MRI) cervix BT with interstitial component treatments. Practical aspects and improvements to be implemented into the TPSs are discussed. This review is based on the clinical expertise of a group of radiation oncologists and medical physicists and on interactive demos provided by the software manufacturers. The TPS versions considered include all the new tools currently in development for future commercial releases. The specialists from the supplier companies were asked to propose solutions to some of the challenges often encountered in a clinical environment through a questionnaire. The results include not only such answers but also comments by the authors that, in their opinion, could help solve the challenges covered in these questions. This study summarizes the possibilities offered nowadays by commercial TPSs, highlighting the absence of some useful tools that would notably improve the planning of MR-based interstitial component cervix brachytherapy.

Technical Note: Commissioning of an ultrasound-compatible surrogate vaginal cylinder for transvaginal ultrasound-based gynecologic high-dose-rate brachytherapy

PURPOSE

To provide a comprehensive set of commissioning tests for clinical implementation of three-dimensional transvaginal ultrasound (3D TVUS) as a replacement of computed tomography (CT) for applicator reconstruction in gynecologic intracavitary high-dose-rate brachytherapy with a multi-channel vaginal cylinder.

METHODS

We introduce an ultrasound-compatible "surrogate" vaginal cylinder (SVC) for reconstruction of Elekta's CT-MR Multi Channel Applicator (MCVC) in 3D TVUS. The MCVC is digitized over the SVC in 3DUS using digital library model overlay. Consulting guidelines from various sources (CPQR, GEC-ESTRO, AAPM), we identify and describe three tests specific to commissioning the SVC: 1) verification of SVC outer dimensions, 2) source position accuracy of MCVC digitization over the SVC in 3D TVUS, and 3) MRI/US registration error.

RESULTS

The SVC outer dimensions (diameter and A-D marker locations) were well matched to the MCVC, however a 0.6 mm discrepancy in length between cylinder tips was observed. Source position accuracy was within 1 mm (tolerance recommended by CPQR) when reconstructing the MCVC in 3D TVUS. Dice similarity coefficients and target registration error for MRI/3D TVUS registration was similar to MRI/CT registration, which is the clinical standard.

CONCLUSIONS

These commissioning tests are performed using institutional equipment but provide the framework for any practitioners to repeat in their own setup, to demonstrate safe adoption of the 3D TVUS system for patient treatments. We demonstrate that MRI/US-based workflow achieves similar source position accuracy and image registration error as standard MRI/CT, which is consistent with standard tolerances. This is a critical step towards replacement of CT with US in gynecologic high-dose-rate brachytherapy treatments with the MCVC. This article is protected by copyright. All rights reserved.

Diagnostic Application and Systematic Evaluation of Image Registration Software in External Radiotherapy

Purpose: Analyze the clinical application of MIM maestro in cancer radiotherapy and evaluate the advantage of the software compare to the clinical applied tools. Materials and Methods: Potentially relevant studies published were identified through a pubmed and web of science search using words “MIM Maestro,” “Atlas,” “image registration,” “dose accumulation,” “irradiation.” Combinations of words were also searched as were bibliographies of downloaded papers in order to avoid missing relevant publications. Results: In many patients with cancer radiotherapy, multiple types of images are demanded, MIM Maestro is a multi-modality image information processing system for radiotherapy. Contour atlas and image registration among dose accumulation and individual fractions is beneficial for radiotherapy. Overall 34 papers were enrolled for analysis. The MIM appears to provide excellent clinical applications such as the function of contour altas, image fusion and registration, dose accumulation in radiotherapy compared to the other software. Conclusions: The regular optimization of radiotherapy technology and the development of image technology, improve the clinical efficiency. The current paper give a systematic review of MIM Maestro multi-modality image processing software.

Dosimetric impact of target definition in brachytherapy for cervical cancer - Computed tomography and trans rectal ultrasound versus magnetic resonance imaging

Background and Purpose

Magnetic Resonance Imaging (MRI) based target definition in cervix brachytherapy is limited by its availability, logistics and financial implications, therefore, use of computed tomography (CT) and Trans Rectal UltraSonography (TRUS) has been explored. The current study evaluated the dosimetric impact of CT + TRUS based target volumes as compared to gold standard MRI.

Methods and Materials

Images of patients (n = 21) who underwent TRUS followed by MRI and CT, were delineated with High-Risk Clinical Target Volume in CT (CTVHR-CT) and in MRI (CTVHR-MR). CTVHR-CT was drawn on CT images with TRUS assistance. For each patient, two treatment plans were made, on MRI and CT, followed by fusion and transfer of CTVHR-MR to the CT images, referred as CTVHR-MRonCT. The agreement between CTVHR-MRonCT and CTVHR-CT was evaluated for dosimetric parameters (D90, D98 and D50; Dose received by 90%, 98% and 50% of the volumes) using Bland-Altman plots, linear regression, and Pearson correlation.

Results

No statistically significant systematic difference was found between MRI and CT. Mean difference (±1.96 SD) of D90, D98 and D50 between CTVHR-MRonCT and CTVHR-CT was 2.0, 1.2 and 5.6 Gy respectively. The number of patients who have met the dose constraints of D90 > 85 Gy were 90% and 80% in MR and in CT respectively, others were in the borderline, with a minimum dose of 80 Gy. The mean ± SD dose-difference between MR and CT plans for bladder was significant (5 ± 13 Gy; p = 0.12) for D0.1cm3, while others were statistically insignificant.

Conclusion

CT + TRUS based delineation of CTVHR appear promising, provide useful information to optimally utilize for brachytherapy planning, however, MRI remains the gold standard.

Emerging technologies in brachytherapy

Brachytherapy is a mature treatment modality. The literature is abundant in terms of review articles and comprehensive books on the latest established as well as evolving clinical practices. The intent of this article is to part ways and look beyond the current state-of-the-art and review emerging technologies that are noteworthy and perhaps may drive the future innovations in the field. There are plenty of candidate topics that deserve a deeper look, of course, but with practical limits in this communicative platform, we explore four topics that perhaps is worthwhile to review in detail at this time. First, intensity modulated brachytherapy (IMBT) is reviewed. The IMBT takes advantage ofanisotropicradiation profile generated through intelligent high-density shielding designs incorporated onto sources and applicators such to achieve high quality plans. Second, emerging applications of 3D printing (i.e. additive manufacturing) in brachytherapy are reviewed. With the advent of 3D printing, interest in this technology in brachytherapy has been immense and translation swift due to their potential to tailor applicators and treatments customizable to each individual patient. This is followed by, in third, innovations in treatment planning concerning catheter placement and dwell times where new modelling approaches, solution algorithms, and technological advances are reviewed. And, fourth and lastly, applications of a new machine learning technique, called deep learning, which has the potential to improve and automate all aspects of brachytherapy workflow, are reviewed. We do not expect that all ideas and innovations reviewed in this article will ultimately reach clinic but, nonetheless, this review provides a decent glimpse of what is to come. It would be exciting to monitor as IMBT, 3D printing, novel optimization algorithms, and deep learning technologies evolve over time and translate into pilot testing and sensibly phased clinical trials, and ultimately make a difference for cancer patients. Today's fancy is tomorrow's reality. The future is bright for brachytherapy.

Dose Summation Strategies for External Beam Radiation Therapy and Brachytherapy in Gynecologic Malignancy: A Review from the NRG Oncology and NCTN Medical Physics Subcommittees

Definitive, nonsurgical management of gynecologic malignancies involves external beam radiation therapy (EBRT) and/or brachytherapy (BT). Summation of the cumulative dose is critical to assess the total biologic effective dose to targets and organs at risk. Cumulative dose calculation from EBRT and BT can be performed with or without image registration (IR) and biologic dose summation. Among these dose summation strategies, linear addition of dose-volume histogram (DVH) parameters without IR is the global standard for composite dose reporting. This approach stems from an era without image guidance and simple external beam and brachytherapy treatment approaches. With technological advances, EBRT and high-dose-rate BT have evolved to allow for volume-based treatment planning and delivery. Modern conformal therapeutic radiation involves volumetric or intensity modulated EBRT, capable of simultaneously treating multiple targets at different specified dose levels. Therefore, given the complexity of modern radiation treatment, the linear addition of DVH parameters from EBRT and high-dose-rate BT is challenging to represent the combined dose distribution. Deformable image registration (DIR) between EBRT and image guided brachytherapy (IGBT) data sets may provide a more nuanced calculation of multimodal dose accumulation. However, DIR is still nascent in this regard, and needs further development for accuracy and efficiency for clinical use. Biologic dose summation can combine physical dose maps from EBRT and each IGBT fraction, thereby generating a composite DVH from the biologic effective dose. However, accurate radiobiologic parameters are tissue-dependent and not well characterized. A combination of voxel-based DIR and biologic weighted dose maps may be the best approximation of dose accumulation but remains invalidated. The purpose of this report is to review dose summation strategies for EBRT and BT, including conventional equivalent dose in 2-Gy fractions dose summation without image registration, physical dose summation using 3-dimensional rigid IR and DIR, and biologic dose summation. We also provide general clinical workflows for IGBT with a focus on cervical cancer.

Interfraction dose deviation and catheter position in cervical interstitial and intracavitary image guided HDR brachytherapy

Interstitial and intracavitary gynecological HDR brachytherapy involve precise, localized delivery to targets with high dose gradients, sparing adjacent organs at risk (OAR). Due to the proximity of the rectum, bowel and bladder to the target, deviations in the applicator or catheter with respect to patient anatomy can significantly increase dose to OAR. The magnitude and direction of applicator and catheter migration at each fraction was assessed for template interstitial and tandem and ring (T&R) cohorts. The cohort included twelve gynecological patients with intact cervical lesions treated with external beam and brachytherapy. Pre-treatment CT images were registered to the simulation CT with respect to the target. Treatment catheter positions transformed into the planning CT coordinate system to evaluate localized catheter displacement and dose distributions calculated at each fraction. Dose was evaluated on the planning CT with planning contours and dwell locations at treatment position. Absolute deviation, depth and deflection angle for all patients were 4.6 ± 4.2 mm, -1.4 ± 4.0 mm, and 3.1 ± 2.3° respectively (n = 516 catheter positions for all treatment fractions and patients, mean ± SD). Absolute catheter deviation and deflection magnitude for interstitial treatments increased overall with each subsequent fraction with an overall increase of catheter retraction at each fraction (p < 0.005, n = 492 catheters, Kruskal-Wallis). A target EQD2 D90 reduction of 10 ± 10% and 7.7 ± 8.7% of the planned dose for interstitial and T&R cohorts respectively. There was an overall increase in bladder and rectal doses at each fraction. Catheter tracking in interstitial and intracavitary gynecological treatments with CT imaging revealed significant changes in catheter positioning with respect to the target volume. Overall deviations increased in magnitude with each subsequent fraction in the interstitial treatments. This caused patient dosimetry deviations, including target dose reduction and adjacent OAR doses changes.

Feasibility of fusing three-dimensional transabdominal and transrectal ultrasound images for comprehensive intraoperative visualization of gynecologic brachytherapy applicators

PURPOSE

In this study, we propose combining three-dimensional (3D) transrectal ultrasound (TRUS) and 3D transabdominal ultrasound (TAUS) images of gynecologic brachytherapy applicators to leverage the advantages of each imaging perspective, providing a broader field-of-view and allowing previously obscured features to be recovered. The aim of this study was to evaluate the feasibility of fusing these 3D ultrasound (US) perspectives based on the applicator geometry in a phantom prior to clinical implementation.

METHODS

In proof-of-concept experiments, 3D US images of application-specific multimodality pelvic phantoms were acquired with tandem-and-ring and tandem-and-ovoids applicators using previously validated imaging systems. Two TRUS images were acquired at different insertion depths and manually fused based on the position of the ring/ovoids to broaden the TRUS field-of-view. The phantom design allowed "abdominal thickness" to be modified to represent different body habitus and TAUS images were acquired at three thicknesses for each applicator. The merged TRUS images were then combined with TAUS images by rigidly aligning applicator components and manually refining the registration using the positions of source channels and known tandem length, as well as the ring diameter for the tandem-and-ring applicator. Combined 3D US images were manually, rigidly registered to images from a second modality (magnetic resonance (MR) imaging for the tandem-and-ring applicator and X-ray computed tomography (CT) for the tandem-and-ovoids applicator (based on applicator compatibility)) to assess alignment. Four spherical fiducials were used to calculate target registration errors (TREs), providing a metric for validating registrations, where TREs were computed using root-mean-square distances to describe the alignment of manually identified corresponding fiducials. An analysis of variance (ANOVA) was used to identify statistically significant differences (p < 0.05) between the TREs for the three abdominal thicknesses for each applicator type. As an additional indicator of geometric accuracy, the bladder was segmented in the 3D US and corresponding MR/CT images, and volumetric differences and Dice similarity coefficients (DSCs) were calculated.

RESULTS

For both applicator types, the combination of 3D TRUS with 3D TAUS images allowed image information obscured by the shadowing artifacts under single imaging perspectives to be recovered. For the tandem-and-ring applicator, the mean ± one standard deviation (SD) TREs from the images with increasing thicknesses were 1.37 ± 1.35 mm, 1.84 ± 1.22 mm, and 1.60 ± 1.00 mm. Similarly, for the tandem-and-ovoids applicator, the mean ± SD TREs from the images with increasing thicknesses were 1.37 ± 0.35 mm, 1.95 ± 0.90 mm, and 1.61 ± 0.76 mm. No statistically significant difference was detected in the TREs for the three thicknesses for either applicator type. The mean volume differences for the bladder segmentations were 3.14% and 2.33% and mean DSCs were 87.8% and 87.7% for the tandem-and-ring and tandem-and-ovoids applicators, respectively.

CONCLUSIONS

In this proof-of-concept study, we demonstrated the feasibility of fusing 3D TRUS and 3D TAUS images based on the geometry of tandem-and-ring and tandem-and-ovoids applicators. This represents a step toward an accessible and low-cost 3D imaging method for gynecologic brachytherapy, with the potential to extend this approach to other intracavitary configurations and hybrid applicators.

Radiomics Analysis of 3D Dose Distributions to Predict Toxicity of Radiotherapy for Cervical Cancer

Standard treatment for locally advanced cervical cancer (LACC) is chemoradiotherapy followed by brachytherapy. Despite radiation therapy advances, the toxicity rate remains significant. In this study, we compared the prediction of toxicity events after radiotherapy for locally advanced cervical cancer (LACC), based on either dose-volume histogram (DVH) parameters or the use of a radiomics approach applied to dose maps at the voxel level. Toxicity scores using the Common Terminology Criteria for Adverse Events (CTCAE v4), spatial dose distributions, and usual clinical predictors for the toxicity of 102 patients treated with chemoradiotherapy followed by brachytherapy for LACC were used in this study. In addition to usual DVH parameters, 91 radiomic features were extracted from rectum, bladder and vaginal 3D dose distributions, after discretization into a fixed bin width of 1 Gy. They were evaluated for predictive modelling of rectal, genitourinary (GU) and vaginal toxicities (grade ≥ 2). Logistic Normal Tissue Complication Probability (NTCP) models were derived using clinical parameters only or combinations of clinical, DVH and radiomics. For rectal acute/late toxicities, the area under the curve (AUC) using clinical parameters was 0.53/0.65, which increased to 0.66/0.63, and 0.76/0.87, with the addition of DVH or radiomics parameters, respectively. For GU acute/late toxicities, the AUC increased from 0.55/0.56 (clinical only) to 0.84/0.90 (+DVH) and 0.83/0.96 (clinical + DVH + radiomics). For vaginal acute/late toxicities, the AUC increased from 0.51/0.57 (clinical only) to 0.58/0.72 (+DVH) and 0.82/0.89 (clinical + DVH + radiomics). The predictive performance of NTCP models based on radiomics features was higher than the commonly used clinical and DVH parameters. Dosimetric radiomics analysis is a promising tool for NTCP modelling in radiotherapy.

Interstitial brachytherapy for gynecologic malignancies: Complications, toxicities, and management

From both a disease and management perspective, locally advanced gynecologic cancers present a significant challenge. Dose escalation with brachytherapy serves as a key treatment, providing conformal radiation while sparing at-risk organs. Intracavitary brachytherapy techniques have been shown to be effective, with improving tumor control and toxicity profiles with the advent of three-dimensional image planning. Despite this, the variations in tumor size, location, and pelvic anatomy may lead to suboptimal dosimetry with standard intracavitary applicators in some clinical scenarios. The addition of interstitial needles (interstitial brachytherapy (ISBT)) can improve the conformality of brachytherapy treatments by adding needles to peripheral (and central) regions of the target volume, improving the ability to escalate doses in these undercovered regions while sparing organs at risk. Interstitial brachytherapy can be delivered by intracavitary and interstitial hybrid applicators (ICBT/ISBT), perineal template (P-ISBT), or by free-hand technique. ISBT has however yet to be widely available because of concerns of complications and toxicities from this specialized treatment. However, with the increasing use of three-dimensional image-guided brachytherapy, there is an opportunity to increase the level of expertise in the gynecologic radiation oncology community with an improved understanding of the potential complications and morbidity. In this article, we review the acute and long-term toxicity in both ICBT/ISBT and P-ISBT using image-guided brachytherapy.

Quality Assurance in Modern Gynecological HDR-Brachytherapy (Interventional Radiotherapy): Clinical Considerations and Comments

Simple Summary

This is a focused review discussing quality assurance during interventional brachytherapy in gynecological cancers. This topic is very large and is usually addressed from the technical and physical sides, therefore, we decided to select “hot-spots” under this large title and discuss them from the point of view of clinicians. We hope that this concise and focused review will help clinicians in improving their quality assurance protocols and draw attention to the discussed issues.

Abstract

The use of brachytherapy (interventional radiotherapy) in the treatment of gynecological cancers is a crucial element in both definitive and adjuvant settings. The recent developments in high-dose rate remote afterloaders, modern applicators, treatment-planning software, image guidance, and dose monitoring systems have led to improvement in the local control rates and in some cases improved the survival rates. The development of these highly advanced and complicated treatment modalities has been accompanied by challenges, which have made the existence of quality assurance protocols a must to ensure the integrity of the treatment process. Quality assurance aims at standardizing the technical and clinical procedures involved in the treatment of patients, which could eventually decrease the source of uncertainties whether technical (source/equipment related) or clinical. This commentary review sheds light (from a clinical point of view) on some potential sources of uncertainties associated with the use of modern brachytherapy in the treatment of gynecological cancers.

Modeling Complex Deformations of the Sigmoid Colon Between External Beam Radiation Therapy and Brachytherapy Images of Cervical Cancer

PURPOSE

In this study, we investigated registration methods for estimating the large interfractional sigmoid deformations that occur between external beam radiation therapy (EBRT) and brachytherapy (BT) for cervical cancer.

METHODS AND MATERIALS

Sixty-three patients were retrospectively analyzed. The sigmoid colon was delineated on 2 computed tomography images acquired during EBRT (without applicator) and BT (with applicator) for each patient. Five registration approaches were compared to propagate the contour of the sigmoid from BT to EBRT anatomies: rigid registration, commercial hybrid (ANAtomically CONstrained Deformation Algorithm), controlling ROI surface projection of RayStation, and the classical and constrained symmetrical thin-plate spline robust point matching (sTPS-RPM) methods. Deformation of the sigmoid due to insertion of the BT applicator was reported. Registration performance was compared by using the Dice similarity coefficient (DSC), distance to agreement, and Hausdorff distance. The 2 sTPS-RPM methods were compared by using surface triangle quality criteria between deformed surfaces. Using the deformable approaches, the BT dose of the sigmoid was deformed toward the EBRT anatomy. The displacement and discrepancy between the deformable methods to propagate the planned D1cm³ and D2cm³ of the sigmoid from BT to EBRT anatomies were reported for 55 patients.

RESULTS

Large and complex deformations of the sigmoid were observed for each patient. Rigid registration resulted in poor sigmoid alignment with a mean DSC of 0.26. Using the contour to drive the deformation, ANAtomically CONstrained Deformation Algorithm was able to slightly improve the alignment of the sigmoid with a mean DSC of 0.57. Using only the sigmoid surface as controlling ROI, the mean DSC was improved to 0.79. The classical and constrained sTPS-RPM methods provided mean DSCs of 0.95 and 0.96, respectively, with an average inverse consistency error <1 mm. The constrained sTPS-RPM provided more realistic deformations and better surface topology of the deformed sigmoids. The planned mean (range) D1cm³ and D2cm³ of the sigmoid were 13.4 Gy (1-24.1) and 12.2 Gy (1-21.5) on the BT anatomy, respectively. Using the constrained sTPS-RPM to deform the sigmoid from BT to EBRT anatomies, these hotspots had a mean (range) displacement of 27.1 mm (6.8-81).

CONCLUSIONS

Large deformations of the sigmoid were observed between the EBRT and BT anatomies, suggesting that the D1cm³ and D2cm³ of the sigmoid would unlikely to be at the same position throughout treatment. The proposed constrained sTPS-RPM seems to be the preferred approach to manage the large deformation due to BT applicator insertion. Such an approach could be used to map the EBRT dose to the BT anatomy for personalized BT planning optimization.