Abstract
Although evidence exists for a causal association between 25-hydroxyvitamin D (25(OH)D) serum levels, and multiple sclerosis (MS), the role of variation in vitamin D receptor (VDR) binding in MS is unknown. Here, we leveraged previously identified variants associated with allele imbalance in VDR binding (VDR-binding variant; VDR-BV) in ChIP-exo data from calcitriol-stimulated lymphoblastoid cell lines and 25(OH)D serum levels from genome-wide association studies to construct genetic instrumental variables (GIVs). GIVs are composed of one or more genetic variants that serve as proxies for exposures of interest. Here, GIVs for both VDR-BVs and 25(OH)D were used in a two-sample Mendelian Randomization study to investigate the relationship between VDR binding at a locus, 25(OH)D serum levels, and MS risk. Data for 13,598 MS cases and 38,887 controls of European ancestry from Kaiser Permanente Northern California, Swedish MS studies, and the UK Biobank were included. We estimated the association between each VDR-BV GIV and MS. Significant interaction between a VDR-BV GIV and a GIV for serum 25OH(D) was evidence for a causal association between VDR-BVs and MS unbiased by pleiotropy. We observed evidence for associations between two VDR-BVs (rs2881514, rs2531804) and MS after correction for multiple tests. There was evidence of interaction between rs2881514 and a 25(OH)D GIV, providing evidence of a causal association between rs2881514 and MS. This study is the first to demonstrate evidence that variation in VDR binding at a locus contributes to MS risk. Our results are relevant to other autoimmune diseases in which vitamin D plays a role.