 Research
 Open Access
 Published:
Prediction of gastrointestinal toxicity after external beam radiotherapy for localized prostate cancer
Radiation Oncologyvolume 10, Article number: 80 (2015)
Abstract
Background
Gastrointestinal (GI) toxicity is a common effect following radiation therapy (RT) for prostate cancer. Purpose of the present work is to compare two Normal Tissue Complication Probability (NTCP) modelling approaches for prediction of late radioinduced GI toxicity after prostate external beam radiotherapy.
Methods
The study includes 84 prostate cancer patients evaluated for late rectal toxicity after 3D conformal radiotherapy. Median age was 72 years (range 5385). All patients received a total dose of 76 Gy to the prostate gland with daily fractions of 2 Gy. The acute and late radioinduced GI complications were classified according to the RTOG/EORTC scoring system. Rectum dosevolume histograms were extracted for LymanKutcherBurman (LKB) NTCP model fitting using Maximum Likelihood Estimation. The bootstrap method was employed to test the fit robustness. The area under the receiver operating characteristic curve (AUC) was used to evaluate the predictive power of the LKB and to compare it with a multivariate logistic NTCP model previously determined.
Results
At a median followup of 36 months, 42% (35/84) of patients experienced grade 12 (G12) acute GI events while 25% (21/84) of patients developed G12 late GI events. The bestestimate of fitting parameters for LKB NTCP model for mild\moderate GI toxicity resulted to be: D _{ 50 } = 87.3 Gy, m = 0.37 and n = 0.10. Bootstrap result showed that the parameter fit was robust. The AUC values for the LKB and for the multivariate logistic models were 0.60 and 0.75, respectively.
Conclusions
We derived the parameters of the LKB model for mild\moderate GI toxicity prediction and we compared its performance with that of a datadriven multivariate model. Compared to LKB, the multivariate model confirmed a higher predictive power as showed by the AUC values.
Background
To predict the normal tissue complication probability after a radiotherapy (RT) treatment one of the most wellknown and used method is the LymanKutcherBurman (LKB) model [1] based on dosimetric data available from dosevolume histograms.
Several authors have published studies about fitting of LKB model for radiation induced rectal toxicity and many of them considered rectal bleeding (G12) as toxicity endpoint. Burman et al. were the first authors who computed the model parameters suggesting the values D _{ 50 } = 80.0 Gy, m = 0.15, n = 0.12 for late rectal bleeding [2]. More recently, Rancati et al. [3] have found D _{ 50 } = 81.9 Gy, m = 0.19, n = 0.23 and QUANTEC recommended the values D _{ 50 } = 76.9 Gy, m = 0.13, n = 0.09 to apply LKB model for the same endpoint [46]. While many studies have been published about NTCP models for late rectal bleeding, so far very few authors have reported results about moderate/mild toxicity as high stool frequency, loose stools and rectal urgency as endpoints [710].
The traditional NTCP models use only the dose information to predict toxicity. In the last few years several studies have been published concerning the probability of rectal injury depending on many individual factors extracted from clinical information: the drug prescription (antihypertensives and/or anticoagulants), smoking history, previous abdominal surgery, pretreatment morbidities (hypertension, cardiovascular history), diabetes mellitus, presence of acute gastrointestinal toxicity [8,1115]. Defraene et al. demonstrated the benefit of including clinical factors in the predictive power of different NTCP models for all the endpoints considered, i.e. rectal bleeding, high stool frequency and fecal incontinence [8].
In a previous work, we derived a multivariate NTCP model for late GI toxicity from a set of 57 patients treated with RT for localized prostate cancer [15] and we compared its performance with the LKB model using the parameters suggested in literature [7]. In the present work we extended the sample size with an additional dataset up to 84 patients in order to: 1) fit the LKB model deriving the parameters for specific rectal toxicity; 2) validate the multivariate model on an extended cohort of patients and compare its predictive power with the LKB model.
Methods
Clinical and dosimetric data
We retrospectively analyzed the treatment plans of 99 consecutive patients affected by localized prostate adenocarcinoma treated with radiotherapy at the radiation oncology department of the University “Federico II” of Naples. Data of radiationinduced rectal toxicity and treatment characteristics were obtained from clinical reports and matched with the dosimetric data. All participants gave written informed consent and the patient data were analyzed anonymously. This retrospective study was approved by the local Ethics Committee. For 14 out of 99 patients physics and/or followup data were not available while one patient experienced GI disorders before treatment. These patients were excluded from further evaluation. Clinical data included cardiac comorbidities, smoking history, hormonal therapy, drugs prescription, prostate specific antigen (PSA) at diagnosis. Dosimetric data were extracted from the cumulative dosevolume histograms (DVHs).
All patients were treated with full threedimensional radiation treatment planning with a total dose of 76 Gy in 38 daily fractions of 2 Gy. RT was administered with 20 MV photon beams from linear accelerator by conformal radiation technique (CRT) with six field arrangement or by conformal dynamic arc radiation technique (ART). XIO (Elekta CMS) and ERGO (3DLine Medical System) treatment planning systems were used. The above RT techniques are described more in detail in [15].
Treatment planning was based on computed tomography (CT) performed with empty rectum, comfortably filled bladder, and with the patient in prone position using vacuumlocked mattress. Fivemillimeter increment CT slices of the pelvis extending from L4L5 to 2 cm caudal to the bottom of ischial tuberosities were acquired. CT images were electronically transferred to the CT simulation software (Focal Ease 4.2, Elekta CMS) for target and critical organs contouring. Clinical target volume (CTV) included the prostate gland or the prostate gland plus the seminal vesicles. A 1 cm margin was 3D automatically added around the CTV to define the planning target volume (PTV), except at the boundary between the anterior rectal wall and the prostate, where a 0.5 cm margin was used. The rectum delineation was performed on purpose by the same radiation oncologist (M.C.) according to the male Radiation Therapy Oncology Group (RTOG) Normal Pelvis Atlas [16]. The prescription dose was specified at the center of the PTV. Field weightings were adjusted to achieve the 95% of prescription dose to 95% of the PTV.
End point
The definition of acute and late GI toxicity was evaluated according to Radiation Therapy Oncology Group/European Organization for Research and Treatment of Cancer (RTOG/EORTC) criteria [17]. GI treatment toxicity was classified as acute if it occurred during radiotherapy and in the first 3 months thereafter, while it was classified as late GI toxicity if present after 3 months from the completion of treatment. The monitoring of patients for GI toxicity was part of clinical routine, followup visits were planned every 3 months for the first year, then every 6 months for the next 3 years, and yearly thereafter.
Statistical analysis and modelling
Univariate logistic analysis was performed using the Spearman’s rank correlation (R_{s}) coefficient to assess correlation of patient clinical data with late GI toxicity.
Normal tissue complication probability models
In this study, we have analyzed two NTCP models: LymanKutcherBurman (LKB) model and multivariate logistic NTCP model.
According to LKB model with generalized equivalent uniform dose (gEUD) formulation, NTCP is expressed by the following equations [18]:
D _{ 50 } is the value of the dose corresponding to the 50% probability to induce normal tissue complication, the parameter m is inversely proportional to the slope of doseresponse curve, the parameter n can assume values in the range 01 and accounts for volume effect of the organ; v _{ i } is the relative volume that receives the dose d _{ i }, the sum is over all the bins of DVH.
The Maximum Likelihood (ML) method was used to find the bestfit values of the parameters D _{ 50 }, m and n by maximizing the logarithm of the likelihood (LLH):
The sum is over all the patients with different outcome y(i) = 1 and y(i) = 0, i.e. with and without GI toxicity respectively.
The LLH function was numerically maximized by the NelderMead Simplex Method (Matlab implementation: FMINSEARCH function) using an inhouse developed library for Matlab [19].
Ninetyfive percent confidence intervals for parameters estimates were obtained using the profile likelihood method [20]. Following this method, each parameter belonging to the set (D _{ 50 }, m, n) was varied around its ML estimate (optimum LLH) while the other two parameters were fixed at their ML estimate.
The 95% confidence bounds were determined reducing the maximum LLH by one half of the χ^{2} inverse cumulative distribution function associated to 95%, so as to obtain the isolikelihood contours in each Cartesian plane of the parameters space (D _{ 50 } , m, n). In correspondence to the parameters values belonging to the isolikelihood contours, a bundle of NTCP curves was calculated and the 95% confidence region for the model fit was thus estimated [21].
In order to perform an internal validation of the fitting results and to test the fit robustness, the bootstrap method was here employed to determine the spread in ML estimation of NTCP parameters. The bootstrap resampling method works by refitting the NTCP model using the ML estimation to many pseudodatasets which are created by subsampling the input data set (20000 bootstrap resample runs with a number of folds of 80).
The logistic regression model is based on the sigmoidal relationship between dose and response endpoint. The normal tissue complication probability is given by:
where x represents a vector of input variables and g(x) is given by the following equation:
The model order is defined as s, while α and β_{i} are the corresponding set of model coefficients determined by maximizing the probability that the data gave rise to the observed outcomes.
Model’s evaluation and comparison
In a previous analysis of the GI toxicity on a subset (57 patients) of the present dataset we developed a 3variable logistic regression model consisting of the percentage of rectum volume which receives at least a dose of 65 Gy (V65), the use of antihypertensive and/or anticoagulant drugs (AH/AC) and the manifestation of acute GI toxicity. The expression of g(x) with the values of the coefficients is given by:
To assess the correlation of the LKB and the above multivariate model with GI toxicity the R_{s} coefficient was analyzed. To evaluate the discriminating ability of model fits, the receiveroperating characteristic (ROC) analysis was performed and the area under the receiver operating characteristic curve (AUC) was calculated (SPPS Inc., Chicago IL, vs. 18). The discrimination value on the ROC curve, i.e. the cutoff point optimally classifying patients in a binary prediction problem, was determined by Youden’s J statistic [22]. ROC curve comparison was performed by a Z test.
Results
Clinical and dosimetric data were collected for eightyfour patients of which 42% (35/84) had developed acute GI toxicity while 25% (21/84) had developed late GI toxicity. Among cases of acute rectal morbidity, 74% (26/35) were G1 and 26% (9/35) were G2. Of note, no cases of late rectal toxicity greater than grade 1 were reported.
The univariate logistic regression analysis showed a negative relation between antihypertensive/anticoagulant treatments and GI late toxicity, while a positive relation was found between acute GI toxicity and late GI toxicity (Table 1).
With regard to the dosimetric evaluation, in Figure 1 the mean cumulative rectum DVHs for GI late toxicity patients and complicationfree patients were compared: on average, the rectal volume irradiated in the dose range of 2570 Gy for GI late toxicity patients was greater.
LKB fitting
The optimal NTCP parameters values for LKB model resulted to be D _{ 50 } = 87.3 Gy (95% CI 75.9102.2 Gy), m = 0.37 (95% CI 0.260.64), n = 0.10 (95% CI 0.020.26) and the corresponding value of loglikelihood is LLH = 46.3. Figures 2ac illustrate the isolikelihood contours in each Cartesian plane of the parameters space (D _{ 50 }, m, n). In Figure 2d the bundle of NTCP curves corresponding to the 95% confidence interval region for the model fit is plotted.
The optimal parameters were used to calculate the NTCP values as a function of D_{x} doses, being D_{x} the minimum dose to the x volume. The NTCP values and the fractional rectal volume were then correlated using the R_{s} coefficients (Figure 3).
To test the fit robustness we performed a bootstrap method. The mean and the standard deviation of LKB NTCP model parameters obtained for bootstrap samples are D _{ 50 } = 87 Gy (SD = 6 Gy), m = 0.37 (SD = 0.08), n = 0.10 (SD = 0.03). The mean values of m and n parameters are close to the exact fit to the whole patient cohort.
Comparison of the predictive capability of LKB and logistic NTCP models
In a previous work we suggested a threevariable logistic NTCP model, including clinical patientspecific factors [15]. Accordingly, the risk of G12 late GI toxicity increased as V65 increased, it was higher for patients experiencing previous acute toxicity and lower for patients taking antihypertensive and/or anticoagulant drugs. The model exhibited a good predictive performance (AUC = 0.79). When applied to the present extended dataset, the logistic NTCP performance is still good with an AUC value of 0.75 (95% CI 0.6130.891).
For comparison, the AUC and R_{s} values of the multivariate logistic and LKB models are reported in Table 2 and the ROC curves are shown in Figure 4.
The discrimination value on the ROC curve for LKBNTCP model is 26% and for the logistic NTCP model is 39%. According to LKB model, GI toxicity incidence was higher in patients with NTCP ≥ 26% than in those with NTCP < 26% (37.1% vs. 16.3%) (Figures 5.a and b). According to the logistic model GI toxicity occurred more frequently in the group with NTCP ≥ 39% than in those with NTCP < 39% (15/27, 55.6% vs. 6/57, 10.5%).
The ROC curve comparison showed a significant difference in the prediction capability of the two models: the logistic NTCP resulted in being significantly more predictive when compared to LKB NTCP (z > 1.96, P < 0.05).
Discussion
In recent years, technological advances in radiation therapy have allowed to deliver higher prescribed dose for localized prostate cancer [23] and to reduce the risk of severe adverse effects. Consequently, we are observing a change in the toxicity profiles of external beam RT and the focus of toxicity analysis is changing from rectal bleeding to quality of life assessment [24]. Alterations of intestinal motility and peristalsis such as high stool frequency, loose stools and rectal urgency can greatly affect the quality of life of longsurviving patients. As a consequence it is extremely important to model different and specific aspects of rectal toxicity in a robust way.
Many studies report on LKB normal tissue complication probability parameters for rectal bleeding [25,26] while very few data are available on late mild/moderate radioinduced toxicity for the rectum. Gulliford et al. [7] have found the LKB parameters for specific rectal complications observed in clinical practice including stool frequency, loose stools and rectal urgency.
In the present study, we retrospectively reviewed 84 patients undergoing radiotherapy for prostate cancer. The most frequently observed symptoms were high stool frequency, loose stools and rectal urgency while no rectal bleeding was recorded.
Using the mild/moderate rectal late toxicity as endpoint, two different modelling approaches were compared: the LymanKutcherBurman (LKB) and the datadriven multivariate logistic NTCP models. We compared the predictive power of the two models to understand the benefits of a datadriven approach to NTCP modeling.
Indeed in a previous study [15], we have analyzed dosimetric and clinical data of a subset of patients (57) with localized prostate cancer treated with radiation therapy and we have obtained a multivariate normal tissue complication model for the rectum. We compared its predictive performance with the LKB model using a set of parameters proposed in the literature [7]. In order to be fair in comparing two different model philosophies, in the present study based on an extended cohort, we first identified the best set of parameters for the LKB model on our dataset and then we compared its performance with the logistic NTCP model. Figure 4 shows the better performance of multivariate NTCP model (R_{s} = 0.378, p < .001, AUC = 0.751) compared to the LKB NTCP model (R_{s} = 0.143, p < .001, AUC = 0.595). The obtained results confirm the importance of including, besides the dose, clinical factors such as the use of anticoagulant and/or antihypertensive drugs and the appearance of acute toxicity to obtain a robust prediction of the late toxicity risk. Models that take into account relationships among different patientrelated and dosimetric factors may offer a powerful approach to the optimization of risk ascertainment in order to establish tailored strategies for a patient adapted RT.
Using our cohort of patients, we estimated the best fit values for the parameters of the traditional LKB model by the Maximum Likelihood (ML). The LKB fitted parameter values estimated from our cohort resulted to be D _{ 50 } = 87.3 Gy, m = 0.37, and n = 0.10. We compared this result with that of Gulliford et al. for similar endpoints [7]: the m parameter resulted comparable, while we found different values for D _{ 50 } and n values. Indeed, Gulliford et al. report D _{ 50 } value in the range of 54.162.6 Gy, m value in the range of 0.340.6 and n value in the range 0.360.40.
Very recently, Ospina et al. [5] found lower n values (0.003 to 0.06) for all types of late rectal toxicity (grade ≥ 2) showing the serial behavior of rectum. The slight differences in NTCP parameters may depend on several factors, such as the DVH shape, the conformity of the technique, the dosevolume constraints, the dose calculation algorithm and the rectal volume definition [5]. A volumeexponent of n = 0.10 reflects a remarkable influence of the highdose region, which is mostly determined by the shape of the rectum near the prostate. Several studies explore the effects of the setup error and organ motion on the deviations in the dose delivery and the predicted complications uncertainty [2730]. Of note, in this study no motioninclusive dose distribution has been used to take into account the variation between the planned dose distribution to the rectum and the treated dose distribution. A recent study reports the differences in associations for tenesmus and rectal bleeding using the planned over the motioninclusive dose distributions by a simulations of random and/or systematic motion to the planned dose distributions [27]. Except for the uncertainties arising from setup errors, volume organ variability and contouring method, from our results the rectum seems to act as a more “seriallike” organ. The values of D _{ 50 } = 87.3 Gy (CI 95%, 75.9102.2 Gy) are similar to literature values for G12 rectal bleeding [2,13,31,32] and for stool frequency and fecal incontinence of more severity grade [9].
We expected that the dose response for mild grade complications be shifted to lower dose while the dose response for severity grade of toxicity be shifted to higher doses. However, this expectation was not confirmed by the results. Similarly, Michalski et al. [25] found that the LKB parameters relative to grade ≥ 3 late rectal bleeding were broadly similar to those relative to grade ≥ 2 late rectal bleeding. They attributed this observation to the deviation of daily rectal position from the position of the simulation, which would lead to higher cumulative rectal DVHs than planned.
In this framework, AAPM Task Group 166 suggested great caution to use phenomenological models as representative of radiobiological model of the rectum due to the complex and different clinical situations from which they were derived.
An increased dependence on the highest dose for moderate rectal toxicity of grade 12 is in accordance with a smaller n value [33]. This effect is evidenced in Figure 2.b. Our results suggest a serially constructed organ when moderate late toxicity endpoints are considered. As shown in Figure 3, the correlation coefficient for the LKB NTCP values decrease as the fractional volume increases, implying that the highdose parts of the patient’s histograms corresponds to a better representation of the patients’ actual risk. A similar result was already reported by Dale et al. [34] who observed an increased correlation coefficient as n decreases.
Conclusions
In this study, we demonstrated that the multivariate NTCP model deriving from a previous work outperformed the LKB model derived from the same set of patients. In future perspective both models have to be tested on dose distributions including the changes in rectum shape based on multiple CT scans. A dosevolume effect analysis, with and without including rectal motion, will eventually disclose whether the rectal motion model improves the prediction of GI toxicity by the logistic model.
Abbreviations
 GI:

