Dose–volume response analyses of late rectal bleeding after radiotherapy for prostate cancer

A comparative study of different parameter estimation methods for predictive models of Normal Tissue Complication Probability (NTCP) of radiation-induced temporal lobe injury following intensity-modulated radiotherapy in nasopharyngeal carcinoma

BACKGROUND

Normal Tissue Complication Probability (NTCP) models predict temporal lobe injury risk post-intensity-modulated radiotherapy in nasopharyngeal carcinoma patients. Optimal parameter estimation methods for NTCP models need refinement.

PURPOSE

To identify optimal method for parameter estimation in Normal Tissue Complication Probability models for temporal lobe injury following intensity-modulated radiotherapy in nasopharyngeal carcinoma patients.

MATERIALS AND METHODS

In this study, all patients underwent curative intensity-modulated radiation therapy at two research centers. Data of temporal lobes from three cohorts [Data-A, n = 278(training set); Data-B, n = 119(external validation set); Data-C, n = 215(internal validation set)]. Five NTCP models were considered, including the Serial Reconstruction Unit (SRU) model, Poisson model, Lyman model, Logit model and Logistic model. Three parameter estimation methods, namely Bayesian estimation (BE), Least Squares Estimation (LSE) and Maximum Likelihood Estimation (MLE), were applied to calibrate the five NTCP models. Area Under Curve (AUC), confusion matrices, dose-response curves were used to compare the performance of the models.

RESULTS

Six hundred twelve patients were enrolled, with 278 in the Data-A; 119 in the Data-B; 215 in the Data-C. The Poisson-NTCP model was evaluated using AUC and R2 values across three parameter estimation methods (BE, LSE, and MLE) on three datasets. The results were as follows: Data-A: BE (AUC: 0.938, R2: 0.953), LSE (0.942, 0.986), MLE (0.940, 0.843); Data-B: BE (0.744, 0.958), LSE (0.743, 0.697), MLE (0.745, 0.857); Data-C: BE (0.867, 0.915), LSE (0.862, 0.916), MLE (0.865, 0.896). Compared with the remaining models, the Poisson-NTCP model based on BE had also better performance of fitting dose-response curve and recall rate, accuracy and specificity of confusion matrix.

CONCLUSION

Bayesian Estimation (BE) is the best parameter estimation method among the three parameter estimation methods. The Poisson-NTCP model based on BE exhibited the best fit to the data in predicting post-IMRT incidence of TLI in NPC.

Combining dosiomics and machine learning methods for predicting severe cardiac diseases in childhood cancer survivors: the French Childhood Cancer Survivor Study

Background

Cardiac disease (CD) is a primary long-term diagnosed pathology among childhood cancer survivors. Dosiomics (radiomics extracted from the dose distribution) have received attention in the past few years to assess better the induced risk of radiotherapy (RT) than standard dosimetric features such as dose-volume indicators. Hence, using the spatial information contained in the dosiomics features with machine learning methods may improve the prediction of CD.

Methods

We considered the 7670 5-year survivors of the French Childhood Cancer Survivors Study (FCCSS). Dose-volume and dosiomics features are extracted from the radiation dose distribution of 3943 patients treated with RT. Survival analysis is performed considering several groups of features and several models [Cox Proportional Hazard with Lasso penalty, Cox with Bootstrap Lasso selection, Random Survival Forests (RSF)]. We establish the performance of dosiomics compared to baseline models by estimating C-index and Integrated Brier Score (IBS) metrics with 5-fold stratified cross-validation and compare their time-dependent error curves.

Results

An RSF model adjusted on the first-order dosiomics predictors extracted from the whole heart performed best regarding the C-index (0.792 ± 0.049), and an RSF model adjusted on the first-order dosiomics predictors extracted from the heart's subparts performed best regarding the IBS (0.069 ± 0.05). However, the difference is not statistically significant with the standard models (C-index of Cox PH adjusted on dose-volume indicators: 0.791 ± 0.044; IBS of Cox PH adjusted on the mean dose to the heart: 0.074 ± 0.056).

Conclusion

In this study, dosiomics models have slightly better performance metrics but they do not outperform the standard models significantly. Quantiles of the dose distribution may contain enough information to estimate the risk of late radio-induced high-grade CD in childhood cancer survivors.

In Silico Comparison of Three Different Beam Arrangements for Intensity-Modulated Proton Therapy for Postoperative Whole Pelvic Irradiation of Prostate Cancer

Background and purpose: Proton therapy has been shown to provide dosimetric benefits in comparison with IMRT when treating prostate cancer with whole pelvis radiation; however, the optimal proton beam arrangement has yet to be established. The aim of this study was to evaluate three different intensity-modulated proton therapy (IMPT) beam arrangements when treating the prostate bed and pelvis in the postoperative setting. Materials and Methods: Twenty-three post-prostatectomy patients were planned using three different beam arrangements: two-field (IMPT2B) (opposed laterals), three-field (IMPT3B) (opposed laterals inferiorly matched to a posterior-anterior beam superiorly), and four-field (IMPT4B) (opposed laterals inferiorly matched to two posterior oblique beams superiorly) arrangements. The prescription was 50 Gy radiobiological equivalent (GyE) to the pelvis and 70 GyE to the prostate bed. Comparisons were made using paired two-sided Wilcoxon signed-rank tests. Results: CTV coverages were met for all IMPT plans, with 99% of CTVs receiving ≥ 100% of prescription doses. All organ at risk (OAR) objectives were met with IMPT3B and IMPT4B plans, while several rectum objectives were exceeded by IMPT2B plans. IMPT4B provided the lowest doses to OARs for the majority of analyzed outcomes, with significantly lower doses than IMPT2B +/- IMPT3B for bladder V30-V50 and mean dose; bowel V15-V45 and mean dose; sigmoid maximum dose; rectum V40-V72.1, maximum dose, and mean dose; femoral head V37-40 and maximum dose; bone V40 and mean dose; penile bulb mean dose; and skin maximum dose. Conclusion: This study is the first to compare proton beam arrangements when treating the prostate bed and pelvis. four-field plans provided better sparing of the bladder, bowel, and rectum than 2- and three-field plans. The data presented herein may help inform the future delivery of whole pelvis IMPT for prostate cancer.

Refinement & validation of rectal wall dose volume objectives for prostate hypofractionation in 20 fractions

•

LRB was correlated to irradiation technique and several % rectal wall cut points.

•

The predictive role of dosimetric variables relates to the irradiation technique.

•

Dose volume objectives for patients treated with IMRT/VMAT are reported.

Background and purpose

Dose-volume objectives for the rectum have been proposed to limit long term toxicity after moderately hypofractionated radiotherapy (MHRT) for localized prostate cancer. The purpose of the present study is to validate and possibly refine dose volume objective for the rectal wall after 20-fraction MHRT.

Materials and methods

All patients treated by 20-fraction MHRT at a single Institution were identified and relative rectal wall (%RW) DVH retrieved. The endpoint of the study is the development of grade 2+ late rectal bleeding (LRB) according to a modified RTOG scale. Clinical and dosimetric predictors of LRB were investigated at both uni- and multi-variable analysis.

Results

293 patients were identified and analyzed. Of them, 35 (12%) developed the endpoint. At univariable analysis, antithrombotic drug usage (yes vs no), technique (3DCRT vs IMRT/VMAT) and several %RW DVH cut-points were significantly correlated with LRB. However, within patients treated by 3DCRT (N = 106), a bi-variable model including anti-thrombotic drug usage and selected %RW dose/volume metrics failed to identify independent dosimetric predictors of LRB. Conversely, within patients treated with intensity modulation (N = 187), the same model showed a progressively higher impact of the percent of RW receiving doses above 40 Gy. Based on this model, we were able to confirm (V32), refine (V60) and identify a novel (V50) cut-point for the %RW.

Conclusion

We recommend the following dose volume objectives for the %RW in order to minimize the risk of LRB after 20-fraction MHRT: V32 ≤ 50%; V50 ≤ 25.8% and V60 ≤ 10%.

Perspectives on making big data analytics work for oncology

