Email updates

Keep up to date with the latest news and content from Breast Cancer Research and BioMed Central.

Open Access Highly Accessed Research article

Plasma 25-hydroxyvitamin D and risk of breast cancer in the Nurses' Health Study II

A Heather Eliassen12*, Donna Spiegelman23, Bruce W Hollis5, Ronald L Horst6, Walter C Willett124 and Susan E Hankinson12

Author Affiliations

1 Channing Laboratory, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, 181 Longwood Avenue, Boston, MA 02115, USA

2 Department of Epidemiology, Harvard School of Public Health, 677 Huntington Avenue, Boston, MA 02115, USA

3 Department of Biostatistics, Harvard School of Public Health, 677 Huntington Avenue, Boston, MA 02115, USA

4 Department of Nutrition, Harvard School of Public Health, 677 Huntington Avenue, Boston, MA 02115, USA

5 Department of Pediatrics and Biochemistry and Molecular Biology, Medical University of South Carolina, 171 Ashley Avenue, Charleston, SC 29425, USA

6 Heartland Assays, Inc., 2325 North Loop Drive, Suite 6300, Ames, IA 50010, USA

For all author emails, please log on.

Breast Cancer Research 2011, 13:R50  doi:10.1186/bcr2880


The electronic version of this article is the complete one and can be found online at: http://breast-cancer-research.com/content/13/3/R50


Received:21 January 2011
Revisions received:1 April 2011
Accepted:11 May 2011
Published:11 May 2011

© 2011 Eliassen et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Introduction

Experimental evidence indicates vitamin D may play an important role in breast cancer etiology but epidemiologic evidence to date is inconsistent. Vitamin D comes from dietary intake and sun exposure and plasma levels of 25-hydroxyvitamin D (25(OH)D) are considered the best measure of vitamin D status.

Methods

We conducted a prospective nested case-control study within the Nurses' Health Study II (NHSII). Plasma samples collected in 1996 to 1999 were assayed for 25(OH)D in 613 cases, diagnosed after blood collection and before 1 June 2007, and in 1,218 matched controls. Multivariate relative risks (RR) and 95% confidence intervals (CI) were calculated by conditional logistic regression, adjusting for several breast cancer risk factors.

Results

No significant association was observed between plasma 25(OH)D levels and breast cancer risk (top vs. bottom quartile multivariate RR = 1.20, 95% CI (0.88 to 1.63), P-value, test for trend = 0.32). Results were similar when season-specific quartile cut points were used. Results did not change when restricted to women who were premenopausal at blood collection or premenopausal at diagnosis. Results were similar between estrogen receptor (ER)+/progesterone receptor (PR)+ and ER-/PR- tumors (P-value, test for heterogeneity = 0.51). The association did not vary by age at blood collection or season of blood collection, but did vary when stratified by body mass index (P-value, test for heterogeneity = 0.01).

Conclusions

Circulating 25(OH)D levels were not significantly associated with breast cancer risk in this predominantly premenopausal population.

Introduction

Vitamin D is an essential component of the endocrine system, responsible for calcium level maintenance, and is hypothesized to be involved in many other body systems. There are two circulating forms of vitamin D: 25-hydroxyvitamin D (25(OH)D) and 1,25-dihydroxyvitamin D3 (1,25(OH)2D). 25(OH)D is the main form of circulating vitamin D that is converted, by the 1α hydroxylase enzyme, into 1,25(OH)2D, which is considered the biologically active metabolite, binding to nuclear vitamin D receptors (VDR). However, although 1,25(OH)2D binds to VDR with 1,000 times the affinity of 25(OH)D, the latter is more abundant suggesting its interaction with VDR may not be trivial [1]. Although the conversion of 25(OH)D primarily takes place in the kidney, 1α hydroxylase also is present in breast and other tissues, suggesting that some portion of 1,25(OH)2D is produced and used locally [2-6]. 1,25(OH)2D is tightly regulated to maintain calcium homeostasis, while 25(OH)D serves as a pool for biologically active vitamin D and thus is a marker of overall vitamin D status [7,8].

Experimental evidence supports an important role of vitamin D in breast cancer etiology with 1,25(OH)2D reducing proliferation and promoting differentiation and apoptosis in breast cancer cell lines [9-11]. In animal models, vitamin D inhibited growth of both estrogen receptor (ER)+ and ER- tumors [12-14]. In a recent study in breast cancer cells, 1,25(OH)2D decreased the expression of aromatase [15]. Although aromatase in adipose tissue plays an integral role in estrogen levels in postmenopausal women, aromatase inhibitors may increase ovarian estradiol production in premenopausal women [16-19].

Vitamin D, which is hydroxylated in the liver into 25(OH)D, comes from dietary and supplement intake as well as sun exposure. Although several studies have assessed sun exposure and dietary intake separately, plasma levels of 25(OH)D are considered a better marker of vitamin D status because they incorporate all sources of exposure as well as skin color, weight, and other factors that influence an individual's vitamin D levels [20]. The association between circulating 25(OH)D and breast cancer has been assessed in seven prospective studies with conflicting results [21-27]. Most studies have included primarily or exclusively postmenopausal women, with the exception of a small Finnish study (n = 100 cases) in which plasma 25(OH)D levels measured during pregnancy were not significantly associated with subsequent breast cancer (RR = 1.4, 95% CI (0.6 to 3.4)) [25]. The inconsistent results overall and limited data in premenopausal women, coupled with the potential public health importance of vitamin D as a chemopreventive agent, prompted our assessment of this association in the prospective Nurses' Health Study II (NHSII) cohort.

Materials and methods

Study population

The NHSII was established in 1989, when 116,430 female registered nurses, aged 25 to 42 years, completed and returned a questionnaire. The cohort has been followed biennially by questionnaire to update exposures and ascertain newly diagnosed disease.

Between 1996 and 1999, 29,611 cohort members, who were cancer-free and between the ages of 32 and 54 years, provided blood and urine samples. These women were similar to the overall cohort with respect to lifestyle factors, such as body mass index, parity, age at menarche, past oral contraceptive use, and only differed slightly in the prevalence of family history of breast cancer (19% vs. 15% in the overall cohort) [28]. Of the 29,611 women who gave blood, 18,521 were premenopausal participants (that is, still having menstrual periods) who provided two blood samples and one urine sample timed within the menstrual cycle (one follicular sample collected on the third to fifth day and one luteal sample collected seven to nine days before the anticipated start of their next cycle). Participants were sent a short questionnaire and a blood and urine collection kit containing the necessary supplies to have blood samples drawn by a local laboratory or a colleague. Follicular plasma was aliquotted by the participants 8 to 24 hours after collection and stored in their home freezer until the luteal collection. The day of the luteal collection, follicular and luteal blood samples and luteal urine samples were shipped, via overnight courier with an ice-pack, to our laboratory where the luteal blood sample was processed and separated into plasma, and red blood cell and white blood cell components. Samples have been stored in liquid nitrogen freezers (<-130°C) since collection. The study was approved by the Committee on the Use of Human Subjects in Research at Harvard School of Public Health and Brigham and Women's Hospital. Informed consent was implied by receipt of completed questionnaires and blood samples.

Cases

Breast cancer cases were identified on biennial questionnaires; the National Death Index was searched for non-responders. Cases had no previously reported cancer diagnosis and were diagnosed with breast cancer after blood collection but before 1 June 2007. Overall, 613 cases of breast cancer (n = 415 invasive) were reported on biennial questionnaires and confirmed by medical record review (n = 582) or verbal confirmation by the nurse (n = 31). Given the 99% confirmation rate on medical record review, these latter cases were included. Information on invasiveness and hormone receptor status was abstracted from the medical record. Mean time from blood draw to diagnosis was 57 months (range = 1 to 127). Two controls were matched to each case (n = 1,218) on age (± 2 years); menopausal status at diagnosis; month/year of collection (± 2 months); ethnicity (African-American, Asian, Hispanic, Caucasian, Other); luteal day ((date of next period-date of luteal blood draw) ± 1 day); and for each blood collection, time of day (± 2 hours) and fasting status (< 2 hours, 2 to 4, 5 to 7, 8 to 11, 12+).

Laboratory assays

