Abstract
Background/Aim: We examined the prognostic value of intraprostatic gross tumour volume (GTV) as measured by multiparametric MRI (mpMRI) in patients with prostate cancer following (primary) external beam radiation therapy (EBRT). Patients and Methods: In a retrospective monocentric study, we analysed patients with prostate cancer (PCa) after EBRT. GTV was delineated in pre-treatment mpMRI (GTV-MRI) using T2-weighted images. Cox-regression analyses were performed considering biochemical failure recurrence-free survival (BRFS) as outcome variable. Results: Among 131 patients, after a median follow-up of 57 months, biochemical failure occurred in 27 (21%). GTV-MRI was not correlated with % of positive biopsy cores, Gleason score and initial PSA (all r<0.2) and only moderately correlated with cT stage (r=0.32). In univariate analysis, cT stage, Gleason score and GTV-MRI were higher in subjects with shorter BRFS (p<0.05). GTV-MRI remained a significant predictor for BRFS in multivariate analyses, independent of Gleason score and cT stage. Conclusion: GTV, defined using mpMRI, provides incremental prognostic value for BRFS, independent of established risk factors. This supports the implementation of imaging-based GTV for risk-stratification, although further validation is needed.
For patients with primary prostate cancer (PCa), an accurate risk stratification is mandatory to select a suitable treatment. The National Comprehensive Cancer Network risk-classification system for primary PCa patients lists clinical T stage (cT stage), Gleason score, the percentage of positive biopsy cores (PPBC) and initial prostate-specific antigen (iPSA) levels as important factors (1). Each has been shown to be an independent risk factor of biochemical failure after definitive radiotherapy (RT) (2). However, this model does not always predict the final outcome and its predictions may be very heterogeneous (3). Consequently, more robust risk-factors are urgently needed.
The implementation of multiparametric magnetic resonance imaging (mpMRI) shows promising results in the delineation of the intraprostatic gross tumour volume (GTV) and is considered as the standard of care diagnostic imaging modality for patients with primary PCa (1). Three studies examined whether local recurrences of PCa after primary RT occur at the location of the primary lesion using pre- and post-treatment MRI in a limited number of patients. Each study concluded that local recurrence after RT appears mostly at the primary tumour's location (4-6). Similarly, two studies postulated that even if mpMRI may not detect the entire PCa, it is able to show a clinically sufficient amount of the tumour (7, 8). Additionally, a recent study postulated that the RT dose distribution in the MR-defined GTV correlates with outcome after definitive external beam radiotherapy (EBRT) in PCa patients (9).
By calculating the percentage PPBC (10, 11) or by measuring the PCa volume in a surgical specimen (12, 13), several studies have proposed that the amount of malignant tissue within the prostate is an independent prognostic factor for patients with PCa. However, it is still not known whether the volume of imaging-defined PCa is a prognostic factor of the outcome after primary RT. First studies proposed that the intraprostatic tumour extension defined by MRI may be used as a surrogate marker for outcome in patients treated with primary RT (14, 15). However, in these studies, only tumour extension in MR images was estimated. In this study, we examined whether the clinically relevant tumour mass can be quantified using GTV that are manually outlined on mpMR images. Consequently, we tested GTV correlation with established prognostic risk-factors and whether GTV may influence the outcome of PCa patients after definitive EBRT.
Patients and Methods
Patients. In this retrospective, single institution analysis, all patients with a localized and histologically proven PCa who received EBRT with or without androgen deprivation therapy (ADT) were screened from February 2008 to October 2016. The availability of mpMRI data at the maximum of 6 months prior to EBRT was obligatory. Patients were excluded from the analysis in case of cN1 or cM1 disease, EBRT of the pelvic lymph nodes, and initial PSA serum values above 50 ng/ml. The administration of ADT over 1 month prior to conduction of MRI scans was also an exclusion criterion. Finally, 131 patients were included in the analysis.
Magnetic resonance imaging. In 116 patients (88.5%), MR images were acquired on various 1.5 Tesla systems (Siemens Healthcare, Berlin and Munich, Germany), and for 15 patients (11.5%) a 3 Tesla system (Tim Trio, Siemens Healthcare) was used. For signal reception, a body array surface coil in combination with an integrated spine array coil was used, and no endo-rectal coil was applied. The mpMRI imaging protocol contained a T2-weighted fast spin echo (T2w-TSE) acquisition, diffusion weighted imaging (DWI) and dynamic contrast-enhanced (DCE) perfusion MRI - a detailed list of all acquisition parameters of the mpMRI protocol is given in (16). In case of multiple mpMRI scans before RT, the last pre-treatment scan was selected.
Treatment protocol. EBRT was delivered in 39 (30%) and 92 (70%) patients with 3D-conformal and intensity-modu lated RT (IMRT), respectively. All patients had image-guided RT (IGRT) using daily 2D/2D imaging and at least one cone-beam CT per week. Intraprostatic fiducial markers were implanted in 123 (94%) patients prior to EBRT. Using the cone-beam CT information, organs at risk as well as target volumes were adapted during IMRT. According to the patient risk stratification the clinical target volume (CTV) was defined. In low-risk patients the prostatic gland and the base of the seminal vesicles were defined as CTV. In intermediate- and high-risk patients the prostatic gland and at least half of the seminal vesicles were expanded by 4-5 mm (excluding bladder and rectum) to create the CTV. The respective CTVs were expanded with 6 mm to create the planning target volume (PTV).
The aimed and the delivered median prescription dose to the PTV were 76 and 74 Gy (range=66-78 Gy), respectively. Ten patients had single fractions of 1.8 Gy, and 121 patients had single fractions of 2 Gy. A total of 71 patients (54%) received ADT parallel to EBRT for a median duration of 7 months (range=3-24 months). In the follow-up, patients were seen every 3-6 months for the first two years, and every 6-12 months thereafter for physical examination, PSA measurements and radiological examination, if necessary.
Image co-registration and generation of contours. Axial T2w sequences and planning CT scans were matched in a RT planning system (Eclipse v13.5, Varian, Palo Alto, CA, USA) using mutual information registration. In case visual assessment showed anatomical mismatch, a manual adjustment was performed. The usage of an axial T2w sequence and at least one DWI- or DCE sequence was obligatory. From the DWI data, the calculated apparent diffusion coefficient (ADC) maps were registered, whereas from the DCE time series a post-injection time frame was manually selected for having a maximum contrast in the prostate. For further alignment between T2w images and the DWI- or DCE sequences the respective pre-set (hardware based) registrations were used.
The prostatic gland (PG) was contoured using CT and T2w-MRI information. The radiology reports of all MRI scans were available during contouring. T2W images were available for all 131 patients, 81 (61.8%) had DWI/ADC images, and for 50 patients (38.2%) DCE information was additionally available. Considering the PI-RADs v2 criteria (17), a GTV was delineated based on the mpMRI data (GTV-MRI) by two readers individually (Figure 1): CZ (GTV-MRI_exp) and BO (GTV-MRI_res). CZ is a board licensed radiation oncologist with over 5 years of experience in prostate MR image interpretation. Additionally, CZ had training according to Rischke et al. (18). BO is a second-year resident at the department of radiology without any additional training in prostate MR image interpretation.
Data analysis and statistics. According to the Phoenix criteria, biochemical failure after EBRT (19) was defined as the study endpoint. The GTV volumes below and above 6 ml were considered as continuous and dichotomous variables - the cut-off point of 6 ml was chosen according to literature (20, 21). Univariate and multivariate Cox regression analyses were performed considering bio chemical failure recurrence-free survival (BRFS).
For graphical representation, the Kaplan-Meier survival curves for BRFS were used. The Kaplan-Meier curves for the two patient groups (GTV ≤6 ml vs. GTV >6 ml) were compared by a log-rank test. Correlation analyses were performed with nonparametric Spearman rho test. All tests were considered to be statistically significant at p<0.05. Statistical analysis was conducted with SPSS v25 (IBM, Armonk, NY, USA).
Results
Patient and treatment characteristics. Detailed characteristics of the study cohort are listed in Table I. The majority of patients had intermediate-risk disease according to NCNN classification. The median volume of GTV-MRI_exp and GTV-MRI_res was 2.6 ml (range=0-38 ml) and 2.5 ml (range=0.2-60.4 ml), respectively (p=0.039 in Wilcoxon matched-pairs test). Among 131 patients (median age=74; range=56-85 years), after a median follow-up of 57 months (range=14-138 months) following EBRT, biochemical failure occurred in 27 (21%). Median PSA level at the time of recurrence was 4.5 ng/ml (range=2.3-14.9 ng/ml). At the time of last evaluation 2 patients died of PCa. Univariate Cox-regression showed that the prescription dose, the RT technique (IMRT vs. 3D) and the usage of ADT had no impact on BRFS (p>0.05).
Correlation between GTV-MRI and established prognostic factors for BRFS. For the GTV contours of both readers a very weak correlation was observed with PPBC, Gleason score and initial PSA (Spearman r<20, p>0.05). A weak correlation with the cT stage was found for GTV-MRI_exp (Spearman r=0.31, p<0.01) and GTV-MRI_res (Spearman r=0.33, p<0.01), respectively. The correlation between GTV-MRI_exp volumes and clinical parameters is presented in Figure 2.
Impact of clinical parameters on BRFS. The Kaplan-Meier curves (Figure 3) demonstrate the impact of absolute GTV-MRI volumes (cut-off: 6 ml) on BRFS. Table II shows the results of Cox regression analyses to assess the impact of different patient-related factors on BRFS. Univariate analyses revealed that 3 out of 10 factors were statistically significant (p<0.05): volume of GTV-MRI, Gleason score and cT stage. In multi variate analysis, including all significant parameters from the univariate analysis, the volume of GTV-MRI_exp (cut-off: 6ml and continuous variable), the volume of GTV-MRI_res (cut-off: 6ml) and the Gleason score were independent risk factors for BRFS.
Discussion
The NCCN risk-stratification system is widely used for patients with primary PCa, although it may not always predict the outcome correctly (3). To improve risk-stratification systems, several studies showed that the PCa amount has an impact on the outcome of patients treated with prostatectomy or EBRT for primary PCa (10-12, 22). However, in most patients treated with primary RT, the size of the PCa was indirectly assessed by PPBC (23) or by disease localization based on MRI imaging information (24-26). Kaufmann et al. reported that the largest axial tumour diameter in MR imaging was strongly associated with biochemical outcomes after RT for primary PCa (26). In this study, we postulated that mpMRI may define the significant PCa amount. Likewise, we hypothesized that the volume of the manually outlined GTV-MRI should predict the outcome in terms of biochemical failure after primary EBRT.
With multivariate Cox regression analyses, we demonstrated that a tumour volume of more than 6 ml is an independent prognostic factor for BRFS, suggesting an up to 3-fold increase in the likelihood of biochemical tumour recurrence. Additionally, when defined by an expert, the absolute volume can be used as a continuous variable where each 1 ml-increase of GTV volume leads to an about 10% increase in BRFS risk. These findings may be explained by the increased number of clonogenic PCa cells (27) or by a possible increase in the number of hypoxic regions in larger tumours (28). Furthermore, McNeal et al. postulated that the PCa volume in surgery specimens strongly correlates with nodal metastases (29). Kim et al. (12) used a threshold of 5 ml to categorize extensive PCa in prostatectomy specimens, which is close to our 6 ml-threshold. The cut-off value of 6 ml was also suggested as a predictive factor for worse outcomes by the prostatectomy studies of Stamney et al. (20, 21).
Several previous studies had a huge interobserver variability for manually delineated GTVs based on MR images (18, 30). In this study, two observers with different levels of experience were contouring the GTVs. No statistically significant differences in absolute GTV-volumes were observed, and the cut-off value of 6 ml was a significant predictor of BRFS for both readers. However, only the expert contour had a significant impact on BRFS as a continuous variable, indicating that the variations in the contours drawn by the novice may lead to fluctuations that mask the correlation between GTV volume and outcome. Future studies should examine whether implementation of automatic GTV segmentation methods based on deep convolutional neuronal networks may improve the robustness of GTV segmentation (31).
In our study, no correlation was found between iPSA serum values and GTV-MRI, which may be explained by the non-specific PSA expression of non-malignant cells in the prostate and by other non-tumour specific factors that can influence PSA expression (32). Interestingly, only a weak correlation between cT stage and GTV-MRI was observed: patients with cT3b disease had lower median GTV-MRI values compared to cT3a PCa patients.
GTV-MRI did not correlate with Gleason score, which is consistent with Donato et al. who did not find a correlation between PCa volume in ADC-MRI and Gleason score (33). Furthermore, no correlation between PPBC and GTV-MRI was detected, and some of the patients with PPBC >80% had low GTV-MRI volumes (<2 ml). It should be mentioned that an exact reconstruction of the intraprostatic tumour mass is not assured in trans-rectal biopsies (34), which may explain the missing correlation between PPBC and GTV volume. Only 4% of the patients received a trans perineal TRUS-fusion biopsy (guided by mpMRI or PSMA PET/CT). Future studies should analyse whether the volume of GTV correlates better with PPBC and Gleason score after transperineal TRUS-fusion biopsies. Additionally, the usage of ADT between biopsy and pre-RT images may serve as an explanation for the lack of correlation between PPBC and GTV-MRI. ADT administration might shrink the volume of the lesion (35), moving some patients under the threshold even if the original volume of the GTV was ≥6 mm. In our study, 25 of the 71 patients (35%) with ADT received mpMRI scans during ADT. Future studies are needed to evaluate whether GTV for risk-stratification should be derived before or after ADT. In summary, our findings suggest that GTV-MRI volume may offer complementary information to the NCCN risk-stratification system as an independent risk factor. It would be interesting to further investigate the correlation of other parameters extracted from GTV-MRI and the corresponding prognosis. Two studies proposed the usage of diffusion-weighted MRI as an imaging biomarker for monitoring response of PCa to primary RT (36, 37). Turchan et al. showed that PCa lesions with PI-RADs score of 5 increase the risk of failure after RT for primary PCa (38).
This study has several limitations due to its retrospective character: treatment protocols (e.g. ADT duration, RT technique) were not identical for all patients. However, RT technique, RT dose and the usage of ADT had no significant impact on BRFS in our study. Furthermore, MRI systems and acquisition protocols were not identical for the entire group, but imaging contrasts such as T2, diffusion coefficient and perfusion-induced contrast change were used that are robust against different field strengths and imaging hardware. Despite the usage of different imaging methods and imaging protocols, we showed that image-defined GTV had a significant impact on BRFS. Interestingly, we reproduced the cut-off value of 6 ml for GTV, which was already described by prostatectomy studies. Additionally, the usage of different MRI systems and protocols reflects the current clinical situation, and every new generalizable biomarker for risk prediction has to be robust in various clinical settings: a mandatory condition for the implementation of imaging-derived GTV for risk stratification systems must be a maximum independence of the obtained results from the utilised imaging system. Another limitation of our study is the relatively short follow-up time (approx. 5 years), which did not enable the evaluation of other relevant endpoints like PCa-specific survival or overall survival. However, other similar studies had even shorter follow-up times of 43-44 months (24, 39). Although mpMRI imaging is the standard of care for primary PCa staging, a higher sensitivity for PSMA PET/CT in GTV detection was reported (40, 41). Future studies should address whether PSMA-based GTV may also be used for risk prediction of patients with primary PCa. To account for these issues, our findings need external validation, preferably from prospective studies with longer follow-up.
This study showed that GTV defined on mpMRI has a weak correlation with known prognostic factors for risk stratification in patients with primary PCa. GTV derived from pre-therapeutic MR images is an independent predictive factor for biochemical failure after EBRT. These findings support the implementation of GTV delineation based on modern imaging tools for risk-stratification systems in patients with primary PCa.
Footnotes
Authors' Contributions
BO acquired clinical follow-up data, delineated GTV in mpMRI images, carried out statistical analysis and was a major contributor to the final manuscript. MVB acquired primary, clinical follow-up data and patients. CK acquired primary and clinical follow-up data, delineated GTV in mpMRI images. SK acquired primary and clinical follow-up data. NN acquired primary and clinical follow-up data and was a major contributor to the final manuscript. TS acquired primary and clinical follow-up data. SS acquired primary and clinical follow-up data. DDG carried out statistical analysis. LH extracted imaging information and processed delineated GTV. LB extracted imaging information and processed delineated GTV. MBo carried out statistical analysis and was a major contributor to the final manuscript. AG acquired primary and clinical follow-up data and was a major contributor to the final manuscript. FB supervised imaging acquisition and was a major contributor to the final manuscript. MBe supervised imaging acquisition and was a major contributor to the final manuscript. CZ supervised imaging acquisition acquired primary and follow-up clinical data, delineated GTV in mpMRI images, carried out statistical analysis and was a major contributor to the final manuscript.
This article is freely accessible online.
Conflicts of Interest
All Authors disclose that no financial and personal relationships with other people or organisations that could inappropriately influence the work exist.
- Received July 15, 2020.
- Revision received July 29, 2020.
- Accepted August 5, 2020.
- Copyright© 2020, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved