Genetic variation in the vitamin D receptor (VDR) and the vitamin D binding protein (GC) and risk of colorectal cancer: Results from the Colon Cancer Family Registry

doi:10.1158/1055-9965.EPI-09-0662

Cancer Epidemiol Biomarkers Prev. Author manuscript; available in PMC 2011 February 1.

Published in final edited form as:

Cancer Epidemiol Biomarkers Prev. 2010 February; 19(2): 525.

Published online 2010 January 19. doi: 10.1158/1055-9965.EPI-09-0662

PMCID: PMC2819604

NIHMSID: NIHMS160891

Genetic variation in the vitamin D receptor (VDR) and the vitamin D binding protein (GC) and risk of colorectal cancer: Results from the Colon Cancer Family Registry

Jenny N. Poynter,^1,² Elizabeth T. Jacobs,³ Jane C. Figueiredo,¹ Won H. Lee,¹ David V. Conti,¹ Peter T. Campbell,^4,⁵ A. Joan Levine,¹ Paul Limburg,⁶ Loic Le Marchand,⁷ Michelle Cotterchio,⁸ Polly A. Newcomb,⁵ John D. Potter,⁵ Mark A. Jenkins,⁹ John L. Hopper,⁹ David J. Duggan,¹⁰ John A. Baron,¹¹ and Robert W. Haile¹

¹Department of Preventive Medicine, University of Southern California, Los Angeles, CA

²Division of Pediatric Epidemiology and Clinical Research, University of Minnesota, Minneapolis, MN

³Arizona Cancer Center and Mel and Enid Zuckerman College of Public Health, University of Arizona, Tuscon, AZ

⁴Department of Epidemiology, American Cancer Society, Atlanta, GA

⁵Cancer Prevention Program, Fred Hutchinson Cancer Research Center, Seattle, WA

⁶Mayo Clinic, Rochester, MN

⁷Cancer Research Center of Hawaii, University of Hawaii, Honolulu, HI

⁸Cancer Care Ontario, Toronto, ON, Canada

⁹Centre for Molecular, Environmental, Genetic and Analytic Epidemiology, University of Melbourne, Melbourne, Australia

¹⁰Translational Genomics Research Institute (TGen), Phoenix, AZ

¹¹Department of Medicine, Dartmouth Medical School, Lebanon, NH

Requests for Reprints: Jenny N. Poynter, Ph.D., Division of Pediatric Epidemiology and Clinical Research, University of Minnesota, 420 Delaware St SE, Moos 1-117, Minneapolis, MN 55455, Tel: 612-625-4232, Fax: 612-624-7147, Email: poynt006/at/umn.edu

The publisher's final edited version of this article is available free at Cancer Epidemiol Biomarkers Prev

See other articles in PMC that cite the published article.

Abstract

Epidemiologic evidence supports a role for vitamin D in colorectal cancer (CRC) risk. Variants in vitamin D-related genes might modify the association between vitamin D levels and CRC risk. In this analysis, we performed a comprehensive evaluation of common variants in the vitamin D receptor (VDR) and the vitamin D binding protein (GC, group-specific component) genes using a population-based case-unaffected sibling control design that included 1,750 sibships recruited into the Colon Cancer Family Registry (Colon CFR). We also evaluated whether any associations differed by calcium supplement use, family history of CRC, or tumor characteristics. Heterogeneity by calcium and vitamin D intake was evaluated for a subset of 585 cases and 837 sibling controls who completed a detailed food frequency questionnaire (FFQ). Age- and sex-adjusted associations were estimated using conditional logistic regression. Overall, we did not find evidence for an association between any SNP in VDR or GC and risk of CRC (range of unadjusted p-values 0.01—0.98 for VDR and 0.07—0.95 for GC). None of these associations was significant after adjustment for multiple comparisons. We also found no evidence that calcium or vitamin D intake (food and supplement) from the FFQ modified the association estimates between VDR and GC SNPs and CRC. We did observe associations between SNPs in GC and microsatellite unstable CRC, although these results should be confirmed in additional studies. Overall, our results do not provide evidence for a role of common genetic variants in VDR or GC in susceptibility to CRC.

Keywords: vitamin D, vitamin D receptor (VDR), vitamin D binding protein (GC, DBP), colorectal cancer

Introduction

Sunlight and vitamin D were first hypothesized to reduce risk of colorectal cancer in 1980 (1). The epidemiologic evidence has supported this relationship as both cohort (2-5) and case-control studies (6-8) have reported an inverse relationship between vitamin D intake and CRC, although the association did not reach statistical significance in all studies. An inverse association between serum levels of 25-OHD, a biomarker of vitamin D status, and CRC risk has been also demonstrated in epidemiologic studies (5, 9-15). In addition, high plasma 25-OHD was associated with decreased proliferation in colonic tissue from subjects at risk for colorectal neoplasms because they had a personal history of adenomas (16).

