PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Cancer Epidemiol Biomarkers Prev. Author manuscript; available in PMC 2010 January 1.
Published in final edited form as:
PMCID: PMC2651680
NIHMSID: NIHMS88000

The Association of Bone Mineral Density with Prostate Cancer Risk in the Osteoporotic Fractures in Men (MrOS) Study

Abstract

We investigated the association of bone mineral density (BMD) measures with prostate cancer (PCa) risk in older men enrolled in the Osteoporotic Fractures in Men (MrOS) Study. We hypothesized that men with higher BMD, a marker of exposure to endogenous sex hormones, would have an increased incidence of PCa.

The cohort included 4597 men (89% white, 65 years or older) with no prior history of PCa. Baseline total body, total hip, and spine BMD were assessed using dual energy x-ray absorptiometry. Prostate cancer was confirmed by review of medical records. Cox regression was used to assess the association of BMD quartiles with incident PCa, adjusting for age, BMI, and other covariates.

During an average follow-up of 5.2 years, 5.6% (N=255) of men developed PCa. Total body BMD was inversely associated with incident PCa, with a significant trend for decreasing PCa risk with increasing BMD quartiles (p-trend=0.007). Men in the highest total body BMD quartile had a 41% reduced risk for prostate cancer (HR=0.59, 95% CI: 0.40–0.86), compared to men in the lowest quartile. Total hip and spine BMD did not exhibit significant relationships with PCa. Associations of BMD measures differed for low-grade (Gleason sum 2–6) vs. high-grade tumors (Gleason sum ≥7). Significant inverse relationships with high-grade disease were noted at the total body and total hip sites. However, no associations were observed with low-grade disease.

Our results provide support for an inverse association between BMD and prostate cancer risk. Possible pathophyisological mechanisms linking BMD and prostate cancer should be elucidated.

Keywords: Bone mineral density, prostate cancer, MrOS Study

BACKGROUND

Prostate cancer (PCa) is the most common malignancy in men. According to the American Cancer Society, it is estimated to result in 218,890 new cases and 27,050 deaths in 2007 (1). Despite its public health burden, the etiology of PCa remains poorly understood. Apart from the established risk factors for PCa including older age, race, and family history of the disease, only a few other factors have been suggested to contribute to the development of this malignancy.

Androgens are thought to stimulate cell proliferation of the prostate epithelium and have been related to increased risk of PCa in some but not most prospective epidemiologic studies (2, 3). Higher levels of insulin-like growth factor 1 (IGF-1) also have been implicated in prostate carcinogenesis (4). Additionally, high calcium intake has been associated with elevated PCa risk (5). However, evidence for the role of these factors is still inconclusive as conflicting results were reported by other studies (6, 7). A common limitation of these studies was the use of a single assessment of androgens and growth factor levels. Owing to the intra-individual variations in the levels of these factors, a single measurement at one time point may not accurately reflect average or cumulative exposure.

Bone mineral density (BMD) has been regarded as a possible surrogate marker for lifetime exposure to endogenous sex hormones, IGF-1, and calcium intake (8, 9). Based on this, the association of BMD with prostate cancer was examined in a few epidemiologic studies (1012). Results from the Tobago Prostate Survey indicated a cross-sectional association between higher BMD and increased prostate cancer prevalence in Afro-Caribbean men aged 60–79 years (10). Similarly, in the Framingham Study, a trend for higher prostate cancer risk with increasing BMD quartiles was noted in Caucasian men (12). However, conflicting findings were reported in the NHANES I Epidemiologic Follow-up Survey where a non-significant decline in prostate cancer risk was observed with higher BMD (11).

The majority of these studies have not utilized state of the art assessments of bone density (11, 12). While the Tobago Prostate Survey used dual-energy X-ray absorptiometry (DXA), it was cross-sectional in nature, and therefore did not account for the possible confounding effect of PCa on bone density (10). To our knowledge, no study has prospectively assessed the association of DXA-determined BMD measures with incident prostate cancer in men with no history of the disease. Additionally, it is not known whether this association varies by tumor grade.

The aim of the present study was to investigate the association of BMD measures with the subsequent development of prostate cancer in older men participating in the Osteoporotic Fractures in Men (MrOS) study. Our hypothesis was that men with higher BMD would have an increased risk of prostate cancer.

