Search

Background: Electromagnetic Tracking (EMT) technology has been integrated in a prototype high-dose-rate brachytherapy (HDR-BT) afterloading device. Its potential for dwell position (DP) monitoring has earlier been demonstrated in prostate phantoms. However, its performance for prostate BT in the clinical setting remains to be assessed., Aim: Assess the reliability and value of EMT measurements in transrectal ultrasound-based (TRUS-based) and computed tomography-based (CT-based) prostate HDR-BT., Methods: EMT measurements were conducted on 20 patients undergoing dual-fraction prostate HDR-BT monotherapy. In each treatment fraction an individual TRUS-based or CT-based treatment plan was generated. The measurements were compared to DPs of manually reconstructed needles in those TRUS-based or CT-based treatment plans. An internal reference sensor was also placed in one needle to assess internal movement levels and its potential for movement correction., Results: For TRUS-based treatments, median Euclidean distances (ED) of 1.00 mm were observed between EMT measurements and manual DP determination. Reference sensor movement was minimal at a median of 0.18 mm. For DPs measured in the CT-room and treatment room, median EDs of 1.60 mm and 2.24 mm compared to CT-based DP determination respectively were observed, indicating the system's ability to detect changes in implant geometry over time and after patient repositioning. Median reference sensor movement of 0.97 mm was observed. Implementing reference sensor-based movement correction led to a significant but small decrease in ED for CT-based treatments., Conclusion: EMT is suitable for TRUS-based prostate HDR-BT quality assurance and error detection. While EMT can identify changes in implant geometry in CT-based prostate HDR-BT treatments, it showed lower accuracy in this setting., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier B.V. All rights reserved.)

Purpose: In evaluating second primary cancers (SPCs) following External Beam Radiotherapy (EBRT), the role of lifestyle factors is frequently not considered due to data limitations. We investigated the association between smoking, comorbidities, and SPC risks within EBRT-treated patients for localized prostate cancer (PCa)., Patients & Methods: The study included 1,883 PCa survivors aged 50-79, treated between 2006 and 2013, with intensity-modulated radiotherapy (IMRT) or three-dimensional conformal radiotherapy (3D-CRT). Clinical data were combined with SPC and survival data from the Netherlands Cancer Registry with a 12-month latency period. Standardized Incidence Ratios (SIRs) were calculated comparing the EBRT cohort with the general Dutch population. To explore the effect of patient and treatment characteristics on SPCs we conducted a Cox regression analysis. Lastly, we estimated cumulative incidences of developing solid SPC, pelvis SPC, and non-pelvis SPC using a competing risk analysis., Results: Significantly increased SIRs were observed for all SPC (SIR = 1.21, 95% confidence interval [CI]: 1.08-1.34), pelvis SPC (SIR = 1.46, 95% CI: 1.18-1.78), and non-pelvis SPC (SIR = 1.18, 95% CI [1.04-1.34]). Smoking status was significantly associated with pelvic and non-pelvic SPCs. Charlson comorbidity index (CCI) ≥ 1 (Hazard Ratio [HR] = 1.45, 95% CI: 1.10-1.91), cardiovascular disease (HR = 1.41, 95% CI: 1.05-1.88), and chronic obstructive pulmonary disease (COPD) (HR = 1.91, 95% CI: 1.30-2.79) were significantly associated with non-pelvis SPC. The proportion of active smoking numbers in the cohort was similar to the general population., Interpretation: We conclude that the presence of comorbidities in the EBRT population might be a relevant factor in observed excess non-pelvis SPC risk, but not for excess pelvis SPC risk.

Background: Modelling studies suggest that advanced intensity-modulated radiotherapy may increase second primary cancer (SPC) risks, due to increased radiation exposure of tissues located outside the treatment fields. In the current study we investigated the association between SPC risks and characteristics of applied external beam radiotherapy (EBRT) protocols for localized prostate cancer (PCa)., Methods: We collected EBRT protocol characteristics (2000-2016) from five Dutch RT institutes for the 3D-CRT and advanced EBRT era (N = 7908). From the Netherlands Cancer Registry we obtained patient/tumour characteristics, SPC data, and survival information. Standardized incidence ratios (SIR) were calculated for pelvis and non-pelvis SPC. Nationwide SIRs were calculated as a reference, using calendar period as a proxy to label 3D-CRT/advanced EBRT., Results: From 2000-2006, 3D-CRT with 68-78 Gy in 2 Gy fractions, delivered with 10-23 MV and weekly portal imaging was the most dominant protocol. By the year 2010 all institutes routinely used advanced EBRT (IMRT, VMAT, tomotherapy), mainly delivering 78 Gy in 2 Gy fractions, using various kV/MV imaging protocols. Sixteen percent (N = 1268) developed ≥ 1 SPC. SIRs for pelvis and non-pelvis SPC (all institutes, advanced EBRT vs 3D-CRT) were 1.17 (1.00-1.36) vs 1.39 (1.21-1.59), and 1.01 (0.89-1.07) vs 1.03 (0.94-1.13), respectively. Nationwide non-pelvis SIR was 1.07 (1.01-1.13) vs 1.02 (0.98-1.07). Other RT protocol characteristics did not correlate with SPC endpoints., Conclusion: None of the studied RT characteristics of advanced EBRT was associated with increased out-of-field SPC risks. With constantly evolving EBRT protocols, evaluation of associated SPC risks remains important., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