Gastrointestinal
 RT:

Radiation therapy
 NTCP:

Normal tissue complication probability
 LKB:

LymanKutcherBurman
 AUC:

Area under the curve
 G12:

Grade 12
 PSA:

Prostate specific antigen
 DVH:

Dose volume histogram
 CRT:

Conformal radiation technique
 ART:

Dynamic arc radiation technique
 CT:

Computed tomography
 CTV:

Clinical target volume
 PTV:

Planning target volume
 RTOG:

Radiation Therapy Oncology Group
 EORTC:

European Organization for Research and Treatment of Cancer
 R_{s} :

Spearman’s rank correlation
 gEUD:

generalized equivalent uniform dose
 ML:

Maximum Likelihood
 LLH:

Logarithm of the likelihood
 AH/AC:

Antihypertensive and/or anticoagulant
 ROC:

Receiver operator characteristic
 SD:

Standard deviation
References
 1.
Kutcher GJ, Burman C. Calculation of complication probability factors for nonuniform normal tissue irradiation: the effective volume method. Int J Radiat Oncol Biol Phys. 1989;16:1623–30.
 2.
Burman C, Kutcher GJ, Emami B, Goitein M. Fitting of normal tissue tolerance data to an analytic function. Int J Radiat Oncol Biol Phys. 1991;21:123–35.
 3.
