eISSN: 2081-2841
ISSN: 1689-832X
Journal of Contemporary Brachytherapy
Current Issue Archive Supplements Articles in Press Journal Information Aims and Scope Editorial Office Editorial Board Register as Author Register as Reviewer Instructions for Authors Abstracting and indexing Subscription Advertising Information Links
SCImago Journal & Country Rank

4/2019
vol. 11
 
Share:
Share:
more
 
 
Original paper

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

Bin Tang
,
Xiangyu Liu
,
Xianliang Wang
,
Shengwei Kang
,
Pei Wang
,
Jie Li
,
Lucia Clara Orlandini

J Contemp Brachytherapy 2019; 11, 4: 379–383
Online publish date: 2019/08/22
Article file
Get citation
ENW
EndNote
BIB
JabRef, Mendeley
RIS
Papers, Reference Manager, RefWorks, Zotero
AMA
APA
Chicago
Harvard
MLA
Vancouver
 
 

Purpose

Conventional radiotherapy for cervical cancer combines external beam radiotherapy (EBRT) with a boost of high-dose-rate (HDR) brachytherapy. The efficacy of a brachytherapy boost has been well recognized [1,2,3,4,5], contributing significantly to local tumor control, with higher doses accumulated within the central lesion while sparing surrounding normal tissues.

During the optimization process in intracavitary brachytherapy (ICBT) planning, once the afterloading catheters are positioned in site (i.e., the Fletcher applicator), the dose delivered can be controlled by changing the source positions within the catheters and dwell times. There is a large number of combinations of dwell positions, dwell times, and corresponding weights that lead to different target dose coverage and dose to organs at risk (OARs). The identification of the best combination is one of the challenges of the treatment planning process. Graphical optimization (GO) and inverse planning simulated annealing (IPSA) are widely available treatment planning optimization techniques used in patients undergoing HDR 3D brachytherapy treatment.

In GO, the isodose lines on the transverse, sagittal, and coronal computed tomography (CT) images are adjusted manually. Then, the dwell positions and weights of the radioactive sources along the catheters are calculated to achieve the new dose distribution. GO is quite a fast method, and it is able to improve the coverage of the target while decreasing hot spots. Moreover, the optimization is not based on the anatomy.

The inverse planning approach [6] can be defined as a method of brachytherapy treatment planning, where the clinical objectives and treatment parameters are set up and adjusted along the optimization process to achieve the desired dose distribution. CT or magnetic resonance imaging (MRI) images are used not only to define the anatomy for visual assessment and dose calculation, but also to optimize the dose distribution [7]. Therefore, they provide to a physician additional flexibility and control to shape the dose distribution. In inverse planning, the anatomical features along with the dose constraints, constitute the starting point of the dose optimization process. In IPSA optimization, we have to define dose objectives for the target and dose constraints for OARs not only at the surface, but also inside the regions. A cost function with an increasing penalty for those failing to satisfy the predefined objectives and constraints is used. Then, the cost function is continuously decreased through an iterative process with an automatic adjustment of the sources and dwell weights [8], until an optimal solution is found. The main benefit of the inverse planning approach is that all clinical requirements (e.g., dose coverage and normal tissue protection) are simultaneously and automatically taken into account in the planning process.

In brachytherapy, the application of IPSA started in prostate cancer [9,10,11] and then extended to gynecological cancer [12,13,14]. However, several studies [15,16,17] have shown that IPSA gives promising results in target coverage, OARs sparing, and planning time in cervical cancer brachytherapy, either different applicators or the absence of interstitial implants might lead to different dosimetric results. This is the case in our study, which aims to compare the Fletcher applicator-based cervical brachytherapy plans optimized with GO and IPSA in terms of target and OARs dose-volume histogram (DVH) and irradiation times, when no interstitial implants are used.

Material and methods