Objective . In radiotherapy, the internal movement of organs between treatment sessions causes errors in the final radiation dose delivery. To assess the need for adaptation, motion models can be used to simulate dominant motion patterns and assess anatomical robustness before delivery. Traditionally, such models are based on principal component analysis (PCA) and are either patient-specific (requiring several scans per patient) or population-based, applying the same set of deformations to all patients. We present a hybrid approach which, based on population data, allows to predict patient-specific inter-fraction variations for an individual patient. Approach . We propose a deep learning probabilistic framework that generates deformation vector fields warping a patient's planning computed tomography (CT) into possible patient-specific anatomies. This daily anatomy model (DAM) uses few random variables capturing groups of correlated movements. Given a new planning CT, DAM estimates the joint distribution over the variables, with each sample from the distribution corresponding to a different deformation. We train our model using dataset of 312 CT pairs with prostate, bladder, and rectum delineations from 38 prostate cancer patients. For 2 additional patients (22 CTs), we compute the contour overlap between real and generated images, and compare the sampled and 'ground truth' distributions of volume and center of mass changes. Results . With a DICE score of 0.86 ± 0.05 and a distance between prostate contours of 1.09 ± 0.93 mm, DAM matches and improves upon previously published PCA-based models, using as few as 8 latent variables. The overlap between distributions further indicates that DAM's sampled movements match the range and frequency of clinically observed daily changes on repeat CTs. Significance . Conditioned only on planning CT values and organ contours of a new patient without any pre-processing, DAM can accurately deformations seen during following treatment sessions, enabling anatomically robust treatment planning and robustness evaluation against inter-fraction anatomical changes., (Creative Commons Attribution license.)

A review of studies on seminal vesicle motion was performed to improve the understanding of these treatment uncertainties. This will aid planning target volume margin reduction, which is necessary for hypofractionation of high-risk prostate cancer. Embase, Medline, Web of science Core collection, Cochrane CENTRAL register of trials and Google scholar were searched for publications including 3D information on seminal vesicle motion. In total 646 publications were found of which 22 publications were eligible for inclusion. The mean, systematic and random error of inter- and intra-fraction translations are reported, as well as rotations. The translations of the seminal vesicles is smallest in the left-right direction, whereas the rotation was largest around this axis. Although rectal and bladder filling status were the main cause for seminal vesicle motion, no apparent effect on magnitude of motion was seen when different bladder and rectal preparation protocols were used. Inter- and intra-fraction motion of the seminal vesicles is significant. In the studies, systematic and random errors range between 1-7 mm and 1-5 mm respectively, and are largely uncorrelated to prostate motion. The maximum correlation between seminal vesicle and prostate motion was reported with an R 2 of 0.7, while 3 other studies report lower and/or non-significant correlations. Five studies report a planning target volume margin of approximately 8 mm. This margin is in line with the results of four relevant dosimetric studies. Mitigating the inter- and intra-fraction motion of the seminal vesicles, including prostate tracking, has the potential to reduce planning target volume margins., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

Background/purpose: Intensity-modulated proton therapy (IMPT) is highly sensitive to anatomical variations which can cause inadequate target coverage during treatment. Available mitigation techniques include robust treatment planning and online-adaptive IMPT. This study compares a robust planning strategy to two online-adaptive IMPT strategies to determine the benefit of online adaptation., Materials/methods: We derived the robustness settings and safety margins needed to yield adequate target coverage (V 95% ≥98%) for >90% of 11 patients in a prostate cancer cohort (88 repeat CTs). For each patient, we also adapted a non-robust prior plan using a simple restoration and a full adaptation method. The restoration uses energy-adaptation followed by a fast spot-intensity re-optimization. The full adaptation uses energy-adaptation followed by the addition of new spots and a range-robust spot-intensity optimization. Dose was prescribed as 55 Gy(RBE) to the low-dose target (lymph nodes and seminal vesicles) with a boost to 74 Gy(RBE) to the high-dose target (prostate). Daily patient set-up was simulated using implanted intra-prostatic markers., Results: Margins of 4 and 8 mm around the high- and low-dose target regions, a 6 mm setup error and a 3% range error were found to obtain adequate target coverage for all repeat CTs of 10/11 patients (94.3% of all 88 repeat CTs). Both online-adaptive strategies yielded V 95% ≥98% and better OAR sparing in 11/11 patients. Median OAR improvements up to 11%-point and 16%-point were observed when moving from robust planning to respectively restoration and full adaption., Conclusion: Both full plan adaptation and simple dose restoration can increase OAR sparing besides better conforming to the target criteria compared to robust treatment planning., (Copyright © 2020 The Authors. Published by Elsevier B.V. All rights reserved.)

Purpose: To develop and validate a robust and accurate registration pipeline for automatic contour propagation for online adaptive Intensity-Modulated Proton Therapy (IMPT) of prostate cancer using elastix software and deep learning., Methods: A three-dimensional (3D) Convolutional Neural Network was trained for automatic bladder segmentation of the computed tomography (CT) scans. The automatic bladder segmentation alongside the computed tomography (CT) scan is jointly optimized to add explicit knowledge about the underlying anatomy to the registration algorithm. We included three datasets from different institutes and CT manufacturers. The first was used for training and testing the ConvNet, where the second and the third were used for evaluation of the proposed pipeline. The system performance was quantified geometrically using the dice similarity coefficient (DSC), the mean surface distance (MSD), and the 95% Hausdorff distance (HD). The propagated contours were validated clinically through generating the associated IMPT plans and compare it with the IMPT plans based on the manual delineations. Propagated contours were considered clinically acceptable if their treatment plans met the dosimetric coverage constraints on the manual contours., Results: The bladder segmentation network achieved a DSC of 88% and 82% on the test datasets. The proposed registration pipeline achieved a MSD of 1.29 ± 0.39, 1.48 ± 1.16, and 1.49 ± 0.44 mm for the prostate, seminal vesicles, and lymph nodes, respectively, on the second dataset and a MSD of 2.31 ± 1.92 and 1.76 ± 1.39 mm for the prostate and seminal vesicles on the third dataset. The automatically propagated contours met the dose coverage constraints in 86%, 91%, and 99% of the cases for the prostate, seminal vesicles, and lymph nodes, respectively. A Conservative Success Rate (CSR) of 80% was obtained, compared to 65% when only using intensity-based registration., Conclusion: The proposed registration pipeline obtained highly promising results for generating treatment plans adapted to the daily anatomy. With 80% of the automatically generated treatment plans directly usable without manual correction, a substantial improvement in system robustness was reached compared to a previous approach. The proposed method therefore facilitates more precise proton therapy of prostate cancer, potentially leading to fewer treatment-related adverse side effects., (© 2019 The Authors. Medical Physics published by Wiley Periodicals, Inc. on behalf of American Association of Physicists in Medicine.)