Rancati T, Fiorino C, Gagliardi G, Cattaneo GM, Sanguineti G, Casanova Borca V, et al. Fitting late rectal bleeding data using different NTCP models: results from an Italian multicentric study (AIROPROS0101). Radiother Oncol. 2004;73:21–32.
 4.
Liu M, Moiseenko V, Agranovich A, Karvat A, Kwan W, Saleh ZH, et al. Normal Tissue Complication Probability (NTCP) modeling of late rectal bleeding following external beam radiotherapy for prostate cancer: A Test of the QUANTECrecommended NTCP model. Acta Oncol. 2010;49:1040–4.
 5.
Ospina JD, Zhu J, Chira C, Bossi A, Delobel JB, Beckendorf V, et al. Random forests to predict rectal toxicity following prostate cancer radiation therapy. Int J Radiat Oncol Biol Phys. 2014;89:1024–31.
 6.
Benadjaoud MA, Blanchard P, Shwartz B, Champoudry J, Bouaita R, Lefkopoulos D, et al. Functional Data Analysis in NTCP Modeling: A New Method to Explore the Radiation DoseVolume effects. Int J Radiation Oncol Biol Phys. 2014;90:654–63.
 7.
Gulliford SL, Partridge M, Sydes MR, Webb S, Evans PM, Dearnaley DP. Parameters for the Lyman Kutcher Burman (LKB) model of Normal Tissue Complication Probability (NTCP) for specific rectal complications observed in clinical practise. Radiother Oncol. 2012;102:347–51.
 8.
