Dose-surface maps identifying local dose-effects for acute gastrointestinal toxicity after radiotherapy for prostate cancer

The use of dose surface maps as a tool to investigate spatial dose delivery accuracy for the rectum during prostate radiotherapy

PURPOSE

This study aims to address the lack of spatial dose comparisons of planned and delivered rectal doses during prostate radiotherapy by using dose-surface maps (DSMs) to analyze dose delivery accuracy and comparing these results to those derived using DVHs.

METHODS

Two independent cohorts were used in this study: twenty patients treated with 36.25 Gy in five fractions (SBRT) and 20 treated with 60 Gy in 20 fractions (IMRT). Daily delivered rectum doses for each patient were retrospectively calculated using daily CBCT images. For each cohort, planned and average-delivered DVHs were generated and compared, as were planned and accumulated DSMs. Permutation testing was used to identify DVH metrics and DSM regions where significant dose differences occurred. Changes in rectal volume and position between planning and delivery were also evaluated to determine possible correlation to dosimetric changes.

RESULTS

For both cohorts, DVHs and DSMs reported conflicting findings on how planned and delivered rectum doses differed from each other. DVH analysis determined average-delivered DVHs were on average 7.1% ± 7.6% (p ≤ 0.002) and 5.0 ± 7.4% (p ≤ 0.021) higher than planned for the IMRT and SBRT cohorts, respectively. Meanwhile, DSM analysis found average delivered posterior rectal wall dose was 3.8 ± 0.6 Gy (p = 0.014) lower than planned in the IMRT cohort and no significant dose differences in the SBRT cohort. Observed dose differences were moderately correlated with anterior-posterior rectal wall motion, as well as PTV superior-inferior motion in the IMRT cohort. Evidence of both these relationships were discernable in DSMs.

CONCLUSION

DSMs enabled spatial investigations of planned and delivered doses can uncover associations with interfraction motion that are otherwise masked in DVHs. Investigations of dose delivery accuracy in radiotherapy may benefit from using DSMs over DVHs for certain organs such as the rectum.

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.

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.

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.

Normal tissue complication probability models for prospectively scored late rectal and urinary morbidity after proton therapy of prostate cancer

Background and purpose

Photons and protons have fundamentally different properties, i.e. protons have a reduced dose bath but a higher relative biological effectiveness. Photon-based normal tissue complication probability (NTCP) models may therefore not immediately be applicable to proton therapy (PT). The aim was to derive parameters of the Lyman-Kutcher-Burman (LKB) NTCP model using prospectively recorded late morbidity data from PT, focusing on rectal morbidity and prostate cancer.

Materials and methods

Prospectively collected data were available for 1151 prostate cancer patients treated with passive scattering PT and prescribed target doses of 78–82 Gy (RBE = 1.1) in 2 Gy fractions. Morbidity data (CTCAE v3.0) consisted of two alternative late grade 2 rectal bleeding endpoints: Medical Grade2A (GR2A) and procedural Grade2B (GR2B), as well as late grade 3 + urinary morbidity. GR2A + 2B were observed in 156/1047 patients (15%), GR2B in 45/1047 patients (4%), and urinary grade 3 + in 51/1151 patients (4%). LKB NTCP model parameters (D50, m, and n) were derived by maximum likelihood estimation.

Results

For the rectum/rectal wall the volume parameter n was low (0.07–0.14) for both GR2A + 2B and GR2B, as was the m parameter (range: 0.16–0.20). For the bladder/bladder wall both parameters were high (n-range: 0.20–0.36; m-range: 0.32–0.36). D50 parameters were higher for GR2B of the rectum/rectal wall (95.9–98.0 Gy) and bladder/bladder wall (118.1–119.9 Gy), but lower for GR2A2B (71.7–73.6 Gy).

Conclusion

PT specific LKB NTCP model parameters were derived from a population of more than 1000 patients. The D50 parameter differed for all structures and endpoints and deviated from typical photon-based LKB model values.

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.

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.

Patient-reported acute GI symptoms in locally advanced cervical cancer patients correlate with rectal dose

BACKGROUND AND PURPOSE

To investigate relationships between patient-reported acute gastro-intestinal symptoms in a locally advanced cervical cancer (LACC) prospective cohort and clinical and dosimetric parameters, while also taking spatial dose into account.

MATERIAL AND METHODS

A total of 103 patients was included, receiving radiotherapy based on a plan-library-based plan-of-the-day protocol, combined either with concurrent chemotherapy or with neo-adjuvant chemotherapy and concomitant hyperthermia. Toxicity endpoints were extracted from questionnaires sent out weekly during treatment and regularly in the acute phase after treatment. Endpoints were defined for symptoms concerning obstipation, diarrhea, fecal leakage, bowel cramps and rectal bleeding. Dose surface maps were constructed for the rectum. Clinical parameters and dosimetric parameters of the bowel bag and rectum were collected for all patients.

RESULTS

The use of concomitant chemotherapy and an increase in Planning Target Volume (PTV) resulted in a significant increase in reported diarrhea. The dose-volume parameters V_5Gy-V_25Gy of the rectum were found to be significant, unlike dose-volume parameters of the bowel bag. Additionally, a significantly higher dose to the inferior part of the rectum was found for patients reporting diarrhea. No significance was reached for fecal leakage and bowel cramps.

CONCLUSION

The significance of results for patients reporting diarrhea symptoms found for PTV volume indicates a potential benefit for a plan-of-the-day protocol. Additionally, the results suggest that a reduction of inferior rectum dose could decrease patient-reported diarrhea symptoms, while the administration of concomitant chemotherapy appears to lead to radiosensitizing effects that increase these symptoms.

Local Dose Effects for Late Gastrointestinal Toxicity After Hypofractionated and Conventionally Fractionated Modern Radiotherapy for Prostate Cancer in the HYPRO Trial