25(OH)D was assayed at Heartland Assays, Inc. (Ames, IA, USA) using a radioimmunoassay with radioiodinated tracers after acetonitrile extraction. Case-control sets were assayed together and samples were ordered randomly and labeled to mask case-control status. Samples were assayed in two batches: cases (and matched controls) diagnosed between 1999 and 2005 and those diagnosed from 2005 to 2007. The overall coefficients of variation (CVs) from masked replicate quality control samples included in each batch were 10.7% and 6.0%. Due to technical difficulties with the assay, three case-control sets failed and were removed, resulting in our final study population of 613 cases and 1,218 controls.

Covariate data

We obtained breast cancer risk factor information from a questionnaire completed at the time of blood collection and from the biennial NHSII questionnaires. Age at menarche, weight at age 18, and height were queried in 1989. Age at first birth, parity, and diagnosis of benign breast disease were assessed biennially. Family history of breast cancer in the participant's mother and sisters was queried in 1989 and 1997. Weight was assessed at the time of blood collection.

Statistical analysis

Quartile cut points were defined in two ways, based on the distribution in controls. The first set of cut points was calculated among all controls. Next we calculated season-specific cut points, using November to April blood collection dates to define "winter" and May to October blood collection dates to define "summer". We identified statistical outliers using the generalized extreme studentized deviate many-outlier detection approach [29]. Although we identified seven participants who had outlier values on the low end of the 25(OH)D distribution (≤3.6 ng/ml (8.99 nmol/L)), we affirmed that these were likely truly low levels given that most of these blood samples were collected in winter months and the participants had very low vitamin D dietary intake (for example, 50 to 100 IU/day). Therefore, we retained these values in the analysis.

We used conditional logistic regression to estimate relative risks (RRs) and 95% confidence intervals (CIs) to adjust for the matched design of this study. Multivariate models adjusted for body mass index (BMI) at age 18 and at the time of blood collection, ages at menarche and first birth, parity, family history of breast cancer, and history of benign breast disease. Given that outdoor physical activity may serve as a proxy for sun exposure, an important source of vitamin D, our main model did not include physical activity but we adjusted for it in a sensitivity analysis. In stratified analyses, we used unconditional logistic regression, adjusting for matching factors, since results from multivariate unconditional and conditional logistic regression models were essentially identical. To test whether the association differed by ER and PR status of the tumor, we used polychotomous logistic regression [30] with three endpoints (for example, ER+/PR+, ER-/PR-, and no breast cancer). We used a likelihood ratio test to compare a model with separate 25(OH)D slopes in each case group with a model with a common slope. Wald tests for interaction between stratification variables and hormones were used to compare the slope of the quartile medians between groups. Tests for trend were conducted by modeling quartile median concentrations and calculating the Wald statistic. The shape of the dose-response curves and tests for non-linearity were assessed using restricted cubic spline models [31,32]. We used previously published reproducibility data [33] to correct for random within-person and laboratory error [34]. All P-values were based on two-sided tests and were considered statistically significant if ≤0.05. All analyses were conducted using SAS software version 9 (SAS Institute, Cary, NC, USA).

Results

Cases had slightly lower BMI at age 18 and at blood collection and had gained less weight since age 18 than controls (Table 1). Cases were more likely than controls to be nulliparous and have a family history of breast cancer and a history of benign breast disease.

Table 1. Characteristics of breast cancer cases and matched controls, NHSII; mean ± SD or %.

In the simple conditional logistic regression model accounting for matching factors only, women in the highest quartile of 25(OH)D were not at significantly different risk of breast cancer compared with women in the lowest quartile (RR = 1.26, 95% CI (0.94 to 1.69), P-value, test for trend = 0.14) (Table 2). After adjustment for breast cancer risk factors, the RR was slightly attenuated (RR = 1.20, 95% CI (0.88 to 1.63), P-value, test for trend = 0.32), with adjustment for BMI at blood collection accounting for the majority of the difference between the simple and multivariate estimates. To assess the impact of physical activity, which could be both a confounder and a source of vitamin D, we added it to our multivariate model; results still were not significant (top vs. bottom quartile RR = 1.30, 95% CI (0.95 to 1.78), P-value, test for trend = 0.14). Analyses using season-specific cut points were similarly null (comparable RR = 1.06, 95% CI (0.78 to 1.44), P-value, test for trend = 0.79). Using restricted cubic spline models, no significant linear or non-linear associations were detected.