Proton therapy plans are very sensitive to anatomical changes such as density changes along the pencil-beam paths and changes in organ shape and location. Previously, we developed a restoration method which compensates for density changes along the pencil-beam paths but which is unable to adapt for anatomical changes. This study's purpose is to develop and evaluate an automated method for adaptation of IMPT plans in near real-time to the anatomy of the day. We developed an automated treatment plan adaptation method using (1) a restoration of spot positions (Bragg peaks) by adapting the energies to the new water equivalent path lengths; and (2) a spot addition to fully cover the target of the day, followed by a fast reference point method optimization of the spot weights resulting in a Pareto optimal plan for the daily anatomy. The method was developed and evaluated using 8-10 repeat CT scans of 11 prostate cancer patients, prescribing 55 Gy(RBE) (seminal vesicles and lymph nodes) with a boost to 74 Gy(RBE) (prostate). Applying the automated adaptation method resulted in a clinically acceptable target coverage (V 95% [Formula: see text] 98% and V 107% [Formula: see text] 2%) for 96% of the scans after a single iteration of adding 2500 spots. The other scans obtained target coverages with V 95% [Formula: see text] 98% and 2  <  V 107% [Formula: see text] 5%. When using two spot-addition iterations, all scans obtained clinically acceptable results. Compared to the restoration method the adaptation lowered the mean dose to rectum and bladder with median values of 6.2 Gy(RBE) and 4.7 Gy(RBE) respectively. The largest improvements were obtained for V 45Gy(RBE) for both rectum and bladder, with median differences of 10.3%-point and 10.8%-point respectively, and maximum differences up to 22%-point. The two adaptation steps took on average 7.3 s and 1.7 min respectively. No user interaction was needed, making this fast and fully automated method a first step towards online adaptive proton therapy.

Proton therapy is very sensitive to daily density changes along the pencil beam paths. The purpose of this study is to develop and evaluate an automated method for adaptation of IMPT plans to compensate for these daily tissue density variations. A two-step restoration method for 'densities-of-the-day' was created: (1) restoration of spot positions (Bragg peaks) by adapting the energy of each pencil beam to the new water equivalent path length; and (2) re-optimization of pencil beam weights by minimizing the dosimetric difference with the planned dose distribution, using a fast and exact quadratic solver. The method was developed and evaluated using 8-10 repeat CT scans of 10 prostate cancer patients. Experiments demonstrated that giving a high weight to the PTV in the re-optimization resulted in clinically acceptable restorations. For all scans we obtained V 95%   ⩾  98% and V 107%   ⩽  2%. For the bladder, the differences between the restored and the intended treatment plan were below  +2 Gy and  +2%-point. The rectum differences were below  +2 Gy and  +2%-point for 90% of the scans. In the remaining scans the rectum was filled with air, which partly overlapped with the PTV. The air cavity distorted the Bragg peak resulting in less favorable rectum doses.