Purpose: Late gastrointestinal (GI) toxicity after radiotherapy for prostate cancer may have significant impact on the cancer survivor's quality of life. To date, little is known about local dose-effects after modern radiotherapy including hypofractionation. In the current study we related the local spatial distribution of radiation dose in the rectum to late patient-reported gastrointestinal (GI) toxicities for conventionally fractionated (CF) and hypofractionated (HF) modern radiotherapy in the randomized HYPRO trial. Material and Methods: Patients treated to 78 Gy in 2 Gy fractions (n = 298) or 64.6 Gy in 3.4 Gy fractions (n = 295) with available late toxicity questionnaires (n ≥ 2 within 1-5 years post-treatment) and available 3D planning data were eligible for this study. The majority received intensity modulated radiotherapy (IMRT). We calculated two types of dose surface maps: (1) the total delineated rectum with its central axis scaled to unity, and (2) the delineated rectum with a length of 7 cm along its central axis aligned on the prostate's half-height point (prostate-half). For each patient-reported GI symptom, dose difference maps were constructed by subtracting average co-registered EQD2 (equivalent dose in 2 Gy) dose maps of patients with and without the symptom of interest, separately for HF and CF. P-values were derived from permutation tests. We evaluated patient-reported moderate to severe GI symptoms. Results: Observed incidences of rectal bleeding and increased stool frequency were significantly higher in the HF group. For rectal bleeding (p = 0.016), mucus discharge (p = 0.015), and fecal incontinence (p = 0.001), significant local dose-effects were observed in HF patients but not in CF patients. For rectal pain, similar local dose-effects (p < 0.05) were observed in both groups. No significant local dose-effects were observed for increased stool frequency. Total rectum mapping vs. prostate-half mapping showed similar results. Conclusion: We demonstrated significant local dose-effect relationships for patient-reported late GI toxicity in patients treated with modern RT. HF patients were at higher risk for increased stool frequency and rectal bleeding, and showed the most pronounced local dose-effects in intermediate-high dose regions. These findings suggest that improvement of current treatment optimization protocols could lead to clinical benefit, in particular for HF treatment.

Dosimetric impact of organ at risk daily variation during prostate stereotactic ablative radiotherapy

OBJECTIVE

Prostate stereotactic ablative radiotherapy (SABR) delivers large doses using a fast dose rate. This amplifies the effect geometric uncertainties have on normal tissue dose. The aim of this study was to determine whether the treatment dose-volume histogram (DVH) agrees with the planned dose to organs at risk (OAR).

METHODS

41 low-intermediate risk prostate cancer patients were treated with SABR using a linac based technique. Dose prescribed was 35 Gy in five fractions delivered on alternate days, planned using volumetric modulated arc therapy (VMAT) with 10X flattening filter free (FFF). On treatment, prostate was matched to fiducial markers on cone beam CT (CBCT). OAR were retrospectively delineated on 205 pre-treatment CBCT images. Daily CBCT contours were overlaid on the planning CT for dosimetric analysis. Verification plan used to evaluate the daily DVH for each structure. The daily doses received by OAR were recorded using the D%.

RESULTS

The median rectum and bladder volumes at planning were 67.1 cm3 (interquartile range 56.4-78.2) and 164.4 cm3 (interquartile range 120.3-213.4) respectively. There was no statistically significant difference in median rectal volume at each of the five treatment scans compared to the planning scan (p = 0.99). This was also the case for median bladder volume (p = 0.79). The median dose received by rectum and bladder at each fraction was higher than planned, at the majority of dose levels. For rectum the increase ranged from 0.78-1.64Gy and for bladder 0.14-1.07Gy. The percentage of patients failing for rectum D35% < 18 Gy (p = 0.016), D10% < 28 Gy (p = 0.004), D5% < 32 Gy (p = 0.0001), D1% < 35 Gy (p = 0.0001) and bladder D1% < 35 Gy (p = 0.001) at treatment were all statistically significant.

CONCLUSION

In this cohort of prostate SABR patients, we estimate the OAR treatment DVH was higher than planned. This was due to rectal and bladder organ variation.

ADVANCES IN KNOWLEDGE

OAR variation in prostate SABR using a FFF technique, may cause the treatment DVH to be higher than planned.

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.

Evaluating the predictive value of quantec rectum tolerance dose suggestions on acute rectal toxicity in prostate carcinoma patients treated with IMRT

Aim

To investigate the predictive value of convenience of rectum dosimetry with Quantitative Analysis of Normal Tissue Effects in the Clinic (QUANTEC) dose limits, maximum rectum dose (Dmax), total rectal volume (TVrectum), rectal volume included in PTV (VrectumPTV) on Grade 2-3 acute rectal toxicity for utilization in clinical practice.

Background

Numerous previous data have reported frequent acute proctitis after external-beam RT of prostate cancer. Predicting toxicity limited with dose information is inadequate in clinical practice due to comorbidities and medications used.

Materials and Method

Sixty-four non-metastatic prostate cancer patients treated with IMRT were enrolled. Patients were treated to a total dose of 70-76 Gy. Rectal dose volume histograms (DVH) of all patients were evaluated retrospectively, and a QUANTEC Score between 0 and 5 was calculated for each patient. The correlation between the rectal DVH data, QUANTEC score, TVrectum, VrectumPTV, rectum Dmax and Grade 2-3 rectal toxicity was investigated.

Results

In the whole group grade 1, 2 and 3 acute rectal toxicities were 25%, 18.8% and 3.1%, respectively. In the DVH data, rectum doses of all patients were under RTOG dose limits. Statistically significant correlation was found between grade 2-3 rectal toxicity and TVrectum (p = 0,043); however. It was not correlated with QUANTEC score, VrectumPTV and Dmax.

Conclusion

Our results were not able to show any significant correlation between increasing convenience with QUANTEC limits and lower rectal toxicity. Conclusively, new dosimetric definitions are warranted to predict acute rectal toxicity more accurately in prostate cancer patients during IMRT treatment.

Assessing localized dosimetric effects due to unplanned gas cavities during pelvic MR-guided radiotherapy using Monte Carlo simulations

PURPOSE

It has been proposed that beam modulation and opposing beam configurations can cancel effects of the Electron Return Effect (ERE) during MR-guided radiotherapy (MRgRT). However, this may not always be the case for unplanned gas cavities outside of the target in the pelvic region. We evaluate dosimetric effects, including effects in the rectal wall, due to unplanned spherical air cavities during MRgRT.

METHODS

Nine virtual cuboid water phantoms containing spherical air cavities (0.5-7.5 cm diameter) and a reference phantom without air were created. Monte Carlo dose calculations of 7 MV photons under the influence of a 1.5 T transverse magnetic field were produced using Monaco 5.19.02 Treatment Planning System (TPS) (Elekta AB, Stockholm, Sweden). Cavities in the path of a single and multiple beam plans were considered. Dose distributions of phantoms with and without air cavities were compared (ΔD% ) using a spherical coordinate system originating in the center of the cavity. Effects in the rectal wall were quantified by comparing dose volume histogram (DVH) parameters for solid and gaseous filling from simulated rectal wall structures.

RESULTS

Max(ΔD% ) of ~70% and 20% were observed around large cavities in the path of a single and multiple beam plans, respectively. Approximately 45 cm3 of phantom surrounding the largest cavity in a single beam received dose changes of >10%. Dmean in the rectal wall was unchanged when comparing gaseous and solid filling in the path of a single beam; however, D1cc and Dmax increased by up to ~45% and ~63%, respectively.

CONCLUSIONS

Unplanned gas cavities in the path of a single beam during pelvic MRgRT with a 1.5 T transverse magnetic field cause dose changes which may impact toxicity in the rectal wall, depending on local dose and fractionation. Effects are reduced but not eliminated with a five-beam plan.

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.

Dose-volume analysis of predictors for acute anal toxicity after radiotherapy in prostate cancer patients

Background

This study aimed to evaluate the clinical and dosimetric factors predictive of acute anal toxicity (AAT) after radiotherapy in prostate cancer (PCa) patients with or without hemorrhoids.

Methods

We analyzed data from 347 PCa patients (248 cases treated from July 2013 to November 2017 for training cohort and 99 cases treated in 2018 for validation cohort) treated with pelvic radiotherapy at a single institution. Anal canal dose–volume histogram was used to determine the prescribed dose. Univariate and multivariate analyses were used to evaluate the risk of AAT as a function of clinical and dosimetric factors.

Results

Totally, 39.5% (98/248) and 31.3% (31/99) of the PCa patients developed AAT in training and validation cohorts, respectively. The incidence of AAT was much higher in patients with hemorrhoids than in those without hemorrhoids in both training and validation cohorts. Hemorrhoids and volume received more than 20 Gy (V20) were valuated as independent factors for predicting AAT in training cohort. Similar results were also observed in our validation cohort. The combination of hemorrhoids and high anal canal V20 (> 74.93% as determined by ROC curves) showed the highest specificity and positive predictive values for predicting AAT in both training and validation cohorts.

Conclusions

AAT occurs commonly in PCa patients with hemorrhoids during and after pelvic radiotherapy. Hemorrhoids and anal canal V20 are independent predictors of AAT. These factors should be carefully considered during treatment planning to minimize the incidence of AAT.

Sildenafil Protects Endothelial Cells From Radiation-Induced Oxidative Stress

Introduction:

The etiology of radiation-induced erectile dysfunction (ED) is complex and multifactorial, and it appears to be mainly atherogenic.

Aim:

To focus on vascular aspects of radiation-induced ED and to elucidate whether the protective effects of sildenafil are mediated by attenuation of oxidative stress and apoptosis in the endothelial cells.

Methods:

Bovine aortic endothelial cells (BAECs), with or without pretreatment of sildenafil (5 μM at 5 minutes before radiation), were used to test endothelial dysfunction in response to external beam radiation at 10e15 Gy. Generation of reactive oxygen species (ROS) was studied. Extracellular hydrogen peroxide (H₂O₂) was measured using the Amplex Red assay and intracellular H₂O₂ using a fluorescent sensor. In addition, ROS superoxide (O₂•-) was measured using a O₂•- chemiluminescence enhancer. Both H₂O₂ and O2•- are known to reduce the bioavailability of nitric oxide, which is the most significant chemical mediator of penile erection. Generation of cellular peroxynitrite (ONOO⁻) was measured using a chemiluminescence assay with the PNCL probe. Subsequently, we measured the activation of acid sphingomyelinase (ASMase) enzyme by radioenzymatic assay using [¹⁴C-methylcholine] sphingomyelin as substrate, and the generation of the proapoptotic C₁₆-ceramide was assessed using the diacylglycerol kinase assay. Endothelial cells apoptosis was measured as a readout of these cells’ dysfunction.

Main Outcome Measures:

Single high-dose radiation therapy induced NADPH oxidases (NOXs) activation and ROS generation via the proapoptotic ASMase/ceramide pathway. The radio-protective effect of sildenafil on BAECs was due to inhibition of this pathway.

Results:

Here, we demonstrate for the first time that radiation activated NOXs and induced generation of ROS in BAECs. In addition, we showed that sildenafil significantly reduced radiation-induced O₂•- and as a result there was reduction in the generation of peroxynitrite in these cells. Subsequently, sildena fil protected the endothelial cells from radiation therapy-induced apoptosis.

Strengths and Limitations:

This is the first study demonstrating that single high-dose radiation therapy induced NOXs activation, resulting in the generation of O₂•- and peroxynitrite in endothelial cells. Sildenafil reduced ROS generation by inhibiting the ASMase/ceramide pathway. These studies should be followed in an animal model of ED.

Conclusions:

