First application of a pixel-wise analysis on bladder dose-surface maps in prostate cancer radiotherapy

More than one way to skin a dose volume: the impact of dose-surface map calculation approach on study reproducibility

Objective.Dose-surface maps (DSMs) provide spatial representations of the radiation dose to organ surfaces during radiotherapy and are a valuable tool for identifying dose deposition patterns that are predictive of radiation toxicities. Over the years, many different DSM calculation approaches have been introduced and used in dose-outcome studies. However, little consideration has been given to how these calculation approaches may be impacting the reproducibility of studies in the field. Therefore, we conducted an investigation to determine the level of equivalence of DSMs calculated with different approaches and their subsequent impact on study results.Approach.Rectum and bladder DSMs were calculated for 20 prostate radiotherapy patients using combinations of the most common slice orientation and spacing styles in the literature. Equivalence of differently calculated DSMs was evaluated using pixel-wise comparisons and DSM features (rectum only). Finally, mock cohort comparison studies were conducted with DSMs calculated using each approach to determine the level of dosimetric study reproducibility between calculation approaches.Main results.We found that rectum DSMs calculated using the planar and non-coplanar orientation styles were non-equivalent in the posterior rectal region and that equivalence of DSMs calculated with different slice spacing styles was conditional on the choice of inter-slice distance used. DSM features were highly sensitive to choice of slice orientation style and DSM sampling resolution. Finally, while general result trends were consistent between the comparison studies performed using different DSMs, statisitically significant subregions and features could vary greatly in position and magnitude.Significance.We have determined that DSMs calculated with different calculation approaches are frequently non-equivalent and can lead to differing conclusions between studies performed using the same dataset. We recommend that the DSM research community work to establish consensus calculation approaches to ensure reproducibility within the field.

Patient-reported persistent symptoms after radiotherapy and association with quality of life for prostate cancer survivors

PURPOSE

To evaluate the persistence of symptoms after radiotherapy (RT) for localised prostate cancer (PCa) and the association with quality of life (QOL).

MATERIALS AND METHODS

Prospective patient-reported outcome (PRO) from a multi-institutional study on PCa treated with radical RT (2010-2014) was analysed. Data was collected at baseline (BL) and follow-ups (FUPs) up to 5 years. Patients with BL and ≥3 late FUPs (≥6 months) were analysed. PRO was scored by means of the IPSS and ICIQ-SF (urinary), LENT-SOMA (gastrointestinal [GI]), and EORTC-C30 (pain, insomnia, fatigue, and QOL) questionnaires. Symptoms were defined 'persistent' if the median score over FUPs was ≥3 (urinary) or ≥2 (GI, pain, insomnia, and fatigue), and worse than BL. Different thresholds were chosen to have enough events for each symptom. QOL was linearly transformed on a continuous scale (0-100). Linear-mixed models were used to identify significant differences between groups with and without persistent symptoms including age, smoking status, previous abdominal surgery, and diabetes as confounders. Mean QOL differences between groups were evaluated longitudinally over FUPs.

RESULTS

The analysis included 293 patients. Persistent urinary symptoms ranged from 2% (straining) to 12% (weak stream, and nocturia). Gastrointestinal symptoms ranged from 7% (rectal pain, and incontinence) to 30% (urgency). Proportions of pain, insomnia, and fatigue were 6, 13, and 18%. Significant QOL differences of small-to-medium clinical relevance were found for urinary incontinence, frequency, urgency, and nocturia. Among GI symptoms, rectal pain and incontinence showed small-to-medium differences. Fatigue was associated with the largest differences.

CONCLUSIONS

The analysis showed that symptoms after RT for PCa occur with different persistence and their association with QOL varies in magnitude. A number of persistent urinary and GI symptoms showed differences in a comparable range. Urinary incontinence and frequency, rectal pain, and faecal incontinence more often had significant associations. Fatigue was also prevalent and associated with largely deteriorated QOL.

Selective sparing of bladder and rectum sub-regions in radiotherapy of prostate cancer combining knowledge-based automatic planning and multicriteria optimization

Background and Purpose

The association between dose to selected bladder and rectum symptom-related sub-regions (SRS) and late toxicity after prostate cancer radiotherapy has been evidenced by voxel-wise analyses. The aim of the current study was to explore the feasibility of combining knowledge-based (KB) and multi-criteria optimization (MCO) to spare SRSs without compromising planning target volume (PTV) dose delivery, including pelvic-node irradiation.

Materials and Methods

Forty-five previously treated patients (74.2 Gy/28fr) were selected and SRSs (in the bladder, associated with late dysuria/hematuria/retention; in the rectum, associated with bleeding) were generated using deformable registration. A KB model was used to obtain clinically suitable plans (KB-plan). KB-plans were further optimized using MCO, aiming to reduce dose to the SRSs while safeguarding target dose coverage, homogeneity and avoiding worsening dose volume histograms of the whole bladder, rectum and other organs at risk. The resulting MCO-generated plans were examined to identify the best-compromise plan (KB + MCO-plan).

Results

The mean SRS dose decreased in almost all patients for each SRS. D1% also decreased in the large majority, less frequently for dysuria/bleeding SRS. Mean differences were statistically significant (p < 0.05) and ranged between 1.3 and 2.2 Gy with maximum reduction of mean dose up to 3-5 Gy for the four SRSs. The better sparing of SRSs was obtained without compromising PTVs coverage.

Conclusions

Selectively sparing SRSs without compromising PTV coverage is feasible and has the potential to reduce toxicities in prostate cancer radiotherapy. Further investigation to better quantify the expected risk reduction of late toxicities is warranted.

Experimental validation of a novel method of dose accumulation for the rectum

BACKGROUND

Dose-surface maps (DSMs) are an increasingly popular tool to evaluate spatial dose-outcome relationships for the rectum. Recently, DSM addition has been proposed as an alternative method of dose accumulation from deformable registration-based techniques. In this study, we performed the first experimental investigation of the accuracy at which DSM accumulation can capture the total dose delivered to a rectum's surface in the presence of inter-fraction motion.

MATERIAL AND METHODS

A custom PVC rectum phantom capable of representing typical rectum inter-fraction motion and filling variations was constructed for this project. The phantom allowed for the placement of EBT3 film sheets on the representative rectum surface to measure rectum surface dose. A multi-fraction prostate VMAT treatment was designed and delivered to the phantom in a water tank for a variety of inter-fraction motion scenarios. DSMs for each fraction were calculated in two ways using CBCT images acquired during delivery and summed to produce accumulated DSMs. Accumulated DSMs were then compared to film measurements using gamma analysis (3%/2 mm criteria). Similarity of isodose clusters between films and DSMs was also investigated.

RESULTS

Baseline agreement between film measurements and accumulated DSMs for a stationary rectum was 95.6%. Agreement between film and accumulated DSMs in the presence of different types of inter.-fraction motion was ≥92%, and isodose cluster mean distance to agreement was within 1.5 mm for most scenarios. Overall, DSM accumulation performed the best when using DSMs that accounted for changes in rectum path orientation.

CONCLUSION

Dose accumulation performed with DSMs was found to accurately replicate total delivered dose to a rectum phantom in the presence of inter-fraction motion.

