MRI-Based Nomogram of Prostate Maximum Sectional Area and Its Zone Area for Prediction of Prostate Cancer

Objective: To reduce unnecessary prostate biopsies, we designed a magnetic resonance imaging (MRI)-based nomogram prediction model of prostate maximum sectional area (PA) and investigated its zone area for diagnosing prostate cancer (PCa).

Methods: MRI was administered to 691 consecutive patients before prostate biopsies from January 2012 to January 2020. PA, central gland sectional area (CGA), and peripheral zone sectional area (PZA) were measured on axial T2-weighted prostate MRI. Multivariate logistic regression analysis and area under the receiver operating characteristic (ROC) curve were performed to evaluate and integrate the predictors of PCa. Based on multivariate logistic regression coefficients after excluding combinations of collinear variables, three models and nomograms were generated and intercompared by Delong test, calibration curve, and decision curve analysis (DCA).

Results: The positive rate of PCa was 46.74% (323/691). Multivariate analysis revealed that age, PSA, MRI, transCGA, coroPZA, transPA, and transPAI (transverse PZA-to-CGA ratio) were independent predictors of PCa. Compared with no PCa patients, transCGA (AUC = 0.801) was significantly lower and transPAI (AUC = 0.749) was significantly higher in PCa patients. Both of them have a significantly higher AUC than PSA (AUC = 0.714) and PV (AUC = 0.725). Our best predictive model included the factors age, PSA, MRI, transCGA, and coroPZA with the AUC of 0.918 for predicting PCa status. Based on this predictive model, a novel nomogram for predicting PCa was conducted and internally validated (C-index = 0.913).

Conclusions: We found the potential clinical utility of transCGA and transPAI in predicting PCa. Then, we firstly built the nomogram based on PA and its zone area to evaluate its diagnostic efficacy for PCa, which could reduce unnecessary prostate biopsies.

Introduction

Prostate cancer (PCa) is the most common cancer among men in the Western world, and it has an increasing prevalence (1). There is an international consensus that early detection and treatment of PCa can improve the survival rate of PCa patients. Prostate-specific antigen (PSA) is the most widely used screening marker to detect PCa at an early stage. The larger clinical trial found that patients having undergone PSA screening had 25% lower PCa death rates than those who did not (2). After tests reveal an elevated serum PSA level, most patients require puncture biopsy of the prostate, because the prostate biopsy remains the gold standard method for diagnosing PCa. However, we have to face a clinical problem that the prostate biopsy is an invasive operation. It not only brings pain and fear to the patients, but also may cause medical complications such as infection and hemorrhage (3). Because prostate biopsy always has the probability of missing tumor tissue, it is not able to make a 100% diagnosis of PCa. The rate of negative prostate biopsies was substantially high (58.51%–69.30%) especially in cases with only elevated PSA levels, thus greatly affecting patients’ quality of life (4, 5). Therefore, it is rational to avoid the biopsy on patients who are ultimately proved to be negative cases.

In order to overcome the limitations of PSA test, Benson et al. proposed the concept of PSA density (PSAD, PSA value divided by prostate volume), which was considered to increase the accuracy of PSA test for diagnosing the PCa (6). The main principle is that PCa tissues can release more PSA per unit volume to blood serum than enlarged or normal prostate tissues do. Recent research had also shown that PSA density could overweigh PSA in distinguishing clinically significant PCa and intraprostatic inflammation before prostate biopsy (7). However, it has been reported that the prostate volume was frequently roughly calculated using the prolate ellipsoid formula before operation, in which there is 10%–20% error compared with prostatectomy specimens in the clinical situation (8, 9). So, this kind of prostate volume should be further improved for assessing the exact risk of PCa. Therefore, we considered finding new MRI-based predictors to enhance the role of roughly calculated prostate volume for predicting PCa.

In MRI images of prostate zonal anatomy, the prostate comprises the peripheral zone, transition zone, central zone, and anterior fibromuscular stroma (10). A related study showed that approximately 75%–85% of PCa cases are located in the peripheral zone, rather than the central gland. The central gland is the typical site of BPH, which includes the transition zone, central zone, and anterior fibromuscular stroma (11). MRI images can clearly identify different anatomical areas of the prostate, which is beneficial to improve the detection rate of PCa. We first propose the concept of the prostate maximum sectional area (PA) for predicting PCa by MRI images, which includes both central gland sectional area (CGA) and peripheral zone sectional area (PZA). Comparatively, MRI is regarded as the most precise noninvasive method, as it can assess PA with high reproducibility and accuracy compared with rough prostate volume calculated by the common formula (12).