Oncology, with its unique combination of clinical, physical, technological, and biological data provides an ideal case study for applying big data analytics to improve cancer treatment safety and outcomes. An oncology treatment course such as chemoradiotherapy can generate a large pool of information carrying the 5Vs hallmarks of big data. This data is comprised of a heterogeneous mixture of patient demographics, radiation/chemo dosimetry, multimodality imaging features, and biological markers generated over a treatment period that can span few days to several weeks. Efforts using commercial and in-house tools are underway to facilitate data aggregation, ontology creation, sharing, visualization and varying analytics in a secure environment. However, open questions related to proper data structure representation and effective analytics tools to support oncology decision-making need to be addressed. It is recognized that oncology data constitutes a mix of structured (tabulated) and unstructured (electronic documents) that need to be processed to facilitate searching and subsequent knowledge discovery from relational or NoSQL databases. In this context, methods based on advanced analytics and image feature extraction for oncology applications will be discussed. On the other hand, the classical p (variables)≫n (samples) inference problem of statistical learning is challenged in the Big data realm and this is particularly true for oncology applications where p-omics is witnessing exponential growth while the number of cancer incidences has generally plateaued over the past 5-years leading to a quasi-linear growth in samples per patient. Within the Big data paradigm, this kind of phenomenon may yield undesirable effects such as echo chamber anomalies, Yule-Simpson reversal paradox, or misleading ghost analytics. In this work, we will present these effects as they pertain to oncology and engage small thinking methodologies to counter these effects ranging from incorporating prior knowledge, using information-theoretic techniques to modern ensemble machine learning approaches or combination of these. We will particularly discuss the pros and cons of different approaches to improve mining of big data in oncology.

An image-based method to quantify biomechanical properties of the rectum in radiotherapy of prostate cancer

BACKGROUND

Gastrointestinal morbidity after radiotherapy (RT) for prostate cancer may be related to the biomechanical properties of the rectum. In this study we present a magnetic resonance imaging (MRI)-based method to quantitate the thickness and elasticity of the rectal wall in prostate cancer patients treated with RT.

MATERIAL AND METHODS

Four patients previously treated with RT for prostate cancer underwent an MRI session with stepwise rectal bag deflation (from a maximum tolerable volume to 0 ml, in 50 ml steps), with a probe inserted inside the bag to monitor the internal rectal pressure. MRIs were acquired using Dixon sequences (4 mm axial slice thickness) at each deflation step. Rectal walls were defined from the recto-sigmoid junction to 3 cm above the anal canal as the space between the inner and outer wall surfaces. The wall thickness was determined and biomechanical properties (strain and stress) were calculated from the pressure measurements and the MRI-segmented rectal walls.

RESULTS

The integral rectal pressure varied for the maximum tolerable volume (range 150-250 ml) across patients and ranged from 1.3 to 4.0 kPa (SD = 1.2 kPa). Wall thickness was found to vary between patients and also across different rectum segments, with a mean (SD) thickness for the different segments at the 50 ml distension volume of 1.8-4.0 (0.6) mm. Stress showed larger variation than strain, with mean (SD) values for the different segments ranging between 1.5 and 7.0 (1.5) kPa.

CONCLUSION

We have developed a method to quantify biomechanical properties of the rectal wall. The resulting rectal wall thickness, strain and stress differed between patients, as well as across different rectal wall sections. These findings could provide guidance in future predictive outcome modelling in order to better understand the rectal dose-volume response relationship.

Clinical effect of multileaf collimator width on the incidence of late rectal bleeding after high-dose intensity-modulated radiotherapy for localized prostate carcinoma

BACKGROUND

Several studies have confirmed a dosimetric advantage associated with use of a smaller leaf in intensity-modulated radiation therapy (IMRT). However, no studies have identified any clinical benefits. We investigated the effect of a smaller multileaf collimator (MLC) width on the onset of late rectal bleeding after high-dose prostate IMRT.

MATERIALS AND METHODS

Two hundred and five prostate cancer patients were treated with a total dose of 78 Gy in 39 fractions by use of a dynamic MLC technique; however, two different MLC were used: a 10-mm-wide device and a 5-mm-wide device. Gastrointestinal toxicity and several clinical factors were assessed.

RESULTS

The 5-year actuarial risk of grade 2 or higher rectal bleeding was 6.9 % for the 10-mm-wide group (n = 132) and 1.8 % for the 5-mm-wide group (n = 73) (p = 0.04). The median estimated rectal doses for the two groups were 55.1 and 50.6 Gy (p < 0.001), respectively. Univariate analysis showed that acute toxicity, rectal V30-60, median rectal dose, normal tissue complication probability (NTCP), and MLC type were significant predictive factors for late rectal toxicity. In multivariate analysis, acute toxicity and NTCP remained significant.

CONCLUSION

In our planning approach for prostate IMRT, a decrease in MLC width from 10 to 5 mm contributed to further rectal dose reduction, which was the most important predictor of late rectal toxicity.

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

PURPOSE

This study compared normal tissue complication probability (NTCP) modeling of chronic gastrointestinal toxicities following prostate cancer treatment for 2 treatment modalities. Possible factors causing discrepancies in optimal NTCP model parameters between 3-dimensional conformal radiation therapy (3D-CRT) and intensity modulated RT (IMRT) were analyzed and discussed, including the impact of patient characteristics, image guidance, toxicity scoring bias, and NTCP model limitations.

METHODS AND MATERIALS

Rectal wall dose-volume histograms of 1115 patients treated for prostate cancer under an adaptive radiation therapy protocol were used to model gastrointestinal toxicity grade ≥2 (according to Common Terminology Criteria for Adverse Events). A total of 457 patients were treated with 3D-CRT and 658 with IMRT. 3D-CRT patients were matched to IMRT patients based on various patient characteristics, using a propensity score-based algorithm. Parameters of the Lyman equivalent uniform dose and cut-off dose logistic regression NTCP models were estimated for the 2 matched treatment modalities and the combined group.

RESULTS

After they were matched, the 3D-CRT and IMRT groups contained 275 and 550 patients with a large discrepancy of 28.7% versus 7.8% toxicities, respectively (P<.001). For both NTCP models, optimal parameters found for the 3D-CRT groups did not fit the IMRT patients well and vice versa. Models developed for the combined data overestimated NTCP for the IMRT patients and underestimated NTCP for the 3D-CRT group.

CONCLUSIONS

Our analysis did not reveal a single definitive cause for discrepancies of model parameters between 3D-CRT and IMRT. Patient characteristics and bias in toxicity scoring, as well as image guidance alone, are unlikely causes of the large discrepancy of toxicities. Whether the cause was inherent to the specific NTCP models used in this study needs to be verified by future investigations. Because IMRT is increasingly used clinically, it is important that appropriate NTCP model parameters are determined for this treatment modality.

Reducing rectal injury during external beam radiotherapy for prostate cancer

Rectal bleeding and faecal incontinence are serious injuries that men with prostate cancer who receive radiotherapy can experience. Although technical advances--including the use of intensity-modulated radiotherapy coupled with image-guided radiotherapy--have enabled the delivery of dose distributions that conform to the shape of the tumour target with steep dose gradients that reduce the dose given to surrounding tissues, radiotherapy-associated toxicity can not be avoided completely. Many large-scale prospective studies have analysed the correlations of patient-related and treatment-related parameters with acute and late toxicity to optimize patient selection and treatment planning. The careful application of dose-volume constraints and the tuning of these constraints to the individual patient's characteristics are now considered the most effective ways of reducing rectal morbidity. Additionally, the use of endorectal balloons (to reduce the margins between the clinical target volume and planning target volume) and the insertion of tissue spacers into the region between the prostate and anterior rectal wall have been investigated as means to further reduce late rectal injury. Finally, some drugs and other compounds are also being considered to help protect healthy tissue. Overall, a number of approaches exist that must be fully explored in large prospective trials to address the important issue of rectal toxicity in prostate cancer radiotherapy.

Dosimetric assessment of prostate cancer patients through principal component analysis (PCA)