Predicting cancer relapse following lung stereotactic radiotherapy: an external validation study using real-world evidence

Purpose

For patients receiving lung stereotactic ablative radiotherapy (SABR), evidence suggests that high peritumor density predicts an increased risk of microscopic disease (MDE) and local-regional failure, but only if there is low or heterogenous incidental dose surrounding the tumor (GTV). A data-mining method (Cox-per-radius) has been developed to investigate this dose-density interaction. We apply the method to predict local relapse (LR) and regional failure (RF) in patients with non-small cell lung cancer.

Methods

199 patients treated in a routine setting were collated from a single institution for training, and 76 patients from an external institution for validation. Three density metrics (mean, 90th percentile, standard deviation (SD)) were studied in 1mm annuli between 0.5cm inside and 2cm outside the GTV boundary. Dose SD and fraction of volume receiving less than 30Gy were studied in annuli 0.5-2cm outside the GTV to describe incidental MDE dosage. Heat-maps were created that correlate with changes in LR and RF rates due to the interaction between dose heterogeneity and density at each distance combination. Regions of significant improvement were studied in Cox proportional hazards models, and explored with and without re-fitting in external data. Correlations between the dose component of the interaction and common dose metrics were reported.

Results

Local relapse occurred at a rate of 6.5% in the training cohort, and 18% in the validation cohort, which included larger and more centrally located tumors. High peritumor density in combination with high dose variability (0.5 - 1.6cm) predicts LR. No interactions predicted RF. The LR interaction improved the predictive ability compared to using clinical variables alone (optimism-adjusted C-index; 0.82 vs 0.76). Re-fitting model coefficients in external data confirmed the importance of this interaction (C-index; 0.86 vs 0.76). Dose variability in the 0.5-1.6 cm annular region strongly correlates with heterogeneity inside the target volume (SD; ρ = 0.53 training, ρ = 0.65 validation).

Conclusion

In these real-world cohorts, the combination of relatively high peritumor density and high dose variability predicts increase in LR, but not RF, following lung SABR. This external validation justifies potential use of the model to increase low-dose CTV margins for high-risk patients.

Technical note: rtdsm-An open-source software for radiotherapy dose-surface map generation and analysis

BACKGROUND

Dose-outcome studies in radiation oncology have historically excluded spatial information due to dose-volume histograms being the most dominant source of dosimetric information. In recent years, dose-surface maps (DSMs) have become increasingly popular for characterization of spatial dose distributions and identification of radiosensitive subregions for hollow organs. However, methodological variations and lack of open-source, publicly offered code-sharing between research groups have limited reproducibility and wider adoption.

PURPOSE

This paper presents rtdsm, an open-source software for DSM calculation with the intent to improve the reproducibility of and the access to DSM-based research in medical physics and radiation oncology.

METHODS

A literature review was conducted to identify essential functionalities and prevailing calculation approaches to guide development. The described software has been designed to calculate DSMs from DICOM data with a high degree of user customizability and to facilitate DSM feature analysis. Core functionalities include DSM calculation, equivalent dose conversions, common DSM feature extraction, and simple DSM accumulation.

RESULTS

A number of use cases were used to qualitatively and quantitatively demonstrate the use and usefulness of rtdsm. Specifically, two DSM slicing methods, planar and noncoplanar, were implemented and tested, and the effects of method choice on output DSMs were demonstrated. An example comparison of DSMs from two different treatments was used to highlight the use cases of various built-in analysis functions for equivalent dose conversion and DSM feature extraction.

CONCLUSIONS

We developed and implemented rtdsm as a standalone software that provides all essential functionalities required to perform a DSM-based study. It has been made freely accessible under an open-source license on Github to encourage collaboration and community use.

High weekly integral dose and larger fraction size increase risk of fatigue and worsening of functional outcomes following radiotherapy for localized prostate cancer

Introduction

We hypothesized that increasing the pelvic integral dose (ID) and a higher dose per fraction correlate with worsening fatigue and functional outcomes in localized prostate cancer (PCa) patients treated with external beam radiotherapy (EBRT).

Methods

The study design was a retrospective analysis of two prospective observational cohorts, REQUITE (development, n=543) and DUE-01 (validation, n=228). Data were available for comorbidities, medication, androgen deprivation therapy, previous surgeries, smoking, age, and body mass index. The ID was calculated as the product of the mean body dose and body volume. The weekly ID accounted for differences in fractionation. The worsening (end of radiotherapy versus baseline) of European Organisation for Research and Treatment of Cancer EORTC) Quality of Life Questionnaire (QLQ)-C30 scores in physical/role/social functioning and fatigue symptom scales were evaluated, and two outcome measures were defined as worsening in ≥2 (WS2) or ≥3 (WS3) scales, respectively. The weekly ID and clinical risk factors were tested in multivariable logistic regression analysis.

Results

In REQUITE, WS2 was seen in 28% and WS3 in 16% of patients. The median weekly ID was 13.1 L·Gy/week [interquartile (IQ) range 10.2-19.3]. The weekly ID, diabetes, the use of intensity-modulated radiotherapy, and the dose per fraction were significantly associated with WS2 [AUC (area under the receiver operating characteristics curve) =0.59; 95% CI 0.55-0.63] and WS3 (AUC=0.60; 95% CI 0.55-0.64). The prevalence of WS2 (15.3%) and WS3 (6.1%) was lower in DUE-01, but the median weekly ID was higher (15.8 L·Gy/week; IQ range 13.2-19.3). The model for WS2 was validated with reduced discrimination (AUC=0.52 95% CI 0.47-0.61), The AUC for WS3 was 0.58.

Conclusion

Increasing the weekly ID and the dose per fraction lead to the worsening of fatigue and functional outcomes in patients with localized PCa treated with EBRT.

Radial Data Mining to Identify Density-Dose Interactions That Predict Distant Failure Following SABR

Purpose

Lower dose outside the planned treatment area in lung stereotactic radiotherapy has been linked to increased risk of distant metastasis (DM) possibly due to underdosage of microscopic disease (MDE). Independently, tumour density on pretreatment computed tomography (CT) has been linked to risk of MDE. No studies have investigated the interaction between imaging biomarkers and incidental dose. The interaction would showcase whether the impact of dose on outcome is dependent on imaging and, hence, if imaging could inform which patients require dose escalation outside the gross tumour volume (GTV). We propose an image-based data mining methodology to investigate density-dose interactions radially from the GTV to predict DM with no a priori assumption on location.

Methods

Dose and density were quantified in 1-mm annuli around the GTV for 199 patients with early-stage lung cancer treated with 60 Gy in 5 fractions. Each annulus was summarised by three density and three dose parameters. For parameter combinations, Cox regressions were performed including a dose-density interaction in independent annuli. Heatmaps were created that described improvement in DM prediction due to the interaction. Regions of significant improvement were identified and studied in overall outcome models.

Results

Dose-density interactions were identified that significantly improved prediction for over 50% of bootstrap resamples. Dose and density parameters were not significant when the interaction was omitted. Tumour density variance and high peritumour density were associated with DM for patients with more cold spots (less than 30-Gy EQD2) and non-uniform dose about 3 cm outside of the GTV. Associations identified were independent of the mean GTV dose.

Conclusions

