ORIGINAL RESEARCH article
Influence of Urethra Sparing on Tumor Control Probability and Normal Tissue Complication Probability in Focal Dose Escalated Hypofractionated Radiotherapy: A Planning Study Based on Histopathology Reference
- 1Department of Radiation Oncology, Medical Center – University of Freiburg, Faculty of Medicine, University of Freiburg, Freiburg, Germany
- 2German Cancer Consortium (DKTK). Partner Site Freiburg, Freiburg, Germany
- 3Berta-Ottenstein-Programme, Faculty of Medicine, University of Freiburg, Freiburg, Germany
- 4Division of Medical Physics, Department of Radiation Oncology, Medical Center - University of Freiburg, Faculty of Medicine, University of Freiburg, Freiburg, Germany
- 5Department of Urology, Medical Center – University of Freiburg, Faculty of Medicine, University of Freiburg, Freiburg, Germany
- 6Institute for Surgical Pathology, Medical Center – University of Freiburg, Faculty of Medicine, University of Freiburg, Freiburg, Germany
- 7Department of Nuclear Medicine, Medical Center – University of Freiburg, Faculty of Medicine, University of Freiburg, Freiburg, Germany
- 8Department of Radiology, Medical Center – University of Freiburg, Faculty of Medicine, University of Freiburg, Freiburg, Germany
Purpose: Multiparametric magnetic resonance tomography (mpMRI) and prostate specific membrane antigen positron emission tomography (PSMA-PET/CT) are used to guide focal radiotherapy (RT) dose escalation concepts. Besides improvements of treatment effectiveness, maintenance of a good quality of life is essential. Therefore, this planning study investigates whether urethral sparing in moderately hypofractionated RT with focal RT dose escalation influences tumour control probability (TCP) and normal tissue complication probability (NTCP).
Patients and Methods: 10 patients with primary prostate cancer (PCa), who underwent 68Ga PSMA-PET/CT and mpMRI followed by radical prostatectomy were enrolled. Intraprostatic tumour volumes (gross tumor volume, GTV) based on both imaging techniques (GTV-MRI and -PET) were contoured manually using validated contouring techniques and GTV-Union was created by summing both. For each patient three IMRT plans were generated with 60 Gy to the whole prostate and a simultaneous integrated boost up to 70 Gy to GTV-Union in 20 fractions by (Plan 1) not respecting and (Plan 2) respecting dose constraints for urethra as well as (Plan 3) respecting dose constraints for planning organ at risk volume for urethra (PRV = urethra + 2mm expansion). NTCP for urethra was calculated applying a Lyman-Kutcher-Burman model. TCP-Histo was calculated based on PCa distribution in co-registered histology (GTV-Histo). Complication free tumour control probability (P+) was calculated. Furthermore, the intrafractional movement was considered.
Results: Median overlap of GTV-Union and PRV-Urethra was 1.6% (IQR 0-7%). Median minimum distance of GTV-Histo to urethra was 3.6 mm (IQR 2 – 7 mm) and of GTV-Union to urethra was 1.8 mm (IQR 0.0 – 5.0 mm). The respective prescription doses and dose constraints were reached in all plans. Urethra-sparing in Plans 2 and 3 reached significantly lower NTCP-Urethra (p = 0.002) without significantly affecting TCP-GTV-Histo (p = p > 0.28), NTCP-Bladder (p > 0.85) or NTCP-Rectum (p = 0.85), resulting in better P+ (p = 0.006). Simulation of intrafractional movement yielded even higher P+ values for Plans 2 and 3 compared to Plan 1.
Conclusion: Urethral sparing may increase the therapeutic ratio and should be implemented in focal RT dose escalation concepts.
Radiotherapy (RT) of primary Prostate cancer (PCa) is currently experiencing an individualization, utilizing modern imaging techniques for staging and definition of intraprostatic gross tumor volume (GTV). Since an increase in RT dose improves tumor control rates (1), concepts of focal dose escalation have developed to deliver higher doses to the tumor and thereby improving rates of biochemical recurrence (2, 3) without risking higher toxicities by respecting OAR restrictions. Recently the long-term result of the phase III FLAME trial demonstrated that mpMRI-defined focal dose escalation significantly improves biochemical disease free survival (4). Earlier publications from this trial demonstrated the feasibility and reported no significant increase in acute and late toxicities (5). These results are encouraging, that unfavorable intermediate- and high-risk PCa patients, who’s proportion is on the rise (6), benefit from these advanced treatments. Besides multiparametric magnetic resonance tomography (mpMRI) being the gold standard for diagnostics in PCa (7), prostate specific membrane antigen positron emission tomography (PSMA-PET) has emerged as a diagnostic tool of high quality (8–13). Recently, the superiority of PSMA-PET for initial staging compared to conventional imaging was prospectively proved, which led to therapy management change in 28% of cases (14). Regarding depiction of the intraprostatic GTV PSMA-PET/CT reveals GTVs more concordant with biopsy reference (9), whereas mpMRI underestimates the true tumor and misses significant tumour lesions (15–17). Previously conducted planning studies from our group and from Goodman et al. demonstrated that despite putative limitations for focal therapy approaches due to larger volumes, boosting of PSMA-PET/CT delineated GTVs is technically feasible (18–20). Bettermann et al. and Eiber et al. could clearly demonstrate that the combined use of mpMRI and PSMA-PET (GTV-Union) significantly improved sensitivity (9, 11). A planning study by Zamboglou et al. revealed significantly increased tissue control probabilities (TCP) for GTV-Union based focal dose escalation compared to GTV-PET or GTV-MRI-based dose escalation (20). Prospective trials will evaluate whether these advances in imaging and diagnostic accuracy can be translated into improved clinical outcomes. A modern approach includes moderately hypofractionated RT (MHRT) to the whole prostate with simultaneously integrated dose escalation to mpMRI- and PSMA-PET/CT-defined GTVs. Although the impact of accountable structures such as bladder, bladder trigone and urethra stay vague, the urethra as a serial organ may be of particular importance in this setting. This planning study aims to investigate whether urethral sparing in MHRT with focal dose escalation delivered to mpMRI and PSMA-PET/CT defined GTVs, influences tumor control probability (TCP) and normal tissue complication probability (NTCP). NTCP was calculated based on the Lyman-Kutcher-Burman (LKB) model with parameters defined by Panettiere et al. (21), TCP was calculated based on 3D dose distribution in co-registered histopathology as standard of reference. Furthermore, the influence of intrafractional movement was assessed (22).
Material and Methods
The utilized study cohort consisted of ten (10) patients with primary PCa, who underwent 68Ga-HBED-CC-PSMA (68Ga-PSMA-PET) and mpMRI followed by radical prostatectomy. Patient characteristics are listed in Table 1. A written informed consent was obtained from each patient and the institutional review board of the Albert-Ludwigs-University of Freiburg approved the study (No.: 469/14).
PET/CT and MRI Imaging
Diagnostic images were acquired using a diagnostic setup.
PET/CT scans using the ligand 68Ga-HBED-CC-PSMA (23) were performed in 9 patients with a 64-slice GEMINI TF PET/CT and in 1 patient with a Vereos PET/CT (both Philips Healthcare, USA). The imaging systems were cross-calibrated to ensure the comparability of the quantitative measurements and both scanners fulfilled the requirements indicated in the European Association of Nuclear Medicine (EANM) imaging guidelines and obtained EANM Research Ltd. (EARL) accreditation during acquisition. The spatial resolution in the transverse direction near the centre is 4.8 mm for GEMINI TF (24) and 4.2mm for Vereos (25) Patients underwent the whole-body PET scan starting 1 h after injection andwere asked to urinate prior PET imaging. The uptake of 68Ga-PSMA-HBED-CC was quantified by standardized uptake values (SUV). A detailed description of the used 68Ga-HBED-CC-PSMA PET/CT imaging protocol is described in (26).
MR images were acquired on a 3 Tesla system (5 patients on TrioTim, 1 patient on Magnetom Vida, 1 patient on Skyra, all Siemens, Germany) and on a 1.5 Tesla system (3 patients on Aera, Siemens, Germany). The MR imaging systems were equipped with a surface phased array (Body Matrix) in combination with an integrated spine array coil. No endo-rectal coil was used. Not additional cross-calibration was performed. Essentially, T2-weighted fast spin echo (T2W-TSE) images, diffusion weighted images (DWI) and dynamic contrast-enhanced (DCE) perfusion images were acquired. Apparent diffusion coefficient (ADC) maps were calculated from the DWIs using information from all measured b-values. ADC maps were generated with a monoexponential model as implemented in syngo.via (syngo.via ADC & b-value tool, Siemens Healthcare, Germany). Extrapolated high b-value images (b = 1400 s/mm2) were calculated with syngo.via using information from all measured b-values. These extrapolated images were considered the high b-value DWIs for prostate MRI reading according to PI-RADS v2.0 (27). MR protocols were heterogeneous in terms of slice thickness, gap between slices and b-values. A detailed description of the used T2w, DWI and DCE MRI imaging protocol can be found in (28).
Intraprostatic Tumour Mass
GTV-PET was contoured manually using a validated scaling of SUVmin-max: 0-5 (29) within the prostate using Eclipse™ Treatment Planning System (Varian, USA). GTV-MRI was contoured manually based on MRI T2-w and ADC images, applying imaging criteria PI-RADSv2.0 and considering lesions with a PI-RADS score of ≥ 3 as relevant (27). Final GTVs were the respective consensus contour between two readers with >4 years experience in PET and MRI interpretation. Subsequently careful manual co-registration of in-vivo CT and in-vivo MRI was performed to transfer GTV-MRI to the corresponding in-vivo-CT image and to create GTV-Union composed of the sum of GTV-PET and GTV-MRI. GTV-Union was used based on the benefit in terms of higher sensitivity and complementary information of both techniques (9, 13).
Organs at Risk (OAR), Clinical Target Volume (CTV) and Planning Target Volume (PTV)
Bladder, rectum, femoral head as OAR were contoured based on the planning-CT scan according to RTOG guidelines (30). Urethra was contoured based on the co-registered MRI. Planning organ at risk volume (PRV)-urethra was created from applying 2 mm isotropic extension of urethra according to the hypo-Flame trial (31).
The CTV was created by following the ESTRO-ACROP guidelines (32). CTV1 was defined as the prostate including extracapsular PCa + 3mm isotropic extension (excluding rectum and bladder). CTV_SV was defined as the proximal 1.4 cm and 2.2 cm of the seminal vesicle (SV) in unfavorable intermediate risk and high-risk patients accordingly. In case of tumor infiltration of the SV the respective regions were included in CTV-SV. CTV2 was defined as the prostate and the base of the SV including parts of the SV with visible tumor burden. PTV 1 was created from isotropic 6 mm-extension of CTV1 and 8 mm of CTV_SV, followed by merging both volumes. PTV2 was created from isotropic 6 mm-expansion of CTV2. PTV3 was created from isotropic 2mm-extension of GTV-Union and consequent remove of existing overlaps with organ at risks (OAR) contours. For analysis purposes three different PTV3 were generated: PTV3_1 was defined as the GTV-Union isotropically expanded by 2mm. PTV3_2 was created from the subtraction of urethra from PTV3_1 and finally PTV3_3 was created from the subtraction of PRV-Urethra from PTV_1. See Figure 1 for illustration of volumes.
Figure 1 Shows GTV-MRI (blue), GTV-PET (green), urethra (orange) PRV-urethra (yellow), prostate (red), PTV3_3 (boost volume minus PRV-urethra, purple) and PTV1 (pink).
PCa lesions in whole mount histopathology were used as standard of reference as previously conducted by our group (9, 33). After fixation, the resected prostate was fixed in a customized localizer with a 4 mm grid and an ex-vivo CT scan (16-channel Brilliance Big Bore, Phillips, Germany) was performed. Subsequently, whole-mount step sections were cut every 4 mm using an in–house cutting device to guarantee equal cutting angles between histological and corresponding ex-vivo CT slices. Following paraffin embedding, specimens were cut using a Leica microtome. Haematoxylin and eosin staining were performed following routine protocols. A board-certified experienced pathologist marked PCa lesions. Subsequently, histopathological information was digitalized via intermediate registration to ex-vivo CT using MITK software (MITK Workbench 2015.5.2). Automatic interpolation was performed to create GTV-Histo (GTV based on histopathology). Images were transferred to Eclipse™ Treatment Planning System v.15.6 (Varian Medical Systems, USA). Ex-vivo CT and in vivo-CT (from PSMA-PET/CT scans) were carefully manual co-registered allowing non-rigid deformation and considering the 4mm grid and anatomical landmarks such as urethra and cyst and prostate capsule in particular. Hence, this registration workflow takes into account non-linear shrinkage and distortion of the prostate gland after resection.
Distances to Urethra
Minimum distance of GTV-Histo to urethra was evaluated on each hematoxylin and eosin stained (H&E) slice of the respective patient. Accordingly, minimum distance of GTV-Union to urethra was evaluated on the corresponding CT-slice on the in vivo CT.
IMRT plans were created in Eclipse™ Treatment Planning System v15.1 (Varian, USA) with a calculation grid size of 1.5 mm. Dose prescription protocols were the following: PTV1 45 Gy in 15 fractions and PTV2 15 Gy in 5 fractions, resulting in 60 Gy for PTV2. A simultaneous integrated boost (SIB) up to 70 Gy for PTV3 for all 20 fractions was prescribed. Adapted from the DELINEATE trial (34) and based on findings from Martinez et al. (1), our dose concept aimed for boost doses near 100 Gy (EQD2, a/b=1.6). For PTV2 D98% was ≥ 58.8 Gy and D2% ≤ 70 Gy, for PTV3 D98% was ≥ 68.6 Gy and D2% ≤ 71.4 Gy. Three different plans were created using three different boost volumes for the simultaneous integrated boost (SIB): The SIB volumes were PTV3_1, PTV3_2 and PTV3_3 for plan 1, 2 and 3 respectively. Dose constraints for organs at risk were considered according to CHHiP-, FLAME- and DELINEATE-trial (5, 34–36). Dose constraints for Urethra and PRV-Urethra were 62.4 Gy for D2%. Details of RT planning prescription doses and OAR constraints can be found in Supplementary Material 1.
To evaluate the impact of urethral sparing three different IMRT plans were calculated: (i) Plan 1 without any dose constraints for urethra, (ii) Plan 2 considered the D2% dose constraint for urethra and (iii) Plan 3 considered the D2% dose constraints for PRV-Urethra.
TCP and NTCP Modeling
Structure sets and calculated 3D-dose matrices of the radiotherapy plans were exported as DICOM files. Furthermore, using a Varian ESAPI script (https://varianapis.github.io/), dose matrix voxels for each structure were exported (https://github.com/isachpaz/ESAPICommander). TCP was calculated based on the linear quadratic (LQ) Poisson model (37–41):
Where ρcl is the homogeneous clonogenic cell density (# cells/cm3) in the tumor of volume V. EQD0 is the equi-effective dose for 0Gy fractionation given by Eq. 2, and α is the coefficient of LQ-model defining the linear-term of cell killing.
d is the dose per fraction, and D is the total dose delivered in N fractions, D=N*d, where α/β is the ratio of linear to quadratic cell killing probability according to the LQ-model.
In the present study, the tumor cell density was set ρcl=2.8*108 cells/cm (42–44). α/β value of 1.6 Gy was assumed, based on the recent meta-analysis results by Vogelius et al., which included studies with mildly- and ultrahypofractionated radiotherapy (45). To account for diversity of published α/β values we performed a robustness analysis for TCPGTV-Histo with three different parameter sets encompassing the range for α/β described by Vogelius et al. (45). α was each time fitted (Table 2), so that 70% TCPGTV-Histo would be reached in our patient cohort with a conventional dose of 60 Gy in 3 Gy fractions, as we have described in our previous publication (19).
NTCP for bladder and rectum (NTCPBladder, NTCPRectum) were calculated based on the relative seriality model as described by Bostel et al. (46). For bladder a D50 of 80.0 Gy as EQD2 for symptomatic contracture and volume loss, a relative seriality parameter value s of 1.3 and γ = 2.59 were used (47). For rectum a D50 of 80.0 Gy as EQD2 for severe proctitis/necrosis/stenosis/fistula (2, 47–50), s = 0.75 and γ = 1.79 were considered (47). An α/β value of 3.0 Gy for bladder and rectum was assumed (34). For NTCPUrethra the Lyman-Kutcher-Burman (LKB) model was applied for the endpoint urethral stricture as published by Panitierri et al. (21): D50 = 116.7 Gy, m = 0.23, n = 0.3, and α/β = 5.0 Gy. We additionally performed NTCPUrethra calculations for the 68% confident intervals (CI) with a step of 1.0 Gy for D50 and a step of 0.01 for m. In total 364 combinations of D50 and m were evaluated.
Complication Free Tumour Control Probability P+
In order to account for the injuries or risk for complications to each of the healthy organs (OARs) involved in a given clinical case, the following expression is usually applied for the total probability of injury PI:
where NTCPj is the probability of injuring the normal tissue (OAR), wj is a weighting factor expressing the relative clinical importance of each endpoint, and NOARs is the total number of healthy organs involved in the clinical case. The effectiveness of a given dose distribution can be evaluated by the comparison of its advantages in terms of tumour control (benefit B) against its disadvantages considering normal tissues complications (injury I). The probability of complication free tumour control P+, is defined as
where PB is the probability of getting benefit from the treatment (tumour control, Eq. 1) and PI is the probability of causing injury to normal tissues (Eq. 3). For the case of complete independency of response of tumor and OARs, P+ becomes:
P+ is an overall parameter for evaluation of complex dose distributions and treatment localisations and is suggested to support decision on treatment plan selection and treatment adaptation (46, 51–56).
As previously performed by our group, TCP and NTCP calculations were calculated with and without movement (22). This was achieved by changing the relative positioning between structure matrix and dose matrix implementing Gaussian filtering of the dose matrix. Based on results of Langen et al. (57) the standard deviation of a three-dimensional Gaussian kernel, was set to 0.92 mm, 1.59 mm and 1.54 mm for left-right, anterior-posterior and cranio-caudal, respectively.
The Sörensen-dice coefficient was calculated for spatial overlap of GTV-Histo with GTV-Union, GTV-PET and GTV-MRI and for spatial overlap of PTV3_1, PTV3_2, PTV3_3 and GTV-Histo.
Statistical analysis of volumes was performed with GraphPad Prism v8.4.2 (GraphPad Sofware, USA). Data normality was tested using the Shapiro-Wilk test. For not normally distributed variables, Friedman test and uncorrected Dunn’s test was used for comparison of more than two variables and two-sided Wilcoxon matched-pairs signed rank test was used for comparison of two variables (both at a significance level of 0.05). For normally distributed variables, repeated measures one-way ANOVA with the Geisser-Greenhouse correction and Fisher’s LSD was used for comparison of more than two variables and two-sided paired t test was used for comparison of two variables (both at a significance level of 0.05). For statistical analysis of unpaired and not normally distributed data (minimum distance to urethra on H&E slices and CT images) Mann-Whitney test at a significance level of 0.05 was used.
Exploratory statistical analysis of TCPs, NTCPs and dosimetric analysis was performed with R (version 3.6.2) (58). Wilcoxon matched pairs signed-rank test was used with a significance level of 0.05.
Volumes and Distance of GTVs to Urethra
Median volume for GTV-Histo was 4.5 ml (IQR 1.8 – 6.9ml) and for GTV-Union 5.7 ml (IQR 2.9 – 13.3 ml). Median intersection volume of PRV-Urethra with GTV-Histo was 0.05 ml (IQR 0.00 - 0.25 ml) and with GTV-Union 0.1 ml (IQR 0.00 – 0.88 ml) respectively. Expressed in percentage of the PRV-urethra volume, intersection of GTV-Union with PRV-Urethra was median 1.6% (IQR 0.0 – 6.5%) and maximum 8.5% in patient 10. Please see Supplementary Table 1 for details.
Median volumes for PTV3_1 was 13.5 ml (IQR 7.0 - 22.6 ml), for PTV3_2 13.2 ml (IQR 6.9 – 22.0 ml) and for PTV3_3 12.8 ml (IQR 6.6 – 20.6 ml), respectively. PTV3_3 was not statistically significantly smaller than PTV3_2 (p = 0.053) but significantly smaller than PTV3_1 (p = 0.031) (Supplementary Table 2).
The median intersection volume of GTV-Histo with PTV3_1, PTV3_2 and PTV3_3 was 2.7 ml (IQR 1.5 – 6.3), 2.7 ml (IQR 1.5 – 6.2) and 2.7 ml (1.4 – 5.9), respectively. There was no statistically significant difference between the DSCs for GTV-Histo and the three PTVs (p > 0.96) (Supplementary Tables 2 and 3).
Median coverage of GTV-Histo by GTV-Union, GTV-PET and GTV-MRI was 79% (IQR 55 – 97%), 76% (IQR 37 – 83%) and 53% (IQR 13 – 74%). Coverage by GTV-Union was significantly higher than by GTV-PET (p = 0.014) and GTV-MRI (p = 0.004), whereas there was no significant difference between GTV-PET and GTV-MRI (p = 0.058). Median coverage of GTV Histo by PTV3_1, PTV3_2 and PTV3_3 was 90% (IQR 70 – 92%), 89% (IQR 70 – 91%) and 85% (IQR 65 – 88%). Coverage by PTV3_3 was significantly lower than by PTV3_1 (p=0.016) (Supplementary Tables 3 and 4).
In 3 patients contact between GTV-Histo and urethra could be observed on H&E slices. In 6 patients contact between GTV-Union and urethra could be observed on in-vivo CT slices. Discrepancies between patients with detected contact on slices but without intersection volumes were manually verified. In all cases intersection volumes were present but too small to be quantified in Eclipse™ Treatment Planning System. The median minimum distance of GTV-Histo to urethra on each slice was 3.6 mm (IQR 2.2 – 7.3 mm) and median minimum distance of GTV-Union to urethra was 1.8 mm (IQR 0.0 – 5.0 mm). Distance of GTV-Union to urethra was statistically significantly smaller (p = 0.02). Median minimum distance of GTV-Histo to urethra per patient was 1.9 mm (IQR 0.0 – 3.6 mm) and median minimum distance of GTV-Union to urethra per patient was 0.0 mm (IQR 0.0 – 1.5mm). Again, distance of GTV-Union to urethra was statistically significantly smaller (p = 0.02).
Doses Distribution in Target Volumes
Median D98%, D50% and D2% doses for, PTV3_1-3 (boost volume), GTV-Histo, urethra and PRV-Urethra for plan 1-3 are shown in Table 3 with and without consideration of the intrafractionary movement, respectively. Without consideration of intrafractionary movement following doses were statistically significant different: For PTV3 D98% of plans 2 and 3 were significantly smaller than for plan 1, whereas D50% showed no significant difference between the three plans. D2% was significantly higher in plans 2 and 3 than in plan 1. For GTV-Histo, only D98% was significantly lower for plan 2 and 3 compared to plan 1. With consideration of intrafractionary movement following doses were statistically significantly different:
Table 3 Dose volume parameter values for different volumes without and with consideration of prostate intrafractional movement.
For PTV3, D98% and D50% of plan 3 were slightly but significantly smaller compared to plan 1, whereas dose parameters of plan 2 showed no statistical significance to dose parameters of plan 3. D2% was significantly higher in plans 2 and 3 than in plan 1. For GTV-Histo D98% was significantly smaller and D2% significantly higher in plans 2 and 3 compared to plan 1.
Doses for urethra and PRV-urethra were significantly lower in both plan 2 and plan 3 compared to plan 1 in all cases. Furthermore, all doses were significantly lower in plan 3 compared to plan 2 except for D98% with movement. For details about p-values see Supplementary Table 6. Figure 2 shows cumulative dose-volume-histograms for boost volumes, urethra, bladder and rectum without and with movement.
Figure 2 Shows cumulative dose-volume-histograms for boost volumes (PTV3_1-3), urethra, bladder and rectum without movement (A) and with movement (B), respectively.
All plans complied with the constraints for bladder and rectum.
Without consideration of intrafractional movement, in plan 1 (no dose constraints for urethra considered in optimization) constraints for D2% for urethra were not reached in the majority of the planed cases, 8 out of 10. In plan 2 (respecting dose constraints for urethra), constraints for D2% for PRV-urethra were not reached in 7 patients.
When intrafractional movement is considered, in plan 1 constraints for D2% for urethra were not reached again in in 8 patients. In plan 2, constraints for D2% for urethra were not reached in 4 patients and D2% for PRV-urethra in 7 patients.
In plan 3 urethra- and PRV-Urethra constraints were reached in all patients without and with movement consideration.
TCP/NTCP/P+ Without Intrafractional Movement
Please see Table 4 for median P+, TCPGTV-Histo, NTCPUrethra, NTCPBladder and NTCPRectum as well as p-values. Urethra-sparing resulted in significantly lower NTCPUrethra without significantly affecting TCPGTV-Histo or NTCPBladder and NTCPRectum. Consequently, P+ was statically significantly better for plans respecting urethral sparing. Radiobiological modeling was also performed by assuming α/β values of 1.2 Gy and 2.7 Gy for tumor tissue (see Supplementary Tables 7 and 8). Summarized P+ shows the same behavior for α/β = 1.2 Gy, whereas for α/β = 2.7 Gy no significant differences between all three plans could be observed. For the calculation of total probability of injury PI, required for the complication free tumour control P+ (Eq. 5), all three weighting factors wj in Eq. 3 for the OARs are set to 1.0 (equal clinical importance). Analysis on patient level revealed, that P+ was higher in plan 2 and plan 3 compared to plan 1 in all patients.
Table 4 Median P+, TCP- and NTCP-values for plans 1-3, as well as p-values for comparison of plan 1 vs 2 and 3, respectively, not considering intrafractional movement with α/β 1.6 Gy for tumor tissue and 3 Gy for bladder and rectum.
TCP/NTCP/P+ With Intrafractional Movement
Implementation of intrafractional movement into the model yielded in even slightly higher P+ for both urethral sparring plans (Table 5). For median P+, TCPGTV-Histo, NTCPUrethra, NTCPBladder and NTCPRectum as well as p-values considering intrafractional movement see Table 5. Radiobiological modeling was performed by assuming α/β values of 1.2 Gy and 2.7 Gy for tumor tissue (see Supplementary Tables 9 and 10), and summarized P+ shows again the same behavior for α/β = 1.2 Gy, whereas for α/β = 2.7 Gy no significant differences between all three plans could be observed. Analysis on patient level revealed, that P+ was higher in plan 2 and plan 3 compared to plan 1 in all patients.
Table 5 Median P+, TCP- and NTCP-values for plans 1-3, as well as p-values for comparison of plan 1 vs 2 and 3, respectively, considering intrafractional movement with α/β 1.6 Gy for tumor tissue and 3 Gy for bladder and rectum.
Re-run of NTCPUrethra, calculation in order to consider uncertainties, showed no deviation from initial outcomes.
In the context of focal escalation, the results of our planning study demonstrate that boosting of PSMA-PET/CT and mpMRI defined GTVs using MHRT is technically feasible and prescription doses as well as dose constraints are achieved even when considering organ movement. Furthermore, urethral sparing achieves significantly lower NTCPs for urethral toxicities without affecting TCPs and NTCPs for bladder and rectum, consequently results in a better therapeutic ratio in terms of P+ and should be implemented in focal RT dose escalation concepts. We discuss the different aspects in the following sections in detail.
Urethra sparing is performed in SBRT and brachytherapy, since higher urethral doses are associated with higher GU toxicities (59, 60). The recently published toxicity reports of the hypoFlame trial suggest that prioritization of OAR constraints yields acceptable toxicities for focal dose escalation using SBRT. At the time of publication, the only published toxicity reports of trials investigating moderately-hyofractionated dose escalation and urethral sparing was the DELINEATE trial, which used MRI-defined boost volumes with dose escalation up to 67 Gy and showed slightly higher, but comparable acute and late GI and GU toxicity rates to dose escalation with conventional fractionation (34). In the cohort receiving focal dose escalated MHRT cumulative grade 2 or worse GU and GI toxicities after 3 years were 22.1% and 14.0%, respectively. The dose regimen chosen in our study utilized a higher prescription dose of 70 Gy for the boost volume defined by validated contouring approaches for GTV definition of PET and mpMRI imaging modalities (29, 61). Not surprisingly, volumes for GTV-Union (median 5.7 ml, IQR 2.9 – 13.3 ml) were significantly larger than GTV-Histo (median 4.5 ml, IQR 1.8 – 6.9 ml, p = 0.01). Consequently, prioritization of maintaining standard toxicity rates is pivotal when boosting larger volumes and therefore we conducted this planning study to evaluate the effect of urethral sparing on NTCPs and TCP to evaluate its potential in MHRT with focal dose escalation with this novel boost volume definition.
Constraints and prescription doses were achieved for all patients as intended in the respective plans. This also applies when implementing organ movement into the plan evaluation. These results suggest that putative negative consequences in terms of under- or overdosing were not relevant and consequently organ movement did not affect the highly conformal IMRT plans. Furthermore, we confirmed the feasibility of the applied dose and constraint prescription.
To evaluate the impact of urethra-sparing we chose as endpoint for NTCPUrethra stricture requiring urethrectomy within 4 years based on the LKB model by Panettiere et al. (21). Based on recently published α/β values from Vogelius et al. we used an α/β of 1.6 Gy to calculate the TCP. Considering the published range for α/β values, we performed the same analysis with an α/β 1.2 Gy and 2.7 Gy. Application of 1.2 Gy yielded similar results with a significant better P+ for plan 3 (0.924 and 0.9285 without and with movement, respectively, p = <0.01), whereas application of 2.7 Gy resulted in no significant improvement of P+ (see Supplementary Tables 7–10). However, the mentioned meta-analysis suggests that α/β of 2.7 Gy is likely to be too high, particularly in a setting of hypofractionation. Therefore, we refer on the results derived from α/β of 1.6 in the following. Urethral sparing in IMRT planning significantly reduced the median NTCPUrethra from 7.2% up to 4.2% (p=0.002) with only minimal and statically not significant reduction of median TCP from 99.7% up to 99.2% (p = 0.105). Noteworthy, NTCPBladder and NTCPRectum were not affected by urethral sparing, precluding the possibility of improving NTCPUrethra at the costs of other toxicities. Consequently, urethral sparing resulted in significantly better P+ value (88.8% vs up to 91.9%, p = 0.006). Considering that the urethra is a serial organ the minimum distance of urethra to PCa is of particular clinical relevance. Evaluation of minimum distances of GTV-Union to urethra was significantly smaller than GTV-Histo and in 60% of patients, contact of GTV-Union with urethra could be determined, supporting the rationale of urethra sparing. The minimal impact on TCP can be attributed to the small intersection of PCa tissue and urethra. Even applying a margin to the urethra resulted in intersection of GTV-Histo and PRV-Urethra of median 0.8% and maximum 6.4%, intersection of GTV-Union with PRV-Urethra was median 1.6% and maximum 8.5%. Considering volume analysis subtraction of PRV-Urethra from PTV3 resulted in slightly, but significantly smaller PTV-volumes, as well as slightly but significantly lower coverage of GTV-Histo. Nevertheless, coverage was still very high (85%) and even though doses partly showed significant differences between IMRT plans, the mentioned differences had no significant consequences on TCPGTV-Histo. Additionally, coverage of GTV-Histo by GTV-Union was statistically significantly higher than by GTV-PET or GTV-MRI, supporting the rationale to implement both imaging modalities in boost volume definition. Overall results of volumetric analysis and TCP/NTCP calculation suggest that boosting of GTV-Union is compatible with sufficient urethra sparing in most cases. A study by Leibovich et al., which found that the mean distance from the urethra to the nearest cancer was 3mm (62). In our study median distance of GTV-Histo to urethra was comparable with 3.6 mm. Even though contact with urethra was detectable in the majority of cases, intersection volumes were very small, supporting the estimation of little consequences of intersection of urethra and GTV.
Comparison of our results with different planning studies is hampered due to lack of data. One other study evaluated NTCPUrethra and showed extremely high NTCP-values >60% by using TD50 of 70.7 Gy (63). The results of our study are still higher than clinical reported urethral stricture rates after external beam RT (EBRT) ranging between 2-3%, nevertheless applied doses were lower and these studies did not use focal dose escalation (64–66). Considering this aspect our results represent realistic estimations and should be compared with eagerly awaited long term results of clinical trials investigating focal dose escalated EBRT.
Additionally, we simulated intrafractional organ movement in order to evaluate its consequences on IMRT plans and TCP/NTCP calculation. This implementation had slightly positive effects on P+. Furthermore, urethral sparing did still significantly reduce NTCPUrethra (from 7.2% to 4.3%) without significantly affecting TCPHisto, NTCPBladder or NTCPRectum. Remarkable NTCPBladder and NTCPRectum were even slightly better. These results suggest that intrafractional movement potentially influences positively P+ and the used margins were adequate to compensate intrafractional movements. This complies with previously reported results by Thomann et al., which demonstrated this effect in cases where boost volumes do not fully comply with GTV-Histo (22). These results are encouraging that urethral sparing might significantly reduce GU toxicities without significantly affecting tumour control, in particular since all patients in our study benefited from urethral sparing in terms of improved P+-values in plan 2 or 3 compared to plan 1. In order to evaluate the clinical benefit of this approach it should be evaluated in clinical trials. This enables to evaluate, whether specific patient subgroups don’t benefit from this approach. Likely, in patients with high tumour burden or niches with radio-resistant PCa cells (67) surrounding the urethra, sparing might not be an advantage. Whether a threshold in terms of absolute or relative volume of intersection between PRV-urethra and boost-PTV exist, from which on positive effects are reversed, should be evaluated in further studies and larger cohorts. In this context, individual radiosensitivity might be another important aspect, possibly causing a reduced tumour control with urethral sparing in patients with low radiosensitivity and a negligible impact of urethra sparing in patients with high radiosensitivity due to sufficient dose delivery. We estimated an equal radiosensitivity for all patients in our planning study. To the best of our knowledge validated surrogate parameters to determine radiosensitivity are missing and our data don’t allow do draw conclusions in this regard. Future research might enable to consider this aspect. However, urethra sparing offers another tool for individualizing radiotherapy and acknowledging patients’ preferences, for instance a high demand for safety vs. tumour control.
Furthermore, adherence to high quality of image acquisition (68), image co-registration (69) as well as accuracy of delineation (29, 70) and radiation delivery (71) is a prerequisite for implementation of this individualized radiotherapy approach. In the context of delineation, progress in diagnostics has to be considered. Regarding PSMA-PET/CT, Fluorine-18-labeled tracers like 18F-PSMA-1007 have been implemented in nuclear medicine practise. Current research shows, that 18F-PSMA-1007 shows very high sensitivities and high specificities (70, 72). Since accurate delineation of boost volumes for focal therapy approaches depends on the applied windowing (73), usage of validated contouring approaches is necessary (13). Whether usage of 18F-PSMA-1007 affects TCP calculation compared to 68Ga-PSMA should be evaluated in future studies.
We acknowledge the limitations of our study. Firstly, it should be mentioned that the NTCP for the urethra was modeled based on a previous publication of Panettieri et al. (21). The analysis was based on 258 which received EBRT and brachytherapy. Thus, without loss of the generality our analysis was based on a parameter-set, which has been modeled with combined treatment. Secondly, we enrolled a relatively small number of patients, which is a result of the elaborate co-registration pathway of the histopathologic specimens. Thirdly, the co-registration between histopathologic 3D-volumes and cross-sectional images bears risks of uncertainty due to non-linear shrinkage of the prostate after prostatectomy and co-registration mismatch susceptibilities. Consequently, coverage of GTV-Histo boost volumes might lack precision. Fourthly, used models for TCP and NTCP calculation could not be validated with the institutions own experiences, since the follow-up database of patients treated with mildly hypofractionated EBRT was not sufficient. Fifthly, co-registered images were acquired in a diagnostic setup, potentially affecting image registration and dose calculation. Therefore, the included patients, which are part of a larger cohort, were selected in terms of bowel and bladder preparation and positioning enable BRT planning. However, our experiences for image co-registration are in line with a recently published study, demonstrating no significant differences in MRI acquisition in diagnostic and radiotherapy setups (74). Furthermore, different PET/CT and MRI scanners were used. This limitation was considered by cross-calibration of the PET scanners and a reasonable and recommended slice thickness of 3 mm was acquired in all patients.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
The studies involving human participants were reviewed and approved by Institutional review board of the Albert-Ludwigs-University of Freiburg. The patients/participants provided their written informed consent to participate in this study.
SS, IS, and CZ contributed to conception and design of the study. RW and BT created IMRT plans. MB, IS, and DB implemented intrafractional movement modeling. SS and IS conducted the statistical analysis. AS was involved in surgery indication and prostatectomy. PB conducted the histopathological processing and marking of tumor lesions. JR was responsible for conduction and reporting of the PSMA-PET/CTs. MB was responsible for conduction and reporting of the mpMRIs. TS, NN, DB, and AG supervised the study. SS and IS wrote the first draft of the manuscript. CZ wrote sections of the manuscript. All authors contributed to the article and approved the submitted version.
Ilias Sachpazidis and Constantinos Zamboglou received funding from the “Dekade gegen Krebs” initiative from the BMBF. SS received funding from the EraPERMED Call 2019 (PersoRad, BMBF). Material and publication fee is funded by the EraPERMED Call 2019 (PersoRad, BMBF).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.652678/full#supplementary-material
1. Martinez AA, Gonzalez J, Ye H, Ghilezan M, Shetty S, Kernen K, et al. Dose Escalation Improves Cancer-Related Events At 10 Years for Intermediate- and High-Risk Prostate Cancer Patients Treated With Hypofractionated High-Dose-Rate Boost and External Beam Radiotherapy. Int J Radiat Oncol Biol Phys (2011) 79(2):363–70. doi: 10.1016/j.ijrobp.2009.10.035
2. Michalski JM, Moughan J, Purdy J, Bosch W, Bruner DW, Bahary JP, et al. Effect of Standard vs Dose-Escalated Radiation Therapy for Patients With Intermediate-Risk Prostate Cancer: The NRG Oncology RTOG 0126 Randomized Clinical Trial. JAMA Oncol (2018) 4(6):e180039. doi: 10.1001/jamaoncol.2018.0039
3. Morris WJ, Tyldesley S, Rodda S, Halperin R, Pai H, McKenzie M, et al. Androgen Suppression Combined With Elective Nodal and Dose Escalated Radiation Therapy (the ASCENDE-RT Trial): An Analysis of Survival Endpoints for a Randomized Trial Comparing a Low-Dose-Rate Brachytherapy Boost to a Dose-Escalated External Beam Boost for High- and Intermediate-risk Prostate Cancer. Int J Radiat Oncol Biol Phys (2017) 98(2):275–85. doi: 10.1016/j.ijrobp.2016.11.026
4. Kerkmeijer LGW, Groen VH, Pos FJ, Haustermans K, Monninkhof EM, Smeenk RJ, et al. Focal Boost to the Intraprostatic Tumor in External Beam Radiotherapy for Patients With Localized Prostate Cancer: Results From the FLAME Randomized Phase III Trial. J Clin Oncol (2021) 39(7):787–96. doi: 10.1200/JCO.20.02873
5. Lips IM, van der Heide UA, Haustermans K, van Lin EN, Pos F, Franken SP, et al. Single Blind Randomized Phase III Trial to Investigate the Benefit of a Focal Lesion Ablative Microboost in Prostate Cancer (Flame-Trial): Study Protocol for a Randomized Controlled Trial. Trials (2011) 12:255. doi: 10.1186/1745-6215-12-255
6. Ahlering T, Huynh LM, Kaler KS, Williams S, Osann K, Joseph J, et al. Unintended Consequences of Decreased PSA-Based Prostate Cancer Screening. World J Urol (2019) 37(3):489–96. doi: 10.1007/s00345-018-2407-3
7. Kasivisvanathan V, Rannikko AS, Borghi M, Panebianco V, Mynderse LA, Vaarala MH, et al. Mri-Targeted or Standard Biopsy for Prostate-Cancer Diagnosis. New Engl J Med (2018) 378(19):1767–77. doi: 10.1056/NEJMoa1801993
8. Berger I, Annabattula C, Lewis J, Shetty DV, Kam J, Maclean F, et al. (68)Ga-Psma PET/CT vs. mpMRI for Locoregional Prostate Cancer Staging: Correlation With Final Histopathology. Prostate Cancer Prostatic Dis (2018) 21(2):204–11. doi: 10.1038/s41391-018-0048-7
9. Bettermann AS, Zamboglou C, Kiefer S, Jilg CA, Spohn S, Kranz-Rudolph J, et al. [(68)Ga-]Psma-11 PET/CT and Multiparametric MRI for Gross Tumor Volume Delineation in a Slice by Slice Analysis With Whole Mount Histopathology as a Reference Standard - Implications for Focal Radiotherapy Planning in Primary Prostate Cancer. Radiother Oncol J Eur Soc Ther Radiol Oncol (2019) 141:214–9. doi: 10.1016/j.radonc.2019.07.005
10. Chen M, Zhang Q, Zhang C, Zhao X, Marra G, Gao J, et al. Combination of (68)Ga-PSMA PET/CT and Multiparameter Mri Improves the Detection of Clinically Significant Prostate Cancer: A Lesion by Lesion Analysis. J Nucl Med (2019) 60(7):944–9. doi: 10.2967/jnumed.118.221010
11. Eiber M, Weirich G, Holzapfel K, Souvatzoglou M, Haller B, Rauscher I, et al. Simultaneous 68ga-Psma HBED-CC Pet/Mri Improves the Localization of Primary Prostate Cancer. Eur Urol (2016) 70(5):829–36. doi: 10.1016/j.eururo.2015.12.053
12. Scheltema MJ, Chang JI, Stricker PD, van Leeuwen PJ, Nguyen QA, Ho B, et al. Diagnostic Accuracy of (68) Ga-prostate-specific Membrane Antigen (Psma) Positron-Emission Tomography (PET) and Multiparametric (Mp)MRI to Detect Intermediate-Grade Intra-Prostatic Prostate Cancer Using Whole-Mount Pathology: Impact of the Addition of (68) Ga-PSMA PET to Mpmri. BJU Int (2019) 124(Suppl 1):42–9. doi: 10.1111/bju.14794
13. Spohn S, Jaegle C, Fassbender TF, Sprave T, Gkika E, Nicolay NH, et al. Intraindividual Comparison Between (68)Ga-PSMA-PET/CT and mpMRI for Intraprostatic Tumor Delineation in Patients With Primary Prostate Cancer: A Retrospective Analysis in 101 Patients. Eur J Nucl Med Mol Imaging (2020) 47:2796–803. doi: 10.1007/s00259-020-04827-6
14. Hofman MS, Lawrentschuk N, Francis RJ, Tang C, Vela I, Thomas P, et al. Prostate-Specific Membrane Antigen PET-CT in Patients With High-Risk Prostate Cancer Before Curative-Intent Surgery or Radiotherapy (Propsma): A Prospective, Randomised, Multicentre Study. Lancet (London England) (2020) 395(10231):1208–16. doi: 10.1016/S0140-6736(20)30314-7
15. Johnson DC, Raman SS, Mirak SA, Kwan L, Bajgiran AM, Hsu W, et al. Detection of Individual Prostate Cancer Foci Via Multiparametric Magnetic Resonance Imaging. Eur Urol (2019) 75(5):712–20. doi: 10.1016/j.eururo.2018.11.031
16. Priester A, Natarajan S, Khoshnoodi P, Margolis DJ, Raman SS, Reiter RE, et al. Magnetic Resonance Imaging Underestimation of Prostate Cancer Geometry: Use of Patient Specific Molds to Correlate Images With Whole Mount Pathology. J Urol (2017) 197(2):320–6. doi: 10.1016/j.juro.2016.07.084
17. Sachpazidis I, Mavroidis P, Zamboglou C, Klein CM, Grosu A-L, Baltas D. Prostate Cancer Tumour Control Probability Modelling for External Beam Radiotherapy Based on Multi-Parametric Mri-GTV Definition. Radiat Oncol (2020) 15(1):242. doi: 10.1186/s13014-020-01683-4
18. Goodman CD, Fakir H, Pautler S, Chin J, Bauman GS. Dosimetric Evaluation of PSMA Pet-Delineated Dominant Intraprostatic Lesion Simultaneous Infield Boosts. Adv Radiat Oncol (2020) 5(2):212–20. doi: 10.1016/j.adro.2019.09.004
19. Zamboglou C, Sachpazidis I, Koubar K, Drendel V, Wiehle R, Kirste S, et al. Evaluation of Intensity Modulated Radiation Therapy Dose Painting for Localized Prostate Cancer Using (68)Ga-HBED-CC Psma-Pet/Ct: A Planning Study Based on Histopathology Reference. Radiother Oncol (2017) 123(3):472–7. doi: 10.1016/j.radonc.2017.04.021
20. Zamboglou C, Thomann B, Koubar K, Bronsert P, Krauss T, Rischke HC, et al. Focal Dose Escalation for Prostate Cancer Using (68)Ga-Hbed-Cc PSMA PET/CT and MRI: A Planning Study Based on Histology Reference. Radiat Oncol (2018) 13(1):81. doi: 10.1186/s13014-018-1036-8
21. Panettieri V, Rancati T, Onjukka E, Ebert MA, Joseph DJ, Denham JW, et al. External Validation of a Predictive Model of Urethral Strictures for Prostate Patients Treated With HDR Brachytherapy Boost. Front Oncol (2020) 10:910. doi: 10.3389/fonc.2020.00910
22. Thomann B, Sachpazidis I, Koubar K, Zamboglou C, Mavroidis P, Wiehle R, et al. Influence of Inhomogeneous Radiosensitivity Distributions and Intrafractional Organ Movement on the Tumour Control Probability of Focused IMRT in Prostate Cancer. Radiother Oncol (2018) 127(1):62–7. doi: 10.1016/j.radonc.2018.02.006
23. Eder M, Neels O, Müller M, Bauder-Wüst U, Remde Y, Schäfer M, et al. Novel Preclinical and Radiopharmaceutical Aspects of [68Ga]Ga-PSMA-HBED-CC: A New PET Tracer for Imaging of Prostate Cancer. Pharmaceuticals (Basel Switzerland) (2014) 7(7):779–96. doi: 10.3390/ph7070779
24. Surti S, Kuhn A, Werner ME, Perkins AE, Kolthammer J, Karp JS. Performance of Philips Gemini TF Pet/Ct Scanner With Special Consideration for its Time-of-Flight Imaging Capabilities. J Nucl Med Off Publication Soc Nucl Med (2007) 48(3):471–80.
25. Rausch I, Ruiz A, Valverde-Pascual I, Cal-González J, Beyer T, Carrio I. Performance Evaluation of the Vereos Pet/Ct System According to the NEMA Nu2-2012 Standard. J Nucl Med (2019) 60(4):561–7. doi: 10.2967/jnumed.118.215541
26. Zamboglou C, Carles M, Fechter T, Kiefer S, Reichel K, Fassbender TF, et al. Radiomic Features From PSMA PET for non-Invasive Intraprostatic Tumor Discrimination and Characterization in Patients With Intermediate- and High-Risk Prostate Cancer - a Comparison Study With Histology Reference. Theranostics (2019) 9(9):2595–605. doi: 10.7150/thno.32376
27. Weinreb JC, Barentsz JO, Choyke PL, Cornud F, Haider MA, Macura KJ, et al. Pi-Rads Prostate Imaging - Reporting and Data System: 2015, Version 2. Eur Urol (2016) 69(1):16–40. doi: 10.1016/j.eururo.2015.08.052
28. Zamboglou C, Wieser G, Hennies S, Rempel I, Kirste S, Soschynski M, et al. Mri Versus 6 8Ga-Psma PET/CT for Gross Tumour Volume Delineation in Radiation Treatment Planning of Primary Prostate Cancer. Eur J Nucl Med Mol Imaging (2016) 43(5):889–97. doi: 10.1007/s00259-015-3257-5
29. Zamboglou C, Fassbender TF, Steffan L, Schiller F, Fechter T, Carles M, et al. Validation of Different PSMA-PET/CT-Based Contouring Techniques for Intraprostatic Tumor Definition Using Histopathology as Standard of Reference. Radiother Oncol (2019) 141:208–13. doi: 10.1016/j.radonc.2019.07.002
30. Gay HA, Barthold HJ, O’Meara E, Bosch WR, El Naqa I, Al-Lozi R, et al. Pelvic Normal Tissue Contouring Guidelines for Radiation Therapy: A Radiation Therapy Oncology Group Consensus Panel Atlas. Int J Radiat Oncol Biol Phys (2012) 83(3):e353–62. doi: 10.1016/j.ijrobp.2012.01.023
31. Draulans C, van der Heide UA, Haustermans K, Pos FJ, van der Voort van Zyp J, De Boer H, et al. Primary Endpoint Analysis of the Multicentre Phase II Hypo-FLAME Trial for Intermediate and High Risk Prostate Cancer. Radiother Oncol J Eur Soc Ther Radiol Oncol (2020) 147:92–8. doi: 10.1016/j.radonc.2020.03.015
32. Salembier C, Villeirs G, De Bari B, Hoskin P, Pieters BR, Van Vulpen M, et al. Estro ACROP Consensus Guideline on CT- and MRI-based Target Volume Delineation for Primary Radiation Therapy of Localized Prostate Cancer. Radiother Oncol J Eur Soc Ther Radiol Oncol (2018) 127(1):49–61. doi: 10.1016/j.radonc.2018.01.014
33. Zamboglou C, Drendel V, Jilg CA, Rischke HC, Beck TI, Schultze-Seemann W, et al. Comparison of (68)Ga-HBED-CC PSMA-PET/CT and Multiparametric MRI for Gross Tumour Volume Detection in Patients With Primary Prostate Cancer Based on Slice by Slice Comparison With Histopathology. Theranostics (2017) 7(1):228–37. doi: 10.7150/thno.16638
34. Murray JR, Tree AC, Alexander EJ, Sohaib A, Hazell S, Thomas K, et al. Standard and Hypofractionated Dose Escalation to Intraprostatic Tumor Nodules in Localized Prostate Cancer: Efficacy and Toxicity in the DELINEATE Trial. Int J Radiat Oncol (2020) 106(4):715–24. doi: 10.1016/j.ijrobp.2019.11.402
35. Dearnaley D, Syndikus I, Mossop H, Khoo V, Birtle A, Bloomfield D, et al. Conventional Versus Hypofractionated High-Dose Intensity-Modulated Radiotherapy for Prostate Cancer: 5-Year Outcomes of the Randomised, non-Inferiority, Phase 3 CHHiP Trial. Lancet Oncol (2016) 17(8):1047–60. doi: 10.1016/S1470-2045(16)30102-4
36. Wilkins A, Naismith O, Brand D, Fernandez K, Hall E, Dearnaley D, et al. Derivation of Dose/Volume Constraints for the Anorectum From Clinician- and Patient-Reported Outcomes in the CHHiP Trial of Radiation Therapy Fractionation. Int J Radiat Oncol Biol Phys (2020) 106(5):928–38. doi: 10.1016/j.ijrobp.2020.01.003
39. Lind BK, Mavroidis P, Hyödynmaa S, Kappas C. Optimization of the Dose Level for a Given Treatment Plan to Maximize the Complication-Free Tumor Cure. Acta Oncol (Stockholm Sweden) (1999) 38(6):787–98. doi: 10.1080/028418699432950
40. Wheldon TE, Deehan C, Wheldon EG, Barrett A. The Linear-Quadratic Transformation of Dose–Volume Histograms in Fractionated Radiotherapy. Radiother Oncol (1998) 46(3):285–95. doi: 10.1016/S0167-8140(97)00162-X
42. van Lin EN, Fütterer JJ, Heijmink SW, van der Vight LP, Hoffmann AL, van Kollenburg P, et al. Imrt Boost Dose Planning on Dominant Intraprostatic Lesions: Gold Marker-Based Three-Dimensional Fusion of CT With Dynamic Contrast-Enhanced and 1H-Spectroscopic Mri. Int J Radiat Oncol Biol Phys (2006) 65(1):291–303. doi: 10.1016/j.ijrobp.2005.12.046
43. Casares-Magaz O, van der Heide UA, Rørvik J, Steenbergen P, Muren LP. A Tumour Control Probability Model for Radiotherapy of Prostate Cancer Using Magnetic Resonance Imaging-Based Apparent Diffusion Coefficient Maps. Radiother Oncol J Eur Soc Ther Radiol Oncol (2016) 119(1):111–6. doi: 10.1016/j.radonc.2016.02.030
44. Ghobadi G, de Jong J, Hollmann BG, van Triest B, van der Poel HG, Vens C, et al. Histopathology-Derived Modeling of Prostate Cancer Tumor Control Probability: Implications for the Dose to the Tumor and the Gland. Radiother Oncol J Eur Soc Ther Radiol Oncol (2016) 119(1):97–103. doi: 10.1016/j.radonc.2016.02.015
45. Vogelius IR, Bentzen SM. Diminishing Returns From Ultrahypofractionated Radiation Therapy for Prostate Cancer. Int J Radiat Oncol Biol Phys (2020) 107(2):299–304. doi: 10.1016/j.ijrobp.2020.01.010
46. Bostel T, Sachpazidis I, Splinter M, Bougatf N, Fechter T, Zamboglou C, et al. Dosimetric Impact of Interfractional Variations in Prostate Cancer Radiotherapy—Implications for Imaging Frequency and Treatment Adaptation. Front Oncol (2019) 9:(940). doi: 10.3389/fonc.2019.00940
47. Takam R, Bezak E, Yeoh EE, Marcu L. Assessment of Normal Tissue Complications Following Prostate Cancer Irradiation: Comparison of Radiation Treatment Modalities Using NTCP Models. Med Phys (2010) 37(9):5126–37. doi: 10.1118/1.3481514
48. Kuang Y, Wu L, Hirata E, Miyazaki K, Sato M, Kwee SA. Volumetric Modulated Arc Therapy Planning for Primary Prostate Cancer With Selective Intraprostatic Boost Determined by 18F-Choline Pet/Ct. Int J Radiat Oncol Biol Phys (2015) 91(5):1017–25. doi: 10.1016/j.ijrobp.2014.12.052
49. Liu M, Moiseenko V, Agranovich A, Karvat A, Kwan W, Saleh ZH, et al. Normal Tissue Complication Probability (NTCP) Modeling of Late Rectal Bleeding Following External Beam Radiotherapy for Prostate Cancer: A Test of the QUANTEC-recommended NTCP Model. Acta Oncol (Stockholm Sweden) (2010) 49(7):1040–4. doi: 10.3109/0284186X.2010.509736
50. Peeters ST, Hoogeman MS, Heemsbergen WD, Hart AA, Koper PC, Lebesque JV. Rectal Bleeding, Fecal Incontinence, and High Stool Frequency After Conformal Radiotherapy for Prostate Cancer: Normal Tissue Complication Probability Modeling. Int J Radiat Oncol Biol Phys (2006) 66(1):11–9. doi: 10.1016/j.ijrobp.2006.03.034
51. Sachpazidis I, Hense J, Mavroidis P, Gainey M, Baltas D. Investigating the Role of Constrained CVT and CVT in HIPO Inverse Planning for HDR Brachytherapy of Prostate Cancer. Med Phys (2019) 46(7):2955–68. doi: 10.1002/mp.13564
52. Splinter M, Bostel T, Sachpazidis I, Fechter T, Zamboglou C, Jäkel O, et al. Dosimetric Impact of Interfractional Variations for Post-prostatectomy Radiotherapy to the Prostatic Fossa—Relevance for the Frequency of Position Verification Imaging and Treatment Adaptation. Fron Oncol (2019) 9(1191). doi: 10.3389/fonc.2019.01191
53. Mavroidis P, Ferreira BC, Papanikolaou N, Lopes Mdo C. Analysis of Fractionation Correction Methodologies for Multiple Phase Treatment Plans in Radiation Therapy. Med Phys (2013) 40(3):031715. doi: 10.1118/1.4792636
54. Roland T, Mavroidis P, Gutierrez A, Goytia V, Papanikolaou N. A Radiobiological Analysis of the Effect of 3D Versus 4d Image-Based Planning in Lung Cancer Radiotherapy. Phys Med Biol (2009) 54(18):5509–23. doi: 10.1088/0031-9155/54/18/011
55. Deb P, Fielding A. Radiobiological Model Comparison of 3D Conformal Radiotherapy and IMRT Plans for the Treatment of Prostate Cancer. Australas Phys Eng Sci Med (2009) 32(2):51. doi: 10.1007/BF03178629
57. Langen KM, Willoughby TR, Meeks SL, Santhanam A, Cunningham A, Levine L, et al. Observations on Real-Time Prostate Gland Motion Using Electromagnetic Tracking. Int J Radiat Oncol Biol Phys (2008) 71(4):1084–90. doi: 10.1016/j.ijrobp.2007.11.054
58. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing (2020). Available at: https://www.R-project.org/.
59. Yamada Y, Rogers L, Demanes DJ, Morton G, Prestidge BR, Pouliot J, et al. American Brachytherapy Society Consensus Guidelines for High-Dose-Rate Prostate Brachytherapy. Brachytherapy (2012) 11(1):20–32. doi: 10.1016/j.brachy.2011.09.008
60. Zilli T, Jorcano S, Bral S, Rubio C, Bruynzeel AME, Oliveira A, et al. Once-a-Week or Every-Other-Day Urethra-Sparing Prostate Cancer Stereotactic Body Radiotherapy, a Randomized Phase II Trial: 18 Months Follow-Up Results. Cancer Med (2020) 9(9):3097–106. doi: 10.1002/cam4.2966
61. Turkbey B, Rosenkrantz AB, Haider MA, Padhani AR, Villeirs G, Macura KJ, et al. Prostate Imaging Reporting and Data System Version 2.1: 2019 Update of Prostate Imaging Reporting and Data System Version 2. Eur Urol (2019) 76(3):340–51. doi: 10.1016/j.eururo.2019.02.033
62. Leibovich BC, Blute ML, Bostwick DG, Wilson TM, Pisansky TM, Davis BJ, et al. Proximity of Prostate Cancer to the Urethra: Implications for Minimally Invasive Ablative Therapies. Urology (2000) 56(5):726–9. doi: 10.1016/S0090-4295(00)00792-5
63. Wang T, Zhou J, Tian S, Wang Y, Patel P, Jani AB, et al. A Planning Study of Focal Dose Escalations to Multiparametric MRI-Defined Dominant Intraprostatic Lesions in Prostate Proton Radiation Therapy. Br J Radiol (2020) 93(1107):20190845. doi: 10.1259/bjr.20190845
64. Mohammed N, Kestin L, Ghilezan M, Krauss D, Vicini F, Brabbins D, et al. Comparison of Acute and Late Toxicities for Three Modern High-Dose Radiation Treatment Techniques for Localized Prostate Cancer. Int J Radiat Oncol Biol Phys (2012) 82(1):204–12. doi: 10.1016/j.ijrobp.2010.10.009
66. King CR, Brooks JD, Gill H, Presti JC. Long-Term Outcomes From a Prospective Trial of Stereotactic Body Radiotherapy for Low-Risk Prostate Cancer. Int J Radiat Oncol Biol Phys (2012) 82(2):877–82. doi: 10.1016/j.ijrobp.2010.11.054
68. Ullrich T, Quentin M, Oelers C, Dietzel F, Sawicki LM, Arsov C, et al. Magnetic Resonance Imaging of the Prostate At 1.5 Versus 3.0t: A Prospective Comparison Study of Image Quality. Eur J Radiol (2017) 90:192–7. doi: 10.1016/j.ejrad.2017.02.044
69. Parker CC, Damyanovich A, Haycocks T, Haider M, Bayley A, Catton CN. Magnetic Resonance Imaging in the Radiation Treatment Planning of Localized Prostate Cancer Using Intra-Prostatic Fiducial Markers for Computed Tomography Co-Registration. Radiother Oncol (2003) 66(2):217–24. doi: 10.1016/S0167-8140(02)00407-3
70. Spohn SKB, Kramer M, Kiefer S, Bronsert P, Sigle A, Schultze-Seemann W, et al. Comparison of Manual and Semi-Automatic [18f]Psma-1007 PET Based Contouring Techniques for Intraprostatic Tumor Delineation in Patients With Primary Prostate Cancer and Validation With Histopathology as Standard of Reference. Front Oncol (2020) 10:2725. doi: 10.3389/fonc.2020.600690
71. Wegener D, Zips D, Thorwarth D, Weiß J, Othman AE, Grosse U, et al. Precision of T2 Tse MRI-CT-image Fusions Based on Gold Fiducials and Repetitive T2 Tse Mri-MRI-fusions for Adaptive IGRT of Prostate Cancer by Using Phantom and Patient Data. Acta Oncol (2019) 58(1):88–94. doi: 10.1080/0284186X.2018.1518594
72. Kuten J, Fahoum I, Savin Z, Shamni O, Gitstein G, Hershkovitz D, et al. Head-to-Head Comparison of (68)Ga-PSMA-11 With (18)F-Psma-1007 PET/CT in Staging Prostate Cancer Using Histopathology and Immunohistochemical Analysis as a Reference Standard. J Nucl Med Off Publication Soc Nucl Med (2020) 61(4):527–32. doi: 10.2967/jnumed.119.234187
73. Draulans C, De Roover R, van der Heide UA, Kerkmeijer L, Smeenk RJ, Pos F, et al. Optimal 68Ga-PSMA and 18F-PSMA PET Window Levelling for Gross Tumour Volume Delineation in Primary Prostate Cancer. Eur J Nucl Med Mol Imaging (2021) 48:1211–8. doi: 10.1007/s00259-020-05059-4
74. Sabater S, Pastor-Juan MR, Andres I, López-Martinez L, Lopez-Honrubia V, Tercero-Azorin MI, et al. Mri Prostate Contouring is Not Impaired by the Use of a Radiotherapy Image Acquisition Set-Up. An Intra- and Inter-Observer Paired Comparative Analysis With Diagnostic Set-Up Images. Cancer/Radiothérapie (2021) 25(2):107–13. doi: 10.1016/j.canrad.2020.05.024
Keywords: hypofractionated radiotherapy, PSMA - prostate specific membrane antigen, focal dose escalation, tumor control probability (TCP), NTCP (normal tissue complication probability) model, mpMRI, primary prostate cancer, histopathology
Citation: Spohn SKB, Sachpazidis I, Wiehle R, Thomann B, Sigle A, Bronsert P, Ruf J, Benndorf M, Nicolay NH, Sprave T, Grosu AL, Baltas D and Zamboglou C (2021) Influence of Urethra Sparing on Tumor Control Probability and Normal Tissue Complication Probability in Focal Dose Escalated Hypofractionated Radiotherapy: A Planning Study Based on Histopathology Reference. Front. Oncol. 11:652678. doi: 10.3389/fonc.2021.652678
Received: 12 January 2021; Accepted: 08 April 2021;
Published: 14 May 2021.
Edited by:Tone Frost Bathen, Norwegian University of Science and Technology, Norway
Reviewed by:Torgrim Tandstad, St. Olavs University Hospital, Norway
Ingerid Skjei Knudtsen, Norwegian University of Science and Technology, Norway
Kathrine Røe Redalen, Norwegian University of Science and Technology, Norway
Copyright © 2021 Spohn, Sachpazidis, Wiehle, Thomann, Sigle, Bronsert, Ruf, Benndorf, Nicolay, Sprave, Grosu, Baltas and Zamboglou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Simon K. B. Spohn, Simon.Spohn@uniklinik-freiburg.de
†These authors have contributed equally to this work and share first authorship