Purpose: Bi-objective treatment planning for high-dose-rate prostate brachytherapy is a novel treatment planning method with two separate objectives that represent target coverage and organ-at-risk sparing. In this study, we investigated the feasibility and plan quality of this method by means of a retrospective observer study. Methods and Materials: Current planning sessions were recorded to configure a bi-objective optimization model and to assess its applicability to our clinical practice. Optimization software, GOMEA, was then used to automatically generate a large set of plans with different trade-offs in the two objectives for each of 18 patients treated with high-dose-rate prostate brachytherapy. From this set, five plans per patient were selected for comparison to the clinical plan in terms of satisfaction of planning criteria and in a retrospective observer study. Three brachytherapists were asked to evaluate the blinded plans and select the preferred one. Results: Recordings demonstrated applicability of the bi-objective optimization model to our clinical practice. For 14/18 patients, GOMEA plans satisfied all planning criteria, compared with 4/18 clinical plans. In the observer study, in 53/54 cases, a GOMEA plan was preferred over the clinical plan. When asked for consensus among observers, this ratio was 17/18 patients. Observers highly appreciated the insight gained from comparing multiple plans with different trade-offs simultaneously. Conclusions: The bi-objective optimization model adapted well to our clinical practice. GOMEA plans were considered equal or superior to the clinical plans. In addition, presenting multiple high-quality plans provided novel insight into patient-specific trade-offs.

Purpose: To investigate the variation in computed dose-volume (DV) indices for high-dose-rate (HDR) prostate brachytherapy that can result from typical differences in computation settings in treatment planning systems (TPSs). Methods: Five factors were taken into account: number of dose-calculation points, radioactive source description, interpolation between delineated contours, intersections between delineated organ contours, and organ shape at the top and bottom contour using either full or partial slice thickness. Using in-house developed software, the DV indices of the treatment plans of 26 patients were calculated with different settings, and compared to a baseline setting that closely followed the default settings of the TPS used in our medical center. Studied organs were prostate and seminal vesicles, denoted as targets, and bladder, rectum, and urethra, denoted as organs at risk (OARs), which were delineated on MRI scans with a 3.3 mm slice thickness. Results: When sampling a fixed number of points in each organ, in order to achieve a width of the 95% confidence interval over all patients of the DV indices of 1% or less, only 32,000 points had to be sampled per target, but 256,000 points had to be sampled per OAR. For the remaining factors, DV indices changed up to 0.4% for rectum, 1.3% for urethra, and 2.6% for prostate. DV indices of the bladder changed especially if the high-dose-region was (partly) located at the most caudal contour, up to 8.5%, and DV indices of the vesicles changed especially if there were few delineated contours, up to 9.8%, both due to the use of full slice thickness for the top and bottom contour. Conclusions: The values of DV indices used in prostate HDR brachytherapy treatment planning are influenced by the computation settings in a TPS, especially at the most caudal part of the bladder, as well as in the seminal vesicles.