Defraene G, Van den Bergh L, AlMamgani A, Haustermans K, Heemsbergen W, Van Den Heuvel F, et al. The benefits of including clinical factors in rectal normal tissue complication probability modeling after radiotherapy for prostate cancer. Int J Radiat Oncol Biol Phys. 2012;82:1233–42.
 9.
Peeters ST, Hoogeman MS, Heemsbergen WD, Hart AAM, Koper PCM, Lebesque JV. Rectal bleeding, fecal incontinence, and high stool frequency after conformal radiotherapy for prostate cancer: normal tissue complication probability modeling. Int J Radiat Oncol Biol Phys. 2006;66:11–9.
 10.
Fiorino C, Cozzarini C, Vavassori V, Sanguineti G, Bianchi C, Cattaneo GM, et al. Relationships between DVHs and late rectal bleeding after radiotherapy for prostate cancer: analysis of a large group of patients pooled from three institutions. Radiother Oncol. 2002;64:1–12.
 11.
Fellin G, Fiorino C, Rancati T, Vavassori V, Baccolini M, Bianchi C, et al. Clinical and dosimetric predictors of late rectal toxicity after conformal radiation for localized prostate cancer: results of a large multicenter observational study. Radiother Oncol. 2009;93:197–202.
 12.
Fiorino C, Fellin G, Rancati T, Vavassori V, Bianchi C, Borca VC, et al. Clinical and dosimetric predictors of late rectal syndrome after 3DCRT for localized prostate cancer: preliminary results of a multicenter prospective study. Int J Radiat Oncol Biol Phys. 2008;70:1130–7.
 13.