The vitamin D receptor (VDR) is a member of the steroid superfamily of nuclear receptors and plays a key role in regulating the transcriptional activity of the vitamin D metabolite, 1,25-(OH)₂D₃ (17). Polymorphisms in VDR have been studied extensively in epidemiologic studies of CRC (18-26). These studies have mostly focused on a few selected variants, including the FokI (rs2228570, previously rs10735810), TaqI (rs731236), BsmI (rs1544410), and ApaI (rs7975232) restriction sites, a polymorphism in the CDX2 binding element in VDR (rs11568820), and a poly-A microsatellite. Overall, the results have been mixed for studies of individual variants in CRC (18-22, 24-28). Several studies have suggested that haplotypes of these selected variants are associated with risk of CRC (22, 23, 29) and that the association with these variants may be modified by dietary factors (18, 20, 21, 30, 31), physical activity (32), and BMI (32).

The gene for the vitamin D binding protein (GC) is also a logical candidate in the vitamin D pathway because its major function is to transport vitamin D metabolites in the blood. Variants in this gene have been shown to alter plasma concentrations of 25-OHD (33, 34). Selected variants in GC have been evaluated in studies of breast (35, 36) and prostate cancers (37); however, to date, no studies have been reported with colorectal neoplasia.

Despite the large number of studies that have evaluated the association between selected genetic variants in VDR and CRC, a comprehensive analysis of VDR and GC has not yet been published. In this analysis, we used data from a multi-site family-based case-control study conducted by the Colon Cancer Family Registry (Colon CFR) to evaluate associations between common single nucleotide polymorphisms (SNPs) in VDR and GC and CRC. In addition, we evaluated heterogeneity by calcium supplement use, dietary vitamin D and calcium intake, family history of CRC, tumor subsite in the colon or rectum, and presence/absence of tumor microsatellite instability.

Materials and Methods

Study Population

Participants were recruited for the Colon CFR from six registry centers: the University of Hawaii (Honolulu), Fred Hutchinson Cancer Research Center (Seattle, WA), Mayo Clinic (Rochester, MN), University of Southern California Consortium (Los Angeles, CA), Cancer Care Ontario (Toronto, Canada), and the University of Melbourne (Victoria, Australia). Families were ascertained through population-based cancer registries (population-based) and high risk clinics (clinic-based). Some centers recruited all incident cases of CRC whereas others over-sampled cases with a family history of CRC or young age at diagnosis of CRC. Standardized procedures were used to collect epidemiologic data and blood samples. Tumor blocks and pathology reports from cases were obtained from the Colon CFR Jeremy Jass Memorial Pathology Bank. Detailed information about the Colon CFR can be found online¹ and is summarized in Newcomb et al. (38).

In this analysis, we used a case-unaffected sibling control design including only population-based families. Cases were defined as probands who were diagnosed with invasive CRC from 1997-2005 and affected relatives of probands with a diagnosis of invasive CRC. All cases were interviewed within 5 years of diagnosis (75% within 2 years). Controls were defined as siblings without a diagnosis of CRC. A total of 4,704 population-based cases (n=1,815) and unaffected sibling controls (n=2,889) were genotyped for this analysis. The majority of discordant sibships had only one case (96.3%), and the remaining sibships included one or more affected relatives. In addition, we also genotyped a random set of unrelated population-based controls (n=447) from one of the Colon CFR sites (i.e., FHCRC).

We obtained informed consent from all participants. The study was approved by the Institutional Review Board(s) at each Colon CFR site.

Data Collection

A core risk-factor questionnaire, administered to all participants at the time of recruitment, was used to collect information on personal and family history of polyps, colorectal cancer, and other cancers, as well as other risk factors, including demographic information, medication use, reproductive history (females only), physical activity, alcohol consumption, tobacco use, and a brief section on diet. In addition, a well-validated food frequency questionnaire (FFQ), developed at the Cancer Research Center of Hawaii (39), was administered to all participants at baseline for three of the Colon CFR sites (USC consortium, Ontario, and Hawaii) representing 32% of the overall study population. The FFQ included questions on both dietary and supplemental intake of calcium and vitamin D. Calcium and vitamin D intake from food and supplements were evaluated as nutrient densities (per 1000 KCAL/day). Individuals were asked to report exposures at a reference time of two years before diagnosis for cases and two years before recruitment for unaffected siblings.

Genotyping

TagSNPs were selected using Haploview Tagger (40) with the following criteria: minor allele frequency (MAF)>5%; pairwise r²>0.95; and distance from closest SNP >60 bp. Linkage disequilibrium (LD) blocks were determined using data from HapMap data release #16c.1, June 2005, on NCBI B34 assembly, dbSNP b124. For each gene, we extended examination of the 5′- and 3′-UTR regions to include the 5′- and 3′-most SNP within the LD block (approximately 10kb upstream and 5kb downstream). In regions with no or low LD, SNPs (MAF>5%) were selected from either HapMap or dbSNP at a density of approximately 1 per kb. Finally, all non-synonymous SNPs and SNPs previously reported in the literature were included, regardless of MAF.

