|Year : 2010 | Volume
| Issue : 2 | Page : 113-119
Dosimetric analysis and comparison of IMRT and HDR brachytherapy in treatment of localized prostate cancer
V Murali, P. G. G. Kurup, P Mahadev, S Mahalakshmi
Department of Radiation Oncology, Apollo Speciality Hospital, 320, Anna Salai, Chennai-600 035, India
|Date of Submission||17-Aug-2009|
|Date of Decision||25-Sep-2009|
|Date of Acceptance||30-Dec-2009|
|Date of Web Publication||3-Apr-2010|
Apollo Speciality Hospital, 320, Padma Complex, Chennai - 600 035
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Radical radiotherapy is one of the options for the management of prostate cancer. In external beam therapy, 3D conformal radiotherapy (3DCRT) and intensity modulated radiotherapy (IMRT) are the options for delivery of increased radiation dose, as vital organs are very close to the prostate and a higher dose to these structures leads to an increased toxicity. In brachytherapy, low dose rate brachytherapy with permanent implant of radioactive seeds and high dose rate brachytherapy (HDR) with remote after loaders are available. A dosimetric analysis has been made on IMRT and HDR brachytherapy plans. Ten cases from each IMRT and HDR brachytherapy have been taken for the study. The analysis includes comparison of conformity and homogeneity indices, D100, D95, D90, D80, D50, D10 and D5 of the target. For the organs at risk (OAR), namely rectum and bladder, V100, V90 and V50 are compared. In HDR brachytherapy, the doses to 1 cc and 0.1 cc of urethra have also been studied. Since a very high dose surrounds the source, the 300% dose volumes in the target and within the catheters are also studied in two plans, to estimate the actual volume of target receiving dose over 300%. This study shows that the prescribed dose covers 93 and 92% of the target volume in IMRT and HDR brachytherapy respectively. HDR brachytherapy delivers a much lesser dose to OAR, compared to the IMRT. For rectum, the V50 in IMRT is 34.0cc whilst it is 7.5cc in HDR brachytherapy. With the graphic optimization tool in HDR brachytherapy planning, the dose to urethra could be kept within 120% of the target dose. Hence it is concluded that HDR brachytherapy may be the choice of treatment for cancer of prostate in the early stage.
Keywords: Brachytherapy, conformity, intensity modulated radiotherapy, prostate
|How to cite this article:|
Murali V, Kurup P, Mahadev P, Mahalakshmi S. Dosimetric analysis and comparison of IMRT and HDR brachytherapy in treatment of localized prostate cancer. J Med Phys 2010;35:113-9
|How to cite this URL:|
Murali V, Kurup P, Mahadev P, Mahalakshmi S. Dosimetric analysis and comparison of IMRT and HDR brachytherapy in treatment of localized prostate cancer. J Med Phys [serial online] 2010 [cited 2019 Apr 18];35:113-9. Available from: http://www.jmp.org.in/text.asp?2010/35/2/113/62201
| Introduction|| |
External beam radiotherapy (EBRT) for prostate cancer, using conformal technique, was practiced over many decades. According to Brenner and Hall,  for prostate cancer, the tumor and the surrounding late responding normal tissues are likely to have similar a/b values and thus similar sensitivities to changes in fractionation. The result is that healthy tissues can only be spared by the reduction of dose. On the other hand, increasing the dose to the tumor improves local control and survival. The only way to increase the dose to the tumor, without delivering more doses to the healthy tissues, is to make the dose distribution conformal. After the advent of IMRT, EBRT has become the choice of treatment for cancer of prostate. The IMRT has its own advantages; one of which is the treatment is non-invasive. The patient doesn't need hospitalization during the course of treatment. The dose distribution could be optimized to cover the target volume and minimize the dose to rectum. Hence the target dose could be escalated.
Brachytherapy to prostate became available as early as 1911  using Radium needles. Transrectal ultrasound based (TRUS) brachytherapy using Iodine-125 (125 I) seeds came into practice in the 1980s.  After the invention of high dose rate (HDR) remote after-loading machines, in the early 80's, removable implants were done.  There are advantages with removable implants with HDR machines, like no radiation exposure to the staff involved and short course of treatment for the patient. The catheters are implanted under TRUS guidance, using template. The source positions (dwell positions) and dwell time could be optimized for better target coverage while minimizing the dose to urethra and rectum. The catheters can be positioned outside the capsule to treat the extra capsular disease also. In EBRT, the patient is immobilized for treatment set-up reproducibility and daily setup verification is done by matching either the portal images with the DRR generated from the planning system, or the cone beam CT images (CBCT) with the CT images used for planning. The intra-fraction movement of the prostate could, however, be hardly controlled. With the removable implants, this issue doesn't arise at all. The prostate is held in position by the catheters and the catheters are held in position by the template, which is sutured to the perineum.
The objective of this study is to make a dosimetric comparison based on physical dose distribution, between IMRT and HDR brachytherapy (given as mono therapy), for the treatment of cancer of the prostate
| Materials and Methods|| |
Intensity modulated radiation therapy
Ten treatment plans from each modality of treatment are selected for comparison. The treatment plans for IMRT were based on CT images acquired at 3 mm intervals. The OAR considered are, bladder, rectum and femoral heads. All the IMRT plans are done with seven coplanar beams of 6MV x-rays using the Plato Sun Rise Planning System (Nucletron, The Netherlands). The beams are placed almost at equal intervals of the gantry angles. The system uses gradient search method for the optimization, based on the user defined dose constraints and weights. The IMRT plans are MLC based, with step and shoot method. The width of the leaf at isocenter is 1 cm and the beam-let size is 1cm x 1cm. The target is the planning target volume (PTV), which includes the prostate and a 5mm margin around it. The margin is reduced to 3mm near the bladder and rectum. A dose of 76 Gy in 38 fractions is prescribed to the target. The total number of segments in a plan is approximately 70. The dose distribution of a typical plan is shown in [Figure 1].
Remote after-loading brachytherapy
A TRUS is done two days prior to the implant to map the prostate and to estimate the approximate number of catheters required. Nucletron prostate template is used for implant. The template has pre-drilled holes at 0.5 cm spacing with a square grid arrangement. Transperineal insertion of the catheters is done under TRUS. Post implant treatment plans are done on CT images acquired at 3mm intervals. The CT images are acquired roughly an hour after implantation is over. Rectum, bladder and urethra are taken as the OAR. While contouring the bladder and rectum, the outer most mucosa is contoured. For the urethra the outer surface of the Foley's is contoured. Plans are done with Iridium-192 (192 Ir) stepping source, using Plato Sun Rise Planning System. The step size for the dwell position used is 2.5mm. Dose points are defined on the periphery of the target volume at 8mm interval and optimization is made on these target dose points with a gradient factor of 0.5. A dose of 30 Gy in three fractions is prescribed to the target. One fraction is given on the day of the implant and the remaining two fractions on the following day.
Optimization in the brachytherapy module is different from the one in the inverse planning of external beam therapy module. The optimization plays with the dwell times and dwell positions. Only the dose points described on the target volume are considered for the optimization. No consideration is given for the hot spot inside the target. The objective is to achieve maximum coverage of target to the prescribed dose. No heterogeneity correction is applied and calculations are done for homogeneous tissue medium. Studies have been reported on the use of inverse planning with simulated annealing (IPSA) for HDR brachytherapy. [4,5] After the dose point optimization, graphic optimization is performed by dragging the isodose lines in such a way that the 125% and 150% isodose lines are mostly towards the boundary of the capsule and the dose to urethra doesn't exceed 125% . Sometimes 'manual dwell weights' for individual source position are adjusted to obtain the required dose distributions and make sure that dwell time is not less than one second. This is especially done when the source activity is between 9 and 10 Ci. A typical dose distribution is shown in [Figure 2].
Comparison of dosimetry
For the comparison of dosimetry or dose distribution of plans, the following criteria are considered:
For target dose comparison, the best quantitative evaluation can be done with the use of dose volume histogram (DVH). The DVH is a good tool to evaluate a complex dose distribution;  but unfortunately it does not provide information with regard to the location of the global maximum dose within the target volume. In our study, we compared D 100 , D 98 , D 90 , D 80 , D 50 , D 10 and D 5 of the target [Table 1] and [Table 2]. Since the absolute dose prescribed for the IMRT and brachytherapy are different, the prescription dose (reference dose) to target volume is normalized to 100% and the relative dose values are compared. [Figure 3] shows the comparative DVH of the target. The target dose uniformity and conformity are also evaluated and compared. The conformity index (CI), as defined in ICRU  is,
- Conformity of the prescribed dose to the target
- Homogeneity of dose within the target volume
- Dose to the organs at risk
CI= 1.00 is for an ideal case.
The conformity index is a scalar quantity and by itself cannot represent the true conformity.  Several indices have been proposed for the evaluation of dose distribution of brachytherapy implants. , Meertens et al,  have mentioned that the clinical significance and the correlations of these indices were not yet clear. Since this is a comparative study, it is decided not to use two different methods to evaluate the homogeneity of IMRT and brachytherapy plans. The Homogeneity Index (HI)  is calculated as follows
where Dp is the prescribed dose to the target. The conformity and homogeneity indices are shown in the [Table 3] and [Table 4] respectively.
In radiation therapy, we are not only concerned with the dose received by the specified target volume, but also that received by the organs around the target, which are at risk. Although DVH is a useful tool to represent the statistics of dose distribution, its applicability to hollow organs such as esophagus and rectum remains unresolved. ,,,, Moreover, in the present study, the brachytherapy calculation algorithm considers the hollow structures also as unit density tissue medium. Some authors have compared the DVH, the dose-surface histogram (DSH) and dose-wall histogram (DWH) for cylindrical and spherical models of rectum and bladder respectively. ,,, In this study, the actual volumes of rectum and bladder receiving the 100%(V 100 ), 90%(V 90 ) and 50%(V 50 ) of the target dose are compared [Table 5] and [Table 6].The DVHs of rectum and bladder are illustrated in [Figure 4] and [Figure 5]. The dose to urethra [Table 7] is studied in brachytherapy for two reasons. One, the dose per fraction is 10Gy and the other one is, in the immediate vicinity of sources, the dose is very high, which is the very nature of brachytherapy. In IMRT, as the prescribed dose per fraction is only 2Gy to the target and the maximum dose in the target doesn't exceed 115% [Table 1], the urethral dose is not of great concern.
Even though the radiobiological aspect is beyond the scope of this article, the biological equivalence between IMRT and brachytherapy is found out purely for academic interest. The biologically effective dose (BED) is calculated for both the IMRT and brachytherapy using the equation, suggested by Fowler 
where N - the no. of fractions, d - the dose per fraction in Gy and a/b the tissue specific parameter. a/b value isassumed to be 3.1Gy for the tumor.  The above equation does not take into account the tumor cell proliferation factor. The BED values calculated for IMRT (76 Gy =2 Gy X 38#) and brachytherapy (30 Gy=10.0 Gy x 3 #) are 125.0Gy and 126.7Gy respectively. This shows that the dose regimens followed are radio-biologically equivalent.
| Results and Discussion|| |
The P-values mentioned in the tables are estimated using student t-score method. From [Table 1] and [Table 2] it is seen that the target volume ranges from 70 -113 cc with a mean volume of 91 cc in IMRT. In brachytherapy the target volume ranges from 34-110cc with a mean of 71 cc. It is seen that the mean value of D 90 is almost the same in both IMRT and brachytherapy. In brachytherapy the D 100 and D 98 are slightly lower compared to IMRT. One of the reasons for this is the inability to insert catheters especially on the lateral and anterior aspects due to narrowing of the pubic arch. As the catheters are inserted through the template, there is very limited scope to manipulate the catheters to reach the target.
In this study, the mean value of D 5 in IMRT is 112%, whereas in brachytherapy it is as high as 300%. A comparative DVH of the target is shown in [Figure 3]. In this analysis we have considered D 5 and not D 2 for the maximum dose within the target; for two reasons - one, in brachytherapy, total volume of the catheters constitute about 4-5% of the volume of the target. Of course, this value changes with the number of catheters used and the volume of the prostate. The other, within the catheters the dose would be immensely high to be considered as high dose within the tissue. To establish this fact, the 300% dose within the catheters, are calculated for study nos. 7 and 10, who were implanted with 17 and 18 catheters, respectively. Those parts of the catheters, which are in the target volume alone, are marked on CT images and considered for the analysis. In study no.7, volumes of the 300% dose within the target and within the catheters are 3.12 cc and 1.53 cc respectively. Hence the net volume of target tissue is 1.59 cc. This is only 1.95% for the target volume of 81.4cc, whereas the volume of catheters that receive 300% dose is 1.87%. For the study no.10, it is found that the volume of the target tissue is 2.84% and that of the catheters is 1.63% of the target volume. In both these cases, the volume of the catheters receiving 300% dose is close to 2% of the target volume and hence considering D 2 as the maximum target dose may not be correct.
The following formula could be used to find the actual volume of the target that receives the high doses in the region of 300 to 350%, which are formed in and around the catheters and especially when the catheters are very close.
V ACT - Actual volume of target that receives a particular % dose
V PTV - Total volume of target that receives the particular % dose
ΣCTH - Total volume of catheters that receive the particular % dose
For illustration, the details are given for study no.10 only, in the [Table 8]. The study no.10 listed in [Table 1] and [Table 2] underwent brachytherapy only. For the sake of comparison and analysis, the same CT images [Figure 1] and [Figure 2] with structure set are used to plan IMRT with the usual protocol. On comparing the DVHs of OAR [Figure 4] and [Figure 5], it is seen that the dose below 20%, is more in brachytherapy than in IMRT. The reason is very clear - that in brachytherapy optimization, the dose to OAR is not considered. In graphical optimization, only the high dose regions are adjusted.
| Conclusions|| |
- The difference in target volume coverage between IMRT and brachytherapy plans is not very significant as optimization is done for the target volume in both the plans.
- The mean conformity index 0.933 in IMRT is slightly better compared to the value of 0.922 in brachytherapy. Even though the dose distribution in the brachytherapy could be optimized and improved, it still largely depends on the implant
- The homogeneity index is very good in IMRT. For the brachytherapy plans, probably new indices could be used to evaluate the homogeneity considering the very high dose within the target volume.
The temporary implant with HDR brachytherapy provides an acceptable dose distribution to the prostate in terms of coverage and a much-reduced dose to rectum and bladder, compared to the IMRT. The only concern is the homogeneity of dose within the target. By performing optimization in the HDR brachytherapy planning, the dose to urethra could also be kept within 120% of the target dose. From the dosimetry point of view, HDR brachytherapy would be a better choice for the treatment of early cancer of the prostate.
| References|| |
|1.||Brenner Hall EJ. Fractionation and protraction for radiotherapy of prostate carcinoma. Int J Radiat Oncol Biol Phys 1999;43:1095-101. |
|2.||Porter AT. Prostate brachytherapy: An overview. In: Martinez AA, Orton CG, Mould RF, editors. Brachytherapy HDR and LDR, proceedings brachytherapy meeting, Remote after loading: State of the art. Nucletron: 1989. p. 267-80. |
|3.||Poter AT, Scrimgour JW, Pocha JS. Remote interstitial after loading in cancer of the prostate: Preliminary experience with microSelectron. Int J Radiat Oncol Biol Phys 1988;14:571-5. |
|4.||Lachance B, Nadeau D, Lessard E. Early clinical experience with anatomy based inverse planning optimization for HDR boost of the prostate. Int J Radiat Oncol Biol Phys 2002;54:86-100. |
|5.||Alterovitz R, Lessard E, Pouliot J, Hsu IC, O'Brien JF, Goldberg K. Optimisation of high-dose rate brachytherapy dose distribution using linear programming with penalty constraints. Med Phys 2006;33:4012-9. |
|6.||Yoon M, Park SY, Shin D, Lee SB, Pyo HR, Kim DY, et al. A new homogeneity index based on statistical analysis of the dose-volume histogram. J Appl Clin Med Phys 2007;8:9-16. [PUBMED] |
|7.||International Commission on Radiation Units and Measurements (ICRU). Prescribing, recording, and reporting photon beam therapy. Supplement to ICRU Report 50, ICRU Report 62. Bethesda (MD): ICRU, 1999. |
|8.||Feuvret L, Noel G, Mazeron JJ, Bey P. Conformity index: A review. Int J Radiat Oncol Biol Phys 2006;64:333-42. |
|9.||Saw CB, Sunthralingam N. Quantitative assessment of interstitial implant. Int J Radiat Oncol Biol Phys 1990;20:135-9. |
|10.||Meertens H, Borger J, Steggerda M, Blom A. Evaluation and optimization of interstitial brachytherapy dose distribution. In: Mould RF, Battermann JJ, Martinez AA, Spesier BL. editors. Brachytherapy from radium to optimization. Nucletron; 1994. p. 300-6. |
|11.||Li S, Boyer A, Lu Y, Chen GT. Analysis of the dose-surface histogram and dose-wall histogram for rectum and bladder. Med Phys 1997;24:1107-16. [PUBMED] |
|12.||Fiorino C, Gianolini S, Nahum AE. A cylindrical model of the rectum: Comparing dose-volume, dose-surface and dose-wall histograms in the radiotherapy of prostate cancer. Phys Med Biol 2003;48:2603-16. [PUBMED] [FULLTEXT] |
|13.||Zhou SM, Marks LB, Tracton GS, Sibley GS, Light KL, Maguire PD, et al. A new three dimensional dose distribution reduction scheme for tubular organs. Med Phys 2000;27:1727-31. [PUBMED] |
|14.||Ting JY, Wu X, Fiedler JA, Yang C, Watzichand ML, Markoe A. Dose-volume histograms for bladder and rectum. Int J Radiat Oncol Biol Phys 1997;38:1105-11 . |
|15.||Lu Y, Li S, Spelbring D, Song P, Vijayakumar S, Pelizzari C, et al. Dose-surface histogram as treatment planning tool for prostate conformal therapy. Med Phys 1995;22:279-84. [PUBMED] |
|16.||John CR, Danny R, Spellbring D, Vijayakumar S, Jeffrey DF, Chen GT. Calculation of rectal dose surface histograms in the presence of time varying deformations. Int J Radiat Oncol Biol Phys 1996;34:391. |
|17.||Fowler JF. The linear quadratic formula and progress in fractionated radiotherapy: A review. Br J Radiol 1989;62:679-94. [PUBMED] |
|18.||Wang JZ, Li XA, Yu CX, Dibiase SJ. The low alpha/ beta ratio for prostate cancer: What does the clinical outcome of HDR brachytherapy tell us? Int J Radiat Oncol Biol Phys 2003;57:1101-8. [PUBMED] [FULLTEXT] |
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5]
[Table 1], [Table 2], [Table 3], [Table 4], [Table 5], [Table 6], [Table 7], [Table 8]
|This article has been cited by|
||Analysis of high–dose rate brachytherapy dose distribution resemblance in CyberKnife hypofractionated treatment plans of localized prostate cancer
| ||H. Sudahar,P.G.G. Kurup,V. Murali,P. Mahadev,J. Velmurugan |
| ||Medical Dosimetry. 2013; 38(4): 385 |
|[Pubmed] | [DOI]|
||Dosimetry of two new interstitial brachytherapy sources
| ||Saidi, P. and Sadeghi, M. |
| ||Polish Journal of Medical Physics and Engineering. 2012; 17(3): 129-141 |
||Equivalent normalized total dose estimates in cyberknife radiotherapy dose delivery in prostate cancer hypofractionation regimens
| ||Sudahar, H. and Kurup, P.G.G. and Murali, V. and Mahadev, P. and Velmurugan, J. |
| ||Journal of Medical Physics. 2012; 37(2): 90-96 |
||Dosimetry analysis of two new design 103Pd interstitial brachytherapy sources
| ||Saidi, P. and Sadeghi, M. and Shirazi, A. and Mahdavi, S.R. and Mirzaii, M. |
| ||International Journal of Low Radiation. 2011; 8(5-6): 363-373 |
||Dosimetry of two new interstitial brachytherapy sources
| ||Pooneh Saidi,Mahdi Sadeghi |
| ||Polish Journal of Medical Physics And Engineering. 2011; 17(3) |
|[Pubmed] | [DOI]|