Rancati T, Fiorino C, Fellin G, Vavassori V, Cagna E, Borca VC, et al. Inclusion of clinical risk factors into NTCP modelling of late rectal toxicity after high dose radiotherapy for prostate cancer. Radiother Oncol. 2011;100:124–30.
 14.
Hamstra DA, Stenmark MH, Ritter T, Litzenberg D, Jackson W, Johnson S, et al. Age and comorbid illness are associated with late rectal toxicity following doseescalated radiation therapy for prostate cancer. Int J Radiat Oncol Biol Phys. 2012;85:1246–53.
 15.
Cella L, D'Avino V, Liuzzi R, Conson M, Doria F, Faiella A, et al. Multivariate normal tissue complication probability modeling of gastrointestinal toxicity after external beam radiotherapy for localized prostate cancer. Radiat Oncol. 2013;8:221.
 16.
Gay HA, Barthold HJ, O'Meara E, Bosch WR, El Naqa I, AlLozi R, et al. Pelvic normal tissue contouring guidelines for radiation therapy: a Radiation Therapy Oncology Group consensus panel atlas. Int J Radiat Oncol Biol Phys. 2012;83:e353–62.
 17.
Cox JD, Stetz J, Pajak TF. Toxicity criteria of the Radiation Therapy Oncology Group (RTOG) and the European Organization for Research and Treatment of Cancer (EORTC). Int J Radiat Oncol Biol Phys. 1995;31:1341–6.
 18.