Twenty-one cervical cancer patients treated in our center were included in this retrospective study. Patient median age was 51 (range, 42-71), and they all presented with cervical squamous cell carcinoma with FIGO stages IIB (10 patients) and IIIB (11 patients). The patients were treated with EBRT with a prescription dose of 45 Gy in 25 fractions, and were subsequently scheduled with 5/6 fractions (6 Gy/fraction) for ICBT, based on MRI scans and gynecological examination; no interstitial needles were used due to patients’ positive responses after EBRT.

For ICBT, the patients were arranged in a supine position on a vacuum bag with a thermoplastic mask modeled on their bodies. CT images with a slice thickness of 3 mm were acquired with the Philips Brilliance Big Bore (Philips, Eindhoven, The Netherlands), from L5 to 2 cm below the vulva with a CT/MRI Fletcher applicator (Elekta AB, Stockholm, Sweden) in situ. In the same configuration, T2-weighted MRI images of 3 mm thickness were acquired on a 1.5T Archivia MRI scanner (Philips Medical Systems B.V., Eindhoven, The Netherlands) at the same day. CT and MRI data sets were imported into the Oncentra Brachy treatment planning system, version 4.3 (Elekta AB, Stockholm, Sweden) and manually registered. Following the international study on MRI-guided brachytherapy in locally advanced cervical cancer (EMBRACE study), the radiation oncologist identified and delineated the clinical target volume (CTV) as well as the OARs (bladder and rectum) on the MRI T2 images; then, they were copied to CT, based on the registration between CT and MRI images. No margins were added to the CTV and OARs, since the setup uncertainties and internal motions were minimal. The CTV prescribed dose was 6 Gy, whereas for the rectum and bladder, we strove to reach a dose as low as achievable. Considering that the OARs were exposed to the full EBRT dose of 45 Gy, physical dose constraints per fraction of 5.4 Gy for the bladder and 4.2 Gy for the rectum were calculated with the EQD2 model. An experienced medical physicist dedicated to brachytherapy treatment plans (department workload 30 to 35 treatments per day) optimized the ICBT plan, and GO technique was used as first. During GO optimization, the global-local cursor was set to the middle, and 2.5 mm spacing of dwell points of a microSelectron HDR 192Ir v2r source was applied. A second plan was then obtained using IPSA optimization, using the same parameters (i.e., contours, catheters, or location of dwell points), while setting the dwell time deviation constraint to 0.4. In Table 1, the objectives and constraints for the CTV and OARs used in the clinical practice of our center are shown. For the target dose coverage comparison, the volume receiving 100%, 150%, and 200% of prescription dose (V100,V150, and V200, respectively); the conformity index (CI) and the homogeneity index (HI) were evaluated. CI and HI were defined as:

Table 1

Dose objectives set for inverse planning simulated annealing (IPSA) optimization

ROIUsageMarginSurfaceVolume
WeightMinMaxWeight
CTVRef. target1006.0
BladderOrgan4.080
RectumOrgan4.080
SigmoidOrgan4.080

[i] CTV – clinical target volume

(1)
HI=V100V150V100 /f/fulltexts/JCB/37409/JCB-11-37409-eq1_min.jpg
(2)
CI=VCTV,refVCTVVCTV,refVref /f/fulltexts/JCB/37409/JCB-11-37409-eq2_min.jpg

where VCTV,ref is the target volume covered by the prescription dose, VCTV is the target volume, and Vref is the total volume covered by the prescription dose. For the OARs, the maximum doses received by 1 and 2 cm3 (D1cc and D2cc, respectively) and the volume receiving 75% of the prescription dose (V75) for the bladder and rectum were extracted from the DVH and used for the comparison; α/β = 3 was used in the determination of EQD2 [18]. The paired t-test was performed using SPSS19.0 (SPSS Inc., Chicago, USA), and p < 0.05 was considered statistically significant.

Results

The DVH dosimetric results obtained for the CTV and OARs for both optimizations are reported in Table 2. The datasets of the 21 patients are shown as the mean ± standard deviation (x ±s); the mean volume of CTV is 35 (range, 23-60 cc).

Table 2

CTV and OARs DVH parameters of 21 patients, obtained with GO and IPSA optimizations techniques

ParameterGOIPSAtp value
CTVV100 (%)87.46 ±3.9786.81 ±4.520.6400.529
V150 (%)
V200 (%)
53.71 ±3.04
32.27 ±3.64
54.61 ±4.13
33.61 ±3.51
0.795
1.093
0.436
0.287
HI0.38 ±0.050.37 ±0.061.0020.328
CI0.53 ±0.100.57 ±0.111.3840.182
BladderD1cc (Gy)3.35 ±1.453.12 ±1.403.5960.002*
D2cc (Gy)3.02 ±1.362.81 ±1.313.4900.002*
V75 (cc)1.04 ±1.461.56 ±3.050.8240.420
RectumD1cc (Gy)3.38 ±0.903.02 ±0.943.6820.002*
D2cc (Gy)2.98 ±0.842.61 ±0.863.7510.002*
V75 (cc)0.51 ±0.520.52 ±0.720.0920.928

* p ≤ 0.05

CTV – clinical target volume, OARs – organs at risk, DVH – dose-volume histogram, GO – graphical optimization, IPSA – inverse planning simulated annealing, HI – homogeneity index, CI – conformity index

CTV V100, V150, and V200 for both optimization techniques gave similar mean values without a statistical difference (p > 0.05); the same behavior (p > 0.05) was registered for HI and CI. The comparison of CTV V100, V150, and V200 for each patient and optimization technique is shown in Figure 1.

Fig. 1

Clinical target volume (CTV) V100, V150, and V200 obtained with graphical optimization (GO) and inverse planning simulated annealing (IPSA) optimization techniques

/f/fulltexts/JCB/37409/JCB-11-37409-g001_min.jpg

When considering the OARs, no statistical difference was observed for the bladder and rectum V75 for both optimizations, whereas IPSA plans registered the bladder and rectum D1cc and D2cc values significantly lower than those obtained with GO (p = 0.002), appearing to have a better ability to control high doses. Figure 2 shows the isodose distribution on sagittal view of both IPSA and GO plans.

Fig. 2

Isodose distribution on the sagittal view of plans optimized with graphical optimization (GO) and inverse planning simulated annealing (IPSA)

/f/fulltexts/JCB/37409/JCB-11-37409-g002_min.jpg

After normalization to a source strength of 38.830 cGy. cm2.h-1, the mean irradiation time over 21 patients for GO and IPSA was 417.2 ±75.1 s and 407.9 ±62.4 s, respectively, with no statistical difference (t = 0.632, p = 0.534); the results obtained for each patient are reported in Figure 3A. However, the dwell point weight (DPW) distribution varied for GO and IPSA, as shown in Figure 3B, for a particular patient. IPSA demonstrated a greater range of relative dwell weight than GO. In IPSA, only 17 out of 44 dwell positions were effectively used and some specific points had relatively large dwell times, whereas in GO, a smoother distribution of DPW was shown and none of the dwell points were completely turned off.

Fig. 3

The comparison of irradiation time (A) between graphical optimization (GO) and inverse planning simulated annealing (IPSA), and the dwell time distribution for a particular patient (B)

/f/fulltexts/JCB/37409/JCB-11-37409-g003_min.jpg

Discussion

To date, the Groupe Européen de Curiethérapie in the European Society for Therapeutic Radiology and Oncology (GEC-ESTRO) and the American Brachytherapy Society (ABS) have not recommended an optimization method for ICBT of cervical cancer. Two previous dosimetric studies [19,20] conducted in our center showed that ICBT optimized by GO and IPSA, respectively, can both provide good target coverage and maintain low-dose to OARs when combined with EBRT. Through further comparison, this study shows that GO and IPSA plans do not have a significant difference in target volume dose distribution, homogeneity, and conformity for ICBT of cervical cancer. Although Hsu et al. [21] compared plans of prostate cancer optimized by GO or IPSA and they found that IPSA plan was superior to GO plan in target (V100), HI, and CI, and the dose of bladder and rectum was lower in IPSA plan. The main reason for this discrepancy between Hsu’s study and ours is due to limited number of Fletcher applicator catheters used, resulting in a small optimization space of dwell position and weight. Therefore, the potential and flexibility of IPSA in our study was not as obvious as in multichannel implanting brachytherapy. Another limitation of this study is that only surface dose points were set as the objectives in IPSA optimization process, particularly Dmin for the CTV and Dmax for the OARs, and no DVH parameters were used. However, these results revealed that IPSA plans are still better at controlling the maximum dose points of the bladder and rectum.