The aims of this study were twofold: first, to determine the impact of variance in dose-volume histograms (DVH) on patient-specific toxicity after 2 high-dose fractions in a sample of 22 men with prostate cancer; and second, to compare the effectiveness of traditional DVH analysis and principal component analysis (PCA) in predicting rectum and urethra toxicity. A series of 22 patients diagnosed with prostate adenocarcinoma was treated with 45 Gy external beam and 20 Gy dose rate brachytherapy. Principal component analysis was applied to model the shapes of the rectum and urethra dose-volume histograms. We used logistic regression to measure the correlations between the principal components and the incidence of rectal bleeding and urethra stricture. We also calculated the equivalent uniform dose (EUD) and normal tissue complication probability (NTCP) for the urethra and rectum, and tumor control probability (TCP) for the prostate using BioSuite software. We evaluated their correlations with rectal and urethra toxicity. The rectum DVHs are well described by one principal component (PC1), which accounts for 93.5% of the variance in their shapes. The urethra DVHs are described by two principal components, PC1 and PC2, which account for 94.98% and 3.15% of the variance, respectively. Multivariate exact logistic regression suggests that urethra PC2 is a good predictor of stricture, with Nagelkerke's R2 estimated at 0.798 and a Wald criterion of 5.421 (p &lt; 0.021). The average NTCPs were 0.06% ± 0.04% and 1.25% ± 0.22% for the rectum and urethra, respectively. The average TCP was 85.29% ± 2.28%. This study suggests that principal component analysis can be used to identify the shape variation in dose-volume histograms, and that the principal components can be correlated with the toxicity of a treatment plan based on multivariate analysis. The principal components are also correlated with traditional dosimetric parameters.

Metabolic imaging biomarkers of postradiotherapy xerostomia

PURPOSE

Xerostomia is a major complication of head and neck radiotherapy (RT). Available xerostomia measures remain flawed. [(18)F]fluorodeoxyglucose-labeled positron emission tomography-computed tomography (FDG-PET-CT) is routinely used for staging and response assessment of head and neck cancer. We investigated quantitative measurement of parotid gland FDG uptake as a potential biomarker for post-RT xerostomia.

METHODS AND MATERIALS

Ninety-eight locally advanced head and neck cancer patients receiving definitive RT underwent baseline and post-RT FDG-PET-CT on a prospective imaging trial. A separate validation cohort of 14 patients underwent identical imaging while prospectively enrolled in a second trial collecting sialometry and patient-reported outcomes. Radiation dose and pre- and post-RT standard uptake values (SUVs) for all voxels contained within parotid gland ROI were deformably registered.

RESULTS

Average whole-gland or voxel-by-voxel models incorporating parotid D(Met) (defined as the pretreatment parotid SUV weighted by dose) accurately predicted posttreatment changes in parotid FDG uptake (e.g., fractional parotid SUV). Fractional loss of parotid FDG uptake closely paralleled early parotid toxicity defined by posttreatment salivary output (p < 0.01) and Radiation Therapy Oncology Group/European Organization for Research and Treatment of Cancer xerostomia scores (p < 0.01).

CONCLUSIONS

In this pilot series, loss of parotid FDG uptake was strongly associated with acute clinical post-RT parotid toxicity. D(Met) may potentially be used to guide function-sparing treatment planning. Prospective validation of FDG-PET-CT as a convenient, quantifiable imaging biomarker of parotid function is warranted and ongoing.

Monte Carlo role in radiobiological modelling of radiotherapy outcomes

Radiobiological models are essential components of modern radiotherapy. They are increasingly applied to optimize and evaluate the quality of different treatment planning modalities. They are frequently used in designing new radiotherapy clinical trials by estimating the expected therapeutic ratio of new protocols. In radiobiology, the therapeutic ratio is estimated from the expected gain in tumour control probability (TCP) to the risk of normal tissue complication probability (NTCP). However, estimates of TCP/NTCP are currently based on the deterministic and simplistic linear-quadratic formalism with limited prediction power when applied prospectively. Given the complex and stochastic nature of the physical, chemical and biological interactions associated with spatial and temporal radiation induced effects in living tissues, it is conjectured that methods based on Monte Carlo (MC) analysis may provide better estimates of TCP/NTCP for radiotherapy treatment planning and trial design. Indeed, over the past few decades, methods based on MC have demonstrated superior performance for accurate simulation of radiation transport, tumour growth and particle track structures; however, successful application of modelling radiobiological response and outcomes in radiotherapy is still hampered with several challenges. In this review, we provide an overview of some of the main techniques used in radiobiological modelling for radiotherapy, with focus on the MC role as a promising computational vehicle. We highlight the current challenges, issues and future potentials of the MC approach towards a comprehensive systems-based framework in radiobiological modelling for radiotherapy.

Rectal dose constraints for intensity modulated radiation therapy of the prostate

Prostate dose escalation appears to have resulted in increased cancer control. Such escalation has been made possible by the ability to deliver more conformal treatment that spares normal tissue from the higher radiation doses. The supposition is that this has enabled higher doses, but without an increase in toxicity. The most disabling toxicity in prostate cancer radiotherapy is rectal. We evaluated the current status of conformal radiation and late rectal toxicity with the goal of determining whether reasonable rectal dose and volume constraints can be determined. Although the literature is inexact, we believe that some generalized constraints can be recommended and show that those recommendations are consistent with what is being used at experienced centers.

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

PURPOSE

To determine the predictive risk factors for Grade 2 or worse rectal bleeding after high-dose-rate brachytherapy (HDR-BT) combined with hypofractionated external-beam radiotherapy (EBRT) for prostate cancer using dose-volume histogram analysis.

METHODS AND MATERIALS

The records of 216 patients treated with HDR-BT combined with EBRT were analyzed. The treatment protocols for HDR-BT were 5 Gy × five times in 3 days or 7 Gy × three, 10.5 Gy × two, or 9 Gy × two in 2 days. The EBRT doses ranged from 45 to 51 Gy with a fractional dose of 3 Gy.

RESULTS

In 20 patients Grade 2 or worse rectal bleeding developed, and the cumulative incidence rate was 9% at 5 years. By converting the HDR-BT and EBRT radiation doses into biologic effective doses (BED), the BED(3) at rectal volumes of 5% and 10% in the patients who experienced bleeding were significantly higher than those in the remaining 196 patients. Univariate analysis showed that a higher rectal BED(3-5%) and the use of fewer needles in brachytherapy were correlated with the incidence of bleeding, but BED(3-5%) was found to be the only significant factor on multivariate analysis.

CONCLUSIONS

The radiation dose delivered to small rectal lesions as 5% is important for predicting Grade 2 or worse rectal bleeding after HDR-BT combined with EBRT for prostate cancer.

Comparison of dose response models for predicting normal tissue complications from cancer radiotherapy: application in rat spinal cord

Seven different radiobiological dose-response models have been compared with regard to their ability to describe experimental data. The first four models, namely the critical volume, the relative seriality, the inverse tumor and the critical element models are mainly based on cell survival biology. The other three models: the Lyman (Gaussian distribution), the parallel architecture and the Weibull distribution models are semi-empirical and rather based on statistical distributions. The maximum likelihood estimation was used to fit the models to experimental data and the χ2-distribution, AIC criterion and F-test were applied to compare the goodness-of-fit of the models. The comparison was performed using experimental data for rat spinal cord injury. Both the shape of the dose-response curve and the ability of handling the volume dependence were separately compared for each model. All the models were found to be acceptable in describing the present experimental dataset (p > 0.05). For the white matter necrosis dataset, the Weibull and Lyman models were clearly superior to the other models, whereas for the vascular damage case, the Relative Seriality model seems to have the best performance although the Critical volume, Inverse tumor, Critical element and Parallel architecture models gave similar results. Although the differences between many of the investigated models are rather small, they still may be of importance in indicating the advantages and limitations of each particular model. It appears that most of the models have favorable properties for describing dose-response data, which indicates that they may be suitable to be used in biologically optimized intensity modulated radiation therapy planning, provided a proper estimation of their radiobiological parameters had been performed for every tissue and clinical endpoint.

Dose-volume constraints to reduce rectal side effects from prostate radiotherapy: evidence from MRC RT01 Trial ISRCTN 47772397

