Introduction

Radical prostatectomy (RP) is a standard first-line treatment option for localized prostate cancer (PCa), with nerve-sparing surgery restricted to preoperatively potent patients without high-grade prostatic cancer (HGPC) [1,2,3]. Less invasive treatment options (e.g., active surveillance or high-intensity focalized ultrasound (HIFU) ablation) are available but appropriate only for less aggressive cancers but with a recurrence rate that remains high at approximately 30% [4]. This may be due, in part, to the difficulty of selecting the most appropriate patients for each treatment option, which remains problematic even with the benefit of nuclear magnetic resonance imaging [5]. To improve outcomes, patient selection needs to be further refined; in this context, androgen hormonal status merits further consideration.

The impact of hypogonadism or testosterone deficiency on localized PCa remains controversial [6, 7]. Many studies have methodological limitations, e.g., small sample size or suboptimal assay protocols. Furthermore, in most studies, hypo-gonadism is based solely on total serum testosterone (TT) levels, without considering bioactive (bioavailable) testosterone (BT) which accounts for approximately 50% of TT [6]. Similarly, the complexities of the androgen pathway leading to testosterone (and allied hormonal levels) are rarely fully considered. Practical issues, such as the recommendation to perform blood collection in the morning between 7 and 11 a.m., are also often overlooked [8, 9], while radioimmunoassay (RIA) has largely been superseded by gas chromatography-mass spectrometry (GC-MS) [10, 11].

The aim of this prospective, multicenter study was to correlate aggressiveness of the PCa, as defined by the predominant Gleason pattern (prdGP) with androgen levels and behavioral/quality of life (QoL) assessments [13, 14].

Methods

Study Design and Study Population

This was a prospective cohort study, drawn from four large French metropolitan hospitals between June 2013 and June 2016, recruiting consecutive patients scheduled to undergo RP (in accordance with both the EAU and AUA recommendations for clinically localized PCa (T1-T2c)) [2, 3]. We excluded patients with previous or ongoing local treatment (e.g., radiotherapy, phototherapy, thermotherapy, HIFU) or systemic treatment likely to interfere with hormonal status (e.g., androgen receptor blockers, LHRH agonists/analogues, testosterone supplementation) or with comorbidities likely to affect gonadal status. No age or PSA level limitations were set.

The study protocol was approved by the institutional review boards of all centers. The trial was registered (NCT02235142) and conducted in accordance with the Declaration of Helsinki. Written informed consent was obtained from all patients.

Assessments

Following admission for RP, demographic and anthropometric data and routine laboratory samples (including metabolic syndrome parameters) were collected on all patients. For QoL assessment, patients completed two different questionnaires; the aging male symptom (AMS) scale and the International Index of Erectile Function 5-item (IIEF-5) questionnaire [12, 13].

Assessment of Androgen Levels and PSA

Immediately prior to surgery, blood samples (30 mL) were taken between 7 and 10 a.m. following an overnight fast, as recommended in the Endocrine Society guidelines [8] and stored at − 20 °C before assay by a single laboratory. PSA, follicle-stimulating hormone (FSH), luteinizing hormone (LH), and sex hormone–binding globulin (SHBG) were measured by RIA, and TT, dihydrotestosterone (DHT), dehydroepiandrosterone (DHEA), Δ5 (androstenediol, D5), ∆4 (androstenedione, D4), estrone (E1), estradiol (E2), and DHEA sulfate (DHEA-S) were measured by GC-MS [11]. BT and free testosterone (FT) were calculated through a standardized formula using specific association constants of testosterone for SHBG and albumin. The threshold for hypogonadism was set at TT < 300 ng/dL (< 10.4 nmol/L), with a corresponding BT level set at the lower level of physiologic BT (< 80 ng/dL [< 2.8 nmol/L]).

Prostate Histological Study

Preoperative prostate core biopsies (12 from each patient) and RP specimens were processed and analyzed according to standard procedures [14]. As PCa aggressiveness is not sufficiently represented by a single Gleason score—which does not distinguish between (prognostically different) 3 + 4 and 4 + 3 scores—we used prdGP grading, with prdGP4 defining HGPC or ISUP ≥ 3, an approach proposed by Epstein et al. [15]. To ensure reporting consistency, core biopsies and RP specimens were assessed by a single pathologist at each center, then blindly centrally reviewed by a single reference uropathologist, with any differences resolved by consensus.

Data Analysis

Based on previous experience [16], we anticipated that prdGP4 would be identified in approximately one third of surgical resections. Using a margin of error of < 2.5% with 95% confidence intervals (CI), we calculated a required sample size of approximately 1350 patients; accounting for an estimated attrition rate of 20%, the recruitment target was set at 1620.

For continuous variables, data are presented as median (25th–75th percentile); for categorical variables, counts and percentages are presented. The cohort was classified into two groups according to the prdGP3 or prdGP4 grading on RP specimens. For comparisons between different levels of hypogonadism, patients were stratified into four groups: (i) “T+, B+” (TT ≥ 300 ng/dL; BT ≥ 80 ng/dL); (ii) “T+, B−” (TT ≥ 300 ng/dL; BT < 80 ng/dL); (iii) “T−, B+” (TT < 300 ng/dL; BT ≥ 80 ng/dL); (iv) “T−, B−” (TT < 300 ng/dL; BT < 80 ng/dL).

For continuous variables, differences between groups were analyzed using the non-parametric Kruskal-Wallis test, and when global p values were significant, pairwise Dunn tests were performed. The Kruskal-Wallis test was also used to compare lower grade (prdGP3) and higher grade (prdGP4) PCa. For categorical variables, Fisher’s exact test (2 × 2 or r × c) was used, and when a global significance was present, 2 × 2 Fisher’s exact test was used for pairwise comparisons. For multiple pairwise comparisons, the Bonferroni-Holms correction was used.

Multivariable logistic regression analysis was used to identify factors associated with the likelihood of upgrading from prdGP3 on biopsy to prdGP4 on RP specimens. For patients with prdGP3 on biopsy, a total of 110 logistic regressions were performed using all continuous and categorical parameter values available prior to surgery; for each regression, a random sample of 55% of cases was drawn to construct the model, with the remaining 45% used for model validation [17]. Prediction of upgrading to prdG4 on RP was at least as good as chance in 24 models, and those parameters that featured in at least two of these models were retained for the final logistic regression. Receiver operator characteristic (ROC) curve analysis was used to determine the predictive accuracy of this multivariate model.

Statistical analysis was conducted using NCSS, version 11 (NCSS Ltd., Kaysville, UT, USA) and R (version 3.3.2; The R Foundation, Vienna, Austria). A two-sided p value < 0.05 was used for statistical significance.

Results

Study Population

During the enrollment period, 2073 RPs were performed. Of these, 1611 patients were enrolled and 1343 patients were evaluable. Twenty patients (1.2%) who did not undergo RP, 147 patients (9.1%) with protocol deviations, 24 patients (1.5%) with unevaluable hormonal blood samples, and 80 patients (5.0%) who withdrew consent were excluded from the final analysis (Fig. 1).

Fig. 1
figure 1

Study flow

Demographic and Clinic-Pathological Characteristics

Of the 1343 patients, 68% (n = 912) had prdGP3 (i.e., ISUP ≤ 2) tumors and 32% with prdGP4 (i.e., ISUP ≥ 3) tumors (n = 431). Characteristics of the overall cohort and stratified by prdGP grading are shown in Table 1. Compared to patients with prdGP3 tumors, patients with prdGP4 PCa were older (by approximately 1.5 years), with significantly higher PSA levels (by approximately 2.5 units), with a lower proportion of tumors at T1c stage (44% vs. 57%) and a higher proportion at clinical/nodal stages (pT ≥ 3a, pN ≥ 1). No differences in medical comorbidities and allied biochemistry were apparent, and no differences in AMS scale and IIEF-5 questionnaires were observed.

Table 1 Population characteristics: stratification between prdGP3 and prdGP4 on radical prostatectomy specimen

Androgen Characteristics

Clear differences in androgen pathway hormone levels were seen with different tumor grades (Table 1). Compared with patients with prdGP3 tumors, significantly more patients with prdGP4 cancers had demonstrable hypogonadism, characterized either by BT levels (< 80 ng/dL; 17.4% vs. 10.7%, p < 0.001) or by TT levels (< 300 ng/dL; 14.2% vs. 9.7%, p = 0.02). BT levels were also lower in those patients with prdGP4 tumors compared with those in prdGP3 patients (119 ng/dL [25th–75th percentile, 92–150 ng/dL] vs. 124 ng/dL [25th–75th percentile, 99–156 ng/dL] respectively, p = 0.008), as were FT values (and conversely substantially higher SHBG levels). Median LH and TT values did not differ significantly between prdGP3 and prdGP4 patients, with no significant difference in other steroid levels, regardless of an adrenal (i.e., DHEA, ∆5) or testicular origin (i.e., ∆4, E2, E1, or DHT).

Compared with eugonadal patients (i.e., “T+, B+”), both groups of hypogonadal patients (TT < 300 ng/dL) had significantly higher values for the main markers of the metabolic syndrome and were more likely to have cardiovascular disease, although there was no significant difference in the prevalence of diabetes (Table 2). While no significant differences in PSA values or the incidence of prdGP4 on presurgery biopsies between hypogonadal and eugonadal groups were apparent, significantly more RP specimens were graded as prdGP4 in fully hypogonadal patients (T−B−) compared to fully eugonadal (T+B+) subjects (47.5% vs. 29.9%, p = 0.004). The correlation between PSA and BT or TT was extremely low (Spearman’s coefficient of 0.018 and 0.017, respectively). The proportion of patients with prdGP4 tumors in both groups with low BT (< 80 ng/dL) was about 10% greater than in groups with higher BT, and higher tumor clinical/nodal stages were also seen in hypogonadal patients (Table 2). Significant differences in SHBG concentrations were seen for every pairwise comparison, and SHBG concentrations were higher in both TT ≥ 300 ng/dL groups. The global p values for differences in the androgen pathway hormone levels (DHEA, DHEA−S, ∆4, ∆5, TT, BT, FT, DHT, E1, E2) were all highly significant. Higher AMS sexual and IIEF-5 scores were observed in eugonadal patients compared to that in hypogonadal patients.

Table 2 Patient characteristics according to their gonadal status

Multivariate logistic regression was performed on all presurgical parameters and replicated to select parameters that retained in at least 2 out of 24 models with suitable predictive value for prdGP4 on RP resection. Increases in age, PSA, triglycerides, and FSH; decreased LH and BT levels; non-Caucasian ethnicity; and coexisting cardiovascular disease each independently increase the risk of upgrading from prdGP3 on biopsy to prdGP4 on RP specimen and were included in the final model (Table 3). The corresponding integrated ROC curve for the final model is shown in Fig. 2. The model provided a predictive accuracy of 61%.

Table 3 Multivariable logistic regression assessing independent predictors of prdGP4 on radical prostatectomy specimen in patients with prdGP3 on biopsy
Fig. 2
figure 2

ROC curve analysis for model as a predictor of upgrading from prdGP3 on biopsy to prdGP4 on resection specimen. Model incorporates all independent predictors of prdGP4 identified in multivariate logistic regression. Model for Logit(prdG4): X = − 3.22 + 0.04 × age + 0.04 × PSA + 0.18 × triglycerides + 0.03 × FSH − 0.08 × LH − 0.11 × BT + 0.17 × non-Caucasian ethnic group + 0.69 × presence of cardiovascular disorder, with Prob(Y = grade 4) = 1/(1 + Exp(−X))

Discussion

The impact of testosterone deficiency on PCa remains controversial, with conflicting results due in part to variations in study design/methodological limitations [6]. In the present study, we applied current recommended standards of androgen evaluation, including morning androgen measurements, centralized steroid determination by GC-MS in a single laboratory [11], and a fuller evaluation of the androgen cascade, measuring not only TT but also its precursors and metabolites and in particular the biologically active form of testosterone (BT) and androgen plasma–binding globulin (SHBG).

As most circulating testosterone is unavailable to prostatic cells, being tightly bound to SHBG, BT typically only represents 50% of TT [6]. Although prostatic tissue is sensitive to changes in serum TT at low concentrations, a saturation limit exists at higher androgen levels [18]. Arguably, BT is a better indicator of circulating androgens available to target organs, and so any possible impact of hypogonadism on PCa should be evaluated using both TT and also BT levels. However, to date, few studies have suggested a correlation between low BT and high-grade PCA. Indeed, most studies do not specifically assay BT; usually, FT is calculated from a RIA measurement of TT; then, BT values are calculated based on standard formulae. The role of BT remains infrequently reported. In their pooled analysis of 18 studies recruiting 3886 men with incident PCa, Roddam and colleagues did not report BT levels [19], although more recently, Garcia-Cruz et al. reported an almost fivefold increased risk of detection of PCa on biopsy in patients with low BT [20].

In the present study, we found that BT appeared twice in a random sample of 24 models that had good predictive quality for HGPC or ISUP ≥ 3 (i.e., prdGP4) making it a candidate for the subset of variables that were to be included in such a model; in contrast, TT was not present in any of these models. Hypogonadic patients had a higher rate of prdGP4, with a lower risk in patients with physiological TT levels (Table 2). No significant differences in adrenal (i.e., DHEA, ∆5) or testicular (i.e., ∆4, E2, E1, TT, or DHT) steroid levels in patients with different tumor grades were observed. However, in prdGP4 patients, BT and FT levels were substantially and significantly lower (Table 1), a direct consequence of higher SHBG levels in patients with prdGP4 disease.

We suggest that using both TT and BT to define hypo-gonadism is a highly relevant refinement in terms of pathology, androgen cascade, and clinical parameters; testosterone deficiency should not be defined using TT in isolation. In this study, testosterone deficiency based on TT alone was found in 11.2% (14.2% in prdGP4), while adding cases based on BT increased this percentage by 8.3% (10.7% in prdGP4) leading to a total of 19.5% cases (24.9% in prdGP4); data similar to that was reported elsewhere [7]. Our data also support the growing evidence that testosterone deficiency is independently associated with higher tumor aggressiveness. Our results have implications for patient selection for nerve-sparing surgery, as EAU guidelines highlight that predominant Gleason grade 4 jeopardizes outcomes of nerve-sparing RP [1, 2].

We found weak concordance between Gleason scores of biopsy samples and full RP specimens; 21.1% of our 1045 patients with a prdGP3 biopsy had prdGP4 cancer on surgical specimens. More surprisingly, in 75 (26.8%) of the 279 patients with prdGP4 on biopsy, no prdGP4 could be detected on the surgical specimen. Although other studies report concordance ranging between 30 and 67% [21], the moderate concordance we found was despite systematic sampling (12 core biopsies) and centralized reviewing of biopsies and prostatectomy specimens by a single reference uropathologist. Discordance between biopsy and RP specimen grading is not necessarily problematic, as long as this is considered in clinical decision-making. In this context, we and others have previously reported that Gleason scores are more often under-evaluated by prostate biopsies in hypogonadal patients [16, 22], which further supports the view that hormonal status should be fully considered in PCa case assessment.

In a previous study, we found that obesity and hypo-gonadism, while not interdependent, can coexist and together may be associated with more aggressive disease [23]. Even if obesity is not the primary cause of hypogonadism, there is a large overlap between obesity and hypogonadism, and both impact PCa status. In the present study, we found that low preoperative serum TT levels, obesity, and metabolic syndrome parameters were not independently associated with aggressive features of PCa (such as extracapsular extension) nor with tumor grading and that weight, BMI, and body fat are not predictive factors of HGPC. We also found a relationship between obesity and low TT levels (but not low BT levels). This, and the fact that BT, not TT, played a role in some of the models developed to predict prdGP4 in patients with prdGP3 on biopsy, suggests a dichotomy in the role of TT and BT; namely, TT relates to adiposity and BT to PCa aggressiveness. To our knowledge, this dichotomy has not previously been suggested. One possible explanation is that BMI and metabolic syndrome parameters are only surrogate indicators of low levels of circulating androgens in PCa patients. Indeed, metabolic syndrome has been shown to result in decreased FT and BT levels in PCa, although other biological mechanisms may influence the association of metabolic disorders and aggressive PCa [24]. Higher SHBG levels were also associated with a decrease in BT, with TT levels essentially remaining normal. This may explain why obesity and diabetes can lead to hypogonadism in aging males.

The role of SHBG as a predictor of PCa extracapsular extension in men undergoing RP has previously been reported [25], although that group subsequently reported that SHBG levels are not multivariate predictors of high-risk PCa (as defined by the NCCN Practice guidelines) [26, 27], highlighting the importance of using consistent, well-defined terminology across studies evaluating PCa. Our results provide some support for predictive benefit of SHBG, with higher levels in patients with prdGP4 tumors compared to prdGP3 patients.

Finally, while we found no statistically significant differences in E1 and E2 levels between prdGP3 and prdGP4 patient groups, the absolute levels (25.6 to 26.1 pg/mL for E2), determined by GC-MS, were lower than levels (determined by RIA) reported previously by Salonia et al., where a significant association between the rate of HGPC and E2 level ≥ 50 pg/mL was seen in a small patient subset [28]. In the present study, using GC-MS, we found relatively few patients (1%) with E2 ≥ 50 pg/mL and suggest this cutoff level is not relevant for clinical practice.

Our study has strengths and limitations. The prospective multicenter aspect allowed consecutive enrolment of a com-paratively large number of patients (none of whom had received testosterone therapy or had comorbidity-associated hypo-gonadism) with data collected in a consistent manner and little if any selection bias. All RP specimens were graded by a single experienced specialist pathologist (using the most recent criteria.) Another distinguishing feature (which may help explain the discrepancy with previous study findings) was our use of a single central laboratory with analyses by GC-MS [11] and sampling according to recent recommendations [8]. Nevertheless, limitations exist; to fully assess differences influenced by any BT-TT interaction, a very large trial is required. Although total patient numbers in this study are comparatively high, less than 20% of our patients were hypogonadal, as assessed by either TT and/or BT, and so our results should be viewed in this context.

In conclusion, we found that testosterone deficiency (defined by TT and/or BT levels) was independently associated with higher PCa aggressiveness. Furthermore, BT is a key hormonal marker, and accurate assessment of BT levels is necessary when evaluating hypogonadism in patients with localized PCa. A previously undescribed dichotomy was observed, TT having a stronger link with adiposity and BT a stronger link to cancer aggressiveness. SHBG levels should also be considered in treatment decisions.