This study demonstrated that sildenafil protects BAECs from radiation-induced oxidative stress by reducing NOX-induced ROS generation, thus resulting in decreased endothelial dysfunction. Therefore, it provides a potential mechanism to better understand the atherogenic etiology of postradiation ED.

[Prediction of rectal toxicity of radiotherapy for prostate cancer based on multi-modality feature and multi-classifiers]

OBJECTIVE

To evaluate rectal toxicity of radiotherapy for prostate cancer using a novel predictive model based on multi-modality and multi-classifier fusion.

METHODS

We retrospectively collected the clinical data from 44 prostate cancer patients receiving external beam radiation (EBRT), including the treatment data, clinical parameters, planning CT data and the treatment plans. The clinical parameter features and dosimetric features were extracted as two different modality features, and a subset of features was selected to train the 5 base classifiers (SVM, Decision Tree, K-nearest-neighbor, Random forests and XGBoost). To establish the multi-modality and multi-classifier fusion model, a multi-criteria decision-making based weight assignment algorithm was used to assign weights for each base classifier under the same modality. A repeat 5-fold cross-validation and the 4 indexes including the area under ROC curve (AUC), accuracy, sensitivity and specificity were used to evaluate the proposed model. In addition, the proposed model was compared quantitatively with different feature selection methods, different weight allocation algorithms, the model based on single mode single classifier, and two integrated models using other fusion methods.

RESULTS

Repeated (5 times) 5-fold cross validation of the proposed model showed an accuracy of 0.78 for distinguishing toxicity from non-toxicity with an AUC of 0.83, a specificity of 0.79 and a sensitivity of 0.76.

CONCLUSIONS

Compared with the models based on a single mode or a single classifier and other fusion models, the proposed model can more accurately predict rectal toxicity of radiotherapy for prostate cancer.

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.

Ano-rectal wall dose-surface maps localize the dosimetric benefit of hydrogel rectum spacers in prostate cancer radiotherapy

BACKGROUND AND PURPOSE

To evaluate spatial differences in dose distributions of the ano-rectal wall (ARW) using dose-surface maps (DSMs) between prostate cancer patients receiving intensity-modulated radiation therapy with and without implantable rectum spacer (IMRT+IRS; IMRT-IRS, respectively), and to correlate this with late gastro-intestinal (GI) toxicities using validated spatial and non-spatial normal-tissue complication probability (NTCP) models.

MATERIALS AND METHODS

For 26 patients DSMs of the ARW were generated. From the DSMs various shape-based dose measures were calculated at different dose levels: lateral extent, longitudinal extent, and eccentricity. The contiguity of the ARW dose distribution was assessed by the contiguous-DSH (cDSH). Predicted complication rates between IMRT+IRS and IMRT-IRS plans were assessed using a spatial NTCP model and compared against a non-spatial NTCP model.

RESULTS

Dose surface maps are generated for prostate radiotherapy using an IRS. Lateral extent, longitudinal extent and cDSH were significantly lower in IMRT+IRS than for IMRT-IRS at high-dose levels. Largest significant differences were observed for cDSH at dose levels >50 Gy, followed by lateral extent at doses >57 Gy, and longitudinal extent in anterior and superior-inferior directions. Significant decreases (p = 0.01) in median rectal and anal NTCPs (respectively, Gr 2 late rectal bleeding and subjective sphincter control) were predicted when using an IRS.

CONCLUSIONS

Local-dose effects are predicted to be significantly reduced by an IRS. The spatial NTCP model predicts a significant decrease in Gr 2 late rectal bleeding and subjective sphincter control. Dose constraints can be improved for current clinical treatment planning.

An Investigation of Dosimetric Correlates of Acute Toxicity in Prostate Stereotactic Body Radiotherapy: Dose to Urinary Trigone is Associated with Acute Urinary Toxicity

AIMS

There are limited data on dosimetric correlates of toxicity in stereotactic body radiotherapy (SBRT) for prostate cancer. We aimed to identify potential relationships between dose and toxicity using conventional dose-volume histograms (DVHs) and dose-surface maps (DSMs).

MATERIALS AND METHODS

Urinary bladder trigone and rectum DSMs were produced for a single-institution service evaluation cohort of 50 patients receiving SBRT for localised prostate cancer, together with conventional DVHs for bladder and rectum. Patients had been prospectively recruited to this cohort and treated according to a pre-defined protocol to a dose of 36.25 Gy in five fractions. Radiation Therapy Oncology Group (RTOG) and International Prostate Symptom Score (IPSS) toxicity data were recorded prospectively. Logistic regression was used to identify dosimetric predictors of acute IPSS+10 (rise of 10 points or more above baseline) and grade 2+ RTOG toxicity.

RESULTS

On univariate analysis, trigone area receiving 40 Gy and trigone Dmax were associated with IPSS+10 (odds ratio 1.06 [1.02-1.11], P = 0.007 and odds ratio 1.54 [1.06-2.25], P = 0.024, respectively). These two variables were highly correlated. In a multivariate model, including all baseline variables, trigone Dmax remained associated with IPSS+10 (odds ratio 1.91 [1.13-3.22], P = 0.016). These findings were not significant with Holm-Bonferroni correction for multiple testing (corrected P value threshold 0.006). No associations were seen between rectal toxicity and DVH or DSM parameters.

CONCLUSIONS

Our study suggests a potential relationship between high doses to the urinary bladder trigone and patient-reported urinary toxicity in prostate SBRT, and is consistent with previous studies in conventionally fractionated radiotherapy, justifying further evaluation in larger cohorts.

Reporting Erectile Function Outcomes After Radiation Therapy for Prostate Cancer: Challenges in Data Interpretation

BACKGROUND

Choice of prostate cancer treatment is frequently influenced by the expected chance of treatment-induced side effects such as erectile dysfunction (ED). However, great discrepancy in cited ED rates exists in the contemporary radiation therapy literature.