PURPOSE

Radical radiotherapy for prostate cancer is effective but dose limited because of the proximity of normal tissues. Comprehensive dose-volume analysis of the incidence of clinically relevant late rectal toxicities could indicate how the dose to the rectum should be constrained. Previous emphasis has been on constraining the mid-to-high dose range (>/=50 Gy). Evidence is emerging that lower doses could also be important.

METHODS AND MATERIALS

Data from a large multicenter randomized trial were used to investigate the correlation between seven clinically relevant rectal toxicity endpoints (including patient- and clinician-reported outcomes) and an absolute 5% increase in the volume of rectum receiving the specified doses. The results were quantified using odds ratios. Rectal dose-volume constraints were applied retrospectively to investigate the association of constraints with the incidence of late rectal toxicity.

RESULTS

A statistically significant dose-volume response was observed for six of the seven endpoints for at least one of the dose levels tested in the range of 30-70 Gy. Statistically significant reductions in the incidence of these late rectal toxicities were observed for the group of patients whose treatment plans met specific proposed dose-volume constraints. The incidence of moderate/severe toxicity (any endpoint) decreased incrementally for patients whose treatment plans met increasing numbers of dose-volume constraints from the set of V30

CONCLUSION

Considering the entire dose distribution to the rectum by applying dose-volume constraints such as those tested here in the present will reduce the incidence of late rectal toxicity.

Collapse

Key Words Collapse

MESH Headings

• Dose-Response Relationship, Radiation
• Gastrointestinal Hemorrhage/etiology
• Humans
• Male
• Odds Ratio
• Prostatic Neoplasms/radiotherapy
• Radiation Injuries/complications
• Radiation Injuries/prevention & control
• Radiotherapy Dosage
• Rectal Diseases/etiology
• Rectum/pathology
• Rectum/radiation effects

Collapse

Grants

• G0501019 Medical Research Council
• MC_U122861330 Medical Research Council
• C46/A2131 Cancer Research UK

Collapse

Affiliation(s)

Sarah L Gulliford Joint Department of Physics, Institute of Cancer Research and Royal Marsden National Health Service Foundation Trust, Sutton, United Kingdom
Kerwyn Foo
Rachel C Morgan
Edwin G Aird
A Margaret Bidmead
Helen Critchley
Philip M Evans
Stefano Gianolini
W Philip Mayles
A Rollo Moore
Beatriz Sánchez-Nieto
Mike Partridge
Matthew R Sydes
Steve Webb
David P Dearnaley

Collapse

Collaborators Collapse

Assessing correlations between the spatial distribution of the dose to the rectal wall and late rectal toxicity after prostate radiotherapy: an analysis of data from the MRC RT01 trial (ISRCTN 47772397)

Many studies have been performed to assess correlations between measures derived from dose-volume histograms and late rectal toxicities for radiotherapy of prostate cancer. The purpose of this study was to quantify correlations between measures describing the shape and location of the dose distribution and different outcomes. The dose to the rectal wall was projected on a two-dimensional map. In order to characterize the dose distribution, its centre of mass, longitudinal and lateral extent, and eccentricity were calculated at different dose levels. Furthermore, the dose-surface histogram (DSH) was determined. Correlations between these measures and seven clinically relevant rectal-toxicity endpoints were quantified by maximally selected standardized Wilcoxon rank statistics. The analysis was performed using data from the RT01 prostate radiotherapy trial. For some endpoints, the shape of the dose distribution is more strongly correlated with the outcome than simple DSHs. Rectal bleeding was most strongly correlated with the lateral extent of the dose distribution. For loose stools, the strongest correlations were found for longitudinal extent; proctitis was most strongly correlated with DSH. For the other endpoints no statistically significant correlations could be found. The strengths of the correlations between the shape of the dose distribution and outcome differed considerably between the different endpoints. Due to these significant correlations, it is desirable to use shape-based tools in order to assess the quality of a dose distribution.

Volumetric modulated arc therapy (VMAT) vs. serial tomotherapy, step-and-shoot IMRT and 3D-conformal RT for treatment of prostate cancer

INTRODUCTION

Volumetric modulated arc therapy (VMAT), a complex treatment strategy for intensity-modulated radiation therapy, may increase treatment efficiency and has recently been established clinically. This analysis compares VMAT against established IMRT and 3D-conformal radiation therapy (3D-CRT) delivery techniques.

METHODS

Based on CT datasets of 9 patients treated for prostate cancer step-and-shoot IMRT, serial tomotherapy (MIMiC), 3D-CRT and VMAT were compared with regard to plan quality and treatment efficiency. Two VMAT approaches (one rotation (VMAT1x) and one rotation plus a second 200 degrees rotation (VMAT2x)) were calculated for the plan comparison. Plan quality was assessed by calculating homogeneity and conformity index (HI and CI), dose to normal tissue (non-target) and D(95%) (dose encompassing 95% of the target volume). For plan efficiency evaluation, treatment time and number of monitor units (MU) were considered.

RESULTS

For MIMiC/IMRT(MLC)/VMAT2x/VMAT1x/3D-CRT, mean CI was 1.5/1.23/1.45/1.51/1.46 and HI was 1.19/1.1/1.09/1.11/1.04. For a prescribed dose of 76 Gy, mean doses to organs-at-risk (OAR) were 50.69 Gy/53.99 Gy/60.29 Gy/61.59 Gy/66.33 Gy for the anterior half of the rectum and 31.85 Gy/34.89 Gy/38.75 Gy/38.57 Gy/55.43 Gy for the posterior rectum. Volumes of non-target normal tissue receiving > or =70% of prescribed dose (53 Gy) were 337 ml/284 ml/482 ml/505 ml/414 ml, for > or =50% (38 Gy) 869 ml/933 ml/1155 ml/1231 ml/1993 ml and for > or =30% (23 Gy) 2819 ml/3414 ml/3340 ml/3438 ml /3061 ml. D(95%) was 69.79 Gy/70.51 Gy/71,7 Gy/71.59 Gy/73.42 Gy. Mean treatment time was 12 min/6 min/3.7 min/1.8 min/2.5 min.

CONCLUSION

All approaches yield treatment plans of improved quality when compared to 3D-conformal treatments, with serial tomotherapy providing best OAR sparing and VMAT being the most efficient treatment option in our comparison. Plans which were calculated with 3D-CRT provided good target coverage but resulted in higher dose to the rectum.

A comparison of HDR brachytherapy and IMRT techniques for dose escalation in prostate cancer: a radiobiological modeling study

A course of one to three large fractions of high dose rate (HDR) interstitial brachytherapy is an attractive alternative to intensity modulated radiation therapy (IMRT) for delivering boost doses to the prostate in combination with additional external beam irradiation for intermediate risk disease. The purpose of this work is to quantitatively compare single-fraction HDR boosts to biologically equivalent fractionated IMRT boosts, assuming idealized image guided delivery (igIMRT) and conventional delivery (cIMRT). For nine prostate patients, both seven-field IMRT and HDR boosts were planned. The linear-quadratic model was used to compute biologically equivalent dose prescriptions. The cIMRT plan was evaluated as a static plan and with simulated random and setup errors. The authors conclude that HDR delivery produces a therapeutic ratio which is significantly better than the conventional IMRT and comparable to or better than the igIMRT delivery. For the HDR, the rectal gBEUD analysis is strongly influenced by high dose DVH tails. A saturation BED, beyond which no further injury can occur, must be assumed. Modeling of organ motion uncertainties yields mean outcomes similar to static plan outcomes.

Modeling of alpha/beta for late rectal toxicity from a randomized phase II study: conventional versus hypofractionated scheme for localized prostate cancer

BACKGROUND

Recently, the use of hypo-fractionated treatment schemes for the prostate cancer has been encouraged due to the fact that alpha/beta ratio for prostate cancer should be low. However a major concern on the use of hypofractionation is the late rectal toxicity, it is important to be able to predict the risk of toxicity for alternative treatment schemes, with the best accuracy. The main purpose of this study is to evaluate the response of rectum wall to changes in fractionation and to quantify the alpha/beta ratio for late rectal toxicity