Table 2. Relative risks (95% CI) of breast cancer according to quartile of plasma 25(OH)D, NHSII.

Analyses by tumor subtype did not show any significant differences (Table 2). When we restricted the analyses to invasive cases (N = 415), results were similar to the overall (top vs. bottom RR = 1.29, 95% CI (0.92 to 1.81), P-value, test for trend = 0.14). RRs also were similar across hormone receptor subtypes (ER+ vs. ER- P-value, test for heterogeneity = 0.41; ER+/PR+ vs. ER-/PR- P-value, test for heterogeneity = 0.51).

To assess whether the association between 25(OH)D levels and breast cancer risk differed by other factors, we conducted secondary analyses restricted to premenopausal women and women diagnosed more than two years after blood collection. Results of analyses restricted to women who were premenopausal at blood collection (N = 470 cases) were still null (top vs. bottom quartile RR = 1.23, 95% CI (0.87 to 1.73), P-value, test for trend = 0.29), as were results of analyses restricted to women premenopausal at diagnosis (N = 294 cases) (comparable RR = 1.40, 95% CI (0.91 to 2.15), P-value, test for trend = 0.16). Results also were unchanged when we excluded cases diagnosed in the first two years (RR = 1.04, 95% CI (0.74 to 1.47), P-value, test for trend = 0.91).

We also stratified by age at blood collection (< 45 vs. ≥ 45), BMI at blood collection (< 25 vs. ≥ 25 kg/m2), and season of blood collection (winter vs. summer) (Table 3). When we stratified by BMI, we observed no significant association among leaner women (top vs. bottom quartile RR = 0.90, 95% CI (0.60 to 1.33), P-value, test for trend = 0.45) and a significant positive association among overweight and obese women (comparable RR = 1.90, 95% CI (1.19 to 3.03), P-value, test for trend = 0.005). The interaction between plasma 25(OH)D and BMI was statistically significant (P-value, test for heterogeneity = 0.006) and a significant linear association was observed in restricted cubic spline models for BMI ≥ 25 (P = 0.007) but no significant association for BMI < 25 (P = 0.09). No non-linearity of 25(OH)D in either BMI strata was evident. Results were unchanged when stratified by age at blood collection (P-value, test for heterogeneity = 0.29).

Table 3. Stratified relative risks (95% CI) of breast cancer according to quartile of plasma 25(OH)D, NHSII.

Additional stratified analyses included family history of breast cancer and among blood samples collected from premenopausal women in the luteal phase of the menstrual cycle. Although there was a significant interaction by family history of breast cancer (P-value, test for heterogeneity < 0.001), the numbers are small and the effects in each strata of family history are not clear (for example, quartiles 2, 3, 4 vs. quartile 1 RR = 1.15, 1.14, 1.19, P-value, test for trend = 0.33 for no family history and RR = 0.54, 0.32, 1.49, P-value, test for trend = 0.31 for positive family history). Results among women with luteal phase samples were similar to overall results (data not shown).

The association between 25(OH)D and breast cancer did not vary by season of blood collection using overall cut points (top vs. bottom quartile summer RR = 1.17, 95% CI (0.74 to 1.88), P-value, test for trend = 0.39, winter RR = 1.23, 95% CI (0.83 to 1.84), P-value, test for trend = 0.44; P-value, test for heterogeneity = 0.33) or season-specific cut points (top vs. bottom quartile summer RR = 0.97, 95% CI (0.63 to 1.52), P-value, test for trend = 0.95, winter RR = 1.22, 95% CI (0.81 to 1.82), P-value, test for trend = 0.52) (Table 3).

When we corrected for random within-person and laboratory error using the reproducibility data (intraclass correlation coefficient (ICC) = 0.72), results were largely unchanged (for example, median of top vs. bottom quartile RR = 1.09, 95% CI (0.86 to 1.36) uncorrected and RR = 1.11, 95% CI (0.83 to 1.48) corrected).