Nomogram is a simple intuitive graph of a complex mathematical formula (13). It is widely used for cancer prediction, primarily because of their ability to use biologic and clinical variables building a graphically depictive statistical predictive model that is tailored to an individual patient (14). User-friendly graphical interfaces for generating these estimates facilitate the use of nomograms to aid in clinical decision-making.

The purpose of the current study was to establish a new nomogram about PA and its associated zone area such as CGA and PZA on axial T2 fat-saturated MRI for diagnosing PCa. To the best of our knowledge, no previous literature has employed MRI-based PA and the associated zone area for the prediction of PCa before prostate biopsy among Chinese population.

Methods

Study Population

The study was designed as a retrospective cohort study that was conducted in the Laboratory of Urology and the Department of Urology of Fujian Medical University Union Hospital (Fuzhou, China). We enrolled 691 consecutive patients who underwent multiparametric magnetic resonance image (mp-MRI) before initial transrectal ultrasound (TRUS)-guided prostate biopsy from January 2012 to January 2020, followed by data anonymization. Eligible patients who matched the selection criteria were identified by the following criteria: elevated PSA levels (≥10 ng·ml⁻¹), suspected cancer on digital rectal examination (DRE), hyperechoic or hypoechoic TRUS, or abnormal MRI findings. For PSA values between 4 and 10 ng·ml⁻¹, the biopsy criterion was ratio of free to total PSA < 16%. The exclusion criteria were as follows: previous prostate biopsy, history of prostate surgery, pathological examination revealing tumors other than adenocarcinoma, and incomplete mp-MRI information or imaging artifacts (Figure 1). The study was ethically approved by the Institutional Review Board of Fujian Medical University Union Hospital with an approval number of 2020KY059. Written informed consent was obtained from patients before the study commenced. Details of patients’ identity had to be omitted. Our work complies with the Code of Ethics of the World Medical Association (Declaration of Helsinki, revised in 2013).

FIGURE 1

Figure 1 Flow chart of patient selection.

Clinical Data and Variable Definitions

Clinical characteristics including age, body mass index (BMI), PSA, free PSA (FPSA), free-to-total PSA (FTPSA), prostate volume (PV), PSA density (PSAD), MRI, transverse prostate maximum sectional area (transPA), coronal prostate maximum sectional area (coroPA), sagittal prostate maximum sectional area (sagiPA), transverse peripheral zone sectional area (transPZA), coronal peripheral zone sectional area (coroPZA), sagittal peripheral zone sectional area (sagiPZA), transverse central gland sectional area (transCGA), coronal central gland sectional area (coroCGA), sagittal central gland sectional area (sagiCGA), alkaline phosphatase (ALP), and lactate dehydrogenase (LDH) were collected before prostate biopsy. Subsequently, transverse PSA-to-PA ratio (transPSAPA), coronal PSA-to-PA ratio (coroPSAPA), sagittal PSA-to-PA ratio (sagiPSAPA), transverse PSA-to-PZA ratio (transPSAPZA), coronal PSA-to-PZA ratio (coroPSAPZA), sagittal PSA-to-PZA ratio (sagiPSAPZA), transverse PSA-to-CGA ratio (transPSACGA), coronal PSA-to-CGA ratio (coroPSACGA), sagittal PSA-to-CGA ratio (sagiPSACGA), transverse PZA-to-CGA ratio (transverse prostate area index, transPAI), coronal PZA-to-CGA ratio (coronal prostate area index, coroPAI), and sagittal PZA-to-CGA ratio (sagittal prostate area index, sagiPAI) were calculated. The prostate maximum sectional area on prostate T2WI MRI had the following definition: In the transverse plane, when the bilateral prostate lobes are basically symmetrical, and the quasi-circular internal urethral sphincter can be seen in the middle of the prostate, the maximum section is the one for which the sectional area becomes smaller when scanning upward or downward. In the coronal plane, when the bilateral prostate lobes are basically symmetrical, and the strip-type internal urethral sphincter can be seen in the middle of the prostate, the maximum section is the one for which the sectional area becomes smaller when scanning upward or downward. In the sagittal plane, when the strip-type internal urethral sphincter can be seen in the middle of the prostate, the maximum section is the one for which the sectional area becomes smaller when scanning upward or downward.

Image Acquisition and Interpretation