AIM

To analyze the reported ED rates and cause of discrepancies and explore the strengths and limitations in the literature on radiation-induced ED.

METHODS

We performed a PubMed literature search and reviewed the literature on ED rates associated with external-beam radiotherapy and brachytherapy from the past 10 years. Eighteen studies were eligible for inclusion and subsequently reviewed.

OUTCOMES

Variables required for interpretation of erectile function outcomes, including patient demographics, treatment characteristics, and sexual function outcomes.

RESULTS

A large variety in the reported incidence of ED was found among studies. In part, these differences resulted from large variations in (i) study populations, (ii) patient characteristics, (iii) treatment characteristics, (iv) prescription of androgen deprivation therapy, (v) means of data acquisition, (vi) definitions of ED, (vii) temporal considerations, and (viii) erectile aid use. Relevant data required for adequate appraisal of sexual function outcomes were not always reported.

CLINICAL IMPLICATIONS

Based on the present findings, we present general recommendations for reporting of erectile function outcomes after radiotherapy for prostate cancer. These should improve future reports.

STRENGTHS AND LIMITATIONS

This is the first report that presents general requirements on reporting erectile function outcomes in the setting of radiotherapy for prostate cancer. We did not conduct a formal meta-analysis because we focused on concepts of research design; this might be considered a limitation.

CONCLUSION

In this review, we have highlighted the strengths and deficiencies of the current literature on ED after external-beam radiotherapy and brachytherapy for prostate cancer. We have made general recommendations to achieve some degree of standardization among reports and improve clinical interpretability. Wortel RC, Incrocci L, Muhall JP. Reporting Erectile Function Outcomes After Radiation Therapy for Prostate Cancer: Challenges in Data Interpretation. J Sex Med 2017;14:1260-1269.

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.

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-institutional analysis demonstrates the importance of lower than previously anticipated dose regions to prevent late rectal bleeding following prostate radiotherapy

PURPOSE

To investigate whether inter-institutional cohort analysis uncovers more reliable dose-response relationships exemplified for late rectal bleeding (LRB) following prostate radiotherapy.

MATERIAL AND METHODS

Data from five institutions were used. Rectal dose-volume histograms (DVHs) for 989 patients treated with 3DCRT or IMRT to 70-86.4 Gy@1.8-2.0 Gy/fraction were obtained, and corrected for fractionation effects (α/β = 3 Gy). Cohorts with best-fit Lyman-Kutcher-Burman volume-effect parameter a were pooled after calibration adjustments of the available LRB definitions. In the pooled cohort, dose-response modeling (incorporating rectal dose and geometry, and patient characteristics) was conducted on a training cohort (70%) followed by final testing on the remaining 30%. Multivariate logistic regression was performed to build models with bootstrap stability.

RESULTS

Two cohorts with low bleeding rates (2%) were judged to be inconsistent with the remaining data, and were excluded. In the remaining pooled cohorts (n = 690; LRB rate = 12%), an optimal model was generated for 3DCRT using the minimum rectal dose and the absolute rectal volume receiving less than 55 Gy (AUC = 0.67; p = 0.0002; Hosmer-Lemeshow p-value, pHL = 0.59). The model performed nearly as well in the hold-out testing data (AUC = 0.71; p < 0.0001; pHL = 0.63), indicating a logistically shaped dose-response.

CONCLUSION

We have demonstrated the importance of integrating datasets from multiple institutions, thereby reducing the impact of intra-institutional dose-volume parameters explicitly correlated with prescription dose levels. This uncovered an unexpected emphasis on sparing of the low to intermediate rectal dose range in the etiology of late rectal bleeding following prostate radiotherapy.

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.

Predictors of urinary and rectal toxicity after external conformed radiation therapy in prostate cancer: Correlation between clinical, tumour and dosimetric parameters and radical and postoperative radiation therapy

OBJECTIVE

To determine rectal and urinary toxicity after external beam radiation therapy (EBRT), assessing the results of patients who undergo radical or postoperative therapy for prostate cancer (pancreatic cancer) and their correlation with potential risk factors.

METHOD

A total of 333 patients were treated with EBRT. Of these, 285 underwent radical therapy and 48 underwent postoperative therapy (39 cases of rescue and 9 of adjuvant therapy). We collected clinical, tumour and dosimetric variable to correlate with toxicity parameters. We developed decision trees based on the degree of statistical significance.

RESULTS

The rate of severe acute toxicity, both urinary and rectal, was 5.4% and 1.5%, respectively. The rate of chronic toxicity was 4.5% and 2.7%, respectively. Twenty-seven patients presented haematuria, and 9 presented haemorrhagic rectitis. Twenty-five patients (7.5%) presented permanent limiting sequela. The patients with lower urinary tract symptoms prior to the radiation therapy presented poorer tolerance, with greater acute bladder toxicity (P=0.041). In terms of acute rectal toxicity, 63% of the patients with mean rectal doses >45Gy and anticoagulant/antiplatelet therapy developed mild toxicity compared with 37% of the patients with mean rectal doses <45 Gy and without anticoagulant therapy. We were unable to establish predictors of chronic toxicity in the multivariate analysis. The long-term sequelae were greater in the patients who underwent urological operations prior to the radiation therapy and who were undergoing anticoagulant therapy.

CONCLUSIONS

The tolerance to EBRT was good, and severe toxicity was uncommon. Baseline urinary symptoms constitute the predictor that most influenced the acute urinary toxicity. Rectal toxicity is related to the mean rectal dose and with anticoagulant/antiplatelet therapy. There were no significant differences in severe toxicity between radical versus postoperative radiation therapy.

Spatial rectal dose/volume metrics predict patient-reported gastro-intestinal symptoms after radiotherapy for prostate cancer

BACKGROUND