METHODS

Study Population

Participants were enrolled in the MrOS study, a multi-center longitudinal study evaluating risk factors and sequelae of vertebral and non-vertebral fractures in older men. The cohort included 5995 community dwelling, ambulatory men aged 65 years or older. Men were recruited from March 2000 through April 2002 at six geographic regions of the United States: Birmingham, Alabama; Minneapolis, Minnesota; Palo Alto, California; Pittsburgh, Pennsylvania; Portland, Oregon; and San Diego, California. Men were not eligible for inclusion if they: (1) were unable to walk without the assistance of another person or aide, (2) had bilateral hip replacements, (3) were unable to provide self-reported data, (4) were not expected to reside near a clinical site for the duration of the study, (5) had a medical condition that (in the judgment of the investigator) would result in imminent death, or (6) were unable to understand and sign an informed consent. To qualify as an enrollee, the participant had to answer the self-administered questionnaire, attend the clinic visit, and complete at least the anthropometric, DXA, and vertebral X-ray procedures. The Institutional Review Board at each recruitment site approved the study protocol, and written informed consent was obtained from all participants (13, 14).

The current analysis included 4597 participants. We excluded men with a self-reported history of prostate cancer (N= 709) or any other cancer (except for non-melanoma skin cancer) (N= 384), and men who have received osteoporosis medications (N= 90), androgen (N= 22) or antiandrogen (N= 150) therapy, or testosterone injections (N= 55).

Prostate Cancer Diagnosis

Prostate cancer cases occurring between the baseline visit and February 2007 (range of follow-up= 0.0–6.8 years, mean= 5.2 years, median= 5.4 years) were identified through self-report using a mailed tri-annual follow-up questionnaire. For participants who did not return the questionnaire, information about events was elicited through in-person or telephone interviews. For each reported event, medical records were requested from the hospital or clinic including: pathology reports for initial diagnosis of prostate cancer, PSA lab reports prior to diagnosis, clinical notes ordering biopsy, post-diagnosis studies reports, and post-diagnosis clinic notes. Medical records were reviewed and events were adjudicated centrally at the MrOS coordinating center (University of California, San Francisco) and California Pacific Medical Center Research Institute without knowledge of BMD or other risk factors. Key prognostic characteristics of the tumor (stage and Gleason histologic scores) and type of treatment were also collected.

BMD Measurement

Areal BMD (g/cm2) measures of the total body, total hip, and spine were measured at the baseline visit using DXA (QDR 4500 W scanner, Hologic Inc., Waltham, MA, USA). Scans were performed at each study center by certified technicians using a standardized protocol for participant positioning and scan analysis. For quality assurance, the MrOS coordinating center reviewed a random subset of scans, scans with exceptionally high or low BMD, and problematic scans identified by technicians at the clinics. Cross-calibration studies were performed prior to the baseline MrOS visit. No linear differences across scanners were observed and the inter-scanner coefficient of variation was 0.9% for the hip and 0.6% for the spine.

Covariates

At the baseline visit, participants completed a self-administered questionnaire and were interviewed and examined by trained and certified clinical staff. Demographic characteristics included age, race/ethnicity (Caucasian/White, African American/Black, Asian, Hispanic and Other) and educational level. Lifestyle risk factors included alcohol consumption (current drinking vs. not), smoking (current, past, never), and physical activity as reported on the Physical Activity Scale for the Elderly (PASE) (15). Personal history of specific medical conditions (e.g., cancer, diabetes mellitus, osteoporosis, and prostatitis) was assessed. Participants were asked to bring in current prescription medications to the clinic visit. Specific classes of medications (e.g. thiazide diuretics, oral corticosteroids, statins, osteoporosis medications, etc) were coded by trained staff using a computerized database. The intakes of dietary calcium and vitamin D from foods and supplements were estimated using a modified Block semiquantitative food frequency questionnaire developed specifically for MrOS by Block Dietary Systems (Berkeley, CA) (16). Total intakes of calcium and vitamin D were calculated by summing dietary calcium intake (milligrams per day) and daily dosage of calcium supplements (milligrams per day). Anthropometric measures such as weight and height were measured using standard equipment, including a Harpenden stadiometer and a balanced beam or digital scale. Body mass index (BMI) was calculated as weight in kilograms divided by the square of height in meters.