DNA was extracted from blood samples for genotyping (38). Genotyping was performed using the Illumina GoldenGate genotyping assay (41). We implemented a series of quality-control checks based on Illumina metrics, and SNPs were excluded from analysis based on the following standards: GenTrain Score <0.4; 10%GC Score <0.25; AB T Dev >0.1239; call rate < 0.95; more than 2 P-P-C errors. Inter- and intra-plate replicates were included at a rate of ~5% on all plates, and SNPs were excluded from the analysis if there were greater than two errors on replicate genotypes. In addition, genotype data from 30 CEPH trios (Coriell Cell Repository, Camden, NJ) were used to confirm reliability and reproducibility of the genotyping. SNPs were excluded from the analysis if more than 3 discordant genotypes were discovered on comparison with genotypes from the International HapMap Project (42).

We performed additional genotyping using Sequenom's iPLEX Gold for tagSNPs that were not successfully genotyped on the Illumina platform and for additional SNPs selected to ensure adequate coverage based on updated HapMap data (v.21). These additional SNPs were selected using Haploview Snagger (r²>0.95, MAF> 0.05) (43). Polymerase chain reaction (PCR) and extension primers for these SNPs were designed using the MassARRAY Assay Design 3.0 software (Sequenom, Inc) and are available upon request. PCR amplification and single base extension reactions were performed according to the manufacturer's instructions. Extension product sizes were determined by mass spectrometry using Sequenom's Compact MALDI-TOF mass spectrometer. The resulting mass-spectra were converted to genotype data using SpectroTYPER-RT software.

As a QC measure, the frequency of discordant genotypes was estimated: 2 of 281 (0.7%) blinded replicates were discordant; these samples were excluded from all analyses. In addition, 3 samples were excluded because their call rate was less than 0.85 for SNPs genotyped on the Illumina platform.

In this analysis, we report results for two genes in the vitamin D pathway: VDR and GC. We genotyped 50 tagSNPs in VDR and 26 tagSNPs in GC. Seven SNPs in VDR were excluded due to discordant genotypes compared with HapMap data (rs3847987) or poor quality by Illumina QC standards (rs1540339, rs4237855, rs739837, rs7975232, rs7136534, rs8179174). We excluded one SNP in GC from the analysis because of discordant genotypes compared with HapMap data for the 30 CEPH trios (rs6837549).

Microsatellite instability status

Microsatellite instability (MSI) was evaluated for all cases with available tumor tissue using a panel of 10 markers (BAT25, BAT26, BAT40, MYCL, D5S346, D17S250, ACTC, D18S55, D10S197, and BAT34C4) using standard techniques (44). Results were required for at least four markers to determine MSI status. Tumors were deemed MSI-H if instability was observed at ≥30% of markers; MSI-L if >0 and <30% of markers were unstable; and MSS if all markers were stable.

Tumor Location

Tumors were classified by location in the colon using International Classification of Diseases for Oncology, third edition (ICD-O-3) codes (45). Tumors located in the cecum, ascending colon, hepatic flexure, transverse colon, and splenic flexure (ICD-O-3 codes C180, C182, C183, C184, and C185) were classified as proximal colon. Tumors located in the descending colon and sigmoid colon (ICD-O-3 codes C186 and C187) were classified as distal colon. Rectal tumors included those of the rectosigmoid junction and rectum (ICD-O-3 codes C199 and C209).

Statistical analysis

Statistical analyses were conducted in R (Version 2.7.1) (46). Minor allele frequency was estimated using genotype data collected on unrelated population based controls. We determined pairwise linkage disequilibrium between SNPs within a gene using the square of the correlation coefficient between markers (r²). We also evaluated Hardy Weinberg Equilibrium for each SNP.

We estimated associations between variants and risk of CRC using multivariable conditional logistic regression with sibship as the matching factor. All analyses were adjusted for age (continuous) and sex. We evaluated calcium supplement use (yes/no), multivitamin use (yes/no), and physical activity (average weekly MET hours of physical activity throughout adulthood) as potential confounders. Potential confounders that changed the odds ratio (OR) by more than 10% were included in the final model. We tested the associations between tagSNPs and CRC using a log additive model, in which the estimates of effect represent the risk due to one additional variant. Because we do not know which, if any, of these tagSNPs are causal variants, we used a robust variance estimator to prevent biased estimates from testing association in the presence of linkage (47). P-values were adjusted for multiple testing, taking into account correlated tagSNPs using the p_ACT test of Conneely and Boehnke (48) with a modification to allow for inclusion of covariates.

We estimated stratum-specific ORs to evaluate heterogeneity by calcium supplement use, total vitamin D and calcium intake from the FFQ, and family history of CRC in a first-degree relative as reported by the proband. We evaluated differences in the associations by MSI status and by tumor location, by stratifying the matched sets on the tumor characteristics of the case. We assigned the unaffected sibling to the same MSI or tumor-site category as the case and included interaction terms in the conditional logistic regression models to estimate these stratum-specific ORs. A likelihood ratio test was used to assess heterogeneity by comparing models with the interaction terms to a model that included only the main effects of genotype. We used a Bonferroni correction to determine statistical significance of the stratum-specific ORs.

Results

