ORIGINAL RESEARCH article

Front. Oncol., 18 May 2026

Sec. Genitourinary Oncology

Volume 16 - 2026 | https://doi.org/10.3389/fonc.2026.1835850

A multivariable prediction model combining 18F-PSMA PET/CT and mpMRI for clinically significant prostate cancer: development and validation

  • 1. Department of Urology, Peking University First Hospital, Beijing, China

  • 2. Institute of Urology, Peking University, Beijing, China

  • 3. National Urological Cancer Center, Beijing, China

Abstract

Background:

Multiparametric MRI (mpMRI) and prostate-specific membrane antigen (PSMA) PET/CT provide complementary information for prostate cancer diagnosis, but optimal integration strategies remain unclear. We aimed to develop and validate a multivariable model combining clinical and imaging parameters to predict clinically significant prostate cancer (csPCa).

Methods:

This study included 1305 consecutive patients with suspected prostate cancer who underwent both mpMRI and 18F-PSMA PET/CT with biopsy at Peking University First Hospital. The development and temporal validation cohorts were temporally divided; the development cohort was split into training (70%) and internal test (30%) sets. csPCa was defined as Gleason Grade Group ≥2. Predictors were selected using LASSO regression (1-standard error criterion) followed by forward stepwise logistic regression (AIC, p < 0.05). Model performance was assessed by AUC, calibration plots, and decision curve analysis.

Results:

The final model retained 3 predictors: PRIMARY score, PI-RADS score, and PSAD. The training AUC was 0.916. At the Youden Index cutoff (≥84%), sensitivity was 79.3% and specificity 76.6%; at the recommended screening cutoff (≥46%), sensitivity reached 96.0%. The AUC was 0.914 (95% CI: 0.882–0.941) in the internal test set and 0.837 (95% CI: 0.778–0.891) in temporal validation. The model showed good calibration (Brier score: 0.096) and superior clinical utility over individual imaging parameters.

Conclusions:

The multivariable model integrating 18F-PSMA PET/CT and mpMRI parameters provides accurate risk stratification for csPCa and may help optimize biopsy decisions.

Introduction

Prostate cancer (PCa) is the most commonly diagnosed malignancy in men worldwide (). While early detection through prostate-specific antigen (PSA) screening has increased the diagnosis of localized disease, it has also led to substantial overdiagnosis and overtreatment of indolent cancers that would never cause clinical harm (). This has driven a paradigm shift toward identifying clinically significant prostate cancer (csPCa), which requires active treatment, while avoiding unnecessary interventions for low-risk disease. However, the definition of csPCa has varied across studies, creating challenges for clinical decision-making and research comparisons. Based on the 2014 International Society of Urological Pathology (ISUP) consensus conference and subsequent clinical guidelines, csPCa is now widely defined as Gleason Grade Group (GG) ≥2 (Gleason score ≥3 + 4), a threshold that identifies tumors with metastatic potential warranting definitive treatment ().

Traditional screening with prostate-specific antigen (PSA) testing has been limited by its inability to distinguish csPCa from indolent disease, leading to overdiagnosis and overtreatment (). Multiparametric magnetic resonance imaging (mpMRI) has emerged as a valuable tool for pre-biopsy risk stratification. The Prostate Imaging Reporting and Data System (PI-RADS) provides a standardized framework for lesion assessment, with PI-RADS ≥3 triggering targeted biopsy consideration (). Landmark trials including PROMIS and PRECISION have established mpMRI as a gatekeeper for prostate biopsy (, ). However, mpMRI has notable limitations, including moderate interobserver variability and suboptimal specificity, particularly for PI-RADS 3 lesions (). PSA density (PSAD), calculated as serum PSA divided by prostate volume, has also been shown to improve diagnostic specificity ().

Prostate-specific membrane antigen (PSMA) positron emission tomography/computed tomography (PET/CT) has gained attention as a molecular imaging modality with high tumor specificity. PSMA is a type II transmembrane glycoprotein overexpressed in prostate cancer cells, making it an ideal target for molecular imaging (). While 68Ga-labeled PSMA tracers were initially developed, 18F-labeled PSMA ligands offer several advantages including longer half-life, higher image resolution, and cyclotron-based production enabling wider availability (). Our institution has developed 18F-PSMA-137, an 18F-labeled PSMA ligand with favorable pharmacokinetic properties and reduced urinary excretion, improving pelvic visualization (). The PRIMARY scoring system provides a standardized qualitative assessment framework for intraprostatic PSMA uptake ().

The integration of mpMRI and PSMA PET may provide complementary diagnostic information, combining anatomical detail with molecular characteristics. However, the optimal method for combining these modalities remains unclear. Simple rule-based approaches may sacrifice specificity, while requiring concordant positive findings may miss cancers.

Therefore, this study aimed to develop and validate a multivariable prediction model integrating clinical parameters, mpMRI findings, and 18F-PSMA PET/CT data to improve pre-biopsy risk stratification for csPCa.

Materials and methods

Study design and population

This diagnostic accuracy study employed a combined retrospective-prospective design with consecutive patient enrollment. Male patients aged ≥18 years with suspected PCa who underwent both mpMRI and 18F-PSMA PET/CT with prostate biopsy at Peking University First Hospital were enrolled retrospectively from August 2021 to August 2025 and prospectively from September 2025 to December 2025. All patients meeting the inclusion criteria during the respective time periods were consecutively included to minimize selection bias.

The indications for biopsy included: (1) suspicious nodules on digital rectal examination (DRE); (2) elevated serum PSA level (>4 ng/mL); or (3) suspicious lesions identified on imaging.

Patients were excluded based on the following criteria: (1) time interval exceeding 3 months between imaging modalities and biopsy; (2) history of prostate surgery, including transurethral resection of the prostate (TURP); (3) prior treatment for PCa, including androgen deprivation therapy (ADT), radiotherapy, or focal therapy; (4) incomplete imaging or pathological data.

The study was approved by the Institutional Ethics Committee (Approval No. 2025R0380-0002), with waiver of informed consent for the retrospective cohort and written informed consent obtained from all participants in the prospective cohort.

Clinical and pathological data

Baseline clinical characteristics including age at the time of biopsy, serum PSA level (ng/mL), and prostate volume (PV, in mL) measured on imaging were collected. PSA density (PSAD) was calculated as PSA/PV (). The reference standard was histopathological diagnosis from systematic (9 or 12-core) and/or MRI-targeted prostate biopsy. The reference standard for all patients was histopathological diagnosis via a prostate biopsy. Biopsy results were graded per ISUP 2014 consensus (). The primary endpoint of the study was csPCa defined as ISUP GG ≥ 2 (Gleason score ≥ 3 + 4 = 7). The patients who had GG 1 (Gleason score of 3 + 3 = 6) or benign pathology were deemed not to have csPCa.

Image acquisition and interpretation

All mpMRI examinations were performed on 1.5T or 3.0T scanners, with T2-weighted imaging (T2WI) and diffusion-weighted imaging (DWI) sequences acquired according to PI-RADS v2.1 guidelines (). 18F-PSMA PET/CT was performed 60 minutes after intravenous injection of 18F-PSMA-137. All images were reviewed by two independent readers (board-certified radiologists with over 10 years of experience in prostate MRI, and nuclear medicine physicians with over 10 years of experience in PSMA PET/CT) blinded to pathology results.

For mpMRI, the PI-RADS v2.1 scoring system was applied (). T stage was determined based on MRI findings according to the TNM classification criteria, assessing for prostate-confined lesion (T2), extracapsular extension (T3a), seminal vesicle invasion (T3b), and invasion of adjacent structures (T4). It should be noted that this MRI-based T staging differs from clinical T staging defined in international guidelines, which is based on DRE findings; MRI-based staging was used in this study as it provides more detailed anatomical assessment for pre-biopsy risk stratification.

For 18F-PSMA PET/CT, the PRIMARY scoring system was used for qualitative assessment (detailed scoring criteria provided in Supplementary Table 1) (). Briefly, a PRIMARY score using a combination of pattern information and SUVmax was assigned to each patient: Score 1 indicated no pattern and low-grade activity; Score 2 indicated diffuse transition zone (TZ) or symmetric central zone (CZ) activity without focal uptake; Score 3 indicated focal TZ activity visually greater than twice the background TZ activity; Score 4 indicated focal peripheral zone (PZ) activity; Score 5 indicated any pattern with SUVmax ≥12. The maximum standardized uptake value (SUVmax) of the dominant lesion was also recorded for quantitative analysis. Discrepancies were resolved by consensus.

Data splitting strategy

The study utilized a internal-temporal validation framework. The development cohort and temporal validation cohort were temporally divided with February 2025 as the cutoff point. The development cohort was randomly split using stratified sampling (maintaining outcome proportions) into: (1) a training set (70%) for model construction and coefficient estimation, and (2) an internal test set (30%) for initial performance evaluation on unseen data from the same time period. The temporally separated temporal validation cohort served as independent validation to assess model generalizability across different time periods.

This design allows assessment of: (1) reproducibility on contemporaneous unseen data in the internal test set, and (2) temporal transportability in the temporal validation cohort.

Statistical analysis

Continuous variables were expressed as medians with interquartile ranges (IQR) and compared using Mann-Whitney U test. Categorical variables were reported as frequencies with percentages and compared using chi-square test.

Continuous variables, including age (cutoff 60 years), PSAD (cutoff 0.15 ng/mL2), and SUVmax (cutoffs 4 and 8), were categorized based on clinically meaningful thresholds (). To address potential computational collinearity, PSAD was included as a candidate variable while PSA and prostate volume were excluded (). Candidate variables from the training set were entered into least absolute shrinkage and selection operator (LASSO) regression for initial variable selection (), with the optimal regularization parameter (λ) determined by 10-fold cross-validation using the 1-standard error (1se) criterion. The selected variables were then forwarded to a stepwise logistic regression model based on the Akaike Information Criterion (AIC) and a multivariable significance threshold of p < 0.05 for final model construction. The optimal cutoff of the model was determined by calculating the Youden Index, while a recommended screening threshold was also identified for optimal clinical prioritization. At these two specific thresholds, the sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were calculated.

Model assessment included internal testing and temporal validation. To address potential overfitting and assess model stability, internal validation was also performed using 1000 bootstrap resamples to calculate the optimism-corrected area under the curve (AUC). The probabilities calculated by the models were used to draw receiver operating characteristic (ROC) curves (). Model discrimination was evaluated using the area under the curve (AUC) with 95% confidence intervals (CI). Model calibration was assessed using calibration plots and the Brier score (). Decision curve analysis (DCA) evaluated clinical utility across a range of threshold probabilities ().

The study flowchart was shown in Figure 1. Statistical analyses were performed using Python 3.9 with scikit-learn and statsmodels packages (). Two-sided p < 0.05 was considered statistically significant.

Figure 1

Results

Patient characteristics

A total of 1305 consecutive patients were included: 894 in the development cohort (training set: n=625; internal test set: n=269) and 411 in the temporal validation cohort. The baseline characteristics are summarized in Table 1.

Table 1

CharacteristicDevelopment (n=894)Validation (n=411)p-value
Age, year68.0 [62.0-73.0]69.0 [63.0-73.0]0.129
PSA, ng/mL10.1 [6.9-14.0]10.6 [7.5-16.0]<0.001
PV, cmL40.5 [31.0-57.5]36.6 [27.4-49.1]<0.001
PSAD, ng/mL20.2 [0.1-0.3]0.3 [0.2-0.5]<0.001
SUVmax9.6 [5.5-17.6]10.3 [5.7-17.7]0.198
T stage<0.001
 ≤ T2605 (67.7%)310 (75.4%)
 T3a152 (17.0%)36 (8.8%)
 T3b119 (13.3%)51 (12.4%)
 T418 (2.0%)14 (3.4%)
PI-RADS Score0.182
 1-279 (8.8%)30 (7.3%)
 3153 (17.1%)56 (13.6%)
 4352 (39.4%)162 (39.4%)
 5310 (34.7%)163 (39.7%)
PRIMARY Score0.082
 1-2166 (18.6%)54 (13.1%)
 383 (9.3%)46 (11.2%)
 4275 (30.8%)136 (33.1%)
 5370 (41.4%)175 (42.6%)
Biopsy
 Total cores11.0 ± 3.010.4 ± 3.0<0.001
 Positive cores4.4 ± 2.95.2 ± 2.8<0.001
 Positive rate, %42.69 ± 28.7253.07 ± 29.0<0.001
ISUP Grade Group0.003
 Benign112 (12.5%)28 (6.8%)
 182 (9.2%)36 (8.8%)
 2285 (31.9%)152 (37.0%)
 3201 (22.5%)103 (25.1%)
 488 (9.8%)52 (12.7%)
 5126 (14.1%)40 (9.7%)
csPCa (ISUP≥2)700 (78.3%)347 (84.4%)0.012

Baseline patient characteristics (development and validation cohorts).

PSA, prostate specific antigen; PV, prostate volume; PSAD, PSA density; SUVmax, maximum standardized uptake value; PI-RADS, Prostate Imaging Reporting and Data System; ISUP, International Society of Urological Pathology; csPCa, clinically significant prostate cancer.

Continuous variables are presented as median [IQR]. Categorical variables are presented as count (%). P values from Mann-Whitney U test (continuous) or chi-square test (categorical).

csPCa was detected in 700 patients (78.3%) in the development cohort and 347 patients (84.4%) in the temporal validation cohort. The temporal validation cohort showed higher median PSA, PSAD, and SUVmax compared to the development cohort (all p < 0.05), reflecting temporal changes in referral patterns.

Predictor selection

Candidate variables from the training set (Table 2) were entered into LASSO regression. LASSO regression with 10-fold cross-validation using the 1se criterion identified 5 macro-variable groups: PSAD, MRI-based T stage, SUVmax, PI-RADS score, and PRIMARY score (Figure 2). These candidate variables were then subjected to forward stepwise regression based on AIC and multivariable p < 0.05. The stepwise process refined the predictors, ultimately retaining three core variables: PSAD, PI-RADS score, and PRIMARY score. (Table 3, Figure 3).

Table 2

VariablecsPCa (+) (n=489)csPCa (−) (n=136)p-value
Age, year68.0 [63.0-73.0]66.0 [61.0-72.0]0.010
< 6062 (12.7%)29 (21.3%)0.017
 ≥ 60427 (87.3%)107 (78.7%)
PSAD, ng/mL²0.2 [0.2-0.4]0.1 [0.1-0.2]<0.001
 < 0.1587 (17.8%)61 (44.9%)<0.001
 ≥ 0.15402 (82.2%)75 (55.1%)
T stage<0.001
 ≤ T2294 (60.1%)129 (94.9%)
 > T2195 (39.9%)7 (5.1%)
SUVmax11.6 [6.6-19.6]4.8 [4.0-6.2]<0.001
 < 417 (3.5%)33 (24.3%)<0.001
 4-8142 (29.0%)81 (59.6%)
 ≥ 8330 (67.5%)22 (16.2%)
PI-RADS Score<0.001
 1-28 (1.6%)47 (34.6%)
 363 (12.9%)47 (34.6%)
 4215 (44.0%)35 (25.7%)
 5203 (41.5%)7 (5.1%)
PRIMARY Score<0.001
 1-224 (4.9%)90 (66.2%)
 349 (10.0%)7 (5.1%)
 4176 (36.0%)27 (19.9%)
 5240 (49.1%)12 (8.8%)

Candidate variables characteristics of training set.

PSA, prostate specific antigen; PV, prostate volume; PSAD, PSA density; SUVmax, maximum standardized uptake value; PI-RADS, Prostate Imaging Reporting and Data System; ISUP, International Society of Urological Pathology; csPCa, clinically significant prostate cancer.

Continuous variables are presented as median [IQR]. Categorical variables are presented as count (%). P values from Mann-Whitney U test (continuous) or chi-square test (categorical).

Figure 2

Table 3

VariableOR95% CIp-value
PRIMARY Score3.052.36-3.94<0.001
PI-RADS Score3.682.66-5.11<0.001
PSAD0.007
< 0.15Reference
≥ 0.152.231.25-3.980.007

Variables of multivariable logistic regression model.

OR, odds ratio; CI, confidence interval; PSAD, prostate specific antigen density; PI-RADS, Prostate Imaging Reporting and Data System.

Forward stepwise regression analysis was used to select variables.

Figure 3

Among the retained imaging parameters, PI-RADS score showed the strongest association with csPCa (OR = 3.68, 95% CI: 2.66-5.11, p < 0.001). PRIMARY score also demonstrated strong independent predictive value (OR = 3.05, 95% CI: 2.36-3.94, p < 0.001). Elevated PSAD (≥0.15 ng/mL2) was a significantly independent clinical predictor (OR = 2.23, 95% CI: 1.25-3.98, p = 0.007).

Model development and performance

The final multivariable model strictly included three predictors: PSAD, PI-RADS score, and PRIMARY score. A point-based nomogram incorporating the three selected predictors was constructed to facilitate clinical application (Figure 4). In this nomogram, each predictor is assigned a specific point value on the top axis, and the sum of these points maps directly to the estimated probability of csPCa on the bottom axis (). The optimal threshold determined by the Youden index was ≥84% probability, which yielded a sensitivity of 79.3% and specificity of 76.6%. A clinically recommended threshold of ≥46% probability was also provided to maximize csPCa detection, achieving a sensitivity of 96.0% and specificity of 46.6% (Figure 5, Table 4).

Figure 4

Figure 5

Table 4

Cut-off values of predictive probabilitiesSensitivity (ratio) 95% CISpecificty (ratio) 95% CIPPV (ratio) 95% CINPV (ratio) 95% CI
Youden
index≥84%
79.3% (275/347)
74.7%-83.2%
76.6% (49/64)
64.9%-85.3%
94.8% (275/290)
91.6%-96.8%
40.5% (49/121)
32.2%-49.4%
recommended
value≥46%
96.0% (333/347)
93.3%-97.6%
46.9% (30/64)
35.2%-58.9%
90.7% (333/367)
87.3%-93.3%
68.2% (30/44)
53.4%-80.0%

Clinical Application cut-offs values and model diagnostic performance.

CI, confidence interval; PPV, positive predictive value; NPV, negative predictive value.

The Youden index cut-off values of predicting probabilities were obtained at the maximum Youden index (sensitivity + specificity − 1).

We recommended a cut-off value of 46% for the csPCa-predicting nomogram to achieve the sensitivity of 96.0%.

The nomogram allows calculation of individualized probability of csPCa based on patient-specific clinical and imaging characteristics.

Model test and validation

To rigorously evaluate potential overfitting, internal validation using 1000 bootstrap resamples was performed. The apparent AUC was 0.9156 with a low calculated optimism of 0.0026, yielding an optimism-corrected AUC of 0.9130. In the internal test set, the model achieved an AUC of 0.914 (95% CI: 0.882-0.941) (Figure 6A). In the temporally separated temporal validation cohort, the model maintained an AUC of 0.837 (95% CI: 0.778-0.891) (Figure 6B). Furthermore, a head-to-head comparison in the temporal validation cohort demonstrated that our proposed model (AUC 0.837) possessed superior discriminative ability compared to the established PCPT 2.0 (AUC 0.623) and ERSPC RC3 (AUC 0.741) risk calculators (Figure 7). The calibration plots demonstrated good agreement between predicted probabilities and observed outcomes in both the internal test set and temporal validation cohort (Figure 8), with Brier scores of 0.076 and 0.096, respectively.

Figure 6

Figure 7

Figure 8

Decision curve analysis showed that the combined model provided superior net benefit across a wide range of threshold probabilities (20% to 90%) in both internal and external cohorts (Figure 9), supporting its robust clinical utility.

Figure 9

Discussion

We developed and validated a multivariable prediction model for csPCa integrating clinical variables, mpMRI parameters, and 18F-PSMA PET/CT findings. In the training set, the model achieved an apparent AUC of 0.916. It demonstrated excellent discriminative ability in the internal test set (AUC 0.914, 95% CI: 0.882-0.941) and maintained robust transportability in the temporal validation cohort (AUC 0.837, 95% CI: 0.778-0.891).

While mpMRI represents the current standard for pre-biopsy risk stratification, it is primarily limited by moderate interobserver variability and suboptimal specificity. The integration of 18F-PSMA PET/CT provides complementary molecular insight, specifically mapping the overexpression of prostate-specific membrane antigen in clinically significant lesions (). Recent evidence, such as the PRIMARY2 trial, highlights the critical role of PSMA PET strategies in safely reducing unnecessary biopsies. Furthermore, there is an evolving clinical landscape toward utilizing quantitative PET parameters, such as SUVmax, for robust risk stratification and direct therapeutic decision-making (). Our integrated approach further demonstrates that combining PSMA PET/CT, mpMRI, and PSAD optimizes risk stratification beyond what either modality can achieve alone. In our rigorous dual-stage variable screening, PRIMARY score demonstrated profound independent prognostic value, whereas quantitative SUVmax and MRI-based T-stage were eliminated during multivariable stepwise adjustment. This exclusion underscores a critical principle of statistical parsimony: the variance explained by SUVmax and T-stage was highly collinear with the combined qualitative assessment provided by PI-RADS and PRIMARY scores. Consequently, our 3-variable model avoids the redundancy of including overlapping imaging metrics, offering a streamlined yet powerful tool for clinical decision-making. It is worth noting that, based on our nomogram, the isolated predictive probability of a PRIMARY score of 5 appears somewhat lower than the >90% detection rate reported in the original PRIMARY publication (). This apparent reduction is an inherent characteristic of our multivariable model design. Because the overall predictive weight is systematically distributed across three distinct clinical and imaging variables (PRIMARY score, PI-RADS score, and PSAD), the maximum predictive threshold can only be reached when these complementary parameters are combined. Therefore, a high PRIMARY score in our model should not be interpreted in isolation, but rather as a synergistic component of the integrated approach, which ultimately yields a more comprehensive risk stratification.

Furthermore, our final model strictly retained prostate-specific antigen density (PSAD) as the sole clinical biomarker, actively excluding raw serum PSA and prostate volume to circumvent computational collinearity. PSAD has been extensively validated in contemporary literature as a superior metric to PSA, effectively neutralizing the confounding effect of benign prostatic hyperplasia and significantly improving diagnostic specificity for csPCa (, ). Our stepwise AIC selection objectively confirmed that PSAD provides optimal, non-redundant predictive weight alongside advanced imaging parameters.

To maximize clinical applicability, we evaluated our nomogram across two distinct probability thresholds. The Youden Index identified an optimal cutoff of ≥84%, which balances sensitivity and specificity for confident csPCa diagnosis. Conversely, selecting a lower threshold of ≥46% drastically elevates sensitivity to 96.0%. This lower cutoff is clinically recommended for screening purposes, aligning with the primary objective of pre-biopsy risk stratification: safely avoiding unnecessary biopsies without missing high-risk cancers, analogous to the triage thresholds employed by contemporary risk calculators ().

The observed reduction in AUC from 0.914 in the internal test set to 0.837 in the temporal validation cohort represents a natural calibration drift commonly encountered when prediction models are strictly transported to temporally or geographically distinct populations (, ). Contributing factors likely include differences in baseline cohort characteristics—such as variations in csPCa prevalence and median PSAD—as well as evolving referral patterns over time. Additionally, this study has limitations, including its retrospective nature which may introduce selection bias, and the relatively small sample size in the validation cohort. Future prospective multicenter trials are warranted to further refine these thresholds and establish broad generalizability.

Currently, several pre-biopsy risk calculators are utilized in clinical practice to assist decision-making. For instance, the Prostate Cancer Prevention Trial (PCPT) risk calculator, incorporating clinical variables such as age and PSA, yields an area under the curve (AUC) of approximately 0.70 to 0.75 (). Similarly, the European Randomized Study of Screening for Prostate Cancer (ERSPC) risk calculator demonstrates an AUC of 0.75 to 0.80 (). More advanced tools, such as the Stockholm3 model, which integrates clinical variables with protein markers and genetic polymorphisms, can achieve AUCs of 0.80 to 0.85, albeit at a significantly higher diagnostic cost (). Compared to these existing risk models, a key innovation of our study is the formal integration of quantitative and qualitative PSMA PET parameters alongside standard clinical and MRI factors, ultimately yielding robust diagnostic discrimination (temporal validation AUC 0.837) without the requisite for complex genetic profiling.

Several limitations should be acknowledged. First, despite consecutive enrollment, the retrospective component may be subject to selection bias. Second, the reference standard was systematic biopsy rather than whole-mount prostatectomy pathology, which may underestimate true tumor burden (). Third, the prevalence of csPCa in our cohort is remarkably high (78.3% and 84.4%). This reflects a tertiary-center spectrum bias, as our population primarily consists of patients referred with highly elevated PSA levels or suspicious prior imaging. This bias may overestimate the positive predictive value and limits the direct applicability of our model to primary screening populations; the nomogram is primarily intended to aid decision-making in high-risk populations where multimodal imaging is already being considered. Fourth, MRI-based T staging rather than DRE-based clinical T staging was used (). Fifth, despite temporal validation, further validation in larger, multi-centric, and prospective settings is necessary prior to wide-scale clinical adoption. Furthermore, the combined utilization of mpMRI and 18F-PSMA PET/CT imposes a significant financial burden and relies on advanced medical infrastructure. We acknowledge that this integrated approach may be cost-prohibitive and less accessible in developing countries or resource-limited settings. Consequently, while our model demonstrates excellent theoretical diagnostic performance, its immediate practical application may be restricted to well-resourced tertiary centers rather than routine primary screening. Finally, while all scans were reviewed by highly experienced readers, we did not formally calculate inter-observer agreement metrics (such as Cohen’s Kappa) at the time of study design, which represents a potential source of bias ().

Future studies should focus on prospective multicenter validation, evaluation in populations with lower csPCa prevalence, and incorporation of emerging biomarkers. The model may serve as a framework for integrating multimodal imaging data in prostate cancer diagnosis ().

Conclusion

We developed and validated a multivariable prediction model for csPCa integrating clinical variables, mpMRI parameters, and 18F-PSMA PET/CT findings. Compared to existing risk models, a key innovation of our study is the formal integration of quantitative and qualitative PSMA PET parameters alongside standard clinical and MRI factors. The model achieved an AUC of 0.916 in the training set, with maintained discriminative ability in both internal test (AUC = 0.914) and temporal validation (AUC = 0.837) cohorts.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by Institutional Ethics Committee of Peking University First Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

CY: Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft. ZZ: Data curation, Investigation, Validation, Writing – original draft. PT: Data curation, Investigation, Writing – original draft. JZ: Data curation, Formal analysis, Software, Writing – original draft. LC: Data curation, Formal analysis, Writing – original draft. JQ: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing. KG: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the National High Level Hospital Clinical Research Funding (Scientific and Technological Achievements Transformation Incubation Guidance Fund Project of Peking University First Hospital, 2024CX23; High Quality Clinical Research Project of Peking University First Hospital, 2022CR75), National Natural Science Foundation of China (No. 82141103; 82103153; 82172617; 82172665; 82103153), Beijing Natural Science Foundation (QY23068; No. 7232176), and Capital’s Funds for Health Improvement and Research (2022-2-4074).

Acknowledgments

We would like to give our sincere appreciation to all the participants in our cohort and their families for their support.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2026.1835850/full#supplementary-material

References

  • 1

    SiegelRLKratzerTBGiaquintoANSungHJemalA. Cancer statistics, 2025. Ca. (2025) 75:10. doi: 10.3322/caac.21871. PMID:

  • 2

    AndrioleGLCrawfordEDGrubbRL3rdBuysSSChiaDChurchTRet al. Mortality results from a randomized prostate-cancer screening trial. N Engl J Med. (2009) 360:1310–9. doi: 10.1056/nejmoa0810696. PMID:

  • 3

    EpsteinJIEgevadLAminMBDelahuntBSrigleyJRHumphreyPAet al. The 2014 International Society of Urological Pathology (ISUP) consensus conference on Gleason grading of prostatic carcinoma: definition of grading patterns and proposal for a new grading system. Am J Surg Pathol. (2016) 40:244–52. doi: 10.1097/pas.0000000000000530. PMID:

  • 4

    LoebSBjurlinMANicholsonJTammelaTPensonDFCarterHBet al. Overdiagnosis and overtreatment of prostate cancer. Eur Urol. (2014) 65:1046–55. doi: 10.1016/j.eururo.2013.12.062. PMID:

  • 5

    TurkbeyBRosenkrantzABHaiderMAPadhaniARVilleirsGMacuraKJet al. Prostate imaging reporting and data system version 2.1: 2019 update of prostate imaging reporting and data system version 2. Eur Urol. (2019) 76:340–51. doi: 10.1016/j.eururo.2019.02.033. PMID:

  • 6

    AhmedHUEl-Shater BosailyABrownLCGabeRKaplanRParmarMKet al. Diagnostic accuracy of multi-parametric MRI and TRUS biopsy in prostate cancer (PROMIS): a paired validating confirmatory study. Lancet. (2017) 389:815–22. doi: 10.1016/s0140-6736(16)32401-1. PMID:

  • 7

    KasivisvanathanVTakwoingiYEmbertonMMooreCM. MRI-targeted or standard biopsy for prostate-cancer diagnosis. N Engl J Med. (2018) 378:1767–77. doi: 10.1056/nejmoa1801993. PMID:

  • 8

    PadhaniARBarentszJVilleirsGRosenkrantzABMargolisDJTurkbeyBet al. PI-RADS steering committee: the PI-RADS multiparametric MRI and MRI-directed biopsy pathway. Radiology. (2019) 292:464–74. doi: 10.1148/radiol.2019182946. PMID:

  • 9

    BensonMCWhangISPantuckARingKKaplanSAOlssonCA. Prostate specific antigen density: a means of distinguishing benign prostatic hypertrophy and prostate cancer. J Urol. (1992) 147:815–6. doi: 10.1016/s0022-5347(17)37393-7. PMID:

  • 10

    FendlerWPEiberMBeheshtiMBomanjiJCeciFChoSet al. 68Ga-PSMA PET/CT: joint EANM and SNMMI procedure guideline for prostate cancer imaging: version 1.0. Eur J Nucl Med Mol Imaging. (2017) 44:1014–24. doi: 10.1007/s00259-017-3670-z. PMID:

  • 11

    GieselFLHadaschikBCardinaleJRadtkeJVinsensiaMLehnertWet al. F-18 labelled PSMA-1007: biodistribution, radiation dosimetry and histopathological validation of tumor lesions in prostate cancer patients. Eur J Nucl Med Mol Imaging. (2017) 44:678–88. doi: 10.1007/s00259-016-3573-4. PMID:

  • 12

    RenYLiuCLiuTDuanXZhangQLiuJet al. Preclinical evaluation and first in human study of Al18F radiolabeled ODAP-urea-based PSMA targeting ligand for PET imaging of prostate cancer. Front Oncol. (2022) 12:1030187. doi: 10.3389/fonc.2022.1030187. PMID:

  • 13

    EmmettLPapaNButeauJHoBLiuVRobertsMet al. The PRIMARY score: using intraprostatic 68Ga-PSMA PET/CT patterns to optimize prostate cancer diagnosis. J Nucl Med. (2022) 63:1644–50. doi: 10.2967/jnumed.121.263448. PMID:

  • 14

    CornfordPvan den BerghRCNBriersEVan den BroeckTBrunckhorstODarraughJet al. EAU-EANM-ESTRO-ESUR-ISUP-SIOG guidelines on prostate cancer—2024 update. Part I: screening, diagnosis, and local treatment with curative intent. Eur Urol. (2024) 86:148–63. doi: 10.1016/j.eururo.2024.03.027. PMID:

  • 15

    TibshiraniR. Regression shrinkage and selection via the lasso. J R Stat Soc Ser B Stat Method. (1996) 58:267–88. doi: 10.1111/j.2517-6161.1996.tb02080.x. PMID:

  • 16

    MetzCE. Basic principles of ROC analysis. Semin Nucl Med. (1978) 8(4):283–98. doi: 10.1016/s0001-2998(78)80014-2

  • 17

    Van CalsterBMcLernonDJvan SmedenMWynantsLSteyerbergEWSTRATOS initiativeet al. Calibration: the Achilles heel of predictive analytics. BMC Med. (2019) 17:230. doi: 10.1186/s12916-019-1466-7. PMID:

  • 18

    VickersAJElkinEB. Decision curve analysis: a novel method for evaluating prediction models. Med Decision Making. (2006) 26:565–74. doi: 10.1177/0272989x06295361. PMID:

  • 19

    McKinneyW. Data structures for statistical computing in Python. Scipy. (2010) 445:51–6. doi: 10.25080/majora-92bf1922-00a

  • 20

    ChengCLiuJYiXYinHQiuDZhangJet al. Prediction of clinically significant prostate cancer using a novel 68Ga-PSMA PET-CT and multiparametric MRI-based model. Trans Andrology Urol. (2023) 12:1115. doi: 10.21037/tau-22-832. PMID:

  • 21

    EiberMHerrmannKCalaisJHadaschikBGieselFLHartenbachMet al. Prostate cancer molecular imaging standardized evaluation (PROMISE): proposed miTNM classification for the interpretation of PSMA-ligand PET/CT. J Nucl Med. (2018) 59:469–78. doi: 10.2967/jnumed.117.198119. PMID:

  • 22

    RodríguezDMAlmazán TreviñoLLamm WiechersLGonzález CosióRGalicia BelauzaranLFHerrera MuñozJAet al. Radical prostatectomy without prior biopsy in patients with positive PSMA PET and/or mpMRI: A practiced but poorly documented approach. Retrospective multicenter cohort study. Actas Urológicas Españolas (English Edition). (2026), 501921. doi: 10.1016/j.acuroe.2026.501921

  • 23

    SiddiquiMMRais-BahramiSTurkbeyBGeorgeAKRothwaxJShakirNet al. Comparison of MR/ultrasound fusion–guided biopsy with ultrasound-guided biopsy for the diagnosis of prostate cancer. Jama. (2015) 313:390–7. doi: 10.1001/jama.2014.17942. PMID:

  • 24

    SteyerbergEWHarrellFE. Prediction models need appropriate internal, internal-external, and external validation. J Clin Epidemiol. (2015) 69:245–55. doi: 10.1016/j.jclinepi.2015.04.005. PMID:

  • 25

    ThompsonIMAnkerstDPChiCGoodmanPJTangenCMLuciaMSet al. Assessing prostate cancer risk: results from the Prostate Cancer Prevention Trial. J Natl Cancer Institute. (2006) 98:529–34. doi: 10.1093/jnci/djj131. PMID:

  • 26

    RoobolMJSteyerbergEWKranseRWoltersTvan den BerghRCBangmaCHet al. A risk-based strategy improves prostate-specific antigen–driven detection of prostate cancer. Eur Urol. (2010) 57:7985. doi: 10.1016/j.eururo.2009.08.025. PMID:

  • 27

    GrönbergHAdolfssonJAlyMNordströmTWiklundPBrandbergYet al. Prostate cancer screening in men aged 50–69 years (STHLM3): a prospective population-based diagnostic study. Lancet Oncol. (2015) 16:1667–76. doi: 10.1016/S1470-2045(15)00361-7

  • 28

    AhdootMWilburARReeseSELebastchiAHMehralivandSGomellaPTet al. MRI-targeted, systematic, and combined biopsy for prostate cancer diagnosis. N Engl J Med. (2020) 382:917–28. doi: 10.1056/nejmoa1910038. PMID:

  • 29

    TayaMBehrSCWestphalenAC. Perspectives on technology: Prostate Imaging‐Reporting and Data System (PI‐RADS) interobserver variability. BJU Int. (2024) 134:510–8. doi: 10.1111/bju.16452

  • 30

    KalaparaAANzenzaTPanHYCBallokZRamdaveSO'SullivanRet al. Detection and localisation of primary prostate cancer using 68gallium prostate‐specific membrane antigen positron emission tomography/computed tomography compared with multiparametric magnetic resonance imaging and radical prostatectomy specimen pathology. BJU Int. (2020) 126:8390. doi: 10.1111/bju.14858. PMID:

Summary

Keywords

clinically significant prostate cancer, multiparametric MRI, prediction model, prostate cancer, PSMA PET/CT

Citation

Yu C, Zhao Z, Tian P, Zhou J, Cai L, Qiu J and Gong K (2026) A multivariable prediction model combining 18F-PSMA PET/CT and mpMRI for clinically significant prostate cancer: development and validation. Front. Oncol. 16:1835850. doi: 10.3389/fonc.2026.1835850

Received

21 March 2026

Revised

23 April 2026

Accepted

29 April 2026

Published

18 May 2026

Volume

16 - 2026

Edited by

Carmelo Caldarella, Fondazione Policlinico Universitario A. Gemelli IRCCS, Italy

Reviewed by

Jorge Gustavo Morales Montor, Hospital Gea Gonzalez, Mexico

Michelangelo Cobangbang, St. Luke’s Medical Center, Philippines

Updates

Copyright

*Correspondence: Jianhui Qiu,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics