Inverse planning simulated annealing for magnetic resonance imaging-based intracavitary high-dose-rate brachytherapy for cervical cancer

Comparison of two inverse planning algorithms for cervical cancer brachytherapy

PURPOSE

To compare two inverse planning algorithms, the hybrid inverse planning optimization (HIPO) algorithm and the inverse planning simulated annealing (IPSA) algorithm, for cervical cancer brachytherapy and provide suggestions for their usage.

MATERIAL AND METHODS

This study consisted of 24 cervical cancer patients treated with CT image-based high-dose-rate brachytherapy using various combinations of tandem/ovoid applicator and interstitial needles. For fixed catheter configurations, plans were retrospectively optimized with two methods: IPSA and HIPO. The dosimetric parameters with respect to target coverage, localization of high dose volume (LHDV), conformal index (COIN), and sparing of organs at risk (OARs) were evaluated. A plan assessment method which combines a graphical analysis and a scoring index was used to compare the quality of two plans for each case. The characteristics of dwell time distributions of the two plans were also analyzed in detail.

RESULTS

Both IPSA and HIPO can produce clinically acceptable treatment plans. The rectum D_2cc was slightly lower for HIPO as compared to IPSA (P = 0.002). All other dosimetric parameters for targets and OARs were not significantly different between the two algorithms. The generated radar plots and scores intuitively presented the plan properties and enabled to reflect the clinical priorities for the treatment plans. Significant different characteristics were observed between the dwell time distributions generated by IPSA and HIPO.

CONCLUSIONS

Both algorithms could generate high-quality treatment plans, but their performances were slightly different in terms of each specific patient. The clinical decision on the optimal plan for each patient can be made quickly and consistently with the help of the plan assessment method. Besides, the characteristics of dwell time distribution were suggested to be taken into account during plan selection. Compared to IPSA, the dwell time distributions generated by HIPO may be closer to clinical preference.

Inverse treatment planning for an electronic brachytherapy system delivering anisotropic radiation therapy

An inverse radiation treatment planning algorithm for Sensus Healthcare's Sculptura^TM electronic brachytherapy system has been designed. The algorithm makes use of simulated annealing to optimize the conformation number (CN) of the treatment plan. The highly anisotropic dose distributions produced by the Sculptura^TM x-ray source empower the inverse treatment planning algorithm to achieve highly conformal treatment plans for a wide range of prescribed planning target volumes. Over a set of 10 datasets the algorithm achieved an average CN of 0.79 ± 0.08 and an average gamma passing rate of 0.90 ± 0.10 at 5%/5 mm. A regularization term that encouraged short treatment plans was used, and it was found that the total treatment time could be reduced by 20% with only a nominal reduction in the CN and gamma passing rate. It was also found that downsampling the voxelized volume (from 320³ to 64³ voxels) prior to optimization resulted in a 150× speedup in the optimization time (from 2 + minutes to < 1 s) without affecting the quality of the treatment plan.

Automatic optimization of treatment dosimetry to improve visual outcomes in episcleral plaque brachytherapy

PURPOSE

The aim of this study was to develop an automated dose optimization algorithm, powered by simulated annealing, for inverse planning of ocular plaque brachytherapy treatments, to improve tumor coverage and critical structure sparing for reduced treatment-related morbidity.

METHODS AND MATERIALS

A simulated annealing based algorithm was built and evaluated by finding variable seed strengths that maximized dose uniformity across tumor base for model cases of various tumor heights and plaque sizes. This problem assessed its potential to find the minima in the energy function using differential loading. The algorithm was further developed to decrease doses to critical structures while maintaining desired tumor coverage. Doses to the optic disc and fovea were compared to those using uniform seed strengths for various model cases. Finally, 10 retrospective patient cases treated with uniform seed strengths were replanned with the developed algorithm to evaluate effectiveness for clinical application.

RESULTS

The developed algorithm achieved dose reductions of up to 37.3% and 39.6% to the optic disc and fovea, respectively, compared to those using uniform seed strengths, depending on the tumor and plaque size used for the model cases. Applying the algorithm to 10 clinical scenarios resulted in dose reductions of 14% to the optic disc and fovea relative to clinical treatments performed with uniform seed strengths.