Gastro-intestinal (GI) toxicity after radiotherapy (RT) for prostate cancer reduces patient's quality of life. In this study, we explored associations between spatial rectal dose/volume metrics and patient-reported GI symptoms after RT for localized prostate cancer, and compared these with those of dose-surface/volume histogram (DSH/DVH) metrics.

MATERIAL AND METHODS

Dose distributions and six GI symptoms (defecation urgency/emptying difficulties/fecal leakage, ≥Grade 2, median follow-up: 3.6 y) were extracted for 200 patients treated with image-guided RT in 2005-2007. Three hundred and nine metrics assessed from 2D rectal dose maps or DSHs/DVHs were subject to 50-times iterated five-fold cross-validated univariate and multivariate logistic regression analysis (UVA, MVA). Performance of the most frequently selected MVA models was evaluated by the area under the receiving-operating characteristics curve (AUC).

RESULTS

The AUC increased for dose-map compared to DSH/DVH-based models (mean SD: 0.64 ± 0.03 vs. 0.61 ± 0.01), and significant relations were found for six versus four symptoms. Defecation urgency and faecal leakage were explained by high doses at the central/upper and central areas, respectively; while emptying difficulties were explained by longitudinal extensions of intermediate doses.

CONCLUSIONS

Predictability of patient-reported GI toxicity increased using spatial metrics compared to DSH/DVH metrics. Novel associations were particularly identified for emptying difficulties using both approaches in which intermediate doses were emphasized.

Deep convolutional neural network with transfer learning for rectum toxicity prediction in cervical cancer radiotherapy: a feasibility study

Better understanding of the dose-toxicity relationship is critical for safe dose escalation to improve local control in late-stage cervical cancer radiotherapy. In this study, we introduced a convolutional neural network (CNN) model to analyze rectum dose distribution and predict rectum toxicity. Forty-two cervical cancer patients treated with combined external beam radiotherapy (EBRT) and brachytherapy (BT) were retrospectively collected, including twelve toxicity patients and thirty non-toxicity patients. We adopted a transfer learning strategy to overcome the limited patient data issue. A 16-layers CNN developed by the visual geometry group (VGG-16) of the University of Oxford was pre-trained on a large-scale natural image database, ImageNet, and fine-tuned with patient rectum surface dose maps (RSDMs), which were accumulated EBRT  +  BT doses on the unfolded rectum surface. We used the adaptive synthetic sampling approach and the data augmentation method to address the two challenges, data imbalance and data scarcity. The gradient-weighted class activation maps (Grad-CAM) were also generated to highlight the discriminative regions on the RSDM along with the prediction model. We compare different CNN coefficients fine-tuning strategies, and compare the predictive performance using the traditional dose volume parameters, e.g. D _0.1/1/2cc, and the texture features extracted from the RSDM. Satisfactory prediction performance was achieved with the proposed scheme, and we found that the mean Grad-CAM over the toxicity patient group has geometric consistence of distribution with the statistical analysis result, which indicates possible rectum toxicity location. The evaluation results have demonstrated the feasibility of building a CNN-based rectum dose-toxicity prediction model with transfer learning for cervical cancer radiotherapy.

Local Protocol Variations for Image Guided Radiation Therapy in the Multicenter Dutch Hypofractionation (HYPRO) Trial: Impact of Rectal Balloon and MRI Delineation on Anorectal Dose and Gastrointestinal Toxicity Levels

PURPOSE

The phase 3 HYpofractionated irradiation for PROstate cancer (HYPRO) trial randomized patients with intermediate- to high-risk localized prostate cancer to conventionally fractionated (78 Gy in 39 fractions) or hypofractionated (64.6 Gy in 19 fractions) radiation therapy. Differences in techniques and treatment protocols were present between participating centers. This study aimed to compare dose parameters and patient-reported gastrointestinal symptoms between these centers.

METHODS AND MATERIALS

From the trial population, we selected patients (N=572) from 4 treatment centers who received image guided (IG) intensity modulated radiation therapy (IMRT). Center A (n=242) applied planning target volume (PTV) margins of 5 to 6 mm and was considered the reference center. In center B (n=170, 7-mm margins), magnetic resonance imaging (MRI) was integrated in treatment planning. An endorectal balloon (ERB) was applied in center C (n=85, 7-mm margins). Center D (n=75) applied the largest PTV margins of 8 mm. The study protocol provided identical anorectal dose constraints, and local protocols were applied for further treatment optimization. Anorectal dose-surface histograms were compared by applying t tests. Rectal complaints during follow-up (6 months to 4 years) were compared in a generalized linear model, adjusting for age, follow-up, treatment arm, and hormone therapy.

RESULTS

Favorable anorectal dose distributions were found for center B (MRI delineation) and center C (ERB application) as compared with centers A and D. These were associated with significantly lower incidences of patient-reported complaints of rectal incontinence, use of incontinence pads, and rectal discomfort in these centers. Furthermore, lower incidences of increased stool frequency (≥4 per day) and mucous loss were observed for center C.

CONCLUSIONS

Despite comparable IG-IMRT techniques and predefined dose constraints, pronounced differences in dose distributions and toxicity rates were observed. MRI delineation and ERB application were associated with favorable rectal dose parameters and toxicity profiles, whereas a 2- to 3-mm difference in PTV margins did not translate into observed differences. We conclude that choices for treatment optimization of IG-IMRT are important and clinically relevant for patients since these affect symptoms experienced in daily life.

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.

Delivered dose can be a better predictor of rectal toxicity than planned dose in prostate radiotherapy

Background and purpose

For the first time, delivered dose to the rectum has been calculated and accumulated throughout the course of prostate radiotherapy using megavoltage computed tomography (MVCT) image guidance scans. Dosimetric parameters were linked with toxicity to test the hypothesis that delivered dose is a stronger predictor of toxicity than planned dose.

Material and methods

