Introduction

Prostate-specific membrane antigen (PSMA), a type II transmembrane protein encoded by the FOLH1 gene, is highly expressed in prostate cancer (PCa) and has become a widely used target for molecular theranostics of PCa. The advent of PSMA-PET imaging and 177Lu-PSMA-617 radioligand therapy have profoundly reshaped the landscape of clinical practice for the management of lethal metastatic castration-resistant PCa (mCRPC)1,2. However, a subset of patients derive limited benefit from PSMA-targeted therapies, presenting a significant clinical challenge3,4. This may be in part attributed to the heterogeneous or absent expression of PSMA leading to suboptimal responses to PSMA-targeted strategy in some cases5. The heterogeneity of PSMA expression in both primary hormone-sensitive prostate cancer (HSPC) and mCRPC remains insufficiently explored in the current literature.

Several studies have reported membranous-cytoplasmic co-staining of PSMA in immunohistochemistry (IHC), with some research suggesting that membranous expression is the clinically actionable target, while cytoplasmic expression is considered as a potential interference5. In contrast, other studies propose that cytoplasmic expression reflects the internalization activity of membrane receptors6and integrating both membranous and cytoplasmic expression may offer a better prediction of prognosis and therapeutic efficacy7.

Many studies relied on nuclear medicine techniques to evaluate PSMA expression8which, while informative and mainstream, do not provide detailed information on cellular or subcellular localization of PSMA expression. In comparison, IHC offers both qualitative and subcellular-quantitative data, and is commonly applied in clinical practice. While discrepancies between nuclear imaging and IHC results have been reported9low PSMA expression in IHC often corresponds to limited uptake of PSMA-targeted radiopharmaceuticals5.

Precision treatment strategies targeting cell surface tumor-associated antigens (TAAs) have become prominent approaches in oncology, with antibody-drug conjugates (ADCs) leading the field in development. ADCs combine specific antibodies with cytotoxic agents10,11. The response to ADCs can be influenced by various factors, including antibody specificity, cytotoxic payload, linker chemistry, and conjugation methods, with the expression level of TAAs being of paramount importance12,13,14. Although controversial, higher TAA expression levels are generally believed to enhance ADC efficacy theoretically11. As of early 2024, fifteen ADCs targeting eleven TAAs have been approved for both solid and hematologic tumors. Among these, four targets—human epidermal growth factor receptor 2 (HER2), nectin cell adhesion molecule 4 (NECTIN4), tumor-associated calcium signal transducer 2 (TROP2), and tissue factor (TF)—have been FDA-approved for use in solid tumors11. In the field of PCa, relatively well-established ADC targets further include six-transmembrane epithelial antigen of the prostate 1 (STEAP1) and B7-homolog 3/CD276 (B7H3)15,16.

The aim of this study was to systematically characterize the heterogeneity of PSMA expression in both HSPC and mCRPC, and to head-to-head explore the expression profiles of these alternative well-established TAAs in PSMA-negative cases.

Results

A total of 127 primary lesions (PL) from HSPC cohort and 76 bone metastatic lesions (BML) from mCRPC cohort, including 27 pairs of matched PL-BML samples were finally included, Supplementary Fig. 1 showed the flowchart of sample selection process. Of note, although the mCRPC cohort predominantly consisted of castration-resistant cases (69/76, 90.8%), it also included a small proportion of metastatic HSPC (mHSPC) cases (7/76, 9.2%).

The baseline characteristics of patents included were summarized in Table 1 and Supplementary Table 1.

Table 1 Baseline characteristics of the patients with the matched PL-BML samples.

The PSMA mambranous H-score (MHscore) of PL samples did not show a significant correlation with age or tPSA levels (Fig. 1A, B), the PSMA MHscore was significantly higher in the ISUP = 2 HSPC group compared to the ISUP = 5 HSPC group (p = 0.026), while no significant differences were observed across the remaining subgroups (Fig. 1C and Supplementary Table 2). Figure 1D illustrates both the inter-patient and intra-tumoral heterogeneity of Membranous PSMA expression. The Shannon’s diversity index (SDI) reflects the heterogeneity abundance within each sample. SDI greater than 1, indicating relatively high intra-tumoral heterogeneity of PSMA expression, was observed in 86 (67.7%) PL samples. Twenty (15.7%) PL samples had SDI between 0.5 and 1, indicating moderate intra-tumoral heterogeneity. PSMA-negative was identified in 19 cases(15.0%) within the HSPC cohort. In the PSMA-negative HSPC cases, B7H3 and TROP2 demonstrated significantly better complementarity to PSMA and exhibited higher MHscore levels compared to other TAAs (Fig. 1E, F and Supplementary Table 3). HER2 was consistently negative in the PSMA-negative HSPC cases.

Fig. 1
figure 1

Membranous PSMA expression of PL is heterogeneous in clinical samples obtained from patients with HSPC at initial diagnosis of PCa. (A) Scatter plot of PSMA MHscore and age, showing no correlation was observed. (B) Scatter plot of PSMA MHscore and log(tPSA + 1). Due to the wide distribution of tPSA, a mathematical transformation was applied. However, no significant correlation was observed before or after the transformation. (C) Box plot of PSMA MHscore distribution across ISUP groups. Significant differences across subgroups are marked in the plot. (D) Membranous PSMA expression at the initial HSPC status. Expression of membranous PSMA quantified by H-score and presented in order of increasing MHscore. Degree of PSMA heterogeneity was measured by Shannon’s diversity index (SDI) and depicted as heat map ranging from low heterogeneity (light green) to high heterogeneity (dark green). PSMA-negative cases are indicated with a dashed line in the plot with the definition as MHscore ≤ 20. (E) Head-to-head comparison bar plots showing the expression profiles of other well-established TAAs in PSMA-negative HSPC cases. (F) Box plot showing the comparison of MHscore for well-established TAAs in PSMA-negative HSPC cases. Significant differences across subgroups are marked in the plot, *represents p < 0.05, **represents p < 0.01, ***represents p < 0.001.

The heterogeneity of membranous PSMA expression in mCRPC cohort was exhibited in Fig. 2A. SDI greater than 1 was observed in 23 (30.1%) BML samples, twenty-five (32.9%) BML had SDI between 0.5 and 1. PSMA-negative was identified in 28 cases (36.8%) within the mCRPC cohort. PSMA MHscore did not show any differences across the various bone-metastatic sites (Fig. 2B). Both STEAP1 and B7H3 exhibited significantly higher complementarily to PSMA, with higher MHscore levels in PSMA-negative mCRPC cases. Additionally, HER2 MHscore levels greater than 50 were observed in 28.6% of cases (Fig. 2C, D and Supplementary Table 4). TF emerged as the lowest expressing TAA in PSMA-negative mCRPC cases. The normalized membrane ratio (NMR) analysis revealed that STEAP1 and B7H3 displayed consistently high NMR expression patterns, while HER2 showed relatively heterogeneous NMR pattern, with 21.4% of PSMA-negative mCRPC cases showing a cytoplasmic-dominant HER2 expression pattern (Fig. 2E).

Fourteen mCRPC cases were identified in which neuroendocrine differentiation (NE) markers (Syn/CgA/NSE/CD56) had been assessed from the clinical pathology system of our center. Among these, three cases (21.4%) were classified as neuroendocrine prostate cancer (NEPC). Details of NE marker staining and PSMA expression level data of these 14 cases were presented in Supplementary Table 5. All the 3 NEPC cases were PSMA-negative cases. Using the Wilcoxon rank-sum test, a trend toward lower PSMA membrane expression levels in NEPC compared to non-NEPC BML was observed (p = 0.088). The NMR of PSMA expression was found to be significantly lower in NEPC than in non-NEPC BML (p = 0.035), indicating that PSMA may be predominantly localized in the cytoplasm to some extent rather than on the membrane in NEPC.

A total of 55 mCRPC cases had complete clinicopathological records of AR IHC staining. Of these 55 cases, 6 were classified as AR-negative, 9 as low expression, 11 as intermediate expression, and 29 as high expression. Membranous PSMA expression levels differed significantly across AR expression groups (overall p < 0.001). A positive correlation was observed between membranous PSMA expression and AR status, with the AR-negative group exhibiting significantly lower PSMA MHscores than other AR groups (p < 0.05), as shown in Supplementary Fig. 2 and Supplementary Table 6. These findings suggest that negative or low AR expression may be associated with reduced PSMA expression in BML. All six AR-negative mCRPC cases showed PSMA- negative expression, and the expression profiles of other actionable TAAs in these BML samples are detailed in Supplementary Table 7.

Fig. 2
figure 2

Membranous PSMA expression of BML is heterogeneous in clinical samples obtained from patients with mCRPC. (A) Membranous PSMA expression in mCRPC status. Expression of membranous PSMA quantified by H-score and presented in order of increasing MHscore. Degree of heterogeneity in mPSMA was measured by SDI and depicted as heat map ranging from low heterogeneity (light green) to high heterogeneity (dark green). PSMA-negative cases are indicated with a dashed line in the plot with the definition as MHscore ≤ 20. (B) Box plot of PSMA MHscore distribution across bone metastatic sites, with no significant difference was observed. (C) Head-to-head comparison bar plots showing the expression profiles of other well-established TAAs in PSMA-negative mCRPC cases. (D) Box plot showing the comparison of MHscore for well-established TAAs in PSMA-negative mCRPC cases. Significant differences across subgroups are marked in the plot, *represents p < 0.05, **represents p < 0.01, ***represents p < 0.001. (E) Heatmap of NMR classification across TAAs in PSMA-negative mCRPC cases. NMR > 50% indicates a membranous dominant expression pattern, while NMR < 50% indicates a cytoplasmic dominant expression pattern.

Overall, membranous PSMA expression in PL is significantly higher than in BML (p<0.001) (Fig. 3A). Meanwhile, the NMR levels in PL are also significantly higher than in BML(p<0.001) (Fig. 3B). In matched samples, although some patients exhibited increase in membranous PSMA expression of BML compared to the initial PL, 63.0% of patients showed decrease in membranous PSMA expression of BML. The treatment history of these patients during the interval between the two status is shown in Fig. 3C. There were no significant differences of the PSMA alteration (dPSMA) across the different categories of SGARA, chemotherapy, or radiation in matched samples (Supplementary Table 8). Figure 3D-G illustrates the different patterns of intra-patient PSMA expression heterogeneity.

Fig. 3
figure 3

Comparison of PSMA expression profile between HSPC-PL and mCRPC-BML status. (A) Box plot showing the comparison of PSMA MHscore between PL and BML, PSMA expression level in BML was significant lower than PL. (B) Violin plot showing the comparison of PSMA NMR between PL and BML, PSMA NMR in BML was also significant lower than PL. (C) Waterfall plot showing the alteration of membranous PSMA expression from PL to BML, along with the corresponding treatment history heatmap (second generation androgen receptor antagonist (SGARA), chemotherapy and radiotherapy), no significant correlation was observed between PSMA alteration and any of the various treatment history. (D, E,F, G) Different patterns of intra-patient PSMA expression heterogeneity in matched samples. (D) Heterogeneously low expression in different regions of the same HSPC sample (: within same regions, even adjacent tumor tissues can exhibit heterogeneous PSMA expression level, the tumor tissue in the red dashed area shows absent PSMA expression;completely absent PSMA expression tumor region), with high and uniform PSMA expression in matched mCRPC. (E) High PSMA expression in HSPC, low expression in matched mCRPC. (F) Heterogeneously low expression in both HSPC and matched mCRPC. (G) High expression in primary HSPC, heterogeneously low expression in both matched synchronously mHSPC and metachronously mCRPC. All the black scale bars in the IHC image represent 100 μm. Brown scale bars in the IHC image represent 500 μm.

Discussion

PSMA-targeted theranostic strategies, including PSMA-PET imaging and 177Lu-PSMA-617 radioligand therapy, have significantly revolutionized the clinical management of PCa. However, a subset of patients derived limited benefit from these approaches, in part due to heterogeneous or absent PSMA expression1,5,17,18emphasizing the urgent need for alternative theranostic targets. This study systematically evaluated the spatiotemporal inter- and intra-patient heterogeneity of PSMA expression in HSPC and mCRPC, revealing the PSMA-negative expression rate, and exploring potential alternative TAA expression in PSMA-negative cases. To our knowledge, this may be the first report to comprehensively head-to-head assess the expression of a panel of clinically actionable TAAs in PSMA-negative PCa within the same cohort, using real-world clinical samples of PCa, it may provide some insights for the clinical management of PCa in future.

Although the heterogeneity of PSMA expression in PCa has been documented, the detailed patterns of the heterogeneity remain elusive5,9,18,19. This study attempts to comprehensively characterize PSMA heterogeneity at three levels: inter-patient, intra-patient, and intra-tumoral. Consistent with previous literature, inter-patient heterogeneity of PSMA expression is evident in both HSPC-PL and mCRPC-BML in our study. Inter-patient variability in PSMA expression was observed in both PL and BML, with H-score ranging from 1.42 to 197.16. Notably, we observed significantly lower PSMA expression levels in BML compared to PL, which aligns with prior findings indicating that a considerable proportion of mCRPC patients exhibit PSMA expression loss19. Our study further revealed spatiotemporally intra-patient heterogeneity patterns of PSMA expression in 27 patients with matched PCa samples, with only 2 cases consistently showing negative PSMA expression, while 6 cases exhibiting HSPC-/mCRPC + and 7 cases exhibiting HSPC+/mCRPC- PSMA expression patterns in HSPC and mCRPC status, indicating that PSMA may display certain degrees of spatiotemporally dynamic expression. To further explore intra-tumoral diversity, PSMA expression was quantified through SDI analysis. We found that over 60% of PCa cases exhibited at least moderately intra-tumoral heterogeneity, different levels of PSMA expression even observed across adjacent tumoral sites within the same tissue sample in IHC images, reflecting the complexity of tumor subgroups within the tumor microenvironment.

Tumor components with low or absent PSMA expression may not be effectively targeted, resulting in residual disease. To address the clinical challenge posed by the loss of PSMA expression in these cases, we further investigated the expression of alternative clinically actionable TAAs as potential theranostic targets. B7H3 and STEAP1 were both first described more than two decades ago and have been identified as highly expressed cell surface antigens in PCa15,20. Our findings support the high expression levels of these two TAAs in PSMA-negative mCRPC, along with relatively homogeneous expression and predominant membrane localization (Fig. 2D), suggesting their potential as effectively complementary targets to PSMA. Future multi-target strategies combining targeted radionuclide therapy with ADCs, chimeric antigen receptor T-cell therapy (CART) or bi-specific antibody drugs may significantly improve comprehensive tumor targeting and killing in lethal mCRPC.

Although PSMA is primarily a membrane antigen, many cases in our study exhibited both membranous and cytoplasmic staining, with notably higher cytoplasmic levels in mCRPC compared to primary HSPC. This finding is consistent with recent studies in PCa, which also demonstrated membranous-cytoplasmic co-staining in IHC15,21. Paschalis et al. highlighted the importance of distinguishing these staining patterns, suggesting that membranous PSMA may represent more clinically actionable target, with cytoplasmic expression affecting membranous staining assessment5. Conversely, Pereira et al. believed that endocytosis and recycling of membranous receptors may explain this co-expression phenomenon6. These dynamic processes can influence TAA density on the cell surface and contribute to cytoplasmic staining. Garassino MC et al. found that in lung cancer, the NMR index of TROP2, incorporating both membranous and cytoplasmic expression, was better at predicting prognosis and therapeutic responses to TROP2-ADCs than membranous expression alone7. Therefore, we suggest that cytoplasmic PSMA staining should not be regarded merely as a limitation but integrated into the systematic assessment of PSMA expression. Another noteworthy observation from the literature is that Watanabe et al. reported PSMA-expressing tumor endothelial cells in approximately 12% of prostate cancer patients, with PSMA expression in endothelial cells being predominantly cytoplasmic. Their findings suggest a potential role for PSMA in promoting angiogenesis in vascular endothelial cells22. However, we would like to clarify that the present study focused exclusively on the expression of targetable markers in tumor cells within prostate cancer tissues. Tumor cells and endothelial cells are morphologically distinct and readily distinguishable. Moreover, the density of endothelial cells within the analyzed tumor regions was extremely low, therefore the endothelial cells were not evaluated in this study. Nonetheless, this phenomenon is intriguing and warrants further investigation in future studies.

Despite the comprehensive profiling of PSMA expression heterogeneity and the exploration of a panel of clinically actionable well-established TAAs in PSMA-negative PCa using large-population real-world clinical samples, this study has certain limitations. Firstly, the single-center, retrospective nature of this study may introduce selection bias. Secondly, although neuroendocrine differentiation has been increasingly reported in advanced prostate cancer, its prevalence in Asian populations appears to be relatively low, and relevant clinical evidence remains limited23,24. In this context, systematic evaluation of neuroendocrine markers was not performed in our study due to the limited availability of serial sections from the clinical samples. Finally, this study primarily focused on the systemic profiling of PSMA expression, it did not correlate target expression with clinical prognosis, nor did it delve into the underlying molecular mechanisms driving the observed phenomena. Future research should address these limitations to further explore the translational potential of our findings.

In conclusion, our study provides a comprehensive analysis of PSMA expression in PCa, highlighting significant heterogeneity in both HSPC and mCRPC, underscoring the importance of considering the dynamic nature of PSMA expression. We also identified alternative theranostic TAAs in PSMA-negative PCa, offering valuable insights into potentially alternative diagnostic and treatment options for future clinical applications in PCa.

Methods

Patients and tissue samples

The experiments involved in this study were approved by the Ethics Committee of Peking University First Hospital (PKUFH) with exemption from informed consent (2023-289-001). All methods were performed in accordance with the protocol approved by the Ethics Committee.

Formalin-fixed paraffin-embedded (FFPE) tissue samples were retrospectively collected from prostate biopsies and bone metastasis surgeries of patients with PCa at PKUFH between 2013 and 2023.

A total of 127 patients had sufficient HSPC-PL samples, including 27 with matched BML samples, the remaining 100 samples were randomly selected from the central database using stratified sampling, with equal weighting assigned to each ISUP group to minimize selection bias, as ISUP group is the most influential clinicopathological factor for PCa. Seventy-six metastatic PCa samples were totally included, with 27 matched BML samples (7 mHSPC and 20 mCRPC cases) and 49 non-matched BML samples (49 mCRPC cases). All samples had sufficient FFPE tissue available for analysis.

Immunohistochemistry

IHC was performed on both cohorts to evaluate PSMA expression. Antigen retrieval was carried out in EDTA buffer (pH 9.0) using microwave heating. Endogenous peroxidase was blocked with 3% H2O2, followed by serum blocking with 3% BSA for 30 min at room temperature. Sections were incubated overnight with primary antibody at 4 °C, then with secondary antibody for 50 min at room temperature. DAB was applied for color development, followed by counterstaining with hematoxylin, differentiation, dehydration, and cover slipping.

Digital images of the slides were captured using Pannoramic DESK, P-MIDI, and P250 scanners (3DHISTECH, Hungary), with the Pannoramic Scanner software, and PSMA expression was quantified using HALO 3.2 software to calculate both the membranous H-score (MHscore) and cytoplasmic H-score (CHscore): H-score = (% weak staining × 1) + (% moderate staining × 2) + (% strong staining × 3). Three random fields per sample were evaluated, and the average score was calculated25,26. The normalized membrane ratio (NMR) was defined as: MHscore/(MHscore + CHscore). PSMA-negative PCa was defined as PSMA MHscore ≤ 205,9. Consecutive sections were derived from the PSMA-negative samples to further assessing the expression of additional TAAs: HER2, NECTIN4, TROP2, TF, STEAP1 and B7H3. Details on the antibodies usage can be found in Supplementary Table 9.

Data collection and statistical analysis

Clinical and pathological data, including age, tPSA, ISUP group and bone metastasis sites were collected. For matched samples, treatment history (ADT, second generation androgen receptor antagonist(SGARA), chemotherapy, radiotherapy) during the sampling interval was also recorded. In the mCRPC cohort, IHC results for AR and neuroendocrine differentiation-related markers, including Syn, CgA, CD56, and NSE, as well as PSA and PSAP, were re-collected and reviewed from the clinical pathology registry system of PKUFH. NEPC was defined based on the following criteria: at least two NE markers positive, negative staining for PSA/PSAP, and morphological features consistent with small-cell carcinoma. The interpretation of AR immunohistochemical (IHC) staining was primarily based on the proportion of AR-positive tumor cells, supplemented by staining intensity: Tumors with < 25% AR-positive cells were classified as AR-negative, those with 25–50% as low expression, 50–75% as intermediate expression, and > 75% as high expression. Linear regression was used to assess the relationship between MHscore and age, tPSA. PSMA MHscore. Other TAA MHscore in PSMA-negative PCa were compared using the Wilcoxon signed-rank test. Intra-sample heterogeneity was quantified by Shannon’s diversity index (SDI). All figures were generated using R 4.3.2 and Adobe Illustrator 27.0. A significance level of 0.05 was used.