Patients with high tumour variance and peritumour density have increased risk of DM if there is a low and non-uniform dose outside the GTV. The dose regions are independent of tumour dose, suggesting that incidental dose may play an important role in controlling occult disease. Understanding such interactions is key to identifying patients who will benefit from dose-escalation. The methodology presented allowed spatial dose-density interactions to be studied at the exploratory stage for the first time. This could accelerate the clinical implementation of imaging biomarkers by demonstrating the impact of incidental dose for tumours of varying characteristics in routine data.

Spatial descriptions of radiotherapy dose: normal tissue complication models and statistical associations

For decades, dose-volume information for segmented anatomy has provided the essential data for correlating radiotherapy dosimetry with treatment-induced complications. Dose-volume information has formed the basis for modelling those associations via normal tissue complication probability (NTCP) models and for driving treatment planning. Limitations to this approach have been identified. Many studies have emerged demonstrating that the incorporation of information describing the spatial nature of the dose distribution, and potentially its correlation with anatomy, can provide more robust associations with toxicity and seed more general NTCP models. Such approaches are culminating in the application of computationally intensive processes such as machine learning and the application of neural networks. The opportunities these approaches have for individualising treatment, predicting toxicity and expanding the solution space for radiation therapy are substantial and have clearly widespread and disruptive potential. Impediments to reaching that potential include issues associated with data collection, model generalisation and validation. This review examines the role of spatial models of complication and summarises relevant published studies. Sources of data for these studies, appropriate statistical methodology frameworks for processing spatial dose information and extracting relevant features are described. Spatial complication modelling is consolidated as a pathway to guiding future developments towards effective, complication-free radiotherapy treatment.

Biologically consistent dose accumulation using daily patient imaging

BACKGROUND

This work addresses a basic inconsistency in the way dose is accumulated in radiotherapy when predicting the biological effect based on the linear quadratic model (LQM). To overcome this inconsistency, we introduce and evaluate the concept of the total biological dose, bEQD_d.

METHODS

Daily computed tomography imaging of nine patients treated for prostate carcinoma with intensity-modulated radiotherapy was used to compute the delivered deformed dose on the basis of deformable image registration (DIR). We compared conventional dose accumulation (DA) with the newly introduced bEQD_d, a new method of accumulating biological dose that considers each fraction dose and tissue radiobiology. We investigated the impact of the applied fractionation scheme (conventional/hypofractionated), uncertainties induced by the DIR and by the assigned α/β-value.

RESULTS

bEQD_d was systematically higher than the conventionally accumulated dose with difference hot spots of 3.3-4.9 Gy detected in six out of nine patients in regions of high dose gradient in the bladder and rectum. For hypofractionation, differences are up to 8.4 Gy. The difference amplitude was found to be in a similar range to worst-case uncertainties induced by DIR and was higher than that induced by α/β.

CONCLUSION

Using bEQD_d for dose accumulation overcomes a potential systematic inaccuracy in biological effect prediction based on accumulated dose. Highest impact is found for serial-type late responding organs at risk in dose gradient regions and for hypofractionation. Although hot spot differences are in the order of several Gray, in dose-volume parameters there is little difference compared with using conventional or biological DA. However, when local dose information is used, e.g. dose surface maps, difference hot spots can potentially change outcomes of dose-response modelling and adaptive treatment strategies.

Vaginal dose-surface maps in cervical cancer brachytherapy: Methodology and preliminary results on correlation with morbidity

PURPOSE

The purpose of the study was to develop a methodology for vaginal dose-surface maps (DSMs) in patients with cervix cancer and to investigate dose-surface histogram metrics as predictors for vaginal stenosis (St) and mucositis (Muc).

METHODS AND MATERIALS

Thirty-one patients with locally advanced cervix cancer with no vaginal St/Muc (CTCAE-v3) G ≥ 2 at baseline were analyzed. Patients were divided in four morbidity groups: 15 with St/Muc G0/1, 6 with St G ≥ 2, 4 with St/Muc G ≥ 2, and 6 with Muc G ≥ 2. Patients received external beam radiotherapy and 4-fraction intracavitary/interstitial high-dose-rate brachytherapy using tandem and ovoids. DSMs were generated from inner/outer vaginal surfaces. DSMs of external beam radiotherapy and brachytherapy (Gy EQD2₃) were added based on a system of homologous points, to generate cumulative DSMs. Dose-surface histogram/dose-volume histogram parameters, location of high/intermediate-dose regions, rectovaginal reference point, vaginal lateral 5 mm point doses, and vagina/implant dimensions were investigated for St and Muc prediction. Average/difference DSMs and one-way analysis of variance were used to compare between groups.

RESULTS

Best predictors of stenosis were D_15-25cm2 and upper-vagina S_65-120Gy(%). Cutoffs of ∼90 Gy EQD2₃ for D_20cm2 and ∼80% for S_65Gy to top 3 cm inner vaginal surface suitably discriminated for stenosis. Spatial dose location on average/difference DSMs showed significantly higher doses (by > 20 Gy, p < 0.001) over longer parts of the dorsolateral vagina and higher rectovaginal reference point doses for any G ≥ 2 morbidity, over the whole circumference of the upper vagina for G ≥ 2 stenosis. Dose-volume histogram parameters were dependent on vaginal wall thickness. An increase of wall thickness from 2 to 4 mm resulted in an increase of D_2cm3 (D_4cm3) of 16% (32%).

CONCLUSIONS

A novel method was developed to generate vaginal DSMs and spatial-dose metrics. DSMs were found to correlate with vaginal stenosis. The findings of this study are promising and should be further validated on a larger patient cohort, treated with different applicators.

Targeting Cytoprotective Autophagy to Enhance Anticancer Therapies

Autophagy is a highly regulated multi-step process that occurs at the basal level in almost all cells. Although the deregulation of the autophagy process has been described in several pathologies, the role of autophagy in cancer as a cytoprotective mechanism is currently well established and supported by experimental and clinical evidence. Our understanding of the molecular mechanism of the autophagy process has largely contributed to defining how we can harness this process to improve the benefit of cancer therapies. While the role of autophagy in tumor resistance to chemotherapy is extensively documented, emerging data point toward autophagy as a mechanism of cancer resistance to radiotherapy, targeted therapy, and immunotherapy. Therefore, manipulating autophagy has emerged as a promising strategy to overcome tumor resistance to various anti-cancer therapies, and autophagy modulators are currently evaluated in combination therapies in several clinical trials. In this review, we will summarize our current knowledge of the impact of genetically and pharmacologically modulating autophagy genes and proteins, involved in the different steps of the autophagy process, on the therapeutic benefit of various cancer therapies. We will also briefly discuss the challenges and limitations to developing potent and selective autophagy inhibitors that could be used in ongoing clinical trials.

Multivariate normal tissue complication probability models for rectal and bladder morbidity in prostate cancer patients treated with proton therapy

BACKGROUND AND PURPOSE

Normal tissue complication probability (NTCP) models applied for model-based patient selection to proton therapy (PT) have usually been derived using dose/volume histogram (DVH) parameters from photon-based radiotherapy. This study aimed to derive PT-specific multivariate NTCP models that also accounted for the spatial dose distribution (rectum only) as well as non-dose/volume related factors.

MATERIALS AND METHODS