Purpose: To shorten delivery times of intensity modulated proton therapy by reducing the number of energy layers in the treatment plan., Methods and Materials: We have developed an energy layer reduction method, which was implemented into our in-house-developed multicriteria treatment planning system "Erasmus-iCycle." The method consisted of 2 components: (1) minimizing the logarithm of the total spot weight per energy layer; and (2) iteratively excluding low-weighted energy layers. The method was benchmarked by comparing a robust "time-efficient plan" (with energy layer reduction) with a robust "standard clinical plan" (without energy layer reduction) for 5 oropharyngeal cases and 5 prostate cases. Both plans of each patient had equal robust plan quality, because the worst-case dose parameters of the standard clinical plan were used as dose constraints for the time-efficient plan. Worst-case robust optimization was performed, accounting for setup errors of 3 mm and range errors of 3% + 1 mm. We evaluated the number of energy layers and the expected delivery time per fraction, assuming 30 seconds per beam direction, 10 ms per spot, and 400 Giga-protons per minute. The energy switching time was varied from 0.1 to 5 seconds., Results: The number of energy layers was on average reduced by 45% (range, 30%-56%) for the oropharyngeal cases and by 28% (range, 25%-32%) for the prostate cases. When assuming 1, 2, or 5 seconds energy switching time, the average delivery time was shortened from 3.9 to 3.0 minutes (25%), 6.0 to 4.2 minutes (32%), or 12.3 to 7.7 minutes (38%) for the oropharyngeal cases, and from 3.4 to 2.9 minutes (16%), 5.2 to 4.2 minutes (20%), or 10.6 to 8.0 minutes (24%) for the prostate cases., Conclusions: Delivery times of intensity modulated proton therapy can be reduced substantially without compromising robust plan quality. Shorter delivery times are likely to reduce treatment uncertainties and costs., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Purpose: To investigate the dosimetric impact of intrafraction prostate motion and the effect of robot correction strategies for hypofractionated CyberKnife treatments with a simultaneously integrated boost., Methods and Materials: A total of 548 real-time prostate motion tracks from 17 patients were available for dosimetric simulations of CyberKnife treatments, in which various correction strategies were included. Fixed time intervals between imaging/correction (15, 60, 180, and 360 seconds) were simulated, as well as adaptive timing (ie, the time interval reduced from 60 to 15 seconds in case prostate motion exceeded 3 mm or 2° in consecutive images). The simulated extent of robot corrections was also varied: no corrections, translational corrections only, and translational corrections combined with rotational corrections up to 5°, 10°, and perfect rotational correction. The correction strategies were evaluated for treatment plans with a 0-mm or 3-mm margin around the clinical target volume (CTV). We recorded CTV coverage (V100%) and dose-volume parameters of the peripheral zone (boost), rectum, bladder, and urethra., Results: Planned dose parameters were increasingly preserved with larger extents of robot corrections. A time interval between corrections of 60 to 180 seconds provided optimal preservation of CTV coverage. To achieve 98% CTV coverage in 98% of the treatments, translational and rotational corrections up to 10° were required for the 0-mm margin plans, whereas translational and rotational corrections up to 5° were required for the 3-mm margin plans. Rectum and bladder were spared considerably better in the 0-mm margin plans. Adaptive timing did not improve delivered dose., Conclusions: Intrafraction prostate motion substantially affected the delivered dose but was compensated for effectively by robot corrections using a time interval of 60 to 180 seconds. A 0-mm margin required larger extents of additional rotational corrections than a 3-mm margin but resulted in lower doses to rectum and bladder., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Background and Purpose: Deformation and correlated target motion remain challenges for margin recipes in radiotherapy (RT). This study presents a statistical deformable motion model for multiple targets and applies it to margin evaluations for locally advanced prostate cancer i.e. RT of the prostate (CTV-p), seminal vesicles (CTV-sv) and pelvic lymph nodes (CTV-ln)., Material and Methods: The 19 patients included in this study, all had 7-10 repeat CT-scans available that were rigidly aligned with the planning CT-scan using intra-prostatic implanted markers, followed by deformable registrations. The displacement vectors from the deformable registrations were used to create patient-specific statistical motion models. The models were applied in treatment simulations to determine probabilities for adequate target coverage, e.g. by establishing distributions of the accumulated dose to 99% of the target volumes (D99) for various CTV-PTV expansions in the planning-CTs., Results: The method allowed for estimation of the expected accumulated dose and its variance of different DVH parameters for each patient. Simulations of inter-fractional motion resulted in 7, 10, and 18 patients with an average D99 >95% of the prescribed dose for CTV-p expansions of 3mm, 4mm and 5mm, respectively. For CTV-sv and CTV-ln, expansions of 3mm, 5mm and 7 mm resulted in 1, 11 and 15 vs. 8, 18 and 18 patients respectively with an average D99 >95% of the prescription., Conclusions: Treatment simulations of target motion revealed large individual differences in accumulated dose mainly for CTV-sv, demanding the largest margins whereas those required for CTV-p and CTV-ln were comparable., (Copyright © 2013 Elsevier Ireland Ltd. All rights reserved.)

Unlabelled: Daily treatment plan selection from a plan library is a major adaptive radiotherapy strategy to account for individual internal anatomy variations. This strategy depends on the initial input images being representative for the variations observed later in the treatment course. Focusing on locally advanced prostate cancer, our aim was to evaluate if residual motion of the prostate (CTV-p) and the elective targets (CTV-sv, CTV-ln) can be prospectively accounted for with a statistical deformable model based on images acquired in the initial part of treatment., Methods: Thirteen patients with locally advanced prostate cancer, each with 9-10 repeat CT scans, were included. Displacement vectors fields (DVF) obtained from contour-based deformable registration of delineations in the repeat- and planning CT scans were used to create patient-specific statistical motion models using principal component analysis (PCA). For each patient and CTV, four PCA-models were created: one with all 9-10 DVF as input in addition to models with only four, five or six DVFs as input. Simulations of target shapes from each PCA-model were used to calculate iso-coverage levels, which were converted to contours. The levels were analyzed for sensitivity and precision., Results: A union of the simulated shapes was able to cover at least 97%, 97% and 95% of the volumes of the evaluated CTV shapes for PCA-models using six, five and four DVFs as input, respectively. There was a decrease in sensitivity with higher iso-coverage levels, with a sharper decline for greater target movements. Apart from having the steepest decline in sensitivity, CTV-sv also displayed the greatest influence on the number of geometries used in the PCA-model., Conclusions: PCA-based simulations of residual motion derived from four to six DVFs as input could account for the majority of the target shapes present during the latter part of the treatment. CTV-sv displayed the greatest range in both sensitivity and precision.

Background: There is growing evidence that prostate cancer (PC) cells are more sensitive to high fraction dose in hypofractionation schemes. High-dose-rate (HDR) brachytherapy as monotherapy is established to be a good treatment option for PC using extremely hypofractionated schemes. This hypofractionation can also be achieved with stereotactic body radiotherapy (SBRT). We report results on toxicity, PSA response, and quality of life (QOL) in patients treated with SBRT for favorable-risk PC., Methods: Over the last 4 years, 50 hormone-naïve patients with low- and intermediate-risk PC were treated with SBRT to a total dose of 38 Gy delivered in four daily fractions of 9.5 Gy. An integrated boost to 11 Gy per fraction was applied to the dominant lesion if visible on MRI. Toxicity and QoL was assessed prospectively using validated questionnaires., Results: Median follow-up was 23 months. The 2-year actuarial biochemical control rate was 100%. Median PSA nadir was 0.6 ng/ml. Median International Prostate Symptoms Score (IPSS) was 9/35 before treatment, with a median increase of 4 at 3 months and remaining stable at 13/35 thereafter. The EORTC/RTOG toxicity scales showed grade 2 and 3 gastrointestinal (GI) acute toxicity in 12% and 2%, respectively. The late grade 2 GI toxicity was 3% during 24 months FU. Genitourinary (GU) grade 2, 3 toxicity was seen in 15%, 8%, in the acute phase and 10%, 6% at 24 months, respectively. The urinary, bowel and sexual domains of the EORTC-PR25 scales recovered over time, showing no significant changes at 24 months post-treatment., Conclusions: SBRT to 38 Gy in 4 daily fractions for low- and intermediate-risk PC patients is feasible with low acute and late genitourinary and gastrointestinal toxicity. Longer follow-up preferably within randomized studies, is required to compare these results with standard fractionation schemes.

