Three-dimensional measurement and analysis of benign prostatic hyperplasia
Introduction
Benign prostatic hyperplasia (BPH) is associated with lower urinary tract symptoms and bladder outlet obstruction. Many non-invasive parameters related to benign prostatic hyperplasia have been widely measured in two-dimensional (2D) planes. The parameters are characterized by non-invasive and convenient comparison, and can be used for BPH diagnosis (1). With the advancement of research, the prostatic urethral has gained much attention. During the hyperplasia process of prostate, the length of prostatic urethra and prostatic urethra angle will be altered (2). Although these changes have been measured in 2D planes (3-5), the degree of which the prostatic urethra is compressed by the prostate and the angle relationship between prostate urethral and bladder have not been investigated. Prostate volume (PV) has attracted considerable attention in 2D planes (6). However, it is difficult to measure intravesical prostatic volume (IPPV), intravesical prostatic thickness (IPPT) and intravesical prostatic protrusion surface area (IPPS).
Three-dimensional (3D) reconstruction can be used to measure the PV (7,8) and the changes of bladder wall thickness (BWT) (9). Its advantages are accurate and intuitive. However, other 2D parameters related to BPH have not been measured in 3D planes as well as other characteristics of prostate, prostate urethral, and bladder. In addition, little attention has been paid to the 3D reconstruction of prostate zonal anatomy.
We measured conventional parameters in 2D plane, as well as other features of BPH using three-dimensional measurement in this study. The purpose of this research is to better understanding the lower urinary tract anatomy of male and to provide references for clinical diagnosis and treatment of BPH. We present the following article in accordance with the STROBE reporting checklist (available at http://dx.doi.org/10.21037/tau-21-142).
Methods
Study population
This retrospective study analyzed clinical database of patients who visited Tianjin first central hospital for lower urinary tract symptoms/BPH from October 2016 to November 2019. A detailed history of 296 patients who underwent T2 MRI was obtained in our study. Exclusion criteria included: history of surgery of lower urinary tract, prostatitis and other non-benign prostatic hyperplasia diseases, patients did not carry a urinary catheter. Patients whose peripheral zone and/or central gland were impossible to identify were also excluded from this study. Finally, 98 patients who carried a urinary catheter with urinary retention which was caused by benign prostatic hyperplasia were eligible for analysis. They were divided into 4 groups by age: 50–60, 60–70, 70–80 and 80–90 (n50-60=14, n60-70=30, n60-70=36, n70-80=18). The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013) and approved by Institutional Review Board (IRB) of the Tianjin First Center Hospital (registration ID2018NO22KY). Individual consent for this retrospective analysis was waived.
3D reconstruction
Patients underwent T2-weighted MRI using a 3-Tesla magnet coil with a slice thickness of 3 mm. The first step, the collected MRI data from patients were imported into Mimics 19.0 software (Materialize, Leuven, Belgium) in DICOM file for 3D reconstruction of prostate, central gland, peripheral zone, intravesical prostatic protrusion, prostatic urethra and bladder based on horizontal, sagittal and coronal of pelvic T2-weighted MRI. To achieve accurate reconstruction of the lower urinary tract anatomy, this study applies corresponding 3D model-based registration for creation of fusion 3D model on account of each horizontal, sagittal and coronal plane. The Second step, the corresponding 3D models were fitted precisely in Geomagic Studio 14.0 software (Geomagic, Rock Hill, SC, USA) based on the same coordinate system. At the same time, the results of parameters related to 2D plane were measured. Finally, 3D measurement and analysis of related parameters were performed in 3-matic software (Materialize, Leuven, Belgium).
Parameter measurements
In 2D plane, the highest diameter of the prostate (HD) was measured on the maximum sagittal section. The longest diameter of the prostate (LD) and the widest diameter of the prostate (WD) were measured on the maximum horizontal section. Then, PV was calculated according to the ellipse volume formula: 0.52 × LD × WD × HD. The measurement of IPP was completed on the sagittal section. A straight line from the left lowest point to the right lowest point where the prostate protrudes into the bladder was drawn. The distance from the highest point of the intravesical prostatic protrusion to this line was defined as IPP.
In 3D plane, the following parameters were measured or analyzed by 3-matic software. IPPV, intravesical prostatic protrusion surface area (IPPS), IPPT, PV, prostate thickness (PT), prostate surface area (PS), central gland volume (CGV), central gland thickness (CGT), central gland surface area (CGS), peripheral zone volume (PZV), peripheral zone thickness (PZT), peripheral zone surface area (PZS), prostatic urethra thickness (PUT) and BWT were analyzed by 3-matic software. Central gland index (CGI) was determined by dividing CGV by PV. Peripheral zone index (PZI) was defined as the ratio PZV/PV. Intravesical prostatic protrusion index (IPPI) was defined as the ratio IPPV/PV. Based on the measurement of PV, they were divided into two groups: <80 mL (n=48) and >80 mL (n=50). LD, WD and HD were determined by 3-matic software from the highest and lowest points of the X-axis, Y-axis, and Z-axis in the world coordinate system, respectively.
The straight line which was used to measure IPP on the 2D plane previously was imported into 3-matic software. The plane in which the intravesical prostatic protrusion lies was established by 3-matic software was based on this straight line. The highest point of intravesical prostatic protrusion was taken as the endmost of intravesical prostatic protrusion tangent to a plane which is parallel to the previously established plane. The distance from the highest point to the plane was defined as IPP. The center of gravity of prostate was analyzed by 3-matic software, and the coronal plane and sagittal plane were built by 3-matic software based on the center of gravity of prostate. The central line of proximal prostatic urethra and distal prostatic urethra were analyzed by 3-matic software. The angle between the two central lines was defined as prostatic urethra angle (PUA). The angle of central line of the proximal prostatic urethra with respect to the coronal plane and sagittal plane were defined as the A-angle and B-angle, respectively. The angle of the central line of distal prostatic urethra with respect to the coronal plane and sagittal plane were defined as the C-angle and D-angle, respectively. The angle between the bladder and prostatic urethra was defined as the E-angle. The highest point and the lowest point of prostatic urethra were analyzed by 3-matic software. Then, the distance from the highest point to the coronal plane and sagittal plane were defined as the A-Distance and B-Distance, respectively. The distance from the lowest point to the coronal plane and sagittal plane were defined as the C-Distance and D-Distance, respectively. The surface distance from proximal endpoint of prostatic urethra to distal endpoint of prostatic urethra was defined as prostatic urethra length (PUL).
Statistical analysis
Values were shown as mean ± standard deviation. Two independent-samples t-tests or Mann-Whitney U test were used to compare the parameters of PV <80 mL group and PV >80 mL group as well as 2-D measurements and 3D measurements. All parameters of proximal prostatic urethra and distal prostatic urethra were compared by the paired t-test or Wilcoxon Rank-Sum test. Comparisons between age groups were made by one-way ANOVA or Kruskall-Wallis test. All statistical methods were chosen according to whether the data conformed to normal distribution. All statistical analysis was analyzed by SPSS software 25.0 (IBM, Armonk, NY, USA). P<0.05 considered statistically significant.
Results
According to the T2-weighted MRI of patients, 3D models of prostate, intravesical prostatic protrusion, central gland, peripheral zone, bladder, and prostatic urethra were successfully reconstructed. These models can be observed from different points of view (Figure 1). Moreover, many parameters of the modes were able to be measured and analyzed by Mimics software, such as angle, diameter line, distance and thickness (Figures 2,3). Table 1 showed the mean measurements of all BPH parameters. Age, PV, IPPV, CGV, PZV, IPPI, CGI, LD, WD, HD, IPP, PUL, A-Distance, A-Angle, PT, IPPT, CGT, PS, IPPS, CGS, PZS of PV >80 mL group were significantly higher than these in PV <80 mL group (P<0.001, 0.001, 0.001, 0.001, 0.01, 0.05, 0.01, 0.001, 0.001, 0.001, 0.01, 0.001, 0.01, 0.05, 0.001, 0.01, 0.001, 0.001, 0.001, 0.001, 0.001, respectively). PZI of PV<80 mL group was significantly higher than it in the PV >80 mL group (P<0.05). No significant difference was noted in B-Distance, C-Distance, D-Distance, PUA, B-Angle, C-Angle, D-Angle, E-Angle, PZT, BWT, PUT between the two groups (Table 2).
Full table
Full table
The WD was significantly higher in 3D measurement than it in 2D measurement (mean 57.79±8.85 vs. 50.42±9.06, P<0.001). Meanwhile, LD (P<0.05). IPP, PV, LD and HD in 2D measurement and 3D measurement were not significantly different (Table 3). A-Distance and A-Angle were significant higher than B-Distance and B-Angle, respectively (both P<0.001). C-Distance and C-Angle were significantly higher than D-Distance and D-Angle, respectively (both P<0.001). A-Distance, B-Distance, A-Angle and B-Angle were significant higher than C-Distance, D-Distance, C-Angle and D-Angle, respectively (all P<0.001) (Table 4).
Full table
Full table
Moreover, PV, IPPV, CGV, PZV, CGI, LD, WD, HD, IPP, PUL, A-Distance, PT, IPPT, CGT, PS, IPPS, CGS and PZS remarkably increased with age (P<0.001, 0.001, 0.001, 0.05, 0.01, 0.001, 0.001, 0.001, 0.01, 0.001, 0.001, 0.001, 0.01, 0.001, 0.001, 0.001, 0.001, 0.001, respectively), while PZI significantly decreased with age (P<0.01). However, there were no significant differences in IPPI, B-Distance, C-Distance, D-Distance, PUA, A-Angle, B-Angle, C-Angle, D-Angle, E-Angle, PZT, BWT and PUT in different age categories (Table 5).
Full table
Discussion
MRI was a robust and reliable method for measuring PV (10). MRI with 3D segmentation software could provide an effective solution for accurate measurements of PV (11). In this study, 3D measurement was used not only to measure various 2D parameters related to BPH but also other characteristics of BPH. This study showed that 3D measurement can be used as a mean to understand the characteristics and pathology progression in BPH.
This study found that 3D plane can provide us with much more accurate measurements than general 2D plane. There was statistical difference in the WD between the 2D measurement and 3D measurement. The WD had important implications in prostatectomy. Based on WD measurement, the anatomic structure of ejaculatory ducts can be determined, which had important implications for retention of ejaculatory function following transurethral resection of prostate (12). Measurement of diameter line in 2D planes is restricted by irregularities in prostate contour and the result is subjectively influenced by the operator. There may be a lack of anatomical features in the three-dimensional model reconstructed by using the horizontal plane alone. To compensate for this defect and achieve accurate reconstruction of the lower urinary tract anatomy, this study applies 3D model-based registration for creation of fusion 3D model on account of each horizontal, sagittal and coronal plane. Moreover, measurement of the diameter line was analyzed by 3-matic software. Therefore, 3D measurement guarantees more reliable and robust results.
Though formula of PV had been commonly explored, measurement errors were inevitable (13,14). Measurement of PV was successfully achieved in this study. Measurement of the PV can serve as a reference standard to use of alpha blockers (15), perform prostatectomy (16), and select prostate biopsy site (17). The standard critical value of PV =80 mL was chosen for the prostatectomy (18). The IPP was the important parameter to assess BPH (19). IPP, IPPT, IPPV, IPPI and IPPS of PV >80 mL group were significantly higher than those in PV <80 mL group in our study. IPP, IPPT and IPPS significantly increased with age. IPPV and PV also strikingly increased with age, while IPPI did not change with age. IPPI was defined as a ratio between IPPV and PV. This indicated that increment of volume was almost the same in the PV and IPPV.
Tian et al. and Kwon et al. had used transrectal ultrasound (TRUS) to measure the PZT. Tian et al. discovered that the mean value of PZT was 13.3±6.3 mm, which was not modulated by age (20). However, Kwon et al. found that the mean value of PZT was 11.10±2.50 mm, which was significantly decreased with age (21). Those results were different from the result of our present study. The mean value of PZT in this study was 13.75±3.77 mm. The differences may result from the inaccuracy of results measured by TRUS. In TRUS, PZT measurement was required to determine the longest line which was draw from the center of the transitional zone to the outer margin of the peripheral zone. Apart from that, peripheral zone was relatively difficult to be differentiated in TRUS. The profile of the peripheral zone was clearer in T2-weighted MRI. Accordingly, peripheral zone was reconstructed precisely. Moreover, PZT was analyzed by 3-matic software. The analysis tools employed to measure PZT was reliable, which guaranteed the quality of this data. PZT had no significant difference in the patients above the age of 50. Our result was consistent with Tian et al.
The shape of peripheral zone is irregular. It was too difficult to measure the PZV and to find a solution to the measurement of PZV. However, 3D measurement was a good way to solve this problem. Currently, a few researches about PZV were controversial. Matsugasumi et al. reported findings from a retrospective analysis that the PZV maintained steadily at about 20 mL in patients with BPH, and the indirect influence of the age on the variation in PZV was small (22). However, findings from a longitudinal study by Williams et al. found that the PZV in the 50–60 years age group was significantly higher than it in <45 years age group, which increased about 5 mL (23). These differences were probably due to variations in study methodology.
Similar to study by Williams et al., we found that with the increased in PV, PZV was increasing as well, but the magnitude of this change is relatively small. Comparing to 50–59 years age group and 60–69 years age group, the PZV for 60–69 years age group and 70–79 years age group increased few. In addition, this study found that PZI significantly decreased with age. This indicated that greater increases in the PV than in the PZV. Subjects in study by Matsugasumi et al. were Japanese and Americans. While in the Williams’s study and our research, the subjects were Americans and Chinese respectively. Therefore, the discrepancy might be caused by the different races.
Comparing to PV <80 mL group, the PZV of PV >80 mL group increased also few. However, the PZI of PV <80 mL group was significantly higher than the PV >80 mL group. PZT showed no significant difference between PV <80 mL group and PV >80 mL group. CGT and PT sharply increased with PV as well as age. This could be the reason why the correlation between BPH and transition zone was significantly higher than peripheral zone (24). PS, CGV, CGS, PZS, LD, WD and HD significantly increased with age as well as PV in our study. During the process of BPH, the prostate got closer to a circle. These results were consistent with the PCAR theory (25).
In 2D planes, PUA had a different definition. PUA was defined as the angle between line from seminal colliculus to bladder neck and rostocaudal axis (4). The vertex angle formed by two rays including both the proximal and distal prostatic urethra was also defined as PUA (26). In contrast to those, considering the characteristics of 3D plane, we used the same criteria to define PUA alteration as study by Ko et al. (2). The PV was increased, leading to prolongation of PUL. The PUL significantly increased with age. However, the PUA decreased with PV although the range was slight. The result was consistent with study by Chen et al. (27). A study by El-Tatawy et al. revealed that PUA did not show a significant difference in age (28), which was concordant with finding of our study. There were no significant differences in PUA between different age groups. With reference to coronal plane and sagittal plane, the angle and distance of which the proximal prostatic urethra was compressed by the prostate were higher than distal prostatic urethra in this study. Therefore, the degrees of the compression of proximal prostatic urethra were larger than distal prostatic urethra. With reference to coronal plane, the angle and distance of the proximal prostatic urethra and distal prostatic urethra were significantly higher than those with reference to sagittal plane. A-Distance and A-Angle of PV >80 mL group were significantly higher than those in PV <80 mL group. A-Distance significantly increased with age.
This indicates that the rearward slope of the prostatic urethra is greater than the left or right offset during the process of BPH. To date, only one study has measured PUT. The PUT was measured as 6 [2–16] mm in a study by Guneyli et al. (29). However, the median (IQR) value of PUT was measured as 11.63 (6.91–17.41) mm in this study. The discrepancy might be due to patients had or had not a urinary catheter. Subjects in our study had a urinary catheter. The fine structure of prostatic urethra will be clearer and the reconstruction will be more accurate if patients had a urinary catheter. The PUT was analyzed by 3-matic software, which avoided the influence of subjective factors. No significant differences were found in PUT between any of the groups. This suggests BPH has little impacts on the PUT. No differences in BWT and E-Angle were also found between any of the groups. BWT was influenced by the volume of urine (30). Besides, patients who were unable to void may have some effect on BWT.
However, this study also has some limitations. First, while we were able to measure multiple features of BPH, our choice of regions of interest was based on anatomical contour, which was time-consuming. Future work is needed to explore methods of automated 3D segmentation. In addition, the BWT and E-Angle results may be affected by the fact that the patients had a urinary catheter.
Conclusions
To sum up, 3D measurement is more accurate than 2D measurement. PZV was increasing, but the magnitude of this change is relatively small and the rearward slope of the prostatic urethra is greater than the left or right offset during the process of BPH. The BPH has a minor impact on the PZT and PUT, and 3D measurement is a great tool for measuring features of BPH, which provides good supports for further clinical and scientific research.
Acknowledgments
Funding: This work was supported by grants from the National Natural Science Foundation of China (31871212, Lidong Zhai), Postdoctoral Science Foundation of China (2016M601272, Lidong Zhai), Chinese Clinical Trial Registry (ChiCTR1800014818, Jianhui Wu).
Footnote
Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at http://dx.doi.org/10.21037/tau-21-142
Peer Review File: Available at http://dx.doi.org/10.21037/tau-21-142
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at http://dx.doi.org/10.21037/tau-21-142). YX serves as an unpaid editorial board member of Translational Andrology and Urology. The other authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The study was approved by Institutional Review Board (IRB) of the Tianjin First Center Hospital (registration ID2018NO22KY) and individual consent for this retrospective analysis was waived.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Malde S, Nambiar AK, Umbach R, et al. Systematic Review of the Performance of Noninvasive Tests in Diagnosing Bladder Outlet Obstruction in Men with Lower Urinary Tract Symptoms. Eur Urol 2017;71:391-402. [Crossref] [PubMed]
- Ko YH, Kim TH, Song PH, et al. Structural Variations of the Prostatic Urethra Within the Prostate Predict the Severities of Obstructive Symptoms: A Prospective Multicenter Observational Study. Urology 2017;104:160-5. [Crossref] [PubMed]
- Ko YH, Song PH. Structural Variation of Prostate Urethra Reflected by the Ratio Between Prostate Volume and Prostatic Urethral Length is Associated with the Degrees of Lower Urinary Tract Symptoms. Low Urin Tract Symptoms 2016;8:113-9. [Crossref] [PubMed]
- Minagawa T, Daimon H, Ogawa N, et al. Morphological and clinical evaluation of prostatic urethra using modified sonourethrography with retrograde jelly injection. Low Urin Tract Symptoms 2019;11:O4-10. [Crossref] [PubMed]
- Agrawal V, Khullar R, Jha AK. Assessment of posterior urethra in benign prostatic hyperplasia and after its surgery. Urol Ann 2020;12:63-8. [Crossref] [PubMed]
- David RA, Badmus TA, Salako AA, et al. Diagnostic performance of transrectal ultrasound for prostate volume estimation in men with benign prostate hyperplasia. Int J Clin Pract 2020;74:e13615 [Crossref] [PubMed]
- Chen L, Feng C, Lv XG, et al. Three-Dimensional Computerized Model Based on the Sonourethrogram: A Novel Technique to Evaluate Anterior Urethral Stricture. J Urol 2018;199:568-75. [Crossref] [PubMed]
- Ni J, Bongers A, Chamoli U, et al. In Vivo 3D MRI Measurement of Tumour Volume in an Orthotopic Mouse Model of Prostate Cancer. Cancer Control 2019;26:1073274819846590 [Crossref] [PubMed]
- Anzia LE, Johnson CJ, Mao L, et al. Comprehensive non-invasive analysis of lower urinary tract anatomy using MRI. Abdom Radiol (NY) 2021;46:1670-6. [Crossref] [PubMed]
- Martins T, Mussi TC. Prostate volume measurement by multiparametric magnetic resonance and transrectal ultrasound: comparison with surgical specimen weight. Einstein 2020;18:eAO4662 [Crossref] [PubMed]
- Bezinque A, Moriarity A, Farrell C, et al. Determination of Prostate Volume: A Comparison of Contemporary Methods. Acad Radiol 2018;25:1582-7. [Crossref] [PubMed]
- Malalasekera AP, Sivasuganthan K, Sarangan S, et al. Morphological variations of the human ejaculatory ducts in relation to the prostatic urethra. Clin Anat 2018;31:456-61. [Crossref] [PubMed]
- Stanzione A, Ponsiglione A, Di Fiore GA, et al. Prostate Volume Estimation on MRI: Accuracy and Effects of Ellipsoid and Bullet-Shaped Measurements on PSA Density. Acad Radiol 2020;S1076-6332(20)30290-7.
- Aprikian S, Luz M, Brimo F, et al. Improving ultrasound-based prostate volume estimation. BMC Urol 2019;19:68. [Crossref] [PubMed]
- El-Adawy MS, Abdelaziz AY, Salem A, et al. Relation of baseline prostate volume to improvement of lower urinary tract symptoms due to tamsulosin monotherapy in benign prostatic hyperplasia: An exploratory, multicenter, prospective study. Urol Ann 2020;12:271-5. [Crossref] [PubMed]
- Wang B, Zhang S, Sun C, et al. Comparison between a transurethral prostate split and transurethral prostate resection for benign prostatic hyperplasia treatment in a small prostate volume: a prospective controlled study. Ann Transl Med 2020;8:1016. [Crossref] [PubMed]
- Stone NN, Crawford ED, Skouteris VM, et al. The Ratio of the Number of Biopsy Specimens to Prostate Volume (Biopsy Density) Greater Than 1.5 Improves the Prostate Cancer Detection Rate in Men Undergoing Transperineal Biopsy of the Prostate. J Urol 2019;202:264-71. [Crossref] [PubMed]
- Gunseren KO, Akdemir S, Çiçek MC, et al. Holmium Laser Enucleation, Laparoscopic Simple Prostatectomy, or Open Prostatectomy: The Role of the Prostate Volume in terms of Operation Time. Urol Int 2021;105:285-90. [PubMed]
- Topazio L, Perugia C, De Nunzio C, et al. Intravescical prostatic protrusion is a predictor of alpha blockers response: results from an observational study. BMC Urol 2018;18:6. [Crossref] [PubMed]
- Tian Y, Liu HM, Yang B, et al. Prostatic anatomical parameters correlate with clinical characteristics suggestive of benign prostatic hyperplasia. Asian J Androl 2021;23:64-8. [Crossref] [PubMed]
- Kwon JK, Han JH, Choi HC, et al. Clinical significance of peripheral zone thickness in men with lower urinary tract symptoms/benign prostatic hyperplasia. BJU Int 2016;117:316-22. [Crossref] [PubMed]
- Matsugasumi T, Fujihara A, Ushijima S, et al. Morphometric analysis of prostate zonal anatomy using magnetic resonance imaging: impact on age-related changes in patients in Japan and the USA. BJU Int 2017;120:497-504. [Crossref] [PubMed]
- Williams AM, Simon I, Landis PK, et al. Prostatic growth rate determined from MRI data: age-related longitudinal changes. J Androl 1999;20:474-80. [PubMed]
- Shi YF, Yu DJ. TRAF6 regulates proliferation of stromal cells in the transition and peripheral zones of benign prostatic hyperplasia via Akt/mTOR signaling. Prostate 2018;78:193-201. [Crossref] [PubMed]
- Watanabe H. New concept of BPH: PCAR theory. Prostate 1998;37:116-25. [Crossref] [PubMed]
- Shim M, Bang WJ, Oh CY. Correlation between prostatic urethral angulation and symptomatic improvement after surgery in patients with lower urinary tract symptoms according to prostate size. World J Urol 2020;38:1997-2003. [Crossref] [PubMed]
- Chen H, Zhang CA, Gill HS. Anatomical study of the prostatic urethra using vinyl polysiloxane casts. Prostate 2020;80:241-6. [Crossref] [PubMed]
- El-Tatawy H, Gameel T, El-enen MA, et al. Clinical significance of prostatic-urethral angulation on the treatment outcome of patients with symptomatic benign prostatic hyperplasia treated with tamsulosin hydrochloride. Arch Ital Urol Androl 2015;87:238-42. [Crossref] [PubMed]
- Guneyli S, Ward E, Peng Y, et al. MRI evaluation of benign prostatic hyperplasia: Correlation with international prostate symptom score. J Magn Reson Imaging 2017;45:917-25. [Crossref] [PubMed]
- Fry CH, McCloskey KD. Spontaneous Activity and the Urinary Bladder. Adv Exp Med Biol 2019;1124:121-47. [Crossref] [PubMed]