Vitamin D-Related Genes, Blood Vitamin D Levels and Colorectal Cancer Risk in Western European Populations.
Fedirko V., Mandle HB., Zhu W., Hughes DJ., Siddiq A., Ferrari P., Romieu I., Riboli E., Bueno-de-Mesquita B., van Duijnhoven FJB., Siersema PD., Tjønneland A., Olsen A., Perduca V., Carbonnel F., Boutron-Ruault M-C., Kühn T., Johnson T., Krasimira A., Trichopoulou A., Makrythanasis P., Thanos D., Panico S., Krogh V., Sacerdote C., Skeie G., Weiderpass E., Colorado-Yohar S., Sala N., Barricarte A., Sanchez M-J., Quirós R., Amiano P., Gylling B., Harlid S., Perez-Cornago A., Heath AK., Tsilidis KK., Aune D., Freisling H., Murphy N., Gunter MJ., Jenab M.
Higher circulating 25-hydroxyvitamin D levels (25(OH)D) have been found to be associated with lower risk for colorectal cancer (CRC) in prospective studies. Whether this association is modified by genetic variation in genes related to vitamin D metabolism and action has not been well studied in humans. We investigated 1307 functional and tagging single-nucleotide polymorphisms (SNPs; individually, and by gene/pathway) in 86 vitamin D-related genes in 1420 incident CRC cases matched to controls from the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort. We also evaluated the association between these SNPs and circulating 25(OH)D in a subset of controls. We confirmed previously reported CRC risk associations between SNPs in the VDR, GC, and CYP27B1 genes. We also identified additional associations with 25(OH)D, as well as CRC risk, and several potentially novel SNPs in genes related to vitamin D transport and action (LRP2, CUBN, NCOA7, and HDAC9). However, none of these SNPs were statistically significant after Benjamini-Hochberg (BH) multiple testing correction. When assessed by a priori defined functional pathways, tumor growth factor β (TGFβ) signaling was associated with CRC risk (P ≤ 0.001), with most statistically significant genes being SMAD7 (PBH = 0.008) and SMAD3 (PBH = 0.008), and 18 SNPs in the vitamin D receptor (VDR) binding sites (P = 0.036). The 25(OH)D-gene pathway analysis suggested that genetic variants in the genes related to VDR complex formation and transcriptional activity are associated with CRC depending on 25(OH)D levels (interaction P = 0.041). Additional studies in large populations and consortia, especially with measured circulating 25(OH)D, are needed to confirm our findings.