Dose–surface maps (DSMs) of the rectal wall were automatically generated from daily MVCT scans for 109 patients within the VoxTox research programme. Accumulated-DSMs, representing total delivered dose, and planned-DSMs, from planning CT data, were parametrised using Equivalent Uniform Dose (EUD) and ‘DSM dose-width’, the lateral dimension of an ellipse fitted to a discrete isodose cluster. Associations with 6 toxicity endpoints were assessed using receiver operator characteristic curve analysis.

Results

For rectal bleeding, the area under the curve (AUC) was greater for accumulated dose than planned dose for DSM dose-widths up to 70 Gy. Accumulated 65 Gy DSM dose-width produced the strongest spatial correlation (AUC 0.664), while accumulated EUD generated the largest AUC overall (0.682). For proctitis, accumulated EUD was the only reportable predictor (AUC 0.673). Accumulated EUD was systematically lower than planned EUD.

Conclusions

Dosimetric parameters extracted from accumulated DSMs have demonstrated stronger correlations with rectal bleeding and proctitis, than planned DSMs.

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.

Comparison of different contouring definitions of the rectum as organ at risk (OAR) and dose-volume parameters predicting rectal inflammation in radiotherapy of prostate cancer: which definition to use?

OBJECTIVE

The objective of this retrospective planning study was to find a contouring definition for the rectum as an organ at risk (OAR) in curative three-dimensional external beam radiotherapy (EBRT) for prostate cancer (PCa) with a predictive correlation between the dose-volume histogram (DVH) and rectal toxicity.

METHODS

In a pre-study, the planning CT scans of 23 patients with PCa receiving definitive EBRT were analyzed. The rectum was contoured according to 13 different definitions, and the dose distribution was correlated with the respective rectal volumes by generating DVH curves. Three definitions were identified to represent the most distinct differences in the shapes of the DVH curves: one anatomical definition recommended by the Radiation Therapy Oncology Group (RTOG) and two functional definitions based on the target volume. In the main study, the correlation between different relative DVH parameters derived from these three contouring definitions and the occurrence of rectal toxicity during and after EBRT was studied in two consecutive collectives. The first cohort consisted of 97 patients receiving primary curative EBRT and the second cohort consisted of 66 patients treated for biochemical recurrence after prostatectomy. Rectal toxicity was investigated by clinical investigation and scored according to the Common Terminology Criteria for Adverse Events. Candidate parameters were the volume of the rectum, mean dose, maximal dose, volume receiving at least 60 Gy (V₆₀), area under the DVH curve up to 25 Gy and area under the DVH curve up to 75 Gy in dependence of each chosen rectum definition. Multivariable logistic regression considered other clinical factors such as pelvine lymphatics vs local target volume, diabetes, prior rectal surgery, anticoagulation or haemorrhoids too.

RESULTS

In Cohort 1 (primary EBRT), the mean rectal volumes for definitions "RTOG", planning target volume "(PTV)-based" and "PTV-linked" were 100 cm³ [standard deviation (SD) 43 cm³], 60 cm³ (SD 26 cm³) and 74 cm³ (SD 31 cm³), respectively (p < 0.01; analysis of variance). The mean rectal doses according to these definitions were 35 Gy (SD 8 Gy), 48 Gy (SD 4 Gy) and 44 Gy (SD 5 Gy) (p < 0.01). In Cohort 2 (salvage EBRT), the mean rectal volumes were 114 cm³ (SD 47 cm³), 64 cm³ (SD 26 cm³) and 81 cm³ (SD 30 cm³) (p < 0.01) and the mean doses received by the rectum were 36 Gy (SD 8 Gy), 49 Gy (SD 5 Gy) and 44 Gy (SD 5 Gy) (p < 0.01). Acute or subacute rectal inflammation occurred in 69 (71.9%) patients in Cohort 1 and in 43 (70.5%) in Cohort 2. We did not find a correlation between all investigated DVH parameters and rectal toxicity, irrespective of the investigated definition. By adding additional variables in multivariate analysis, the predictive ability was substantially improved. Still, there was essentially no difference in the probability of predicting rectal inflammation occurrence between the tested contouring definitions.

CONCLUSION

The RTOG anatomy-based recommendations are questionable in comparison with functional definitions, as they result in higher variances in several relative DVH parameters. Moreover, the anatomy-based definition is no better and no worse in the predictive value concerning clinical end points. Advances in knowledge: Functional definitions for the rectum as OAR are easier to apply, faster to contour, have smaller variances and do not offer less information than the anatomy-based RTOG definition.

Prostate external beam radiotherapy combined with high-dose-rate brachytherapy: dose-volume parameters from deformably-registered plans correlate with late gastrointestinal complications

BACKGROUND

Derivation of dose-volume correlated with toxicity for multi-modal treatments can be difficult due to the perceived need for voxel-by-voxel dose accumulation. With data available for a single-institution cohort with long follow-up, an investigation was undertaken into rectal dose-volume effects for gastrointestinal toxicities after deformably-registering each phase of a combined external beam radiotherapy (EBRT)/high-dose-rate (HDR) brachytherapy prostate treatment.

METHODS

One hundred and eighteen patients received EBRT in 23 fractions of 2 Gy and HDR (TG43 algorithm) in 3 fractions of 6.5 Gy. Results for the Late Effects of Normal Tissues - Subjective, Objective, Management and Analytic toxicity assessments were available with a median follow-up of 72 months. The HDR CT was deformably-registered to the EBRT CT. Doses were corrected for dose fractionation. Rectum dose-volume histogram (DVH) parameters were calculated in two ways. (1) Distribution-adding: parameters were calculated after the EBRT dose distribution was 3D-summed with the registered HDR dose distribution. (2) Parameter-adding: the EBRT DVH parameters were added to HDR DVH parameters. Logistic regressions and Mann-Whitney U-tests were used to correlate parameters with late peak toxicity (dichotomised at grade 1 or 2).

RESULTS

The 48-80, 40-63 and 49-55 Gy dose regions from distribution-adding were significantly correlated with rectal bleeding, urgency/tenesmus and stool frequency respectively. Additionally, urgency/tenesmus and anorectal pain were associated with the 25-26 Gy and 44-48 Gy dose regions from distribution-adding respectively. Parameter-adding also indicated the low-mid dose region was significantly correlated with stool frequency and proctitis.

CONCLUSIONS

This study confirms significant dose-histogram effects for gastrointestinal toxicities after including deformable registration to combine phases of EBRT/HDR prostate cancer treatment. The findings from distribution-adding were in most cases consistent with those from parameter-adding. The mid-high dose range and near maximum doses were important for rectal bleeding. The distribution-adding mid-high dose range was also important for stool frequency and urgency/tenesmus. We encourage additional studies in a variety of institutions using a variety of dose accumulation methods with appropriate inter-fraction motion management.

TRIAL REGISTRATION

NCT NCT00193856 . Retrospectively registered 12 September 2005.

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.

Identification of a rectal subregion highly predictive of rectal bleeding in prostate cancer IMRT

BACKGROUND AND PURPOSE

To identify rectal subregions at risks (SRR) highly predictive of 3-year rectal bleeding (RB) in prostate cancer IMRT.

MATERIALS AND METHODS

Overall, 173 prostate cancer patients treated with IMRT/IGRT were prospectively analyzed, divided into "training" (n=118) and "validation" cohorts (n=53). Dose-volume histograms (DVHs) were calculated in three types of rectal subregions: "geometric", intuitively defined (hemi-rectum,…); "personalized", obtained by non-rigid registration followed by voxel-wise statistical analysis (SRRp); "generic", mapped from SRRps, located within 8×8 rectal subsections (SRRg). DVHs from patients with and without RB were compared and used for toxicity prediction.

RESULTS

Training cohort SRRps were primarily within the inferior anterior hemi-rectum and upper anal canal, with 3.8Gy mean dose increase for Grade⩾1 RB patients. The SRRg, representing 15% of the absolute rectal volume, was located in 10 inferior-anterior rectal subsections. V18-V70 for SRRps and V58-V65 for SRRg were significantly higher for RB patients than non-RB. Maximum areas under the curve (AUCs) for SRRp and SRRg RB prediction were 71% and 64%, respectively. The validation cohort confirmed the predictive value of SRRg for Grade⩾1 RB. The total cohort confirmed the predictive value of SRRg for Grade⩾2 RB. Geometrical subregions were not RB predictors.

CONCLUSION

The inferior-anterior hemi anorectum was highly predictive of RB.

A Voxel-Based Approach to Explore Local Dose Differences Associated With Radiation-Induced Lung Damage

PURPOSE

To apply a voxel-based (VB) approach aimed at exploring local dose differences associated with late radiation-induced lung damage (RILD).

METHODS AND MATERIALS

An interinstitutional database of 98 patients who were Hodgkin lymphoma (HL) survivors treated with postchemotherapy supradiaphragmatic radiation therapy was analyzed in the study. Eighteen patients experienced late RILD, classified according to the Radiation Therapy Oncology Group scoring system. Each patient's computed tomographic (CT) scan was normalized to a single reference case anatomy (common coordinate system, CCS) through a log-diffeomorphic approach. The obtained deformation fields were used to map the dose of each patient into the CCS. The coregistration robustness and the dose mapping accuracy were evaluated by geometric and dose scores. Two different statistical mapping schemes for nonparametric multiple permutation inference on dose maps were applied, and the corresponding P<.05 significance lung subregions were generated. A receiver operating characteristic (ROC)-based test was performed on the mean dose extracted from each subregion.

RESULTS

The coregistration process resulted in a geometrically robust and accurate dose warping. A significantly higher dose was consistently delivered to RILD patients in voxel clusters near the peripheral medial-basal portion of the lungs. The area under the ROC curves (AUC) from the mean dose of the voxel clusters was higher than the corresponding AUC derived from the total lung mean dose.

CONCLUSIONS

We implemented a framework including a robust registration process and a VB approach accounting for the multiple comparison problem in dose-response modeling, and applied it to a cohort of HL survivors to explore a local dose-RILD relationship in the lungs. Patients with RILD received a significantly greater dose in parenchymal regions where low doses (∼6 Gy) were delivered. Interestingly, the relation between differences in the high-dose range and RILD seems to lack a clear spatial signature.

Normal tissue complication probability (NTCP) models for late rectal bleeding, stool frequency and fecal incontinence after radiotherapy in prostate cancer patients

BACKGROUND AND PURPOSE

Curative radiotherapy for prostate cancer may lead to anorectal side effects, including rectal bleeding, fecal incontinence, increased stool frequency and rectal pain. The main objective of this study was to develop multivariable NTCP models for these side effects.

MATERIAL AND METHODS

The study sample was composed of 262 patients with localized or locally advanced prostate cancer (stage T1-3). Anorectal toxicity was prospectively assessed using a standardized follow-up program. Different anatomical subregions within and around the anorectum were delineated. A LASSO logistic regression analysis was used to analyze dose volume effects on toxicity.

RESULTS

In the univariable analysis, rectal bleeding, increase in stool frequency and fecal incontinence were significantly associated with a large number of dosimetric parameters. The collinearity between these predictors was high (VIF>5). In the multivariable model, rectal bleeding was associated with the anorectum (V70) and anticoagulant use, fecal incontinence was associated with the external sphincter (V15) and the iliococcygeal muscle (V55). Finally, increase in stool frequency was associated with the iliococcygeal muscle (V45) and the levator ani (V40). No significant associations were found for rectal pain.

CONCLUSIONS

Different anorectal side effects are associated with different anatomical substructures within and around the anorectum. The dosimetric variables associated with these side effects can be used to optimize radiotherapy treatment planning aiming at prevention of specific side effects and to estimate the benefit of new radiation technologies.