CONCLUSIONS

The developed automatic dose optimization routine was able to achieve significant dose reductions to the critical structures relative to using uniformly loaded plaques both in the model and in the clinical cases.

Commissioning of a 3D image-based treatment planning system for high-dose-rate brachytherapy of cervical cancer

The objective of this work is to present commissioning procedures to clinically implement a three-dimensional (3D), image-based, treatment-planning system (TPS) for high-dose-rate (HDR) brachytherapy (BT) for gynecological (GYN) cancer. The physical dimensions of the GYN applicators and their values in the virtual applicator library were varied by 0.4 mm of their nominal values. Reconstruction uncertainties of the titanium tandem and ovoids (T&O) were less than 0.4 mm on CT phantom studies and on average between 0.8-1.0 mm on MRI when compared with X-rays. In-house software, HDRCalculator, was developed to check HDR plan parameters such as independently verifying active tandem or cylinder probe length and ovoid or cylinder size, source calibration and treatment date, and differences between average Point A dose and prescription dose. Dose-volume histograms were validated using another independent TPS. Comprehensive procedures to commission volume optimization algorithms and process in 3D image-based planning were presented. For the difference between line and volume optimizations, the average absolute differences as a percentage were 1.4% for total reference air KERMA (TRAK) and 1.1% for Point A dose. Volume optimization consistency tests between versions resulted in average absolute differences in 0.2% for TRAK and 0.9 s (0.2%) for total treatment time. The data revealed that the optimizer should run for at least 1 min in order to avoid more than 0.6% dwell time changes. For clinical GYN T&O cases, three different volume optimization techniques (graphical optimization, pure inverse planning, and hybrid inverse optimization) were investigated by comparing them against a conventional Point A technique. End-to-end testing was performed using a T&O phantom to ensure no errors or inconsistencies occurred from imaging through to planning and delivery. The proposed commissioning procedures provide a clinically safe implementation technique for 3D image-based TPS for HDR BT for GYN cancer.

A novel two-step optimization method for tandem and ovoid high-dose-rate brachytherapy treatment for locally advanced cervical cancer

PURPOSE

To present a novel method allowing fast volumetric optimization of tandem and ovoid high-dose-rate treatments and to quantify its benefits.

METHODS AND MATERIALS

Twenty-seven CT-based treatment plans from 6 consecutive cervical cancer patients treated with four to five intracavitary tandem and ovoid insertions were used. Initial single-step optimized plans were manually optimized, approved, and delivered plans created with a goal to cover high-risk clinical target volume (HR-CTV) with D90 >90% and minimize rectum, bladder, and sigmoid D2cc. For the two-step optimized (TSO) plan, each single-step optimized plan was replanned adding a structure created from prescription isodose line to the existent physician delineated HR-CTV, rectum, bladder, and sigmoid. New, more rigorous dose-volume histogram constraints for the critical organs at risks (OARs) were used for the optimization. HR-CTV D90 and OAR D2ccs were evaluated in both plans.

RESULTS

TSO plans had consistently smaller D2ccs for all three OARs while preserving HR-CTV D90. On plans with "excellent" CTV coverage, average D90 of 96% (91-102%), sigmoid, bladder, and rectum D2cc, respectively, reduced on average by 37% (16-73%), 28% (20-47%), and 27% (15-45%). Similar reductions were obtained on plans with "good" coverage, average D90 of 93% (90-99%). For plans with "inferior" coverage, average D90 of 81%, the coverage increased to 87% with concurrent D2cc reductions of 31%, 18%, and 11% for sigmoid, bladder, and rectum, respectively.

CONCLUSIONS

The TSO can be added with minimal planning time increase but with the potential of dramatic and systematic reductions in OAR D2ccs and in some cases with concurrent increase in target dose coverage. These single-fraction modifications would be magnified over the course of four to five intracavitary insertions and may have real clinical implications in terms of decreasing both acute and late toxicities.

Review of clinical brachytherapy uncertainties: analysis guidelines of GEC-ESTRO and the AAPM

Background and purpose