Purpose: To develop a method for margin evaluation accounting for all measured displacements during treatment of prostate cancer., Methods and Materials: For 21 patients treated with stereographic targeting marker-based online translation corrections, dose distributions with varying margins and gradients were created. Sets of possible cumulative delivered dose distributions were simulated by moving voxels and accumulating dose per voxel. Voxel motion was simulated consistent with measured distributions of systematic and random displacements due to stereographic targeting inaccuracies, deformation, rotation, and intrafraction motion. The method of simulation maintained measured correlation of voxel motions due to organ deformation., Results: For the clinical target volume including prostate and seminal vesicles (SV), the probability that some part receives <95% of the prescribed dose, the changes in minimum dose, and volume receiving 95% of prescription dose compared with planning were 80.5% ± 19.2%, 9.0 ± 6.8 Gy, and 3.0% ± 3.7%, respectively, for the smallest studied margins (3 mm prostate, 5 mm SV) and steepest dose gradients. Corresponding values for largest margins (5 mm prostate, 8 mm SV) with a clinical intensity-modulated radiotherapy dose distribution were 46.5% ± 34.7%, 6.7 ± 5.8 Gy, and 1.6% ± 2.3%. For prostate-only clinical target volume, the values were 51.8% ± 17.7%, 3.3 ± 1.6 Gy, and 0.6% ± 0.5% with the smallest margins and 5.2% ± 7.4%, 1.8 ± 0.9 Gy, and 0.1% ± 0.1% for the largest margins. Addition of three-dimensional rotation corrections only improved these values slightly. All rectal planning constraints were met in the actual reconstructed doses for all studied margins., Conclusion: We developed a system for margin validation in the presence of deformations. In our population, a 5-mm margin provided sufficient dosimetric coverage for the prostate. In contrast, an 8-mm SV margin was still insufficient owing to deformations. Addition of three-dimensional rotation corrections was of minor influence., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Purpose: To quantify the residual geometric uncertainties after on-line corrections with intraprostatic fiducial markers, this study analyzed the deformation of the prostate and, in particular, the seminal vesicles relative to such markers., Patients and Methods: A planning computed tomography (CT) scan and three repeat CT scans were obtained for 21 prostate cancer patients who had had three to four cylindrical gold markers placed. The prostate and whole seminal vesicles (clinical target volume [CTV]) were delineated on each scan at a slice thickness of 1.5 mm. Rigid body transformations (translation and rotation) mapping the markers onto the planning scan positions were obtained. The translation only (T(only)) or both translation and rotation were applied to the delineated CTVs. Next, the residue CTV surface displacements were determined using nonrigid registration of the delineated contours. For translation and rotation of the CTV, the residues represented deformation; for T(only), the residues stemmed from deformation and rotation. T(only) represented the residues for most currently applied on-line protocols. The patient and population statistics of the CTV surface displacements were calculated. The intraobserver delineation variation was similarly quantified using repeat delineations for all patients and corrected for., Results: The largest CTV deformations were observed at the anterior and posterior side of the seminal vesicles (population average standard deviation

Purpose: To analyze the correlation between dose-volume parameters of the corpora cavernosa and erectile dysfunction (ED) after external beam radiotherapy (EBRT) for prostate cancer., Methods and Materials: Between June 1997 and February 2003, a randomized dose-escalation trial comparing 68 Gy and 78 Gy was conducted. Patients at our institute were asked to participate in an additional part of the trial evaluating sexual function. After exclusion of patients with less than 2 years of follow-up, ED at baseline, or treatment with hormonal therapy, 96 patients were eligible. The proximal corpora cavernosa (crura), the superiormost 1-cm segment of the crura, and the penile bulb were contoured on the planning computed tomography scan and dose-volume parameters were calculated., Results: Two years after EBRT, 35 of the 96 patients had developed ED. No statistically significant correlations between ED 2 years after EBRT and dose-volume parameters of the crura, the superiormost 1-cm segment of the crura, or the penile bulb were found. The few patients using potency aids typically indicated to have ED., Conclusion: No correlation was found between ED after EBRT for prostate cancer and radiation dose to the crura or penile bulb. The present study is the largest study evaluating the correlation between ED and radiation dose to the corpora cavernosa after EBRT for prostate cancer. Until there is clear evidence that sparing the penile bulb or crura will reduce ED after EBRT, we advise to be careful in sparing these structures, especially when this involves reducing treatment margins.

Purpose: To investigate whether a large rectum filling visible on the planning CT scan was associated with a decrease in freedom from any failure (FFF) and freedom from clinical failure (FFCF) for prostate cancer patients., Methods and Materials: Patients from the Dutch trial (78 Gy vs. 68 Gy) with available acute toxicity data were analyzed (n = 549). A 10-mm margin was applied for the first 68 Gy and 0-5 mm for the 10-Gy boost. The dose in the seminal vesicles (SVs) was prescribed within four treatment groups according to the estimated risk of SV involvement. Two potential risk factors (RFs) for a geometric miss were defined: (1) an anorectal volume > or = 90 cm(3) and > or = 25% of treatment-time diarrhea (RF1); and (2) the mean cross-sectional area of the anorectum (RF2). We tested whether these were significant predictors for FFF and FFCF within each treatment group., Results: Significant results were observed only for patients with a risk of SV involvement > 25% (dose of 68-78 Gy to the SVs, n = 349). We found a decrease in FFF (p = 0.001) and FFCF (p = 0.01) for the 87 patients with RF1 (for RF2, p = 0.02 and p = 0.01, respectively). The estimated decrease in the FFCF rate at 5 years was 15%., Conclusion: Tumor control was significantly decreased in patients with a risk of SV involvement > 25% and at risk of geometric miss. Current image guidance techniques offer several solutions to geometrically optimize the treatment. Additional research is needed to evaluate whether geometric misses can be prevented using these techniques.

Purpose: To determine the planning target volume (PTV) margin for an adaptive radiotherapy procedure that uses five computed tomography (CT) scans to calculate an average prostate position and rectum shape. To evaluate alternative methods to determine an average rectum based on a single delineation., Methods and Materials: Repeat CT scans (8-13) of 19 patients were used. The contoured prostates of the first four scans were matched on the planning CT (pCT) prostate contours. With the resulting translations and rotations the average prostate position was determined. An average rectum was obtained by either averaging the coordinates of corresponding points on the rectal walls or by selecting the "best" rectum or transforming the pCT rectum. Dose distributions were calculated for various expanded average prostates. The remaining CT scans were used to determine the dose received by prostate and rectum during treatment., Results: For the prostate of the pCT scan and a 10-mm margin, all patients received more than 95% of the prescribed dose to 95% of the prostate. For the average prostate, a margin of 7 mm was needed to obtain a similar result (average PTV reduction 30%). The average rectum overestimated the mean dose to the rectum by 0.4 +/- 1.6 Gy, which was better than the pCT rectum (2.1 +/- 3.0 Gy) and the alternative average rectums (1.0 +/- 2.6 Gy and 1.4 +/- 3.2 Gy)., Conclusions: Our adaptive procedure allows for reduction of the PTV margin to 7 mm without decreasing prostate coverage during treatment. For accurate estimation of the rectum dose, rectums need to be delineated and averaged over multiple scans.

Purpose: Late gastrointestinal (GI) toxicity after radiotherapy can be partly explained by late effects of acute toxicity (consequential late damage). We studied whether there is a direct relationship between acute and late GI toxicity., Patients and Methods: A total of 553 evaluable patients from the Dutch dose escalation trial (68 Gy vs. 78 Gy) were included. We defined three outcomes for acute reactions: 1) maximum Radiation Therapy Oncology Group acute toxicity, 2) maximum acute mucous discharge (AMD), and 3) maximum acute proctitis. Within a multivariable model, late endpoints (overall toxicity and five toxicity indicators) were studied as a function of acute toxicity, pretreatment symptoms, and relevant dose parameters., Results: At multivariable analysis, AMD and acute proctitis were strong predictors for overall toxicity, "intermittent bleeding," and "incontinence pads" (p < or = 0.01). For "stools > or =6/day" all three were strong predictors. No significant associations were found for "severe bleeding" and "use of steroids." The predictive power of the dose parameters remained at the same level or became weaker for most late endpoints., Conclusions: Acute GI toxicity is an independent significant predictor of late GI toxicity. This suggests a significant consequential component in the development of late GI toxicity.