Niemierko A. Reporting and analyzing dose distributions: a concept of equivalent uniform dose. Med Phys. 1997;24:103–10.
 19.
Cella L, Palma G, Deasy JO, Oh JH, Liuzzi R, D'Avino V, et al. Complication Probability Models for RadiationInduced Heart Valvular Dysfunction: Do HeartLung Interactions Play a Role? PLoS One. 2014;9:e111753.
 20.
Venzon DJ, Moolgavkar SH. A Method for Computing ProfileLikelihoodBased Confidence Intervals. J R Stat Soc Ser C (Appl Stat). 1988;37:87–94.
 21.
Semenenko VA, Li XA. LymanKutcherBurman NTCP model parameters for radiation pneumonitis and xerostomia based on combined analysis of published clinical data. Phys Med Biol. 2008;53:737–55.
 22.
Steyerberg EW, Vickers AJ, Cook NR, Gerds T, Gonen M, Obuchowski N, et al. Assessing the performance of prediction models: a framework for traditional and novel measures. Epidemiology. 2010;21:128–38.
 23.
Cella L, Lomax A, Miralbell L. New techniques in hadrontherapy: intensity modulated proton beams. Phys Med. 2001;17(Supp 1):100–2.
 24.
Yamazaki H, Nakamura S, Nishimura T, Yoshida K, Yoshioka Y, Koizumi M, et al. Transitioning from conventional radiotherapy to intensitymodulated radiotherapy for localized prostate cancer: changing focus from rectal bleeding to detailed quality of life analysis. J Radiat Res. 2014;55:1033–47.
 25.