The study included rectum and bladder DVHs, 2D rectal dose maps and relevant patient/treatment characteristics from 1151 prostate cancer cases treated with PT. Prospectively scored Grade 2 late rectal bleeding (CTCAE v3.0, also procedural interventions separately) (n = 156 (15%)) and Grade 3+ GU morbidity (n = 51 (4%)) were entered into a multivariate logistic regression analysis. Model evaluation included assessment of the area under the receiver operating characteristic curve (AUC).

RESULTS

Anticoagulant use was a dominant predictor, chosen in four of the six rectum models and in the bladder model. Age was a dominant predictor in all procedural only rectum models while prostate volume, bladder D5% and V75Gy were predictors in the bladder model. The selection frequency of the dose/volume predictors varied widely, where the percentage of the anterior rectum surface receiving >=75 Gy was the most robust. AUC values ranged from 0.58 to 0.70 across all models, with no clear difference between the DVH- and spatial-based models for the rectum.

CONCLUSION

Anticoagulant use and age were the most prominent predictors in the NTCP models. V75Gy of the rectal wall and the bladder was a predictor in the DVH-based models of the rectum and bladder respectively.

Dose surface maps of the heart can identify regions associated with worse survival for lung cancer patients treated with radiotherapy

•

Creation of a modified cylindrical coordinate system to map cardiac surface dose.

•

Validation of mapping using internal anatomical sub-structures.

•

Identification of surface regions were excess dose is associated with poorer survival.

Background and purpose

For lung cancer patients treated with radiotherapy, radiation dose to the heart has been associated with overall survival, with volumetric dose statistics widely presented. However, critical cardiac structures are present on the hearts surface, where this approach may be sub-optimal. In this work we present a methodology for creating cardiac surface dose maps and identify regions where excess dose is associated with in worse overall survival.

Material and methods

A modified cylindrical coordinate system was implemented to map the cardiac surface dose for lung cancer patients. Validation was performed by mapping the cardiac chambers for 55 patients, fitting a point spread function (PSF) to the blurred edge. To account for this uncertainty, dose maps were blurred by a 2D-Gaussian with width described by the PSF. Permutation testing identified regions where excess dose was associated with worse patient survival. The 99th percentile of the max t-value then defined a cardiac surface region to extract dose, from each patient, to be analysed in a multivariable cox-proportional hazards survival model.

Results

Cardiac surface maps were created for 648 lung cancer patients. Cardiac surface dose maps were blurred with a 2D- Gaussian filter of size σ_φ = 4.3° and σ_y = 1.3units to account for mapping uncertainties. Permutation testing identified significant differences across the surface of the right atria, p < 0.001, at all timepoints. The median dose to the region defined by the 99th percentile of the maximum t-value was 18.5 Gy. Multivariable analysis showed the dose to this region was significantly associated with survival, hazard ratio 1.01 Gy⁻¹, p = 0.03, controlling for confounding variables.

Conclusions

Cardiac surface mapping was successfully implemented and identified a region where excess dose was associated with worse patient survival. This region extended over the right atria, potentially suggesting an interaction with the hearts electrical conduction system.

Increased Dose to Organs in Urinary Tract Associates With Measures of Genitourinary Toxicity in Pooled Voxel-Based Analysis of 3 Randomized Phase III Trials

Purpose: Dose information from organ sub-regions has been shown to be more predictive of genitourinary toxicity than whole organ dose volume histogram information. This study aimed to identify anatomically-localized regions where 3D dose is associated with genitourinary toxicities in healthy tissues throughout the pelvic anatomy. Methods and Materials: Dose distributions for up to 656 patients of the Trans-Tasman Radiation Oncology Group 03.04 RADAR trial were deformably registered onto a single exemplar CT dataset. Voxel- based multiple comparison permutation dose difference testing, Cox regression modeling and LASSO feature selection were used to identify regions where 3D dose-increase was associated with late grade ≥ 2 genitourinary dysuria, incontinence and frequency, and late grade ≥ 1 haematuria. This was externally validated by registering dose distributions from the RT01 (up to n = 388) and CHHiP (up to n = 247) trials onto the same exemplar and repeating the voxel-based tests on each of these data sets. All three datasets were then combined, and the tests repeated. Results: Voxel-based Cox regression and multiple comparison permutation dose difference testing revealed regions where increased dose was correlated with genitourinary toxicity. Increased dose in the vicinity of the membranous and spongy urethra was associated with dysuria for all datasets. Haematuria was similarly correlated with increased dose at the membranous and spongy urethra, for the RADAR, CHHiP, and combined datasets. Some evidence was found for the association between incontinence and increased dose at the internal and external urethral sphincter for RADAR and the internal sphincter alone for the combined dataset. Incontinence was also strongly correlated with dose from posterior oblique beams. Patients with fields extending inferiorly and posteriorly to the CTV, adjacent to the membranous and spongy urethra, were found to experience increased frequency. Conclusions: Anatomically-localized dose-toxicity relationships were determined for late genitourinary symptoms in the urethra and urinary sphincters. Low-intermediate doses to the extraprostatic urethra were associated with risk of late dysuria and haematuria, while dose to the urinary sphincters was associated with incontinence.

Associations between voxel-level accumulated dose and rectal toxicity in prostate radiotherapy

Background and Purpose

Associations between dose and rectal toxicity in prostate radiotherapy are generally poorly understood. Evaluating spatial dose distributions to the rectal wall (RW) may lead to improvements in dose-toxicity modelling by incorporating geometric information, masked by dose-volume histograms. Furthermore, predictive power may be strengthened by incorporating the effects of interfraction motion into delivered dose calculations.Here we interrogate 3D dose distributions for patients with and without toxicity to identify rectal subregions at risk (SRR), and compare the discriminatory ability of planned and delivered dose.

Material and Methods

Daily delivered dose to the rectum was calculated using image guidance scans, and accumulated at the voxel level using biomechanical finite element modelling. SRRs were statistically determined for rectal bleeding, proctitis, faecal incontinence and stool frequency from a training set (n = 139), and tested on a validation set (n = 47).

Results

SRR patterns differed per endpoint. Analysing dose to SRRs improved discriminative ability with respect to the full RW for three of four endpoints. Training set AUC and OR analysis produced stronger toxicity associations from accumulated dose than planned dose. For rectal bleeding in particular, accumulated dose to the SRR (AUC 0.76) improved upon dose-toxicity associations derived from planned dose to the RW (AUC 0.63). However, validation results could not be considered significant.

Conclusions

Voxel-level analysis of dose to the RW revealed SRRs associated with rectal toxicity, suggesting non-homogeneous intra-organ radiosensitivity. Incorporating spatial features of accumulated delivered dose improved dose-toxicity associations. This may be an important tool for adaptive radiotherapy in the future.

Voxel-based analysis in radiation oncology: A methodological cookbook

