- Open Access
Prostate external beam radiotherapy combined with high-dose-rate brachytherapy: dose-volume parameters from deformably-registered plans correlate with late gastrointestinal complications
Radiation Oncology volume 11, Article number: 144 (2016)
Derivation of dose-volume correlated with toxicity for multi-modal treatments can be difficult due to the perceived need for voxel-by-voxel dose accumulation. With data available for a single-institution cohort with long follow-up, an investigation was undertaken into rectal dose-volume effects for gastrointestinal toxicities after deformably-registering each phase of a combined external beam radiotherapy (EBRT)/high-dose-rate (HDR) brachytherapy prostate treatment.
One hundred and eighteen patients received EBRT in 23 fractions of 2 Gy and HDR (TG43 algorithm) in 3 fractions of 6.5 Gy. Results for the Late Effects of Normal Tissues — Subjective, Objective, Management and Analytic toxicity assessments were available with a median follow-up of 72 months. The HDR CT was deformably-registered to the EBRT CT. Doses were corrected for dose fractionation. Rectum dose-volume histogram (DVH) parameters were calculated in two ways. (1) Distribution-adding: parameters were calculated after the EBRT dose distribution was 3D-summed with the registered HDR dose distribution. (2) Parameter-adding: the EBRT DVH parameters were added to HDR DVH parameters. Logistic regressions and Mann-Whitney U-tests were used to correlate parameters with late peak toxicity (dichotomised at grade 1 or 2).
The 48–80, 40–63 and 49–55 Gy dose regions from distribution-adding were significantly correlated with rectal bleeding, urgency/tenesmus and stool frequency respectively. Additionally, urgency/tenesmus and anorectal pain were associated with the 25–26 Gy and 44–48 Gy dose regions from distribution-adding respectively. Parameter-adding also indicated the low-mid dose region was significantly correlated with stool frequency and proctitis.
This study confirms significant dose-histogram effects for gastrointestinal toxicities after including deformable registration to combine phases of EBRT/HDR prostate cancer treatment. The findings from distribution-adding were in most cases consistent with those from parameter-adding. The mid-high dose range and near maximum doses were important for rectal bleeding. The distribution-adding mid-high dose range was also important for stool frequency and urgency/tenesmus. We encourage additional studies in a variety of institutions using a variety of dose accumulation methods with appropriate inter-fraction motion management.
NCT NCT00193856. Retrospectively registered 12 September 2005.
External beam radiotherapy (EBRT) with a high-dose-rate brachytherapy (HDR) boost dose is used to treat prostate cancer patients . This treatment and other radiotherapy treatments are planned with consideration of the dose-volume parameters and subsequent constraints associated with acceptable levels of normal tissue toxicity . However, typically the phases of combined EBRT/HDR are planned separately . Constraints on the total planned dose from the two phases would be appropriate for reducing normal tissue toxicity . Constraints could be applied for each phase; however, this is susceptible to anatomical differences between the planning CTs.
When adjustments for anatomical changes are not included, the relevance of plans based on dose-volume constraints depends on how well the planned dose reflects the delivered dose . Hence, studies in other radiotherapy contexts have incorporated dose accumulation [6, 7]. Simple crude addition of the separately planned doses from two modalities is not valid as the anatomy in the CT image study sets may be misaligned due to variations in reference coordinate systems, displacements, deformations and shrinkage . Consequently, a rigid registration is used to align the reference coordinate systems and then deformable image registration (DIR) is applied to adjust for deformations and shrinkage [9, 10]. Additionally, the doses for different fraction schedules should be converted to the equieffective dose given in a reference X Gy per fraction (EQDXα/β Gy) as this adjusts for the biologically non-equivalent fractionation schedules [5, 11, 12].
Adjusting for anatomical differences between planning CTs and subsequently accumulating the phases of planned dose more accurately may allow dose-volume parameters to be more appropriately correlated with toxicity [2, 13]. Studies accumulating the rectum dose from phases of a combined EBRT/HDR prostate treatment by applying deformable registration are lacking. This study uses data from combined EBRT/HDR prostate cancer treatments, which were subject to multicentre trial guidelines, to assess whether the rectum dose-histogram parameters extracted after applying deformable registration are correlated with late gastrointestinal toxicities.
This study included 118 prostate cancer patients (tumour T stage ≥ 2a) who were treated with EBRT followed by HDR at Sir Charles Gairdner Hospital in the period 2004 to 2008. These patients were treated as part of the Trans-Tasman Radiation Oncology Group (TROG) 03.04 Randomized Androgen Deprivation and Radiotherapy (RADAR) trial [14, 15]. The patient criteria and treatment methodology for the RADAR trial have previously been detailed [14, 15]. Aspects of the combined EBRT/HDR treatment process have previously been described . The four-field EBRT plans for a prescription dose of 46 Gy in 23 daily fractions were created in the Elekta XiO treatment planning system (Elekta AB, Stockholm, Sweden). The HDR plans for a prescription dose of 19.5 Gy in 3 fractions across 2 days were created in the BrachyVision treatment planning system (Varian Medical Systems, Palo Alto, US) using the TG43 formalism . Additional file 1 (Supplement A) provides additional patient and treatment details.
The external wall of the rectum was manually delineated by treating-clinicians in the HDR CTs using BrachyVision and in the EBRT CTs using the Elekta Focal contouring system (Elekta AB, Stockholm, Sweden). Author with initials MK reviewed rectum outlines for consistency between patients. For rectum contouring, the inferior-superior limits of the rectum were the rectosigmoid flexure and the last slice where the ischial tuberosities were visible. Patients did not commonly require bowel preparation. Examples of the planning CTs and structures for EBRT and HDR TG43 physical dose plans are provided in Additional file 1 (Supplement B, Figures A1 and A2).
Patients were assessed for various gastrointestinal toxicities at baseline (randomisation) and subsequent time points after randomisation. The median of the most recent patient follow-ups was 72 months (range 12–96 months). The Late Effects of Normal Tissue — Subjective, Objective, Management, and Analytic (LENT SOMA) scales were used to assess rectal bleeding, urgency and tenesmus, stool frequency, diarrhoea, anorectal pain and completeness of evacuation . Proctitis was scored by clinicians according to the Common Toxicity Criteria for Adverse Events (CTCAE version 2) . Additional file 1 (Supplement C, Table A2) provides a summary of the grading systems.
Late peak toxicity was calculated for the period from 3 months after radiation therapy and onwards. Figure 1 provides a summary of the late peak toxicity event rates for the follow-up period. Patients were classified to a toxicity group if the late peak toxicity was at least a certain grade (threshold for dichotomisation). In the interest of modelling a moderate severity of toxicity the threshold was grade 2 for rectal bleeding, stool frequency and completeness of evacuation. The threshold was grade 1 for diarrhoea, anorectal pain, proctitis, urgency and tenesmus due to low toxicity rates for grade ≥ 2 toxicity and/or a lack of significance for grade ≥ 2 toxicity. Alternatively, the chosen thresholds for toxicities are indicated in Fig. 1. The analysis was repeated using the prevalence of toxicity at 36 months post-randomisation. This did not reveal any additional trends for dose-histogram effects and so is reported no further.
The HDR CT was registered to the EBRT CT using rigid registration followed by a B-splines multi-pass DIR in Velocity Advanced Imaging 2.8.1 (Varian Medical Systems, Palo Alto, US) . The registration process has been described in detail previously [20, 21]. Visual inspections for the 118 patients were undertaken by authors (initials CRM, VL and CIT) and did not identify any major registration misalignments (e.g. Additional file 1 [Supplement D, Figure A3]). The registrations were quantitatively evaluated for each patient in this study using the overlap of the EBRT rectum and registered HDR rectum (expressed as a percentage of the volume of the registered HDR rectum). As illustrated in Fig. 2, the median overlap is 80.4 % for alignment of EBRT/registered HDR rectum structure volumes. A general structure overlap of 70 % is considered to be the starting point for satisfactory structure-correspondence in the radiotherapy context [22, 23]. The registrations have also previously been extensively evaluated using structure-correspondence metrics, image similarity metrics and qualitative visual inspection by authors (initials CRM, VL and CIT) .
Obtaining dose-volume histograms
The EBRT and registered HDR 3D-doses were imported into MATLAB™ R2010a (The MathWorks Inc., Massachusetts, US) and the Computational Environment for Radiotherapy Research (CERR, version 4.1) . The voxel doses were converted to equieffective doses given in a reference 2 Gy per fraction using the linear-quadratic model  with an alpha-beta ratio of 3 Gy for the rectum . The analysis was also performed for an alpha-beta ratio of 5.4 Gy to check the sensitivity of results to the upper limit published for the rectum . The EBRT dose was summed voxel-by-voxel with the registered HDR dose (i.e. accumulated). The rectum dose-volume histograms (DVH) in 1 Gy bins from 1 to 80 Gy were extracted for the total registered dose (with EBRT rectum structure), the unregistered EBRT dose (with EBRT rectum structure) and the unregistered HDR dose (with the HDR rectum structure). The parameters extracted from the rectum DVHs were the VX (percentage of the rectal volume receiving at least X Gy after applying an alpha-beta ratio) and DX% (minimum dose to the most irradiated X percent of the rectal volume after applying an alpha-beta ratio). The VX and DX% were calculated using the total registered dose and the EBRT rectum structure (‘distribution-adding’). Additionally, the DX% was alternatively calculated by adding the EBRT DX% to the unregistered HDR DX% using the corresponding rectum structures (‘parameter-adding’).
For each type of toxicity, univariate logistic regression was applied at each VX to obtain an odds ratio (OR) for the increase in toxicity probability per 5 % absolute increase in volume . 95 % confidence intervals (CI) for odds ratios were calculated using bootstrapping with 10,000 resamples from the toxicity and no toxicity groups. Odds ratios were considered significant if the 95 % CIs did not include the value of one. Mann-Whitney U-tests were used to determine whether the median VX (or DX%) values for the toxicity and no toxicity groups were significantly different (p-value < 0.05). This analysis was performed in MATLAB™ R2010a (The MathWorks Inc., Massachusetts, US).
Clinical risk factors were not included in dose-response modelling as a previously published analysis determined that clinical covariates did not significantly influence late toxicities for patients in the RADAR trial . An equivalent analysis for the 118 patients in this study confirmed that clinical covariates did not significantly influence late toxicities. The clinical factors considered were age, tumour T stage, Gleason score, initial PSA, risk group, number of HDR catheters, colorectal disorders, hypertension, diabetes, smoking, use of statins, ACE inhibitors and anti-coagulants.
Unless it is stated otherwise, all figures and tables in this section report distribution-adding dose values in a reference 2 Gy per fraction using an α/β of 3 Gy.
Figure 3 provides the logistic regression odds ratio results for late rectal bleeding, stool frequency, diarrhoea, anorectal pain, urgency and tenesmus respectively. Figure 3 includes an indication for distribution-adding dose levels where the 95 % confidence intervals for the odds ratio did not include a value of one. For completeness of evacuation and proctitis, the odds ratios are not significantly different from one at any dose levels (Additional file 1 [Supplement E, Figure A4]).
Figure 4 provides the distribution-adding VX results for late rectal bleeding, stool frequency, diarrhoea, anorectal pain, urgency and tenesmus. Figure 5 provides the corresponding distribution-adding DX% results for late rectal bleeding, stool frequency, diarrhoea, urgency and tenesmus. Figures 4 and 5 include an indication of dose and volume levels for which there was a significant difference between the toxicity and no toxicity groups (p-value < 0.05). For completeness of evacuation and proctitis, there are no significant differences between the toxicity and no toxicity group results at any distribution-adding VX (Additional file 1 [Supplement E, Figure A5]). For completeness of evacuation, anorectal pain and proctitis, there are no significant differences between the toxicity and no toxicity group results at any distribution-adding DX% (Additional file 1 [Supplement E, Figure A6]).
Table 1 summarises the important distribution-adding doses for the odds ratios and the important distribution-adding doses (volumes) for the VX (DX%) results. Additionally, Table 1 summarises the important volumes for the DX% obtained by parameter-adding (alternatively, see Additional file 1 [Supplement F, Figure A7] for full results for the DX% obtained by parameter-adding). The D8-13% and D5-38% become important for diarrhoea and urgency/tenesmus respectively when distribution-adding is used instead of parameter-adding. The parameter-adding D59-70% and D58-73%, which were significant for proctitis and stool frequency respectively, are not significant when distribution-adding is used. However, the D56-74% becomes important for urgency/tenesmus when distribution-adding is used instead of parameter-adding. Similar trends for significance of DX% are found for all other toxicities regardless of whether distribution-adding or parameter-adding is used.
Table 2 summarises the influence of using an α/β of 5.4 Gy instead of 3 Gy. Alternatively, Additional file 1 (Supplement G, Figures A8, A9 and A10) provide the odds ratio, VX and DX% results when distribution-adding and an α/β of 5.4 Gy are used. Also, Additional file 1 (Supplement H, Figure A11) alternatively provides the DX% results obtained by parameter-adding with an α/β of 5.4 Gy. The odds ratio at EQD2α/β 51 Gy for diarrhoea is no longer significant if an α/β of 5.4 Gy is used instead of 3 Gy. However, similar trends for significance of odds ratios, VX and DX% are found for all other toxicities regardless of whether the α/β is 3 or 5.4 Gy.
It is important to explore dose-toxicity modelling in a variety of registration contexts
This study is the first to use registration-based distribution-addition to obtain accumulated rectum dose-histogram parameters for combined EBRT/HDR prostate cancer treatment and to then correlate the resulting parameters with gastrointestinal toxicities. Studies have estimated the accumulated rectum dose for combined EBRT/HDR prostate cancer treatment without applying deformable registration [3, 7]. However, Kikuchi et al.  acknowledged that deformable image registration should be part of a more accurate method of accumulating the rectum dose. This current study improved upon these studies by applying deformable image registration and then correlating accumulated rectum dose with various gastrointestinal toxicities. This study acknowledges the uncertainties of deformable image registration. Subsequently, it compares the findings for distribution-adding with the findings for parameter-adding. Given the potential uncertainties of deformable registration it is important for registration-based dose-toxicity modelling to be published for a variety of studies. This would allow a multi-institutional comparison of findings to include the confounding factors associated with different registration algorithms, registration circumstances, associated inter-fraction motion constraints and diversity in treatment techniques.
Studies had identified important dose-volume metrics for a variety of prostate radiotherapy techniques
The volume receiving certain doses and the magnitude of dose delivered to volumes have been associated with late gastrointestinal toxicities, typically rectal bleeding, scored after a number of prostate radiotherapy techniques including EBRT only, HDR only, low-dose-rate brachytherapy (LDR) only, combined EBRT/HDR and combined EBRT/LDR [2–4, 26, 28–45]. Table 3 summarises the important dose-response findings for the previously mentioned studies and the findings for this study.
The previously mentioned studies commonly suggested that the incidence of late rectal bleeding following prostate radiotherapy can be reduced by constraining the volume of the rectum receiving high doses (e.g. [4, 28, 42, 43, 46]). Additionally, some of the studies have correlated the mid and low-mid dose regions with late rectal bleeding [28, 30, 32, 33] and stool frequency/urgency/tenesmus [26, 30, 32, 37] respectively. Consequently, the mid-high rectum doses in prostate EBRT are typically managed through constraints on the V40-75 Gy [28, 46] whereas treatments involving prostate brachytherapy (HDR or LDR) should consider the high rectum doses via the V70%-100%, D1cc, D2cc and/or near maximum dose [4, 47–49] due to high-dose hot spots associated with radioactive sources. For treatments involving prostate HDR, the importance of low-dose regions has been explained in terms of considerable inter-patient variation in rectal gas and the distance from the prostate to the anterior rectal wall . The instances where the V80 Gy and V90 Gy have been identified as important for prostate cancer treatments involving EBRT only were related to homogenous irradiation of volumes with hypofractionated doses . The sections to follow will discuss the findings of this study, summarised in Table 3, relative to findings of the previously mentioned dose-response studies, which are also summarised in Table 3.
The findings indicate a serial response for rectal bleeding
In agreement with other studies [2, 28–30] the high-dose metrics for the rectum were significantly correlated with rectal bleeding for both distribution-adding and parameter-adding. The significant correlation between near maximum dose metrics for the rectum and rectal bleeding indicates the dose-volume effects follow a serial response. The confirmation of the expected importance of near maximum doses after registration is important as a previous study without registration did not find any significance for near maximum doses  and the GEC/ESTRO recommendation is to limit the D2cc to 75 Gy . Additionally, the identified seriality is consistent with the suggestion rectal bleeding is associated with epithelial damage and mucositis as a result of exposure of parts of the rectal wall to near maximum doses .
The mid-dose region is important for bleeding/non-bleeding toxicities
Studies have also demonstrated that the mid-dose region (>30 Gy) is important for rectal bleeding [28, 30, 32, 33]. In this study the importance is shifted to relatively higher doses in the mid-high dose range for both distribution-adding and parameter-adding. An influencing factor for the lack of importance of the lower end of the mid-dose range could be that the combined EBRT/HDR treatments were subject to the constraint that the maximum rectum dose from HDR should not exceed 80 % of the 19.5 Gy prescription dose for HDR. Consequently, in the context of the total EBRT/HDR dose this constraint effectively applies more to the lower end of the mid-dose range after adjusting for dose fractionation than it does to the high-dose region. The importance of the high-dose and near maximum dose regions could also be related to the steepness of dose gradients associated with HDR treatments as it has been proposed that a focused high-dose region could aid healing of the vascular sclerosis in high-dose regions via cell migration from the low-dose region . Consequently, it would be useful to determine optimal rectum dose constraints for combined EBRT/HDR based on accumulated dose. A larger sample size containing patients from a variety of institutions would allow for a feasible application of multivariate and cut-point analysis.
The upper end of the mid-high dose range after distribution-adding was important for the non-bleeding toxicities of stool frequency and urgency/tenesmus. This result could support the earlier suggestion that the rectum dose constraint for HDR effectively applies a constraint to the lower end of the mid-dose range when the total EBRT/HDR dose is considered. The dose constraint in one dose region leading to other dose regions becoming important is consistent with a previous study focusing on patients within this trial who received EBRT only . The study indicated the low-mid dose range was important for stool frequency, urgency and tenesmus in the presence of high-dose constraints . More optimised dose constraints for the mid-high dose range based on accumulated dose could be useful for reducing the toxicities associated with these doses. Such constraints could be relatively more important for urgency/tenesmus compared to stool frequency given the higher toxicity rate in this patient sample compared to the other toxicities.
Toxicity is also influenced by low doses and the lower end of the mid-dose range
The association of urgency/tenesmus with distribution-adding doses at the lower end of the mid-dose range is consistent with the finding from another study where violation of the V40 Gy dose constraint was important for urgency . Additionally, the results indicate the lower end of the mid-dose range and the low doses may be associated with anorectal pain and rectal bleeding respectively. The correlation of toxicities with low doses and the lower end of the mid-dose range is possible as it is plausible that a low-dose bath to a large volume will be associated with detriment. However, these findings of association should be considered with respect to toxicity event rates, sample size and the potential of random discovery.
Software developments to improve contour consistency and registration accuracy for the prostate/rectum interface would be of great benefit
Analysis based on contouring and registration is associated with uncertainties. However, the dose regions indicated as being important for toxicity after distribution-adding were in most cases consistent with those indicated as important after parameter-adding. The low-mid dose range for parameter-adding was significantly associated with proctitis and stool frequency. In contrast, these regions after distribution-adding were not identified as important. However, distribution-adding did indicate additional regions as important compared to regions identified by the parameter-adding results. For example, the analysis for parameter-adding did not indicate any significant dose regions for diarrhoea and urgency/tenesmus whereas analysis after distribution-adding indicated the mid-high dose range was important. The alpha-beta ratio is an additional uncertainty for diarrhoea correlations as the mid-high dose range was only important for an alpha-beta ratio of 3 Gy. Further studies for a variety of contouring and registration contexts would be useful for gathering data for the purpose of determining whether registration and distribution-adding reveals correlations which were not identified by parameter-adding.
When considering the distribution-adding findings in isolation it should be noted that errors in contouring and registration accuracy will confound the distribution-adding parameters that have been correlated with toxicity. A median overlap of 80.4 % for the rectum volume correspondence across all patients would indicate the registrations are satisfactory as a general structure overlap of 70 % is considered to be the starting point for satisfactory structure-correspondence in the radiotherapy context [22, 23]. The proximity of the HDR rectum to the HDR catheters makes parameters obtained by registration and distribution-adding sensitive to small localised variations in contouring and/or registration accuracy across the prostate/rectum interface. Given this is the first study to accumulate the rectum dose for combined EBRT/HDR prostate cancer treatments using deformable registration and then correlate the doses with toxicity, we encourage more prostate cancer studies to assess the importance of dose-volume metrics using a variety of registration algorithms. Software developers and treatment planning vendors have the opportunity to greatly improve planning and the reliability of dose-toxicity modelling after registration by improving contouring and registration accuracy for the prostate/rectum interface.
Inter-fraction motion should be considered
A common uncertainty associated with EBRT, HDR and other radiotherapy techniques is inter-fraction motion of patient anatomy [51–54]. In response to this uncertainty, it is becoming more common for institutions to adopt repeat imaging over the course of prostate cancer treatment to correct for inter-fraction motions and improve the correspondence between planned dose and delivered dose [53, 55]. However, many studies including this study do not contain repeat imaging due to the retrospective nature of studies where long follow-up is required to correlate dose with late toxicities. Consequently, studies are constrained by treatments performed in the past with the associated resources and protocols at that time.
A consideration for prostate EBRT in this study is rectum motion and variable rectum contents confounding the accuracy of rectum dose distributions obtained from single static planning CTs [52, 56, 57]. A variety of methods have been used to estimate the impact of inter-fraction motion on rectum dose parameters and dose-response modelling for prostate EBRT [8, 52, 58–61]. To obtain appropriate mean estimates for the difference in EBRT rectum dose between the single CT based values and the motion-corrected values this study analysed the results of another study  that used the same registration software and registration algorithm. Consequently, compared to the motion-corrected values the single CT based values may be conservative estimates by 3.9 % for the D2% and 5.8 % for the equivalent uniform dose .
An important consideration for inter-fraction motion during prostate brachytherapy is the movement of the anterior rectal wall relative to the prostate [54, 62]. A variety of methods have been used to estimate inter-fraction motion in prostate HDR [16, 54] and the subsequent impact on rectum dose parameters for prostate brachytherapy [54, 63–65]. Simnor et at.  calculated that the catheter mean caudal displacements of 7.9 mm and 3.8 mm prior to fractions 2 and 3 were associated with mostly systematic increases to the D2cc of 0.69 Gy (~ 6.6 %) and 0.76 Gy (~ 7.2 %) respectively. For HDR at the institution where patients in this study were treated, the displacement of catheters was checked prior to each of the three fractions using an anterior-posterior radiograph and corrected for using a rigid external holding device as described by Tiong et al. . The catheter mean caudal displacements after this advancement process were reported as 1.7 mm, 1.1 mm and 0.8 mm for fractions 1, 2 and 3 respectively . Consequently, the inter-fraction motion increases to the D2cc reported by Simnor et al.  may be appropriate conservative estimates for the inter-fraction motion of HDR catheters that could be expected for this study.
It is possible that the above mentioned inter-fraction motion could remove significance of dose ranges. However, shifting of dose values identified as significant are likely as the single CT based estimates were mostly identified as being systematically different to the motion-corrected values [52, 54]. The influence of inter-fraction on delivered doses is likely to be important when considering dose constraints recommended by studies where planned dose tends to be less than delivered dose. Consequently, it would be useful to confirm the importance of published dose-volume constraints after registration is applied for repeat daily imaging.
Avenues and recommendations for further analysis
Given this is the first study to apply deformable registration prior to correlating combined EBRT/HDR dose with toxicity, it is important that the model and findings be validated in other contexts with standardised contouring, implanting and planning guidelines for EBRT and HDR.
A larger sample size would make it feasible to explore models that incorporate multiple toxicity events over the follow-up period  or include the persistence of toxicity rather than peak late toxicity .
Image guided radiotherapy or further imaging could improve the reliability of accumulated dose-histogram metrics .
Customised registration algorithms for accurately handling the catheters within the HDR prostate or data for treatments which use plastic HDR catheters are encouraged as prostate and urethra doses are key clinical concerns in the RADAR trial .
It would also be useful to determine whether other aspects of the total registered dose distribution add predictive capability to dose-toxicity modelling e.g. including dose-shape toxicity modelling .
Exploring the association between toxicity and doses to other organs or regions may be useful for further explaining the incidence of toxicity (e.g. doses to the bowel and gastrointestinal tract could be associated with toxicity [46, 68]).
A number of significant dose-histogram effects were revealed for gastrointestinal toxicities after applying deformable registration to adjust for the anatomical differences between planning CTs for each phase of a combine EBRT/HDR prostate cancer treatment. The findings for distribution-adding were in most cases consistent with those for parameter-adding. The mid-high dose range and near maximum doses were important for rectal bleeding. The distribution-adding mid-high dose range was also important for stool frequency and urgency/tenesmus. The anorectum doses which were important for toxicity are reported to guide and encourage future planning of combined EBRT/HDR prostate cancer treatments based on accumulated phases with appropriate inter-fraction motion management. We encourage other studies to report on important dose-histogram effects and spatial aspects of accumulated dose distributions for combined EBRT/HDR.
Computational environment for radiotherapy research
Common toxicity criteria for adverse events
Deformable image registration
- DX% :
Minimum dose to the most irradiated X percent of the rectal volume after applying an alpha-beta ratio
External beam radiotherapy
- EQDXα/β Gy:
Equieffective dose given in a reference X Gy per fraction using α/β
- LENT SOMA:
Late effects of normal tissue — subjective, objective, management and analytic
Randomized androgen deprivation and radiotherapy trial
Trans-Tasman Radiation Oncology Group
- VX :
Percentage of the rectal volume receiving at least X Gy after applying an alpha-beta ratio
Zaorsky NG, Doyle LA, Yamoah K, Andrel JA, Trabulsi EJ, Hurwitz MD, et al. High dose rate brachytherapy boost for prostate cancer: a systematic review. Cancer Treat Rev. 2014; doi:10.1016/j.ctrv.2013.10.006.
Michalski JM, Gay H, Jackson A, Tucker SL, Deasy JO. Radiation dose-volume effects in radiation-induced rectal injury. Int J Radiat Oncol Biol Phys. 2010; doi:10.1016/j.ijrobp.2009.03.078.
Okamoto M, Ishikawa H, Ebara T, Kato H, Tamaki T, Akimoto T, et al. Rectal bleeding after high-dose-rate brachytherapy combined with hypofractionated external-beam radiotherapy for localized prostate cancer: the relationship between dose-volume histogram parameters and the occurrence rate. Int J Radiat Oncol Biol Phys. 2012; doi:10.1016/j.ijrobp.2011.03.041.
Hoskin PJ, Colombo A, Henry A, Niehoff P, Paulsen Hellebust T, Siebert FA, et al. GEC/ESTRO recommendations on high dose rate afterloading brachytherapy for localised prostate cancer: an update. Radiother Oncol. 2013; doi:10.1016/j.radonc.2013.05.002.
Vásquez Osorio EM, Hoogeman MS, Teguh DN, Al-Mamgani A, Kolkman-Deurloo IK, Bondar L, et al. Three-dimensional dose addition of external beam radiotherapy and brachytherapy for oropharyngeal patients using nonrigid registration. Int J Radiat Oncol Biol Phys. 2011; doi:10.1016/j.ijrobp.2010.10.006.
Wortel RC, Incrocci L, Pos FJ, van der Heide UA, Lebesque JV, Aluwini S, et al. Late side effects after image guided intensity modulated radiation therapy compared to 3D-conformal radiation therapy for prostate cancer: results from 2 prospective cohorts. Int J Radiat Oncol Biol Phys. 2016; doi:10.1016/j.ijrobp.2016.01.031.
Tyagi N, Sebastian E, Liang J, Yan D, Ghilezan M, Martinez A. Deformable registration and dose accumulation for image-guided HDR interstitial brachytherapy (IG-BT) boost and external beam pelvic IMRT (EB-IMRT) for intermediate and high-risk prostate cancer patients. 52nd annual meeting of the American Society for Radiation Oncology. Int J Radiat Oncol Biol Phys. 2010; doi:10.1016/j.ijrobp.2010.07.1695.
Thor M, Apte A, Deasy JO, Karlsdóttir À, Moiseenko V, Liu M, et al. Dose/volume-response relations for rectal morbidity using planned and simulated motion-inclusive dose distributions. Radiother Oncol. 2013; doi:10.1016/j.radonc.2013.10.021.
Teo BK, Bonner Millar LP, Ding X, Lin LL. Assessment of cumulative external beam and intracavitary brachytherapy organ doses in gynecologic cancers using deformable dose summation. Radiother Oncol. 2015; doi:10.1016/j.radonc.2015.04.002.
Vestergaard A, Hafeez S, Muren LP, Nill S, Høyer M, Hansen VN, et al. The potential of MRI-guided online adaptive re-optimisation in radiotherapy of urinary bladder cancer. Radiother Oncol. 2016; doi:10.1016/j.radonc.2015.11.003.
Bentzen SM, Dörr W, Gahbauer R, Howell RW, Joiner MC, Jones B, et al. Bioeffect modeling and equieffective dose concepts in radiation oncology – terminology, quantities and units. Radiother Oncol. 2012; doi:10.1016/j.radonc.2012.10.006.
Kikuchi K, Nakamura R, Tanji S, Yamaguchi S, Kakuhara H, Yabuuchi T, et al. Three-dimensional summation of rectal doses in brachytherapy combined with external beam radiotherapy for prostate cancer. Radiother Oncol. 2013; doi:10.1016/j.radonc.2013.03.003.
Troeller A, Yan D, Marina O, Schulze D, Alber M, Parodi K, et al. Comparison and limitations of DVH-based NTCP models derived from 3D-CRT and IMRT data for prediction of gastrointestinal toxicities in prostate cancer patients by using propensity score matched pair analysis. Int J Radiat Oncol Biol Phys. 2015; doi:10.1016/j.ijrobp.2014.09.046.
Denham JW, Joseph D, Lamb DS, Spry NA, Duchesne G, Matthews J, et al. Short-term androgen suppression and radiotherapy versus intermediate-term androgen suppression and radiotherapy, with or without zoledronic acid, in men with locally advanced prostate cancer (TROG 03.04 RADAR): an open-label, randomised, phase 3 factorial trial. Lancet Oncol. 2014; doi:10.1016/S1470-2045(14)70328-6.
Denham JW, Steigler A, Joseph D, Lamb DS, Spry NA, Duchesne G, et al. Radiation dose escalation or longer androgen suppression for locally advanced prostate cancer? Data from the TROG 03.04 RADAR trial. Radiother Oncol. 2015; doi:10.1016/j.radonc.2015.05.016.
Tiong A, Bydder S, Ebert M, Caswell N, Waterhouse D, Spry N, et al. A small tolerance for catheter displacement in high-dose rate prostate brachytherapy is necessary and feasible. Int J Radiat Oncol Biol Phys. 2010; doi:10.1016/j.ijrobp.2009.03.052.
Rivard MJ, Coursey BM, DeWerd LA, Hanson WF, Huq MS, Ibbott GS, et al. Update of AAPM task group no. 43 report: a revised AAPM protocol for brachytherapy dose calculations. Med Phys. 2004; doi:10.1118/1.1646040.
RTOG/EORTC Late Effects Working Group. Late effects on normal tissue consensus conference: lent soma scales for all anatomic sites. Int J Radiat Oncol Biol Phys. 1995; doi:10.1016/0360-3016(95)90159-0.
Trotti A, Byhardt R, Stetz J, Gwede C, Corn B, Fu K, et al. Common toxicity criteria: version 2.0. An improved reference for grading the acute effects of cancer treatment: impact on radiotherapy. Int J Radiat Oncol Biol Phys. 2000; doi:10.1016/S0360-3016(99)00559-3.
Kadoya N, Fujita Y, Katsuta Y, Dobashi S, Takeda K, Kishi K, et al. Evaluation of various deformable image registration algorithms for thoracic images. J Radiat Res. 2014; doi:10.1093/jrr/rrt093.
Moulton CR, House MJ, Lye V, Tang CI, Krawiec M, Joseph DJ, et al. Registering prostate external beam radiotherapy with a boost from high-dose-rate brachytherapy: a comparative evaluation of deformable registration algorithms. Radiat Oncol. 2015; doi:10.1186/s13014-015-0563-9.
Zou KH, Warfield SK, Bharatha A, Tempany CM, Kaus MR, Haker SJ, et al. Statistical validation of image segmentation quality based on a spatial overlap index. Acad Radiol. 2004; doi:10.1016/S1076-6332(03)00671-8.
Acosta O, Dréan G, Ospina JD, Simon A, Haigron P, Lafond C, et al. Voxel-based population analysis for correlating local dose and rectal toxicity in prostate cancer radiotherapy. Phys Med Biol. 2013; doi:10.1088/0031-9155/58/8/2581.
Deasy JO, Blanco AI, Clark VH. CERR: a computational environment for radiotherapy research. Med Phys. 2003; doi:10.1118/1.1568978.
Brenner DJ. Fractionation and late rectal toxicity. Int J Radiat Oncol Biol Phys. 2004; doi:10.1016/j.ijrobp.2004.04.014.
Gulliford SL, Foo K, Morgan RC, Aird EG, Bidmead AM, Critchley H, et al. Dose-volume constraints to reduce rectal side effects from prostate radiotherapy: evidence from MRC RT01 trial ISRCTN 47772397. Int J Radiat Oncol Biol Phys. 2010; doi:10.1016/j.ijrobp.2009.02.025.
Denham JW, Wilcox C, Lamb DS, Spry NA, Duchesne G, Atkinson C, et al. Rectal and urinary dysfunction in the TROG 03.04 RADAR trial for locally advanced prostate cancer. Radiother Oncol. 2012; doi:10.1016/j.radonc.2012.09.018.
Fiorino C, Valdagni R, Rancati T, Sanguineti G. Dose-volume effects for normal tissues in external radiotherapy: pelvis. Radiother Oncol. 2009; doi:10.1016/j.radonc.2009.08.004.
Boersma LJ, van den Brink M, Bruce AM, Shouman T, Gras L, te Velde A, et al. Estimation of the incidence of late bladder and rectum complications after high-dose (70-78 Gy) conformal radiotherapy for prostate cancer, using dose-volume histograms. Int J Radiat Oncol Biol Phys. 1998; doi:10.1016/S0360-3016(98)00037-6.
Ebert MA, Foo K, Haworth A, Gulliford SL, Kennedy A, Joseph DJ, et al. Gastrointestinal dose-histogram effects in the context of dose-volume-constrained prostate radiation therapy: analysis of data from the RADAR prostate radiation therapy trial. Int J Radiat Oncol Biol Phys. 2015; doi:10.1016/j.ijrobp.2014.11.015.
Lalla RV, Bowen J, Barasch A, Elting L, Epstein J, Keefe DM, et al. MASCC/ISOO clinical practice guidelines for the management of mucositis secondary to cancer therapy. Cancer. 2014; doi:10.1002/cncr.28592.
Peeters ST, Lebesque JV, Heemsbergen WD, van Putten WL, Slot A, Dielwart MF, et al. Localized volume effects for late rectal and anal toxicity after radiotherapy for prostate cancer. Int J Radiat Oncol Biol Phys. 2006; doi:10.1016/j.ijrobp.2005.10.002.
Vargas C, Martinez A, Kestin LL, Yan D, Grills I, Brabbins DS, et al. Dose-volume analysis of predictors for chronic rectal toxicity after treatment of prostate cancer with adaptive image-guided radiotherapy. Int J Radiat Oncol Biol Phys. 2005; doi:10.1016/j.ijrobp.2004.12.052.
Konishi K, Yoshioka Y, Isohashi F, Sumida I, Kawaguchi Y, Kotsuma T, et al. Correlation between dosimetric parameters and late rectal and urinary toxicities in patients treated with high-dose-rate brachytherapy used as monotherapy for prostate cancer. Int J Radiat Oncol Biol Phys. 2009; doi:10.1016/j.ijrobp.2008.12.051.
Kim DW, Cho LC, Straka C, Christie A, Lotan Y, Pistenmaa D, et al. Predictors of rectal tolerance observed in a dose-escalated phase 1-2 trial of stereotactic body radiation therapy for prostate cancer. Int J Radiat Oncol Biol Phys. 2014; doi:10.1016/j.ijrobp.2014.03.012.
Pederson AW, Fricano J, Correa D, Pelizzari CA, Liauw SL. Late toxicity after intensity-modulated radiation therapy for localized prostate cancer: an exploration of dose-volume histogram parameters to limit genitourinary and gastrointestinal toxicity. Int J Radiat Oncol Biol Phys. 2012; doi:10.1016/j.ijrobp.2010.09.058.
Stenmark MH, Conlon AS, Johnson S, Daignault S, Litzenberg D, Marsh R, et al. Dose to the inferior rectum is strongly associated with patient reported bowel quality of life after radiation therapy for prostate cancer. Radiother Oncol. 2014; doi:10.1016/j.radonc.2014.01.007.
Yoshioka Y, Nose T, Yoshida K, Oh RJ, Yamada Y, Tanaka E, et al. High-dose-rate brachytherapy as monotherapy for localized prostate cancer: a retrospective analysis with special focus on tolerance and chronic toxicity. Int J Radiat Oncol Biol Phys. 2003; doi:10.1016/S0360-3016(03)00081-6.
Fellin G, Rancati T, Fiorino C, Vavassori V, Antognoni P, Baccolini M, et al. Long term rectal function after high-dose prostatecancer radiotherapy: results from a prospective cohort study. Radiother Oncol. 2014; doi:10.1016/j.radonc.2013.09.028.
Fonteyne V, Ost P, Vanpachtenbeke F, Colman R, Sadeghi S, Villeirs G, et al. Rectal toxicity after intensity modulated radiotherapy for prostate cancer: which rectal dose volume constraints should we use? Radiother Oncol. 2014; doi:10.1016/j.radonc.2014.10.014.
Someya M, Hori M, Tateoka K, Nakata K, Takagi M, Saito M, et al. Results and DVH analysis of late rectal bleeding in patients treated with 3D-CRT or IMRT for localized prostate cancer. J Radiat Res. 2015; doi:10.1093/jrr/rru080.
Price JG, Stone NN, Stock RG. Predictive factors and management of rectal bleeding side effects following prostate cancer brachytherapy. Int J Radiat Oncol Biol Phys. 2013; doi:10.1016/j.ijrobp.2013.04.033.
Sutani S, Ohashi T, Sakayori M, Kaneda T, Yamashita S, Momma T, et al. Comparison of genitourinary and gastrointestinal toxicity among four radiotherapy modalities for prostate cancer: conventional radiotherapy, intensity-modulated radiotherapy, and permanent iodine-125 implantation with or without external beam radiotherapy. Radiother Oncol. 2015; doi:10.1016/j.radonc.2015.08.019.
Akimoto T, Muramatsu H, Takahashi M, Saito J, Kitamoto Y, Harashima K, et al. Rectal bleeding after hypofractionated radiotherapy for prostate cancer: correlation between clinical and dosimetric parameters and the incidence of grade 2 or worse rectal bleeding. Int J Radiat Oncol Biol Phys. 2004; doi:10.1016/j.ijrobp.2004.07.695.
Akimoto T, Katoh H, Kitamoto Y, Tamaki T, Harada K, Shirai K, et al. Rectal bleeding after high-dose-rate brachytherapy combined with hypofractionated external-beam radiotherapy for localized prostate cancer: impact of rectal dose in high-dose-rate brachytherapy on occurrence of grade 2 or worse rectal bleeding. Int J Radiat Oncol Biol Phys. 2006; doi:10.1016/j.ijrobp.2005.12.017.
Landoni V, Fiorino C, Cozzarini C, Sanguineti G, Valdagni R, Rancati T. Predicting toxicity in radiotherapy for prostate cancer. Phys Med. 2016; doi:10.1016/j.ejmp.2016.03.003.
van Nieuwenhuysen J, Waterhouse D, Bydder S, Joseph D, Ebert M, Caswell N. Survey of high-dose-rate prostate brachytherapy practice in Australia and New Zealand, 2010-2011. J Med Imaging Radiat Oncol. 2014; doi:10.1111/1754-9485.12101.
Kirisits C, Goldner G, Berger D, Georg D, Pötter R. Critical discussion of different dose-volume parameters for rectum and urethra in prostate brachytherapy. Brachytherapy. 2009; doi:10.1016/j.brachy.2009.01.003.
Peters M, Hoekstra CJ, van der Voort van Zyp JR, Westendorp H, van de Pol SM, Moerland MA, et al. Rectal dose constraints for salvage iodine-125 prostate brachytherapy. Brachytherapy. 2016; doi:10.1016/j.brachy.2015.10.004.
Munbodh R, Jackson A. Quantifying cell migration distance as a contributing factor to the development of rectal toxicity after prostate radiotherapy. Med Phys. 2014; doi:10.1118/1.4852955.
Kirisits C, Rivard MJ, Baltas D, Ballester F, De Brabandere M, van der Laarse R, et al. Review of clinical brachytherapy uncertainties: analysis guidelines of GEC-ESTRO and the AAPM. Radiother Oncol. 2014; doi:10.1016/j.radonc.2013.11.002.
Akino Y, Yoshioka Y, Fukuda S, Maruoka S, Takahashi Y, Yagi M, et al. Estimation of rectal dose using daily megavoltage cone-beam computed tomography and deformable image registration. Int J Radiat Oncol Biol Phys. 2013; doi:10.1016/j.ijrobp.2013.06.2054.
Thörnqvist S, Hysing LB, Tuomikoski L, Vestergaard A, Tanderup K, Muren LP, et al. Adaptive radiotherapy strategies for pelvic tumors - a systematic review of clinical implementations. Acta Oncol. 2016; doi:10.3109/0284186X.2016.1156738.
Simnor T, Li S, Lowe G, Ostler P, Bryant L, Chapman C, et al. Justification for inter-fraction correction of catheter movement in fractionated high dose-rate brachytherapy treatment of prostate cancer. Radiother Oncol. 2009; doi:10.1016/j.radonc.2009.09.015.
McPartlin AJ, Li XA, Kershaw LE, Heide U, Kerkmeijer L, Lawton C, et al. MRI-guided prostate adaptive radiotherapy - a systematic review. Radiother Oncol. 2016; doi:10.1016/j.radonc.2016.04.014.
Stasi M, Munoz F, Fiorino C, Pasquino M, Baiotto B, Marini P, et al. Emptying the rectum before treatment delivery limits the variations of rectal dose-volume parameters during 3DCRT of prostate cancer. Radiother Oncol. 2006; doi:10.1016/j.radonc.2006.08.007.
Ariyaratne H, Chesham H, Pettingell J, Alonzi R. Image-guided radiotherapy for prostate cancer with cone beam CT: dosimetric effects of imaging frequency and PTV margin. Radiother Oncol. 2016; doi:10.1016/j.radonc.2016.07.018.
Sripadam R, Stratford J, Henry AM, Jackson A, Moore CJ, Price P. Rectal motion can reduce CTV coverage and increase rectal dose during prostate radiotherapy: a daily cone-beam CT study. Radiother Oncol. 2009; doi:10.1016/j.radonc.2008.07.031.
Hatton JA, Greer PB, Tang C, Wright P, Capp A, Gupta S, et al. Does the planning dose-volume histogram represent treatment doses in image-guided prostate radiation therapy? Assessment with cone-beam computerised tomography scans. Radiother Oncol. 2011; doi:10.1016/j.radonc.2011.01.006.
Scaife JE, Thomas SJ, Harrison K, Romanchikova M, Sutcliffe MP, Forman JR, et al. Accumulated dose to the rectum, measured using dose-volume histograms and dose-surface maps, is different from planned dose in all patients treated with radiotherapy for prostate cancer. Br J Radiol. 2015; doi:10.1259/bjr.20150243.
Nassef M, Simon A, Cazoulat G, Duménil A, Blay C, Lafond C, et al. Quantification of dose uncertainties in cumulated dose estimation compared to planned dose in prostate IMRT. Radiother Oncol. 2016; doi:10.1016/j.radonc.2016.03.007.
Merrick GS, Butler WM, Dorsey AT, Dorsey JT. The effect of constipation on rectal dosimetry following prostate brachytherapy. Med Dosim. 2000; doi:10.1016/S0958-3947(00)00047-9.
Foster W, Cunha JA, Hsu IC, Weinberg V, Krishnamurthy D, Pouliot J. Dosimetric impact of interfraction catheter movement in high-dose rate prostate brachytherapy. Int J Radiat Oncol Biol Phys. 2011; doi:10.1016/j.ijrobp.2010.01.016.
Kovalchuk N, Furutani KM, MacDonald OK, Pisansky TM. Dosimetric effect of interfractional needle displacement in prostate high-dose-rate brachytherapy. Brachytherapy. 2012; doi:10.1016/j.brachy.2011.05.006.
Dinkla AM, Pieters BR, Koedooder K, Meijnen P, van Wieringen N, van der Laarse R, et al. Deviations from the planned dose during 48 hours of stepping source prostate brachytherapy caused by anatomical variations. Radiother Oncol. 2013; doi:10.1016/j.radonc.2012.12.011.
Yahya N, Ebert MA, Bulsara M, House MJ, Kennedy A, Joseph DJ, et al. Urinary symptoms following external beam radiotherapy of the prostate: dose-symptom correlates with multiple-event and event-count models. Radiother Oncol. 2015; doi:10.1016/j.radonc.2015.10.003.
Gulliford SL, Partridge M, Sydes MR, Andreyev J, Dearnaley DP. A comparison of dose-volume constraints derived using peak and longitudinal definitions of late rectal toxicity. Radiother Oncol. 2010; doi:10.1016/j.radonc.2009.12.019.
Wortel RC, Witte MG, van der Heide UA, Pos FJ, Lebesque JV, van Herk M, et al. Dose-surface maps identifying local dose-effects for acute gastrointestinal toxicity after radiotherapy for prostate cancer. Radiother Oncol. 2015; doi:10.1016/j.radonc.2015.10.020.
We thank Annette Haworth and radiation oncology staff at Sir Charles Gairdner Hospital for their contributions. We also appreciate the support of the RADAR centres and the Trans-Tasman Radiation Oncology Group.
We acknowledge funding from the National Health and Medical Research Council (300705, 455521, 1006447), the University of Western Australia, an Australian Postgraduate Award, an Ana Africh Scholarship, the Hunter Medical Research Institute, the Health Research Council (New Zealand), Abbott Laboratories and Novartis Pharmaceuticals.
Availability of data and material
CRM, MJH, CIT and MAE have made substantial contributions to design of analysis. CRM, VL, CIT, MK, DJJ, JWD and MAE have made substantial contributions to acquire the data. CRM, MJH, VL, CIT and MAE substantially contributed to analysis and interpretation of data. CRM has been involved in drafting the manuscript. MJH, VL, CIT, DJJ, JWD and MAE revised it critically for important and correct content. JWD and DJJ were involved in the design and coordination of the RADAR trial. All authors read and approved the version to be published.
The authors declare that they have no competing interests.
Consent for publication
The signed patient consent forms for the trial informed patients that their medical information may be used to publish the results of the study. In accordance with the signed patient consent forms, this publication includes only anonymized information and does not include information identifying any patient.
Ethics approval and consent to participate
The TROG 03.04 RADAR Trial is registered with the National Institutes of Health Clinical Trials Registry (number NCT00193856). This trial has approval from the Hunter New England Human Research Ethics Committee (Trial ID. 03/06/11/3.02), the Sir Charles Gairdner Group Human Research Ethics Committee (2003-050) and the University of Western Australia Human Research Ethics Office (RA/4/1/5601). Patients participating in the trial signed consent forms.
Online supplementary material (Supplements A-H.pdf) providing additional treatment details, patient characteristics, visual checks of registrations, results which were not significant, results for α/β = 5.4 Gy and results for parameter-adding. (PDF 2267 kb)
About this article
Cite this article
Moulton, C.R., House, M.J., Lye, V. et al. Prostate external beam radiotherapy combined with high-dose-rate brachytherapy: dose-volume parameters from deformably-registered plans correlate with late gastrointestinal complications. Radiat Oncol 11, 144 (2016). https://doi.org/10.1186/s13014-016-0719-2
- Deformable registration
- Gastrointestinal toxicity