Selected characteristics of the study population are illustrated in Table 1. Nineteen individuals (9 cases and 10 sibling controls) were removed from the analysis because of missing data on age or sex. The cases and unaffected siblings had similar age and sex distributions, and the majority of individuals reported non-Hispanic white ethnicity (Table 1). MSI results were available for 66% of these population-based cases. Cases with missing MSI status were younger and less likely to have a family history of cancer (data not shown). Family history of CRC, defined as having at least one first-degree relative affected with CRC, was reported by 30% of the cases.

Table 1

Selected Characteristics of the Study Population

There was no association between calcium supplement use and CRC (OR=1.07, 95% CI 0.92—1.25, p=0.39 after adjustment for age and sex). In the subset of individuals with FFQ data, total calcium and vitamin D intake were similar in cases and controls (vitamin D: OR=0.86, 95% CI 0.63—1.18 for highest vs. lowest quartile; calcium: OR=0.79, 95% CI 0.56—1.11 for highest vs. lowest quartile).

After adjustment for multiple testing, we observed no statistically significant associations between SNPs in VDR or GC and CRC using log-additive models (Table 2). Models were adjusted for age and sex. Addition of other potential confounders [calcium supplement use (yes/no), multivitamin use (yes/no), and physical activity (average weekly MET hours of physical activity throughout adulthood)] did not appreciably modify the risk estimates (data not shown) and the more parsimonious models are presented. In the three VDR SNPs in our dataset that have previously been evaluated in epidemiologic studies of CRC, we observed no evidence for an association (TaqI: OR=1.04, 95% CI 0.92—1.17; BsmI: OR=1.04, 95% CI 0.92—1.18; CDX2: OR=1.07, 95% CI 0.92—1.25). We repeated the analysis of SNP main effects in the subset of cases with MSI data, and the results did not differ substantially from those shown above (data not shown). We observed no indication that any of the associations with SNPs in VDR or GC were modified by calcium supplement use (Supplementary Table 1).

Table 2

Associations between SNPs in VDR and GC and CRC in population-based sibships

We evaluated associations between VDR and GC variants and risk of CRC after stratification by total vitamin D intake in the subset of sibships from Colon CFR centers that administered the FFQ (N=585 cases and 837 unaffected siblings from 563 sibships). We found an inverse relationship between CRC and VDR (rs2239186) in the stratum with vitamin D intake above the median (OR= 0.72, 95% CI 0.53—0.97), and saw ORs greater than 1 for rs11574143, rs11168267, rs2107301, and rs2239180 in the stratum with vitamin D intake below the median (Table 3). None of these associations was statistically significant after Bonferroni correction for multiple testing. We observed no evidence for modification by total calcium intake measured by the FFQ (supplemental or dietary intake) (data not shown).

Table 3

Associations between VDR and GC variants and CRC by total vitamin D intake

There was no evidence for heterogeneity by tumor location or family history of CRC for either VDR or GC (data not shown). There was also no evidence for heterogeneity by MSI status for SNPs in VDR (data not shown). However, we observed inverse associations that were statistically significant at the 0.05 level for 15 SNPs in GC (rs13117483, rs1491709, rs222014, rs222016, rs222017, rs16847015, rs1352843, rs222029, rs1352844, rs3733359, rs6817912, rs16847039, rs705125, rs16847047, rs843007) and risk of MSI-H tumors (Table 4). These are unlikely to be independent associations as each of these SNPs was highly correlated (r² > 0.80) with at least one other SNP within this list.

Table 4

Association between SNPs in GC and CRC by MSI status

Discussion

In this family-based case-control study of CRC, we did not find evidence to suggest that variation in VDR and GC is associated with risk of CRC overall. Our data also suggested that associations between SNPs in these genes are not modified by calcium supplement use, family history of CRC, or total calcium intake estimated from an FFQ. We found only weak suggestions that dietary vitamin D could modify the association between SNPs in VDR and risk of CRC (Table 3), although these associations were not statistically significant after correction for multiple testing. We observed statistically significant associations between several SNPs in GC and MSI-H CRC; however, these results should be interpreted with caution due to the many comparisons we have made.

There are several mechanisms that may explain the inhibitory action of vitamin D on carcinogenesis, including reduced cellular proliferation (16, 49), induction of apoptosis (49-53), inhibition of inflammation (54) and angiogenesis (55), and protection of colonic epithelium from the damaging effects of bile acids (56). Given that VDR is expressed in a wide variety of tissues, including some with no known role in calcium metabolism (17), it is logical to hypothesize that variation in VDR may influence these potential anti-carcinogenic effects of vitamin D through altered transcriptional control of 1,25-(OH)₂D₃. In addition, two SNPs in GC (rs4588 and rs7041) have been shown to influence 25-OHD levels (33, 34), which suggests that variation in this gene could also influence the actions of vitamin D. We would also expect that the effect of variants in these genes may be most relevant in the context of vitamin D levels.

