Abstract
Background: Epidemiologic evidence indicates that greater intakes of vitamin D may decrease the risk of colorectal cancer. Variants in the vitamin D receptor (VDR) gene have the potential to modify associations between vitamin D intake and colorectal cancer.
Methods: Associations between intakes of vitamin D and colorectal cancer were studied in a large case–control study conducted in central and northeastern Pennsylvania including 1,012 cases with histologically confirmed colorectal cancer and 1,080 population-based controls. Associations between 35 tagSNPs encompassing the VDR gene and risk for colorectal cancer as well as gene–diet associations were also assessed among a subset of the population (770 controls, 710 cases).
Results: No significant trends were observed between vitamin D intake and colorectal cancer risk. After adjustment for multiple comparisons, none of the SNPs or haplotypes within the VDR gene were associated with colorectal cancer. There were also no interactions between dietary factors and variants in the entire VDR gene.
Conclusions: Overall, results from this study suggest that vitamin D intake and variants in the VDR gene have little effect on risk for colorectal cancer.
Impact: Increasing vitamin D intake from the diet may not result in decreasing the incidence of colorectal cancer. Cancer Epidemiol Biomarkers Prev; 24(10); 1635–7. ©2015 AACR.
Introduction
Cancers of the colon and rectum (colorectal cancer) are currently ranked third for both cancer incidence and mortality in the United States. Previous epidemiologic studies have suggested that greater vitamin D intake is associated with reduced colorectal cancer risk, but these results are inconclusive. One possible explanation for the observed inconsistencies between studies is genetic variation within the vitamin D receptor gene (VDR). Individuals with decreased VDR activity could potentially have an increased risk for colorectal cancer compared with those with normal VDR activity, particularly at low vitamin D intakes.
The objectives of the present study were to comprehensively examine associations between colorectal cancer incidence and (1), consumption of supplemental, dietary, and total vitamin D (2), SNPs that tag the common variation in the gene (tagSNPs) in the VDR gene, and (3) the potential interactions between vitamin D intake and variants in the VDR gene, in a population-based case–control study in a high-risk region of central and northeastern Pennsylvania.
Materials and Methods
Detailed descriptions of this study have been previously described (1–3). Briefly, cases were identified from the Pennsylvania State Cancer Registry within 15 months of diagnosis. To be eligible, case participants had to have a first time, histologically confirmed diagnosis of colon, rectal, or colorectal cancer and be English-speaking. Controls were residents of the same region with no history of colorectal cancer and were identified by random digit dialing.
Participants completed a modified version of the Diet History Questionnaire (DHQ), a validated food frequency questionnaire (FFQ) developed by the NCI (4). This DHQ was modified to examine more comprehensively and accurately the distinct dietary intake habits of subjects recruited from the population of central and northeast Pennsylvania as described previously (3).
A total of 35 tagSNPs were selected for the VDR gene by analyzing individuals of northern and western European ancestry collected by the Centre d'Etude du Polymorphisme Humain (CEPH) from HapMap.org according to the following criteria: (i) SNPs were located in the VDR gene or within its 5-kb 5′ or 3′ flanking region, (b) had a minor allele frequency ≥ 0.05, and (c) the other unselected SNPs could be captured by one of the tagging SNPs with a linkage disequilibrium r2 ≥ 0.80 (mean r2 = 0.98). Correction for multiple comparisons was performed using the false-discovery rate (FDR) (5). All analyses were two-sided and considered significant if P < 0.05 for all tests. Statistical analyses were performed using SAS statistical software (version 9.3; SAS Institute). The Institutional Review Boards at the participating institutes, the Northeast Regional Cancer Institute (Scranton, PA) and the Penn State College of Medicine (Hershey, PA), reviewed and approved all study procedures.
Results
Characteristics of the study participants have been described previously (1–3). Overall, controls were younger, better educated, and consumed greater amounts of supplemental and total vitamin D, calcium, and fiber. Controls also had lower body mass index (BMI), were more likely to regularly use NSAIDs, consumed more alcohol, and were less likely to have a family history of colorectal cancer. In multivariate analyses, there was no clear pattern of association between vitamin D intake and risk for colorectal cancer (Table 1). In the age- and sex-adjusted model, subjects consuming >10 μg/d of supplemental vitamin D were inversely associated with colorectal cancer [OR, 0.61; 95% confidence interval (CI), 0.48–078; FDR Ptrend = 0.005]. Similarly, those with the highest intakes of total vitamin D (>15.5 μg/d) were inversely associated with colorectal cancer (OR, 0.68; 95% CI, 0.52–0.88; FDR Ptrend < 0.001); neither of these remained significant in the fully-adjusted model.
. | Range (μg/day) . | Controls . | Cases . | ORa (95% CI) . | FDR P value . | ORb (95% CI) . | FDR P value . |
---|---|---|---|---|---|---|---|
Dietary vitamin D | |||||||
Quartile 1 | ≤2.65 | 270 | 277 | 1.00 | 1.00 | ||
Quartile 2 | 2.66–3.44 | 270 | 215 | 0.71 (0.56–0.92) | 0.018 | 0.76 (0.58–1.01) | 0.079 |
Quartile 3 | 3.45–4.64 | 270 | 257 | 0.81 (0.63–1.04) | 0.139 | 0.88 (0.67–1.15) | 0.384 |
Quartile 4 | >4.65 | 270 | 263 | 0.87 (0.68–1.11) | 0.301 | 1.04 (0.78–1.38) | 0.791 |
P-trend | 0.532 | 0.573 | 0.432 | 0.465 | |||
Supplemental vitamin D | |||||||
No use | 0 (none) | 345 | 349 | 1.00 | 1.00 | ||
Any vs. none | >0.1 | 735 | 663 | 0.85 (0.71–1.03) | 0.088 | 1.14 (0.93–1.41) | 0.226 |
Low vs. none | ≤10 | 450 | 455 | 0.98 (0.80–1.20) | 0.891 | 1.18 (0.95–1.47) | 0.164 |
High vs. none | >10 | 285 | 208 | 0.61 (0.48–0.78) | <0.001 | 0.99 (0.73–1.35) | 0.967 |
P-trend | 0.003 | 0.005 | 0.433 | 0.467 | |||
Total vitamin D | |||||||
Quartile 1 | ≤4.12 | 270 | 260 | 1.00 | 1.00 | ||
Quartile 2 | 4.12–11.27 | 270 | 292 | 1.04 (0.82–1.33) | 0.735 | 1.38 (1.05–1.80) | 0.026 |
Quartile 3 | 11.28–15.48 | 270 | 242 | 0.82 (0.64–1.06) | 0.189 | 1.05 (0.80–1.38) | 0.725 |
Quartile 4 | >15.48 | 270 | 218 | 0.68 (0.52–0.88) | 0.007 | 1.29 (0.90–1.84) | 0.185 |
P-trend | <0.001 | <0.001 | 0.847 | 0.847 |
. | Range (μg/day) . | Controls . | Cases . | ORa (95% CI) . | FDR P value . | ORb (95% CI) . | FDR P value . |
---|---|---|---|---|---|---|---|
Dietary vitamin D | |||||||
Quartile 1 | ≤2.65 | 270 | 277 | 1.00 | 1.00 | ||
Quartile 2 | 2.66–3.44 | 270 | 215 | 0.71 (0.56–0.92) | 0.018 | 0.76 (0.58–1.01) | 0.079 |
Quartile 3 | 3.45–4.64 | 270 | 257 | 0.81 (0.63–1.04) | 0.139 | 0.88 (0.67–1.15) | 0.384 |
Quartile 4 | >4.65 | 270 | 263 | 0.87 (0.68–1.11) | 0.301 | 1.04 (0.78–1.38) | 0.791 |
P-trend | 0.532 | 0.573 | 0.432 | 0.465 | |||
Supplemental vitamin D | |||||||
No use | 0 (none) | 345 | 349 | 1.00 | 1.00 | ||
Any vs. none | >0.1 | 735 | 663 | 0.85 (0.71–1.03) | 0.088 | 1.14 (0.93–1.41) | 0.226 |
Low vs. none | ≤10 | 450 | 455 | 0.98 (0.80–1.20) | 0.891 | 1.18 (0.95–1.47) | 0.164 |
High vs. none | >10 | 285 | 208 | 0.61 (0.48–0.78) | <0.001 | 0.99 (0.73–1.35) | 0.967 |
P-trend | 0.003 | 0.005 | 0.433 | 0.467 | |||
Total vitamin D | |||||||
Quartile 1 | ≤4.12 | 270 | 260 | 1.00 | 1.00 | ||
Quartile 2 | 4.12–11.27 | 270 | 292 | 1.04 (0.82–1.33) | 0.735 | 1.38 (1.05–1.80) | 0.026 |
Quartile 3 | 11.28–15.48 | 270 | 242 | 0.82 (0.64–1.06) | 0.189 | 1.05 (0.80–1.38) | 0.725 |
Quartile 4 | >15.48 | 270 | 218 | 0.68 (0.52–0.88) | 0.007 | 1.29 (0.90–1.84) | 0.185 |
P-trend | <0.001 | <0.001 | 0.847 | 0.847 |
aAdjusted for age and sex.
bAdjusted for age, sex, BMI, family history of colorectal cancer, smoking, education, NSAID use, and intakes of total energy, alcohol, fiber, and dietary calcium.
Because of insufficient genomic DNA obtained from some subjects, there were differences in the sample numbers for the epidemiologic and genetic datasets. There were no major differences between the two datasets. There was no effect of race in the dataset as 98% of the epidemiologic dataset was Caucasian, and only Caucasians were used in the genotyped dataset. Of the 35 tagSNPs evaluated, only one was marginally associated with colorectal cancer; however, none remained statistically significant after FDR correction (Table 2). There were no significant gene–diet interactions, even at the P = 0.1 cutoff. No haplotypes were significant after correction for multiple comparisons.
SNP . | Cases . | Controls . | Allele odds ratio . | 95% Lower CL . | 95% Upper CL . | P value genotype test . | P value for allele test . | P value for trend test . |
---|---|---|---|---|---|---|---|---|
rs3782905 | 705 | 753 | 0.97 | 0.83 | 1.14 | 0.841 | 0.708 | 0.713 |
rs2283342 | 701 | 749 | 0.85 | 0.69 | 1.03 | 0.131 | 0.103 | 0.099 |
rs2239182 | 692 | 733 | 0.93 | 0.80 | 1.08 | 0.571 | 0.345 | 0.342 |
rs11574077 | 709 | 767 | 0.71 | 0.51 | 1.00 | 0.159 | 0.052 | 0.056 |
rs11574027 | 710 | 768 | 1.50 | 0.83 | 2.71 | 0.353 | 0.178 | 0.193 |
rs2238136 | 699 | 732 | 1.22 | 1.03 | 1.44 | 0.054 | 0.021 | 0.020 |
rs11574046 | 687 | 706 | 0.96 | 0.74 | 1.25 | 0.795 | 0.761 | 0.767 |
rs2107301 | 696 | 739 | 1.07 | 0.91 | 1.26 | 0.622 | 0.404 | 0.402 |
rs11168287 | 692 | 739 | 0.96 | 0.83 | 1.12 | 0.596 | 0.633 | 0.641 |
rs3819545 | 697 | 745 | 0.96 | 0.83 | 1.12 | 0.777 | 0.604 | 0.609 |
rs4237855 | 646 | 700 | 0.91 | 0.78 | 1.05 | 0.541 | 0.198 | 0.284 |
rs3847987 | 706 | 750 | 1.01 | 0.81 | 1.26 | 0.809 | 0.912 | 0.912 |
rs2239186 | 701 | 754 | 0.89 | 0.75 | 1.07 | 0.376 | 0.223 | 0.212 |
rs11168275 | 701 | 741 | 1.02 | 0.86 | 1.21 | 0.965 | 0.807 | 0.809 |
rs11574026 | 706 | 762 | 0.85 | 0.67 | 1.06 | 0.353 | 0.150 | 0.149 |
rs2853564 | 695 | 726 | 0.94 | 0.81 | 1.09 | 0.579 | 0.418 | 0.422 |
rs10875695 | 698 | 749 | 1.07 | 0.90 | 1.26 | 0.718 | 0.449 | 0.449 |
rs11568820 | 708 | 758 | 1.05 | 0.88 | 1.26 | 0.721 | 0.587 | 0.601 |
rs7299460 | 701 | 746 | 0.99 | 0.84 | 1.16 | 0.573 | 0.876 | 0.879 |
rs11168267 | 679 | 695 | 1.02 | 0.81 | 1.29 | 0.901 | 0.848 | 0.846 |
rs2239181 | 700 | 757 | 1.01 | 0.80 | 1.28 | 0.619 | 0.941 | 0.940 |
rs12721364 | 697 | 729 | 1.00 | 0.81 | 1.23 | 0.744 | 0.992 | 0.992 |
rs4516035 | 698 | 733 | 0.98 | 0.85 | 1.14 | 0.620 | 0.790 | 0.799 |
rs10875693 | 698 | 735 | 0.95 | 0.81 | 1.11 | 0.418 | 0.503 | 0.511 |
rs2189480 | 697 | 751 | 1.02 | 0.88 | 1.18 | 0.071 | 0.814 | 0.815 |
rs2239179 | 688 | 710 | 0.91 | 0.78 | 1.05 | 0.297 | 0.202 | 0.182 |
rs4760655 | 703 | 750 | 0.93 | 0.80 | 1.09 | 0.664 | 0.376 | 0.378 |
rs4328262 | 699 | 714 | 1.02 | 0.87 | 1.18 | 0.947 | 0.842 | 0.843 |
rs4760648 | 676 | 683 | 1.03 | 0.89 | 1.20 | 0.275 | 0.672 | 0.670 |
rs2254210 | 704 | 748 | 0.92 | 0.79 | 1.08 | 0.402 | 0.312 | 0.304 |
rs12717991 | 703 | 731 | 0.91 | 0.79 | 1.06 | 0.522 | 0.240 | 0.261 |
rs6580642 | 701 | 747 | 1.04 | 0.85 | 1.27 | 0.914 | 0.686 | 0.688 |
rs1544410 | 700 | 742 | 1.12 | 0.96 | 1.30 | 0.244 | 0.149 | 0.147 |
rs4760733 | 687 | 702 | 0.91 | 0.78 | 1.06 | 0.202 | 0.222 | 0.223 |
rs3890733 | 706 | 753 | 1.10 | 0.94 | 1.28 | 0.304 | 0.240 | 0.248 |
SNP . | Cases . | Controls . | Allele odds ratio . | 95% Lower CL . | 95% Upper CL . | P value genotype test . | P value for allele test . | P value for trend test . |
---|---|---|---|---|---|---|---|---|
rs3782905 | 705 | 753 | 0.97 | 0.83 | 1.14 | 0.841 | 0.708 | 0.713 |
rs2283342 | 701 | 749 | 0.85 | 0.69 | 1.03 | 0.131 | 0.103 | 0.099 |
rs2239182 | 692 | 733 | 0.93 | 0.80 | 1.08 | 0.571 | 0.345 | 0.342 |
rs11574077 | 709 | 767 | 0.71 | 0.51 | 1.00 | 0.159 | 0.052 | 0.056 |
rs11574027 | 710 | 768 | 1.50 | 0.83 | 2.71 | 0.353 | 0.178 | 0.193 |
rs2238136 | 699 | 732 | 1.22 | 1.03 | 1.44 | 0.054 | 0.021 | 0.020 |
rs11574046 | 687 | 706 | 0.96 | 0.74 | 1.25 | 0.795 | 0.761 | 0.767 |
rs2107301 | 696 | 739 | 1.07 | 0.91 | 1.26 | 0.622 | 0.404 | 0.402 |
rs11168287 | 692 | 739 | 0.96 | 0.83 | 1.12 | 0.596 | 0.633 | 0.641 |
rs3819545 | 697 | 745 | 0.96 | 0.83 | 1.12 | 0.777 | 0.604 | 0.609 |
rs4237855 | 646 | 700 | 0.91 | 0.78 | 1.05 | 0.541 | 0.198 | 0.284 |
rs3847987 | 706 | 750 | 1.01 | 0.81 | 1.26 | 0.809 | 0.912 | 0.912 |
rs2239186 | 701 | 754 | 0.89 | 0.75 | 1.07 | 0.376 | 0.223 | 0.212 |
rs11168275 | 701 | 741 | 1.02 | 0.86 | 1.21 | 0.965 | 0.807 | 0.809 |
rs11574026 | 706 | 762 | 0.85 | 0.67 | 1.06 | 0.353 | 0.150 | 0.149 |
rs2853564 | 695 | 726 | 0.94 | 0.81 | 1.09 | 0.579 | 0.418 | 0.422 |
rs10875695 | 698 | 749 | 1.07 | 0.90 | 1.26 | 0.718 | 0.449 | 0.449 |
rs11568820 | 708 | 758 | 1.05 | 0.88 | 1.26 | 0.721 | 0.587 | 0.601 |
rs7299460 | 701 | 746 | 0.99 | 0.84 | 1.16 | 0.573 | 0.876 | 0.879 |
rs11168267 | 679 | 695 | 1.02 | 0.81 | 1.29 | 0.901 | 0.848 | 0.846 |
rs2239181 | 700 | 757 | 1.01 | 0.80 | 1.28 | 0.619 | 0.941 | 0.940 |
rs12721364 | 697 | 729 | 1.00 | 0.81 | 1.23 | 0.744 | 0.992 | 0.992 |
rs4516035 | 698 | 733 | 0.98 | 0.85 | 1.14 | 0.620 | 0.790 | 0.799 |
rs10875693 | 698 | 735 | 0.95 | 0.81 | 1.11 | 0.418 | 0.503 | 0.511 |
rs2189480 | 697 | 751 | 1.02 | 0.88 | 1.18 | 0.071 | 0.814 | 0.815 |
rs2239179 | 688 | 710 | 0.91 | 0.78 | 1.05 | 0.297 | 0.202 | 0.182 |
rs4760655 | 703 | 750 | 0.93 | 0.80 | 1.09 | 0.664 | 0.376 | 0.378 |
rs4328262 | 699 | 714 | 1.02 | 0.87 | 1.18 | 0.947 | 0.842 | 0.843 |
rs4760648 | 676 | 683 | 1.03 | 0.89 | 1.20 | 0.275 | 0.672 | 0.670 |
rs2254210 | 704 | 748 | 0.92 | 0.79 | 1.08 | 0.402 | 0.312 | 0.304 |
rs12717991 | 703 | 731 | 0.91 | 0.79 | 1.06 | 0.522 | 0.240 | 0.261 |
rs6580642 | 701 | 747 | 1.04 | 0.85 | 1.27 | 0.914 | 0.686 | 0.688 |
rs1544410 | 700 | 742 | 1.12 | 0.96 | 1.30 | 0.244 | 0.149 | 0.147 |
rs4760733 | 687 | 702 | 0.91 | 0.78 | 1.06 | 0.202 | 0.222 | 0.223 |
rs3890733 | 706 | 753 | 1.10 | 0.94 | 1.28 | 0.304 | 0.240 | 0.248 |
Discussion
In this study, no significant trends were observed between dietary or total vitamin D intake and risk for colorectal cancer after adjustment for known colorectal cancer risk factors (Table 1). There were also no SNPs or haplotypes within the VDR gene associated with colorectal cancer after correction for multiple comparisons. The results from the present study do not support those from a recent meta-analysis showing an inverse association between dietary vitamin D intake and colorectal cancer as well as an inverse association between the BsmI polymorphism (BB vs. bb, rs1544410) and colorectal cancer (6); however, the decreased number of subjects in the present genetic dataset may have limited our power to detect modest associations.
Previous studies that have found inverse associations between dietary vitamin D intake and colorectal cancer have controlled for similar covariates in their multivariate model (6) as was performed in the present study. The addition of calcium to our model has the largest effect on observed associations. When calcium is removed from the model, a nonsignificant inverse association is observed for the upper two quartiles of total vitamin D intake with a near-significant inverse trend (FDR Ptrend = 0.059; data not shown). Our results are similar to those from other studies which have controlled for calcium when evaluating the effect of vitamin D on colorectal cancer (7, 8). Further studies evaluating additional genes or using biomarker data may help elucidate the complex relationship between vitamin D and risk for colorectal cancer.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: C.J. Gallagher, S.M. Lesko, J.E. Muscat, T.J. Hartman, P. Lazarus
Development of methodology: J.E. Muscat
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): S.M. Lesko, J.E. Muscat, T.J. Hartman
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): J.H. Ashmore, C.J. Gallagher, T.J. Hartman, P. Lazarus
Writing, review, and/or revision of the manuscript: J.H. Ashmore, C.J. Gallagher, T.J. Hartman, P. Lazarus
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): J.H. Ashmore, S.M. Lesko
Study supervision: P. Lazarus, S.M. Lesko
Grant Support
This publication was funded by the Pennsylvania Department of Health Grant number 4100038714 (to P. Lazarus).
The costs of publication of this articlewere defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.