3D inverse treatment planning for the tandem and ovoid applicator in cervical cancer

Dosimetric Evaluation of Different Optimization Algorithms Used in Interstitial Brachytherapy of Cervical Carcinoma

BACKGROUND

Conventional optimization techniques are based on the planning approach in which positions and weights are varied to generate the desired dose distribution. Inverse planning simulated annealing (IPSA) is an advanced optimization method developed to automatically determine a suitable combination of positions to design an acceptable plan.

OBJECTIVE

In this study, three optimization techniques namely IPSA, graphical optimization (GROPT), and geometrical optimization (GOPT) methods are compared in high-dose-rate interstitial brachytherapy of cervical carcinoma.

MATERIAL AND METHODS

In this retrospective study, twenty computed tomography (CT) data sets of 10 cervical cancer patients treated with Martinez Universal Perineal Interstitial Template-based interstitial brachytherapy were studied. The treatment plans generated were optimized using the IPSA, and GOPT methods. The prescribed dose was 24 Gy in 4 fractions. Plans produced using IPSA, GrOPT, and GOPT techniques were analyzed for comparison of dosimetric parameters, including target coverage, homogeneity, conformity, and organs at risk (OAR) doses.

RESULTS

V100 values for IPSA, GrOPT and GOPT plans were 95.81±2.33%, 93.12±2.76% and 88.90±4.95%, respectively. The mean D90 values for the IPSA, GrOPT, and GOPT plans were 6.45±0.15 Gy, 6.12±0.21 Gy, and 5.85±0.57 Gy, respectively. Significantly lower doses of OAR were in the IPSA plans that were more homogeneous (HI=0.66). Conformity was comparatively higher in IPSA-based plans (CI=0.75).

CONCLUSION

IPSA plans were superior and resulted in better target coverage, homogeneity, conformity, and minimal OAR doses.

Integration of rotatable tandem applicator to conventional ovoid applicator toward complete framework of intensity modulated brachytherapy (IMBT) for cervical cancer

A new tandem applicator with tungsten shield for Ir-192 radiation source used in intra-cavitary brachytherapy (ICBT) enabled intensity modulated brachytherapy (IMBT) in cervical cancer treatment through fluence-modulation by rotating shield. Our previous work employed group-wise and element-wise sparsity constraints for plan optimization of tandem applicator to minimizes the number of activated angles and source dwell points for delivery efficiency. It, however, did not incorporate the ovoid applicators into the optimizing process, which is generally used to prevent cancer recurrence. To integrate ovoid applicators to the new tandem applicator, this work proposed a comprehensive framework that modifies 1) dose deposition matrix for inverse planning, and 2) plan optimizing algorithm. The dose deposition matrix was newly formulated by the Monte-Carlo simulated dose distribution for 10 positions of ovoid applicators, followed by combining those with tandem-associated dose deposition matrix. The plan optimizing algorithm decomposed entire elements into tandem and ovoid applicators, which were governed by different constraints adaptive to specified plan objectives. The integrated framework was compared against conventional ICBT, and IMBT with tandem only for three patients with asymmetric dose distributions. Integrated IMBT framework resulted in the most optimal plans. Including fluence-modulation by rotating-shield outperformed conventional ICBT in dose sparing to critical organs. Adopting ovoid applicators to the optimization yielded more conformal dose distribution around inferior, laterally expanded region of target volume. The resulting plans reduced D_5cc and D_2cc by 30.9% and 27.8% for critical organs over conventional ICBT, and by 20.6% and 21.5% for target volume over IMBT with tandem only.

Evaluation of dose differences between intracavitary applicators for cervical brachytherapy using knowledge-based models

PURPOSE

Currently, there is a lack of patient-specific tools to guide brachytherapy planning and applicator choice for cervical cancer. The purpose of this study is to evaluate the accuracy of organ-at-risk (OAR) dose predictions using knowledge-based intracavitary models, and the use of these models and clinical data to determine the dosimetric differences of tandem-and-ring (T&R) and tandem-and-ovoids (T&O) applicators.

MATERIALS AND METHODS

Knowledge-based models, which predict organ D_2cc, were trained on 77/75 cases and validated on 32/38 for T&R/T&O applicators. Model performance was quantified using ΔD_2cc=D_2cc,actual-D_{2cc,predicted}, with standard deviation (σ(ΔD_2cc)) representing precision. Model-predicted applicator dose differences were determined by applying T&O models to T&R cases, and vice versa, and compared to clinically-achieved D_2cc differences. Applicator differences were assessed using a Student's t-test (p < 0.05 significant).

RESULTS