Recently, 2D or 3D methods for dose distribution analysis have been proposed as evolutions of the Dose Volume Histogram (DVH) approaches. Those methods, collectively referred to as pixel- or voxel-based (VB) methods, evaluate local dose response patterns and go beyond the organ-based philosophy of Normal Tissue Complication Probability (NTCP) modelling. VB methods have been introduced in the context of radiation oncology in the very last years following the virtuous example of neuroimaging experience. In radiation oncology setting, dose mapping is a suitable scheme to compare spatial patterns of local dose distributions between patients who develop toxicity and who do not. In this critical review, we present the methods that include spatial dose distribution information for evaluating different toxicity endpoints after radiation therapy. The review addresses two main topics. First, the critical aspects in dose map building, namely the spatial normalization of the dose distributions from different patients. Then, the issues related to the actual dose map comparison, i.e. the viable options for a robust VB statistical analysis and the potential pitfalls related to the adopted solutions. To elucidate the different theoretical and technical issues, the covered topics are illustrated in relation to practical applications found in the existing literature. We conclude the overview on the VB philosophy in radiation oncology by introducing new phenomenological approaches to NTCP modelling that accounts for inhomogeneous organ radiosensitivity.

A novel framework for spatial normalization of dose distributions in voxel-based analyses of brain irradiation outcomes

PURPOSE

To devise a novel Spatial Normalization framework for Voxel-based analysis (VBA) in brain radiotherapy. VBAs rely on accurate spatial normalization of different patients' planning CTs on a common coordinate system (CCS). The cerebral anatomy, well characterized by MRI, shows instead poor contrast in CT, resulting in potential inaccuracies in VBAs based on CT alone.

METHODS

We analyzed 50 meningioma patients treated with proton-therapy, undergoing planning CT and T1-weighted (T1w) MRI. The spatial normalization pipeline based on MR and CT images consisted in: intra-patient registration of CT to T1w, inter-patient registration of T1w to MNI space chosen as CCS, doses propagation to MNI. The registration quality was compared with that obtained by Statistical Parametric Mapping software (SPM), used as benchmark. To evaluate the accuracy of dose normalization, the dose organ overlap (DOO) score was computed on gray matter, white matter and cerebrospinal fluid before and after normalization. In addition, the trends in the DOOs distribution were investigated by means of cluster analysis.

RESULTS

The registration quality was higher for the proposed method compared to SPM (p < 0.001). The DOO scores showed a significant improvement after normalization (p < 0.001). The cluster analysis highlighted 2 clusters, with one of them including the majority of data and exhibiting acceptable DOOs.

CONCLUSIONS

Our study presents a robust tool for spatial normalization, specifically tailored for brain dose VBAs. Furthermore, the cluster analysis provides a formal criterion for patient exclusion in case of non-acceptable normalization results. The implemented framework lays the groundwork for future reliable VBAs in brain irradiation studies.

Patterns in ano-rectal dose maps and the risk of late toxicity after prostate IMRT

Purpose: The aim of this work was to determine how the spatial pattern of dose in the ano-rectal wall is related to late gastro-intestinal toxicity for prostate cancer patients treated with mainly IMRT.Patients and methods: Patients from the DUE-01 multicentre study with patient-reported (prospective) follow-up and available dosimetric data were included. Conventionally fractionated patients received 74-80 Gy and hypofractionated patients received 65-75.2 Gy. A large majority of the patients were treated with intensity-modulated radiotherapy (IMRT). Dose-surface maps (DSMs) for the anal canal and rectum as a single structure, and for the anal canal and the rectum separately, were co-registered rigidly in two dimensions and, for the patients with and without toxicity, respectively, the mean value of the dose in each pixel was calculated. A pixel-wise t-test was used to highlight the anatomical areas where there was a significant difference between the 'mean dose maps' of each group. Univariate models were also fitted to a range of spatial parameters. The endpoints considered were a mean grade ≥1 late fecal incontinence and a maximum grade ≥2 late rectal bleeding.Results: Twenty-six out of 213 patients had fecal incontinence, while 21/225 patients had rectal bleeding. Incontinence was associated with a higher dose in the caudal region of the anal canal; the most relevant spatial parameter was the lateral extent of the low and medium isodoses (5-49 Gy in EQD2). Bleeding was associated with high isodoses reaching the posterior rectal wall. The spatial dose parameters with the highest AUC value (.69) were the lateral extent of the 60-70 Gy isodoses.Conclusions: To avoid fecal incontinence it is important to limit the portion of the anal canal irradiated. Our analysis confirms that rectal bleeding is a function of similar spatial dose parameters for patients treated with IMRT, compared to previous studies on patients treated with three-dimensional conformal radiotherapy.

Normal tissue complication probability (NTCP) models for modern radiation therapy

Mathematical models of normal tissue complication probability (NTCP) able to robustly predict radiation-induced morbidities (RIM) play an essential role in the identification of a personalized optimal plan, and represent the key to maximizing the benefits of technological advances in radiation therapy (RT). Most modern RT techniques pose, however, new challenges in estimating the risk of RIM. The aim of this report is to schematically review NTCP models in the framework of advanced radiation therapy techniques. Issues relevant to hypofractionated stereotactic body RT and ion beam therapy are critically reviewed. Reirradiation scenarios for new or recurrent malignances and NTCP are also illustrated. A new phenomenological approach to predict RIM is suggested.

PACE: A Probabilistic Atlas for Normal Tissue Complication Estimation in Radiation Oncology

In radiation oncology, the need for a modern Normal Tissue Complication Probability (NTCP) philosophy to include voxel-based evidence on organ radio-sensitivity (RS) has been acknowledged. Here a new formalism (Probabilistic Atlas for Complication Estimation, PACE) to predict radiation-induced morbidity (RIM) is presented. The adopted strategy basically consists in keeping the structure of a classical, phenomenological NTCP model, such as the Lyman-Kutcher-Burman (LKB), and replacing the dose distribution with a collection of RIM odds, including also significant non-dosimetric covariates, as input of the model framework. The theory was first demonstrated in silico on synthetic dose maps, classified according to synthetic outcomes. PACE was then applied to a clinical dataset of thoracic cancer patients classified for lung fibrosis. LKB models were trained for comparison. Overall, the obtained learning curves showed that the PACE model outperformed the LKB and predicted synthetic outcomes with an accuracy >0.8. On the real patients, PACE performance, evaluated by both discrimination and calibration, was significantly higher than LKB. This trend was confirmed by cross-validation. Furthermore, the capability to infer the spatial pattern of underlying RS map for the analyzed RIM was successfully demonstrated, thus paving the way to new perspectives of NTCP models as learning tools.

Voxel-Based Analysis for Identification of Urethrovesical Subregions Predicting Urinary Toxicity After Prostate Cancer Radiation Therapy

PURPOSE

To apply a voxel-based analysis to identify urethrovesical symptom-related subregions (SRSs) associated with acute and late urinary toxicity in prostate cancer radiation therapy.

METHODS AND MATERIALS

Two hundred seventy-two patients with prostate cancer treated with intensity-modulated radiation therapy/image-guided radiation therapy were analyzed prospectively. Each patient's computed tomography imaging was spatially normalized to a common coordinate system via nonrigid registration. The obtained deformation fields were used to map the dose of each patient to the common coordinate system. A voxel-based statistical analysis was applied to generate 3-dimensional dose-volume maps for different urinary symptoms, allowing the identification of corresponding SRSs with statistically significant dose differences between patients with or without toxicity. Each SRS was propagated back to each individual's native space, and dose-volume histograms (DVHs) for the SRSs and the whole bladder were computed. Logistic and Cox regression were used to estimate the SRS's prediction capability compared with the whole bladder.