A substantial reduction of uncertainties in clinical brachytherapy should result in improved outcome in terms of increased local control and reduced side effects. Types of uncertainties have to be identified, grouped, and quantified.

Methods

A detailed literature review was performed to identify uncertainty components and their relative importance to the combined overall uncertainty.

Results

Very few components (e.g., source strength and afterloader timer) are independent of clinical disease site and location of administered dose. While the influence of medium on dose calculation can be substantial for low energy sources or non-deeply seated implants, the influence of medium is of minor importance for high-energy sources in the pelvic region. The level of uncertainties due to target, organ, applicator, and/or source movement in relation to the geometry assumed for treatment planning is highly dependent on fractionation and the level of image guided adaptive treatment. Most studies to date report the results in a manner that allows no direct reproduction and further comparison with other studies. Often, no distinction is made between variations, uncertainties, and errors or mistakes. The literature review facilitated the drafting of recommendations for uniform uncertainty reporting in clinical BT, which are also provided. The recommended comprehensive uncertainty investigations are key to obtain a general impression of uncertainties, and may help to identify elements of the brachytherapy treatment process that need improvement in terms of diminishing their dosimetric uncertainties. It is recommended to present data on the analyzed parameters (distance shifts, volume changes, source or applicator position, etc.), and also their influence on absorbed dose for clinically-relevant dose parameters (e.g., target parameters such as D₉₀ or OAR doses). Publications on brachytherapy should include a statement of total dose uncertainty for the entire treatment course, taking into account the fractionation schedule and level of image guidance for adaptation.

Conclusions

This report on brachytherapy clinical uncertainties represents a working project developed by the Brachytherapy Physics Quality Assurances System (BRAPHYQS) subcommittee to the Physics Committee within GEC-ESTRO. Further, this report has been reviewed and approved by the American Association of Physicists in Medicine.

Image guided, adaptive, accelerated, high dose brachytherapy as model for advanced small volume radiotherapy

Brachytherapy has consistently provided a very conformal radiation therapy modality. Over the last two decades this has been associated with significant improvements in imaging for brachytherapy applications (prostate, gynecology), resulting in many positive advances in treatment planning, application techniques and clinical outcome. This is emphasized by the increased use of brachytherapy in Europe with gynecology as continuous basis and prostate and breast as more recently growing fields. Image guidance enables exact knowledge of the applicator together with improved visualization of tumor and target volumes as well as of organs at risk providing the basis for very individualized 3D and 4D treatment planning. In this commentary the most important recent developments in prostate, gynecological and breast brachytherapy are reviewed, with a focus on European recent and current research aiming at the definition of areas for important future research. Moreover the positive impact of GEC-ESTRO recommendations and the highlights of brachytherapy physics are discussed what altogether presents a full overview of modern image guided brachytherapy. An overview is finally provided on past and current international brachytherapy publications focusing on "Radiotherapy and Oncology". These data show tremendous increase in almost all research areas over the last three decades strongly influenced recently by translational research in regard to imaging and technology. In order to provide high level clinical evidence for future brachytherapy practice the strong need for comprehensive prospective clinical research addressing brachytherapy issues is high-lighted.

Does help structures play a role in reducing the variation of dwell time in IPSA planning for gynaecological brachytherapy application?

PURPOSE

To report our experience of dosimetric comparison of IPSA and manual plans, with a focus on the use of help structures (HS) during optimization.

MATERIAL AND METHODS

33 patients who underwent MR image-based HDR intracavitary-brachytherapy for cervix cancer based on GYN-ESTRO recommendations were selected for evaluation. Tandem/ovoid (T/O) and Vienna applicators were used. HS of diameter of 5 mm were drawn around the tandem/needles/ovoid and ring. Three plans were generated: manual optimized plan (MOPT), IPSA without help structures (IPSA_woHS) and IPSA with help structures (IPSA_wHS). Dose-volume parameters and the loading pattern were evaluated.

RESULTS