Validation T&O/T&R model precision was 0.65/0.55 Gy, 0.55/0.38 Gy, and 0.43/0.60 Gy for bladder, rectum and sigmoid, respectively, and similar to training. When applying T&O/T&R models to T&R/T&O cases, bladder, rectum and sigmoid D_2cc values in EQD2 were on average 5.69/2.62 Gy, 7.31/6.15 Gy and 3.65/0.69 Gy lower for T&R, with similar HRCTV volume and coverage. Clinical data also showed lower T&R OAR doses, with mean EQD2 D_2cc deviations of 0.61 Gy, 7.96 Gy (p < 0.01) and 5.86 Gy (p < 0.01) for bladder, rectum and sigmoid.

CONCLUSIONS

Accurate knowledge-based dose prediction models were developed for two common intracavitary applicators. These models could be beneficial for standardizing and improving the quality of brachytherapy plans. Both models and clinical data suggest that significant OAR sparing can be achieved with T&R over T&O applicators, particularly for the rectum.

High-dose-rate brachytherapy using inverse planning optimization with tandem and ovoid applicators for locally advanced cervical cancer: a simulation study

The purpose of this study was to evaluate the dosimetric advantage of inverse planning optimization (IPO) in locally advanced cervical cancer using high-dose-rate (HDR) brachytherapy (BT). IPO was compared with point A plan and geometric optimization (GO). The three planning methods were evaluated using doses to the virtual organ-at-risk (OAR) and D90 (the minimum dose covering of 90% of the volume) to the virtual high-risk clinical target volume (HR-CTV) based on quantitative analysis. HR-CTV structures measuring 38.5 cm³ and 59.5 cm³ were created based on the outcomes of a European Group cohort study. The HR-CTVs were located either at the center of the cervix (cervical center model) or shifted toward the left (extensive disease model). The distances from the HR-CTVs to the rectum were set at 0 mm, 3 mm, 5 mm, 8 mm, 10 mm, and 15 mm. We analyzed 336 conditions in this study. All doses were evaluated by conversion to doses in 2-Gy fractions of conventional radiotherapy. The D90 of the HR-CTV using GO and IPO could achieve the desired dose in all conditions, providing better results than that of point A plan. IPO makes it possible to deliver a sufficient dose for the D90 of the HR-CTV while reducing OAR doses with smaller HR-CTV values (< 38.5 cm³) in HDR BT. However, taking into consideration the whole treatment planning time, IPO should be used only when the HR-CTV-to-rectum distance is ≤ 5 mm, and the use of GO cannot spare OARs.

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.

Dosimetric comparison of manual forward planning with uniform dwell times versus volume-based inverse planning in interstitial brachytherapy of cervical malignancies

Aim

Dosimetic comparison of manual forward planning(MFP) with inverse planning(IP) for interstitial brachytherapy(ISBT) in cervical carcinoma.

Background

Brachytherapy planning by MFP is more reliable but time-consuming method, whereas IP has been explored more often for its ease and rapidness. The superiority of either is yet to be established.

Methodology

Two plans were created on data sets of 24 patients of cervical carcinoma who had undergone ISBT, one by MFP with uniform dwell times and another IP on BrachyVision 13.7 planning system with a dose prescription of 600 cGy. Isodose shaper was used for improving conformity & homogeneity. Dosimetric parameters for target and organs at risk (OARs) were recorded. Conformity index (COIN), dose homogeneity index (DHI), overdose index (OI), Coverage index (CI) and dose nonuniformity ratio (DNR) were calculated.

Results

Mean high risk clinical target volume: 73.05(±20.7)cc, D90: 5.51 Gy vs. 5.6 Gy (p = 0.017), V100: 81.77 % vs. 83.74 % (p = 0.002), V150: 21.7 % vs. 24.93 % (p = 0.002), V200: 6.3 % vs. 6.4 % (p=0.75) for IP and MFP, respectively. CI: 0.81(IP) and 0.83(MFP) (p = 0.003); however, COIN was 0.79 for both plans. D2cc of OARs was statistically better with IP (bladder 54.7 % vs. 56.1 %, p = 0.03; rectum 63 % vs. 64.7 %, (p = 0.0008).

Conclusion

Both MFP and IP are equally acceptable dosimetrically. With higher dose achieved to the target, for a similar OAR dose, MFP provides greater user flexibility of dwell positions within the target as well as better optimization. Isodose shaper may be carefully used for fine tuning. Larger sample sizes and clinical correlation will better answer the superiority of one over the other.

Dosimetric comparison of graphical optimization and inverse planning simulated annealing for brachytherapy of cervical cancer

Purpose

Graphical optimization (GO) and inverse planning simulated annealing (IPSA) are the main treatment planning optimization techniques used in patients undergoing 3D brachytherapy treatment. This study aims to compare the dosimetric difference of plans optimized by GO and IPSA in cervical cancer brachytherapy.

Material and methods

21 cervical cancer patients data sets consisted of computed tomography (CT) and magnetic resonance imaging (MRI), acquired with the Fletcher applicator in situ were transferred to the Oncentra brachytherapy planning system. For each patient, the treatment plan was initially optimized with GO to reach a maximal D₉₀ tumor dose (6 Gy/fraction, 5 fractions), while keeping the dose to organs at risk (OARs) as low as possible. A second plan was then optimized with IPSA on the same CT images and data set (i.e., contours, catheters, and location of dwell points). Targets and OARs dose volume histograms and irradiation time were compared; data were analyzed with paired t-test; p value < 0.05 was considered statistically significant.

Results

The plans with both optimizations meet the clinical requirements. The mean D₉₀ of the clinical target volume was comparable for GO and IPSA. Similar values (p > 0.05) of target V₁₀₀, V₁₅₀, V₂₀₀, HI, and CI were registered for GO and IPSA optimizations. Bladder and rectum D_1cc and D_2cc obtained by GO resulted in larger values than those obtained by IPSA (p = 0.002). V₇₅ for bladder and rectum were slightly higher for IPSA, but without statistical difference (p > 0.05). The irradiation time was comparable (p > 0.05).

Conclusions

In 3D brachytherapy of cervical cancer, GO and IPSA optimizations do not present a significant difference in target dose coverage; nevertheless, IPSA may reduce the maximum dose to normal tissue when compared with GO.

Influence of dwell time homogeneity error weight parameter on treatment plan quality in inverse optimized high-dose-rate cervix brachytherapy using SagiPlan

Purpose

Restricting the gradients of dwell times between adjacent dwell positions can potentially be beneficial in reducing the probability of unwanted hot/cold spots occurring, if the planned applicators/anatomy relative positions change before or during treatment. This constraint, however, may degrade plan quality. This study, for the first time, aims to quantify the impact of modulation restriction on plan quality indices in inverse optimization for cervix high-dose-rate (HDR) brachytherapy using the BEBIG SagiPlan treatment planning system.

Material and methods

Ten cervical cancer patient plans were optimized for treatment with a BEBIG SagiNova ⁶⁰Co HDR afterloader using the min/max inverse planning method, with dwell time homogeneity error weight (DTHEW) parameter values of 0 to 10. Dwell time homogeneity and gradients as well as various plan quality indices were analyzed.

Results

For DTHEW = 0, min/max-based optimization yielded higher HR-CTV D₉₀ values than the variance-based option (p < 0.001) and was therefore selected for this study. Averaging over all patients, selecting non-zero DTHEWs resulted in a general increase in dwell time homogeneity and decrease in mean and maximum adjacent dwell time gradients, especially between DTHEWs of 0 and 1. For DTHEW > 1, an increase of this parameter did not always result in more homogeneous dwell times or reduced gradients in individual patients. There was a negative correlation between DTHEW and both HR-CTV D₉₀ and V₁₀₀ (p < 0.001, r = –0.91). Increasing DTHEW also negatively affected conformity index (p < 0.001, r = –0.99). Changes in rectum and sigmoid colon D_2cc were insignificant. There was a strong positive relationship between bladder D_2cc and DTHEW (p < 0.001, r = 0.99).

Conclusions

Assuming a static geometry, statistically significant degradation of plan quality can result from restricting the dwell time homogeneity in min/max-based optimization of cervix HDR brachytherapy plans using SagiPlan. Therefore, setting DTHEW to zero is indicated for the type of patient plans considered in this study.

Error analysis of applicator position for combined internal/external radiation therapy in cervical cancer

The aim of this study was to analyze the error variation in the applicator placement during the first and second radiotherapy session for cervical cancer. We recruited 22 patients with cervical cancer treated with radiotherapy. According to the image output in the first and second CT-Sim inspection, we conducted comparative analysis of image fusion to accurately measure the errors in applicator position in the horizontal (X-), longitudinal (Y-) and vertical (Z)-axes. The calibration processing was implemented in accordance with the data error measured and the location parameters, such as the angle and depth of the applicator. Electronic portal imaging technology (EPID) was used to calibrate posture change amplitude for the extracorporeal irradiation of patients, and dynamic measurement with applicator position was used to describe the error of the parameters. Finally, the data from two measurements in CT-Sim, digital reconstruction radiography (DRR) and EPID were compared. After calibration, the mean value of error of the applicator were significantly smaller. Image registration planning for error parameter calibration of applicator position can effectively reduce the applied horizontal spatial position error in radiotherapy treatment, and improve the accuracy and effectiveness during treatment.

Determining optimal eluter design by modeling physical dose enhancement in brachytherapy

PURPOSE

In situ drug release concurrent with radiation therapy has been proposed to enhance the therapeutic ratio of permanent prostate brachytherapy. Both brachytherapy sources and brachytherapy spacers have been proposed as potential eluters to release compounds, such as nanoparticles or chemotherapeutic agents. The relative effectiveness of the approaches has not been compared yet. This work models the physical dose enhancement of implantable eluters in conjunction with brachytherapy to determine which delivery mechanism provides greatest opportunity to enhance the therapeutic ratio.

MATERIALS AND METHODS

The combined effect of implanted eluters and radioactive sources were modeled in a manner that allowed the comparison of the relative effectiveness of different types of implantable eluters over a range of parameters. Prostate geometry, source, and spacer positions were extracted from treatment plans used for 125 I permanent prostate implants. Compound concentrations were calculated using steady-state solution to the diffusion equation including an elimination term characterized by the diffusion-elimination modulus (ϕb ). Does enhancement was assumed to be dependent on compound concentration up to a saturation concentration (csat ). Equivalent uniform dose (EUD) was used as an objective to determine the optimal configuration of eluters for a range of diffusion-elimination moduli, concentrations, and number of eluters. The compound delivery vehicle that produced the greatest enhanced dose was tallied for points in parameter space mentioned to determine the conditions under whether there are situations where one approach is preferable to the other.

RESULTS

The enhanced effect of implanted eluters was calculated for prostate volumes from 14 to 45 cm3 , ϕb from 0.01 to 4 mm-1 , csat from 0.05 to 7.5 times the steady-state compound concentration released from the surface of the eluter. The number of used eluters (ne ) was simulated from 10 to 60 eluters. For the region of (csat , Φ)-space that results in a large fraction of the gland being maximally sensitized, compound eluting spacers or sources produce equal increase in EUD. In the majority of the remaining (csat , Φ)-space, eluting spacers result in a greater EUD than sources even where sources often produce greater maximal physical dose enhancement. Placing eluting implants in planned locations throughout the prostate results in even greater enhancement than using only source or spacer locations.

CONCLUSIONS

Eluting brachytherapy spacers offer an opportunity to increase EUD during the routine brachytherapy process. Incorporating additional needle placements permits compound eluting spacer placement independent of source placement and thereby allowing a further increase in the therapeutic ratio. Additional work is needed to understand the in vivo spatial distribution of compound around eluters, and to incorporate time dependence of both compound release and radiation dose.

Comparison of graphical optimization or IPSA for improving brachytheraphy plans associated with inadequate target coverage for cervical cancer

Many studies have reported that inverse planning by simulated annealing (IPSA) can improve the quality of brachytherapy plans, and we wanted to examine whether IPSA could improve cervical cancer brachytherapy plans giving D₉₀ < 6 Gy (with 7 Gy per fraction) at our institution. Various IPSA plans involving the tandem and ovoid applicators were developed for 30 consecutive cervical cancer patients on the basis of computed tomography: IPSA1, with a constraint on the maximum dose in the target volume; IPSA1-0, identical to IPSA1 but without a dwell-time deviation constraint; IPSA2, without a constraint on the maximum dose; and IPSA2-0, identical to IPSA2 but without a dwell-time deviation constraint. IPSA2 achieved similar results as graphical optimization, and none of the other IPSA plans was significantly better than graphical optimization. Therefore, other approaches, such as combining interstitial and intracavitary brachytherapy, may be more appropriate for improving the quality of brachytherapy plans associated with inadequate target coverage.

Fast dose optimization for rotating shield brachytherapy

PURPOSE

To provide a fast computational method, based on the proximal graph solver (POGS) - A convex optimization solver using the alternating direction method of multipliers (ADMM), for calculating an optimal treatment plan in rotating shield brachytherapy (RSBT). RSBT treatment planning has more degrees of freedom than conventional high-dose-rate brachytherapy due to the addition of emission direction, and this necessitates a fast optimization technique to enable clinical usage.

METHODS

The multi-helix RSBT (H-RSBT) delivery technique was investigated for five representative cervical cancer patients. Treatment plans were generated for all patients using the POGS method and the commercially available solver IBM ILOG CPLEX. The rectum, bladder, sigmoid colon, high-risk clinical target volume (HR-CTV), and HR-CTV boundary were the structures included in our optimization, which applied an asymmetric dose-volume optimization with smoothness control. Dose calculation resolution was 1 × 1 × 3 mm³ for all cases. The H-RSBT applicator had 6 helices, with 33.3 mm of translation along the applicator per helical rotation and 1.7 mm spacing between dwell positions, yielding 17.5° emission angle spacing per 5 mm along the applicator.

RESULTS

For each patient, HR-CTV D₉₀ , HR-CTV D₁₀₀ , rectum D_2cc , sigmoid D_2cc , and bladder D_2cc matched within 1% for CPLEX and POGS methods. Also, similar EQD2 values between CPLEX and POGS methods were obtained. POGS was around 18 times faster than CPLEX. For all patients, total optimization times were 32.1-65.4 s for CPLEX and 2.1-3.9 s for POGS.

CONCLUSIONS

POGS reduced treatment plan optimization time approximately 18 times for RSBT with similar HR-CTV D₉₀ , organ at risk (OAR) D_2cc values, and EQD2 values compared to CPLEX, which is significant progress toward clinical translation of RSBT.

Model assessment of individual tumor control rate and adverse effects in comparing locally advanced cervical cancer treatment using intracavitary with and without interstitial brachytherapy

Purpose

This study assessed the modeled probability of tumor control and organ at risk toxicities in locally advanced cervical cancer in patients treated by external beam radiation plus brachytherapy using intracavitary combined with interstitial brachytherapy (IC/IS) vs. intracavitary brachytherapy (IC) alone.

Material and methods

Twenty cervical cancer patients with a mean HR-CTV volume of 47.4 cm³ and a mean width of 54 mm were planned with both IC/IS and IC brachytherapy alone. A probit model was utilized to model 3-year (3-yr) local control rate (LC), 3-yr cancer specific survival rate (CSS), and the adverse effect (AE) of the organ at risk by using a modeled data set from multiple institutions. Modeling results were used to estimate the LC, CSS, and AE of the treatments in this study.

Results

Using the IC/IS technique, an EQD₂ increase of 12.3 Gy to D₉₀ (from 76.1 Gy to 88.3 Gy) of HR-CTV is expected to increase 3-yr LC and 3-yr CSS by 12.5%, and 11.0%, respectively. Comparing IC/IS to IC alone, the expected G2⁺ AE were 7.7% vs. 7.9% for the bladder, and 5.9% vs. 6.8% for the rectum.

Conclusions

The IC/IS technique improved dose coverage to the HR-CTV without significantly increasing dose to 2 cm³ of the organ at risk (OAR) surrounding it. With different regimens of EBRT combined with BT, IC/IS can be used to increase the probability of LC and CSS, or decrease the risk of AE.

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.

Comparison of manual and inverse optimisation techniques in high dose rate intracavitary brachytherapy of cervical cancer: A dosimetric study

AIMS AND OBJECTIVES

To compare dosimetrically the manual optimisation with IPSA using dose volume histograms (DVH) among patients treated for carcinoma of cervix with intracavitary brachytherapy.

BACKGROUND

With the advent of advanced imaging modalities, there has been a shift from conventional X-ray based planning to three-dimensional planning. Manual optimisation is widely used across various institutions but it is time consuming and operator dependant. Inverse planning simulated annealing (IPSA) is now available in various brachytherapy planning systems. But there is a paucity of studies comparing manual optimisation and IPSA in treatment of carcinoma cervix with intracavitary brachytherapy and hence this study.

MATERIALS AND METHODS

Fifteen consecutive patients treated between December 2013 and March 2014 with intracavitary brachytherapy for carcinoma of cervix were selected for this study. All patients were initially treated with external beam radiotherapy followed by intracavitary brachytherapy. The DVH was evaluated and compared between manually optimised plans and IPSA in the same set of patients.

RESULTS

There was a significant improvement in the HRCTV coverage, mean V100 of 87.75% and 82.37% (p = 0.001) and conformity index 0.67 and 0.6 (p = 0.007) for plans generated using IPSA and manual optimisation, respectively. Homogeneity index and dose to the OARs remained similar between the two groups.

CONCLUSION

The use of inverse planning in intracavitary brachytherapy of cervix has shown a significant improvement in the target volume coverage when compared with manual planning.

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.

Dosimetric evaluation of manually and inversely optimized treatment planning for high dose rate brachytherapy of cervical cancer

BACKGROUND

To compare five inverse treatment planning methods with the conventional manually optimized planning approach for brachytherapy of cervical cancer with respect to dosimetric parameters.

MATERIAL AND METHODS

Eighteen cervical cancer patients treated with magnetic resonance imaging (MRI)-guided high dose rate (HDR) brachytherapy were included in this study. Six plans were created for each of the 4 HDR brachytherapy fractions for each patient: 1 manually optimized and 5 inversely planned. Three of these were based on inverse planning simulated annealing (IPSA) with and without extra constraints on maximum doses of the target volume, and different constraints on doses to the organs at risk (OARs). In addition there were two plans based on dose to target surface points. The resulting dose-volume histograms were analyzed and compared from the dosimetric point of view by quantifying specific dosimetric parameters, such as clinical target volume (CTV) D₉₀, CTV D₁₀₀, conformal index (COIN), and D₂cm3 for rectum, bladder and the sigmoid colon.

RESULTS

Manual optimization led to a mean target coverage of 78.3% compared to 87.5%, 91.7% and 82.5% with the three IPSA approaches (p < 0.001). Similar COIN values for manual and inverse optimization were found. The manual optimization led to better results with respect to the dose to the OARs expressed as D₂cm3. Overall, the best results were obtained with manual optimization and IPSA plans with volumetric constraints including maximum doses to the target volume.

CONCLUSIONS

Dosimetric evaluation of manual and inverse optimization approaches is indicating the potential of IPSA for brachytherapy of cervical cancer. IPSA with constraints of maximum doses to the target volume is closer related to manual optimization than plans with constraints only to minimum dose to the target volume and maximum doses to OARs. IPSA plans with proper constraints performed better than those based on dose to target surface points and manually optimized plans.

Inverse planning for combination of intracavitary and interstitial brachytherapy for locally advanced cervical cancer

The main purpose of this study was to compare three different treatment plans for locally advanced cervical cancer: (i) the inverse-planning simulated annealing (IPSA) plan for combination brachytherapy (BT) of interstitial and intracavitary brachytherapy, (ii) manual optimization based on the Manchester system for combination-BT, and (iii) the conventional Manchester system using only tandem and ovoids. This was a retrospective study of 25 consecutive implants. The high-risk clinical target volume (HR-CTV) and organs at risk were defined according to the GEC-ESTRO Working Group definitions. A dose of 6 Gy was prescribed. The uniform cost function for dose constraints was applied to all IPSA-generated plans. The coverage of the HR-CTV by IPSA for combination-BT was equivalent to that of manual optimization, and was better than that of the Manchester system using only tandem and ovoids. The mean V100 achieved by IPSA for combination-BT, manual optimization and Manchester was 96 ± 3.7%, 95 ± 5.5% and 80 ± 13.4%, respectively. The mean D100 was 483 ± 80, 487 ± 97 and 335 ± 119 cGy, respectively. The mean D90 was 677 ± 61, 681 ± 88 and 513 ± 150 cGy, respectively. IPSA resulted in significant reductions of the doses to the rectum (IPSA D2cm(3): 408 ± 71 cGy vs manual optimization D2cm(3): 485 ± 105 cGy; P = 0.03) and the bladder (IPSA D2cm(3): 452 ± 60 cGy vs manual optimization D2cm(3): 583 ± 113 cGy; P < 0.0001). In conclusion, combination-BT achieved better tumor coverage, and plans using IPSA provided significant sparing of normal tissues without compromising CTV coverage.

The departmental impact of magnetic resonance imaging in the management of cervical cancer brachytherapy: a discussion paper

This discussion paper will explore the impact of using Magnetic Resonance Imaging (MRI) in the treatment of cervical cancer with brachytherapy. It is written from the perspective of current departmental practice in the UK and aims to highlight the issues associated with using MRI as a tool for image guided brachytherapy planning. To support the discussion, a literature review was undertaken focussing specifically on the use of MRI in brachytherapy treatment planning for cervical cancer. Results from planning case studies and clinical series were analysed, and the literature showed that image guided brachytherapy treatment planning is a promising development. MRI assisted planning could theoretically be implemented in centres that have access to a MRI scanner. Alternative brachytherapy technologies (e.g. Computed Tomography (CT) assisted planning), and alternative radiotherapy modalities (e.g. an external beam radiotherapy boost), were not found to be superior in effect or of implementation. Although MRI shows great promise, the evidence base for MRI in brachytherapy planning for cervical cancer is currently limited and therefore careful implementation and evaluation is required. It is suggested by the authors that new methods of working are devised to ensure consistency and quality in implementation and delivery, and that outcomes are measured and audited to evaluate efficacy.

Preliminary results of MRI-guided brachytherapy in cervical carcinoma: the Chiangmai University experience

This study was performed to evaluate the feasibility of magnetic resonance imaging (MRI) in the treatment planning of image-guided brachytherapy for cervical carcinoma. Seventeen consecutive patients with locally advanced cervical cancer were enrolled in the study. Fifteen patients could be evaluated. When comparing the tumor at diagnosis (GTV-Dx) and the tumor at the first brachytherapy (GTV-BT), 11 of 15 patients showed a tumor regression of more than 80% while only four patients had less than 80% tumor regression. The mean D90 of HR-CTV and the calculated D2cc of the bladder, rectum, and sigmoid were 99.2 ± 11 Gy, 87.7 ± 5.7 Gy, 68.4 ± 5.4 Gy and 70.3 ± 6.8 Gy, respectively. No grade 3-4 acute toxicity was observed. The MRI can be a valuable tool for evaluating tumor response after external beam radiotherapy (EBRT) and is very helpful for prognosis prediction by residual GTV evaluation. Furthermore, MRI-guided brachytherapy allowed us to optimize the dose for both the target volumes and the OARs.

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.

Variation of treatment planning parameters (D90 HR-CTV, D 2cc for OAR) for cervical cancer tandem ring brachytherapy in a multicentre setting: comparison of standard planning and 3D image guided optimisation based on a joint protocol for dose-volume constraints

PURPOSE

To perform a qualitative and quantitative comparison of different treatment planning methods used in different centres for MRI-based brachytherapy (BT) of cervical cancer.

MATERIALS AND METHODS

Two representative patients with advanced cervical cancer (1 "limited volume case"; 1 "extensive volume case") were planned for brachytherapy (BT) with a tandem-ring applicator by six different centres. During a workshop all centres produced an institutional standard plan and an MRI-based adaptive treatment plan for each case. Optimisation was based on the fractionation schedule (HDR, PDR) and method according to the institutional protocol.

RESULTS

The loading pattern, dwell times, shape of the point A isodose varied considerably between institutional standard plans, as did dose-volume parameters for high risk CTV (HR-CTV) and also for the D(2cc) for OAR, violating the dose-volume constraints in many situations. During optimisation, the centres stayed as close as possible to the standard loading pattern and dwell times. The dose distributions and dose-volume parameters between the plans from the different centres became much more comparable after optimisation. The prescribed dose to the HR-CTV could be achieved in the limited volume case by all centres, in the extensive case only if additional needles were applied.

CONCLUSION

Treatment planning for gynaecologic brachytherapy based on different traditions shows less variation in regard to target coverage and OAR dose, when 3D image-based optimisation is performed with a uniform prescription protocol.

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.

Does the entire uterus need to be treated in cancer of the cervix? Role of adaptive brachytherapy

PURPOSE

To evaluate local control and toxicity by use of a method of adaptive cervical brachytherapy (ACB).

METHODS AND MATERIALS

From 1998 to 2008, we identified 65 cervical cancer patients with FIGO (International Federation of Gynecology and Obstetrics) Stage IB1-IVA disease who received definitive external beam radiation therapy and high-dose rate brachytherapy with tandem and ovoid applicators. As tumors regressed, 45 of 65 patients had the tandem source retracted from the uterine fundus at successive brachytherapy insertions, thus decreasing the number of (192)Ir dwell positions. Tests of trend and Fisher's exact test were used to identify the effect of ACB on disease control and toxicity. Kaplan-Meier analyses were performed to evaluate disease control and late complications.

RESULTS

The median follow-up was 24.5 months. Of the patients, 92% received chemotherapy. The 3-year overall survival, 3-year disease-free survival, 3-year distant metastasis-free survival, and local control rates were 67%, 76%, 79%, and 97%, respectively. There was only 1 isolated local failure, and there were no local failures beyond 1 year. Distant failure was involved in 93% of recurrences. No significant trend was identified regarding the extent of retraction of the tandem source start position with either failure or toxicity. Acute and actuarial 3-year late Grade 3 toxicity or greater occurred in 24.6% and 17% of patients, respectively.

CONCLUSIONS

ACB determined by clinical response yielded excellent local control rates. These data indicate that ACB may be useful in decreasing late toxicities from high-dose rate brachytherapy. With the advent of three-dimensional image-guided brachytherapy, additional methods to adapt treatment technique to changes in tumor volume warrant investigation.

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.

Inverse Planning Approach for 3-D MRI-Based Pulse-Dose Rate Intracavitary Brachytherapy in Cervix Cancer

PURPOSE

The purpose of this study was to evaluate the inverse planning simulated annealing (IPSA) software for the optimization of dose distribution in patients with cervix carcinoma treated with MRI-based pulsed-dose rate intracavitary brachytherapy.

METHODS AND MATERIALS

Thirty patients treated with a technique using a customized vaginal mold were selected. Dose-volume parameters obtained using the IPSA method were compared with the classic manual optimization method (MOM). Target volumes and organs at risk were delineated according to the Gynecological Brachytherapy Group/European Society for Therapeutic Radiology and Oncology recommendations. Because the pulsed dose rate program was based on clinical experience with low dose rate, dwell time values were required to be as homogeneous as possible. To achieve this goal, different modifications of the IPSA program were applied.

RESULTS

The first dose distribution calculated by the IPSA algorithm proposed a heterogeneous distribution of dwell time positions. The mean D90, D100, and V100 calculated with both methods did not differ significantly when the constraints were applied. For the bladder, doses calculated at the ICRU reference point derived from the MOM differed significantly from the doses calculated by the IPSA method (mean, 58.4 vs. 55 Gy respectively; p = 0.0001). For the rectum, the doses calculated at the ICRU reference point were also significantly lower with the IPSA method.

CONCLUSIONS

The inverse planning method provided fast and automatic solutions for the optimization of dose distribution. However, the straightforward use of IPSA generated significant heterogeneity in dwell time values. Caution is therefore recommended in the use of inverse optimization tools with clinical relevance study of new dosimetric rules.

Optimization of HDR brachytherapy dose distributions using linear programming with penalty costs

Prostate cancer is increasingly treated with high-dose-rate (HDR) brachytherapy, a type of radiotherapy in which a radioactive source is guided through catheters temporarily implanted in the prostate. Clinicians must set dwell times for the source inside the catheters so the resulting dose distribution minimizes deviation from dose prescriptions that conform to patient-specific anatomy. The primary contribution of this paper is to take the well-established dwell times optimization problem defined by Inverse Planning by Simulated Annealing (IPSA) developed at UCSF and exactly formulate it as a linear programming (LP) problem. Because LP problems can be solved exactly and deterministically, this formulation provides strong performance guarantees: one can rapidly find the dwell times solution that globally minimizes IPSA's objective function for any patient case and clinical criteria parameters. For a sample of 20 prostates with volume ranging from 23 to 103 cc, the new LP method optimized dwell times in less than 15 s per case on a standard PC. The dwell times solutions currently being obtained clinically using simulated annealing (SA), a probabilistic method, were quantitatively compared to the mathematically optimal solutions obtained using the LP method. The LP method resulted in significantly improved objective function values compared to SA (P = 1.54 x 10(-7)), but none of the dosimetric indices indicated a statistically significant difference (P < 0.01). The results indicate that solutions generated by the current version of IPSA are clinically equivalent to the mathematically optimal solutions.

Class solution for inversely planned permanent prostate implants to mimic an experienced dosimetrist

The purpose of this paper is to present a method for the selection of inverse planning parameters and to establish a set of inverse planning parameters (class solution) for the inverse planning included in a commercial permanent prostate implant treatment planning system. The manual planning of more than 750 patients since 1996 led to the establishment of general treatment planning rules. A class solution is tuned to fulfill the treatment planning rules and generate equivalent implants. For ten patients, the inverse planning is compared with manual planning performed by our experienced physicist. The prostate volumes ranged from 17 to 51 cc and are implanted with low activity 1-125 seeds. Dosimetric indices are calculated for comparison. The inverse planning needed about 15 s for each optimization (400 000 iterations on a 2.5 GHz PC). In comparison, the physicist needed about 20 min to perform each manual plan. A class solution is found that consistently produces dosimetric indices equivalent or better than the manual planning. Moreover, even with strict seed placement rules, the inverse planning can produce adequate prostate dose coverage and organ at risk protection. The inverse planning avoids implant with seeds outside of the prostate and too close to the urethra. It also avoids needles with only one seed and needles with three consecutive seeds. This reduces the risk of complication due to seed misplacement and edema. The inverse planning also uses a smaller number of needles, reducing the cause of trauma. The quality of the treatment plans is independent of the gland size and shape. A class solution is established that consistently and rapidly produces equivalent dosimetric indices as manual planning while respecting severe seed placement rules. The class solution can be used as a starting point for every patient, dramatically reducing the time needed to plan individual patient treatments. The class solution works with inverse preplanning, intraoperative inverse preplanning, and intraoperative real-time planning. This technology is not intended to replace the physicist but to accelerate the planning process, making intraoperative treatment planning more effective.