RESULTS

A local dose-effect relationship was found in the bladder and the urethra. SRSs were identified for 5 symptoms: acute incontinence in the urethra, acute retention in the bladder trigone, late retention and dysuria in the posterior part of the bladder, and late hematuria in the superior part of the bladder, with significant dose differences between patients with and without toxicity, ranging from 1.2 to 9.3 Gy. The doses to the SRSs were significantly predictive of toxicity, with maximum areas under the receiver operating characteristic curve of 0.73 for acute incontinence, 0.62 for acute retention, 0.70 for late retention, 0.81 for late dysuria, and 0.67 for late hematuria. The bladder DVH was predictive only for late retention, dysuria, and hematuria (area under the curve, 0.65-0.72).

CONCLUSIONS

The dose delivered to the urethra and the posterior and superior parts of the bladder was predictive of acute incontinence and retention and of late retention, dysuria, and hematuria. The dose to the whole bladder was moderately predictive.

Development of a prediction model for late urinary incontinence, hematuria, pain and voiding frequency among irradiated prostate cancer patients

BACKGROUND AND PURPOSE

Incontinence, hematuria, voiding frequency and pain during voiding are possible side effects of radiotherapy among patients treated for prostate cancer. The objective of this study was to develop multivariable NTCP models for these side effects.

MATERIAL AND METHODS

This prospective cohort study was composed of 243 patients with localized or locally advanced prostate cancer (stage T1-3). Genito-urinary (GU) toxicity was assessed using a standardized follow-up program. The GU toxicity endpoints were scored using the Common Terminology Criteria for Adverse Events version 3.0 (CTCAE 3.0) scoring system. The full bladder and different anatomical subregions within the bladder were delineated. A least absolute shrinkage and selection operator (LASSO) logistic regression analysis was used to analyze dose volume effects on the four individual endpoints.

RESULTS

In the univariable analysis, urinary incontinence was significantly associated with dose distributions in the trigone (V55-V75, mean). Hematuria was significantly associated with the bladder wall dose (V40-V75, mean), bladder dose (V70-V75), cardiovascular disease and anticoagulants use. Pain during urinating was associated with the dose to the trigone (V50-V75, mean) and with trans transurethral resection of the prostate (TURP). In the final multivariable model urinary incontinence was associated with the mean dose of the trigone. Hematuria was associated with bladder wall dose (V75) and cardiovascular disease, while pain during urinating was associated with trigone dose (V75) and TURP. No significant associations were found for increase in voiding frequency.

CONCLUSIONS

Radiation-induced urinary side effects are associated with dose distributions to different organs as risk. Given the dose effect relationships found, decreasing the dose to the trigone and bladder wall may reduce the incidence of incontinence, pain during voiding and hematuria, respectively.

Investigating rectal toxicity associated dosimetric features with deformable accumulated rectal surface dose maps for cervical cancer radiotherapy

BACKGROUND

Better knowledge of the dose-toxicity relationship is essential for safe dose escalation to improve local control in cervical cancer radiotherapy. The conventional dose-toxicity model is based on the dose volume histogram, which is the parameter lacking spatial dose information. To overcome this limit, we explore a comprehensive rectal dose-toxicity model based on both dose volume histogram and dose map features for accurate radiation toxicity prediction.

METHODS

Forty-two cervical cancer patients treated with combined external beam radiotherapy (EBRT) and brachytherapy (BT) were retrospectively studied, including 12 with Grade ≥ 2 rectum toxicity and 30 patients with Grade 0-1 toxicity (non-toxicity patients). The cumulative equivalent 2-Gy rectal surface dose was deformably summed using the deformation vector fields obtained through a recent developed local topology preserved non-rigid point matching algorithm. The cumulative three-dimensional (3D) dose was flattened and mapped to a two-dimensional (2D) plane to obtain the rectum surface dose map (RSDM). The dose volume parameters (DVPs) were calculated from the 3D rectum surface, while the texture features and the dose geometric parameters (DGPs) were extracted from the 2D RSDM. Representative features further computed from DVPs, textures and DGPs by principle component analysis (PCA) and statistical analysis were respectively fed into a support vector machine equipped with a sequential feature selection procedure. The predictive powers of the representative features were compared with the GEC-ESTRO dosimetric parameters D_0.1/1/2cm³.

RESULTS

Satisfactory predictive accuracy of sensitivity 74.75 and 84.75%, specificity 72.67 and 79.87%, and area under the receiver operating characteristic curve (AUC) 0.82 and 0.91 were respectively achieved by the PCA features and statistical significant features, which were superior to the D_0.1/1/2cm³ (AUC 0.71). The relative area in dose levels of 64Gy, 67Gy, 68Gy, 87Gy, 88Gy and 89Gy, perimeters in dose levels of 89Gy, as well as two texture features were ranked as the important factors that were closely correlated with rectal toxicity.

CONCLUSIONS

Our extensive experimental results have demonstrated the feasibility of the proposed scheme. A future large patient cohort study is still needed for model validation.

A case-control study using motion-inclusive spatial dose-volume metrics to account for genito-urinary toxicity following high-precision radiotherapy for prostate cancer

Background and purpose

Material and methods

Results

Conclusions

Image-based Data Mining to Probe Dosimetric Correlates of Radiation-induced Trismus

PURPOSE

To identify imaged regions in which dose is associated with radiation-induced trismus after head and neck cancer radiation therapy (HNRT) using a novel image-based data mining (IBDM) framework.

METHODS AND MATERIALS

A cohort of 86 HNRT patients were analyzed for region identification. Trismus was characterized as a continuous variable by the maximum incisor-to-incisor opening distance (MID) at 6 months after radiation therapy. Patient anatomies and dose distributions were spatially normalized to a common frame of reference using deformable image registration. IBDM was used to identify clusters of voxels associated with MID (P ≤ .05 based on permutation testing). The result was externally tested on a cohort of 35 patients with head and neck cancer. Internally, we also performed a dose-volume histogram-based analysis by comparing the magnitude of the correlation between MID and the mean dose for the IBDM-identified cluster in comparison with 5 delineated masticatory structures.

RESULTS