Whereas the majority of epidemiologic studies of VDR variants and risk of CRC have found associations with the disease, the results for specific variants have not been entirely consistent. Inconsistency has been a common problem in genetic epidemiology studies (57, 58) and is not limited to studies of the VDR gene. Potential explanations for heterogeneity in genetic association studies include inadequate power, bias, and population differences.

The FokI restriction site polymorphism (rs2228570) is the most frequently reported VDR variant in epidemiologic studies of CRC (18-20, 22, 24-26) and adenoma (31, 59, 60). In studies of CRC, three of the seven studies that evaluated this variant found evidence for an association with risk. In contrast, this variant was not associated with adenoma risk in three studies (31, 59, 60), although a statistically significant association for risk of large adenomas was observed in one of these (31). In a recent meta-analysis, the summary OR for this variant was not associated with CRC, although there was significant heterogeneity among studies (27). Unfortunately, we were unable to evaluate associations between the FokI variant and CRC risk in our study as this SNP did not pass quality-control checks.

In three studies, the BsmI variant (rs1544410) was not associated with risk of colorectal adenoma (30, 31, 61); however, data from two of the studies suggested that the variant may modify the relationship between calcium and/or vitamin D intake and risk of adenoma (30, 61). In studies of CRC, this variant was associated with risk of CRC overall in three studies (19, 21, 26) with evidence for modification by dietary intake of calcium and vitamin D (21). However, two additional studies showed no association with CRC (18, 25). In a recent meta-analysis, the summary OR was borderline protective for this variant after one study contributing to heterogeneity in risk estimates was excluded (27). In our analysis, we observed no overall association between the BsmI variant and CRC (Table 2). In addition, we observed no evidence of modification by calcium supplement use or by dietary calcium and vitamin D intake in the subset of individuals with FFQ data.

The TaqI variant (rs731236) was not associated with risk of colorectal adenoma or cancer in three studies (25, 59, 60) but was associated with a borderline statistically significant risk of CRC in one (19). We observed no overall association between this variant and risk of CRC (Table 2). Three previous studies have evaluated associations between the VDR CDX2 variant (rs11568820) and CRC risk, with two studies reporting no association between this variant and risk (18, 23) while the other study reported a statistically significant increased risk in individuals with the AA genotype (24). We observed no association between this SNP and CRC in our study population (Table 2).

To our knowledge, this is the first epidemiologic study to evaluate the association between variation in GC and risk of CRC. Our data suggest that this gene is unlikely to play a major role in genetic susceptibility to this malignancy. We observed limited evidence that variants in GC may be associated with a reduced risk of MSI-H CRC, although these results should be interpreted with caution due to the many comparisons that we have made as well as the lack (at least to date) of a biologic rationale for a role of this gene only in MSI-H CRC.

This study has several strengths including the comprehensive evaluation of both VDR and GC, the availability of systematically collected data on tumor characteristics, and the case-unaffected sibling study design that controls for any potential confounding by ethnicity. This family-based design is more powerful than studies that include unrelated controls for detecting gene-environment interactions (62).

Several limitations must also be considered. While the case-unaffected sibling design is more powerful for detecting gene-environment interactions, this design may have lower power for detecting genetic main effects because the study sample is often overmatched on genotype (62). Because overmatching is an issue of power, rather than bias, evaluating the width of the confidence intervals can provide information regarding the overall power of the study. For the main effects of the genotypes (Table 2), the confidence intervals are quite narrow, especially for the more common variants. This would indicate that the lack of association in our study is not simply an issue of inadequate power. Detailed FFQ data were collected only for a subset of the participants in this study (32%), so we have limited power to detect heterogeneity by calcium and vitamin D intake. In addition, we did not collect information on sun exposure, so we are unable to account for endogenously synthesized vitamin D levels. Cases were interviewed up to five years after diagnosis, which introduces the potential for possible survival bias and increases the potential for exposure misclassification; however, when we restricted our analyses to subjects interviewed within the first two years after diagnosis, we observed very similar results. In addition, we have no strong a priori hypothesis that any of these SNPs would be strongly associated with prognosis.

In conclusion, we did not observe any statistically significant associations between variants in VDR and GC and CRC overall in this family based analysis. After correction for multiple comparisons, we observed no indication that these associations were modified by calcium supplement use and dietary and total intake of calcium and vitamin D. We did observe associations between SNPs in GC and MSI-H CRC, although these results should be confirmed in additional studies.

Supplementary Material

Click here to view.^{(162K, doc)}

Acknowledgments

The authors would like to thank the following individuals who helped prepare the dataset for these analyses: Margreet Luchtenborg, Maj Earle, Barbara Saltzman, Kathleen Kennedy, Darin Taverna, Chris Edlund, Matt Westlake, Paul Mosquin, Darshana Daftary, Douglas Snazel, Allyson Templeton, Terry Teitsch, Helen Chen and Maggie Angelakos. We would also like to thank the participants in the Colon CFR who have generously donated their time for this project.