Michalski JM, Gay H, Jackson A, Tucker SL, Deasy JO. Radiation dosevolume effects in radiationinduced rectal injury. Int J Radiat Oncol Biol Phys. 2010;76:S123–9.
 26.
Trojková D, Judas L, Trojek T. Optimizing the parameters of the Lyman–Kutcher–Burman, Källman, and Logit + EUD models for the rectum – a comparison between normal tissue complication probability and clinical data. Radiat Phys Chem. 2014;104:279–82.
 27.
Thor M, Apte A, Deasy JO, Muren LP. Dose/volumeresponse relations for rectal morbidity using planned and simulated motioninclusive dose distributions. Radiother Oncol. 2013;109:388–93.
 28.
Falco MD, D'Andrea M, Fedele D, Barbarino R, Benassi M, Giudice E, et al. Preliminary experience of a predictive model to define rectal volume and rectal dose during the treatment of prostate cancer. Br J Radiol. 2011;84:819–25.
 29.
Mavroidis P, Axelsson S, Hyodynmaa S, Rajala J, Pitkanen MA, Lind BK, et al. Effects of positioning uncertainty and breathing on dose delivery and radiation pneumonitis prediction in breast cancer. Acta Oncol. 2002;41:471–85.
 30.
Hoogeman MS, van Herk M, Yan D, Boersma LJ, Koper PCM, Lebesque JV. A model to simulate daytoday variations in rectum shape. Int J Radiat Oncol Biol Phys. 2002;54:615–25.
 31.
Sohn M, Yan D, Liang J, Meldolesi E, Vargas C, Alber M. Incidence of late rectal bleeding in highdose conformal radiotherapy of prostate cancer using equivalent uniform dosebased and dosevolumebased normal tissue complication probability models. Int J Radiat Oncol Biol Phys. 2007;67:1066–73.
 32.
Tucker SL, Dong L, Bosch WR, Michalski J, Winter K, Mohan R, et al. Late rectal toxicity on RTOG 9406: analysis using a mixture Lyman model. Int J Radiat Oncol Biol Phys. 2010;78:1253–60.
 33.
Deasy JO. Comments on the use of the LymanKutcherBurman model to describe tissue response to non uniform irradiation. Int J Radiat Oncol Biol Phys. 2000;47:1458–60.
 34.
Dale E, Olsen DR, Fossa SD. Normal tissue complication probabilities correlated with late effects in the rectum after prostate conformal radiotherapy. Int J Radiat Oncol Biol Phys. 1999;43:385–91.
Acknowledgements
This work was partially supported by grants from the Italian Ministry for Education, University and Research (MIUR) in the framework of FIRB (RBFR10Q0PT_001 “DROPS”).
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
LC, and RP conceived and designed the study. FD ,MC reviewed patient clinical and dosimetric data. LC, RL, VDA performed modeling and analysed the data. All authors participated in drafting and revising the manuscript. All authors have given their final approval of the manuscript.
Rights and permissions
About this article
Received
Accepted
Published
DOI
Keywords
 NTCP
 Radiation induced rectum toxicity
 LKB model
 Multivariate model
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.