Discussion

In this analysis among predominantly premenopausal women, we did not observe an inverse association between plasma 25(OH)D levels and breast cancer risk overall. These results did not change when restricted to invasive cases, cases diagnosed more than two years after blood collection, or women who were premenopausal at diagnosis. Results were not substantially different between strata defined by age at blood collection or season of blood collection. Stratification by BMI yielded a significant interaction, with a significant positive association observed among overweight and obese women.

To date, only seven prospective studies of circulating 25(OH)D levels and breast cancer risk have been conducted [21-27]. Results have been inconsistent, though there have been few studies that have been able to examine this association among premenopausal women. In the NHS among mostly postmenopausal women, increasing 25(OH)D was suggestively inversely associated with breast cancer risk (top vs. bottom quintile RR = 0.73, 95% CI (0.49 to 1.07), P-value, test for trend = 0.06) [21]. Results from the Women's Health Initiative among postmenopausal women showed a non-significant inverse association between plasma 25(OH)D levels and breast cancer risk (top vs. bottom quintile RR = 0.82, P-value, test for trend = 0.20) [23]. However, other recent reports among postmenopausal women from the Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial and the Cancer Prevention Study II (CPSII) showed no associations (top vs. bottom quintile RR = 1.04 (P-value, test for trend = 0.81), 1.09 (P-value, test for trend = 0.6), respectively) [22,24]. In the recent Mälmo Diet and Cancer Study, no association with total 25(OH)D was observed overall (top vs. bottom quartile RR = 0.96, P-value, test for trend = 0.78) and a non-significant positive association was observed among the 196 premenopausal cases (top vs. bottom quartile RR = 1.74, 95% CI (0.84 to 3.60), P-value, test for trend = 0.14) [26]. In contrast, the suggested inverse association with serum 25(OH)D in the E3N cohort was stronger among the small numbers of premenopausal women (N = 54) (top vs. bottom tertile OR = 0.37, 95% CI (0.12 to 1.15), P-trend = 0.11) [27]. In a small Finnish study of premenopausal women (n = 100 cases), plasma 25(OH)D levels measured during pregnancy were not significantly associated with breast cancer risk (top vs. bottom quartile RR = 1.4, 95% CI (0.6 to 3.4), P-value, test for trend = 0.4) [25]. Our current study, the largest to date among premenopausal women, shows no significant association between plasma 25(OH)D levels and breast cancer risk.

The significant interaction we observed by BMI was unexpected. Although a significant interaction by BMI also was observed in the NHS study of 25(OH)D and colon cancer, in that study, a stronger inverse association was observed among leaner women [35]. Five prior studies of the association between 25(OH)D levels and breast cancer risk, again among primarily postmenopausal women, reported no significant difference by BMI [21,22,24,26,27]. Upon further exploration, we noted that those with higher 25(OH)D levels in the BMI ≥ 25 group had higher levels of physical activity than those with lower 25(OH)D levels. Although one hypothesis is that this could reflect fewer anovulatory cycles in overweight and obese women with high vitamin D, physical activity was inversely associated with breast cancer risk in premenopausal women in this cohort and there was not a significant difference in this association by BMI [36]. Adjusting for confounding by physical activity in the analysis yielded slightly stronger RRs. Another hypothesis is that overweight women with higher 25(OH)D levels do not benefit from the reduced risk of breast cancer associated with higher BMI in premenopausal women. This is possible given that 25(OH)D levels may inhibit aromatase, which in turn could lead to increased ovarian estrogen production in premenopausal women [15-19]. It also is possible that this finding of a significant interaction is due to chance; this warrants further exploration in other studies.

The prior analysis in the NHS suggested that higher 25(OH)D levels may be associated with a lower risk of ER- or ER-/PR- disease, but no difference by hormone receptor status was observed in the CPSII cohort [21,24]. In this study of predominantly premenopausal women, we did not observe significant differences by hormone receptor status.

This study has several strengths, including the large number of premenopausal cases, extensive information on potential confounders, and blood samples collected before diagnosis. Although this analysis includes only one measure of plasma 25(OH)D to represent longer-term exposure, we have showed good reproducibility of this analyte over a two- to three-year period in the NHS (ICC = 0.72) [33].

Conclusions

While higher plasma 25(OH)D levels have been hypothesized to reduce breast cancer risk, our results suggest such a benefit may not be observed among premenopausal women.

Abbreviations

25(OH)D: 25-hydroxyvitamin D; 1,25(OH)2D: 1,25-dihydroxyvitamin D3; BMI: body mass index; CI: confidence interval; CPSII: Cancer Prevention Study II; CV: coefficient of variation; ER: estrogen receptor; ICC: intraclass correlation coefficient; NHS: Nurses' Health Study; NHSII: Nurses' Health Study II; PR: progesterone receptor; RR: relative risk

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AHE conceived of the study design, acquired and analyzed the data, interpreted the data and wrote the manuscript. DS provided statistical expertise, contributed to the analysis and interpretation of the data, and critically reviewed the manuscript. BWH assisted in the interpretation of the data and critically reviewed the manuscript. RLH performed the 25(OH)D assay and critically reviewed the manuscript. WCW and SEH conceived of the study design, secured funding, acquired the data, contributed to the analysis and interpretation of the data, and critically reviewed the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This study was supported by Research Grants CA67262, CA50385, and P50 CA089393 from the National Cancer Institute.

We would like to thank the participants of the Nurses' Health Study II for their continuing contributions. We also thank the following state cancer registries for their help: AL, AZ, AR, CA, CO, CT, DE, FL, GA, ID, IL, IN, IA, KY, LA, ME, MD, MA, MI, NE, NH, NJ, NY, NC, ND, OH, OK, OR, PA, RI, SC, TN, TX, VA, WA, WY.

References

  1. Ray R: Molecular recognition and structure-activity relations in vitamin D-binding protein and vitamin D receptor. In Vitamin D: Molecular Biology, Physiology, and Clinical Applications. Edited by Holick MF. Totowa, NJ: Humana Press; 1999:147-162. OpenURL

  2. Friedrich M, Rafi L, Mitschele T, Tilgen W, Schmidt W, Reichrath J: Analysis of the vitamin D system in cervical carcinomas, breast cancer and ovarian cancer.

    Recent Results Cancer Res 2003, 164:239-246. PubMed Abstract OpenURL

  3. Welsh J, Wietzke JA, Zinser GM, Byrne B, Smith K, Narvaez CJ: Vitamin D-3 receptor as a target for breast cancer prevention.

    J Nutr 2003, 133:2425S-2433S. PubMed Abstract | Publisher Full Text OpenURL

  4. Chen TC, Wang L, Whitlatch LW, Flanagan JN, Holick MF: Prostatic 25-hydroxyvitamin D-1alpha-hydroxylase and its implication in prostate cancer.

    J Cell Biochem 2003, 88:315-322. PubMed Abstract | Publisher Full Text OpenURL

  5. Tangpricha V, Flanagan JN, Whitlatch LW, Tseng CC, Chen TC, Holt PR, Lipkin MS, Holick MF: 25-hydroxyvitamin D-1alpha-hydroxylase in normal and malignant colon tissue.

    Lancet 2001, 357:1673-1674. PubMed Abstract | Publisher Full Text OpenURL

  6. Holick MF: Vitamin D: importance in the prevention of cancers, type 1 diabetes, heart disease, and osteoporosis.

    Am J Clin Nutr 2004, 79:362-371. PubMed Abstract | Publisher Full Text OpenURL

  7. Reichel H, Koeffler HP, Norman AW: The role of the vitamin D endocrine system in health and disease.

    N Engl J Med 1989, 320:980-991. PubMed Abstract | Publisher Full Text OpenURL

  8. Hollis BW: Assessment of vitamin D nutritional and hormonal status: what to measure and how to do it.

    Calcif Tissue Int 1996, 58:4-5. PubMed Abstract | Publisher Full Text OpenURL

  9. Welsh J, Wietzke JA, Zinser GM, Smyczek S, Romu S, Tribble E, Welsh JC, Byrne B, Narvaez CJ: Impact of the Vitamin D3 receptor on growth-regulatory pathways in mammary gland and breast cancer.

    J Steroid Biochem Mol Biol 2002, 83:85-92. PubMed Abstract | Publisher Full Text OpenURL

  10. Jensen SS, Madsen MW, Lukas J, Binderup L, Bartek J: Inhibitory effects of 1alpha,25-dihydroxyvitamin D(3) on the G(1)-S phase-controlling machinery.

    Mol Endocrinol 2001, 15:1370-1380. PubMed Abstract | Publisher Full Text OpenURL

  11. Colston KW, Berger U, Coombes RC: Possible role for vitamin D in controlling breast cancer cell proliferation.

    Lancet 1989, 1:188-191. PubMed Abstract OpenURL

  12. Colston KW, Pirianov G, Bramm E, Hamberg KJ, Binderup L: Effects of Seocalcitol (EB1089) on nitrosomethyl urea-induced rat mammary tumors.

    Breast Cancer Res Treat 2003, 80:303-311. PubMed Abstract | Publisher Full Text OpenURL

  13. Nolan E, Donepudi M, VanWeelden K, Flanagan L, Welsh J: Dissociation of vitamin D3 and anti-estrogen mediated growth regulation in MCF-7 breast cancer cells.

    Mol Cell Biochem 1998, 188:13-20. PubMed Abstract | Publisher Full Text OpenURL

  14. VanWeelden K, Flanagan L, Binderup L, Tenniswood M, Welsh J: Apoptotic regression of MCF-7 xenografts in nude mice treated with the vitamin D3 analog, EB1089.

    Endocrinology 1998, 139:2102-2110. PubMed Abstract | Publisher Full Text OpenURL

  15. Krishnan AV, Swami S, Feldman D: Vitamin D and breast cancer: inhibition of estrogen synthesis and signaling.

    J Steroid Biochem Mol Biol 2010, 121:343-348. PubMed Abstract | Publisher Full Text OpenURL

  16. Grodin JM, Siiteri PK, MacDonald PC: Source of estrogen production in postmenopausal women.

    J Clin Endocrinol Metab 1973, 36:207-214. PubMed Abstract | Publisher Full Text OpenURL

  17. Siiteri PK: Adipose tissue as a source of hormones.

    Am J Clin Nutr 1987, 45:277-282. PubMed Abstract | Publisher Full Text OpenURL

  18. Swaby RF, Sharma CG, Jordan VC: SERMs for the treatment and prevention of breast cancer.

    Rev Endocr Metab Disord 2007, 8:229-239. PubMed Abstract | Publisher Full Text OpenURL

  19. de Ziegler D, Mattenberger C, Luyet C, Romoscanu I, Irion NF, Bianchi-Demicheli F: Clinical use of aromatase inhibitors (AI) in premenopausal women.

    J Steroid Biochem Mol Biol 2005, 95:121-127. PubMed Abstract | Publisher Full Text OpenURL

  20. Giovannucci E, Liu Y, Rimm EB, Hollis BW, Fuchs CS, Stampfer MJ, Willett WC: Prospective study of predictors of vitamin D status and cancer incidence and mortality in men.

    J Natl Cancer Inst 2006, 98:451-459. PubMed Abstract | Publisher Full Text OpenURL

  21. Bertone-Johnson ER, Chen WY, Holick MF, Hollis BW, Colditz GA, Willett WC, Hankinson SE: Plasma 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D and risk of breast cancer.

    Cancer Epidemiol Biomarkers Prev 2005, 14:1991-1997. PubMed Abstract | Publisher Full Text OpenURL

  22. Freedman DM, Chang SC, Falk RT, Purdue MP, Huang WY, McCarty CA, Hollis BW, Graubard BI, Berg CD, Ziegler RG: Serum levels of vitamin D metabolites and breast cancer risk in the prostate, lung, colorectal, and ovarian cancer screening trial.

    Cancer Epidemiol Biomarkers Prev 2008, 17:889-894. PubMed Abstract | Publisher Full Text OpenURL

  23. Chlebowski RT, Johnson KC, Kooperberg C, Pettinger M, Wactawski-Wende J, Rohan T, Rossouw J, Lane D, O'Sullivan MJ, Yasmeen S, Hiatt RA, Shikany JM, Vitolins M, Khandekar J, Hubbell FA, Women's Health Initiative Investigators: Calcium plus vitamin D supplementation and the risk of breast cancer.

    J Natl Cancer Inst 2008, 100:1581-1591. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. McCullough ML, Stevens VL, Patel R, Jacobs EJ, Bain EB, Horst RL, Gapstur SM, Thun MJ, Calle EE: Serum 25-hydroxyvitamin D concentrations and postmenopausal breast cancer risk: a nested case control study in the Cancer Prevention Study-II Nutrition Cohort.

    Breast Cancer Res 2009, 11:R64. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  25. Agborsangaya CB, Surcel H-M, Toriola AT, Pukkala E, Parkkila S, Tuohimaa P, Lukanova A, Lehtinen M: Serum 25-hydroxyvitamin D at pregnancy and risk of breast cancer in a prospective study.

    Eur J Cancer 2010, 46:467-470. PubMed Abstract | Publisher Full Text OpenURL

  26. Almquist M, Bondeson A-G, Bondeson L, Malm J, Manjer J: Serum levels of vitamin D, PTH and calcium and breast cancer risk-a prospective nested case-control study.

    Int J Cancer 2010, 127:2159-2168. PubMed Abstract | Publisher Full Text OpenURL

  27. Engel P, Fagherazzi G, Boutten A, Dupré T, Mesrine S, Boutron-Ruault M-C, Clavel-Chapelon F: Serum 25(OH) vitamin D and risk of breast cancer: a nested case-control study from the French E3N cohort.

    Cancer Epidemiol Biomarkers Prev 2010, 19:2341-2350. PubMed Abstract | Publisher Full Text OpenURL

  28. Tworoger SS, Sluss P, Hankinson SE: Association between plasma prolactin concentrations and risk of breast cancer among predominately premenopausal women.

    Cancer Res 2006, 66:2476-2482. PubMed Abstract | Publisher Full Text OpenURL

  29. Rosner B: Percentage points for a generalized ESD many-outlier procedure.

    Technometrics 1983, 25:165-172. Publisher Full Text OpenURL

  30. Marshall RJ, Chisholm EM: Hypothesis testing in the polychotomous logistic model with an application to detecting gastrointestinal cancer.

    Stat Med 1985, 4:337-344. PubMed Abstract | Publisher Full Text OpenURL

  31. Durrleman S, Simon R: Flexible regression models with cubic splines.

    Stat Med 1989, 8:551-561. PubMed Abstract | Publisher Full Text OpenURL

  32. Govindarajulu US, Spiegelman D, Thurston SW, Ganguli B, Eisen EA: Comparing smoothing techniques in Cox models for exposure-response relationships.

    Stat Med 2007, 26:3735-3752. PubMed Abstract | Publisher Full Text OpenURL

  33. Kotsopoulos J, Tworoger SS, Campos H, Chung F-L, Clevenger CV, Franke AA, Mantzoros CS, Ricchiuti V, Willett WC, Hankinson SE, Eliassen AH: Reproducibility of plasma and urine biomarkers among premenopausal and postmenopausal women from the Nurses' Health Studies.

    Cancer Epidemiol Biomarkers Prev 2010, 19:938-946. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Rosner B, Spiegelman D, Willett WC: Correction of logistic regression relative risk estimates and confidence intervals for random within-person measurement error.

    Am J Epidemiol 1992, 136:1400-1413. PubMed Abstract OpenURL

  35. Wu K, Feskanich D, Fuchs CS, Willett WC, Hollis BW, Giovannucci EL: A nested case control study of plasma 25-hydroxyvitamin D concentrations and risk of colorectal cancer.

    J Natl Cancer Inst 2007, 99:1120-1129. PubMed Abstract | Publisher Full Text OpenURL

  36. Maruti SS, Willett WC, Feskanich D, Rosner B, Colditz GA: A prospective study of age-specific physical activity and premenopausal breast cancer.

    J Natl Cancer Inst 2008, 100:728-737. PubMed Abstract | Publisher Full Text OpenURL