Variants in the vitamin D pathway, serum levels of vitamin D, and estrogen receptor negative breast cancer among African-American women: a case-control study
- Song Yao1Email author,
- Gary Zirpoli1,
- Dana H Bovbjerg2,
- Lina Jandorf3,
- Chi Chen Hong1,
- Hua Zhao1,
- Lara E Sucheston1,
- Li Tang1,
- Michelle Roberts1,
- Gregory Ciupak1,
- Warren Davis1,
- Helena Hwang4,
- Candace S Johnson5,
- Donald L Trump6,
- Susan E McCann1,
- Foluso Ademuyiwa6,
- Karen S Pawlish7,
- Elisa V Bandera8 and
- Christine B Ambrosone1
© Yao et al.; licensee BioMed Central Ltd. 2012
Received: 15 November 2011
Accepted: 4 April 2012
Published: 4 April 2012
American women of African ancestry (AA) are more likely than European Americans (EA) to have estrogen receptor (ER)-negative breast cancer. 25-hydroxyvitamin D (25OHD) is low in AAs, and was associated with ER-negative tumors in EAs. We hypothesized that racial differences in 25OHD levels, as well as in inherited genetic variations, may contribute, in part, to the differences in tumor characteristics.
In a case (n = 928)-control (n = 843) study of breast cancer in AA and EA women, we measured serum 25OHD levels in controls and tested associations between risk and tag single nucleotide polymorphisms (SNPs) in VDR, CYP24A1 and CYP27B1, particularly by ER status.
More AAs had severe vitamin D deficiency (< 10 ng/ml) than EAs (34.3% vs 5.9%), with lowest levels among those with the highest African ancestry. Associations for SNPs differed by race. Among AAs, VDR SNP rs2239186, associated with higher serum levels of 25OHD, decreased risk after correction for multiple testing (OR = 0.53, 95% CI = 0.31-0.79, p by permutation = 0.03), but had no effect in EAs. The majority of associations were for ER-negative breast cancer, with seven differential associations between AA and EA women for CYP24A1 (p for interaction < 0.10). SNP rs27622941 was associated with a > twofold increased risk of ER-negative breast cancer among AAs (OR = 2.62, 95% CI = 1.38-4.98), but had no effect in EAs. rs2209314 decreased risk among EAs (OR = 0.38, 95% CI = 0.20-0.73), with no associations in AAs. The increased risk of ER-negative breast cancer in AAs compared to EAs was reduced and became non-significant (OR = 1.20, 95% CI = 0.80-1.79) after adjusting for these two CYP24A1 SNPs.
These data suggest that genetic variants in the vitamin D pathway may be related to the higher prevalence of ER-negative breast cancer in AA women.
American women of African ancestry (AA) are more likely to develop breast cancer at a younger age than those with European ancestry (EA) and are more likely to have tumors with aggressive characteristics, including high histological grade, negative estrogen receptor (ER) status, and basal-like - ER- and/or progesterone receptor (PR)-, HER2-, and cytokeratin 5/6+ and/or HER1+ -features [1, 2]. The reasons for these racial disparities are unknown.
It is clear that, among geographically diverse populations, certain genotypic and phenotypic characteristics may be selected for in response to local environmental pressures . Skin pigmentation, the primary factor that provides protection from ultraviolet (UV) radiation, is correlated with latitude, and dark skin pigmentation is likely to be the original ancestral trait in humans. Migrations to Europe and Asia eventually gave rise to decreased pigmentation and lighter skin [4–6]. As much as 90% of vitamin D is derived from sun exposure, but high skin melanin concentration prevents penetration of UVB light and compromises synthesis efficiency by 10 to 50 times . Although high pigmentation would reduce absorption of vitamin D, intense sun exposure in sub-Saharan Africa would compensate. However, in high-latitude areas where UVB intensity is low and where more time may be spent indoors (particularly in winter), vitamin D deficiency may result among individuals with higher skin pigmentation. Indeed, in the US, the prevalence of 25-hydroxyvitamin D (25OHD) of less than 15 ng/mL is almost 10 times higher in AA than in EA women , and the prevalence of severe vitamin D deficiency (< 10 ng/mL) among AAs was 29% in 2001 to 2004 . In contrast, in Guinea-Bissau, the average 25OHD levels in healthy Africans were 34 ng/mL, and the prevalence of severe vitamin D deficiency was as low as 1% .
Endogenous 25OHD may also be affected by variability in metabolic pathways, and synthesis and metabolism catalyzed by two major enzymes, 1α-hydroxylase and 24-hydroxylase, which are encoded by CYP27B1 and CYP24A1, respectively. Binding of vitamin D to the vitamin D receptor (VDR) activates or suppresses gene transcription, depending on the type of response elements , and genetic variability in the above genes, known to differ by ancestry , is likely to affect vitamin D signaling.
Laboratory, preclinical, and clinical findings support the hypothesis that low levels of vitamin D are related to breast cancer risk. In the human mammary gland, VDR is expressed in all cell types , and vitamin D treatment inhibits breast cancer cell proliferation, induces cell apoptosis, and prevents carcinogenesis in rodent models [14, 15]. However, epidemiologic evidence for associations between vitamin D and breast cancer risk is considered 'limited' , and one randomized trial showed little impact of vitamin D supplementation on breast cancer incidence . These inconclusive findings could be due to tumor heterogeneity, which implies that the effects of vitamin D may only present in specific breast cancer subtypes. In fact, Vdr knockout mice were more likely than their wild-type littermates to develop Er/Pr- tumors . Consistent with results from these preclinical studies, some epidemiologic studies have also indicated that effects of vitamin D may be strongest for breast cancers with poor prognostic characteristics and that lower serum 25OHD levels are found among women with ER- compared with ER+ tumors [19–21]. Recently, we reported lower levels of serum 25OHD in women with high- versus low-grade breast tumors and in women with triple-negative versus luminal A breast tumors .
Here, we examined levels of 25OHD in AA and EA women without breast cancer in relation to self-reported race as well as ancestry, which was estimated by using ancestry informative markers (AIMs). 25OHD levels in women with breast cancer could be a result of disease processes, and some samples were obtained after chemotherapy was initiated; thus, we did not compare serum levels of 25OHD between cases and controls. Instead, we evaluated variants in vitamin D activity and major metabolism (VDR, CYP27B1, and CYP24A1) in relation to breast cancer risk, particularly in relation to self-reported race and estrogen receptor status. We also tested whether vitamin D-related genetic variants could explain, in part, the higher prevalence of ER- breast cancer among AA women.
Materials and methods
The Women's Circle of Health Study (WCHS) is an ongoing study designed specifically to examine the role of genetic and non-genetic factors in early/aggressive breast cancer in AA and EA women. Study design, enrollment, and collection of data and biospecimens have been described in detail previously . Briefly, women with diagnosed incident breast cancer were identified through both hospital-based case ascertainment in targeted hospitals that had large referral patterns of AAs in four boroughs of the metropolitan New York City area and population-based case ascertainment in seven counties in New Jersey through the New Jersey State Cancer Registry. The eligibility criteria for cases were the following: self-identified AA and EA women, 20 to 75 years of age at diagnosis, no previous history of cancer other than non-melanoma skin cancer, recent diagnosis of primary, histologically confirmed breast cancer, and English-speaking. EA women with breast cancer, more prevalent in the catchment area than AA cases, were randomly selected for recruitment and were matched by age and county to AA cases. Controls who did not have a history of diagnosis of any cancer other than non-melanoma skin cancer and who were living in the same area as cases were identified through random digit dialing and were matched to cases by self-reported race and 5-year age categories. After agreement to participate was obtained, in-person interviews were conducted to complete informed consent and to query participants on a number of potential risk factors, including medical history, family history of cancer, diet, physical activity, and other lifestyle factors. Anthropometric measures were taken, and biospecimens were collected. Blood samples were initially collected, but owing to logistical and cost constraints, we transitioned to saliva samples after the enrollment of approximately 850 participants. Permission to obtain pathology data, including ER status, and tumor tissue blocks was included in the informed consent form. The participation rates were 80.2% and 53.4% for AA cases and controls, respectively, and 80.0% and 48.9% for EA cases and controls, respectively. This study was approved by the institutional review boards of the Roswell Park Cancer Institute, the Cancer Institute of New Jersey, the Mount Sinai School of Medicine, and the participating hospitals in New York City. At the time of the genotyping (April 2010), DNA and data were available for 553 AA cases and 466 AA controls and 383 EA cases and 382 EA controls from the WCHS.
Selection of multi-population tag single-nucleotide polymorphisms
We used a two-step approach to select a set of multi-population tag single-nucleotide polymorphisms (SNPs) that represent common genetic variations - minor allele frequency of at least 0.05 - in VDR in both AA and EA populations. First, 122 SNPs in the VDR region plus 15-kb regions from both 3' and 5' ends were selected from HapMap  and other resequencing projects by using the Genome Variation Server at SeattleSNPs . These SNPs were then genotyped in 60 AA and 60 EA controls. With the TAGster program , 49 multi-population tag SNPs were subsequently selected for genotyping in the WCHS. Also among those VDR variants finally selected were commonly studied SNPs, including Cdx2 (rs11568820), Fok1 (rs2228570), Bsm1 (rs1544410), Apa1 (rs7975232), and Taq1 (rs731236). For CYP24A1, 15 multi-population tag SNPs were selected by using publicly available CEU (Utah residents with ancestry from Northern and Western Europe) and YRI (Yoruba in Ibadan) genotype data from HapMap. For CYP27B1, only one SNP had a minor allele frequency of at least 0.05 in the dbSNP database and was thus included. To control for potential bias due to population admixture and to examine serum 25OHD levels in relation to ancestry, a panel of 108 AIMs that have been shown to be effective in correcting for admixture in case-control studies  was chosen.
DNA was extracted from blood samples by using FlexiGene™ DNA kits (Qiagen, Valencia, CA, USA) in accordance with the instructions of the manufacturer and from saliva collected in Oragene™ kits (DNA Genotek Inc., Kanata, ON, Canada). Genomic DNA was evaluated and quantified by a Nanodrop UV spectrometer (Thermo Fisher Scientific Inc., Wilmington, DE, USA) and a PicoGreen-based fluorometric assay (Molecular Probes, now part of Invitrogen Corporation, Carlsbad, CA, USA) and stored at -80°C until analysis. Selected tag SNPs and AIMs were genotyped by an Illumina GoldenGate assay (Illumina Inc., San Diego, CA, USA) at the Genomics Core Facility at the Roswell Park Cancer Institute. Five percent duplicates and two sets of in-house trio samples were included for quality control purposes. The average successful genotyping rate for each sample and each SNP was at least 99%, and no SNPs violated Hardy-Weinberg equilibrium in controls or mendelian inheritance. Clustering plots of SNPs that were significant in the statistical analysis were manually re-inspected post hoc to ensure that the calls were robust.
Measurement of serum levels of 25OHD in WCHS controls
Serum samples were available from 242 AA and 187 EA women in the control group and were used to measure levels of 25OHD by immunochemiluminometric assay. The assay coefficient of variation was 10.3%.
Estimation of ancestry and comparison of serum levels of 25OHD
Individual ancestral proportions for EA and AA were estimated with the Bayesian Markov chain Monte Carlo clustering algorithm implemented in STRUCTURE 2.3 . We included publicly available genotypes from the YRI and CEU ancestral populations, and the program was run multiple times assuming K = 2 underlying ancestries. Women with more than 85% of genomic race other than the self-identified race were excluded from the analysis (n = 13); therefore, 547 AA cases and 461 AA controls and 381 EA cases and 382 EA controls were included in the final analyses for SNPs and breast cancer risk. There were no exclusions due to conflicting self-report and marker ancestry in the serum analyses performed in controls only. To compare serum levels of 25OHD between AA and EA women, least squares means and standard errors were calculated with adjustment for age, body mass index (BMI), and season of blood collection (four seasons). To test genotype-phenotype correlations between SNPs and 25OHD levels in controls, the Pearson correlation test was used and genotypes were coded as 0, 1, and 2 to reflect the number of copies of the minor allele.
Associations between breast cancer risk and genotypes and haplotypes
Descriptive variables were analyzed by Student t test or chi-squared test. All genotype analyses were performed for AA and EA populations separately. A genotypic (co-dominant) model was assumed for SNP effects. When genotype frequencies of the rare homozygote were not more than 5% in both populations, categories were collapsed (homozygote rare and heterozygote) for power considerations. To test whether there was a linear dose effect of the variant alleles (log-additive genetic model for trend test), SNPs were coded as 0, 1, and 2 as described above. Univariate single SNP analysis was performed, and resulting P values were plotted after logarithm transformation with accompanying linkage disequilibrium (LD) map by using the snp.plotter R package . Covariates, including age at diagnosis, BMI, European ancestry, family history of breast cancer, and education, were then adjusted in multivariate logistic regression models to derive odds ratios (ORs) and 95% confidence intervals (CIs).
For VDR and CYP24A1, haplotype structure was determined by using the method of Gabriel and colleagues , and each haplotype was tested in the regression model in comparison with all other haplotypes. For both single SNP and haplotype analyses, we controlled the family-wise error rate by using permutation testing (n = 10,000) as implemented in PLINK .
Potential modification of associations by race, menopausal status, and estrogen receptor status
To examine whether the associations of SNPs with breast cancer risk differed between AA and EA women, interaction with race was tested by including a race*SNP term in the logistic regression model without an estimate of ancestry. A similar approach was used to test modification effects by menopausal status.
To test whether selected SNPs contributed to differential risk of ER- breast cancer between AA and EA women, OR of ER- breast cancer by race was first estimated from a base model containing race and other covariates. SNPs that were differentially associated with ER- breast cancer risk by race (P for interaction was less than 0.10) were then entered in the base model and subjected to backward selection. A substantial reduction of OR (≥ 10%) for AA versus EA race by adding selected SNPs would indicate that those SNPs explained, in part, the higher risk of ER- breast cancer in AA than in EA women.
Descriptive characteristics of African-American and European-American women by case-control status in the Women's Circle of Health Study
(n = 461)
(n = 381)
(n = 382)
Body mass index, kg/m2
Percentage with European ancestry
Less than high school
College or above
Hormone replacement therapy
Serum levels of 25OHD
Associations between genetic variants and breast cancer risk by self-reported race
Single-nucleotide polymorphisms in VDR and CYP24A1 and differential associations with breast cancer risk between African-American and European-American women in the Women's Circle of Health Study
Among EA women, the VDR 'at-risk' G allele for rs2239186 was more common in EA women but was not associated with breast cancer risk (OR = 0.85, 95% CI = 0.62 to 1.17), nor were VDR rs12721364 SNPs. There were increases in risk by the VDR SNP rs11608702 and significant decreases in risk by two CYP24A1 variants: rs912505 and rs3787555; however, these did not remain significant after correction for multiple testing.
Results from haplotype analysis were consistent with those from single SNP analysis for VDR rs2239186. Among AA women, a G-G-G haplotype consisting of this SNP and two neighboring variants was associated with a decreased risk of breast cancer after adjustment for multiple testing (OR = 0.55, 95% CI = 0.38 to 0.81, P = 0.04) (Table S2 of Additional file 1). Among EA women, similar results were found for haplotypes containing rs11608702 in VDR and haplotypes containing rs3787555 in CYP24A1. The commonly studied haplotype in the 3' untranslated region of VDR consisting of Taq1, Apa1, and Bsm1 was not associated with breast cancer risk in AA women, but a modest decreased risk was observed in EA women and significance was marginal (OR = 0.82, 95% CI = 0.67 to 1.02).
Estrogen receptor-negative breast cancer and CYP24A1variants
Single-nucleotide polymorphisms in VDR and CYP24A1 and differential association with estrogen receptor-specific breast cancer risk among African-American and European-American women
OR (95% CI)
OR (95% CI)
Changes in risk of estrogen receptor-negative breast cancer by race with inclusion of single-nucleotide polymorphisms in CYP24A1
Adjusted OR (95% CI)
Race (AA versus EA)
Base model + SNPs
Race (AA versus EA)
rs2209314 (AG/GG versus AA)
rs2762941 (AG versus AA)
rs2762941 (GG versus AA)
Lastly, there were significant interactions for two SNPs in VDR with menopausal status (Table S5 of Additional file 1). The increased risk associated with rs886441 in AA women was restricted to pre-menopausal women (OR = 2.27, 95% CI = 1.32 to 3.90), and the increased risk associated with rs7975232 (Apa1) in EA women was restricted to post-menopausal women (OR = 2.24, 95% CI = 1.19 to 4.21).
In this study, we found that relationships between breast cancer risk and variants in genes associated with vitamin D activity and metabolism, VDR and CYP24A1, differed depending upon self-reported race and that associations were most notable for risk of ER- breast cancer in both AA and EA women. Importantly, we found that rs2209314 and rs2762941 in CYP24A1 contributed significantly to the higher risk of ER- breast cancer in AA than EA women. Among controls in the WCHS, serum levels of 25OHD were notably lower in AA women than EA women, the lowest levels were among women with the greatest African ancestry estimated by AIMs, and VDR expression levels, as estimated from published data on cultured lymphoblastoid cells , were also lower in AA women. In a previous study among EA women, we found that low 25OHD levels were associated with increased risk of ER- breast cancer, both in comparison with controls and with women with ER+ breast cancer . Given all of these data, it is possible that low 25OHD levels in AA women, coupled with unique 'at-risk' genetic variants, contribute, in part, to the higher prevalence of ER- breast cancer among AA women. If these potential associations were to be consistently observed in future studies, our results would support a public health effort for vitamin D supplementation to reduce risk of aggressive breast cancer among AA women.
The finding of an inverse association between African ancestry estimated by AIMs and blood 25OHD levels is consistent with a recent community cohort study of AA men and women . Our findings of extensive racial differences in allele frequencies and LD patterns for SNPs in VDR and CYP24A1 are also consistent with those from an earlier study . Previous studies on VDR polymorphisms and breast cancer risk have focused on only a few SNPs. However, we did not find any relationship with Fok1, Bsm1, or Taq1 or a three-SNP haplotype consisting of Bsm1, Apa1, and Taq1 in either AA or EA women. The variant homozygote of Apa1 was associated with increased risk of breast cancer in EA women, but the effect was limited to post-menopausal women. Increased risk of breast cancer was also reported in a previous study for Apa1 ; however, results in the literature are conflicting [35, 36]. The G allele of Cdx2 was associated with lower risk of breast cancer in AA women in our study and this finding was in contrast to the speculated functional alteration that the variant G allele resulted in lower binding of the Cdx2 protein and thus lower transcriptional activity of VDR [37, 38]. Similar to a study among women in Germany , our study found no association of Cdx2 with breast cancer risk in EA women.
To date, three studies have examined selected SNPs in VDR with breast cancer risk in both AA and EA women. Two of them did not find associations for Fok1, Bsm1, Bgl1 (rs739837), or the 3' untranslated region poly(A) microsatellite in either AA or EA women [40, 41], and a third study found increased risk by Bsm1 variant in EA but not AA women . However, none of the above studies examined the associations by ER status or used the systematic approach we employed to capture variation across the genes.
We found that four SNPs in VDR and two SNPs in CYP24A1 had differential associations with breast cancer by race (P for interaction was not more than 0.10). The fact that the associations were not consistent in AA and EA populations corroborates the differences in blood levels of vitamin D and frequency and LD pattern of vitamin D-related genetic variants, implying that the race-specific associations might be the result of gene-environment interactions. In further analyses stratified by ER status, one SNP in VDR and seven SNPs in CYP24A1 were specifically associated with ER- but not ER+ cancer risk, and the associations differed between AA and EA women. Controlling for the two SNPs in CYP24A1 in a multivariate model substantially reduced the increased ER- breast cancer risk associated with AA race and made the association no longer significant. This provides the first evidence supporting the contribution of vitamin D-related genetic variants to higher risk of more aggressive breast cancer in AA women.
We found significant associations between breast cancer risk and a number of tag SNPs in VDR without previously known functionality. However, in our analyses, SNP rs2239186 was associated with increased serum levels of 25OHD in AA women without breast cancer. This SNP and a haplotype containing it were also significantly associated with reduced breast cancer risk in AA women, irrespectively of ER status. This SNP resides in an intronic region and thus is unlikely to be the causal variant. However, it may be a marker for a causal SNP outside of the VDR gene. This SNP has not been implicated in other breast cancer studies but has been shown to be associated with reduced risk of colorectal cancer in individuals with low vitamin D intake  and was also implicated in type I diabetes . This SNP may warrant future replication and fine-mapping studies.
The two SNPs, rs2209314 and rs2762941, in CYP24A1 shown to be associated with racial differences in ER- breast cancer risk are intronic. Although these SNPs have not been implicated previously in breast cancer, elevated expression of CYP24A1 was found in breast cancer tissues , indicating a potential role in breast cancer etiology. We did not observe associations of these two SNPs in CYP24A1 with serum 25OHD levels in either AA or EA populations, indicating that these two SNPs themselves or linked causal variants may affect ER- breast cancer not through altering circulating 25OHD levels but availability of vitamin D in local mammary tissues.
One limitation of our study is the lack of validation for the significant findings. The number of patients with breast cancer and controls was relatively limited, especially after stratification by race and ER status. None of the associations with SNPs, except for rs2239186, remained significant after correction for multiple comparisons. We thus could not exclude the possibility of false-positive findings in our data. However, the fact that rs2239186 was associated with higher serum 25OHD levels in AA women as well as reduced breast cancer risk in this population is biologically coherent and reduces the likelihood of spurious findings for this VDR SNP. Another limitation is that only three genes in vitamin D-related pathways were included in this study. Although VDR, CYP27B1, and CYP24A1 are the three key genes in this pathway, genes encoding for some other vitamin D metabolizing enzymes, particularly GC encoding for vitamin D binding protein (which has been related to circulating vitamin D levels), may also be related to breast cancer risk and warrant further studies. AA women are more likely to develop breast cancer at a younger age than EAs; we enrolled all eligible AA women but randomly selected eligible EA women, frequency-matching by 5-year age categories. We also initially limited eligibility to women 65 years or younger because of low participation of older women without breast cancer to case-control studies. Thus, the overall study population is younger than that of some other studies. Although we found no evidence of modification effects by menopausal status for any but two SNPs (Table S5 of Additional file 1), the high proportion of pre-menopausal women in this study needs to be considered in relation to generalizability.
We found notable differences in blood levels of 25OHD and genetic variants in VDR and CYP24A1 between AA and EA women. Specifically, we found race-specific associations with breast cancer risk and these associations may be due to distinct genetic background and differences in 25OHD levels between the two populations. Our study provides evidence that variants in vitamin D-related genes may contribute to higher risk of ER- breast cancer in AA than EA women. Future studies are warranted to validate our findings and to investigate whether pre-diagnostic blood levels of 25OHD may also be related to racial differences in risk of ER- breast cancer.
ancestry informative marker
body mass index
Utah residents with ancestry from Northern and Western Europe
lymphoblastoid cell line
vitamin D receptor
Women's Circle of Health Study
Yoruba in Ibadan.
This work was supported by grants from the US Army Medical Research and Material Command (DAMD-17-01-1-0334), the National Cancer Institute (R01 CA100598 and P01 CA151135), and the Breast Cancer Research Foundation and a gift from the Philip L Hubbell family. SY was the recipient of a pre-doctoral training award from the Department of Defense Breast Cancer Research Program (W81XWH-08-1-0223). EVB was the recipient of a career transition award from the National Cancer Institute (K22 CA138563). Samples were stored and managed by the Roswell Park Cancer Institute (RPCI) DataBank and BioRepository, and genotyping was performed in the RPCI Genomics Core Facility; both are Cancer Center Support Grant-shared resources supported by P30 CA016056-32. The New Jersey State Cancer Registry (NJSCR) is a participant in the Centers for Disease Control and Prevention's National Program of Cancer Registries and is a National Cancer Institute Surveillance, Epidemiology, and End Results (SEER) Expansion Registry. The NJSCR is supported by the Centers for Disease Control and Prevention under cooperative agreement DP07-703 awarded to the New Jersey Department of Health & Senior Services. The collection of State of New Jersey cancer incidence data is also supported by the National Cancer Institute's SEER Program under contract N01-PC-95001-20. The funding agents played no role in design, in the collection, analysis, and interpretation of data, in the writing of the manuscript, or in the decision to submit the manuscript for publication.
We are grateful to the women who participated in this study and to colleagues, physicians, and clinical staff at participating hospitals in New York who facilitated identification and enrollment of cases into the study: Kandace Amend (i3 Drug Safety), Helena Furberg (Memorial Sloan-Kettering Cancer Center), Thomas Rohan and Joseph Sparano (Albert Einstein College of Medicine), Kitwaw Demissie (University of Medicine and Dentistry of New Jersey), Paul Tartter and Alison Estabrook (St. Luke's Roosevelt Hospital), James Reilly (Kings County Hospital Center), Benjamin Pace, George Raptis, and Christina Weltz (Mount Sinai School of Medicine), Maria Castaldi (Jacob Medical Center), Sheldon Feldman (New York-Presbyterian), and Margaret Kemeny (Queens Hospital Center). We also thank Dr. Eileen M. Dolan for the permission for using the expression data of vitamin D receptor in HapMap lymphoblastoid cell lines.
- Amend K, Hicks D, Ambrosone CB: Breast cancer in African-American women: differences in tumor biology from European-American women. Cancer Res. 2006, 66: 8327-8330. 10.1158/0008-5472.CAN-06-1927.PubMedView ArticleGoogle Scholar
- Ooi LL, Zhou H, Kalak R, Zheng Y, Conigrave AD, Seibel MJ, Dunstan CR: Vitamin D deficiency promotes human breast cancer growth in a murine model of bone metastasis. Cancer Res. 2010, 70: 1835-1844. 10.1158/0008-5472.CAN-09-3194.PubMedView ArticleGoogle Scholar
- Fumagalli M, Sironi M, Pozzoli U, Ferrer-Admettla A, Pattini L, Nielsen R: Signatures of environmental genetic adaptation pinpoint pathogens as the main selective pressure through human evolution. PLoS Genet. 2011, 7: e100235-View ArticleGoogle Scholar
- Yuen AW, Jablonski NG: Vitamin D: in the evolution of human skin colour. Med Hypotheses. 2010, 74: 39-44. 10.1016/j.mehy.2009.08.007.PubMedView ArticleGoogle Scholar
- Jablonski NG: The evolution of human skin and skin color. Ann Rev Anthropol. 2004, 33: 585-623. 10.1146/annurev.anthro.33.070203.143955.View ArticleGoogle Scholar
- Chaplin G, Jablonski NG: Vitamin D and the evolution of human depigmentation. Am J Phys Anthropol. 2009, 139: 451-461. 10.1002/ajpa.21079.PubMedView ArticleGoogle Scholar
- Clemens TL, Adams JS, Henderson SL, Holick MF: Increased skin pigment reduces the capacity of skin to synthesise vitamin D3. Lancet. 1982, 1: 74-76.PubMedView ArticleGoogle Scholar
- Nesby-O'Dell S, Scanlon KS, Cogswell ME, Gillespie C, Hollis BW, Looker AC, Allen C, Doughertly C, Gunter EW, Bowman BA: Hypovitaminosis D prevalence and determinants among African American and white women of reproductive age: third National Health and Nutrition Examination Survey, 1988-1994. Am J Clin Nutr. 2002, 76: 187-192.PubMedGoogle Scholar
- Ginde AA, Liu MC, Camargo CA: Demographic differences and trends of vitamin D insufficiency in the US population, 1988-2004. Arch Intern Med. 2009, 169: 626-632. 10.1001/archinternmed.2008.604.PubMedPubMed CentralView ArticleGoogle Scholar
- Wejse C, Olesen R, Rabna P, Kaestel P, Gustafson P, Aaby P, Andersen PL, Glerup H, Sodemann M: Serum 25-hydroxyvitamin D in a West African population of tuberculosis patients and unmatched healthy controls. Am J Clin Nutr. 2007, 86: 1376-1383.PubMedGoogle Scholar
- Deeb KK, Trump DL, Johnson CS: Vitamin D signalling pathways in cancer: potential for anticancer therapeutics. Nat Rev Cancer. 2007, 7: 684-700. 10.1038/nrc2196.PubMedView ArticleGoogle Scholar
- Nejentsev S, Godfrey L, Snook H, Rance H, Nutland S, Walker NM, Lam AC, Guja C, Ionescu-Tirgoviste C, Undlien DE, Ronningen KS, Tuomilehto-Wolf E, Tuomilehto J, Newport MJ, Clayton DG, Todd JA: Comparative high-resolution analysis of linkage disequilibrium and tag single nucleotide polymorphisms between populations in the vitamin D receptor gene. Hum Mol Genet. 2004, 13: 1633-1639. 10.1093/hmg/ddh169.PubMedView ArticleGoogle Scholar
- Zinser G, Packman K, Welsh J: Vitamin D(3) receptor ablation alters mammary gland morphogenesis. Development. 2002, 129: 3067-3076.PubMedGoogle Scholar
- Welsh J: Vitamin D and breast cancer: insights from animal models. Am J Clin Nutr. 2004, 80: 1721S-1724S.PubMedGoogle Scholar
- Welsh J: Vitamin D and prevention of breast cancer. Acta Pharmacol Sin. 2007, 28: 1373-1382. 10.1111/j.1745-7254.2007.00700.x.PubMedView ArticleGoogle Scholar
- International Agency for Research on Cancer: Vitamin D and cancer. IACR Working Group Reports. 2008, Lyon, France: International Agency for Research on Cancer, 5:Google Scholar
- 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: Calcium plus vitamin D supplementation and the risk of breast cancer. J Natl Cancer Inst. 2008, 100: 1581-1591. 10.1093/jnci/djn360.PubMedPubMed CentralView ArticleGoogle Scholar
- Zinser GM, Suckow M, Welsh J: Vitamin D receptor (VDR) ablation alters carcinogen-induced tumorigenesis in mammary gland, epidermis and lymphoid tissues. J Steroid Biochem Mol Biol. 2005, 97: 153-164. 10.1016/j.jsbmb.2005.06.024.PubMedView ArticleGoogle Scholar
- Neuhouser ML, Sorensen B, Hollis BW, Ambs A, Ulrich CM, McTiernan A, Bernstein L, Wayne S, Gilliland F, Baumgartner K, Baumgartner R, Ballard-Barbash R: Vitamin D insufficiency in a multiethnic cohort of breast cancer survivors. Am J Clin Nutr. 2008, 88: 133-139.PubMedPubMed CentralGoogle Scholar
- Abbas S, Chang-Claude J, Linseisen J: Plasma 25-hydroxyvitamin D and premenopausal breast cancer risk in a German case-control study. Int J Cancer. 2009, 124: 250-255. 10.1002/ijc.23904.PubMedView ArticleGoogle Scholar
- Abbas S, Linseisen J, Chang-Claude J: Dietary vitamin D and calcium intake and premenopausal breast cancer risk in a German case-control study. Nutr Cancer. 2007, 59: 54-61. 10.1080/01635580701390223.PubMedView ArticleGoogle Scholar
- Yao S, Sucheston LE, Millen AE, Johnson CS, Trump DL, Nesline MK, Davis W, Hong CC, McCann SE, Hwang H, Kulkarni S, Edge SB, O'Connor TL, Ambrosone CB: Pretreatment serum concentrations of 25-hydroxyvitamin D and breast cancer prognostic characteristics: a case-control and a case-series study. PLoS One. 2011, 6: e17251-10.1371/journal.pone.0017251.PubMedPubMed CentralView ArticleGoogle Scholar
- Ambrosone CB, Ciupak GL, Bandera EV, Jandorf L, Bovbjerg DH, Zirpoli G, Pawlish K, Godbold J, Furberg H, Fatone A, Valdimarsdottir H, Yao S, Li Y, Hwang H, Davis W, Roberts M, Sucheston L, Demissie K, Amend KL, Tartter P, Reilly J, Pace BW, Rohan T, Sparano J, Raptis G, Castaldi M, Estabrook A, Feldman S, Weltz C, Kemeny M: Conducting molecular epidemiological research in the age of HIPAA: a multi-institutional case-control study of breast cancer in African-American and European-American women. J Oncol. 2009, 2009: 871250-PubMedPubMed CentralView ArticleGoogle Scholar
- The International HapMap Consortium: The International HapMap Project. Nature. 2003, 426: 789-796. 10.1038/nature02168.View ArticleGoogle Scholar
- SeattleSNPs Genome Variation Server. [http://gvs.gs.washington.edu/GVS/]
- Xu Z, Kaplan NL, Taylor JA: TAGster: efficient selection of LD tag SNPs in single or multiple populations. Bioinformatics. 2007, 23: 3254-3255. 10.1093/bioinformatics/btm426.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsai HJ, Choudhry S, Naqvi M, Rodriguez-Cintron W, Burchard EG, Ziv E: Comparison of three methods to estimate genetic ancestry and control for stratification in genetic association studies among admixed populations. Hum Genet. 2005, 118: 424-433. 10.1007/s00439-005-0067-z.PubMedView ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMedPubMed CentralGoogle Scholar
- Luna A, Nicodemus KK: snp.plotter: an R-based SNP/haplotype association and linkage disequilibrium plotting package. Bioinformatics. 2007, 23: 774-776. 10.1093/bioinformatics/btl657.PubMedView ArticleGoogle Scholar
- Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science. 2002, 296: 2225-2229. 10.1126/science.1069424.PubMedView ArticleGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang W, Duan S, Kistner EO, Bleibel WK, Huang RS, Clark TA, Chen TX, Schweitzer AC, Blume JE, Cox NJ, Dolan ME: Evaluation of genetic variation contributing to differences in gene expression between populations. Am J Hum Genet. 2008, 82: 631-640. 10.1016/j.ajhg.2007.12.015.PubMedPubMed CentralView ArticleGoogle Scholar
- Signorello LB, Williams SM, Zheng W, Smith JR, Long J, Cai Q, Hargreaves MK, Hollis BW, Blot WJ: Blood vitamin d levels in relation to genetic estimation of African ancestry. Cancer Epidemiol Biomarkers Prev. 2010, 19: 2325-2331. 10.1158/1055-9965.EPI-10-0482.PubMedPubMed CentralView ArticleGoogle Scholar
- Curran JE, Vaughan T, Lea RA, Weinstein SR, Morrison NA, Griffiths LR: Association of A vitamin D receptor polymorphism with sporadic breast cancer development. Int J Cancer. 1999, 83: 723-726. 10.1002/(SICI)1097-0215(19991210)83:6<723::AID-IJC4>3.0.CO;2-3.PubMedView ArticleGoogle Scholar
- Sillanpaa P, Hirvonen A, Kataja V, Eskelinen M, Kosma VM, Uusitupa M, Vainio H, Mitrunen K: Vitamin D receptor gene polymorphism as an important modifier of positive family history related breast cancer risk. Pharmacogenetics. 2004, 14: 239-245. 10.1097/00008571-200404000-00003.PubMedView ArticleGoogle Scholar
- Hou MF, Tien YC, Lin GT, Chen CJ, Liu CS, Lin SY, Huang TJ: Association of vitamin D receptor gene polymorphism with sporadic breast cancer in Taiwanese patients. Breast Cancer Res Treat. 2002, 74: 1-7. 10.1023/A:1016048900049.PubMedView ArticleGoogle Scholar
- Arai H, Miyamoto KI, Yoshida M, Yamamoto H, Taketani Y, Morita K, Kubota M, Yoshida S, Ikeda M, Watabe F, Kanemasa Y, Takeda E: The polymorphism in the caudal-related homeodomain protein Cdx-2 binding element in the human vitamin D receptor gene. J Bone Miner Res. 2001, 16: 1256-1264. 10.1359/jbmr.2001.16.7.1256.PubMedView ArticleGoogle Scholar
- Fang Y, van Meurs JB, Bergink AP, Hofman A, van Duijn CM, van Leeuwen JP, Pols HA, Uitterlinden AG: Cdx-2 polymorphism in the promoter region of the human vitamin D receptor gene determines susceptibility to fracture in the elderly. J Bone Miner Res. 2003, 18: 1632-1641. 10.1359/jbmr.2003.18.9.1632.PubMedView ArticleGoogle Scholar
- Abbas S, Nieters A, Linseisen J, Slanger T, Kropp S, Mutschelknauss EJ, Flesch-Janys D, Chang-Claude J: Vitamin D receptor gene polymorphisms and haplotypes and postmenopausal breast cancer risk. Breast Cancer Res. 2008, 10: R31-10.1186/bcr1994.PubMedPubMed CentralView ArticleGoogle Scholar
- John EM, Schwartz GG, Koo J, Wang W, Ingles SA: Sun exposure, vitamin D receptor gene polymorphisms, and breast cancer risk in a multiethnic population. Am J Epidemiol. 2007, 166: 1409-1419. 10.1093/aje/kwm259.PubMedView ArticleGoogle Scholar
- VandeVord PJ, Wooley PH, Darga LL, Severson RK, Wu B, Nelson DA: Genetic determinants of bone mass do not relate with breast cancer risk in US white and African-American women. Breast Cancer Res Treat. 2006, 100: 103-107. 10.1007/s10549-006-9228-6.PubMedView ArticleGoogle Scholar
- Trabert B, Malone KE, Daling JR, Doody DR, Bernstein L, Ursin G, Marchbanks PA, Strom BL, Humphrey MC, Ostrander EA: Vitamin D receptor polymorphisms and breast cancer risk in a large population-based case-control study of Caucasian and African-American women. Breast Cancer Res. 2007, 9: R84-10.1186/bcr1833.PubMedPubMed CentralView ArticleGoogle Scholar
- Poynter JN, Jacobs ET, Figueiredo JC, Lee WH, Conti DV, Campbell PT, Levine AJ, Limburg P, Le Marchand L, Cotterchio M, Newcomb PA, Potter JD, Jenkins MA, Hopper JL, Duggan DJ, Baron JA, Haile RW: Genetic variation in the vitamin D receptor (VDR) and the vitamin D-binding protein (GC) and risk for colorectal cancer: results from the Colon Cancer Family Registry. Cancer Epidemiol Biomarkers Prev. 2010, 19: 525-536. 10.1158/1055-9965.EPI-09-0662.PubMedPubMed CentralView ArticleGoogle Scholar
- Kahles H, Morahan G, Todd JA, Badenhoop K: Association analyses of the vitamin D receptor gene in 1654 families with type I diabetes. Genes Immun. 2009, 10 (Suppl 1): S60-63.PubMedPubMed CentralView ArticleGoogle Scholar
- Townsend K, Banwell CM, Guy M, Colston KW, Mansi JL, Stewart PM, Campbell MJ, Hewison M: Autocrine metabolism of vitamin D in normal and malignant breast tissue. Clin Cancer Res. 2005, 11: 3579-3586. 10.1158/1078-0432.CCR-04-2359.PubMedView ArticleGoogle Scholar
This article is published under license to 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.