Variations in Rectal Volumes and Dosimetry Values Including NTCP due to Interfractional Variability When Administering 2D-Based IG-IMRT for Prostate Cancer
We estimated variations in rectal volumes and dosimetry values including NTCP with interfractional motion during prostate IG-IMRT. Rectal volumes, DVH parameters, and NTCPs of 20 patients were analyzed. For this patient population, the median (range) volume on the initial plan for the rectum was 45.6 cc (31.3–82.0), showing on-treatment spread around the initial prediction based on the initial plan. DVH parameters of on-treatment CBCT analyses showed systematic regularity shift from the prediction based on the initial plan. Using the Lyman-Kutcher-Burman model, NTCPs of predicted late rectal bleeding toxicity of rectal grade ≥ 2 (RTOG) and the QUANTEC update rectal toxicity for the prediction based on the initial plan were 0.09% (0.02–0.24) and 0.02% (0.00–0.07), respectively, with NTCPs from on-treatment CBCT analyses being 0.35% (0.01–6.16) and 0.12% (0.00–4.11), respectively. Using the relative seriality model, for grade ≥ 2 bleeding rectal toxicity, NTCP of the prediction based on the initial plan was 0.64% (0.15–1.22) versus 1.48% (0.18–7.66) for on-treatment CBCT analysis. Interfraction variations in rectal volumes occur in all patients due to physiological changes. Thus, rectal assessment during 2D-based IG-IMRT using NTCP models has the potential to provide useful and practical dosimetric verification.
Intensity modulated radiotherapy (IMRT) allows successful delivery of escalated dose radiation to the tumor target with high precision, while simultaneously sparing sensitive adjacent tissues. This irradiation technique is adequate for treatment of prostate cancer, as increasing the dose to the prostate achieves better local control, but there is a risk of increased late toxicity . The success of IMRT is related to the ability to accurately target the prostate during treatment. Several factors influence this ability, such as the accuracy of organ delineation, adequate patient immobilization, and intrafraction and interfraction organ motions. Because the treatment plan is normally based on treatment planning computed tomography (CT) images, these errors could result in the dose distribution delivered being different from that planned . This may result in insufficient dose coverage to the prostate tumor volume, that is, less than that prescribed, and thereby decrease the tumor control rate. On the other hand, if the organ at risk is not at the same position as in the planning CT image, there is a potential of delivering an excessive dose possibly resulting in more serious complications.
Delivered doses can only be fully assessed using CT scans acquired at the time of treatment. Accurate interfractional patient repositioning before prostate radiotherapy has become possible in this era of image-guided radiotherapy with several in-room systems [6–9]. The KV/MV cone beam CT (CBCT) technique is one of these methods, allowing on-treatment acquisition of 3D images with excellent bony and reasonable soft-tissue definition, at acceptably low radiation doses [8–10]. Information can therefore be acquired for treatment set-up verification, including organ positioning, such that the volume treated and the daily delivered dose may be accurately calculated [11, 12].
Late rectal bleeding is one of the most important late toxicities of radiotherapy for prostate cancer as it can have a major impact on the patient’s quality of life. Thus, separate doses should be delivered to the rectum during the treatment period. Changes to rectal dimensions due to rectal motions may have a substantial impact on excessive normal tissue toxicity in the rectum. Hellebust et al. reported variations in rectal doses given to patients receiving fractionated high dose rate gynecological brachytherapy using CT-based 3D treatment planning and dose volume histograms (DVH) . For linear accelerator treatments, Sripadam et al. acquired CBCT scans for 15 patients with bony anatomy based set-ups and noted differences in rectal doses from the treatment plan, with rectal volumes being both smaller and larger than calculated averages . However, it is difficult to devise an index, represented by doses or DVH parameters, allowing intuitive prediction of complications.
The possibility of combining this information in one parameter reliably expressing the probability of late rectal bleeding, as with the normal tissue complication probability (NTCP), is very attractive . In this study, variations in rectal volumes and dosimetry values including NTCP of interfractional motion during prostate image-guided IMRT (IG-IMRT) using CBCT were estimated to assess their usefulness for evaluating late complications. On-treatment verification CBCT images were used to study the effects of rectal motion during 2D-based IG-IMRT prostate radiotherapy.
2. Materials and Methods
2.1. Patient Data
Twenty patients with clinically localized stage T1c-T3b prostate carcinoma were recruited, without preselection, between January 2010 and March 2014. The median patient age was 72 years, and the median prostate specific antigen level and Gleason score were 8.5 ng·mL−1 and 7.5, respectively. All patients completed their planned treatment course. There were no preparations for these patients after those of the bladder and bowel protocol, aimed at an empty rectum and full bladder, either before the initial planning scan or at any time during the treatment. For treatment planning, patients were scanned on a large bore fan-beam CT scanner with a 2.0 mm slice thickness in the supine position.
2.2. On-Treatment CBCT Acquisition
At each treatment session, the patient was placed in the appropriate position using simulation markers and alignment lasers, and orthogonal kV images were then acquired. These images were used online to match the pelvic bone positions to the initial plan reference images and the patient was moved to the aligned position. No tolerance for shifting was applied, such that all daily measured set-up errors were corrected with online bony 2D matching. Patients remained supine on a flat carbon fiber couch top with vacuum cushions. CBCTs were scheduled after the initial patient set-up to confirm the interfractional motions of the prostate. Acquired images were subsequently analyzed on a weekly basis during the treatment course (8 CBCTs/patient).
2.3. Dose Calculation and Assessment
The rectal volume was contoured from anal verge to the rectosigmoid junction. The planning treatment volume (PTV) was determined using the XiO planning system (Elekta, Stockholm, Sweden) for the scan from the initial planning CT and each on-treatment CBCT. All scans used a margin of 6 mm posterior and 0.7–1 cm in other directions around the clinical target volume to determine the PTV. The PTV was defined as the prostate alone to which a dose of 74 Gy in 37 fractions had been prescribed for 3 patients and the prostate including seminal vesicles to which a dose of 78 Gy in 39 fractions had been prescribed for 10 patients and 80 Gy in 40 fractions for 7 patients. Energy equivalent to an X-ray dose of 10 MV was used for treatment, normalized at the PTV for 95% isodose coverage. This initial plan was then applied to all CBCTs for each patient to calculate the daily dose delivered. Doses calculated using Hounsfield units (HU) from CBCTs have errors due to inaccuracies inherent to the HU values. Since CBCT images suffer from an increased contribution of scatter, resulting in reconstruction errors in these images, thereby leading to dose errors, comparisons were made by not adopting heterogeneity corrections of dose calculations for either the initial planning or the CBCT images, allowing assessment of the actual effects of rectal motions.
For each patient, the on-treatment DVHs for the rectum were exported and compared with the initial plan DVH. Rectal volume and DVH parameters, that is, , , , , and , were calculated and analyzed. Here, and are the maximum and mean dose deposited in rectum, and is the minimal dose deposited in % or cc of the rectal volume.
NTCPs were analyzed using the Lyman-Kutcher-Burman (LKB) model and the relative seriality (RS) model. To account for variations in dose per fraction in different subvolumes of the rectum with changes in fractionation schedules, the total physical dose corresponding to each DVH bin, , was converted into an isoeffective dose in 2 Gy fractions using the following equation : where is the number of fractions. According to the literature, radiation responses are generally determined for 2 Gy fractions, making it good practice to utilize LQED2 rather than the physical dose in calculations of NTCP . The LKB model is defined as follows : where is the dose that, if given uniformly to the entire volume, will lead to the same NTCP as the actual nonuniform dose distribution, TD50 is the uniform dose given to the entire organ that results in a 50% complication risk, is a measure of the slope of the sigmoid curve, is the volume effect parameter, and is the fractional organ volume receiving a dose . The RS model is based on Poisson statistics and accounts for the architecture of the organ through the parameter. The RS model is derived from the ratio of serial subunits to all subunits in the organ. For a heterogeneous dose distribution the complication probability is given by  where () is the complication rate derived from the Poisson model. In-house software was used to estimate rectal NTCP using the LKB model and the RS model, and the model parameters used for the calculations for each rectum are summarized in Table 1 [3–5].
3.1. Rectum Volumes
A list of volumetric values of the rectum for prostate radiation treatment with IG-IMRT by online bony 2D matching techniques is presented in Table 2. These values from the initial planning CT and analyses from on-treatment CBCTs are compared. Across the patient population, the median (range) volume on the planning CT for the rectum was 45.6 cc (31.3–82.0). Rectal volumes show the spread of on-treatment plan values around the prediction based on the initial plan. In several patients, the on-treatment rectal volume is more than twice that of the initial plan.
3.2. DVH Parameters
Differences in DVH parameters, that is, , , , , and , are summarized in Table 3. In this study, for comparison of two values, results are given in terms of relative differences. The relative differences in () between the prediction based on the initial plan () and on-treatment CBCT analyses () are given as follows: Differences in DVH parameters show systematic regularity shift from the initial plan, since 17 patients had median based on the on-treatment CBCT analyses higher than the prediction based on the initial plan, while 2 patients had lower . One patient (patient 13) had median and from on-treatment CBCT analyses higher than the prediction based on the initial plan, while other , that is, , , and , were lower. To illustrate individual patient data, Figure 1 shows the spread of for on-treatment values around the prediction based on the initial plan, as an example. None of the patients show any recognizable variation pattern.
3.3. Rectal NTCP
Table 4 shows NTCP of late rectal bleeding as predicted from the initial plan and on-treatment CBCT analyses in 2D-based IG-IMRT. Calculations, made with the LKB model and the RS model, indicate the risks of rectal complications based on the on-treatment CBCT analyses to be higher than those predicted based on the initial plan. Figures 2 to 4 show the details of NTCP from on-treatment CBCT analyses, spread around the prediction based on the initial plan, in individual patients. None of the patients showed any recognizable pattern in variation or in , with some having a broad range of values, while in others the spread of NTCPs was within a narrow range.
It is well known that, compared to kV CT images, CBCT images suffer from increased scatter . In general, X-ray scatter reduces image contrast, increases image noise, and may introduce reconstruction errors into CBCT images. In addition, the contribution of scatter is patient geometry dependent. Consequently, the CT numbers of a CBCT scan cannot simply be converted to electron density and directly used for dose recomputation, as this may lead to dose errors [17, 18]. Yang et al. have shown differences in the maximum critical structure dose between the results of the Monte Carlo method and the radiotherapy planning system without heterogeneity correction being negligible . For these reasons, we considered using methods without heterogeneity correction for dose calculations to be appropriate for making comparisons between initial planning CT images and on-treatment CBCT images.
None of our patients showed consistent differences in rectal volumes (Table 2); however, the majority of patients showed systematic regularity shift in DVH parameters between the prediction based on the initial plan and on-treatment CBCT analyses (Table 3). No apparent patterns were discernible when making comparisons among patients or within individual patients, in data values (Figure 1). Some patients did exhibit larger variations than others, but, overall, the variations were essentially random. In a previous study, Hatton et al. found that rectal DVH results tended to underestimate the on-treatment rectal dose, with 65% of all on-treatment plans showing higher rectal doses than predicted . Their study demonstrated the generally limited effectiveness of the rectal preparation protocols used, based on CBCT scans. The study by Kupelian et al. included a rectal enema for the initial CT scan, with no specific instructions for bowel preparations prior to treatment, and found similar significant variations in rectal DVH, leading the authors to conclude that the planned dose was a poor predictor of actual doses delivered to the rectum . Pawlowski et al., who proposed that variations in rectal filling are reflected in dosimetric variations, specified no rectal preparation . In their study, no preparation protocol for the rectum was prescribed, resulting in randomness in the rectal filling variation, and the doses delivered to the rectum were unpredictable. Variation in rectal volumes occurs in all patients due to physiological changes. Whatever the factors underlying these rectal volume differences, the current study provides a quantitative assessment of the dosimetric impact of such variations.
Because of the daily variation in rectal filling, the actual delivered radiation doses did not match the treatment plan in the large majority of fractions, leading to the spread of NTCPs from on-treatment CBCT analyses around the prediction based on the initial plan (Figures 2, 3, and 4). In this study, 2 different NTCP models were used to assess endpoints for 3 types of late rectal bleeding. Results for the prediction of NTCP values based on the initial plan showed a tendency for underestimation of the on-treatment NTCP values, with all NTCP models, indicating that, regardless of which models or parameters are used, the same extents and tendencies are observed for all results (Table 4). Reviewing the NTCP results for individual patients, severe late rectal bleeding was rare. It should be noted that such NCTPs are assumed to represent radiation dose deliveries given to patients in a total of 37, 39, or 40 fractions, such that there is a possibility of overestimating the NTCP by on-treatment CBCT analyses. However, overestimating NTCP involving the organs at risk requires adequate attention to radiation safety, and, in addition, it is helpful to clinicians since the NTCP reliably expresses the probability of late rectal bleeding. The possibility of complications can, in fact, more easily be assessed intuitively than with physical doses or DVH parameters. Additionally, the extent of the range of values for patients is not systematically shown by mapping the range of NTCPs in any currently available model. Spread of generalized equivalent uniform dose (gEUD) for on-treatment values around the prediction based on the initial plan is shown in Figure 5 for comparison, indicating similar extents and tendencies to the result of [15, 22]. These findings confirm that it is difficult to predict late toxicity (i.e., rectal bleeding in this study) using only a physical dose index, like one of the DVH parameters and gEUD, as a predictive factor, especially for organs at risk or having parallel components. In a clinical setting and with on-treatment scans acquired during treatment, conformal plans calculated based on the on-treatment scans have the potential to provide accurate predictions of the overall physical doses delivered and to indicate NTCP. Developing the methods allowing feedback of such information to the clinician for implementation of treatment planning recalculation or for quality assurance would be useful . However, the aspect of adding the process of NTCP calculations into adaptive workflows or into the evaluation phase for resimulation and replanning will need further investigation.
The 20 case studies presented herein are large enough sample to reflect the outcomes in our entire patient population. In addition, detailed examination of individual patient data allows recognition of outcomes which could easily be missed in grouped data. As our findings show, this case series was large enough to identify common issues relevant to the accuracy of the initial treatment plan dose prediction and NTCP for the rectum as compared to that actually experienced throughout the treatment course. These observations warrant further investigation.
Employing an online matching protocol for pelvic bone, we assessed the influences of variations in interfraction motions on rectal volumes and dosimetry values, such as DVH parameters and NTCP, and detected no recognizable variation patterns. This approach to assessing interfractional variability in the rectum receiving 2D-based IG-IMRT during radiation treatment for prostate cancer using NTCP models could provide a useful and practical means of dosimetric verification and overcome limitations of the physical dose index.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors are grateful to Mrs. Kazuko Ogawa at Keio University School of Medicine for her support and assistance.
M. J. Zelefsky, E. J. Levin, M. Hunt et al., “Incidence of late rectal and urinary toxicities after three-dimensional conformal radiotherapy and intensity-modulated radiotherapy for localized prostate cancer,” International Journal of Radiation Oncology Biology Physics, vol. 70, no. 4, pp. 1124–1129, 2008.View at: Publisher Site | Google Scholar
R. Varadhan, S. K. Hui, S. Way, and K. Nisi, “Assessing prostate, bladder and rectal doses during image guided radiation therapy—need for plan adaptation?” Journal of applied clinical medical physics, vol. 10, article 2883, no. 3, 2009.View at: Google Scholar
T. R. Willoughby, P. A. Kupelian, J. Pouliot et al., “Target localization and real-time tracking using the Calypso 4D localization system in patients with localized prostate cancer,” International Journal of Radiation Oncology Biology Physics, vol. 65, no. 2, pp. 528–534, 2006.View at: Publisher Site | Google Scholar
J. Boda-Heggemann, F. M. Köhler, B. Küpper et al., “Accuracy of ultrasound-based (BAT) prostate-repositioning: a three-dimensional on-line fiducial-based assessment with cone-beam computed tomography,” International Journal of Radiation Oncology, Biology, Physics, vol. 70, no. 4, pp. 1247–1255, 2008.View at: Publisher Site | Google Scholar
M. Palombarini, S. Mengoli, P. Fantazzini, C. Cadioli, C. Degli Esposti, and G. P. Frezza, “Analysis of inter-fraction setup errors and organ motion by daily kilovoltage cone beam computed tomography in intensity modulated radiotherapy of prostate cancer,” Radiation Oncology, vol. 7, no. 1, article 56, 2012.View at: Publisher Site | Google Scholar
R. Sripadam, J. Stratford, A. M. Henry, A. Jackson, C. J. Moore, and P. Price, “Rectal motion can reduce CTV coverage and increase rectal dose during prostate radiotherapy: a daily cone-beam CT study,” Radiotherapy & Oncology, vol. 90, no. 3, pp. 312–317, 2009.View at: Publisher Site | Google Scholar
J. A. Hatton, P. B. Greer, C. Tang et al., “Does the planning dose-volume histogram represent treatment doses in image-guided prostate radiation therapy? Assessment with cone-beam computerised tomography scans,” Radiotherapy and Oncology, vol. 98, no. 2, pp. 162–168, 2011.View at: Publisher Site | Google Scholar
T. P. Hellebust, E. Dale, A. Skjønsberg, and D. R. Olsen, “Inter fraction variations in rectum and bladder volumes and dose distributions during high dose rate brachytherapy treatment of the uterine cervix investigated by repetitive CT-examinations,” Radiotherapy and Oncology, vol. 60, no. 3, pp. 273–280, 2001.View at: Publisher Site | Google Scholar
J. M. Pawlowski, E. S. Yang, A. W. Malcolm, C. W. Coffey, and G. X. Ding, “Reduction of dose delivered to organs at risk in prostate cancer patients via image-guided radiation therapy,” International Journal of Radiation Oncology Biology Physics, vol. 76, no. 3, pp. 924–934, 2010.View at: Publisher Site | Google Scholar
A. Niemierko, “A generalized concept of equivalent uniform dose (EUD),” Medical Physics, vol. 26, p. 1100, 1999.View at: Google Scholar