METHODS

162 patients with localized prostate cancer, treated with conformal radiotherapy, were enrolled in a phase II randomized trial. The patients were randomly assigned to 80 Gy in 40 fractions over 8 weeks (arm A) or 62 Gy in 20 fractions over 5 weeks (arm B). The median follow-up was 30 months. The late rectal toxicity was evaluated using the Radiation Therapy Oncology Group (RTOG) scale. It was assumed >or= Grade 2 (G2) toxicity incidence as primary end point. Fit of toxicity incidence by the Lyman-Burman-Kutcher (LKB) model was performed.

RESULTS

The crude incidence of late rectal toxicity >or= G2 was 14% and 12% for the standard arm and the hypofractionated arm, respectively. The crude incidence of late rectal toxicity >or= G2 was 14.0% and 12.3% for the arm A and B, respectively. For the arm A, volumes receiving >or= 50 Gy (V50) and 70 Gy (V70) were 38.3 +/- 7.5% and 23.4 +/- 5.5%; for arm B, V38 and V54 were 40.9 +/- 6.8% and 24.5 +/- 4.4%. An alpha/beta ratio for late rectal toxicity very close to 3 Gy was found.

CONCLUSION

The >or= G2 late toxicities in both arms were comparable, indicating the feasibility of hypofractionated regimes in prostate cancer. An alpha/beta ratio for late rectal toxicity very close to 3 Gy was found.

Comparison of dose contribution to normal pelvic tissues among conventional, conformal and intensity-modulated radiotherapy techniques in prostate cancer

High-energy external radiotherapy has become one of the most common treatment in localized prostate cancer. We compared the difference of dose distribution, mainly at the 5-30 Gy dose level, in the irradiated pelvic volume among three modalities of radiotherapy for patients with prostate cancer: conventional, conformal and intensity-modulated radiotherapy (IMRT). We selected six patients with prostate cancer treated by conformal radiotherapy at the doses of 46 Gy to PTVN (prostate and seminal vesicles), and 70 Gy to PTV-T (prostate). The conventional technique": an 8-field arrangement was used; the conformal technique 4 fields with a boost through 6 fields. For IMRT, a five-beam arrangement was used. Dose-volume histograms (DVH) were analyzed and compared among the three techniques. The IMRT technique significantly increased the pelvic volume covered by the isodose surfaces below 15 Gy as compared with the conventional and conformal techniques. The mean absolute increase for the pelvic volume included between 5-30 Gy for the IMRT technique, was about 2 900 ml as compared with the conventional technique. However, IMRT significantly reduced the irradiated volume of the rectum in the dose range of 5 to 40 Gy, also significantly reduced the irradiated volume of bladder and femoral heads, and obtained a similar or improved isodose distribution in the PTVs. In addition, the use of IMRT slightly increased the relative dose delivered to the body volume outside the pelvis, as estimated by the use of specific software. A long-term follow-up will be needed to evaluate potential late treatment complications related to the use of IMRT and the low or moderate irradiation dose level obtained in the pelvis and in the whole body.

Radiation-induced rib fractures after hypofractionated stereotactic body radiation therapy of non-small cell lung cancer: a dose- and volume-response analysis

BACKGROUND AND PURPOSE

The aim of this study is to analyse the dose-response and the volume-response of radiation-induced rib fractures after hypofractionated stereotactic body radiation therapy (SBRT).

MATERIALS AND METHODS

During the period 1998-2005, 68 patients with medically inoperable stage I non-small cell lung cancer (NSCLC) were treated with hypofractionated SBRT to 45 Gy in 3 fractions. Among the 33 patients with complete treatment records and radiographic follow-up exceeding 15 months (median: 29 months), 13 fractures were found in seven patients. Identifying all ribs receiving at least 21 Gy, 81 ribs (13 with and 68 without fracture) in 26 patients were separately contoured and their dose-volume histograms (DVHs) were obtained. The DVHs were assessed with the mean dose and cut-off models. Maximum likelihood estimation was used to fit dose-response and volume-response curves to each model.

RESULTS

It was possible to quantify the risk of radiation-induced rib fracture using response curves and information contained in the DVHs. Absolute volumes provided better fits than relative volumes and dose-response curves were more suitable than volume-response curves. For the dose given by the 2 cm(3) cut-off volume, D(2 cm(3)), the logistic dose-response curve for three fractions was parameterised by D(50)=49.8 Gy and gamma(50)=2.05. Consequently, for a median follow-up of 29 months, if D(2 cm(3))<3 x 7.0 Gy the risk is close to 0, and the 5% and 50% risks are given by D(2 cm(3))=3 x 9.1 Gy and 3 x 16.6 Gy, respectively.

CONCLUSIONS

In this group of patients, the risk for radiation-induced rib fracture following hypofractionated SBRT was related to the dose to 2 cm(3) of the rib.

Dosimetric and anatomic indicators of late rectal toxicity after high-dose intensity modulated radiation therapy for prostate cancer

We seek to identify dosimetric and anatomic indicators of late rectal toxicity in prostate cancer patients treated with intensity modulated radiation therapy (IMRT). Data from 49 patients sampled from 698 patients treated for clinically localized prostate cancer at the Memorial Sloan-Kettering Cancer Center with IMRT to a dose of 81 Gy were analyzed. The end point of the study was late Grade 2 or worse rectal toxicity within 30 months of treatment. Dosimetric analysis was performed on the rectum surface in three dimensions and on two-dimensional dose maps obtained by flattening the rectum surface using a conformal mapping procedure. Several parameters including the percentage and absolute surface area of the rectum irradiated, mean dose as a function of location on the rectum, planning target volume (PTV) size and rectum size were analyzed for correlation to toxicity. Significance was set at p < 0.05 for a two-sided t-test. Correlation between absolute areas irradiated and toxicity was observed on both the rectum surface and flattened rectum. Patients with toxicity also received a significantly higher mean dose to the superior 25% of the rectum surface and 15% of the flattened rectum. PTV volume, PTV height, rectum surface area and average cross-sectional area were significantly larger in patients with toxicity. The conformal mapping procedure has potential utility for evaluating dose to the rectum and risk of toxicity. Late rectal toxicity was related to the irradiation of the upper part of the rectum and also to the absolute area irradiated, PTV size, and rectum size on the planning computed tomography (CT) scan.

Principal Component Analysis-Based Pattern Analysis of Dose–Volume Histograms and Influence on Rectal Toxicity

PURPOSE

The variability of dose-volume histogram (DVH) shapes in a patient population can be quantified using principal component analysis (PCA). We applied this to rectal DVHs of prostate cancer patients and investigated the correlation of the PCA parameters with late bleeding.

METHODS AND MATERIALS

PCA was applied to the rectal wall DVHs of 262 patients, who had been treated with a four-field box, conformal adaptive radiotherapy technique. The correlated changes in the DVH pattern were revealed as "eigenmodes," which were ordered by their importance to represent data set variability. Each DVH is uniquely characterized by its principal components (PCs). The correlation of the first three PCs and chronic rectal bleeding of Grade 2 or greater was investigated with uni- and multivariate logistic regression analyses.

RESULTS

Rectal wall DVHs in four-field conformal RT can primarily be represented by the first two or three PCs, which describe approximately 94% or 96% of the DVH shape variability, respectively. The first eigenmode models the total irradiated rectal volume; thus, PC1 correlates to the mean dose. Mode 2 describes the interpatient differences of the relative rectal volume in the two- or four-field overlap region. Mode 3 reveals correlations of volumes with intermediate doses ( approximately 40-45 Gy) and volumes with doses >70 Gy; thus, PC3 is associated with the maximal dose. According to univariate logistic regression analysis, only PC2 correlated significantly with toxicity. However, multivariate logistic regression analysis with the first two or three PCs revealed an increased probability of bleeding for DVHs with more than one large PC.

CONCLUSIONS

PCA can reveal the correlation structure of DVHs for a patient population as imposed by the treatment technique and provide information about its relationship to toxicity. It proves useful for augmenting normal tissue complication probability modeling approaches.

Long-term results of the M. D. Anderson randomized dose-escalation trial for prostate cancer