For T/O, the use of HS did not make significant impact in the dose-volume parameters and in the loading of tandem and ovoids, however steep variation was found in the individual dwell time. In case of Vienna applicator, inclusion of HS in the optimization made a significant impact in loading of needles. The percentage ratio of total time of needles to the tandem (T(N/T%)) was found to be 14±2.5, 53±9, 22±6 for MOPT, IPSA_woHS and IPSA_wHS, respectively, which implies that in IPSA_woHS the dwell time in needles were half of the dwell time in the tandem, while in MOPT the needles were loaded only in 14%, and in IPSA_wHS it was 22% of the dwell time of tandem. Inclusion of HS in the optimization has reduced the contribution of dwell time of needle in IPSA_wHS. The individual variation of dwell time was also reduced in IPSA_wHS, however drawing of HS is a time consuming procedure and may not be practical for a routine practice.

CONCLUSION

The role of HS was evaluated for IPSA for T/O and Vienna-applicator, the use of HS may be beneficial in case of combined intracavitary - interstitial approach.

A detailed dosimetric comparison between manual and inverse plans in HDR intracavitary/interstitial cervical cancer brachytherapy

PURPOSE

The purpose of this study was to compare two inverse planning algorithms for cervical cancer brachytherapy and a conventional manual treatment planning according to the MUW (Medical University of Vienna) protocol.

MATERIAL AND METHODS

For 20 patients, manually optimized, and, inversely optimized treatment plans with Hybrid Inverse treatment Planning and Optimization (HIPO) and with Inverse Planning Simulated Annealing (IPSA) were created. Dosimetric parameters, absolute volumes of normal tissue receiving reference doses, absolute loading times of tandem, ring and interstitial needles, Paddick and COIN conformity indices were evaluated.

RESULTS

HIPO was able to achieve a similar dose distribution to manual planning with the restriction of high dose regions. It reduced the loading time of needles and the overall treatment time. The values of both conformity indices were the lowest. IPSA was able to achieve acceptable dosimetric results. However, it overloaded the needles. This resulted in high dose regions located in the normal tissue. The Paddick index for the volume of two times prescribed dose was outstandingly low.

CONCLUSIONS

HIPO can produce clinically acceptable treatment plans with the elimination of high dose regions in normal tissue. Compared to IPSA, it is an inverse optimization method which takes into account current clinical experience gained from manual treatment planning.

Investigation of geometric distortions on magnetic resonance and cone beam computed tomography images used for planning and verification of high-dose rate brachytherapy cervical cancer treatment

PURPOSE

To measure the amount of geometric distortions present in the three-dimensional imaging modalities--cone beam computed tomography (CBCT) and magnetic resonance imaging (MRI)--used at University of California, San Francisco, CA, for gynecologic high dose rate brachytherapy.

METHODS AND MATERIALS

An MRI- and CT-compatible water phantom with two different sets of support structures was designed and built for this study. The support structures were used to precisely position catheters that were filled with either an MRI contrast agent or a string of radio-opaque markers. The first support structure without anatomy was built to test system-based distortions. A second structure included two types of gynecologic applicators as well as several anatomical structures, including bones and rectum to test object-induced distortions. Images were acquired with CT (for reference), kilovoltage CBCT, and MRI (1.5 T with T1- and T2-weighted images). The difference in catheter positions between the images and the CT images was analyzed.

RESULTS

For CBCT, the mean of the absolute deviations was below 1mm in all directions. The inherent uncertainty in the measurement of distortion was less than 0.5mm. MRI presented mean absolute system-based distortions between 0.6 and 1.1mm in the central region of the image and between 0.7 and 2.3mm in the outer region. Images with the applicator and anatomy in place created mean absolute distortions of 0.4, 0.8, and 0.8mm or less for CBCT, MR-T1, and MR-T2 images, respectively.

CONCLUSIONS

The distortions measured in the presence of applicators are small enough to validate the use of CBCT and 1.5 T MRI for GYN brachytherapy treatment planning and verification.

New inverse planning technology for image-guided cervical cancer brachytherapy: description and evaluation within a clinical frame

PURPOSE

To test the feasibility of a new inverse planning technology based on the Hybrid Inverse treatment Planning and Optimisation (HIPO) algorithm for image-guided cervical cancer brachytherapy in comparison to conventional manual optimisation as applied in recent clinical practice based on long-term intracavitary cervical cancer brachytherapy experience.