Data Analysis

Baseline characteristics and BMD measures of men with or without incident PCa were compared using chi-square test for categorical variables and either 2-sample t-test or Wilcoxon rank-sum test for continuous data. The incidence rate of prostate cancer was calculated as the number of cases divided by the person-years of follow-up. Follow-up time was calculated from the date of study entry to the date of prostate cancer diagnosis, death, or last contact with the participant. Cox proportional hazards regression was used to estimate the hazard ratio (HR) of incident PCa by BMD quartiles. Unadjusted and adjusted models were fitted for each BMD variable separately. Variables were selected for entry into the regression models if they were significantly associated with prostate cancer in univariate analysis. Race and body mass index were adjusted for regardless of their statistical significance. Linear trend in the risk of PCa across BMD quartiles was tested by including the median values for BMD quartiles as a single continuous variable. Additional analyses were performed to model the risk of PCa per one SD decrease in BMD (calculated as the deviation from the mean BMD divided by the standard deviation of the BMD measure). The proportional hazards assumption was checked by testing the significance of interaction terms of BMD variables with time. In secondary analyses, separate models were performed to estimate the risk of high-grade (defined as Gleason sum ≥ 7) vs. low-grade (defined as Gleason sum 2–6) tumors by BMD quartiles. The level of significance was set at 0.05. Data was analyzed using SAS version 8.01 (SAS Institute Inc, Cary, NC, USA).

RESULTS

During an average follow-up of 5.2 years (range of follow-up= 0.0–6.8 years, median= 5.4 years), 5.5% (255 out of 4597) of the men were diagnosed with prostate cancer. The median age at diagnosis was 74.7 years. Fifty one percent of the cases had a high histologic grade (defined as Gleason sum ≥7). Based on their baseline characteristics, men with incident prostate cancer were younger, had a higher level of physical activity and daily calcium intake, and were more likely to have a family history of prostate cancer, compared to men without prostate cancer (Table 1).

Table 1
Baseline characteristics and bone mineral density measures (mean ± SD, N (%)) by prostate cancer status in the MrOS cohort

For total body BMD, the incidence rates of prostate cancer (per 1000 person-years) were found to decrease by increasing BMD quartiles (12.3 in the lowest quartile, 11.7 in the second, 11.4 in the third, and 7.7 in the highest quartile). A significant trend for decreasing prostate cancer risk with increasing BMD quartiles was observed in unadjusted (p-value for trend= 0.02), age-adjusted (p-value for trend= 0.01), and multivariate models (p-value for trend= 0.007). Men in the highest total body BMD quartile had a 41% reduced risk for prostate cancer (multivariate-adjusted HR= 0.59, 95% CI: 0.40–0.86), compared to men in the lowest quartile. However, men in the second and third quartiles of BMD did not have a significant reduction in PCa risk as compared to men in the lowest quartile. When total body BMD was entered in the model as a continuous variable, each SD increase in BMD was associated with a 14% reduced risk of prostate cancer (multivariate-adjusted HR= 0.86, 95% CI: 0.75–0.98) (Table 2).

Table 2
Risk of prostate cancer by baseline total body BMD quartiles in the MrOS cohort

Total hip BMD was not associated with the incidence of prostate cancer in unadjusted or adjusted analyses. Compared to the lowest quartile of total hip BMD, the adjusted relative risk of prostate cancer was 0.97 (95% CI: 0.98–1.38) in the second quartile, 1.03 (95% CI: 0.72–1.46) in the third quartile, and 0.91 (95% CI: 0.62–1.33) in the highest quartile (Table 3).

Table 3
Risk of prostate cancer by baseline total hip BMD quartiles in the MrOS cohort

Similarly, in unadjusted and adjusted Cox regression models, spine BMD was not significantly associated with the risk of prostate cancer (Table 4).

Table 4
Risk of prostate cancer by baseline spine BMD quartiles in the MrOS cohort

Associations of BMD measures differed for low-grade vs. high-grade tumors. A significant inverse relationship was observed for total body BMD with high-grade PCa. A significant trend for lower high-grade PCa risk was observed with increasing BMD quartiles (p for trend= 0.01), with men in the highest quartile having a 57% reduced risk, compared to men in the lowest quartile (HR= 0.43, 95% CI 0.25–0.74). For total hip BMD, men in the second quartile had a 46% reduced risk of high-grade disease, as compared to men in the first quartile (HR= 0.54, 95% CI 0.32–0.92). On the other hand, none of the BMD measures was significantly associated with the risk of developing low-grade PCa (Table 5).