PURPOSE

To report the long-term results of a randomized radiotherapy dose escalation trial for prostate cancer.

METHODS AND MATERIALS

From 1993 to 1998, a total of 301 patients with stage T1b to T3 prostate cancer were accrued to a randomized external beam dose escalation trial using 70 Gy versus 78 Gy. The median follow-up is now 8.7 years. Kaplan-Meier analysis was used to compute rates of prostate-specific antigen (PSA) failure (nadir + 2), clinical failure, distant metastasis, disease-specific, and overall survival as well as complication rates at 8 years post-treatment.

RESULTS

For all patients, freedom from biochemical or clinical failure (FFF) was superior for the 78-Gy arm, 78%, as compared with 59% for the 70-Gy arm (p = 0.004, and an even greater benefit was seen in patients with initial PSA >10 ng/ml (78% vs. 39%, p = 0.001). The clinical failure rate was significantly reduced in the 78-Gy arm as well (7% vs. 15%, p = 0.014). Twice as many patients either died of prostate cancer or are currently alive with cancer in the 70-Gy arm. Gastrointestinal toxicity of grade 2 or greater occurred twice as often in the high dose patients (26% vs. 13%), although genitourinary toxicity of grade 2 or greater was less (13% vs. 8%) and not statistically significantly different. Dose-volume histogram analysis showed that the complication rate could be significantly decreased by reducing the amount of treated rectum.

CONCLUSIONS

Modest escalation in radiotherapy dose improved freedom from biochemical and clinical progression with the largest benefit in prostate cancer patients with PSA >10 ng/ml.

Self-consistent tumor control probability and normal tissue complication probability models based on generalized EUDa)

Traditional methods to compute the tumor control probability (TCP) or normal tissue complication probability (NTCP) typically require a heterogeneous radiation dose distribution to be converted into a simple uniform dose distribution with an equivalent biological effect. Several power-law type dose-volume-histogram reduction schemes, particularly Niemierko's generalized equivalent uniform dose model [Med. Phys. 26, 1000 (1999)], have been proposed to achieve this goal. In this study, we carefully examine the mathematical outcome of these schemes. We demonstrate that (1) for tumors, with each tumor cell independently responding to local radiation dose, a closed-form analytical solution for tumor survival fraction and TCP can be obtained; (2) for serial structured normal tissues, an exponential power-law form relating survival to functional sub-unit (FSU) radiation is required, and a closed-form analytical solution for the related NTCP is provided; (3) in the case of a parallel structured normal tissue, when NTCP is determined solely by the number of the surviving FSUs, a mathematical solution is available only when there is a non-zero threshold dose and/or a finite critical dose defining the radiotherapy response. Some discussion is offered for the partial irradiation effect on normal tissues in this category; (4) for normal tissues with alternative architectures, where the radiation response of FSU is inhomogeneous, there is no exact global mathematical solution for SF or NTCP within the available schemes. Finally, numerical fits of our models to some experimental data are also presented.

Investigation of bladder dose and volume factors influencing late urinary toxicity after external beam radiotherapy for prostate cancer

BACKGROUND

We sought to identify the bladder dose-volume factors associated with an increased risk of late urinary toxicity among prostate cancer patients treated with radiotherapy.

METHODS AND MATERIALS

This retrospective analysis included data from 128 prostate cancer patients treated on protocol with 2 Gy/fraction to 46 Gy followed by a boost to 78 Gy. The endpoint for this analysis was Grade 1 or greater late genitourinary (GU) toxicity occurring within two years of treatment. The Lyman-Kutcher-Burman, mean dose, threshold dose, and hottest volume models were fitted to the toxicity data using the maximum likelihood method.

RESULTS

Model fits based on dose-volume histograms tended to fit the toxicity data better than models based on dose-wall histograms. The hottest volume (hotspot) model was found to be the best-fitting model investigated. The best fit was for the hottest 2.9% of bladder (95% CI, 1.1-6.8%). This model has an area under the receiver operating characteristic curve of 0.74. The hotspot model separated the patients into clinically meaningful subgroups with approximately 25% of the patients who received <78 Gy to the hottest 2.9% of bladder had GU toxicity at eight years compared with approximately 50% when the dose was > or =78 Gy (p = 0.002).

CONCLUSION

This provides the first evidence supporting that bladder "hotspots" are related to GU toxicity within two years after external beam radiotherapy for prostate cancer. Confirming data are needed from other investigators. Particular attention should be given to hotspots higher than 78 Gy in bladder in radiation treatment planning.

The effect of positional realignment on dose delivery to the prostate and organs-at-risk for 3DCRT

In this study, we evaluate the impact of daily image-guided patient repositioning on dose delivery to prostate and sensitive organs in the treatment of prostate carcinoma with 3-dimensional conformal radiation therapy (3DCRT). Five patients with substantial ultrasound-documented interfractional prostate motion during their 3DCRT treatment course were selected. Starting with the original treatment plan, 2 additional plans were retrospectively generated for each patient. In one set, organ contours were moved for each fraction, thus simulating positioning with misalignment caused by organ motion if ultrasound guidance were not used. In a second set of plans, the isocenter was shifted, as were the organ contours, simulating realignment based on the ultrasound image. In all cases, the number of planned monitor units was set to those of the original plan. For a given patient, isodose distributions, dose-volume histograms (DVHs), equivalent uniform dose (EUD) for prostate, and generalized equivalent uniform dose (gEUDs) for bladder and rectum were calculated for each fraction and then combined for each shift condition. In all reconstructed plans, the results show no substantial changes in dose coverage of the prostate <0.21% change in EUD) compared to the original plan. However, in some cases with no realignment, a larger volume of the bladder or rectum gets higher dose, with the consequent gEUD for each organ significantly greater compared to the original plan.

Incidence of late rectal bleeding in high-dose conformal radiotherapy of prostate cancer using equivalent uniform dose-based and dose-volume-based normal tissue complication probability models

PURPOSE

Accurate modeling of rectal complications based on dose-volume histogram (DVH) data are necessary to allow safe dose escalation in radiotherapy of prostate cancer. We applied different equivalent uniform dose (EUD)-based and dose-volume-based normal tissue complication probability (NTCP) models to rectal wall DVHs and follow-up data for 319 prostate cancer patients to identify the dosimetric factors most predictive for Grade > or = 2 rectal bleeding.

METHODS AND MATERIALS

Data for 319 patients treated at the William Beaumont Hospital with three-dimensional conformal radiotherapy (3D-CRT) under an adaptive radiotherapy protocol were used for this study. The following models were considered: (1) Lyman model and (2) logit-formula with DVH reduced to generalized EUD, (3) serial reconstruction unit (RU) model, (4) Poisson-EUD model, and (5) mean dose- and (6) cutoff dose-logistic regression model. The parameters and their confidence intervals were determined using maximum likelihood estimation.

RESULTS

Of the patients, 51 (16.0%) showed Grade 2 or higher bleeding. As assessed qualitatively and quantitatively, the Lyman- and Logit-EUD, serial RU, and Poisson-EUD model fitted the data very well. Rectal wall mean dose did not correlate to Grade 2 or higher bleeding. For the cutoff dose model, the volume receiving > 73.7 Gy showed most significant correlation to bleeding. However, this model fitted the data more poorly than the EUD-based models.

CONCLUSIONS

Our study clearly confirms a volume effect for late rectal bleeding. This can be described very well by the EUD-like models, of which the serial RU- and Poisson-EUD model can describe the data with only two parameters. Dose-volume-based cutoff-dose models performed worse.

The atlas of complication incidence: a proposal for a new standard for reporting the results of radiotherapy protocols

We present a new method of reporting the results of radiotherapy protocols. The dose-volume atlas of complication incidence is a comprehensive and unbiased summary of the dose-volume exposures and complications occurring in patients after treatment. This new tool provides clear and systematic information about the safety of regions of dose-volume exposure previously treated that can be used when considering new treatments. Actuarial and model-dependent versions of the atlas are described. By using the raw data in the appropriate forms of the atlas, logistic regression, Kaplan-Meier, and Cox proportional hazards analysis can be performed, allowing for the independent calculation of dose-volume response. The data required are simple enough that provided compatible definitions of dose, volume, and complications are used, atlases from different protocols are potentially additive, facilitating the meta-analysis of inter-interinstitutional data. If this method were adopted as a standard for reporting the outcome of treatment protocols, a potentially synergistic increase in the utility of each protocol could result.

Association of DNA repair and steroid metabolism gene polymorphisms with clinical late toxicity in patients treated with conformal radiotherapy for prostate cancer

OBJECTIVE

To explore the possible relationship between single nucleotide polymorphisms (SNP) in candidate genes encoding DNA damage recognition/repair/response and steroid metabolism proteins with respect to clinical radiation toxicity in a retrospective cohort of patients previously treated with three-dimensional conformal radiotherapy (3-DCRT) for prostate cancer.

EXPERIMENTAL DESIGN

One hundred twenty-four patients with prostate cancer underwent 3-DCRT at our institution between September 1996 and December 2000. Of these, 83 consented for follow-up of blood sampling and SNP analysis. Twenty-eight patients were documented as having experienced grade >/=2 late bladder or rectal toxicity (scoring system of Radiation Therapy Oncology Group) on at least one follow-up visit. We analyzed 49 SNPs in BRCA1, BRCA2, ESR1, XRCC1, XRCC2, XRCC3, NBN, RAD51, RAD52, LIG4, ATM, BCL2, TGFB1, MSH6, ERCC2, XPF, NR3C1, CYP1A1, CYP2C9, CYP2C19, CYP3A5, CYP2D6, CYP11B2, and CYP17A1 genes using the Pyrosequencing technique.

RESULTS