Financial Support: This work was supported by the National Cancer Institute, National Institutes of Health under RFA # CA-95-011 and through cooperative agreements with the Australasian Colorectal Cancer Family Registry (U01 CA097735), the USC Familial Colorectal Neoplasia Collaborative Group (U01 CA074799), the Mayo Clinic Cooperative Family Registry for Colon Cancer Studies (U01 CA074800), the Ontario Registry for Studies of Familial Colorectal Cancer (U01 CA074783), the Seattle Colorectal Cancer Family Registry (U01 CA074794), and the University of Hawaii Colorectal Cancer Family Registry (U01 CA074806) as well as NCI T32 CA009142 (JNP), NCI R01 CA112237 (RWH), and NCI 1K07CA10629-01A1 (ETJ). P.T.C. and J.C.F. were supported in part by National Cancer Institute of Canada post-PhD Fellowships (#18735 and #17602).

Footnotes

Disclosure: Dr. Limburg has a consulting agreement with Genomic Health, Inc.

¹http://epi.grants.cancer.gov/CFR/

References

1. Garland CF, Garland FC. Do sunlight and vitamin D reduce the likelihood of colon cancer? Int J Epidemiol. 1980;9:227–31. [PubMed]

2. Bostick RM, Potter JD, Sellers TA, McKenzie DR, Kushi LH, Folsom AR. Relation of calcium, vitamin D, and dairy food intake to incidence of colon cancer among older women. The Iowa Women's Health Study Am J Epidemiol. 1993;137:1302–17. [PubMed]

3. Martinez ME, Marshall JR, Sampliner R, Wilkinson J, Alberts DS. Calcium, vitamin D, and risk of adenoma recurrence (United States) Cancer Causes Control. 2002;13:213–20. [PubMed]

4. Martinez ME, Giovannucci EL, Colditz GA, et al. Calcium, vitamin D, and the occurrence of colorectal cancer among women. J Natl Cancer Inst. 1996;88:1375–82. [PubMed]

5. Giovannucci E, Liu Y, Rimm EB, et al. Prospective study of predictors of vitamin D status and cancer incidence and mortality in men. J Natl Cancer Inst. 2006;98:451–9. [PubMed]

6. Pritchard RS, Baron JA, Gerhardsson de Verdier M. Dietary calcium, vitamin D, and the risk of colorectal cancer in Stockholm, Sweden. Cancer Epidemiol Biomarkers Prev. 1996;5:897–900. [PubMed]

7. Marcus PM, Newcomb PA. The association of calcium and vitamin D, and colon and rectal cancer in Wisconsin women. Int J Epidemiol. 1998;27:788–93. [PubMed]

8. Kampman E, Slattery ML, Caan B, Potter JD. Calcium, vitamin D, sunshine exposure, dairy products and colon cancer risk (United States) Cancer Causes Control. 2000;11:459–66. [PubMed]

9. Garland CF, Comstock GW, Garland FC, Helsing KJ, Shaw EK, Gorham ED. Serum 25-hydroxyvitamin D and colon cancer: eight-year prospective study. Lancet. 1989;2:1176–8. [PubMed]

10. Braun MM, Helzlsouer KJ, Hollis BW, Comstock GW. Colon cancer and serum vitamin D metabolite levels 10-17 years prior to diagnosis. Am J Epidemiol. 1995;142:608–11. [PubMed]

11. Tangrea J, Helzlsouer K, Pietinen P, et al. Serum levels of vitamin D metabolites and the subsequent risk of colon and rectal cancer in Finnish men. Cancer Causes Control. 1997;8:615–25. [PubMed]

12. Feskanich D, Ma J, Fuchs CS, et al. Plasma vitamin D metabolites and risk of colorectal cancer in women. Cancer Epidemiol Biomarkers Prev. 2004;13:1502–8. [PubMed]

13. Wactawski-Wende J, Kotchen JM, Anderson GL, et al. Calcium plus vitamin D supplementation and the risk of colorectal cancer. N Engl J Med. 2006;354:684–96. [PubMed]

14. 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–9. [PubMed]

15. Otani T, Iwasaki M, Sasazuki S, Inoue M, Tsugane S. Plasma vitamin D and risk of colorectal cancer: the Japan Public Health Center-Based Prospective Study. Br J Cancer. 2007;97:446–51. [PMC free article] [PubMed]

16. Holt PR, Arber N, Halmos B, et al. Colonic epithelial cell proliferation decreases with increasing levels of serum 25-hydroxy vitamin D. Cancer Epidemiol Biomarkers Prev. 2002;11:113–9. [PubMed]

17. Haussler MR, Whitfield GK, Haussler CA, et al. The nuclear vitamin D receptor: biological and molecular regulatory properties revealed. J Bone Miner Res. 1998;13:325–49. [PubMed]

18. Theodoratou E, Farrington SM, Tenesa A, et al. Modification of the inverse association between dietary vitamin D intake and colorectal cancer risk by a FokI variant supports a chemoprotective action of Vitamin D intake mediated through VDR binding. Int J Cancer. 2008;123:2170–9. [PubMed]

19. Slattery ML, Yakumo K, Hoffman M, Neuhausen S. Variants of the VDR gene and risk of colon cancer (United States) Cancer Causes Control. 2001;12:359–64. [PubMed]

