Visual Abstract
Abstract
177Lu-PSMA-617 is an effective therapeutic option in metastasized castration-resistant prostate cancer (mCRPC). However, some patients progress under treatment. We hypothesized that the tracer kinetics within the metastases may influence the therapy effectiveness and tested this hypothesis by analyzing uptake parameters on 2 consecutive posttherapy SPECT/CT scans. Methods: mCRPC patients treated with 177Lu-PSMA-617 and with available posttherapy SPECT/CT imaging (24 and 48 h after the first treatment) were enrolled retrospectively. Volumes of interest were defined on lymph node metastasis (LNM) and bone metastasis (BM) on both SPECT/CT scans. The reduction of the percentage injected dose (%IDred) between the 2 SPECT/CT scans was computed. We compared %IDred of responders (prostate-specific antigen drop ≥ 50% after 2 cycles of 177Lu-PSMA-617) and nonresponders. We tested the association of %IDred with progression-free survival and overall survival (OS) using a univariate Kaplan–Meier (KM) analysis and a multivariate Cox regression model. Results: Fifty-five patients (median age, 73 y; range, 54–87 y) were included. %IDred in LNM and BM was greater in nonresponders than in responders (for LNM, 36% in nonresponders [interquartile range (IQR), 26%–47%] vs. 24% in responders [IQR, 12%–33%] [P = 0.003]; for BM, 35% in nonresponders [IQR, 27%–52%] vs. 18% in responders [IQR, 15%–29%] [P = 0.002]). For progression-free survival, in KM analysis, greater %IDred in LNM (P = 0.008) and BM (P = 0.001) was associated with shorter survival, whereas in multivariate analysis, only %IDred in LNM was retained (P = 0.03). In univariate KM analysis of OS, greater %IDred in BM was associated with shorter survival (P = 0.002). In multivariate OS analysis, BM %IDred (P = 0.009) was retained. Conclusion: The 177Lu-PSMA-617 clearance rate from mCRPC metastases appears to be a relevant prognosticator of response and survival, with faster clearing possibly signaling a shorter radiopharmaceutical residence time and absorbed dose. Dual-time-point analysis appears to be a feasible and readily available approach to estimate the likelihood of response and patients’ survival.
Lutetium-177-PSMA-617 radioligand therapy (RLT) is a treatment option for prostate-specific membrane antigen (PSMA)–positive metastasized castration-resistant prostate cancer (mCRPC), which showed a survival rate superior to supportive treatment and a higher response rate with fewer adverse events than second-line chemotherapy (1,2). This treatment can induce a satisfactory response (prostate-specific antigen [PSA] decrease > 50%) in more than half of the treated subjects and hinder progression in 30% of them (2–6). It is safe and well tolerated, with toxicity-related discontinuation being rare (7).
However, some patients progress during the treatment, despite having shown adequate PSMA uptake on the pretherapy PET (8), which should have implied adequate radiopharmaceutical binding and tumor radiation exposure. Nonetheless, the tumor dose depends on the area under the curve of uptake intensity and 177Lu-PSMA-617 residence time within the lesions; those with a fast radiopharmaceutical turnover could therefore receive an insufficient radiation dose. Moreover, the 177Lu-PSMA-617 kinetics could vary across metastasis types (lymph node metastasis [LNM], bone metastasis [BM], and organ metastasis localizations). Estimating the radiopharmaceutical residence time from PET images, using multiple acquisitions, is hardly feasible, given the short half-life of the current diagnostic radionuclides 68Ga and 18F. In opposition, this parameter could be estimated on the post-RLT scans using the γ-radiation component of 177Lu-PSMA-617 (9).
A recent study confirmed that the dose to the metastases is pivotal in inducing response: patients receiving less than 10 Gy to the tumor localizations were unlikely to achieve a response (10). A full dosimetry approach is time-consuming and requires the patients to undergo a series of diagnostic acquisitions up to 7 d after 177Lu-PSMA-617 administration (10). Therefore, it is not yet used in the clinical routine. However, in theory, the PSMA kinetics within the lesions could be estimated with as few as 2 imaging time points, which may allow identification of the rate of change between the scans. In this study, we measured the percentage injected dose (%ID) within the mCRPC secondary lesions on 2 posttherapy SPECT/CT scans; the objective was to test whether the percentual change of this parameter correlates with the patients’ response, as well as with progression-free survival (PFS) and overall survival (OS).
MATERIALS AND METHODS
Patients and Ethics
The local database was retrospectively searched for mCRPC patients treated with 177Lu-PSMA-617 between April 2016 and September 2020. Criteria of eligibility for the treatment were histologically confirmed prostate cancer with progression to mCRPC status, adequate PSMA expression on a recent (<60 d) PET scan (18F-PSMA-1007 or 68Ga-PSMA-11), and approval by a registered multidisciplinary tumor board. Criteria of inclusion in the study were availability of dual-time-point SPECT/CT imaging (24 and 48 h after injection) at the time of the first RLT and at least 2 cycles of 177Lu-PSMA-617. There was no set maximum number of RLT cycles; treatment was carried on until complete response, disease progression, or intolerable therapy-related toxicity.
Exclusion criteria were unavailability or insufficient quality of SPECT/CT imaging, neuroendocrine histology, presence of PSMA-negative or 18F-FDG–positive disease (11), and later onset of PSMA-negative progression. Approval for this retrospective study was obtained from the local ethics committee of the University Hospital of Tübingen (decision 672/2019BO2), and the need to obtain specific consent was waived.
177Lu-PSMA-617 SPECT/CT Acquisition and Reconstruction
Images were acquired 24 and 48 h after 177Lu-PSMA-617 administration on a hybrid, dual-head γ-camera (NM/CT 670pro; GE Healthcare) equipped with medium-energy collimators. Whole-body planar acquisition (bed speed, 20 cm/min) and 2 or 3 SPECT/CT scans, to cover the field of view between the vertex and the mid thigh, were acquired. SPECT data (obtained with a 180° range, 6° steps, and a step time of 20 s) were reconstructed using a 3-dimensional, ordered-subset expectation maximization algorithm (with CT-based attenuation correction, scatter correction, and resolution recovery). CT had low-dose settings (30 mA, 120 kV). The SPECT and CT images were then coregistered (Xeleris V.4; GE Healthcare).
Image Analysis
SPECT/CT images acquired after the first 177Lu-PSMA-617 administration were analyzed using a commercial application (Volumetrix MI; GE Healthcare). On loading the first SPECT/CT scan, the injected and residual activity, as well as calibration and administration times, were inserted. Volumes of interest (VOIs) were constructed semiautomatically on the LNM, BM, and organ metastasis localization by first encasing the target lesions in a large VOI and then applying automatic thresholding (41% of the maximum voxel value); finally, nontarget areas were manually excluded. The minimum VOI was set at 5 cm3. The sum of all VOIs was labeled tumor burden.
For each VOI, the system calculated the mean and maximum activity concentrations and %ID (that is, the percentage of the initial injected dose that was retained within a specific volume in the 2 time points). We computed this parameter for the LNM, BM, and organ metastasis localizations separately, as well as a single score including all disease localizations (%ID total). For the main purposes of our analysis, we considered this last parameter only because of its robustness (12). Concentration indices were tested only in the differentiation between responders and nonresponders; we did not compute weight-normalized indices, such as SUVmean or SUVmax, because they can underperform in the evaluation of tissue-specific tracers (13). All VOIs were then exported onto the second dataset, and positioning adjustments were made whenever necessary. Finally, we calculated the relative reduction of %ID (%IDred) from the first to the second SPECT/CT scan in each district (LNM, BM, and organ metastasis).
Study Objectives
The first objective was to test whether there are differences in %IDred between responders (PSA drop ≥ 50% after 2 RLTs) and nonresponders. In a subanalysis, patients were stratified into good responders (PSA drop ≥ 50%), partial responders (PSA drop, 1%–50%), and nonresponders (no PSA drop or progression). The PSA variation was calculated as the percent difference between the baseline value (on the day of the first RLT, before 177Lu-PSMA-617 administration) and the measurement 1 mo after the second RLT cycle.
The second objective was to identify whether %IDred or baseline %ID in any of the considered disease localizations correlated with PFS or OS. All patients had routinely undergone a clinical or laboratory data evaluation 1 mo after each treatment, and imaging with 68Ga-PSMA-11 or 18F-PSMA-1007 PET/CT had been performed 1 mo after every second RLT. Parameters related to the patients’ clinical history, laboratory values, and imaging-related indices were tested; these indices were chosen on the basis of clinical reasoning and available literature data (14).
PFS was defined as the time elapsed between RLT start and progression; that is, 2 consecutive increases in PSA levels totaling at least 25% of the nadir (15), appearance of new metastatic lesions on any imaging, worsening of disease-related symptoms requiring an adjustment in supportive therapy, deterioration of general conditions, or any instance in which an mCRPC therapy change was required. OS was defined as the time between RLT start and confirmed disease-related death of the patient.
Statistical Analysis
Data are presented as median and interquartile range (IQR) unless otherwise specified. Comparisons between continuous variables between groups were made using the unpaired Student t test (2 groups) or 1-way ANOVA with least significance difference post hoc comparison (3 groups) for normally distributed data and using the Mann–Whitney test (2 groups) or Kruskal–Wallis test (3 groups) for nonnormally distributed data. Categoric variables were compared using the χ2 or Fisher exact test, as appropriate. Correlations between continuous variables were made by means of bivariate analysis, using the Pearson coefficient. Predictors of PFS and OS were tested using the Kaplan–Meier estimator (dichotomized or categoric variables) and with a univariable Cox regression model. A multivariable model was then constructed using a combination of statistically significant image-derived and clinical parameters. The multivariable analysis was performed with the Cox proportional hazard model, and the proportional hazard assumptions were tested using Schoenfeld residuals. Multicollinearity was checked by matrix correlation. A P value of less than 0.05 was considered statistically significant.
RESULTS
Patients’ Cohort Characteristics
Initially, 125 patients treated with 177Lu-PSMA-617 RLT were identified. Sixty subjects were excluded because of unavailable 2–time point imaging. Five more patients were excluded who, because of mixed 68Ga/18F-PSMA and 18F-FDG–positive disease or neuroendocrine histology, had received 177Lu-PSMA-617 as a disease control attempt due to the lack of viable alternatives. Five patients were excluded because of treatment discontinuation unrelated to disease progression. Thus, 55 patients were ultimately included (Fig. 1).
Table 1 shows an overview of the 55 included patients. The subjects had a long disease history, with a median of 8.3 y between diagnosis and RLT (IQR, 3.7–14.2 y). The median time between diagnosis and mCRPC progression was 3.1 y (IQR, 1.4–9.1 y); afterward, the patients were treated with standard approaches for a median time of 3.1 y (IQR, 2.1–4.9 y) before 177Lu-PSMA-617.
On RLT start, 43 (78%), 43 (78%), or 14 (26%) subjects had LNM, BM, or organ metastasis, respectively. The latter category included liver (n = 6, 43%), lung (n = 7, 50%), and meningeal (n = 1, 7%) localizations. Subjects received a median of 3 cycles of 177Lu-PSMA-617 RLT (IQR, 2–8). After the first 2 treatment cycles, 25 (45%), 12 (22%), and 18 (33%) patients were classified as responders, partial responders, and nonresponders, respectively. Tumor volume and %ID were not significantly different between responders and nonresponders.
Tracer Kinetics Across Responders and Nonresponders
Median values of %IDred within LNM were significantly higher in nonresponders (36%; IQR, 26%–47%) than in responders (24%; IQR, 12%–33%; P = 0.003; Fig. 2). When examining the variation across good, partial, and nonresponders, an increasing trend of %IDred was noted, indicating rapid tracer clearance in partial and nonresponders (partial responders, 30% [IQR, 26%–38%]; nonresponders, 42% [IQR, 33%–57%]; P = 0.001 between nonresponders and good responders; Fig. 2).
A similar pattern was detected in BM. Nonresponders had a more pronounced median %IDred drop than did responders (nonresponders, 35% [IQR, 27%–52%]; responders, 18% [IQR, 15%–29%]; P = 0.002). Partial responders had a large median %IDred overlap with nonresponders (partial responders, 28% [IQR, 23%–55%]; nonresponders, 37% [IQR, 29%–46%]; P = 0.002; Fig. 2). Figure 3 shows visual representations of the activity kinetics in BM.
In the patients’ subgroup with organ metastasis, a visual and nonsignificant trend toward a difference in %IDred between responders and nonresponders could be observed (Fig. 2).
Compared with %ID, mean radioactivity concentration appeared to have at least noninferior performance in telling apart responders from nonresponders; maximum radioactivity concentration seemed to perform less well (Supplemental Table 1 [supplemental materials are available at http://jnm.snmjournals.org]).
Finally, there was a direct correlation between %IDred in LNM and %IDred in BM (R = 0.494, P = 0.005), as well as between BM and organ metastasis (R = 0.763, P = 0.006). Conversely, no correlation was found between LNM and organ metastasis (R = 0.037, P = 0.908).
Prognostic Factors of PFS
Progression occurred in all patients, with a median time to the event of 7.5 mo (95% CI, 4.7–10.3 mo). An %IDred value above the median in both LNM and BM was associated with shorter survival (Fig. 3). Moreover, higher %ID within the LNM correlated with longer survival (Fig. 4). The kinetics of 177Lu-PSMA-617 within the organ metastasis had no prognostic implications (data not shown).
Supplemental Table 2 and Table 2 present an overview of all tested PFS predictors in the Kaplan–Meier and simple proportional Cox regression analyses, respectively. In the multivariable Cox regression analysis, which consisted of all factors showing a significant association with PFS in the univariable tests, only %IDred in the LNM retained a significant association with PFS (Table 3).
Prognostic Factors of OS
Forty-three patients died of disease. OS from the RLT start was 16.3 mo (95% CI, 11.1–21.6). %IDred in BM predicted OS, but neither LNM nor BM %ID was associated with OS (Fig. 5).
All tested factors in Kaplan–Meier and univariable Cox regression are presented in Supplemental Table 3 and Table 4, respectively.
In the multivariable Cox regression analysis, indices of tumor burden, as well as the degree of %IDred in BM, were associated with a worse outcome; in contrast, baseline activity within the BM lesion reduced the risk (Table 5). Similar results were obtained using the %ID total instead of the BM-specific one (Supplemental Table 4).
DISCUSSION
This study presents evidence of a prognostic role of a kinetics analysis of posttherapy images in patients treated with 177Lu-PSMA-617 RLT. In particular, our data show that assessment of 2 time points of radiopharmaceutical distribution can reveal diverging patterns of uptake and washout from prostate cancer cells. Patients who will respond to the treatment show 177Lu-PSMA-617 activity remaining stable or undergoing a moderate drop from the first to the second observation point. Conversely, nonresponders show more significant activity reduction within the metastases. Finally, the activity drop between the 2 time points appears to be a strong and independent prognostic marker of both PFS and OS.
The kinetics analysis suggests that a more significant drop in activity could be linked with faster tracer washout from the cells and therefore with a shorter radiopharmaceutical residence time, which could entail a smaller delivered dose to the target lesions. Conversely, stable radioactivity concentration indices could signal a longer residence time of 177Lu-PSMA, resulting in a higher delivered dose and better RLT effectiveness. This hypothesis is supported by the correlation between the uptake variation within the lesion and the long-term outcome and, more specifically, by the difference in the determinants of PFS and OS: the kinetics within the LNM was predictive of disease progression, whereas the uptake patterns within the BM metastases correlated with OS. mCRPC LNM disease can quickly escalate if left unchecked or undertreated (16,17). However, LNM can only rarely cause morbidity or mortality per se. Conversely, the overall burden of BM has a proven impact on skeletal function and survival (18–20), and failure in controlling the metastatic infiltration of the skeleton can effectively reduce patients’ OS. Our data confirmed these notions, showing that tracer uptake, washout rate, and tumor load were relevant to OS.
Existing data on 177Lu-PSMA-617 RLT dosimetry indicate that the radiation dose delivered to the target is key in ensuring response (21,22). However, in a recent systematic review, a relevant variability in tumor dose (range, 7.5–77.6 Gy) was highlighted across and within populations (23). Technical issues and differences in tracer uptake intensity could partly account for this variability. However, because the administered activity per cycle is fixed (7,400 MBq), and given that all lesions show adequate PSMA expression (24), it could be hypothesized that tracer uptake and washout kinetics could be the main determinants of dose variability.
In this study, we tested the efficacy of a dual-point tracer concentration measurement approach as a sort of dosimetry proxy. Such an approach could be implemented as a quick clinical tool to identify the kinetics pattern shortly after the first treatment administration, potentially enabling salvage in cases of accelerated 177Lu-PSMA-617 washout. However, this information can be obtained only after the first treatment cycle has been administered. To overcome this limitation, further studies could be designed to transfer this concept to the pre-RLT imaging setting. Using longer–half-life PSMA tracers, such as 64Cu-PSMA, would allow replication of the 177Lu-PSMA-617 setting (25,26). Another potential approach could be represented by scanners with an ultralong axial field of view, which allows very late imaging even when using radiopharmaceuticals with a relatively short half-life, such as 18F-PSMA-1007 (27,28).
Finally, modeling of the tracer kinetics could be used to implement tailored RLT. Specifically, the supposedly reduced residence time of the radiopharmaceutical in patients with rapid 177Lu-PSMA-617 washout could be compensated by increasing the dose delivered per unit of time. This task could be accomplished by increasing the administered activity or using hard-hitting radioisotopes, such as α-emitters, in lieu of 177Lu (29,30).
This study presents some limitations. It represents a single-center, retrospective investigation that is thereby susceptible to selection bias; strict selection criteria were used to avoid distortion. Even though we analyzed organ metastasis, we could not determine whether a difference in the tracer kinetics between responders and nonresponders exists. This limitation could result from the small sample size, as well as the heterogeneity of said metastasis, which were in different organs. Nonetheless, similarities in the washout behavior between organ metastasis and BM could be identified, as highlighted by the correlation analysis, hinting that these hallmarks of advanced disease tend to show comparable 177Lu-PSMA-617 kinetics. Treatment response was always assessed after the second RLT cycle without considering the maximum potential PSA decline; this choice was made to ensure consistency, because the latter information was not always available.
Data from the pre-RLT PET/CT were available, yet the tracer uptake intensity (SUV) was not factored into the univariate and multivariate analyses, given the large variability of time elapsed between imaging and RLT and that the tracer had been changed during the observation period (68Ga-PSMA-11 to 18F-PSMA-1007). Considering later time points in SPECT/CT, for example, 72–120 h after injection, could have improved the prognostic power of the washout index. However, the overall goal of this research line is to identify the kinetics trends as early and with as few time points as possible to make the technique eventually feasible for everyday routine. Finally, many factors can affect OS, especially in patients with a long disease history. To ensure homogeneity and assess the impact of clinical variables on the outcome, we adopted tight selection criteria and tested each of them with a survival analysis model.
CONCLUSION
The 177Lu-PSMA-617 kinetics in LNM and BM appear to be prognostic of treatment response, as well as of survival. In particular, OS appears to be linked with the kinetics parameters of BM, in keeping with the concept that BM is the most threatening pathogenic mechanism of mCRPC. This information, which should be confirmed by prospective trials, is readily obtainable from posttherapy scans and could be used to prognosticate treatment outcomes and design studies aimed to investigate the potential of PET-based prediction, as well as the possibility of patient-adapted therapeutic protocols of RLT.
DISCLOSURE
No potential conflict of interest relevant to this article was reported.
KEY POINTS
QUESTION: Are the 177Lu-PSMA-617 kinetics in mCRPC localizations relevant to treatment effectiveness?
PERTINENT FINDINGS: Greater reduction of 177Lu-PSMA-617 concentration within metastasis over time was associated with reduced response rate and shorter survival.
IMPLICATIONS FOR PATIENT CARE: The kinetics information could be used in studies on RLT tailoring to test the effectiveness of increase of activity or the use of higher-energy isotopes in bolstering the target dose in patients with faster washout.
ACKNOWLEDGMENT
This work was supported by the Deutsche Forschungsgemeinschaft (German Research Foundation, Germany’s Excellence Strategy, EXC2180-390900677).
Footnotes
Published online Jul. 6, 2023.
- © 2023 by the Society of Nuclear Medicine and Molecular Imaging.
REFERENCES
- Received for publication August 7, 2022.
- Revision received April 20, 2023.