A Siemens Magnetom Trio Tim 3.0-T superconducting MRI scanner with an 18-channel phased-array torso coil was used to create all magnetic resonance images [repetition time (TR) 400 ms, echo time (TE) 80 ms, slice thickness = 3 mm, interslice gap = 30%, acquisition four times with fat-suppression technique]. T2-weighted images in the sagittal, coronal, and transverse planes, diffusion-weighted images, apparent diffusion coefficient in the transverse plane, and dynamic contrast-enhanced images were acquired according to the international prostate MRI guidelines (15). Interpretation of the MRI findings was performed by a radiologist and a urologist (with 5 or more years of experience in prostate imaging),who measured PA and CGA on fat-saturated T2WI MRI (Figure S1).

Prostate Biopsy Method

Following local non-infiltrative anesthesia, all prostate biopsies were performed transrectally under TRUS guidance (BK Medical, USA). A standard 13-core systematic prostate biopsy was obtained including transitional, peripheral, and anterior zone from base to apex by an 18-gauge/25-cm biopsy needle (Bard Peripheral Vascular, Inc). All patients underwent standard prostate biopsies, which were performed by an experienced urologist (more than 5 years of experience in prostate biopsy). All biopsy specimens were examined and recorded by two experienced pathologists.

Statistical Analysis

Distributions of variables were compared by the chi-squared test for categorical variables and the Mann–Whitney U test for continuous variables, which was not normally distributed. The values of all continuous variables (age, BMI, PSA, FPSA, FTPSA, transPA, coroPA, sagiPA, transPZA, coroPZA, sagiPZA, transCGA, coroCGA, sagiCGA, transPSAPA, coroPSAPA, sagiPSAPA, transPSAPZA, coroPSAPZA, sagiPSAPZA, transPSACGA, coroPSACGA, sagiPSACGA, ALP, and LDH) were not normally distributed. Variables including BMI, FPSA, FTPSA, transPZA, and ALP were excluded due to lack of statistical significance in univariate logistic regression analysis. We integrated variables with great clinical significance including age, PSA, and MRI into the base model. The remaining variables were reassembled into all kinds of possible combinations through enumeration algorithm. Then, we combined base model and each different combinations together to form our predict models. Correlation analysis was used to detect the multicollinearity between every two variables (Figure S2). Any model that contained two or more multicollinearity variables will be eliminated before the next step. Then, multivariate logistic regression analysis was performed on the rest of these models to identify the independence of each predictor for diagnosing PCa and calculate its variance inflation factor (VIF). Models will also be eliminated when their VIF ≥ 2. The diagnostic efficacy of these models was evaluated by the area under the curve (Figure 2). The first three combinations with the highest AUC were chosen as our final models. The statistical differences among three models and each single predictor were compared by Delong test, respectively. The cutoff value, sensitivity, specificity, and positive and negative likelihood ratios were computed for these variables and prediction models. Nomograms were generated to predict the probability of PCa, based on the multivariate regression coefficients in three models. These models were recalibrated both in the training cohort and the validation cohort to evaluate the nomogram’s discrimination capacity by 1,000 random bootstrap samples with replacement. Calibration slope less than 1 reflects proper fit of the model. The clinical utility of three models was quantified by decision curve analysis (DCA) through summing the benefits (true positives) and subtracting the harms (false positives). Statistical significance was defined as p-value < 0.05. Statistical analysis, nomogram, and calibration plot were generated using R studio (version 4.0.3).

FIGURE 2

Figure 2 Flow chart of statistical analysis.

Results

Clinical Characteristics

A total of 230 (47.23%) of 487 patients in the training cohort and 93 (45.59%) of 204 patients in the validation cohort were diagnosed with PCa (Table 1). Univariate logistic regression analysis showed that all variables were statistically significant predictors of PCa detection except for BMI, FPSA, FTPSA, transPZA, and ALP in the training cohort. No significance was found in variables between the training cohort and validation cohort except for age (Table 2).

TABLE 1

Table 1 Clinical characteristics of patients before the prostate biopsy.

TABLE 2

Table 2 Clinical characteristics of patients in training cohort and validation cohort.

Multivariate Logistic Regression Models

To evaluate the synergistic ability of every single predictor for predicting PCa, we created different models that did not contain multicollinearity. The first three models with the highest AUC were chosen as our final models, which were model 1, model 2, and model 3. Model 1 consists of age, PSA, MRI, transCGA, and coroPZA after excluding sagiPAI, PV, and PSAD. Model 2 consists of age, PSA, MRI, transPAI, coroPZA, and transPA after excluding PV and PSAD. Model 3 consists of age, PSA, MRI, transPAI, and PV after excluding sagiPAI, coroPAI, and PSAD (Table 3). Age, MRI, and PSA were independent predictors of PCa in three models. The statistical significance of PV was not detected in both model 1 and model 2. The statistical significance of PSAD was not detected in all models.