Purpose: To analyze whether inclusion of predisposing clinical features in the Lyman-Kutcher-Burman (LKB) normal tissue complication probability (NTCP) model improves the estimation of late gastrointestinal toxicity., Methods and Materials: This study includes 468 prostate cancer patients participating in a randomized trial comparing 68 with 78 Gy. We fitted the probability of developing late toxicity within 3 years (rectal bleeding, high stool frequency, and fecal incontinence) with the original, and a modified LKB model, in which a clinical feature (e.g., history of abdominal surgery) was taken into account by fitting subset specific TD50s. The ratio of these TD50s is the dose-modifying factor for that clinical feature. Dose distributions of anorectal (bleeding and frequency) and anal wall (fecal incontinence) were used., Results: The modified LKB model gave significantly better fits than the original LKB model. Patients with a history of abdominal surgery had a lower tolerance to radiation than did patients without previous surgery, with a dose-modifying factor of 1.1 for bleeding and of 2.5 for fecal incontinence. The dose-response curve for bleeding was approximately two times steeper than that for frequency and three times steeper than that for fecal incontinence., Conclusions: Inclusion of predisposing clinical features significantly improved the estimation of the NTCP. For patients with a history of abdominal surgery, more severe dose constraints should therefore be used during treatment plan optimization.

Purpose: To identify dosimetric variables predictive of acute gastrointestinal (GI) and genitourinary (GU) toxicity and to determine whether hormonal therapy (HT) is independently associated with acute GI and GU toxicity in prostate cancer patients treated with conformal radiotherapy (RT)., Methods and Materials: This analysis was performed on 336 patients participating in a multicenter (four hospitals) randomized trial comparing 68 Gy and 78 Gy. The clinical target volume consisted of the prostate with or without the seminal vesicles, depending on the risk of seminal vesicle involvement. The margin from the clinical target volume to the planning target volume was 1 cm. For these patients, the treatment plan for a total dose of 68 Gy was used, because nearly all toxicity appeared before the onset of the 10-Gy boost. Acute toxicity (<120 days) was scored according to the Radiation Therapy Oncology Group criteria. The dosimetric parameters were obtained from the relative and absolute dose-volume/surface histograms derived from the rectal wall (rectal wall volume receiving > or =5-65 Gy) and the bladder surface (bladder surface receiving > or =5-65 Gy). Additionally, relative and absolute dose-length histograms of the rectum were created, and the lengths of rectum receiving more than a certain dose over the whole circumference (rectal length receiving > or =5-65 Gy) were computed. The clinical variables taken into account for GI toxicity were neoadjuvant HT, hospital, and dose-volume group; for GU toxicity, the variables pretreatment GU symptoms, neoadjuvant HT, and transurethral resection of the prostate were analyzed. The variable neoadjuvant HT was divided into three categories: no HT, short-term neoadjuvant HT (started < or =3 months before RT), and long-term neoadjuvant HT (started >3 months before RT)., Results: Acute GI toxicity Grade 2 or worse was seen in 46% of the patients. Patients with long-term neoadjuvant HT experienced less Grade 2 or worse toxicity (27%) compared with those receiving short-term neoadjuvant HT (50%) and no HT (50%). The volumes of the prostate and seminal vesicles were significantly smaller in both groups receiving neoadjuvant HT compared with those receiving no HT. In multivariate logistic regression analysis, including the two statistically significant clinical variables neoadjuvant HT and hospital, a volume effect was found for the relative, as well as absolute, rectal wall volumes exposed to intermediate and high doses. Of all the length parameters, the relative rectal length irradiated to doses of > or =5 Gy and > or =30 Gy and absolute lengths receiving > or =5-15 and 30 Gy were significant. Acute GU toxicity Grade 2 or worse was reported in 56% of cases. For patients with pretreatment GU symptoms, the rate was 93%. The use of short-term and long-term neoadjuvant HT resulted in more GU toxicity (73% and 71%) compared with no HT (50%). In multivariate analysis, containing the variables pretreatment symptoms and neoadjuvant HT, only the absolute dose-surface histogram parameters (absolute surface irradiated to > or =40, 45, and 65 Gy) were significantly associated with acute GU toxicity., Conclusion: A volume effect was found for acute GI toxicity for relative, as well as absolute, volumes. With regard to acute GU toxicity, an area effect was found, but only for absolute dose-surface histogram parameters. Neoadjuvant HT appeared to be an independent prognostic factor for acute toxicity, resulting in less acute GI toxicity, but more acute GU toxicity. The presence of pretreatment GU symptoms was the most important prognostic factor for GU symptoms during RT.

The purpose of this study was to quantify to what extent relative and absolute bladder dose-volume and dose-surface histograms of the planning CT scan were representative for the actual treatment. We used data of 17 patients, who each received 11 repeat CT scans and a planning CT scan. The repeat CT scans were matched on the planning CT scan by the bony anatomy. Clinical treatment plans were used to evaluate the impact of bladder filling changes on the four histogram types. The impact was quantified by calculating for this patient group the correlation coefficient between the planning histogram and the treatment histogram. We found that the absolute dose-surface histogram was the most representative one for the actual treatment.

Purpose: To study the correlations between the dose distributions in the anorectal region and late GI symptoms in patients treated for localized prostate carcinoma., Methods and Materials: Data from a randomized study were analyzed. In this trial, patients were treated with either rectangular or conformal fields with a dose of 66 Gy. Data concerning GI symptoms were collected from questionnaires of 197 patients. The distributions of the anorectal region were projected on maps, and the dose parameters were calculated. The incidences of complaints were studied as a function of the dose-area parameters and clinical parameters, using a proportional hazard regression model. Finally, we tested a series of dose parameters originating from different parts of the anorectal region., Results: Analyzing the total region, only a statistically significant dose-area effect relation for bleeding was found (p < 0.01). Defining subareas, we found effect relations for bleeding, soiling, fecal incontinence, and mucus loss. For bleeding and mucus loss, the strongest correlation was found for the dose received by the upper 70-80% of the anorectal region (p < 0.01). For soiling and fecal incontinence, we found the strongest association with the dose to the lower 40-50% (p < 0.05)., Conclusion: We found evidence that complaints originate from specific regions of the irradiated lower GI tract. Bleeding and mucus loss are probably related to irradiation of the upper part of the rectum. Soiling and fecal incontinence are more likely related to the dose to the anal canal and the lower part of the rectum.

Background and Purpose: The goal of this work is to develop and evaluate strategies to reduce the uncertainty in the prostate position and rectum shape that arises in the preparation stage of the radiation treatment of prostate cancer., Patients and Methods: Nineteen prostate cancer patients, who were treated with 3-dimensional conformal radiotherapy, received each a planning CT scan and 8-13 repeat CT scans during the treatment period. We quantified prostate motion relative to the pelvic bone by first matching the repeat CT scans on the planning CT scan using the bony anatomy. Subsequently, each contoured prostate, including seminal vesicles, was matched on the prostate in the planning CT scan to obtain the translations and rotations. The variation in prostate position was determined in terms of the systematic, random and group mean error. We tested the performance of two correction strategies to reduce the systematic error due to prostate motion. The first strategy, the pre-treatment strategy, used only the initial rectum volume in the planning CT scan to adjust the angle of the prostate with respect to the left-right (LR) axis and the shape and position of the rectum. The second strategy, the adaptive strategy, used the data of repeat CT scans to improve the estimate of the prostate position and rectum shape during the treatment., Results: The largest component of prostate motion was a rotation around the LR axis. The systematic error (1 SD) was 5.1 degrees and the random error was 3.6 degrees (1 SD). The average LR-axis rotation between the planning and the repeat CT scans correlated significantly with the rectum volume in the planning CT scan (r=0.86, P<0.0001). Correction of the rotational position on the basis of the planning rectum volume alone reduced the systematic error by 28%. A correction, based on the data of the planning CT scan and 4 repeat CT scans reduced the systematic error over the complete treatment period by a factor of 2. When the correction was carried out later in the treatment (based on the data of more scans) the overall reduction was less. For the rectum, the first strategy performed best at the upper anterior side, where a reduction of the anterior-posterior displacement of 30% could be achieved. The systematic error could be reduced by 43% for the whole rectum by using the data of 4 repeat CT scans and the planning CT scan., Conclusions: Both the pre-treatment as well as the adaptive correction strategy reduced the systematic error in the prostate position and rectum position and shape. A smaller systematic error makes it possible to safely reduce the margin around the clinical tumor volume, so that normal tissues can be spared or the prescription dose can be escalated.