MATERIALS AND METHODS

The clinically applied treatment plans of 10 tandem/ring (T/R) and 10 cases with additional needles (T/R+N) planned with PLATO v14.3 were included. Standard loading patterns were manually optimised to reach an optimal coverage with 7 Gy per fraction to the High Risk CTV and to fulfil dose constraints for organs at risk. For each of these patients an inverse plan was retrospectively created with Oncentra GYN v0.9.14. Anatomy based automatic source activation was based on the topography of target and organs. The HIPO algorithm included individual gradient and modification restrictions for the T/R and needle dwell times to preserve the spatial high-dose distribution as known from the long-term clinical experience in the standard cervical cancer brachytherapy and with manual planning.

RESULTS

HIPO could achieve a better target coverage (V100) for all T/R and 7 T/R+N patients. Changes in the shape of the overdose volume (V200/400) were limited. The D(2 cc) per fraction for bladder, rectum and sigmoid colon was on average lower by 0.2 Gy, 0.4 Gy, 0.2 Gy, respectively, for T/R patients and 0.6 Gy, 0.3 Gy, 0.3 Gy for T/R+N patients (a decrease from 4.5 to 4 Gy per fraction means a total dose reduction of 5 Gy EQD2 for a 4-fraction schedule). In general the dwell times in the additional needles were lower compared to manual planning. The sparing factors were always better for HIPO plans. Additionally, in 7 T/R and 7 T/R+N patients all three D(0.1 cc), D(1 cc) and D(2 cc) for vagina wall were lower and a smaller area of vagina was covered by the reference dose in HIPO plans. Overall loading times in the tandem, the ring and the needles, as well as dose distribution, were largely preserved with adaptations performed due to specific topographical variations, in particular in lateral and caudal directions.

CONCLUSIONS

Inverse planning based on the HIPO algorithm can produce treatment plans for cervical cancer brachytherapy which are comparable to plans based on manual optimisation as applied in clinical practice. It is essential to take into account the spatial dose distribution in addition to the DVH-based constraints. The proposed inverse planning concept is feasible for improving the therapeutic ratio and limiting substantial high-dose regions around needles.

High-dose rate brachytherapy using inverse planning simulated annealing for locoregionally advanced cervical cancer: a clinical report with 2-year follow-up

PURPOSE

We present clinical outcomes of image-guided brachytherapy using inverse planning simulated annealing (IPSA) planned high-dose rate (HDR) brachytherapy boost for locoregionally advanced cervical cancer.

METHODS AND MATERIALS

From February 2004 through December 2006, 51 patients were treated at the University of California, San Francisco with HDR brachytherapy boost as part of definitive radiation for International Federation of Gynecology and Obstetrics Stage IB1 to Stage IVA cervical cancer. Of the patients, 46 received concurrent chemotherapy, 43 with cisplatin alone and 3 with cisplatin/5-fluorouracil. All patients had IPSA-planned HDR brachytherapy boost after whole-pelvis external radiation to a total tumor dose of 85 Gy or greater (for alpha/beta = 10). Toxicities are reported according to National Cancer Institute CTCAE v3.0 (Common Terminology Criteria for Adverse Events version 3.0) guidelines.

RESULTS

At a median follow-up of 24.3 months, there were no toxicities of Grade 4 or greater and the frequencies of Grade 3 acute and late toxicities were 4% and 2%, respectively. The proportion of patients having Grade 1 or 2 gastrointestinal and genitourinary acute toxicities was 48% and 52%, respectively. Low-grade late toxicities included Grade 1 or 2 vaginal, gastrointestinal, and hormonal toxicities in 31%, 18%, and 4% of patients, respectively. During the follow-up period, local recurrence developed in 2 patients, regional recurrence developed in 2, and new distant metastases developed in 15. The rates of locoregional control of disease and overall survival at 24 months were 91% and 86%, respectively.

CONCLUSIONS

Definitive radiation by use of inverse planned HDR brachytherapy boost for locoregionally advanced cervical cancer is well tolerated and achieves excellent local control of disease.