Table 5
Risk of high grade and low grade prostate cancer by baseline total body, total hip, and spine BMD quartiles in the MrOS cohort

DISCUSSION

This prospective analysis evaluated the association of BMD measures with incident prostate cancer in a cohort of older men with no history of PCa. Unexpectedly, we found that higher BMD of the total body was significantly related to reduced risk for prostate cancer. No associations were observed for hip and spine BMD measures with PCa. Additionally, total body BMD was inversely associated with the development of high-grade, but not low-grade disease. A similar but weaker association was observed for total hip BMD with high-grade PCa.

The direction of the association of total body BMD with prostate cancer was contrary to our initial hypothesis. This finding lends support to results from the NHANES I Epidemiologic Follow-up Survey where a decline in prostate cancer risk, although not significant, was observed with higher quartiles of bone density, determined using radiographic absorptiometry. In that study, compared to the lowest BMD quartile, the age, race, and BMI-adjusted rate ratios across BMD quartiles were 0.63 (95% CI: 0.37–1.07) for the second, 0.86 (95% CI: 0.49–1.49) for the third, and 0.72 (95% CI: 0.38–1.38) for the highest quartile (11).

Our results were discordant with those from the Tobago Prostate Survey and the Framingham Study, where higher bone density was associated with an increased prostate cancer risk (10, 12). The Tobago study evaluated Afro-Caribbean men, and reported a significant trend of higher prostate cancer prevalence with increasing total hip BMD quartiles, in men aged 60–79 years. Compared to men in the lowest quartile of BMD, those in the highest quartile had a two-fold higher odds of prostate cancer (OR= 2.12, 95% CI: 1.21–3.71). No such associations were observed in younger men in this cohort (10). Notably, total hip BMD in this population was higher than that observed for African-American men in NHANES III by approximately one standard deviation (17). This may be reflective of a higher exposure to endogenous sex steroid hormones and growth factors. The Framingham study involved 100 cases of prostate cancer, diagnosed at a median age of 75.2 years. It reported a higher risk of PCa in the upper two quartiles of radiogrammetrically-determined metacarpal cortical width, assessed at a mean age of 61 years (12). While these findings were adjusted for important risk factors such as age and BMI, they may be prone to residual confounding by factors such as family history of prostate cancer and calcium and vitamin D intakes, which were not collected in the above studies.

Sex steroids are important regulators of skeletal growth and maintenance of BMD in both men and women (18, 19). In the Swedish arm of the MrOS cohort, free testosterone was found to be a positive predictor of bone density at the total body, total hip, femur trochanter, and arm (20). Androgens are also thought to stimulate cell proliferation of the prostate epithelium and were related to increased risk of PCa in some prospective epidemiologic studies (2, 3). While the direction of association observed in our study does not provide a direct evidence for the role of endogenous sex hormones in the link between BMD and prostate cancer, it suggests a possible involvement of other pathophysiological mechanisms.

Poor vitamin D status may be a common denominator for the inverse association between BMD and prostate cancer risk. Low levels of vitamin D are known to have detrimental effects on bone density (21) and have been implicated in prostate carcinogenesis. In-vivo evidence suggests that calcitriol has anti- proliferative and chemopreventive effects in PCa (22). Additionally, in some epidemiologic studies, low levels of vitamin D metabolites have been assocaited with increased prostate cancer risk (23, 24).

Inflammation may also be involved in the link between BMD and prostate cancer. Pro-inflammatory cytokines have been implicated in osteoporosis and increased fractures risk (25, 26). Interleukin-6 was shown to stimulate osteoclasts, thereby increasing the rates of bone remodeling and bone loss (26). Inflammation is also suggested to play a role in prostate carcinogenesis (27). In the MrOS cohort, self-reported history of prostatitis was found to be positively associated with prevalent prostate cancer (OR= 5.4, 95% CI 4.4–6.6) (28), raising the possibility that chronic inflammation within the prostate may contribute to the pathogenesis of prostate cancer. Owing to the complex multifactorial pathogenesis of prostate cancer and bone mineralization, it is likely that more than one biological mechanism is involved in their link.

Interestingly, in our secondary analysis, the relationship of BMD measures with prostate cancer was limited to high-grade tumors. To our knowledge, our study was the first to investigate the relationship of bone density with tumor grade in PCa. Therefore, no other data are available for direct comparison. Vitamin D insufficiency may be involved in this association. Recent results from the Physician’s Health Study have indicated that low levels of both 25(OH) D and 1,25(OH)2D were related to increased risk of aggressive PCa. However, no such associations were observed for non-aggressive disease (24).

The observed associations were specific to BMD of the total body as well as the total hip, in the case of high grade disease. The lack of association at the spine may be related to the sensitivity of DXA technology to extra-osseous calcification, such as aortic calcification and degenarative osteoarthritic changes, which get incorporated in the region of interest and lead to a falsely increased bone density of the spine.

Interestingly, we observed that men with incident PCa were younger and had a higher level of physical activity, compared to men who did not develop the disease. This may be a reflection of increased awareness to prostate cancer screening and prevention in PCa cases, as triggered by their stronger family history of the disease.

Our results extend previous findings by longitudinally examining the association of prostate cancer and tumor grade with bone density at different skeletal sites. Our study had the benefit of a rigorous adjudication of prostate cancer cases, determination of BMD using a state-of-the-art method, and adjustment of results for a comprehensive set of risk factors for PCa, including calcium intake, family history of PCa, physical activity, and statin use. Limitations of our analysis include the unavailability of serum measurements of vitamin D and sex hormones on the full population, and the generalizability of findings to other populations due to the inclusion of a well-functioning cohort of mainly older white men.

In conclusion, we observed an unexpected inverse association between total body bone density and the risk of prostate cancer in older men who did not have a prior history of the disease. Further research is needed to confirm the direction of the relationship and to elucidate the pathophysiological mechanisms involved in the link between BMD and prostate cancer.

Acknowledgments

Financial Support: The Osteoporotic Fractures in Men (MrOS) Study is supported by National Institutes of Health funding. The following institutes provide support: the National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS), the National Institute on Aging (NIA), the National Center for Research Resources (NCRR), and NIH Roadmap for Medical Research under the following grant numbers: U01 AR45580, U01 AR45614, U01 AR45632, U01 AR45647, U01 AR45654, U01 AR45583, U01 AG18197, U01-AG027810, and UL1 RR024140.

Footnotes

Conflict of Interest: None

Reference List

1. Jemal A, Siegel R, Ward E, et al. Cancer statistics, 2007. CA Cancer J Clin. 2007;57:43–66. [PubMed]
2. Gann PH, Hennekens CH, Ma J, Longcope C, Stampfer MJ. Prospective study of sex hormone levels and risk of prostate cancer. J Natl Cancer Inst. 1996;88:1118–26. [PubMed]
3. Parsons JK, Carter HB, Platz EA, et al. Serum testosterone and the risk of prostate cancer: potential implications for testosterone therapy. Cancer Epidemiol Biomarkers Prev. 2005;14:2257–60. [PubMed]
4. Renehan AG, Zwahlen M, Minder C, et al. Insulin-like growth factor (IGF)-I, IGF binding protein-3, and cancer risk: systematic review and meta-regression analysis. Lancet. 2004;363:1346–53. [PubMed]
5. Chan JM, Giovannucci EL. Dairy products, calcium, and vitamin D and risk of prostate cancer. Epidemiol Rev. 2001;23:87–92. [PubMed]
6. Wiren S, Stocks T, Rinaldi S, et al. Androgens and prostate cancer risk: A prospective study. Prostate. 2007;67:1230–7. [PubMed]
7. Allen NE, Key TJ, Appleby PN, et al. Serum insulin-like growth factor (IGF)-I and IGF-binding protein-3 concentrations and prostate cancer risk: results from the European Prospective Investigation into Cancer and Nutrition. Cancer Epidemiol Biomarkers Prev. 2007;16:1121–7. [PubMed]
8. Zmuda JM, Cauley JA, Ljung BM, et al. Bone mass and breast cancer risk in older women: differences by stage at diagnosis. J Natl Cancer Inst. 2001;93:930–6. [PubMed]
9. Cauley JA, Lucas FL, Kuller LH, et al. Bone mineral density and risk of breast cancer in older women: the study of osteoporotic fractures. Study of Osteoporotic Fractures Research Group. JAMA. 1996;276:1404–8. [PubMed]
10. Bunker CH, Zmuda JM, Patrick AL, et al. High bone density is associated with prostate cancer in older Afro-Caribbean men: Tobago prostate survey. Cancer Causes Control. 2006;17:1083–9. [PubMed]
11. Nelson RL, Turyk M, Kim J, Persky V. Bone mineral density and the subsequent risk of cancer in the NHANES I follow-up cohort. BMC Cancer. 2002;2:22. [PMC free article] [PubMed]
12. Zhang Y, Kiel DP, Ellison RC, et al. Bone mass and the risk of prostate cancer: the Framingham Study. Am J Med. 2002;113:734–9. [PubMed]
13. Orwoll E, Blank JB, Barrett-Connor E, et al. Design and baseline characteristics of the osteoporotic fractures in men (MrOS) study--a large observational study of the determinants of fracture in older men. Contemp Clin Trials. 2005;26:569–85. [PubMed]
14. Blank JB, Cawthon PM, Carrion-Petersen ML, et al. Overview of recruitment for the osteoporotic fractures in men study (MrOS) Contemporary Clinical Trials. 2005;26:557–68. [PubMed]
15. Washburn RA, Smith KW, Jette AM, Janney CA. The Physical Activity Scale for the Elderly (PASE): development and evaluation. J Clin Epidemiol. 1993;46:153–62. [PubMed]
16. Block G, Subar AF. Estimates of nutrient intake from a food frequency questionnaire: the 1987 National Health Interview Survey. J Am Diet Assoc. 1992;92:969–77. [PubMed]
17. Looker AC, Wahner HW, Dunn WL, et al. Updated data on proximal femur bone mineral levels of US adults. Osteoporos Int. 1998;8:468–89. [PubMed]
18. Riggs BL, Khosla S, Melton LJ., III Sex steroids and the construction and conservation of the adult skeleton. Endocr Rev. 2002;23:279–302. [PubMed]
19. Khosla S. Role of hormonal changes in the pathogenesis of osteoporosis in men. Calcif Tissue Int. 2004;75:110–3. [PubMed]
20. Mellstrom D, Johnell O, Ljunggren O, et al. Free testosterone is an independent predictor of BMD and prevalent fractures in elderly men: MrOS Sweden. J Bone Miner Res. 2006;21:529–35. [PubMed]
21. Raisz LG. Pathogenesis of osteoporosis: concepts, conflicts, and prospects. J Clin Invest. 2005;115:3318–25. [PMC free article] [PubMed]
22. Krishnan AV, Moreno J, Nonn L, et al. Novel pathways that contribute to the anti-proliferative and chemopreventive activities of calcitriol in prostate cancer. J Steroid Biochem Mol Biol. 2007;103:694–702. [PubMed]
23. Corder EH, Guess HA, Hulka BS, et al. Vitamin D and prostate cancer: a prediagnostic study with stored sera. Cancer Epidemiol Biomarkers Prev. 1993;2:467–72. [PubMed]
24. Li H, Stampfer MJ, Hollis JB, et al. A prospective study of plasma vitamin D metabolites, vitamin D receptor polymorphisms, and prostate cancer. PLoS Med. 2007;4:e103. [PMC free article] [PubMed]
25. Cauley JA, Danielson ME, Boudreau RM, et al. Inflammatory markers and incident fracture risk in older men and women: the health aging and body composition study. J Bone Miner Res. 2007;22:1088–95. [PubMed]
26. Papanicolaou DA, Wilder RL, Manolagas SC, Chrousos GP. The pathophysiologic roles of interleukin-6 in human disease. Ann Intern Med. 1998;128:127–37. [PubMed]
27. De Marzo AM, Platz EA, Sutcliffe S, et al. Inflammation in prostate carcinogenesis. Nat Rev Cancer. 2007;7:256–69. [PubMed]
28. Daniels NA, Ewing SK, Zmuda JM, Wilt TJ, Bauer DC. Correlates and prevalence of prostatitis in a large community-based cohort of older men. Urology. 2005;66:964–70. [PubMed]