TABLE 3

Table 3 Multivariate logistic regression analysis of predictors associated with PCa before the prostate biopsy.

Comparison of Predictive Accuracy

Furthermore, predictive accuracy of each predictor alone and models were assessed separately with ROC curve analysis. The AUC of model 1 (base model + coroPZA + transCGA) for predicting PCa was the highest among the models of any single predictor alone and the base model combined with any other predictors (Figure 3 and Table 4). Delong test was used to compare the statistical difference of the AUC among three models and single predictors. Compared with model 3 (AUC = 0.907), the AUC of model 1 (AUC = 0.918) and model 2 (AUC = 0.916) had the same higher statistical advantages in the training cohort, while no statistical difference was found among the three models in the validation cohort (Table 5). The AUC of transCGA (0.801) was significantly higher than other single predictors. The AUC of coroPZA (0.635) was lower than other predictors, while there was no significant difference of AUC among transPAI, transPA, PV, and PSA in the training cohort (Table 6).

FIGURE 3

Figure 3 Receiver operating characteristic curves depicting the accuracy of predictors of PCa before the initial biopsy. Base model: age + PSA + MRI. Model 1: Base model + coroPZA + transCGA. Model 2: Base model + transPAI + coroPZA + transPA. Model 3: Base model + transPAI + PV.

TABLE 4

Table 4 The AUC and cutoff values for predicting biopsy outcome and their sensitivity, specificity, and positive and negative likelihood ratios for PCa.

TABLE 5

Table 5 The statistical difference in AUC of predicting PCa among three models.

TABLE 6

Table 6 The statistical difference in AUC of predicting PCa among single predictors in the training cohort.

Nomograms, Calibration Plots, and DCA Curves

Based on the multivariate regression coefficients, the predictive models were visually presented as nomograms (Figure 4 and Figure S3). The nomogram’s discrimination of three models in the training cohort and validation cohort was shown in the calibration plot (Figure 5, Figures S4 and S5). The C-index of model 1 for predicting PCa was 0.918 in the training cohort. The performance of model 2 (0.916) and model 3(0.907) in the calibration plot was not as good as that of model 1 in the training cohort, which demonstrates the superior fit of model 1. DCA curves showed that the nomogram based on model 1 has better net benefit gains in all range of threshold probabilities in the training cohort, while net benefit gains only improved when threshold probabilities are >8% in the validation cohort (Figure 6).

FIGURE 4

Figure 4 Nomogram predicting the probability of PCa at the initial biopsy based on model 1.

FIGURE 5

Figure 5 Calibration plot in training cohort and validation cohort and predictive accuracy for PCa at initial biopsy based on model 1.

FIGURE 6

Figure 6 Decision curve analysis of the effect of the nomogram based on model 1 for predicting prostate cancer in training cohort and validation cohort. Net benefit of nomogram is plotted with threshold probabilities for prostate cancer compared with the strategies of treating all patients or no one. The decision curve illustrated net benefit was improved when threshold probability > 8%.

Discussion

PCa is a malignant form of cancer whose diagnosis depends on the histopathological verification of adenocarcinoma in a prostate biopsy. However, excess of prostate biopsy has led to increased side effects such as bleeding and infection. It also caused the inferior positive rate of 30%–40% (3, 5). A retrospective study including 1,203 patients who underwent prostate biopsy demonstrated that the overall rates of infectious and hemorrhagic complications after prostate biopsy were 8.23% and 15.71%, respectively (16). So, it needs to establish a method to carefully select patients who need prostate biopsy.

Numerous studies have reported the predictive value of prostate volume (PV) and prostate volume-adjusted PSA (i.e., PSAD) for PCa. For example, one study that measured 235 patients’ prostate volume and PSA levels demonstrated that the AUC values of PSAD (0.712) and prostate volume (0.710) were higher than that of PSA (0.517) for diagnosing PCa (17). Our previous research found that the utility of PSAD for performing surveillance in patients at risk of PCa was higher than that of standard variables such as PSA (18). However, a retrospective study found that PSAD and PSA (AUC = 0.620 and 0.530, respectively) failed to outperform prostate volume (AUC = 0.680) for preoperative prediction of PCa (19). The current study confirmed that both PV and PSAD were good predictors of PCa in univariate logistic analysis. However, none of them showed statistical significance in model 1 and model 2. The reason is possibly that the prostate is not a regular geometric solid especially in the malignant growth mode of the tumor. Furthermore, prostate volume is usually estimated by an elliptical sphere formula (PV = 0.52 × length × width × height). Any error on the length, width, or height of prostate may be magnified through the multiplication (20). Previous studies had shown a high magnitude of bias between the calculation of prostate volume by the prolate ellipsoid formula and the actual prostate volume, which casts doubt on the diagnostic efficacy of PV and PSAD in PCa (21). Some studies had confirmed that the bias of calculated prostate volume fluctuates between 10% and 20% (8, 9). As we mentioned above, the increase of PSA level caused by any other reason, except for PCa, may lead to the error of PSAD and reduce its specificity for predicting PCa. Therefore, we do not think that PV and PSAD have the leading advantage of predicting PCa.