GO and IPSA both depend on the experience and skill of planners. In GO optimization, the isodose line is manually dragged to improve dose coverage. Although GO is an intuitive and effective optimization tool, it obviously lacks reproducibility and is largely influenced by experience of a physician. The geometrical shape of the target varies with patients, and personalized adjustment is often necessary; thus, more planning time for GO may be required than for IPSA. In contrast, IPSA is not as visually intuitive as GO. However, it can provide an optimal solution that is derived by minimizing the cost function defined by given dose objectives through an iterative optimization process. Although target and surrounding normal tissues could differ between patients, because the delineations follow some guidelines (i.e., RTOG), there is a similarity in applying template sets of dose objectives as a class solution, which can increase the planning efficiency and is highly feasible in IPSA. As for irradiation time, IPSA and GO show no statistically significant difference, whereas DPW distribution indeed varies with optimization method for a particular case, which agrees with the conclusion from Hsu et al. [21]. The inhomogeneity of DPW, as shown in Figure 3B, is a IPSA problem [17], and those points with large dwell time might result in an overdose to normal tissue around targets where no explicit OARs have been contoured [16]. Setting a higher dwell time deviation constraint and drawing support structures around applicators might be helpful [15].

Conclusions

For ICBT of cervical cancer, plans optimized with GO and IPSA may result in no significant difference in target dose; however, IPSA may reduce the maximum dose to surrounding normal tissues. Nevertheless, in this study, no significant difference was found for V200 of the target between GO and IPSA. There were possibilities that IPSA optimization could result in inhomogeneous dwell time distribution, which would consequently lead to hot spots around the applicator.

Acknowledgements

This research has been financed by the Key Research and Development Project of Sichuan, no. 2018GZ0197.

Disclosure

Authors report no conflict of interest.

References

1 

Eifel PJ, Winter K, Morris M, et al. , authors. Pelvic irradiation with concurrent chemotherapy versus pelvic and para-aortic irradiation for high-risk cervical cancer: an update of radiation therapy oncology group trial (RTOG) 90-01. J Clin Oncol. 2004. 22:p. 872–880

2 

Lindegaard JC, Fokdal LU, Nielsen SK, et al. , authors. MRI-guided adaptive radiotherapy in locally advanced cervical cancer from a Nordic perspective. Acta Oncol. 2013. 52:p. 1510–1519

3 

Nomden CN, de Leeuw AA, Roesink JM, et al. , authors. Clinical outcome and dosimetric parameters of chemo-radiation including MRI guided adaptive brachytherapy with tandem-ovoid applicators for cervical cancer patients: a single institution experience. Radiother Oncol. 2013. 107:p. 69–74

4 

Gupta P, Aich RK, Deb AR , authors. Acute complications following intracavitary high-dose-rate brachytherapy in uterine cancer. J Contemp Brachytherapy. 2014. 6:p. 276–281

5 

Teke F, Yöney A, Teke M, et al. , authors. Evaluation of outcome and prognostic factors in 739 patients with uterine cervix carcinoma: a single institution experience. Contemp Oncol (Pozn). 2015. 19:p. 130–136

6 

Lessard E , author. Development and clinical introduction of an inverse planning dose optimization by simulated annealing (IPSA) for high dose rate brachytherapy. Med Phys. 2004. 31:p. 2935

7 

Tuntipumiamorn L, Lohasammakul S, Dankulchai P, et al. , authors. Comparison of impact of target delineation of computed tomography- and magnetic resonance imaging-guided brachytherapy on dose distribution in cervical cancer. J Contemp Brachytherapy. 2018. 10:p. 418–424

8 

Jamema SV, Sharma S, Mahantshetty U, et al. , authors. Comparison of IPSA with dose-point optimization and manual optimization for interstitial template brachytherapy for gynecologic cancers. Brachytherapy. 2011. 10:p. 306–312

9 

Messing EM, Zhang J, Rubens DJ, et al. , authors. Intraoperative optimized inverse planning for prostate brachytherapy: Early experience. Int J Radiat Oncol Biol Phys. 1999. 44:p. 801–808

10 

Lessard E, Pouliot J , authors. Inverse planning anatomy-based dose optimization for HDR-brachytherapy of the prostate using fast simulated annealing algorithm and dedicated objective function. Med Phys. 2001. 28:p. 773–779

11 

Kolkman-Deurloo IK, Deleye XG, Jansen PP, et al. , authors. Anatomy based inverse planning in HDR prostate brachytherapy. Radiother Oncol. 2004. 73:p. 73–77

12 

Lessard E, Hsu I, Pouliot J , authors. Inverse planning dose optimization of interstitial gyn-template high-dose-rate brachytherapy boost. Int J Radiat Oncol Biol Phys. 2001. 51:p. 329

13 

Lessard E, Hsu IC, Pouliot J , authors. Inverse planning for interstitial gynecologic template brachytherapy: truly anatomy-based planning. Int J Radiat Oncol Biol Phys. 2002. 54:p. 1243–1251

14 

Dewitt KD, Hsu ICJ, Speight J, et al. , authors. 3D inverse treatment planning for the tandem and ovoid applicator in cervical cancer. Int J Radiat Oncol Biol Phys. 2005. 63:p. 1270–1274

15 

Trnková P, Pötter R, Baltas D, et al. , authors. New inverse planning technology for image-guided cervical cancer brachytherapy: description and evaluation within a clinical frame. Radiother Oncol. 2009. 93:p. 331–340

16 

Trnková P, Baltas D, Karabis A, et al. , authors. A detailed dosimetric comparison between manual and inverse plans in HDR intracavitary/interstitial cervical cancer brachytherapy. J Contemp Brachytherapy. 2010. 2:p. 163–170

17 

Krishnamurthy D, Cunha AJ, Hsu IC, et al. , authors. 696 poster Using IPSA algorithm in dose planning for cervix HDR brachytherapy with the Utrecht applicators. first phase results. Radiother Oncol. 2011. 99:p. S278

18 

Grigsby PW, Williamson JF, Clifford Chao KS, et al. , authors. Cervical tumor control evaluated with ICRU 38 reference volumes and integrated reference air kerma. Radiother Oncol. 2001. 58:p. 19–23

19 

Wang X, Li J, Wang P, et al. , authors. Image guided radiation therapy boost in combination with high-dose-rate intracavitary brachytherapy for the treatment of cervical cancer. J Contemp Brachytherapy. 2016. 8:p. 122–127

20 

Yin G, Wang P, Lang J, et al. , authors. Dosimetric study for cervix carcinoma treatment using intensity modulated radiation therapy (IMRT) compensation based on 3D intracavitary brachytherapy technique. J Contemp Brachytherapy. 2016. 8:p. 221–232

21 

Hsu IC, Lessard E, Weinberg V, et al. , authors. Comparison of inverse planning simulated annealing and geometrical optimization for prostate high-dose-rate brachytherapy. Brachytherapy. 2004. 3:p. 147–152

Copyright: © 2019 Termedia Sp. z o. o. This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International (CC BY-NC-SA 4.0) License (http://creativecommons.org/licenses/by-nc-sa/4.0/), allowing third parties to copy and redistribute the material in any medium or format and to remix, transform, and build upon the material, provided the original work is properly cited and states its license.
 
Quick links
© 2019 Termedia Sp. z o.o. All rights reserved.
Developed by Bentus.
PayU - płatności internetowe