Purpose: To identify dose-volume parameters related to late rectal bleeding after radiotherapy for prostate cancer., Materials and Methods: Clinical complication data from a randomized trial were collected and linked to the individual dose-volume data. In this trial, patients with prostate cancer were treated with either conventional (with rectangular fields) or three-dimensional conformal radiotherapy to a dose of 66 Gy. Patient complaints, including rectal blood loss, were collected for 199 patients, using questionnaires. Absolute and relative dose-volume histograms (DVHs) of the rectal wall (with and without the anal region) were calculated with and without rectal filling. A proportional hazard regression (PHR) model was applied to estimate the probability of any rectal blood loss within 3 years, as a function of several DVH parameters. In a multivariable analysis, dose-volume parameters were tested together with patient- and treatment-related parameters (age, smoking, diabetes, cardiovascular disease, tumor stage, neo-adjuvant androgen deprivation, conformal vs. conventional and rectal bleeding during treatment)., Results: The estimated incidence of any and moderate/severe rectal bleeding at 3 years was 33% and 8%, respectively. Differences between the conventional and conformal technique were small and not significant. The analysis of relative DVHs of the rectal wall (with and without the anal region), showed significant (p < 0.01) relations between the irradiated volume and the probability of rectal blood loss within 3 years for dose levels between 25 Gy and 60 Gy. This relationship was shown in subgroups defined by dose-volume cutoff points as well as in the PHR model, in which a continuously rising risk was seen with increasing volumes. For absolute DVHs and DVHs of the rectum including filling, less or no significant results were observed. The most significant volume-effect relation (p = 0.002) was found at 60 Gy for the rectum wall excluding the anal region. The probability of rectal bleeding increased from 10% to 63% when the irradiated rectum volume at 60 Gy increased from 25% to 100%. Other factors. including age, smoking, diabetes, cardiovascular disease, tumor stage, neo-adjuvant androgen deprivation, conformal vs. conventional, rectal bleeding during treatment, rectum length. and whole rectum volume. did not have a significant effect in the multivariable analysis. When controlling for the volumes at 60 Gy, the volumes at lower dose levels (25-55 Gy) were no longer significant (p = 0.5)., Conclusions: For any rectal bleeding within 3 years, an overall incidence of 33% was observed for patients treated to 66 Gy. For this endpoint, a volume-effect relation was found for DVH parameters of the relative rectal wall volume. This relationship appeared to be most significant for the rectum without the anal region and for the higher dose levels (50-60 Gy).

Background and Purpose: To develop a method to project surface elements of a bent tubular organ, e.g. the rectum, in order to create a two-dimensional (2D) map and to use this method to quantify on a local scale shape and position variations of the rectum., Patients and Methods: For this study we used data of 19 patients, who each received a planning CT scan and 9-13 repeat CT scans that were considered representative for the radiotherapy course. We combined maps from multiple CT scans of the same patient to quantify local rectal wall displacements. To make a map we first computed a central axis through the rectum and divided it into segments of equal length assuming that the length of these segments was invariant under rectum shape and position changes. Next, we constructed for each segment a planar cross section through the rectum, which was oriented orthogonally to that segment. The amount of rectal wall tissue was assumed to be constant in all orthogonal cross sections throughout the entire rectum. We unfolded the cross-sected rectal wall at the dorsal side and projected either the associated dose or the coordinates onto the map., Results: The largest variation in the position of the rectal wall during the treatment course occurred at the upper anterior, left and right side (1 SD=5-7 mm). Near the anus the variation was <3 mm (1 SD) and at the posterior side of the rectum <4 mm (1 SD). The anterior-posterior (AP) and left-right displacements between the rectum in the planning CT scan and the mean rectum shape during the treatment were localized between 40 and 80% of the central axis. At the upper anterior, left, and right side the displacements were 5-8 mm (1 SD). These rectal wall displacements correlated with the rectum volume in the planning CT scan. At the upper anterior side the correlation coefficient between the AP displacements and the planning rectum volume was 0.85., Conclusions: We quantified variations in rectum shape and in dose in the rectal wall. The systematic error in rectal wall position was found to be larger than the random shape and position variations. We successfully developed a method to virtually unfold a rectum and to project the dose onto a 2D map. The spatial information of the dose distribution can be used in the analysis of rectum complications.

Purpose: To develop a model that predicts possible rectum configurations that can occur during radiotherapy of prostate cancer on the basis of a planning CT scan and patient group data., Materials and Methods: We used a stochastic shape description model with a limited number of parameters (area, area difference, and curvature) on a slice-by-slice basis to simulate rectum motion. The probability distributions of the chosen parameters were obtained from a group of 9 reference patients, who each received 15-17 repeat CT scans. We used a Monte Carlo technique to generate different rectum configurations from the probability distributions. We verified the model by comparing dose-wall histograms (DWHs) of the originally delineated rectal contours and simulated rectums for a three-field treatment technique with a prescription dose of 78 Gy. The 15-17 sets of rectal contours of each patient are regarded as the golden standard and provide a good estimate of the actual dose received during the treatment. We determined the equivalent uniform dose (EUD) for a quantitative comparison between the actual dose, the dose predicted on the basis of the simulations, and the dose predicted on the basis of a single planning CT scan., Results: The simulated rectum configurations yield a better estimate of the actual dose in the rectal wall than the rectum in the planning CT scan alone. The differences between the EUD based on the planning CT scan and the actual EUD ranged between -1.1 Gy and 2.1 Gy, with respect to a mean actual EUD of 69.8 Gy. This range is smaller for the EUD based on the simulated rectums, namely -0.4 Gy to 0.6 Gy. Furthermore, the simulation generates a set of rectum configurations that provides an estimate of the variation in DWHs during the course of the treatment. This estimate can be used in addition to the DWH of the planning CT scan in the analysis of gastrointestinal toxicity., Conclusions: To simulate rectum shapes, we have developed a model that can be used in addition to the information available in the planning CT scan in the analysis of the received dose to the rectal wall during radiotherapy of prostate cancer.