Purpose:

Intensification of androgen deprivation therapy (ADT) with either docetaxel or androgen receptor axis–targeted therapies (ARAT) are the current standard of care for patients with metastatic castration-sensitive prostate cancer (mCSPC). However, biomarkers guiding treatment selection are lacking. We hypothesized that ADT intensification with ARAT, but not with docetaxel, would be associated with improved outcomes in patients with de novo (dn)-mCSPC harboring SPOP mutations.

Experimental Design:

Patient-level data from a deidentified nationwide (U.S.-based) prostate cancer clinico-genomic database between January 2011 and December 2021 were extracted. Eligibility criteria: diagnosis of metastatic disease within 30 days of original prostate cancer diagnosis, genomic profiling of a tissue biopsy collected within 90 days of original diagnosis, and initiation of ARAT or docetaxel within 120 days of initial diagnosis. The log-rank test and Cox proportional hazards models were used to compare time to castration-resistant prostate cancer (TTCRPC) and overall survival (OS) for patients with and without SPOP mutations undergoing ADT intensification with ARAT or docetaxel.

Results:

In the ARAT cohort, presence of SPOP mutation compared with wild-type was associated with more favorable TTCRPC [not reached (NR) vs. 16.7 months; adjusted HR (aHR), 0.20; 95% confidence interval (CI), 0.06–0.63; P = 0.006] and OS (NR vs. 27.2 months; aHR, 0.19; 95% CI, 0.05–0.79; P = 0.022). In contrast, SPOP mutation status was not associated with TTCRPC or OS in docetaxel-treated cohort.

Conclusions:

In real-world settings, SPOP mutations were associated with improved outcomes to ADT plus ARAT (but not ADT plus docetaxel) in patients with dn-mCSPC. This may serve as a predictive biomarker to guide treatment selection for patients with mCSPC.

Translational Relevance

Current guidelines recommend intensification of androgen deprivation therapy (ADT) with either docetaxel or androgen receptor axis–targeted therapies (ARAT) for patients with metastatic castration-sensitive prostate cancer (mCSPC). However, there are no biomarkers to select treatments. In this hypothesis-driven, real-world study, presence of SPOP mutation (compared with wild-type SPOP) in patients with de novo mCSPC was associated with more favorable survival outcomes (time to castration-resistant prostate cancer and overall survival) with ADT in combination with ARAT compared with docetaxel. The presence of SPOP mutation may be the first biomarker to guide treatment selection for patients with de novo mCSPC.

Over the past 7 years, the therapeutic landscape for advanced prostate cancer has expanded significantly with the introduction of several novel agents and moving of agents typically given in the castration-resistant stage to earlier settings due to significant improvements in survival (1–5). Multiple phase III randomized trials have established the benefit of intensification of androgen deprivation therapy (ADT) with the addition of docetaxel or androgen receptor axis–targeted therapies (ARAT) such as abiraterone, enzalutamide, and apalutamide in improving survival outcomes for patients with metastatic castration-sensitive prostate cancer (mCSPC), and are now the new standard of care (1–7). However, we currently lack biomarkers to identify which patients benefit most from ADT intensification with ARAT or docetaxel.

A lack of prospective clinical trials in prostate cancer focusing on pathogenic alterations other than those belonging to the homologous recombination repair pathway has been recognized as a knowledge gap (8). While prostate cancer is generally viewed as a disease driven by androgen receptor (AR) signaling, multiple reports have identified several unique molecular subtypes (9). One of these subgroups is characterized by recurrent inactivating point mutations in the SPOP gene (9–12). SPOP, “speckle type BTB/POZ protein,” is an adaptor protein of the CUL3-RBX1 E3 ubiquitin ligase complex involved in the binding of multiple oncogenic proteins for ubiquitination and subsequent proteasomal degradation (13). Several studies indicate that loss-of-function SPOP mutations are an early event associated with prostate tumorigenesis, resulting in stabilized AR protein and increased AR signaling. Prior reports suggest that SPOP mutation may be associated with improved survival outcomes in patients with mCSPC treated with ADT, or metastatic castration-resistant prostate cancer treated with an ARAT (14–16). However, the predictive potential of SPOP mutations has not been studied, nor has its ability to influence taxane therapy benefit.

Because of greater dependency of SPOP-altered prostate cancer on androgen signaling, we hypothesized that ADT intensification with ARAT, but not with docetaxel, will be associated with improved outcomes in patients with de novo (dn)-mCSPC harboring SPOP mutations.

Study design and patient selection

The study cohort comprised patients with a confirmed diagnosis of de novo mCSPC included in the U.S.-wide Flatiron Health (FH)-Foundation Medicine (FMI) deidentified mCSPC clinico-genomic database (CGDB) between January 2011 and December 2021 who underwent genomic testing using Foundation Medicine comprehensive genomic profiling (CGP) assays.

The deidentified data originated from approximately 280 U.S. cancer clinics (∼800 sites of care). Retrospective longitudinal clinical data were derived from electronic health record (EHR) data, comprising patient-level structured and unstructured data, curated via technology-enabled abstraction of clinical notes and radiology/pathology reports and were linked to CGP data by deidentified, deterministic matching (17). FH database mortality information is a composite derived from three sources: documents within the EHR, Social Security Death Index, and a commercial death dataset mining data from obituaries and funeral homes. This mortality information has been externally validated in comparison to the National Death Index with greater than 90% accuracy (18). The CGDB has previously supported replications of biomarker associations with survival observed in clinical trials (17, 19, 20). Abstracted clinical data included demographics, clinical features, therapy exposure (with start and stop dates for each therapy line), and survival.

Patient inclusion criteria were as follows: diagnosis of metastatic disease within 30 days of original prostate cancer diagnosis (to ensure de novo metastatic diagnosis), tissue biopsy collected within 90 days of original diagnosis to better reflect phenotype/genotype characteristics at the time of treatment decisions, and start of ARAT or docetaxel therapy within 120 days of initial diagnosis to improve homogeneity of the population (Fig. 1A and B for flowchart).

Figure 1.

Cohort selection. Cohort selection diagram (A) and temporal visualization of the analysis cohort (B) are shown. Dx, diagnosis; met, metastasis; TTNT, time to next treatment; Tx, treatment.

Figure 1.

Cohort selection. Cohort selection diagram (A) and temporal visualization of the analysis cohort (B) are shown. Dx, diagnosis; met, metastasis; TTNT, time to next treatment; Tx, treatment.

Close modal

Comprehensive genomic profiling

The FoundationOne or FoundationOneCDx hybrid capture-based next-generation sequencing assays were performed on tumor biopsies in a Clinical Laboratory Improvement Amendments–certified, College of American Pathologists–accredited laboratory (Foundation Medicine). Samples were evaluated for alterations as described previously (21). SPOP mutations were included if characterized as known pathogenic by existing literature, reflecting the documented hotspot alterations (missense changes, small in-frame insertions/deletions) within the MATH domain (amino acids 38–163). Patient specimens not containing a known pathogenic SPOP mutation are referred to as “SPOP wild-type.”

For larger genomic database analyses in supplementary, tissue biopsies from prostate acinar adenocarcinoma were profiled using FoundationOne (N = 4,200) and FoundationeOneCDx (N = 13,550), and liquid biopsies from patients with prostate cancer were profiled using FoundationOneLiquid CDx (N = 5,629). Approval for this study, including a waiver of informed consent and a Health Insurance Portability and Accountability Act waiver of authorization, was obtained from the WIRB-Copernicus Group Institutional Review Board (protocol no. 20152817).

Statistical analysis and outcome measures

A prespecified statistical analysis plan (SAP) was developed. Consistent with ISPOR guidelines (22), inclusion and exclusion criteria, potential biases, primary outcome measures, exploratory outcome measures, handling of missing data, and all methods described below were defined prior to analysis execution unless noted. The prespecified primary outcome analysis compared the effectiveness of ARAT, and separately docetaxel, in patients with dn-mCSPC with either mutated SPOP (mutSPOP) or wild-type SPOP (wtSPOP), evaluating both time to castration-resistant prostate cancer (TTCRPC) and overall survival (OS) in concert. Exploratory outcomes included effectiveness of ARAT, and separately docetaxel, in patients with dn-mCSPC harboring TMPRSS2-ERG gene fusions; as well as a comparison of ARAT versus docetaxel subgroup analyses adjusted for baseline prognostic features.

χ2 tests and Wilcoxon rank-sum tests were used to assess differences between cohorts of categorical and continuous variables, respectively. Missing values were imputed with the expected values based upon observed covariates using random forests with R package “missForest,” with imputed values treated identically to measured values in subsequent analysis. TTCRPC was calculated from dates of metastatic diagnosis to castration resistance, as assessed by the clinician, or last clinic visit if not yet castration resistant. Patients were right censored if not yet reached CRPC. OS was calculated from metastatic diagnosis to death from any cause, and patients with no record of mortality were right censored at the date of the last clinic visit. Because patients cannot enter the database until a CGP report is delivered, OS risk intervals were left truncated (delayed entry) to the date of the CGP report to account for immortal time (23, 24).

Differences in time-to-event outcomes were assessed with the log-rank test and Cox proportional hazards (PH) models. Multiple comparison adjustments were not performed; P values are reported to quantify the strength of association between biomarker and each respective outcome, not for null hypothesis significance testing. Results were interpreted considering outcome measures collectively (TTCRPC, OS) rather than outcome measures standing on their own.

Biomarker treatment interaction tests (Pinteraction; ref. 25) made use of multivariable Cox PH regression models containing the following variables: drug class (ARAT or docetaxel), biomarker (mutSPOP vs. wtSPOP), and the interaction term between drug class and biomarker, adjusted for risk of CRPC or death.

Pretherapy prognostic features [Eastern Cooperative Oncology Group (ECOG) score, Gleason score, PSA, age, alkaline phosphatase, albumin, hemoglobin] associated with risk of CRPC or OS were combined into multivariable models (Supplementary Table S1) and resulting sum of coefficients used as a risk score for adjustment in all Cox PH models, as described previously (26, 27). Risk adjustment was not declared in SAP.

R version 3.6.3 software was used for all statistical analyses.

Data availability

The data supporting the findings of this study originated from Flatiron Health, Inc. and Foundation Medicine, Inc. These deidentified data may be made available upon request, and are subject to a license agreement with FH and FMI; interested researchers should contact <[email protected]> and <[email protected]> to determine licensing terms.

Characteristics of analysis cohort

Of 4,417 profiled prostate cancer tissue specimens available, 1,497 were obtained within 90 days of original prostate cancer diagnosis, and were thus considered to reflect hormone-sensitive specimens (Fig. 1A). mutSPOP was detected in 152/1,497 (10.2% prevalence; Supplementary Table S2). Of these, 447 unique patients met the eligibility criteria additionally inclusive of metastatic disease diagnosed within 30 days of original diagnosis, and ARAT or docetaxel treatment beginning within 120 days of original diagnosis (Fig. 1A and B), treated with either ARAT (n = 223) or docetaxel (n = 214). mutSPOP was detected in 20/223 (8.6%) and 18/214 (8.4%) patients in the ARAT and docetaxel cohorts, respectively (Fig. 1A; Table 1). Mutations clustered in the MATH domain, with hotspots at residues F102, F133, and W131 (Fig. 2A). Compared with wtSPOP samples, SPOP-mutant samples were mutually exclusive with those harboring TMPRSS2-ERG fusions, were enriched for coalterations in Wnt pathway genes (particularly APC), and were depleted for concurrent alterations in TP53, PTEN, and RB1 (Fig. 2B). As a supplementary analysis, we sought to determine whether similar frequencies of specific MATH domain alterations detected between liquid biopsy and tissue profiling assays were similar. To do this, we made use of the broader Foundation Medicine database (see Materials and Methods) without clinical or outcome associations. Evaluating 17,750 tissue specimens and 5,629 liquid biopsy specimens, the most common hotspots in descending order for tissue profiling and liquid profiling were both F133, F102, and W131 (Supplementary Fig. S1), also mirroring the outcome analysis cohort.

Table 1.

Clinical features of patients by treatment assignment.

ARAT (N = 233)Taxane (N = 214)Total (N = 447)P
SPOP status    0.948 
 Mutation 20 (8.6%) 18 (8.4%) 38 (8.5%)  
 Wild-type 213 (91.4%) 196 (91.6%) 409 (91.5%)  
Months follow-up of OS censored    0.255 
 Median (Q1, Q3) 26.0 (17.2, 33.4) 20.7 (12.5, 34.4) 21.2 (12.8, 34.4)  
Year of Tx initiation    <0.001 
 2014 0 (0.0%) 8 (3.7%) 8 (1.8%)  
 2015 0 (0.0%) 20 (9.3%) 20 (4.5%)  
 2016 0 (0.0%) 27 (12.6%) 27 (6.0%)  
 2017 11 (4.7%) 37 (17.3%) 48 (10.7%)  
 2018 27 (11.6%) 34 (15.9%) 61 (13.6%)  
 2019 63 (27.0%) 35 (16.4%) 98 (21.9%)  
 2020 65 (27.9%) 31 (14.5%) 96 (21.5%)  
 2021 60 (25.8%) 20 (9.3%) 80 (17.9%)  
 2022 7 (3.0%) 2 (0.9%) 9 (2.0%)  
Age    <0.001 
 Median (Q1, Q3) 69.0 (62.0, 77.0) 64.5 (57.0, 70.0) 67.0 (60.0, 74.0)  
ECOG    0.128 
 0 85 (50.3%) 112 (57.7%) 197 (54.3%)  
 1 63 (37.3%) 71 (36.6%) 134 (36.9%)  
 2 18 (10.7%) 9 (4.6%) 27 (7.4%)  
 3 3 (1.8%) 2 (1.0%) 5 (1.4%)  
 Missing observations 64 20 84  
Gleason    0.772 
 6 or less 1 (0.4%) 0 (0.0%) 1 (0.2%)  
 7–8 39 (16.7%) 38 (17.8%) 77 (17.2%)  
 9–10 136 (58.4%) 127 (59.3%) 263 (58.8%)  
 Unknown/Not documented 57 (24.5%) 49 (22.9%) 106 (23.7%)  
PSA at Met Dx    0.475 
 Median (Q1, Q3) 135.0 (34.3, 445.9) 146.5 (37.9, 631.5) 139.1 (35.6, 564.6)  
 Missing observations 10 16  
Alkaline phosphatase    0.024 
 Above ULN 105 (52.5%) 123 (63.7%) 228 (58.0%)  
 Normal 95 (47.5%) 70 (36.3%) 165 (42.0%)  
 Missing observations 33 21 54  
Albumin    0.568 
 Below LLN 15 (7.2%) 12 (5.8%) 27 (6.5%)  
 Normal 193 (92.8%) 194 (94.2%) 387 (93.5%)  
 Missing observations 25 33  
Hemoglobin    0.389 
 Below LLN 113 (53.3%) 119 (57.5%) 232 (55.4%)  
 Normal 99 (46.7%) 88 (42.5%) 187 (44.6%)  
 Missing observations 21 28  
Practice type    0.603 
 Academic 35 (15.0%) 36 (16.8%) 71 (15.9%)  
 Community 198 (85.0%) 178 (83.2%) 376 (84.1%)  
Treatment received    <0.001 
 Abiraterone 143 (61.4%) 0 (0.0%) 143 (32.0%)  
 Apalutamide 19 (8.2%) 0 (0.0%) 19 (4.3%)  
 Cabazitaxel 0 (0.0%) 1 (0.5%) 1 (0.2%)  
 Docetaxel 0 (0.0%) 213 (99.5%) 213 (47.7%)  
 Enzalutamide 71 (30.5%) 0 (0.0%) 71 (15.9%)  
Biopsy site    0.077 
 Bladder 4 (1.7%) 6 (2.8%) 10 (2.2%)  
 Bone 30 (12.9%) 16 (7.5%) 46 (10.3%)  
 Liver 3 (1.3%) 6 (2.8%) 9 (2.0%)  
 Lymph node 25 (10.7%) 28 (13.1%) 53 (11.9%)  
 Other 11 (4.7%) 21 (9.8%) 32 (7.2%)  
 Prostate 160 (68.7%) 137 (64.0%) 297 (66.4%)  
Reached CRPC    <0.001 
 No 125 (53.6%) 54 (25.2%) 179 (40.0%)  
 Yes 108 (46.4%) 160 (74.8%) 268 (60.0%)  
Death events    <0.001 
 Alive 163 (70.0%) 101 (47.2%) 264 (59.1%)  
 Deceased 70 (30.0%) 113 (52.8%) 183 (40.9%)  
ARAT (N = 233)Taxane (N = 214)Total (N = 447)P
SPOP status    0.948 
 Mutation 20 (8.6%) 18 (8.4%) 38 (8.5%)  
 Wild-type 213 (91.4%) 196 (91.6%) 409 (91.5%)  
Months follow-up of OS censored    0.255 
 Median (Q1, Q3) 26.0 (17.2, 33.4) 20.7 (12.5, 34.4) 21.2 (12.8, 34.4)  
Year of Tx initiation    <0.001 
 2014 0 (0.0%) 8 (3.7%) 8 (1.8%)  
 2015 0 (0.0%) 20 (9.3%) 20 (4.5%)  
 2016 0 (0.0%) 27 (12.6%) 27 (6.0%)  
 2017 11 (4.7%) 37 (17.3%) 48 (10.7%)  
 2018 27 (11.6%) 34 (15.9%) 61 (13.6%)  
 2019 63 (27.0%) 35 (16.4%) 98 (21.9%)  
 2020 65 (27.9%) 31 (14.5%) 96 (21.5%)  
 2021 60 (25.8%) 20 (9.3%) 80 (17.9%)  
 2022 7 (3.0%) 2 (0.9%) 9 (2.0%)  
Age    <0.001 
 Median (Q1, Q3) 69.0 (62.0, 77.0) 64.5 (57.0, 70.0) 67.0 (60.0, 74.0)  
ECOG    0.128 
 0 85 (50.3%) 112 (57.7%) 197 (54.3%)  
 1 63 (37.3%) 71 (36.6%) 134 (36.9%)  
 2 18 (10.7%) 9 (4.6%) 27 (7.4%)  
 3 3 (1.8%) 2 (1.0%) 5 (1.4%)  
 Missing observations 64 20 84  
Gleason    0.772 
 6 or less 1 (0.4%) 0 (0.0%) 1 (0.2%)  
 7–8 39 (16.7%) 38 (17.8%) 77 (17.2%)  
 9–10 136 (58.4%) 127 (59.3%) 263 (58.8%)  
 Unknown/Not documented 57 (24.5%) 49 (22.9%) 106 (23.7%)  
PSA at Met Dx    0.475 
 Median (Q1, Q3) 135.0 (34.3, 445.9) 146.5 (37.9, 631.5) 139.1 (35.6, 564.6)  
 Missing observations 10 16  
Alkaline phosphatase    0.024 
 Above ULN 105 (52.5%) 123 (63.7%) 228 (58.0%)  
 Normal 95 (47.5%) 70 (36.3%) 165 (42.0%)  
 Missing observations 33 21 54  
Albumin    0.568 
 Below LLN 15 (7.2%) 12 (5.8%) 27 (6.5%)  
 Normal 193 (92.8%) 194 (94.2%) 387 (93.5%)  
 Missing observations 25 33  
Hemoglobin    0.389 
 Below LLN 113 (53.3%) 119 (57.5%) 232 (55.4%)  
 Normal 99 (46.7%) 88 (42.5%) 187 (44.6%)  
 Missing observations 21 28  
Practice type    0.603 
 Academic 35 (15.0%) 36 (16.8%) 71 (15.9%)  
 Community 198 (85.0%) 178 (83.2%) 376 (84.1%)  
Treatment received    <0.001 
 Abiraterone 143 (61.4%) 0 (0.0%) 143 (32.0%)  
 Apalutamide 19 (8.2%) 0 (0.0%) 19 (4.3%)  
 Cabazitaxel 0 (0.0%) 1 (0.5%) 1 (0.2%)  
 Docetaxel 0 (0.0%) 213 (99.5%) 213 (47.7%)  
 Enzalutamide 71 (30.5%) 0 (0.0%) 71 (15.9%)  
Biopsy site    0.077 
 Bladder 4 (1.7%) 6 (2.8%) 10 (2.2%)  
 Bone 30 (12.9%) 16 (7.5%) 46 (10.3%)  
 Liver 3 (1.3%) 6 (2.8%) 9 (2.0%)  
 Lymph node 25 (10.7%) 28 (13.1%) 53 (11.9%)  
 Other 11 (4.7%) 21 (9.8%) 32 (7.2%)  
 Prostate 160 (68.7%) 137 (64.0%) 297 (66.4%)  
Reached CRPC    <0.001 
 No 125 (53.6%) 54 (25.2%) 179 (40.0%)  
 Yes 108 (46.4%) 160 (74.8%) 268 (60.0%)  
Death events    <0.001 
 Alive 163 (70.0%) 101 (47.2%) 264 (59.1%)  
 Deceased 70 (30.0%) 113 (52.8%) 183 (40.9%)  

Abbreviations: LLN, lower limit of normal; ULN, upper limit of normal.

Figure 2.

SPOP mutation genomic co-occurrence in mCSPC. The location and frequency of mutations in the SPOP gene included in this analysis are shown in lollipop plot (A). The prevalence of co-occurring alterations for SPOP mutations is shown in B.

Figure 2.

SPOP mutation genomic co-occurrence in mCSPC. The location and frequency of mutations in the SPOP gene included in this analysis are shown in lollipop plot (A). The prevalence of co-occurring alterations for SPOP mutations is shown in B.

Close modal

Baseline characteristics of both treatment cohorts (ADT plus ARAT, and ADT plus docetaxel) are presented in Table 1. Patients from both cohorts had a similar ECOG, Gleason score, and PSA at metastatic diagnosis. The subgroup of mutSPOP only had similar clinical associations (Supplementary Table S3). It is important to note that docetaxel was approved prior to ARAT in the mCSPC space, with use by year reflected in Table 1, resulting in longer follow-up for the docetaxel cohort. At the time of data analysis, 46.4% of patients in the ARAT and 74.8% patients in the docetaxel cohort had progressed to the castration-resistant stage. In addition, 70.0% and 47.2% patients were alive in the ARAT and docetaxel cohorts, respectively. In the ARAT cohort, abiraterone (61.4%) was the most frequently used intensification agent followed by enzalutamide (30.5%).

Baseline characteristics of patients with mutSPOP and wtSPOP are presented in Supplementary Table S4. Patients with mutSPOP had a significantly higher median PSA than patients with wtSPOP (232.5 vs. 134.5 ng/mL; P = 0.033). Otherwise, no significant differences were observed with regards to year of treatment initiation, age, ECOG, Gleason sum, alkaline phosphatase, or hemoglobin between both groups.

Patients with SPOP mutations have more favorable TTCRPC and OS on ARAT but not docetaxel

mutSPOP in the ARAT cohort was associated with significantly better TTCRPC compared with wtSPOP [median not reached (NR) vs. 16.7 months, adjusted HR 0.20; 95% confidence interval (CI), 0.06–0.63; P = 0.006; Fig. 3A]. In contrast, more favorable TTCRPC was not observed for patients with mutSPOP versus wtSPOP in the docetaxel cohort [median 18.5 vs. 15.7 months, adjusted HR 0.86 (95% CI, 0.46–1.58; P = 0.62); Fig. 3B].

Figure 3.

Patients with SPOP mutations have more favorable TTCRPC and OS on ARAT but not docetaxel. Outcomes for SPOP mutant versus wt are shown for TTCRPC ARAT (A) and docetaxel (B) and for OS ARAT (C) and docetaxel (D). OS estimates are left truncated to reflect delayed entry to at-risk table (see Materials and Methods).

Figure 3.

Patients with SPOP mutations have more favorable TTCRPC and OS on ARAT but not docetaxel. Outcomes for SPOP mutant versus wt are shown for TTCRPC ARAT (A) and docetaxel (B) and for OS ARAT (C) and docetaxel (D). OS estimates are left truncated to reflect delayed entry to at-risk table (see Materials and Methods).

Close modal

Similarly, patients with mutSPOP in the ARAT cohort had significantly better OS compared with wtSPOP, [median NR vs. 27.2 months, adjusted HR 0.19 (95% CI, 0.05–0.79; P = 0.022); Fig. 3C]. Again, in the docetaxel cohort, OS was similar between patients with either mutSPOP or wtSPOP [median 30.4 vs. 32.3 months, adjusted HR 1.18 (95% CI, 0.57–2.44; P = 0.66); Fig. 3D].

ARAT versus docetaxel subgroup analyses, adjusted for prognostic risk

Treatment interactions were observed between SPOP mutation status (mutant vs. wild-type) and treatment type (ARAT vs. docetaxel) for both TTCRPC (Pinteraction = 0.0181) and OS (Pinteraction = 0.0279; Fig. 4A and B; Supplementary Table S5). Evaluating subgroups of the cohort, patients with mutSPOP had enrichments for outcomes favoring ARAT versus docetaxel for both TTCRPC and OS (Fig. 4A and B).

Figure 4.

ARAT versus docetaxel subgroup analyses, adjusted for prognostic risk. Relative ARAT versus docetaxel outcome analyses are shown for time to CRPC (A) and OS (B). Forest plots contain adjustments for ECOG, Gleason score, PSA, age, alkaline phosphatase, albumin, and hemoglobin.

Figure 4.

ARAT versus docetaxel subgroup analyses, adjusted for prognostic risk. Relative ARAT versus docetaxel outcome analyses are shown for time to CRPC (A) and OS (B). Forest plots contain adjustments for ECOG, Gleason score, PSA, age, alkaline phosphatase, albumin, and hemoglobin.

Close modal

TMPRSS2-ERG gene fusions characterize a separate subtype of prostate cancer and are detected in approximately 30%–50% of patients (9). Because TMPRSS2-ERG fusions are generally mutually exclusive with mutSPOP (9, 11, 28) as observed in this study too (Fig. 2B), we conducted an exploratory analysis in ARAT and docetaxel cohorts as a separate control biomarker. TMPRSS2-ERG fusions were associated with more favorable TTCRPC on ARAT (adjusted HR 0.55; 95% CI, 0.34–0.87; P = 0.011) as well as docetaxel (adjusted HR 0.60; 95% CI, 0.42–0.84; P = 0.003) compared with the absence of such a fusion (Supplementary Fig. S2A and S2B). Directionally consistent but weaker associations were seen with OS (Supplementary Fig. S2C and S2D). In contrast to SPOP, tests for treatment interactions with TMPRSS2-ERG status were not significant.

Our results suggest that mutSPOP is associated with enriched benefit from ARAT but not docetaxel in the castration-sensitive setting. These results are consistent with our hypothesis, and indicate that SPOP mutations may be the first predictive biomarker for prolonged benefit with ARAT (relative to docetaxel) in patients with dn-mCSPC.

We also observed a statistically significant higher baseline PSA in patients with mutSPOP as compared with patients with wtSPOP. In a prior study evaluating 8,559 patients undergoing radical prostatectomy, mutSPOP was associated with a significantly higher preoperative PSA. This finding was observed consistently across multiple cohorts (29). In normal cells, AR and its coactivators BRD4 and TRIM24 are targeted by SPOP for ubiquitination and proteasomal degradation. In contrast, levels of AR protein and its coactivators increase in prostate cancer cells harboring loss-of-function SPOP mutations due to enhanced protein stability (30, 31). In addition, tumors harboring SPOP and FOXA1 mutations have the highest levels of AR-induced transcripts among various molecularly distinct prostate cancer subtypes (9, 10). SPOP mutations have also been reported to result in increased degradation of 17βHSD4, an essential enzyme in steroid metabolism, resulting in increased testosterone/dihydrotesterone production (32). Taken together, these results further emphasize the increased addiction to AR signaling by prostate cancers harboring pathogenic mutations in the SPOP gene. Therefore, it is perhaps not surprising that deeper blockade of AR signaling by the combination of ADT and ARAT is significantly more effective in patients with mutSPOP-positive dn-mCSPC.

We did not observe differential benefit for mutSPOP-positive patients treated with docetaxel compared with wtSPOP. It is possible this may be due to upregulation of stress granules assembly. Stress granule assembly is a survival mechanism employed by cells under stress, and is upregulated in cancer cells. Caprin1 plays a nucleating role in stress granule assembly and is a target for wtSPOP-mediated ubiquitination and degradation. Prostate cancer cell lines harboring mutSPOP exhibit Caprin1 upregulation, increased stress granule assembly and were associated with increased resistance to docetaxel in one study (33). Alternatively, mutSPOP was recently reported to cause aberrant upregulation of global DNA methylation in cell lines and primary prostate cancer specimens (34). Accordingly, whole methylome analysis revealed hypermethylation and silencing of several tumor suppressor genes such as FOXO3, which may lead to decreased efficacy of docetaxel.

While there are no predictive biomarkers approved to guide treatment selection for patients with mCSPC, our findings indicate that patients with dn-mCSPC harboring mutSPOP will likely derive greater benefit from ARAT and not chemotherapy with docetaxel. Furthermore, there is an urgent need to validate these results in larger databases or completed randomized phase III trials to accelerate translation of mutSPOP as a predictive biomarker in this clinical setting.

The limitations of our study include a nonrandomized retrospective analysis. While biases were considered and known imbalances adjusted, unknown imbalances such as access to subsequent therapies and neuroendocrine differentiation likely remain. Quantified disease burden was not available for adjustment, and while PSA level and ECOG were included in adjustments, along with other factors associating with disease burden and severity, these are incomplete proxies for disease severity and burden. Disease volume has been evaluated as a potential biomarker for docetaxel efficacy: prespecified subgroup analyses of the CHAARTED study reported lack of observed benefit from docetaxel in patients with low volume disease (35). However, a subsequent retrospective analysis of the STAMPEDE trial showed no evidence that docetaxel benefit differed on the basis of metastatic burden (36). It is possible that some patients in the docetaxel cohort in our study had low volume disease and therefore may have derived diminished benefit. On the contrary, it is also possible that patients receiving docetaxel in routine practice settings might be more likely to have aggressive disease at presentation, which might have led to poorer outcomes in the docetaxel cohort. In this context, evaluating docetaxel only cohort for TTCRPC and OS based on SPOP status (Fig. 3,B and D) and ARAT versus docetaxel subgroup analyses after adjusting for available prognostic variables such as ECOG, Gleason score, PSA, age, alkaline phosphatase, albumin, and hemoglobin may mitigate disease severity and burden as a confounder. Co-occurring alterations such as Wnt pathway alterations, TP53, PTEN, and RB1 may also affect treatment responses, and our study did not have the power to examine clinical outcomes according to these coalterations (28). An alternate hypothesis is that CHD1 deletion, which is not assessed by the assays utilized but is highly correlated with SPOP mutations, is the more proximal cause for the increased sensitivity to hormonal therapies (16). Finally, an unanswered question is whether SPOP mutation status can help clinicians decide which patients should receive triplet therapy (ADT plus ARAT plus docetaxel) compared with doublet therapy (ADT plus ARAT) in the dn-mCSPC setting, especially as multi-agent regimens begin to enter the armamentarium.

Additional prospective-retrospective analyses of outcomes associated with mutSPOP in completed randomized trials will have reduced possibility of unadjusted confounding, improving the precision of anticipated ARAT versus docetaxel benefit in this patient population, and allow for evaluation of subgroups that additionally received both docetaxel and ARAT.

In summary, in this real-world study of patients with dn-mCSPC, the addition of ARAT (but not docetaxel) was associated with significantly longer TTCRPC and OS in patients harboring mutSPOP compared with wtSPOP. These results indicate that mutSPOP likely represents a treatment-selection biomarker in the mCSPC setting.

U. Swami reports personal fees from Astellas, Exelixis, and Seattle Genetics, as well as grants from Janssen, Exelixis, and Astellas/Seattle Genetics outside the submitted work. R.P. Graf reports personal fees from Foundation Medicine, Inc during the conduct of the study, as well as personal fees from Epic Sciences outside the submitted work; in addition, R.P. Graf has a patent pending. V. Fisher reports personal fees from Roche outside the submitted work. H. Tukachinsky reports personal fees from Foundation Medicine Inc. and other support from Roche during the conduct of the study. A.B. Schrock reports personal fees from Foundation Medicine and Roche during the conduct of the study. G. Li reports personal fees from Foundation Medicine and other support from F. Hoffmann-La Roche outside the submitted work. J.S. Ross reports personal fees and other support from Foundation Medicine outside the submitted work. G.R. Oxnard reports personal fees from Foundation Medicine and Roche during the conduct of the study. N. Agarwal reports (lifetime disclosures) consultancy to Astellas, AstraZeneca, Aveo, Bayer, Bristol Myers Squibb, Calithera, Clovis, Eisai, Eli Lilly, EMD Serono, Exelixis, Foundation Medicine, Genentech, Gilead, Janssen, Merck, MEI Pharma, Nektar, Novartis, Pfizer, Pharmacyclics, and Seattle Genetics, as well as research funding to institution from Astellas, AstraZeneca, Bavarian Nordic, Bayer, Bristol Myers Squibb, Calithera, Celldex, Clovis, Eisai, Eli Lilly, EMD Serono, Exelixis, Genentech, Gilead, GlaxoSmithKline, Immunomedics, Janssen, Medivation, Merck, Nektar, New Link Genetics, Novartis, Pfizer, Prometheus, Rexahn, Roche, Sanofi, Seattle Genetics, Takeda, and Tracon. No disclosures were reported by the other authors.

Foundation Medicine, a wholly owned subsidiary of Roche, is a for-profit company and producer of FDA-regulated molecular diagnostics. Authors employed by Foundation Medicine were involved in the design and conduct of the study, analysis, interpretation of the data, preparation, review, and approval of the article.

U. Swami: Conceptualization, writing–original draft, writing–review and editing, final approval of manuscript. R.P. Graf: Conceptualization, formal analysis, writing–original draft, writing–review and editing, final approval of manuscript. R.H. Nussenzveig: Writing–review and editing, final approval of manuscript. V. Fisher: Formal analysis, writing–review and editing, final approval of manuscript. H. Tukachinsky: Writing–review and editing, final approval of manuscript. A.B. Schrock: Writing–review and editing, final approval of manuscript. G. Li: Writing–review and editing, final approval of manuscript. J.S. Ross: Writing–review and editing, final approval of manuscript. N. Sayegh: Writing–review and editing, final approval of manuscript. N. Tripathi: Writing–review and editing, final approval of manuscript. V. Mathew Thomas: Writing–review and editing, final approval of manuscript. G.R. Oxnard: Administrative support, writing–review and editing, final approval of manuscript. E.S. Antonarakis: Conceptualization, formal analysis, writing–review and editing, final approval of manuscript. N. Agarwal: Conceptualization, formal analysis, writing–review and editing, final approval of manuscript.

We thank the patients whose data made this research possible, the clinical and laboratory staff at Foundation Medicine, and the team at Flatiron Health.

This study was supported by Foundation Medicine.

The publication costs of this article were defrayed in part by the payment of publication fees. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).

1.
Chi
KN
,
Agarwal
N
,
Bjartell
A
,
Chung
BH
,
Pereira de Santana Gomes
AJ
,
Given
R
, et al
.
Apalutamide for metastatic, castration-sensitive prostate cancer
.
N Engl J Med
2019
;
381
:
13
24
.
2.
Davis
ID
,
Martin
AJ
,
Stockler
MR
,
Begbie
S
,
Chi
KN
,
Chowdhury
S
, et al
.
Enzalutamide with standard first-line therapy in metastatic prostate cancer
.
N Engl J Med
2019
;
381
:
121
31
.
3.
Fizazi
K
,
Tran
N
,
Fein
L
,
Matsubara
N
,
Rodriguez-Antolin
A
,
Alekseev
BY
, et al
.
Abiraterone plus prednisone in metastatic, castration-sensitive prostate cancer
.
N Engl J Med
2017
;
377
:
352
60
.
4.
James
ND
,
de Bono
JS
,
Spears
MR
,
Clarke
NW
,
Mason
MD
,
Dearnaley
DP
, et al
.
Abiraterone for prostate cancer not previously treated with hormone therapy
.
N Engl J Med
2017
;
377
:
338
51
.
5.
Sweeney
CJ
,
Chen
Y-H
,
Carducci
M
,
Liu
G
,
Jarrard
DF
,
Eisenberger
M
, et al
.
Chemohormonal therapy in metastatic hormone-sensitive prostate cancer
.
N Engl J Med
2015
;
373
:
737
46
.
6.
Smith
MR
,
Hussain
M
,
Saad
F
,
Fizazi
K
,
Sternberg
CN
,
Crawford
ED
, et al
.
Darolutamide and survival in metastatic, hormone-sensitive prostate cancer
.
N Engl J Med
2022
;
386
:
1132
42
.
7.
Fizazi
K
,
Carles Galceran
J
,
Foulon
S
,
Roubaud
G
,
McDermott
R
,
Fléchon
A
, et al
.
LBA5 A phase III trial with a 2×2 factorial design in men with de novo metastatic castration-sensitive prostate cancer: overall survival with abiraterone acetate plus prednisone in PEACE-1
.
Ann Oncol
2021
;
32
:
S1299
.
8.
Vogl
UM
,
Beer
TM
,
Davis
ID
,
Shore
ND
,
Sweeney
CJ
,
Ost
P
, et al
.
Lack of consensus identifies important areas for future clinical research: Advanced Prostate Cancer Consensus Conference (APCCC) 2019 findings
.
Eur J Cancer
2022
;
160
:
24
60
.
9.
Cancer Genome Atlas Research Network
.
The molecular taxonomy of primary prostate cancer
.
Cell
2015
;
163
:
1011
25
.
10.
Barbieri
CE
,
Baca
SC
,
Lawrence
MS
,
Demichelis
F
,
Blattner
M
,
Theurillat
J-P
, et al
.
Exome sequencing identifies recurrent SPOP, FOXA1 and MED12 mutations in prostate cancer
.
Nat Genet
2012
;
44
:
685
9
.
11.
Chung
JH
,
Dewal
N
,
Sokol
E
,
Mathew
P
,
Whitehead
R
,
Millis
SZ
, et al
.
Prospective comprehensive genomic profiling of primary and metastatic prostate tumors
.
JCO Precis Oncol
2019
;
3
:
PO.18.00283
.
12.
Nguyen
B
,
Fong
C
,
Luthra
A
,
Smith
SA
,
DiNatale
RG
,
Nandakumar
S
, et al
.
Genomic characterization of metastatic patterns from prospective clinical sequencing of 25,000 patients
.
Cell
2022
;
185
:
563
75
.
13.
Wang
Z
,
Song
Y
,
Ye
M
,
Dai
X
,
Zhu
X
,
Wei
W
.
The diverse roles of SPOP in prostate cancer and kidney cancer
.
Nat Rev Urol
2020
;
17
:
339
50
.
14.
Stopsack
KH
,
Nandakumar
S
,
Wibmer
AG
,
Haywood
S
,
Weg
ES
,
Barnett
ES
, et al
.
Oncogenic genomic alterations, clinical phenotypes, and outcomes in metastatic castration-sensitive prostate cancer
.
Clin Cancer Res
2020
;
26
:
3230
8
.
15.
Swami
U
,
Isaacsson Velho
P
,
Nussenzveig
R
,
Chipman
J
,
Sacristan Santos
V
,
Erickson
S
, et al
.
Association of SPOP mutations with outcomes in men with de novo metastatic castration-sensitive prostate cancer
.
Eur Urol
2020
;
78
:
652
6
.
16.
Boysen
G
,
Rodrigues
DN
,
Rescigno
P
,
Seed
G
,
Dolling
D
,
Riisnaes
R
, et al
.
SPOP-mutated/CHD1-deleted lethal prostate cancer and abiraterone sensitivity
.
Clin Cancer Res
2018
;
24
:
5585
93
.
17.
Singal
G
,
Miller
PG
,
Agarwala
V
,
Li
G
,
Kaushik
G
,
Backenroth
D
, et al
.
Association of patient characteristics and tumor genomics with clinical outcomes among patients with non-small cell lung cancer using a clinicogenomic database
.
JAMA
2019
;
321
:
1391
9
.
18.
Zhang
Q
,
Gossai
A
,
Monroe
S
,
Nussbaum
NC
,
Parrinello
CM
.
Validation analysis of a composite real-world mortality endpoint for patients with cancer in the United States
.
Health Serv Res
2021
;
56
:
1281
7
.
19.
Graf
RP
,
Fisher
V
,
Huang
RPS
,
Hamdani
O
,
Gjoerup
OV
,
Stanke
J
, et al
.
Tumor mutational burden as a predictor of first-line immune checkpoint inhibitor versus carboplatin benefit in cisplatin-unfit patients with urothelial carcinoma
.
JCO Precis Oncol
2022
;
e2200121
.
20.
Graf
RP
,
Fisher
V
,
Mateo
J
,
Gjoerup
OV
,
Madison
RW
,
Raskina
K
, et al
.
Predictive genomic biomarkers of hormonal therapy versus chemotherapy benefit in metastatic castration-resistant prostate cancer
.
Eur Urol
2022
;
81
:
37
47
.
21.
Frampton
GM
,
Fichtenholtz
A
,
Otto
GA
,
Wang
K
,
Downing
SR
,
He
J
, et al
.
Development and validation of a clinical cancer genomic profiling test based on massively parallel DNA sequencing
.
Nat Biotechnol
2013
;
31
:
1023
31
.
22.
Berger
ML
,
Mamdani
M
,
Atkins
D
,
Johnson
ML
.
Good research practices for comparative effectiveness research: defining, reporting and interpreting nonrandomized studies of treatment effects using secondary data sources: the ISPOR Good Research Practices for Retrospective Database Analysis Task Force Report–Part I
.
Value Health
2009
;
12
:
1044
52
.
23.
McGough
SF
,
Incerti
D
,
Lyalina
S
,
Copping
R
,
Narasimhan
B
,
Tibshirani
R
.
Penalized regression for left-truncated and right-censored survival data
.
Stat Med
2021
;
40
:
5487
500
.
24.
Brown
S
,
Lavery
JA
,
Shen
R
,
Martin
AS
,
Kehl
KL
,
Sweeney
SM
, et al
.
Implications of selection bias due to delayed study entry in clinical genomic studies
.
JAMA Oncol
2022
;
8
:
287
91
.
25.
Ballman
KV
.
Biomarker: predictive or prognostic?
J Clin Oncol
2015
;
33
:
3968
71
.
26.
Arbogast
PG
,
Ray
WA
.
Performance of disease risk scores, propensity scores, and traditional multivariable outcome regression in the presence of multiple confounders
.
Am J Epidemiol
2011
;
174
:
613
20
.
27.
Graf
RP
,
Hullings
M
,
Barnett
ES
,
Carbone
E
,
Dittamore
R
,
Scher
HI
.
Clinical utility of the nuclear-localized AR-V7 biomarker in circulating tumor cells in improving physician treatment choice in castration-resistant prostate cancer
.
Eur Urol
2020
;
77
:
170
7
.
28.
Nakazawa
M
,
Fang
M
,
Marshall
H.C
,
Lotan
TL
,
Isaacsson Velho
P
,
Antonarakis
ES
.
Clinical and genomic features of SPOP-mutant prostate cancer
.
Prostate
2022
;
82
:
260
8
.
29.
Liu
D
,
Takhar
M
,
Alshalalfa
M
,
Erho
N
,
Shoag
J
,
Jenkins
RB
, et al
.
Impact of the SPOP mutant subtype on the interpretation of clinical parameters in prostate cancer
.
JCO Precis Oncol
2018
;
2018
:
PO.18.00036
.
30.
Janouskova
H
,
El Tekle
G
,
Bellini
E
,
Udeshi
ND
,
Rinaldi
A
,
Ulbricht
A
, et al
.
Opposing effects of cancer-type-specific SPOP mutants on BET protein degradation and sensitivity to BET inhibitors
.
Nat Med
2017
;
23
:
1046
54
.
31.
Theurillat
J-PP
,
Udeshi
ND
,
Errington
WJ
,
Svinkina
T
,
Baca
SC
,
Pop
M
, et al
.
Prostate cancer. Ubiquitylome analysis identifies dysregulation of effector substrates in SPOP-mutant prostate cancer
.
Science
2014
;
346
:
85
9
.
32.
Shi
L
,
Yan
Y
,
He
Y
,
Yan
B
,
Pan
Y
,
Orme
JJ
, et al
.
Mutated SPOP E3 ligase promotes 17βHSD4 protein degradation to drive androgenesis and prostate cancer progression
.
Cancer Res
2021
;
81
:
3593
606
.
33.
Shi
Q
,
Zhu
Y
,
Ma
J
,
Chang
K
,
Ding
D
,
Bai
Y
, et al
.
Prostate cancer-associated SPOP mutations enhance cancer cell survival and docetaxel resistance by upregulating Caprin1-dependent stress granule assembly
.
Mol Cancer
2019
;
18
:
170
.
34.
Zhang
J
,
Gao
K
,
Xie
H
,
Wang
D
,
Zhang
P
,
Wei
T
, et al
.
SPOP mutation induces DNA methylation via stabilizing GLP/G9a
.
Nat Commun
2021
;
12
:
5716
.
35.
Kyriakopoulos
CE
,
Chen
Y-H
,
Carducci
MA
,
Liu
G
,
Jarrard
DF
,
Hahn
NM
, et al
.
Chemohormonal therapy in metastatic hormone-sensitive prostate cancer: long-term survival analysis of the randomized phase III E3805 CHAARTED trial
.
J Clin Oncol
2018
;
36
:
1080
7
.
36.
Clarke
NW
,
Ali
A
,
Ingleby
FC
,
Hoyle
A
,
Amos
CL
,
Attard
G
, et al
.
Addition of docetaxel to hormonal therapy in low- and high-burden metastatic hormone sensitive prostate cancer: long-term survival results from the STAMPEDE trial
.
Ann Oncol
2019
;
30
:
1992
2003
.

Supplementary data