20. Wong HL, Seow A, Arakawa K, Lee HP, Yu MC, Ingles SA. Vitamin D receptor start codon polymorphism and colorectal cancer risk: effect modification by dietary calcium and fat in Singapore Chinese. Carcinogenesis. 2003;24:1091–5. [PubMed]

21. Slattery ML, Neuhausen SL, Hoffman M, et al. Dietary calcium, vitamin D, VDR genotypes and colorectal cancer. Int J Cancer. 2004;111:750–6. [PubMed]

22. Park K, Woo M, Nam J, Kim JC. Start codon polymorphisms in the vitamin D receptor and colorectal cancer risk. Cancer Lett. 2006;237:199–206. [PubMed]

23. Slattery ML, Herrick J, Wolff RK, Caan BJ, Potter JD, Sweeney C. CDX2 VDR polymorphism and colorectal cancer. Cancer Epidemiol Biomarkers Prev. 2007;16:2752–5. [PMC free article] [PubMed]

24. Ochs-Balcom HM, Cicek MS, Thompson CL, et al. Association of vitamin D receptor gene variants, adiposity and colon cancer. Carcinogenesis. 2008;29:1788–93. [PMC free article] [PubMed]

25. Flugge J, Krusekopf S, Goldammer M, et al. Vitamin D receptor haplotypes protect against development of colorectal cancer. Eur J Clin Pharmacol. 2007;63:997–1005. [PubMed]

26. Jenab M, McKay J, Bueno-de-Mesquita HB, et al. Vitamin D receptor and calcium sensing receptor polymorphisms and the risk of colorectal cancer in European populations. Cancer Epidemiol Biomarkers Prev. 2009;18:2485–91. [PubMed]

27. Raimondi S, Johansson H, Maisonneuve P, Gandini S. Review and meta-analysis on vitamin D receptor polymorphisms and cancer risk. Carcinogenesis. 2009;30:1170–80. [PubMed]

28. McCullough ML, Bostick RM, Mayo TL. Vitamin D gene pathway polymorphisms and risk of colorectal, breast, and prostate cancer. Annu Rev Nutr. 2009;29:111–32. [PubMed]

29. Sweeney C, Curtin K, Murtaugh MA, Caan BJ, Potter JD, Slattery ML. Haplotype analysis of common vitamin D receptor variants and colon and rectal cancers. Cancer Epidemiol Biomarkers Prev. 2006;15:744–9. [PubMed]

30. Kim HS, Newcomb PA, Ulrich CM, et al. Vitamin D receptor polymorphism and the risk of colorectal adenomas: evidence of interaction with dietary vitamin D and calcium. Cancer Epidemiol Biomarkers Prev. 2001;10:869–74. [PubMed]

31. Ingles SA, Wang J, Coetzee GA, Lee ER, Frankl HD, Haile RW. Vitamin D receptor polymorphisms and risk of colorectal adenomas (United States) Cancer Causes Control. 2001;12:607–14. [PubMed]

32. Slattery ML, Murtaugh M, Caan B, Ma KN, Wolff R, Samowitz W. Associations between BMI, energy intake, energy expenditure, VDR genotype and colon and rectal cancers (United States) Cancer Causes Control. 2004;15:863–72. [PubMed]

33. Sinotte M, Diorio C, Berube S, Pollak M, Brisson J. Genetic polymorphisms of the vitamin D binding protein and plasma concentrations of 25-hydroxyvitamin D in premenopausal women. Am J Clin Nutr. 2009;89:634–40. [PubMed]

34. Engelman CD, Fingerlin TE, Langefeld CD, et al. Genetic and environmental determinants of 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D levels in Hispanic and African Americans. J Clin Endocrinol Metab. 2008;93:3381–8. [PubMed]

35. McCullough ML, Stevens VL, Diver WR, et al. Vitamin D pathway gene polymorphisms, diet, and risk of postmenopausal breast cancer: a nested case-control study. Breast Cancer Res. 2007;9:R9. [PMC free article] [PubMed]

36. Abbas S, Linseisen J, Slanger T, et al. The Gc2 allele of the vitamin D binding protein is associated with a decreased postmenopausal breast cancer risk, independent of the vitamin D status. Cancer Epidemiol Biomarkers Prev. 2008;17:1339–43. [PubMed]

37. Kidd LC, Paltoo DN, Wang S, et al. Sequence variation within the 5′ regulatory regions of the vitamin D binding protein and receptor genes and prostate cancer risk. Prostate. 2005;64:272–82. [PubMed]

38. Newcomb PA, Baron J, Cotterchio M, et al. Colon cancer family registry: an international resource for studies of the genetic epidemiology of colon cancer. Cancer Epidemiol Biomarkers Prev. 2007;16:2331–43. [PubMed]

39. Stram DO, Hankin JH, Wilkens LR, et al. Calibration of the dietary questionnaire for a multiethnic cohort in Hawaii and Los Angeles. Am J Epidemiol. 2000;151:358–70. [PubMed]

40. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–5. [PubMed]

41. Shen R, Fan JB, Campbell D, et al. High-throughput SNP genotyping on universal bead arrays. Mutat Res. 2005;573:70–82. [PubMed]

42. A haplotype map of the human genome. Nature. 2005;437:1299–320. [PMC free article] [PubMed]

43. Edlund CK, Lee WH, Li D, Van Den Berg DJ, Conti DV. Snagger: a user-friendly program for incorporating additional information for tagSNP selection. BMC Bioinformatics. 2008;9:174. [PMC free article] [PubMed]

44. Lindor NM, Rabe K, Petersen GM, et al. Lower cancer incidence in Amsterdam-I criteria families without mismatch repair deficiency: familial colorectal cancer type X. JAMA. 2005;293:1979–85. [PMC free article] [PubMed]

45. Organization WH. International Classification of Diseases for Oncology. 3rd. Geneva: WHO; 2000.

46. R Core Development Team. R: A language and environment for statistical computing.: R Foundation for Statistical Computing. 2008. Available from: < http://www.r-project.org/>.

47. Siegmund KD, Langholz B, Kraft P, Thomas DC. Testing linkage disequilibrium in sibships. Am J Hum Genet. 2000;67:244–8. [PubMed]

48. Conneely KN, Boehnke M. So Many Correlated Tests, So Little Time! Rapid Adjustment of P Values for Multiple Correlated Tests. Am J Hum Genet. 2007;81 [PubMed]

49. Lamprecht SA, Lipkin M. Chemoprevention of colon cancer by calcium, vitamin D and folate: molecular mechanisms. Nat Rev Cancer. 2003;3:601–14. [PubMed]

50. Palmer HG, Gonzalez-Sancho JM, Espada J, et al. Vitamin D(3) promotes the differentiation of colon carcinoma cells by the induction of E-cadherin and the inhibition of beta-catenin signaling. J Cell Biol. 2001;154:369–87. [PMC free article] [PubMed]

51. Miller EA, Keku TO, Satia JA, Martin CF, Galanko JA, Sandler RS. Calcium, vitamin D, and apoptosis in the rectal epithelium. Cancer Epidemiol Biomarkers Prev. 2005;14:525–8. [PubMed]

52. Diaz GD, Paraskeva C, Thomas MG, Binderup L, Hague A. Apoptosis is induced by the active metabolite of vitamin D3 and its analogue EB1089 in colorectal adenoma and carcinoma cells: possible implications for prevention and therapy. Cancer Res. 2000;60:2304–12. [PubMed]

53. Fedirko V, Bostick RM, Flanders WD, et al. Effects of vitamin D and calcium supplementation on markers of apoptosis in normal colon mucosa: a randomized, double-blind, placebo-controlled clinical trial. Cancer Prev Res (Phila Pa) 2009;2:213–23. [PMC free article] [PubMed]

54. Zhu Y, Mahon BD, Froicu M, Cantorna MT. Calcium and 1 alpha,25-dihydroxyvitamin D3 target the TNF-alpha pathway to suppress experimental inflammatory bowel disease. Eur J Immunol. 2005;35:217–24. [PubMed]

55. Iseki K, Tatsuta M, Uehara H, et al. Inhibition of angiogenesis as a mechanism for inhibition by 1alpha-hydroxyvitamin D3 and 1,25-dihydroxyvitamin D3 of colon carcinogenesis induced by azoxymethane in Wistar rats. Int J Cancer. 1999;81:730–3. [PubMed]

56. Makishima M, Lu TT, Xie W, et al. Vitamin D receptor as an intestinal bile acid sensor. Science. 2002;296:1313–6. [PubMed]

57. Ioannidis JP, Ntzani EE, Trikalinos TA, Contopoulos-Ioannidis DG. Replication validity of genetic association studies. Nat Genet. 2001;29:306–9. [PubMed]

58. Wacholder S, Chanock S, Garcia-Closas M, El Ghormli L, Rothman N. Assessing the probability that a positive report is false: an approach for molecular epidemiology studies. J Natl Cancer Inst. 2004;96:434–42. [PubMed]

59. Peters U, McGlynn KA, Chatterjee N, et al. Vitamin D, calcium, and vitamin D receptor polymorphism in colorectal adenomas. Cancer Epidemiol Biomarkers Prev. 2001;10:1267–74. [PubMed]

60. Grau MV, Baron JA, Sandler RS, et al. Vitamin D, calcium supplementation, and colorectal adenomas: results of a randomized trial. J Natl Cancer Inst. 2003;95:1765–71. [PubMed]

61. Boyapati SM, Bostick RM, McGlynn KA, et al. Calcium, vitamin D, and risk for colorectal adenoma: dependency on vitamin D receptor BsmI polymorphism and nonsteroidal anti-inflammatory drug use? Cancer Epidemiol Biomarkers Prev. 2003;12:631–7. [PubMed]

62. Witte JS, Gauderman WJ, Thomas DC. Asymptotic bias and efficiency in case-control studies of candidate genes and gene-environment interactions: basic family designs. Am J Epidemiol. 1999;149:693–705. [PubMed]