In order to overcome the difficulties mentioned above, we sought to replace the role of prostate volume with the incorporation of more accurate, simple prostatic imaging parameters. We found that PA (prostate maximum sectional area) is a good prostatic imaging parameter for predicting PCa by MRI test in line with the above requirements. MRI has higher spatial resolution and better soft tissue contrast than TRUS, and MRI can provide more accurate PA. So, it can reflect the actual size of the prostate (19). Thus, we used MRI-based PA as an alternative predictor of prostate volume for predicting PCa. As far as we know, we are the first to use prostate maximum sectional areas in sagittal, transverse, and coronal directions to predict PCa. The data of prostate sectional area from three different directions may help to improve its representativeness for irregular prostate. It can possibly find out the shape characteristics of the prostate in different directions. On the other hand, the area and its zone area of prostate were actually measured in MRI segments, which will decrease the systematic error to a great extent compared with the calculated prostate volume by formula. In our research, all sectional area predictors have statistical differences between PCa and no PCa patients in univariate analysis, except for transPZA. It proved that they had great potential in predicting PCa. We found that BPH patients had larger PA and CGA but smaller PZA in three directions compared with those who had PCa. It might result from the fact that BPH contributes to mechanical stress fields by pathological enlargement of the prostate central gland, hence further restraining PCa growth, as PCa mostly originates in the peripheral zone of the prostate (22). Through Delong test among single predictors, we find out the transCGA has the significantly highest AUC (0.801) among all predictors. Compared with no PCa, transCGA is significantly smaller in patients with PCa (p < 0.001). We speculated that this is due to the special behavior pattern of PCa growth on the transverse section. So far, we have not seen any relevant report that needs to be confirmed by further pathological or anatomical studies.

To further explore the potential of prostate sectional area related predictors to predict PCa, we built new predictors like PAI, PSAPA, PSACGA, and PSAPZA (calculated by PZA/CGA, PSA/PA, PSA/CGA, and PSA/PZA) based on the sectional area from three different directions in MRI segments. Each of these predictors had the potential ability to distinguish PCa from no PCa in univariate logistic analysis. Unfortunately, we found that PSAPA, PSACGA, and PSAPZA failed to outperform any predictors when we discovered the different models. So, we did not include these predictors in multivariate logistic analysis. However, PAI (prostate area index) showed favorable predictability in our final model, especially transPAI, which had an AUC (0.749) second only to transCGA. In both model 2 and model 3, transPAI had a certain contribution to the diagnosis of PCa compared with coroPAI and sagiPAI. This also confirms our previous hypothesis about the special behavior pattern of PCa growth on the transverse section. Patients with PCa have higher transPAI. It may be due to PCa often originating in the peripheral zone, which causes the enlargement of PZA. Then, it leads to an increase in PAI (PZA-to-CGA ratio) in turn.

With the changes to people’s living habits, and the progression of population aging, the prevalence of PCa is increasing annually, especially among adults aged over 70 years. This situation has seriously affected the health of older adults (23). Our results also confirm this conclusion. The median age of men with PCa was 75 years, compared with 69 years for men without PCa (p < 0.001). Thus, age can be used as a reference for prostate biopsy.

To reduce unnecessary prostate biopsy and improve the diagnostic accuracy of PCa in clinical practice, nomograms integrating many independent predictors of PCa have been developed and validated. A previous study reported that nomograms could provide more individualized risk estimations of a certain disease, which could help clinicians to make management-related decisions for patients with PCa (24). For example, a nomogram developed on the basis of 1,144 men who underwent TRUS found that the C-index (0.876) was associated with their best model that integrates age, PSA, percentage free PSA, DRE, prostate transition zone volume, and TRUS for predicting PCa (25). Another study integrated age, prostate volume, PSA, FTPSA, TRUS, and DRE as its best model to develop a nomogram for the probability of detecting PCa in all patients, achieving a C-index of 0.853 (26). In the current study, we chose the best prediction model of the base model + transCGA + coroPZA to construct a new nomogram that could provide the risk of PCa for individual Chinese patients. Internal validation showed a predictive accuracy (C-index = 0.913) for PCa, which gained an advantage over some previous nomograms developed by Chinese researchers.

Our study has certain limitations. As with any retrospective study, there was the risk of selection bias in assessing the value of the prostate maximum sectional area on mp-MRI. In addition to the issues surrounding the small sample size, our models were calibrated using internal validation only, with no external validation conducted to ensure their utility. Therefore, further clinical studies that employ long-term follow-up to evaluate our model’s practical applicability are required before it is prospectively applied to patients.

We first found that MRI-based CGA, PZA, and PA in the sagittal, transverse, and coronal section have potential predictive value for diagnosing PCa, especially transCGA. We put forward a new predictor named transPAI. We found that it was significantly higher in PCa patients. Our nomogram model based on age + MRI + PSA + transCGA + coroPZA had great predictive accuracy for PCa. The application of this nomogram model may further decrease the rate of unnecessary biopsies.

Data Availability Statement

Publicly available datasets were analyzed in this study. These data can be found here: https://www.wolai.com/nYw7M5iJDyxn91RDWTmoPj?theme=light.

Author Contributions

SJ conceived the study and carried out the investigation, methodology, and original draft preparation. ZH and BL carried out the investigation, methodology, and formal analysis, and participated in the conceptualization of the manuscript. ZC carried out the data curation, investigation, and methodology. YX, WZ, and YW participated in the formal analysis. ML conceived the study; carried out the supervision, review, and editing of the manuscript; and took responsibility for the integrity of the work as corresponding author. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Startup Fund for Scientific Research, Fujian Medical University (Grant number: 2018QH1044) and Joint Funds for the Innovation of Science and Technology, Fujian province (Grant number: 2017Y9023).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

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.

Acknowledgments

We thank Richard Lipkin, PhD, from Liwen Bianji, Edanz Group China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.

Supplementary Material

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

Supplementary Figure 1 | Prostate location and MR imagine in transverse (A), sagittal (B) and coronal (C) plane. Red: Central gland sectional area, CGA; Blue: Peripheral zone sectional area, PZA.

Supplementary Figure 2 | Correlation analysis among each predictor used to eliminate combinations including collinear predictors.

Supplementary Figure 3 | Nomogram predicting the probability of PCa at the initial biopsy based on model 2 (A) and model 3 (B).

Supplementary Figure 4 | Calibration plot in training cohort and validation cohort and predictive accuracy for PCa at initial biopsy based on model 2.

Supplementary Figure 5 | Calibration plot in training cohort and validation cohort and predictive accuracy for PCa at initial biopsy based on model 3.

Abbreviations

BMI, body mass index; FPSA, free PSA; PSA, prostate-specific antigen; FTPSA, $\frac{\text{free PSA}}{\text{PSA}};$ PV, prostate volume; sagiPSAPA, $\frac{\text{PSA}}{\text{sagi\hspace{0.17em}PA}};$ PSAD, PSA density; sagiPSACGA, $\frac{\text{PSA}}{\text{sagi\hspace{0.17em}CGA}};$ PA, prostate maximum sectional area; sagiPSAPZA, $\frac{\text{PSA}}{\text{sagi}\hspace{0.17em}\text{PZA}};$ CGA, central gland sectional area; transPSAPA, $\frac{\text{PSA}}{\text{trans\hspace{0.17em}PA}};$ PZA, peripheral zone sectional area; transPSACGA, $\frac{\text{PSA}}{\text{trans\hspace{0.17em}CGA}};$ sagiPA, sagittal prostate maximum sectional area; transPSAPZA, $\frac{\text{PSA}}{\text{trans\hspace{0.17em}PZA}};$ sagiCGA, sagittal central gland sectional area; coroPSAPA, $\frac{\text{PSA}}{\text{coro\hspace{0.17em}PA}};$ sagiPZA, sagittal peripheral zone sectional area; coroPSACGA, $\frac{\text{PSA}}{\text{coro\hspace{0.17em}CGA}};$ transPA, transverse prostate maximum sectional area; coroPSAPZA, $\frac{\text{PSA}}{\text{coro\hspace{0.17em}PZA}};$ transCGA, transverse central gland sectional area; sagiPAI, $\frac{\text{sagiPZA}}{\text{sagi}\hspace{0.17em}\text{CGA}};$ transPZA, transverse peripheral zone sectional area; transPAI, $\frac{\text{transPZA}}{\text{trans\hspace{0.17em}CGA}};$ coroPA, coronal prostate maximum sectional area; coroPAI, $\frac{\text{coroPZA}}{\text{coro\hspace{0.17em}CGA}};$ coroCGA, coronal central gland sectional area; ALP, alkaline phosphatase; coroPZA, coronal peripheral zone sectional area; LDH, lactate dehydrogenase

References

1. Henley SJ, Ward EM, Scott S, Ma J, Anderson RN, Firth AU, et al. Annual Report to the Nation on the Status of Cancer, Part I: National Cancer Statistics. Cancer (2020) 126(10):2225–49. doi: 10.1002/cncr.32802

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Cuzick J, Thorat M, Andriole G, Brawley O, Brown P, Culig Z, et al. Prevention and Early Detection of Prostate Cancer. Lancet Oncol (2014) 15(11):e484–92. doi: 10.1016/s1470-2045(14)70211-6

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Smith RA, Andrews KS, Brooks D, Fedewa SA, Manassaram-Baptiste D, Saslow D, et al. Cancer Screening in the United States, 2019: A Review of Current American Cancer Society Guidelines and Current Issues in Cancer Screening. CA: Cancer J Clin (2019) 69(3):184–210. doi: 10.3322/caac.21557

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Mottet N, Bellmunt J, Bolla M, Briers E, Cumberbatch M, De Santis M, et al. EAU-ESTRO-SIOG Guidelines on Prostate Cancer. Part 1: Screening, Diagnosis, and Local Treatment With Curative Intent. Eur Urol (2017) 71(4):618–29. doi: 10.1016/j.eururo.2016.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Shen P, Zhao J, Sun G, Chen N, Zhang X, Gui H, et al. The Roles of Prostate-Specific Antigen (PSA) Density, Prostate Volume, and Their Zone-Adjusted Derivatives in Predicting Prostate Cancer in Patients With PSA Less Than 20.0 Ng/mL. Andrology (2017) 5(3):548–55. doi: 10.1111/andr.12322

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Benson MC, Whang IS, Olsson CA, McMahon DJ, Cooner WH. The Use of Prostate Specific Antigen Density to Enhance the Predictive Value of Intermediate Levels of Serum Prostate Specific Antigen. J Urol (1992) 147:817–21. doi: 10.1016/s0022-5347(17)37394-9

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Bruno SM, Falagario UG, d'Altilia N, Recchia M, Mancini V, Selvaggio O, et al. PSA Density Help to Identify Patients With Elevated PSA Due to Prostate Cancer Rather Than Intraprostatic Inflammation: A Prospective Single Center Study. Front Oncol (2021) 11:693684. doi: 10.3389/fonc.2021.693684

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Zlotta A, Djavan B, Marberger M, Schulman C. Prostate Specific Antigen Density of the Transition Zone: A New Effective Parameter for Prostate Cancer Prediction. J Urol (1997) 157(4):1315–21. doi: 10.1097/00005392-199704000-00041

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Choi YJ, Kim JK, Kim HJ, Cho KS. Interobserver Variability of Transrectal Ultrasound for Prostate Volume Measurement According to Volume and Observer Experience. AJR Am J Roentgenology (2009) 192(2):444–9. doi: 10.2214/ajr.07.3617

CrossRef Full Text | Google Scholar

10. Yacoub JH, Oto A. MR Imaging of Prostate Zonal Anatomy. Radiol Clinics N Am (2018) 56(2):197–209. doi: 10.1016/j.rcl.2017.10.003

CrossRef Full Text | Google Scholar

11. Augustin H, Erbersdobler A, Hammerer PG, Graefen M, Huland H. Prostate Cancers in the Transition Zone: Part 2; Clinical Aspects. BJU Int (2004) 94(9):1226–9. doi: 10.1111/j.1464-410X.2004.05147.x

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Rahmouni A, Yang A, Tempany CM, Frenkel T, Epstein J, Walsh P, et al. Accuracy of in-Vivo Assessment of Prostatic Volume by MRI and Transrectal Ultrasonography. J Comput Assisted Tomography (1992) 16(6):935–40. doi: 10.1097/00004728-199211000-00020

CrossRef Full Text | Google Scholar