A single cluster was identified with the IBDM approach (P < .01), partially overlapping with the ipsilateral masseter. The dose-volume histogram-based analysis confirmed that the IBDM cluster had the strongest association with MID, followed by the ipsilateral masseter and the ipsilateral medial pterygoid (Spearman's rank correlation coefficients: R_s = -0.36, -0.35, -0.32; P = .001, .001, .002, respectively). External validation confirmed an association between mean dose to the IBDM cluster and MID (R_s = -0.45; P = .007).

CONCLUSIONS

IBDM bypasses the common assumption that dose patterns within structures are unimportant. Our novel IBDM approach for continuous outcome variables successfully identified a cluster of voxels that are highly associated with trismus, overlapping partially with the ipsilateral masseter. Tests on an external validation cohort showed an even stronger correlation with trismus. These results support use of the region in HNRT treatment planning to potentially reduce trismus.

Inter-patient image registration algorithms to disentangle regional dose bioeffects

Radiation therapy (RT) technological advances call for a comprehensive reconsideration of the definition of dose features leading to radiation induced morbidity (RIM). In this context, the voxel-based approach (VBA) to dose distribution analysis in RT offers a radically new philosophy to evaluate local dose response patterns, as an alternative to dose-volume-histograms for identifying dose sensitive regions of normal tissue. The VBA relies on mapping patient dose distributions into a single reference case anatomy which serves as anchor for local dosimetric evaluations. The inter-patient elastic image registrations (EIRs) of the planning CTs provide the deformation fields necessary for the actual warp of dose distributions. In this study we assessed the impact of EIR on the VBA results in thoracic patients by identifying two state-of-the-art EIR algorithms (Demons and B-Spline). Our analysis demonstrated that both the EIR algorithms may be successfully used to highlight subregions with dose differences associated with RIM that substantially overlap. Furthermore, the inclusion for the first time of covariates within a dosimetric statistical model that faces the multiple comparison problem expands the potential of VBA, thus paving the way to a reliable voxel-based analysis of RIM in datasets with strong correlation of the outcome with non-dosimetric variables.

Multi-atlas-based segmentation of prostatic urethra from planning CT imaging to quantify dose distribution in prostate cancer radiotherapy

BACKGROUND AND PURPOSE

Segmentation of intra-prostatic urethra for dose assessment from planning CT may help explaining urinary toxicity in prostate cancer radiotherapy. This work sought to: i) propose an automatic method for urethra segmentation in CT, ii) compare it with previously proposed surrogate models and iii) quantify the dose received by the urethra in patients treated with IMRT.

MATERIALS AND METHODS

A weighted multi-atlas-based urethra segmentation method was devised from a training data set of 55 CT scans of patients receiving brachytherapy with visible urinary catheters. Leave-one-out cross validation was performed to quantify the error between the urethra segmentation and the catheter ground truth with two scores: the centerlines distance (CLD) and the percentage of centerline within a certain distance from the catheter (PWR). The segmentation method was then applied to a second test data set of 95 prostate cancer patients having received 78Gy IMRT to quantify dose to the urethra.

RESULTS

Mean CLD was 3.25±1.2mm for the whole urethra and 3.7±1.7mm, 2.52±1.5mm, and 3.01±1.7mm for the top, middle, and bottom thirds, respectively. In average, 53% of the segmented centerlines were within a radius<3.5mm from the centerline ground truth and 83% in a radius<5mm. The proposed method outperformed existing surrogate models. In IMRT, urethra DVH was significantly higher than prostate DVH from V74Gy to V79Gy.

CONCLUSION

A multi-atlas-based segmentation method was proposed enabling assessment of the dose within the prostatic urethra.

Voxel-based analysis unveils regional dose differences associated with radiation-induced morbidity in head and neck cancer patients

The risk of radiation-induced toxicity in patients treated for head and neck (HN) cancer with radiation therapy (RT) is traditionally estimated by condensing the 3D dose distribution into a monodimensional cumulative dose-volume histogram which disregards information on dose localization. We hypothesized that a voxel-based approach would identify correlations between radiation-induced morbidity and local dose release, thus providing a new insight into spatial signature of radiation sensitivity in composite regions like the HN district. This methodology was applied to a cohort of HN cancer patients treated with RT at risk of radiation-induced acute dysphagia (RIAD). We implemented an inter-patient elastic image registration framework that proved robust enough to match even the most elusive HN structures and to provide accurate dose warping. A voxel-based statistical analysis was then performed to test regional dosimetric differences between patients with and without RIAD. We identified a significantly higher dose delivered to RIAD patients in two voxel clusters in correspondence of the cricopharyngeus muscle and cervical esophagus. Our study goes beyond the well-established organ-based philosophy exploring the relationship between radiation-induced morbidity and local dose differences in the HN region. This approach is generally applicable to different HN toxicity endpoints and is not specific to RIAD.

Understanding Urinary Toxicity after Radiotherapy for Prostate Cancer: First Steps Forward

One of the most relevant achievements of Professor Gianni Bonadonna was the implementation of the methodology of controlled clinical trials in medical oncology. It is valid for all cancer types, oncological disciplines and clinical endpoints, both survival and toxicity. This narrative review reports on the status of the current knowledge of the radiation-induced urinary syndrome after external-beam radiotherapy for prostate cancer. In recent years, the syndrome has been the object of large-scale prospective observational trials specifically devoted to investigating the association of patient and treatment features with acute/late urinary toxicity. The first results of these trials allow initial attempts at predictive modeling, which can serve as a basis for the optimization of patient selection and treatment planning.

Associations between volume changes and spatial dose metrics for the urinary bladder during local versus pelvic irradiation for prostate cancer

BACKGROUND

Inter-fractional variation in urinary bladder volumes during the course of radiotherapy (RT) for prostate cancer causes deviations between planned and delivered doses. This study compared planned versus daily cone-beam CT (CBCT)-based spatial bladder dose distributions, for prostate cancer patients receiving local prostate treatment (local treatment) versus prostate including pelvic lymph node irradiation (pelvic treatment).

MATERIAL AND METHODS

Twenty-seven patients (N = 15 local treatment; N = 12 pelvic treatment) were treated using daily image-guided RT (1.8 Gy@43-45 fx), adhering to a full bladder/empty rectum protocol. For each patient, 9-10 CBCTs were registered to the planning CT, using the clinically applied translations. The urinary bladder was manually segmented on each CBCT, 3 mm inner shells were generated, and semi and quadrant sectors were created using axial/coronal cuts. Planned and delivered DVH metrics were compared across patients and between the two groups of treatment (t-test, p < .05; Holm-Bonferroni correction). Associations between bladder volume variations and the dose-volume histograms (DVH) of the bladder and its sectors were evaluated (Spearman's rank correlation coefficient, r_s).

RESULTS

Bladder volumes varied considerably during RT (coefficient of variation: 16-58%). The population-averaged planned and delivered DVH metrics were not significantly different at any dose level. Larger treatment bladder volumes resulted in increased absolute volume of the posterior/inferior bladder sector receiving intermediate-high doses, in both groups. The superior bladder sector received less dose with larger bladder volumes for local treatments (r_s ± SD: -0.47 ± 0.32), but larger doses for pelvic treatments (r_s ± SD: 0.74 ± 0.24).

CONCLUSIONS

Substantial bladder volume changes during the treatment course occurred even though patients were treated under a full bladder/daily image-guided protocol. Larger bladder volumes resulted in less bladder wall spared at the posterior-inferior sector, regardless the treatment received. Contrary, larger bladder volumes meant larger delivered doses to the superior bladder sector for pelvic RT but smaller doses for local treatments.

Spatial features of dose-surface maps from deformably-registered plans correlate with late gastrointestinal complications

This study investigates the associations between spatial distribution of dose to the rectal surface and observed gastrointestinal toxicities after deformably registering each phase of a combined external beam radiotherapy (EBRT)/high-dose-rate brachytherapy (HDRBT) prostate cancer treatment. The study contains data for 118 patients where the HDRBT CT was deformably-registered to the EBRT CT. The EBRT and registered HDRBT TG43 dose distributions in a reference 2 Gy/fraction were 3D-summed. Rectum dose-surface maps (DSMs) were obtained by virtually unfolding the rectum surface slice-by-slice. Associations with late peak gastrointestinal toxicities were investigated using voxel-wise DSM analysis as well as parameterised spatial patterns. The latter were obtained by thresholding DSMs from 1-80 Gy (increment  =  1) and extracting inferior-superior extent, left-right extent, area, perimeter, compactness, circularity and ellipse fit parameters. Logistic regressions and Mann-Whitney U-tests were used to correlate features with toxicities. Rectal bleeding, stool frequency, diarrhoea and urgency/tenesmus were associated with greater lateral and/or longitudinal spread of the high doses near the anterior rectal surface. Rectal bleeding and stool frequency were also influenced by greater low-intermediate doses to the most inferior 20% of the rectum and greater low-intermediate-high doses to 40-80% of the rectum length respectively. Greater low-intermediate doses to the superior 20% and inferior 20% of the rectum length were associated with anorectal pain and urgency/tenesmus respectively. Diarrhoea, completeness of evacuation and proctitis were also related to greater low doses to the posterior side of the rectum. Spatial features for the intermediate-high dose regions such as area, perimeter, compactness, circularity, ellipse eccentricity and confinement to ellipse fits were strongly associated with toxicities other than anorectal pain. Consequently, toxicity is related to the shape of isodoses as well as dose coverage. The findings indicate spatial constraints on doses to certain sections of the rectum may be important for reducing toxicities and optimising dose.

Population model of bladder motion and deformation based on dominant eigenmodes and mixed-effects models in prostate cancer radiotherapy

In radiotherapy for prostate cancer irradiation of neighboring organs at risk may lead to undesirable side-effects. Given this setting, the bladder presents the largest inter-fraction shape variations hampering the computation of the actual delivered dose vs. planned dose. This paper proposes a population model, based on longitudinal data, able to estimate the probability of bladder presence during treatment, using only the planning computed tomography (CT) scan as input information. As in previously-proposed principal component analysis (PCA) population-based models, we have used the data to obtain the dominant eigenmodes that describe bladder geometric variations between fractions. However, we have used a longitudinal analysis along each mode in order to properly characterize patient's variance from the total population variance. We have proposed is a mixed-effects (ME) model in order to separate intra- and inter-patient variability, in an effort to control confounding cohort effects. Other than using PCA, bladder shapes are represented by using spherical harmonics (SPHARM) that additionally enables data compression without information lost. Based on training data from repeated CT scans, the ME model was thus implemented following dimensionality reduction by means of SPHARM and PCA. We have evaluated the model in a leave-one-out cross validation framework on the training data but also using independent data. Probability maps (PMs) were thus generated with several draws from the learnt model as predicted regions where the bladder will likely move and deform. These PMs were compared with the actual regions using metrics based on mutual information distance and misestimated voxels. The prediction was also compared with two previous population PCA-based models. The proposed model was able to reduce the uncertainties in the estimation of the probable region of bladder motion and deformation. This model can thus be used for tailoring radiotherapy treatments.

Modeling Urinary Dysfunction After External Beam Radiation Therapy of the Prostate Using Bladder Dose-Surface Maps: Evidence of Spatially Variable Response of the Bladder Surface

PURPOSE

We assessed the association of the spatial distribution of dose to the bladder surface, described using dose-surface maps, with the risk of urinary dysfunction.

METHODS AND MATERIALS

The bladder dose-surface maps of 754 participants from the TROG 03.04-RADAR trial were generated from the volumetric data by virtually cutting the bladder at the sagittal slice, intersecting the bladder center-of-mass through to the bladder posterior and projecting the dose information on a 2-dimensional plane. Pixelwise dose comparisons were performed between patients with and without symptoms (dysuria, hematuria, incontinence, and an International Prostate Symptom Score increase of ≥10 [ΔIPSS10]). The results with and without permutation-based multiple-comparison adjustments are reported. The pixelwise multivariate analysis findings (peak-event model for dysuria, hematuria, and ΔIPSS10; event-count model for incontinence), with adjustments for clinical factors, are also reported.

RESULTS

The associations of the spatially specific dose measures to urinary dysfunction were dependent on the presence of specific symptoms. The doses received by the anteroinferior and, to lesser extent, posterosuperior surface of the bladder had the strongest relationship with the incidence of dysuria, hematuria, and ΔIPSS10, both with and without adjustment for clinical factors. For the doses to the posteroinferior region corresponding to the area of the trigone, the only symptom with significance was incontinence.

CONCLUSIONS

A spatially variable response of the bladder surface to the dose was found for symptoms of urinary dysfunction. Limiting the dose extending anteriorly might help reduce the risk of urinary dysfunction.

Predicting toxicity in radiotherapy for prostate cancer

This comprehensive review addresses most organs at risk involved in planning optimization for prostate cancer. It can be considered an update of a previous educational review that was published in 2009 (Fiorino et al., 2009). The literature was reviewed based on PubMed and MEDLINE database searches (from January 2009 up to September 2015), including papers in press; for each section/subsection, key title words were used and possibly combined with other more general key-words (such as radiotherapy, dose-volume effects, NTCP, DVH, and predictive model). Publications generally dealing with toxicity without any association with dose-volume effects or correlations with clinical risk factors were disregarded, being outside the aim of the review. A focus was on external beam radiotherapy, including post-prostatectomy, with conventional fractionation or moderate hypofractionation (<4Gy/fraction); extreme hypofractionation is the topic of another paper in this special issue. Gastrointestinal and urinary toxicity are the most investigated endpoints, with quantitative data published in the last 5years suggesting both a dose-response relationship and the existence of a number of clinical/patient related risk factors acting as dose-response modifiers. Some results on erectile dysfunction, bowel toxicity and hematological toxicity are also presented.

Bladder dose-surface maps and urinary toxicity: Robustness with respect to motion in assessing local dose effects

The purpose of this study was to quantify the impact of inter-fraction modifications of bladder during RT of prostate cancer on bladder dose surface maps (DSM). Eighteen patients treated with daily image-guided Tomotherapy and moderate hypofractionation (70-72.8Gy at 2.5-2.6Gy/fr in 28 fractions and full bladder) were considered. Bladder contours were delineated on co-registered daily Megavoltage CT (MVCT) by a single observer and copied on the planning CT to generate dose-volume/surface histograms (DVH/DSH) and bladder DSMs. Discrepancies between planned and daily absorbed doses were analyzed through the average of individual systematic errors, the population systematic errors and the population random errors for the DVH/DSHs and DSMs. In total, 477 DVH/DSH and 472 DSM were available. DSH and DVH showed small population systematic errors of absolute surfaces (<3.4cm(2)) and volumes (<8.4cm(3)) at the highest doses. The dose to the posterior bladder base assessed on DSMs showed a mean systematic error below 1Gy, with population systematic and random errors within 4 and 3Gy, respectively. The region surrounding this area shows higher mean systematic errors (1-3Gy), population systematic (8-11Gy) and random (5-7Gy) errors. In conclusion, DVH/DSH and DSMs are quite stable with respect to inter-fraction variations in the high-dose region, within about 2cm from bladder base. Larger systematic variations occur in the anterior portion and cranially 2.5-3.5cm from the base. Results suggest that dose predictors related to the high dose area (including the trigone dose) are likely to be sufficiently reliable with respect to the expected variations due to variable bladder filling.