Significant univariate associations with late rectal or bladder toxicity (grade >/=2) were found for XRCC3 (A>G 5' untranslated region NT 4541), LIG4 (T>C Asp(568)Asp), MLH1 (C>T, Val(219)Ile), CYP2D6*4 (G>A splicing defect), mean rectal and bladder dose, dose to 30% of rectum or bladder, and age <60 years. On Cox multivariate analysis, significant associations with toxicity were found for LIG4 (T>C, Asp(568)Asp), ERCC2 (G>A, Asp(711)Asp), CYP2D6*4 (G>A, splicing defect), mean bladder dose >60 Gy, and dose to 30% of rectal volume >75 Gy.

CONCLUSIONS

In this study, we identified SNPs in LIG4, ERCC2, and CYP2D6 genes as putative markers to predict individuals at risk for complications arising from radiation therapy in prostate cancer.

Rectal bleeding, fecal incontinence, and high stool frequency after conformal radiotherapy for prostate cancer: normal tissue complication probability modeling

PURPOSE

To analyze whether inclusion of predisposing clinical features in the Lyman-Kutcher-Burman (LKB) normal tissue complication probability (NTCP) model improves the estimation of late gastrointestinal toxicity.

METHODS AND MATERIALS

This study includes 468 prostate cancer patients participating in a randomized trial comparing 68 with 78 Gy. We fitted the probability of developing late toxicity within 3 years (rectal bleeding, high stool frequency, and fecal incontinence) with the original, and a modified LKB model, in which a clinical feature (e.g., history of abdominal surgery) was taken into account by fitting subset specific TD50s. The ratio of these TD50s is the dose-modifying factor for that clinical feature. Dose distributions of anorectal (bleeding and frequency) and anal wall (fecal incontinence) were used.

RESULTS

The modified LKB model gave significantly better fits than the original LKB model. Patients with a history of abdominal surgery had a lower tolerance to radiation than did patients without previous surgery, with a dose-modifying factor of 1.1 for bleeding and of 2.5 for fecal incontinence. The dose-response curve for bleeding was approximately two times steeper than that for frequency and three times steeper than that for fecal incontinence.

CONCLUSIONS

Inclusion of predisposing clinical features significantly improved the estimation of the NTCP. For patients with a history of abdominal surgery, more severe dose constraints should therefore be used during treatment plan optimization.

Cluster model analysis of late rectal bleeding after IMRT of prostate cancer: a case-control study

PURPOSE

Cluster models are newly developed normal-tissue complication probability models in which the spatial aspects of radiation-induced injury are taken into account by considering the size of spatially contiguous aggregates of damaged tissue units. The purpose of this study was to test the validity of a two-dimensional cluster model of late rectal toxicity based on maximum cluster size of damage to rectal surface.

METHODS AND MATERIALS

A paired case-control study was performed in which each of 9 patients experiencing Grade 2 or higher late rectal toxicity after intensity-modulated radiation therapy of localized prostate cancer was paired with a patient having a similar rectal dose-surface histogram but free of rectal toxicity. Numeric simulations were performed to determine the distribution of maximum cluster size on each rectal surface for each of many different choices of possible model parameters.

RESULTS

Model parameters were found for which patients with rectal toxicity were consistently more likely to have a significantly larger mean maximum cluster size than their matched controls. These parameter values correspond to a 50% probability of tissue-unit damage at doses near 30 Gy.

CONCLUSIONS

This study suggests that a cluster model based on maximum cluster size of damage to rectal surface successfully incorporates spatial information beyond that contained in the rectal dose-surface histogram and may therefore provide a useful new tool for predicting rectal normal-tissue complication probability after radiotherapy.

IMRT: Improvement in treatment planning efficiency using NTCP calculation independent of the dose-volume-histogram

The normal tissue complication probability (NTCP) is a predictor of radiobiological effect for organs at risk (OAR). The calculation of the NTCP is based on the dose-volume-histogram (DVH) which is generated by the treatment planning system after calculation of the 3D dose distribution. Including the NTCP in the objective function for intensity modulated radiation therapy (IMRT) plan optimization would make the planning more effective in reducing the postradiation effects. However, doing so would lengthen the total planning time. The purpose of this work is to establish a method for NTCP determination, independent of a DVH calculation, as a quality assurance check and also as a mean of improving the treatment planning efficiency. In the study, the CTs of ten randomly selected prostate patients were used. IMRT optimization was performed with a PINNACLE3 V 6.2b planning system, using planning target volume (PTV) with margins in the range of 2 to 10 mm. The DVH control points of the PTV and OAR were adapted from the prescriptions of Radiation Therapy Oncology Group protocol P-0126 for an escalated prescribed dose of 82 Gy. This paper presents a new model for the determination of the rectal NTCP (R(NTCP)). The method uses a special function, named GVN (from Gy, Volume, NTCP), which describes the R(NTCP) if 1 cm3 of the volume of intersection of the PTV and rectum (R(int)) is irradiated uniformly by a dose of 1 Gy. The function was "geometrically" normalized using a prostate-prostate ratio (PPR) of the patients' prostates. A correction of the R(NTCP) for different prescribed doses, ranging from 70 to 82 Gy, was employed in our model. The argument of the normalized function is the R(int), and parameters are the prescribed dose, prostate volume, PTV margin, and PPR. The R(NTCPs) of another group of patients were calculated by the new method and the resulting difference was < +/- 5% in comparison to the NTCP calculated by the PINNACLE3 software where Kutcher's dose-response model for NTCP calculation is adopted.

Multivariable modeling of radiotherapy outcomes, including dose–volume and clinical factors

PURPOSE

The probability of a specific radiotherapy outcome is typically a complex, unknown function of dosimetric and clinical factors. Current models are usually oversimplified. We describe alternative methods for building multivariable dose-response models.

METHODS

Representative data sets of esophagitis and xerostomia are used. We use a logistic regression framework to approximate the treatment-response function. Bootstrap replications are performed to explore variable selection stability. To guard against under/overfitting, we compare several analytical and data-driven methods for model-order estimation. Spearman's coefficient is used to evaluate performance robustness. Novel graphical displays of variable cross correlations and bootstrap selection are demonstrated.

RESULTS

Bootstrap variable selection techniques improve model building by reducing sample size effects and unveiling variable cross correlations. Inference by resampling and Bayesian approaches produced generally consistent guidance for model order estimation. The optimal esophagitis model consisted of 5 dosimetric/clinical variables. Although the xerostomia model could be improved by combining clinical and dose-volume factors, the improvement would be small.

CONCLUSIONS

Prediction of treatment response can be improved by mixing clinical and dose-volume factors. Graphical tools can mitigate the inherent complexity of multivariable modeling. Bootstrap-based variable selection analysis increases the reliability of reported models. Statistical inference methods combined with Spearman's coefficient provide an efficient approach to estimating optimal model order.

Radiation dosimetry predicts IQ after conformal radiation therapy in pediatric patients with localized ependymoma

PURPOSE

To assess the effects of radiation dose-volume distribution on the trajectory of IQ development after conformal radiation therapy (CRT) in pediatric patients with ependymoma.

METHODS AND MATERIALS

The study included 88 patients (median age, 2.8 years +/- 4.5 years) with localized ependymoma who received CRT (54-59.4 Gy) that used a 1-cm margin on the postoperative tumor bed. Patients were evaluated with tests that included IQ measures at baseline (before CRT) and at 6, 12, 24, 36, 48, and 60 months. Differential dose-volume histograms (DVH) were derived for total-brain, supratentorial-brain, and right and left temporal-lobe volumes. The data were partitioned into three dose intervals and integrated to create variables that represent the fractional volume that received dose over the specified intervals (e.g., V(0-20 Gy), V(20-40 Gy), V(40-65 Gy)) and modeled with clinical variables to develop a regression equation to estimate IQ after CRT.

RESULTS

A total of 327 IQ tests were performed in 66 patients with infratentorial tumors and 20 with supratentorial tumors. The median follow-up was 29.4 months. For all patients, IQ was best estimated by age (years) at CRT; percent volume of the supratentorial brain that received doses between 0 and 20 Gy, 20 and 40 Gy, and 40 and 65 Gy; and time (months) after CRT. Age contributed significantly to the intercept (p > 0.0001), and the dose-volume coefficients were statistically significant (V(0-20 Gy), p = 0.01; V(20-40 Gy), p < 0.001; V(40-65 Gy), p = 0.04). A similar model was developed exclusively for patients with infratentorial tumors but not supratentorial tumors.

CONCLUSION

Radiation dosimetry can be used to predict IQ after CRT in patients with localized ependymoma. The specificity of models may be enhanced by grouping according to tumor location.

Increased risk of biochemical and local failure in patients with distended rectum on the planning CT for prostate cancer radiotherapy

PURPOSE

To retrospectively test the hypothesis that rectal distension on the planning computed tomography (CT) scan is associated with an increased risk of biochemical and local failure among patients irradiated for prostate carcinoma when a daily repositioning technique based on direct prostate-organ localization is not used.

METHODS AND MATERIALS

This study included 127 patients who received definitive three-dimensional conformal radiotherapy for prostate cancer to a total dose of 78 Gy at The University of Texas M. D. Anderson Cancer Center. Rectal distension was assessed by calculation of the average cross-sectional rectal area (CSA; defined as the rectal volume divided by length) and measuring three rectal diameters on the planning CT. The impact of rectal distension on biochemical control, 2-year prostate biopsy results, and incidence of Grade 2 or greater late rectal bleeding was assessed.

RESULTS

The incidence of biochemical failure was significantly higher among patients with distended rectums (CSA >11.2 cm(2)) on the planning CT scan (p = 0.0009, log-rank test). Multivariate analysis indicates that rectal distension and high-risk disease are independent risk factors for biochemical failure, with hazard ratios of 3.89 (95% C.I. 1.58 to 9.56, p = 0.003) and 2.45 (95% C.I. 1.18 to 5.08, p = 0.016), respectively. The probability of residual tumor without evidence of radiation treatment (as scored by the pathologist) increased significantly with rectal distension (p = 0.010, logistic analysis), and a lower incidence of Grade 2 or greater late rectal bleeding within 2 years was simultaneously observed with higher CSA values (p = 0.031, logistic analysis).

CONCLUSIONS

We found strong evidence that rectal distension on the treatment-planning CT scan decreased the probability of biochemical control, local control, and rectal toxicity in patients who were treated without daily image-guided prostate localization, presumably because of geographic misses. Therefore, an empty rectum is warranted at the time of simulation. These results also emphasize the need for image-guided radiotherapy to improve local control in irradiating prostate cancer.

Intensity-modulated radiation therapy: emerging cancer treatment technology

The use of intensity-modulated radiation therapy (IMRT) is rapidly advancing in the field of radiation oncology. Intensity-modulated radiation therapy allows for improved dose conformality, thereby affording the potential to decrease the spectrum of normal tissue toxicities associated with IMRT. Preliminary results with IMRT are quite promising; however, the clinical data is relatively immature and overall patient numbers remain small. High-quality IMRT requires intensive physics support and detailed knowledge of three-dimensional anatomy and patterns of tumour spread. This review focuses on basic principles, and highlights the clinical implementation of IMRT in head and neck and prostate cancer.

IMRT versus conventional 3DCRT on prostate and normal tissue dosimetry using an endorectal balloon for prostate immobilization

This study was undertaken to compare prostate and normal tissue dosimetry in prostate cancer patients treated with intensity-modulated radiation therapy (IMRT) and conventional 3-dimensional conformal radiotherapy (3DCRT) using an endorectal balloon for prostate immobilization. Ten prostate cancer patients were studied using both IMRT and conventional 3DCRT at Houston Veterans Affairs Medical Center. For IMRT, the prescription was 70 Gy at 2 Gy/fraction at the 83.4% isodose line, allowing no more than 15% of the rectum and 33% of the bladder to receive above 68 and 65 Gy, respectively. For conventional 3DCRT, a 6-field arrangement with lateral and oblique fields was used to deliver 76 Gy at 2Gy/fraction, ensuring complete tumor coverage by the 72-Gy isodose line. Mean doses for prostate and seminal vesicles were 75.10 and 75.11 Gy, respectively, for IMRT and 75.40 and 75.02 Gy, respectively, for 3DCRT (p > 0.218). 3DCRT delivered significantly higher doses to 33%, 50%, and 66% volumes of rectum by 3.55, 6.64, and 10.18 Gy, respectively (p < 0.002), and upper rectum by 7.26, 9.86, and 9.16 Gy, respectively (p < 0.007). 3DCRT also delivered higher doses to femur volumes of 33% and 50% by 9.38 and 10.19 Gy, respectively, (p < 0.001). Insignificant differences in tumor control probability (TCP) values between IMRT and 3DCRT were calculated for prostate (p = 0.320) and seminal vesicles (p = 0.289). Compared to 3DCRT, IMRT resulted in significantly reduced normal tissue complication probability (NTCP) only for upper rectum (p = 0.025) and femurs (p = 0.021). This study demonstrates that IMRT achieves superior normal tissue avoidance, especially for rectum and femurs compared to 3DCRT, with comparable target dose escalation.

Implementation and validation of a three-dimensional deformable registration algorithm for targeted prostate cancer radiotherapy

PURPOSE

Daily prostate deformation hinders accurate calculation of dose, especially to intraprostatic targets. We implemented a three-dimensional deformable registration algorithm to aid dose tracking for targeted prostate radiotherapy.

METHODS AND MATERIALS

The algorithm registers two computed tomography (CT) scans by iteratively minimizing their differences in image intensity. For validation, we measured the accuracy in registering (a) a pelvic CT set to its mathematically deformed counterpart, (b) CT scans of a deformable pelvic phantom with and without an endorectal balloon inflated, to simulate intraprostatic targets, 23 CT-opaque seeds were embedded in the prostate, and (c) two pelvic CT scans of a patient obtained on 2 separate days.

RESULTS

The mean (SD) error in registering the pelvic CT set to its transformed set was 0.5 mm (1.5), with correlation coefficient improvement from 0.626 to 0.991. Using the deformable pelvic phantom, the correlation coefficient improved from 0.543 to 0.816 after registration. The mean (SD) error in tracking the intraprostatic seeds was 0.8 mm (0.5). The correlation coefficient improved from 0.610 to 0.944 after registration of the two patient CT sets.

CONCLUSION

The algorithm had an accuracy of about 1 mm. It could be used for optimizing dose calculation and delivery for prostate radiotherapy.