13. Balachandran V, Gonen M, Smith J, DeMatteo R. Nomograms in Oncology: More Than Meets the Eye. Lancet Oncol (2015) 16(4):e173–80. doi: 10.1016/s1470-2045(14)71116-7

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Iasonos A, Schrag D, Raj G, Panageas K. How to Build and Interpret a Nomogram for Cancer Prognosis. J Clin Oncol (2008) 26(8):1364–70. doi: 10.1200/jco.2007.12.9791

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Weinreb J, Barentsz J, Choyke P, Cornud F, Haider M, Macura K, et al. PI-RADS Prostate Imaging - Reporting and Data System: 2015, Version 2. Eur Urol (2016) 69(1):16–40. doi: 10.1016/j.eururo.2015.08.052

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Wu YP, Li XD, Ke ZB, Chen SH, Chen PZ, Wei Y, et al. Risk Factors for Infectious Complications Following Transrectal Ultrasound-Guided Prostate Biopsy. Infect Drug Resist (2018) 11:1491–7. doi: 10.2147/idr.S171162

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Tanaka N, Fujimoto K, Chihara Y, Torimoto M, Hirao Y, Konishi N, et al. Prostatic Volume and Volume-Adjusted Prostate-Specific Antigen as Predictive Parameters for Prostate Cancer Patients With Intermediate PSA Levels. Prostate Cancer Prostatic Dis (2007) 10(3):274–8. doi: 10.1038/sj.pcan.4500957

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Zheng S, Jiang S, Chen Z, Huang Z, Shi W, Liu B, et al. The Roles of MRI-Based Prostate Volume and Associated Zone-Adjusted Prostate-Specific Antigen Concentrations in Predicting Prostate Cancer and High-Risk Prostate Cancer. PLoS One (2019) 14(11):e0218645. doi: 10.1371/journal.pone.0218645

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Peng Y, Shen D, Liao S, Turkbey B, Rais-Bahrami S, Wood B, et al. MRI-Based Prostate Volume-Adjusted Prostate-Specific Antigen in the Diagnosis of Prostate Cancer. J Magn Reson Imaging JMRI (2015) 42(6):1733–9. doi: 10.1002/jmri.24944

CrossRef Full Text | Google Scholar

20. Tang P, Chen H, Uhlman M, Lin Y, Deng X, Wang B, et al. A Nomogram Based on Age, Prostate-Specific Antigen Level, Prostate Volume and Digital Rectal Examination for Predicting Risk of Prostate Cancer. Asian J Androl (2013) 15(1):129–33. doi: 10.1038/aja.2012.111

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Murray NP, Reyes E, Orellana N, Fuentealba C, Duenas R. A Comparative Performance Analysis of Total PSA, Percentage Free PSA, PSA Velocity, and PSA Density Versus the Detection of Primary Circulating Prostate Cells in Predicting Initial Prostate Biopsy Findings in Chilean Men. BioMed Res Int (2014) 2014:676572. doi: 10.1155/2014/676572

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Lorenzo G, Hughes TJR, Dominguez-Frojan P, Reali A, Gomez H. Computer Simulations Suggest That Prostate Enlargement Due to Benign Prostatic Hyperplasia Mechanically Impedes Prostate Cancer Growth. Proc Natl Acad Sci (2019) 116(4):1152–61. doi: 10.1073/pnas.1815735116

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zhang S, Zhao S, Fu X. Intensity Modulated Radiotherapy in Combination With Endocrinotherapy in the Treatment of Middle and Advanced Prostatic Cancer. Pakistan J Med Sci (2019) 35(5):1264–9. doi: 10.12669/pjms.35.5.591

CrossRef Full Text | Google Scholar

24. Shariat SF, Karakiewicz PI, Suardi N, Kattan MW. Comparison of Nomograms With Other Methods for Predicting Outcomes in Prostate Cancer: A Critical Analysis of the Literature. Clin Cancer Res (2008) 14(14):4400–7. doi: 10.1158/1078-0432.Ccr-07-4713

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Wang Y, Xie S, Shangguan X, Pan J, Zhu Y, Xin Z, et al. Prostate Transitional Zone Volume-Based Nomogram for Predicting Prostate Cancer and High Progression Prostate Cancer in a Real-World Population. J Cancer Res Clin Oncol (2017) 143(7):1157–66. doi: 10.1007/s00432-017-2389-3

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Huang Y, Cheng G, Liu B, Shao P, Qin C, Li J, et al. A Prostate Biopsy Strategy Based on a New Clinical Nomogram Reduces the Number of Biopsy Cores Required in High-Risk Patients. BMC Urol (2014) 14:8. doi: 10.1186/1471-2490-14-8

PubMed Abstract | CrossRef Full Text | Google Scholar