Vitamin D and Risk of Multiple Sclerosis: A Mendelian Randomization Study
Many previous studies have shown that
- Genes can reduce (measured) levels of vitamin D
- Genes can also reduce levels of vitamin D which actually get to cells - in addition to measure levels
example: genes associated with Vitamin D Receptor (especially A-A) or Vitamin D Binding Protein, Kidney
- Many non-gene things also reduce the measured and non-measure vitamin D which get to cells
example: Low Magnesium, Low Omega-3,
- Genes can increase the risk of Multiple Sclerosis (as well as other disease)
- Increase in Vitamin D decreases the risk of Multiple Sclerosis (as well as other disease)
This study had found that if a person has the 4 major gene problems, each of which result in low levels of vitamin D, that there is a substantially increased chance of getting Multiple Sclerosis.
The 4 genes considered were DHCR7, CYP2R1,GC, CYP24A1
Only the first two of which are detected by standard vitamin D measurements
Comment by VitaminDWiki
While genes do reduce the measured amounts of vitamin D enough to increase the risk of many diseases, it appears that a person need only take higher doses of vitamin D to compensate.
See Reasons for low response to vitamin D
Note that onlysome of the reasons can be measured by vitamin D tests
- Gene differences can result in 14 ng difference in vitamin D levels– Feb 2014
- GC, CYP2R1 and DHCR7 genes associated with low vitamin D levels in China – 2012, 2013
- Common Vitamin D gene variants and resulting diseases – Jan 2013
- Vitamin D insufficiency was 3.7 X more likely if CYP2R1 gene variation– June 2014
- Continuing Medical Education (CME) for Vitamin D includes genes and questions – July 2014
- Breast Cancer rate reduced 72 percent by vitamin D gene polymorphism CYP24A1 – Nov 2014
Just 1 gene caused a 72% increase in Breast Cancer, vs this study which 4 genes increased MS by 50%
- Overview MS and vitamin D
- Genes (CYP2R1 and GC) which restrict the amount of Vitamin D which gets into bloodstream – June 2014
- All items in category Vitamin D Binding Protein (affected by GC)
- Bioavailable Vitamin D is the same blacks and whites, but measured vit D is not – Oct 2014
- Genes are one of the reasons for low response to vitamin D – Nov 2014
- Gene makes COPD 2.6X more likely unless get more vitamin D – meta-analysis Dec 2014
- similar to study on this page
- Testing misses the Vitamin D tied up with VDBP – June 2014
- Hip Fracture 2.5 times more likely if poor Vitamin D Binding Protein gene - Aug 2015
The articles in both MS and Genetics are:
- Vitamin D genes increase MS relapses in children by 2X – May 2019
- CYP2R1 gene problem increases Multiple Sclerosis risk by 1.4X – Dec 2018
- Multiple Sclerosis more likely if poor vitamin D genes - 22nd study – Aug 2017
- Mendelian proof that low vitamin D (due to 3 genes) increase risk of MS by 20 percent – Nov 2016
- Autoimmune risk gene ZMIZ1 is associated with both MS and Vitamin D – Jan 2017
- Multiple Sclerosis relapse in children is twice as likely having a Vitamin D Gene score of 6 – Oct 2016
- Multiple Sclerosis and obesity share some gene problems (as well as low vitamin D) – June 2016
- Genes make Multiple Sclerosis 2X more likely unless get more vitamin D - Aug 2015
- Multiple Sclerosis is connected to Vitamin D by gene to gene interactions – Aug 2014
- Multiple Sclerosis, gene expression, and vitamin D: Venn diagrams – Aug 2014
- Epigenetics of Multiple Sclerosis – March 2014
- Increased risk of multiple sclerosis risk in African Americans due to genes – June 2013
- 98 pcnt of genes that Vitamin D activates to reduce MS are also activated by Interferon -May 2013
- Transgeneration vitamin D deficiency related to MS was found in mice – Aug 2012
- Epigenetics, vitamin D, and Multiple Sclerosis
- Learning about MS and vitamin D in offspring from mice – Sept 2011
- Vitamin D targets 4 MS genes – May 2011
- Unable to find a gene linking vitamin D and MS – March 2011
- MS and vitamin D may be related by HLA gene – March 2010
- MS due to low level of vitamin D may be due to a specific gene – July 2010
The articles in both MS and Vitamin D Receptor are:
- Multiple Sclerosis 2X-3X more likely if poor Vitamin D Receptor – Meta-analysis Feb 2020
- Risk of Multiple Sclerosis varies with the Vitamin D Receptor – meta-analysis Dec 2019
- Multiple Sclerosis and Vitamin D Receptor super enhancers – March 2019
- Vitamin D genes increase MS relapses in children by 2X – May 2019
- Immunological effects of vitamin D and their relations to autoimmunity – March 2019
- Inflammation and immune responses to Vitamin D (perhaps need to measure active vitamin D) – July 2017
- Multiple Sclerosis more likely if poor vitamin D genes - 22nd study – Aug 2017
- Multiple sclerosis (relapsing-remitting) increases activation of Vitamin D Receptor by 6.6 X – March 2017
- Multiple Sclerosis is more likely if poor Vitamin D Receptor (4X Mexico, 3X Iran)– Feb 2017
- Multiple Sclerosis much more likely if poor Vitamin D Receptor – several studies
- Multiple Sclerosis and the Vitamin D Receptor – meta-analysis July 2014
The articles in both MS and Vitamin D Binding Protein:
- Vitamin D genes increase MS relapses in children by 2X – May 2019
- Mendelian proof that low vitamin D (due to 3 genes) increase risk of MS by 20 percent – Nov 2016
- Genes make Multiple Sclerosis 2X more likely unless get more vitamin D - Aug 2015
- Multiple Sclerosis is associated with about 1.5 X more Vitamin D Binding Protein – Jan 2015
- Late-stage MS associated with protein in spinal cord which blocks vitamin D – Jan 2013
- Cyp2R1 210 hits June 2015
- GC gene 365 hits as of Aug 2015
- CYP24A1 534 hits as of Feb 2015
- DHCR7 70 hits as of Aug 2015
- CPY27A1 279 hits as of Feb 2015
Lauren E. Mokry1, Stephanie Ross1, OmarS. Ahmad1,2, Vincenzo Forgetta1,2, George Davey Smith3, Aaron Leong4,5, Celia M. T. Greenwood6,7,8,9, George Thanassoulis2,10, J. Brent Richards1,2,7,8,email@example.com Centre for Clinical Epidemiology, Department of Epidemiology, Lady Davis Institute for Medical Research, Jewish General Hospital, McGill University, Montreal, Quebec, Canada, 2 Department of Medicine, McGill University, Montreal, Quebec, Canada, 3 MRC Integrative Epidemiology Unit, School of Social and Community Medicine, University of Bristol, Bristol, United Kingdom, 4 Division of General Internal Medicine, Massachusetts General Hospital, Boston, Massachusetts, United States of America, 5 Department of Medicine, Harvard Medical School, Boston, Massachusetts, United States of America, 6 Department of Oncology, McGill University, Montreal, Quebec, Canada, 7 Department of Epidemiology, Biostatistics and Occupational Health, McGill University, Montreal, Quebec, Canada, 8 Department of Human Genetics, McGill University, Montreal, Quebec, Canada, 9 Lady Davis Institute, Jewish General Hospital, Montreal, Quebec, Canada, 10 Preventive and Genomic Cardiology, McGill University Health Center, Montreal, QC,Department of Twin Research and Genetic Epidemiology, King’s College London, United Kingdom
Table of contents
- (PLOS) Editors' Summary
- Association of SUNLIGHT SNPs with Multiple Sclerosis Susceptibility
- Mendelian Randomization Analysis for the Association of 25OHD with Multiple Sclerosis Risk
- Supporting Information
- Additional Information
- Expert reaction to new study on vitamin D deficiency and MS
- Extensive positive comments on study in Neurology Today - Sept 2015
Observational studies have demonstrated an association between decreased vitamin D level and risk of multiple sclerosis (MS); however, it remains unclear whether this relationship is causal. We undertook a Mendelian randomization (MR) study to evaluate whether genetically lowered vitamin D level influences the risk of MS.
Methods and Findings
We identified single nucleotide polymorphisms (SNPs) associated with 25-hydroxyvitamin D (25OHD) level from SUNLIGHT, the largest (n = 33,996) genome-wide association study to date for vitamin D. Four SNPs were genome-wide significant for 25OHD level (p-values ranging from 6 x 10-10 to 2 x 10-109), and all four SNPs lay in, or near, genes strongly implicated in separate mechanisms influencing 25OHD. We then ascertained their effect on 25OHD level in 2,347 participants from a population-based cohort, the Canadian Multicentre Osteoporosis Study, and tested the extent to which the 25OHD-decreasing alleles explained variation in 25OHD level. We found that the count of 25OHD-decreasing alleles across these four SNPs was strongly associated with lower 25OHD level (n = 2,347, F-test statistic = 49.7, p = 2.4 x 10-12). Next, we conducted an MR study to describe the effect of genetically lowered 25OHD on the odds of MS in the International Multiple Sclerosis Genetics Consortium study, the largest genetic association study to date for MS (including up to 14,498 cases and 24,091 healthy controls). Alleles were weighted by their relative effect on 25OHD level, and sensitivity analyses were performed to test MR assumptions. MR analyses found that each genetically determined one-standard-deviation decrease in log- transformed 25OHD level conferred a 2.0-fold increase in the odds of MS (95% CI: 1.7—2.5; p = 7.7 x 10-12; /2 = 63%, 95% CI: 0%-88%). This result persisted in sensitivity analyses excluding SNPs possibly influenced by population stratification or pleiotropy (odds ratio [OR] = 1.7, 95% CI: 1.3-2.2; p = 2.3 x 10-5; /2 = 47%, 95% CI: 0%-85%) and including only SNPs involved in 25OHD synthesis or metabolism (ORsynthesis = 2.1,95% CI: 1.6-2.6, p = 1 x 10-9; ORmetaboiism = 1.9, 95% CI: 1.3-2.7, p = 0.002). While these sensitivity analyses decreased the possibility that pleiotropy may have biased the results, residual pleiotropy is difficult to exclude entirely.
A genetically lowered 25OHD level is strongly associated with increased susceptibility to MS. Whether vitamin D sufficiency can delay, or prevent, MS onset merits further investigation in long-term randomized controlled trials.
Multiple sclerosis (MS) is a debilitating neurological disorder that affects the nerves in the brain and spinal cord and that usually presents between the ages of 20 and 40 years. It is an autoimmune disease—a disease in which the body’s immune system attacks healthy tissues. In the case of MS, the immune system attacks myelin, the fatty tissue that forms an insulating sheath around nerves fibers. Multiple areas of scarring (sclerosis) in the myelin are produced by immune attack, and the damage slows down or blocks the transmission of electrical signals along the nerves to and from the brain and causes symptoms such as loss of vision, muscle stiffness, and daytime fatigue. For most affected individuals, MS initially follows a relapsing-remitting course characterized by periods lasting days or weeks during which symptoms flare up, followed by periods during which symptoms becomes milder or disappear. Half of people with relapsing-remitting MS eventually develop secondary progressive MS, in which their symptoms inexorably worsen. There is currently no cure for MS, but some treatments can relieve its symptoms or reduce the number of relapses.
Why Was This Study Done?
Why the immune system attacks myelin to cause MS is unclear but probably involves both genetic and environmental risk factors. One potential environmental risk factor is reduced levels of vitamin D. Circulating levels of 25-hydroxyvitamin D (25OHD; the clinical determinant of vitamin D status) are determined in part by exposure to sunlight, and MS is more common at higher latitudes, where exposure to sunlight is decreased. Other epidemiological studies (investigations that examine disease patterns in populations) also suggest an association between lower vitamin D level and an increased risk of MS but cannot prove that a decreased vitamin D level actually causes MS. Individuals who develop MS might share another unknown characteristic that increases their risk of MS (confounding), or individuals who have MS might spend less time outdoors and, as a result, have lower circulating vitamin D levels (reverse causation). It is important to know whether a decreased vitamin D level increases the risk of MS because vitamin D insufficiency is becoming increasingly common. Here, the researchers undertake a Mendelian randomization study to determine whether circulating vitamin D level has a causal effect on MS susceptibility. Because gene variants are inherited randomly, they are not prone to confounding. Reverse causation is also prevented since MS does not change genetic variants. So, if vitamin D level actually affects MS risk, genetic variants that affect vitamin D level should be associated with altered susceptibility to MS.
What Did the Researchers Do and Find?
The researchers first ascertained the effect on 25OHD level among participants in the Canadian Multicentre Osteoporosis Study of four single nucleotide polymorphisms (SNPs, a type of genetic variant) that were associated with 25OHD level in a genome-wide association study (SUNLIGHT). Each of the SNPs explained an important proportion of the population-level variance in 25OHD level in the Canadian Multicentre Osteoporosis Study. The researchers then used the SNPs to examine whether there was an association between genetically reduced 25OHD level and susceptibility to MS among participants in the International Multiple Sclerosis Genetics Consortium study, a genome-wide association study involving up to 14,498 people with MS and 24,091 healthy controls. They found that a genetic decrease in the natural-log-transformed 25OHD level by one standard deviation was associated with a 2-fold increased risk of MS. In practical terms, this finding means that increasing an individual’s circulating 25OHD level by approximately 1.5-fold decreases their odds of developing MS by 50%.
What Do These Findings Mean?
These findings show that, among the participants of the International Multiple Sclerosis Genetics Consortium study, all of whom were of European ancestry, genetically lowered 25OHD level was strongly associated with increased susceptibility to MS. Although the Mendelian randomization approach used here largely avoids the possibility of confounding or reverse causation, the reliability of these findings may be limited by some of the assumptions made by the researchers during their analysis. Moreover, although these findings support a role for vitamin D in MS susceptibility, they provide no information about whether vitamin D modulates the course of MS after its onset, and they may not apply to people of non-European ancestry. Nevertheless, these findings provide a strong rationale for undertaking randomized controlled trials to investigate whether vitamin D supplementation can prevent the onset and/or progression of MS.
Multiple sclerosis (MS) is the most common permanent neurological disorder affecting young adults . It is a debilitating autoimmune condition that presents early in life, with a mean age of onset of 28-31 y. Epidemiological studies have indicated that the prevalence of MS varies geographically, such that regions of higher latitude and with decreased levels of sunlight exposure have a higher prevalence of MS [2,3]. Since the circulating level of vitamin D, as measured by the level of 25-hydroxyvitamin D (25OHD, the clinical determinant of vitamin D status), is partially derived from sunlight exposure, it has been suggested that 25OHD deficiency may be the causal risk factor mediating this latitudinal gradient . Further evidence to support the vitamin D hypothesis arose from the Nurses’ Health Study, which reported a protective effect on MS for women who had high levels of daily vitamin D intake . Lower vitamin D level has also been associated with higher rates of MS relapse  and higher MS-specific disease activity and disability . Vitamin D has important effects upon the immune system, and its immune- modulating effects have been observed in multiple cell-culture experiments , providing possible biological mechanisms whereby vitamin D may influence MS risk.
To date, there has been one published meta-analysis investigating the effect of vitamin D supplementation on MS relapse; the meta-analysis included five randomized controlled trials (RCTs) with a total of 254 participants . The authors reported that the effect of high-dose vitamin D treatment on MS relapse was inconclusive (odds ratio [OR] = 0.98, 95% CI: 0.452.16) and that these trials had important methodological limitations, such as small sample size and short duration of vitamin D treatment. In contrast, two non-blinded trials demonstrated improved clinical outcomes with vitamin D therapy; however, disease activity or MRI changes were not the primary outcome of these trials [10,11]. Importantly, all of these trials test whether vitamin D can treat MS, but provide no insight into whether vitamin D can prevent MS.
Consequently, clinical practice guidelines for the treatment of MS  do not include vitamin D therapy. This is at least partially attributable to the possibility of confounding in the above observational studies. Additionally, observational studies are prone to reverse causation, where, for example, individuals with MS may spend less time outdoors and may as a result have lower circulating 25OHD levels. However, if decreased 25OHD levels are causallyassociated with MS, this could have important implications since vitamin D insufficiency, defined as 25OHD level < 50 nmol/l, is common and increasing in prevalence. This was observed in the National Health and Nutrition Examination Survey (NHANES): in 2005,
41.6% of adult Americans were found to be vitamin D insufficient, and mean 25OHD level decreased from 75 nmol/l in 1988 to 50 nmol/l in 2006 [ 13,14].
In the absence of high-quality RCT data, the principles of Mendelian randomization (MR) can be applied to strengthen or refute the causality of biomarkers in disease etiology . MR analysis uses genetic associations to test the effects of biomarkers, such as 25OHD, on the risk of disease. This approach, which is conceptually similar to an RCT, is based on the principle that genetic variants are randomly allocated at meiosis, and consequently these genetic variants are independent of many factors that bias observational studies, such as confounding and reverse causation. MR methods have been used previously to investigate the role of high-density lipoprotein  and C-reactive protein  in predisposition to cardiovascular disease, and have provided strong evidence that PCSK9 inhibition prevents cardiovascular disease . MR methods may be of particular relevance for understanding the etiology of MS since the date of disease onset is often poorly recognized clinically and MR studies assess the effect of lifetime exposures.
Here we adopted an MR design to clarify whether 25OHD level lies in the causal pathway for MS susceptibility. In order to assess whether a reduced level of 25OHD is associated with an increased risk of MS, we selected genome-wide significant single nucleotide polymorphisms (SNPs) as identified by SUNLIGHT (Study of Underlying Genetic Determinants of Vitamin D and Highly Related Traits), the largest genome-wide association study (GWAS) published to date for 25OHD level. Next, we estimated the effect of each of these SNPs upon 25OHD level in the Canadian Multicentre Osteoporosis Study (CaMos) and tested their validity as instrumental variables for MR analyses. Finally, we applied the principles of MR to investigate the association of a lifetime of genetically lowered 25OHD level with MS risk using data from the International Multiple Sclerosis Genetics Consortium (IMSGC).
SNP Selection and Data Sources
Genetic variants associated with 25OHD level at a genome-wide significant level (p < 5 x 10-8) were obtained from SUNLIGHT , a GWAS consisting of 33,996 individuals of European descent from 15 cohorts. 25OHD level in this study was measured by radioimmunoassay, chemiluminescent assay, ELISA, or mass spectrometry. Given that different cohorts used different methods to measure 25OHD level, results were combined across cohorts in SUNLIGHT using Z-score-weighted meta-analysis.
CaMos was used to estimate the effect of each genome-wide significant SNP on 25OHD level, since the effect of each SNP upon 25OHD level could not be used from SUNLIGHT, because of the Z-score meta-analytic approach employed . CaMos is a large population- based cohort and was among the largest included in the replication phase of SUNLIGHT. It includes 2,347 individuals who were genotyped using TaqMan genotyping at the same genome-wide significant vitamin D loci found in SUNLIGHT.
To obtain precise estimates for the association of 25OHD with MS, we tested the effect of each genome-wide significant SNP for vitamin D level in the IMSGC Immunochip study, the largest international genetic study of MS, involving 14,498 MS cases and 24,091 healthy controls . All participants were of European ancestry and were genotyped using the Immuno- chip array, which is a custom array designed to interrogate SNPs with potential immune system effects. Cases were defined as individuals diagnosed by a neurologist according to recognized diagnostic criteria dependent on laboratory and clinical information [22-24]. When data were not available for a specific SNP in the IMSGC Immunochip study, we used data from the second largest MS genetic study, the IMSGC study and IMSGC/Wellcome Trust Case Control Consortium 2 (IMSGC/WTCCC2) study, which included 9,772 cases and 6,332 controls taken from the IMSGC/WTCCC2 common control set .
SNP Validation and Effect Sizes
Linkage disequilibrium assessment. One requirement of MR studies is that the selected SNPs must not be in linkage disequilibrium (LD) since if a selected SNP is highly correlated with other risk factor loci, this may result in confounding . In order to verify that the SNPs in this study met this requirement, we measured LD between all selected SNPs using CEU samples from the 1000 Genomes Project (n = 94) .
Pleiotropy assessment. MR analyses assume that the chosen SNPs do not exert pleiotropic effects on the outcome (in this case, MS) by operating through biological pathways independent of the exposure (in this case, 25OHD level). However, in MR, a SNP may influence the outcome via other factors if the SNP acts upon the other factors through the exposure itself . Previous work has assessed possible pleiotropic actions of the 25OHD-related SNPs used in our analysis by investigating the association between 25OHD-related SNPs and clinical traits in the 1958 British Birth Cohort, which included 6,877 participants of European descent . In this cohort, no associations were found between these SNPs and relevant potential pleiotro- pic pathways, such as sun exposure, time outside, physical activity, oily fish consumption, smoking, alcohol consumption, body mass index (BMI), abdominal obesity, or social class (p > 0.05 for all) . However, we note that some of these factors, such as sun exposure, time outside, BMI, and abdominal obesity, could act at least partially through the vitamin D pathway. Furthermore, SNPs associated with 25OHD level did not associate with other biomarkers (including C-reactive protein, IgE level, von Willebrand factor, tissue plasminogen activator, D-dimer, fibrinogen, triglyceride level, high-density lipoprotein, low-density lipoprotein, total cholesterol, forced expiratory volume, diastolic blood pressure, IGF-1, and HbA1c), and no interactions were observed between the SNPs, biomarkers, and 25OHD level. Additional details are provided in S1 Table.
To further explore sources of pleiotropy, we also conducted a systematic literature search of gene name, gene mutation, and protein name to examine the published literature for possible pleiotropic mechanisms for any of our selected SNPs for effects on MS and autoimmunity using PubMed. Details of this method are described in S1 Methods.
Population stratification assessment. The 1958 British Birth Cohort has previously assessed the potential for population stratification of the 25OHD-associated SNPs, which is a potential source of bias in MR studies since differences in minor allele frequencies between populations may cause the SNP to be associated with both the ancestry and the outcome . In the 1958 British Birth Cohort, each SNP was tested for association with geographic region, which was dichotomized as south and middle UK (southeast England, southwest England, greater London, East Anglia, Midlands, and Wales) versus northern England (north England, northwest England, Yorkshire, and the Humber) and Scotland. We then further assessed potential population stratification of each SNP by testing its association with self-declared ethnicity in the CaMos cohort and tested the association of each SNP with non-European status, defined as exclusion from the European cluster in principal component analysis (PCA).
Effect size estimates of SUNLIGHT SNPs upon 25OHD level. To obtain the effect of each SNP upon 25OHD level, as required for MR analysis, we tested the additive effect of each minor allele on natural-log-transformed 25OHD level in CaMos, while controlling for sex, age, age squared, BMI, and season of 25OHD measurement (using categorical variables for summer [July-September], autumn [October-December], winter [January-March], and spring [April- June]) . Ethnicity was checked by self-report and verified using PCA. To prevent population stratification from confounding our results, individuals who did not cluster with other Europeans were excluded from this analysis and were not used to measure the effect of each SNP upon 25OHD. The number of 25OHD-decreasing alleles was calculated for each participant in the CaMos cohort. This allele count was tested for an association with natural-log- transformed 25OHD level using linear regression, which had been residualized for the above covariates, and the F-statistic for the allele score was reported. The multiply adjusted natural- log-transformed 25OHD level was then assessed for each category of allele count, and a non- parametric trend test across these allele counts was computed.
In order to increase study power and obtain the most precise estimate of the association of 25OHD-associated SNPs with risk of MS, we used summary-level data from the IMSGC Immunochip study, if available (as described above). However, the IMSGC Immunochip geno- typing array used was not genome-wide, so not all SNPs were captured in this experiment. If a SNP was not included in the IMSGC Immunochip study, then summary statistics from the second largest genotyped cohort, the IMSGC/WTCCC2 cohort, were selected. In the event that a SNP was not genotyped in either cohort, summary statistics for a perfect proxy SNP, defined as a surrogate SNP with perfect LD (r2 = 1.0) to the SNP interest, were selected. LD for proxy SNPs was calculated using CEU samples from the 1000 Genomes Project (n = 94) since the IMSGC samples are of the same ancestry . We then assessed whether each SNP was associated with the risk of MS, applying a Bonferroni correction, where statistical significance was declared at p < 0.05/n where n is the number of SNPs associated with 25OHD level from SUNLIGHT.
Mendelian Randomization Estimates
We conducted our MR analysis by assessing the effects of the SNPs upon risk of MS, weighting the effect of each SNP by the magnitude of its effect upon 25OHD level. In this study design, which has been described previously [29-31], the independent SNPs evaluate the association of exposure to genetically lowered 25OHD with MS risk. The individual estimates of the effect were then pooled using statistically efficient estimators formally analogous to those of inverse- variance-weighted meta-analysis . We carried out a meta-analysis of estimates obtained from individual 25OHD-decreasing alleles using both fixed-effects and random-effects models to obtain pooled estimates for the combined effect of the 25OHD SNPs on MS.
Specifically, let x and y denote the centered and scaled natural-log 25OHD and log-odds MS traits, respectively, and suppose these are related by the linear structural equation: y = ax + n. Here, n is a stochastic error term, and in general x and n are correlated because of confounding. The parameter a quantifies the causal effect of x on y, and is thus the parameter we seek to estimate. Let u, denote the allele dosage variable of the ith genetic variant. Let y, and p, denote effect size estimates (derived from GWAS data) of u, on the exposure x (change in natural-log 25OHD level) and outcome y (change in log odds of MS), respectively, and let s(P,) denote the standard error of p,. Then the MR estimate associated with the ith genetic variant is a, = p,//,-, and the variance of this estimate is v, = (s[Pi])/y,)2. Define the precision of the ith MR estimate of a by w, = 1/v,. The inverse-variance-weighted fixed-effects estimate is then
We observe that afixed may also be interpreted as the regression coefficient resulting from the generalized linear regression of the outcome effect sizes p, on the exposure effect sizes y, assuming heteroskedastic errors; in this regression, the ith error term has a variance equal to s (p,)2, and the offset coefficient in the regression is constrained to be zero.
The random-effects estimate arandom and its standard error s(arandom) were also constructed from the individual estimates using standard methods  in which the weights are adjusted to account for the intrinsic variability (or heterogeneity) in the effect size. Heterogeneity may be quantified with the parameter I2, which reports the fraction of the total variance in the meta-analytic estimate that is due to intrinsic variability in the effect size, as distinct from variability arising from measurement error . The random-effects estimate arandom and its standard error s(arandom) are given by equations analogous to those for afixed and s(afixed), in which the weights assigned to individual estimates are adjusted to take into account heterogeneity in the effect size.
For all MR meta-analyses, we report estimates using both fixed-effects and random-effects models. The effect size for each meta-analysis is reported in the main results as the effect of a one-standard-deviation (1-SD) change in natural-log-transformed 25OHD level, since this metric is more interpretable than an arbitrary difference. This measure is given by exp(afixed) for the fixed-effects model and by exp(arandom) for the random-effects model. We also report the I2 as an assessment of heterogeneity.
In order to provide a better clinical interpretation of a 1-SD change in natural-log-transformed 25OHD level, we selected three different clinically relevant 25OHD thresholds for vitamin D status (<25 nmol/l for vitamin D deficiency, <50 nmol/l for vitamin D insufficiency, and >75 nmol/l for vitamin D sufficiency) . These thresholds were converted to the natural log scale because the magnitude of a 1-SD change is not constant on the untransformed scale. For each of these natural-log-transformed 25OHD levels, we then calculated a 1-SD increase in natural-log-transformed 25OHD. To obtain 25OHD levels that correspond to circulating levels in units of nanomoles/liter, we then back-transformed these values.
MR estimates were recalculated after exclusion of SNPs potentially influenced by pleiotropy or population stratification. Since SNPs associated with 25OHD level in SUNLIGHT influence either 25OHD synthesis or 25OHD metabolism , we elected to perform a stratified MR analysis where SNPs involved in either 25OHD synthesis or metabolism were analyzed separately.
Schematic representation of Mendelian randomization analysis. The leftmost box lists SNPs that were genome-wide significant for 25OHD level in SUNLIGHT (n = 33,996). The blue arrow represents the effect of SNPs on multiply adjusted natural-log-transformed 25OHD level using data from CaMos (n = 2,347). The green arrow represents the causal association of decreased 25OHD level with the risk of MS using data from the largest genetic association study to date for MS (the IMSGC Immunochip study, up to 14,498 cases and 24,091 healthy controls).
SNP Selection and Validation
SNP selection. A schematic representation of the MR study design is presented in Fig 1. SUNLIGHT identified four SNPs as genome-wide significant for 25OHD level . These included rs2282679 in GC (association with 25OHD: p = 1.9 x 10-109), rs12785878 near DHCR7 (p = 2.1 x 10-27), rs10741657 near CYP2R1 (p = 3.3 x 10-20), and rs6013897 in CYP24A1 (p = 6.0 x 10-10). We selected these SNPs for our MR study since all are strongly associated with 25OHD level and map to genes implicated in the modulation of 25OHD level through distinct mechanisms . Specifically, GC encodes the vitamin D binding protein (DBP), a group-specific component of serum globulin. DBP acts as the principal protein carrier for 25OHD, transporting 80%-90% of 25OHD to target organs [37-39]. The DHCR7 gene product is known to convert 7-dehydrocholesterol to cholesterol, providing a substrate for vitamin D production. CYP2R1 is a regulator of 25OHD synthesis through 25-hydroxylation of vitamin D in the liver, the first activation step , and, lastly, CYP24A1 inactivates 1a25 (OH)2D, rendering inactive the active form of vitamin D (Fig 2). Therefore, all SNPs used in this study map near genes strongly implicated in vitamin D synthesis, transport, or metabolism. Notably, all four SNPs lie in intergenic or intronic regions, and presently the exact effect of each SNP on these enzymes is unknown. Nevertheless, all SNPs reside near genes strongly implicated in vitamin D synthesis or metabolism .
Linkage disequilibrium and pleiotropy assessment. There was no evidence of LD between any of the four 25OHD SNPs (all pairwise r2 < 0.01) in the 1000 Genomes Project CEU samples. We note that only two of the SNPs, rs10741657 and rs12785878, were located on the same chromosome, which greatly decreases the risk of confounding by LD. As described above, none of the four 25OHD SNPs was associated with relevant pleiotropic pathways in the 1958 British Birth Cohort.
Undertaking a literature review for possible pleiotropic pathways, we found no evidence for pleiotropic mechanisms for the vitamin D metabolism SNPs, rs10741657 (CYP2R1) and rs6013897 (CYP24A1). For rs2282679 (GC), we found that its encoded protein, DBP, has been associated with macrophage activation and may modulate T cell response to vitamin D . Elevated DBP levels are found in the cerebrospinal fluid of patients with Alzheimer disease  and MS , and have been linked to the progression of MS in rats . It has been argued that DBP can act independently of vitamin D to produce clinical phenotypes; therefore, we undertook sensitivity analyses excluding rs2282679 (GC) in our MR analyses. Genetic variation in DHCR7 appears to cause Smith-Lemli-Opitz syndrome, a clinical phenotype relating to cholesterol deficiency. Given that a recent study suggested an interdependence of cholesterol and vitamin D pathways in the etiology of MS , we queried the association of rs12785878 in the largest publically available GWAS results for lipids, those of the Global Lipids Genetics Consortium , and found that this SNP was associated with a minimum p-value of 0.043 across all lipid traits, suggesting that the SNP is not strongly associated with cholesterol.
Population stratification assessment. Previous reports from the 1958 British Birth Cohort demonstrated that rs12785878 (DHCR7) was associated with geographic region . Since rs12785878 is unevenly distributed across geography and the prevalence of MS varies by geographic location, a potential surrogate for local ancestry , we tested whether this SNP was associated with non-European status in CaMos using PCA. The SNP rs12785878 was strongly associated with non-European status in the CaMos cohort (p = 2.7 x 10-13). No other SNP showed any evidence of correlation with non-European status (p > 0.5 for all other SNPs). Given this possible relationship with population stratification, we undertook MR sensitivity analyses excluding the rs12785878 (DHCR7) variant.
Association of SUNLIGHT SNPs with 25OHD Level
Table 1 displays the four SNPs that achieved genome-wide significance for 25OHD level in SUNLIGHT and describes their association with 25OHD . Each of these SNPs explained an important proportion of the population-level variance in 25OHD level, as reflected by the F-statistic.
Vitamin D pathway. In blue are the genes containing, or in proximity to, SNPs that were genome-wide significant for 25OHD level in SUNLIGHT (n = 33,996). The p-values for the association with 25OHD level were 1.9 x 10-109 for GC, 2.1 x 10-27for DHCR7, 3.3 x 10-20for CYP2R1, and 6.0 x 10-10 for CYP24A1. Note that each gene plays an independent role in modulating the level of 25OHD. Kidney and liver images credit: https://openclipart.org/.
Table 1. Characteristics of SNPs used as instrumental variables.
The count of 25OHD-decreasing alleles across these four SNPs was strongly associated lower 25OHD level in the CaMos population, residualized for age, season, sex, and BMI (F-statistic = 49.7, r2 = 2.44%, p for allelic score = 2.4 x 10-12). Fig 3 shows the mean 25OHD levels for individuals with increasing counts of 25OHD-decreasing alleles (non-parametric trend test, p = 3.3x 10-19).
Association of SUNLIGHT SNPs with Multiple Sclerosis Susceptibility
Summary statistics for two of the four 25OHD-associated SNPs (rs10741657 at CYP2R1 and rs12785878 at DHCR7) and their association with MS were taken from the IMSGC Immuno- chip study (Table 1). rs12785878 at DHCR7 was not directly genotyped in the IMSGC Immu- nochip study; however, a perfect proxy for rs12785878, rs4944958, was used (r2 = 1.0 between rs12785878 and rs4944958 in the 1000 Genomes Project CEU samples). Summary statistics for the remaining two SNPs (rs6013897 at CYP24A1 and rs2282679 at GC) were taken from the second largest MS genetic association study, the IMSGC/WTCCC2 study. SNP rs6013897 at CYP24A1 was not present in the IMSGC/WTCCC2 dataset, and therefore a perfect proxy SNP for rs6013897, rs17217119, was used (r2 = 1.0 between rs17217119 andrs6013897 from the 1000 Genomes Project CEU samples).
All four 25OHD-decreasing alleles were associated with an increased risk of MS (Table 1). rs12785878 (DHCR7) achieved genome-wide significance for MS risk, while two 25OHD-decreas- ing alleles (rs10741657 and rs6013897) were moderately associated with MS risk (p = 3.9 x 10-3 andp = 1.7 x 10-3, respectively). The 25OHD-decreasing allele rs2282679 (allele C) (GC) was not significantly associated with MS risk (p = 0.062) (fable 1). However, three of the 25OHD-associ- ated SNPs (rs12785878, rs10741657, and rs6013897) remained associated with MS after a Bonfer- roni correction for the number of independent SNPs (p < 0.05/4 = 0.0125).
Fig 3. 25OHD level by number of 25OHD-decreasing alleles in theCaMos cohort. Here we show the box-plot of natural-log-transformed 25OHD by the count of 25OHD-decreasing alleles in the CaMos population. A count of zero represents individuals with no 25OHD-decreasing alleles (or homozygous at each loci for the 25OHD-increasing allele), and a count of six represents an individual with six 25OHD-decreasing alleles. No individuals with a count of seven or more 25OHD-decreasing alleles were observed in this cohort. The center line and error bars represent the mean level of natural-log-transformed 25OHD and its 95% CI for each respective allele count. Note a negative trend between allele count and mean natural-log-transformed 25OHD. doi:10.1371/journal.pmed.1001866.g003
In order to estimate the association of genetically lowered 25OHD with MS, we used a fixed- effects model in which all four 25OHD-decreasing alleles of the MR set were included. We observed that each 1-SD decrease in natural-log-transformed 25OHD level was associated with an increased risk of MS (OR = 2.02, 95% CI: 1.65-2.46, p = 7.72 x 10-12) (Table 2; Fig 4). Given that the I2 estimate of heterogeneity was somewhat increased (I2 = 63%, 95% CI: 0%-88%), we also undertook random-effects meta-analysis, which generated similar findings (OR = 2.07, 95% CI: 1.45-2.96,p = 5.74 x 10-5) (Table 2; S1 Fig). We note that since our model included only four SNPs, the 95% CIs of the I2 statistic are wide, and consequently heterogeneity cannot be accurately measured using this parameter. In addition, to address the potential effects of population stratification and pleiotropy, we undertook a sensitivity analysis excluding the rs12785878 SNP (DHCR7). Despite removal of this variant, we observed a clear association of genetically lowered 25OHD level with the risk of MS (OR = 1.72, 95% CI: 1.34-2.21, p = 2.28 x 10-5; I2 = 47%, 95% CI: 0%-85%) (Table 3; Fig 5), which remained significant using a random-effects meta-analysis (OR = 1.82, 95% CI: 1.24-2.67, p = 2.13 x 10-3; I2 = 47%, 95% CI: 0%-85%) (Table 3; S2 Fig). Removal of the rs2282679 SNP (GC), which may possibly be influenced by pleiotropy, did not influence the MR results using a fixed-effects or random- effects model (OR = 2.17, 95% CI: 1.73-2.72,p = 1.7 x 10-11; I2 = 67%, 95% CI: 0%-91%; and OR = 2.32, 95% CI: 1.49-3.61, p = 1.8 x 10-4; I2 = 67%, 95% CI: 0%-90%, respectively) (S3 and S4 Figs). To further assess the effect of the independent vitamin D pathways on the risk of MS, we analyzed SNPs near genes implicated in 25OHD synthesis (DHCR7 and CYP2R1) and metabolism (GC and CYP24A1) separately and found that both strongly associated with increased risk of MS (Table 4; S5 and S6 Figs).
The clinical equiva lences of a 1-SD increase in natural-log-transformed 25OHD for the thresholds for vitamin D deficiency (<25 nmol/l), vitamin D insufficiency (<50 nmol/l), and vitamin D sufficiency (>75 nmol/l) are shown in Table 5. We observed that for individuals at the threshold for vitamin D deficiency (25OHD = 25 nmol/l), an increase in 25OHD level to 36.9 nmol/l would be required to decrease the odds of MS by 50%, while for individuals at the thresholds for vitamin D insufficiency (25OHD = 50 nmol/l) and vitamin D sufficiency (25OHD = 75 nmol/l), an increase in 25OHD level to 73.7 nmol/l and 110.6 nmol/l, respectively, would similarly be required.
Using summary-level data for MS and 25OHD levels from large European populations, our study demonstrated that a genetic decrease in natural-log-transformed 25OHD by 1 SD was associated with a 2-fold increase in risk of MS, providing strong evidence in support of a causal role of vitamin D in MS susceptibility. These findings are consistent with evidence from observational studies that have demonstrated that low vitamin D levels influence risk of MS and also reflect findings from functional studies that have implicated vitamin D as an important regulator in the expression of MHC class II genes [47,48]. This evidence provides rationale to further investigate whether vitamin D supplementation may reduce MS susceptibility in those most at risk.
The identification of vitamin D as a causal susceptibility factor for MS may have important public health implications since vitamin D insufficiency is common [13,14], and vitamin D supplementation is both relatively safe and cost-effective . The importance of these findings may be magnified in high-latitude countries, which have disproportionately higher rates of MS and also higher rates of vitamin D insufficiency.
A reasonable first step to understanding the role of vitamin D therapy in delaying the onset or severity of MS would be to treat vitamin D insufficiency in those most at risk of developing MS. MS is often preceded by clinically isolated syndrome, which is a first clinical episode compatible with MS, often accompanied by lesions on magnetic resonance imaging , thereby providing a therapeutic window and rationale for intervening with vitamin D supplementation. Ongoing RCTs are currently assessing vitamin D supplementation for the treatment and prevention of MS [49,50] and may therefore provide needed insights into the role of vitamin D supplementation.
An important difference between MR studies and RCTs is that MR studies describe the association of a lifetime of exposure to vitamin-D-lowering alleles in the general population, whereas RCTs provide insights from supplementation for shorter periods in individuals at risk. Thus, long-term RCTs may be needed to adequately assess the impact of vitamin D supplementation in the prevention or treatment of MS. Lastly, MR may be an ideal study design to understand risk factors for MS, given the long latency period between disease onset and diagnosis, since MR may permit the estimation of lifetime exposure to risk factors.
Table 4. Mendelian randomization estimate of the association of decreased 25OHD with the risk of multiple sclerosis stratified by SNPs near genes involved in 25OHD synthesis versus metabolism using a fixed-effects model.
Our analysis has several strengths. First, by utilizing the random allocation of genetic variants, we were able to overcome potential confounding and reverse causation that may bias estimates from observational studies. Second, using data from the largest genetic consortia for 25OHD level (n = 33,996) and MS risk (up to 14,498 cases and 24,091 controls) has enabled us to more precisely test our study hypothesis than if we had used individual-level data from a small study. Previous work has shown that using estimates from meta-analytic data for uncorrelated genetic variants is similarly efficient to using individual-level data in MR studies . Lastly, the findings from this study represent the association of a lifelong exposure to reduced vitamin D levels with MS in the general European population, and, in the absence of large- scale, long-term RCT data, our findings provide strong evidence in support of a causal role of low vitamin D levels in MS susceptibility.
Our study also has limitations. First, while we have provided evidence supporting a role for vitamin D in MS susceptibility, we cannot conclude that vitamin D plays a role in disease modulation after disease onset. While MR is able to overcome the limitations that may bias observational studies, the possibility of residual pleiotropy could bias estimates in this study. However, in this study the main findings remained robust in multiple sensitivity analyses testing the plei- otropy assumption, thereby decreasing the probability of bias due to pleiotropy. We also note that all four studied SNPs are located in or near 25OHD-associated genes and influence 25OHD levels through known and distinct mechanisms. Additionally, the point estimate for each 25OHD-decreasing allele, as well as the combined 25OHD synthesis and metabolism pathways, was independently associated with increased risk of MS. Therefore, it is unlikely that pleiotropy strongly biased our results. Like in most MR studies, we cannot directly assess whether canalization, which is defined as compensatory feedback interactions, may have influenced our results [15,51,52]. However, since canalization assumes that other physiological mechanisms may attenuate the effect of genetically reduced 25OHD levels, such feedback interactions would tend to bias results toward the null. In contrast, our study has generated results that are very distinct from the null.
Table 5. Clinical equivalence of a 1-SD natural-log increase in 25OHDfor various vitamin D thresholds.
MR analyses using DHCR7, GC, CYP24A1, and CYP2R1 as instruments have been performed in the past [53-57]. We and others recently provided evidence from MR that low vitamin D levels do not increase insulin resistance  or the risk of type 2 diabetes [53,54] or coronary heart disease , but do increase the risk of type 1 diabetes  and possibly blood pressure . Interestingly, MR has shown that 25OHD levels are directly influenced by BMI, and converse effects are likely to be small . Thus, while observational associations between 25OHD and two autoimmune conditions—type 1 diabetes and now MS—have been supported by genetic evidence, associations with cardio-metabolic outcomes have not been supported thus far.
In conclusion, using data from the largest existing genetic consortia, we demonstrate that genetically lowered 25OHD level is associated with an increase in the risk of MS in people of European descent. These findings provide rationale for further investigating the potential therapeutic benefits of vitamin D supplementation in preventing the onset and progression of MS.
S1 Methods. Details of our systematic PubMed search exploring possible sources of pleiot- ropy. (PDF)
Fig. MR estimate of the association of 25OHD level with risk of MS using a random- effects model. (PDF)
Fig. MR estimate of the association of 25OHD level with risk of MS excluding the DHCR7 locus using a random-effects model. (PDF)
Fig. MR estimate of the association of 25OHD level with risk of MS excluding the GC locus using a fixed-effects model. (PDF)
Fig. MR estimate of the association of 25OHD level with risk of MS excluding the GC locus using a random-effects model. (PDF)
Fig. MR estimate of the association of 25OHD level with risk of MS for SNPs involved in 25OHD synthesis using a fixed-effects model. (PDF)
Fig. MR estimate of the association of 25OHD level with risk of MS for SNPs involved in 25OHD metabolism using a fixed-effects model. (PDF)
S1 Table. Summary of previous work assessing the association between 25OHD SNPs and relevant biomarker pathways. (PDF)
We wish to kindly thank the SUNLIGHT Consortium, IMSGC, and the IMSGC/WTCCC2 study for access to their data.
Conceived and designed the experiments: JBR. Performed the experiments: JBR LEM. Analyzed the data: JBR LEM SR. Contributed reagents/materials/analysis tools: JBR. Wrote the first draft of the manuscript: JBR LEM. Contributed to the writing of the manuscript: LEM JBR SR OSA VF AL CMTG GT GDS. Agree with the manuscript’s results and conclusions: LEM JBR SR OSA VF AL CMTG GT GDS. Created all figures: SR. All authors have read, and confirm that they meet, ICMJE criteria for authorship.
- Ramagopalan SV, SadovnickAD. Epidemiology of multiple sclerosis. Neurol Clin. 2011; 29:207-217. doi: 10.1016/i.ncl.2010.12.010 PMID: 21439437
- Simpson S, Blizzard L, Otahal P, Van der Mei I, Taylor B. Latitude is significantly associated with the prevalence of multiple sclerosis: a meta-analysis. J Neurol Neurosurg Psychiatry. 2011; 82:11321141. doi: 10.1136/innp.2011.240432 PMID: 21478203
- Hernan MA, Olek MJ, Ascherio A. Geographic variation of MS incidence in two prospective studies of US women. Neurology. 1999; 53:1711-1718. doi: 10.1212/WNL.53.8.1711 PMID: 10563617
- Ebers GC. Environmental factors and multiple sclerosis. Lancet Neurol. 2008; 7:268-277. doi: 10. 1016/S1474-4422(08)70042-5 PMID: 18275928
- Munger KL, Zhang SM, O’Reilly E, Hernan MA, Olek MJ, Willett WC, et al. Vitamin D intake and incidence of multiple sclerosis. Neurology. 2004; 62:60-65. doi: 10.1212/WNL.63.5.939 PMID: 14718698
- Mowry EM, Krupp LB, Milazzo M, Chabas D, Strober JB, Belman AL, et al. Vitamin D status is associated with relapse rate in pediatric-onset multiple sclerosis. Ann Neurol. 2010; 67:618-624. doi: 10. 1002/ana.21972 PMID: 20437559
- Mowry EM, Waubant E, McCulloch CE, Okuda DT, Evangelista AA, Lincoln RR, etal. Vitamin D status predicts new brain magnetic resonance imaging activity in multiple sclerosis. Ann Neurol. 2012; 72:234-240. doi: 10.1002/ana.23591 PMID: 22926855
- Rosen CJ, Adams JS, Bikle DD, Black DM, Demay MB, Manson JE, etal. The nonskeletal effects of vitamin D: an Endocrine Society scientific statement. Endocr Rev. 2012; 33:456-492. doi: 10.1210/er. 2012-1000 PMID: 22596255
- James E, Dobson R, Kuhle J, Baker D, Giovannoni G, Ramagopalan SV. The effect of vitamin D- related interventions on multiple sclerosis relapses: a meta-analysis. Mult Scler. 2013; 19:1571-1579. doi: 10.1177/1352458513489756 PMID: 23698130
- Burton JM, Kimball S, Vieth R, Bar-Or A, Dosch H-M, Cheung R, etal. A phase I/II dose-escalation trial of vitamin D3 and calcium in multiple sclerosis. Neurology. 2010; 74:1852-1859. doi: 10.1212/WNL. 0b013e3181e1cec2 PMID: 20427749
- Kimball SM, Ursell MR, O’Connor P, Vieth R. Safety of vitamin D3 in adults with multiple sclerosis. Am J Clin Nutr. 2007; 86:645-651. PMID: 17823429
- WrightS, Yadav V, BeverCJ, Bowen J, Bowling A, Weinstock-Guttman B, etal. Summary of evidence- based guideline: complementary and alternative medicine in multiple sclerosis: report of the Guideline Development Subcommittee of the American Academy of Neurology. Neurology. 2014; 83:1484-1486. doi: 10.1212/01.wnl.0000455935.13606.91
- Looker AC, PfeifferCM, Lacher DA, Schleicher RL, Picciano MF, Yetley EA. Serum 25-hydroxyvitamin D status of the US population: 1988-1994 compared with 2000-2004. Am J Clin Nutr. 2008; 88:15191527. doi: 10.3945/aicn.2008.26182 PMID: 19064511
- Yetley EA. Assessing the vitamin D status of the US population. Am J Clin Nutr. 2008; 88:558S-564S. PMID: 18689402
- Lawlor DA, Harbord RM, Sterne JAC, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008; 27:1133-1163. PMID: 17886233
- Voight BF, Peloso GM, Orho-Melander M, Frikke-Schmidt R, Barbalic M, Jensen MK, etal. Plasma HDL cholesterol and risk of myocardial infarction: a Mendelian randomisation study. Lancet. 2012; 380:572-580. doi: 10.1016/S0140-6736(12)60312-2 PMID: 22607825
- Zacho J, Tybjaerg-Hansen A, Jensen JS, Grande P, Sillesen H, Nordestgaard BG. Genetically elevated C-reactive protein and ischemic vascular disease. N Engl J Med. 2008; 359:1897-1908. doi: 10. 1056/NEJMoa0707402 PMID: 18971492
- Cohen JC, Boerwinkle E, Mosley THJ, Hobbs HH. Sequence variations in PCSK9, low LDL, and protection against coronary heart disease. N Engl J Med. 2006; 354:1264-1272. doi: 10.1056/ NEJMoa054013 PMID: 16554528
- WangTJ, Zhang F, Richards JB, Kestenbaum B, Van Meurs JB, Berry D, etal. Common genetic determinants of vitamin D insufficiency: a genome-wide association study. Lancet. 2010; 376:180-188. doi: 10.1016/S0140-6736(10)60588-0 PMID: 20541252
- Faye LL, Sun L, Dimitromanolakis A, Bull SB. A flexible genome-wide bootstrap method that accounts for ranking and threshold-selection bias in GWAS interpretation and replication study design. Stat Med. 2011; 30:1898-1912. doi: 10.1002/sim.4228 PMID: 21538984
- International Multiple Sclerosis Genetic Consortium, Beecham AH, Patsopoulos NA, Xifara DK, Davis MF, Kemppinen A, et al. Analysis of immune-related loci identifies 48 new susceptibility variants for multiple sclerosis. Nat Genet. 2013; 45:1353-1360. doi: 10.1038/ng.2770 PMID: 24076602
- Polman CH, Reingold SC, Edan G, Filippi M, Hartung H-P, Kappos L, etal. Diagnostic criteria for multiple sclerosis: 2005 revisions to the “McDonald Criteria.” Ann Neurol. 2005; 58:840-846. doi: 10.1002/ ana.20703 PMID: 16283615
- McDonald WI, Compston A, Edan G, Goodkin D, Hartung HP, Lublin FD, etal. Recommended diagnostic criteria for multiple sclerosis: guidelines from the international panel on the diagnosis of multiple sclerosis. Ann Neurol. 2001; 50:121-127. PMID: 11456302
- Poser CM, Paty DW, Scheinberg L, McDonald WI, Davis FA, Ebers GC, etal. New diagnostic criteria for multiple sclerosis: guidelines for research protocols. Ann Neurol. 1983; 13:227-231. doi: 10.1002/ ana.410130302 PMID: 6847134
- International Multiple Sclerosis Genetic Consortium, Wellcome Trust Case Control Consortium 2, Saw- cer S, Hellenthal G, Pirinen M, Spencer CCA, etal. Genetic risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis. Nature. 2011; 476:214-219. doi: 10.1038/nature10251 PMID: 21833088
- Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, Handsaker RE, etal. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012; 491:56-65. doi: 10.1038/nature11632 PMID: 23128226
- Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014; 23:R89-R98. doi: 10.1093/hmg/ddu328 PMID: 25064373
- Berry DJ, Vimaleswaran KS, Whittaker JC, Hingorani AD, Hypponen E. Evaluation of genetic markers as instruments for Mendelian randomization studies on vitamin D. PLoSONE. 2012; 7:e37465. doi: 10. 1371/iournal.pone.0037465 PMID: 22629401
- Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013; 37:658-665. doi: 10.1002/gepi.21758 PMID: 24114802
- Ehret GB, Munroe PB, Rice KM, Bochud M, Johnson AD, Chasman DI, etal. Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature. 2011; 478:103-109. doi: 10.1038/nature10405 PMID: 21909115
- Dastani Z, Hivert M-FF, Timpson NJ, Perry JRB, Yuan X, Scott RA., etal. Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals. PLoS Genet. 2012; 8:e1002607. doi: 10.1371/iournal.pgen.1002607 PMID: 22479202
- Patsopoulos NA, Evangelou E, Ioannidis JPA. Sensitivity of between-study heterogeneity in meta-analysis: proposed metrics and empirical evaluation. Int J Epidemiol. 2008; 37:1148-1157. doi: 10.1093/iie/ dyn065 PMID: 18424475
- Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. 2003; 327:557-560. PMID: 12958120
- Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med. 2002; 21:15391558. PMID: 12111919
- Ross AC. Dietary reference intakes for calcium and vitamin D. Washington (District of Columbia): Institute of Medicine. Available: http://www.iom.edu/Reports/2010/Dietary-Reference-Intakes-for-Calcium- and-Vitamin-D.aspx. Accessed 22 July 2015.
- Dastani Z, Li R, Richards B. Genetic regulation of vitamin D levels. Calcif Tissue Int. 2013; 92:106-117. doi: 10.1007/s00223-012-9660-z PMID: 23114382
- Safadi FF, Thornton P, Magiera H, Hollis BW, Gentile M, Haddad JG, et al. Osteopathy and resistance to vitamin D toxicity in mice null for vitamin D binding protein. J Clin Invest. 1999; 103:239-251. PMID: 9916136
- Bikle DD, Siiteri PK, Ryzen E, Haddad JG. Serum protein binding of 1,25-dihydroxyvitamin D: a reevaluation by direct measurement of free metabolite levels. J Clin Endocrinol Metab. 1985; 61:969-975. doi: 10.1210/jcem-61-5-969 PMID: 3840175
- Verboven C, Rabijns A, De Maeyer M, Van Baelen H, Bouillon R, De Ranter C. A structural basis for the unique binding features of the human vitamin D-binding protein. Nat Struct Biol. 2002; 9:131-136. doi: 10.1038/nsb754 PMID: 11799400
- Shinkyo R, Sakaki T, Kamakura M, Ohta M, Inouye K. Metabolism of vitamin D by human microsomal CYP2R1. Biochem Biophys Res Commun. 2004; 324:451-457. doi: 10.1016/i.bbrc.2004.09.073 PMID: 15465040
- KongsbakM, von Essen MR, Levring TB, Schjerling P, Woetmann A, Odum N, etal. Vitamin D-binding protein controls T cell responses to vitamin D. BMC Immunol. 2014; 15:35. doi: 10.1186/s12865-014- 0035-2 PMID: 25230725
- Moon M, Song H, Hong HJ, Nam DW, ChaM-Y, Oh MS, etal. Vitamin D-binding protein interacts with Ap and suppresses Ap-mediated pathology. Cell Death Differ. 2013; 20:630-638. doi: 10.1038/cdd. 2012.161 PMID: 23257976
- Rinaldi AO, Sanseverino I, PurificatoC, Cortese A, Mechelli R, Francisci S, etal. Increased circulating levels of vitamin D binding protein in MS patients. Toxins (Basel). 2015; 7:129-137. doi: 10.3390/ toxins7010129
- Yang M, Qin Z, Zhu Y, Li Y, Qin Y, Jing Y, etal. Vitamin D-binding protein in cerebrospinal fluid is associated with multiple sclerosis progression. Mol Neurobiol. 2013; 47:946-956. doi: 10.1007/s12035-012- 8387-1 PMID: 23339019
- Browne RW, Weinstock-Guttman B, Zivadinov R, Horakova D, BodziakML, Tamano-Blanco M, etal. Serum lipoprotein composition and vitamin D metabolite levels in clinically isolated syndromes: results from a multi-center study. J Steroid Biochem Mol Biol. 2014; 143:424-433. doi: 10.1016/i.isbmb.2014. 06.007 PMID: 24950029
- WillerCJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, etal. Discovery and refinement of loci associated with lipid levels. Nat Genet. 2013; 45:1274-1283. doi: 10.1038/ng.2797 PMID: 24097068
- Van Luijn MM, Kreft KL, JongsmaML, Mes SW, Wierenga-Wolf AF, van Meurs M, etal. Multiple sclerosis-associated CLEC16A controls HLA class II expression via late endosome biogenesis. Brain. 2015; 138:1531-1547. doi: 10.1093/brain/awv080 PMID: 25823473
- Ramagopalan SV, Maugeri NJ, Handunnetthi L, Lincoln MR, Orton S-M, Dyment DA, etal. Expression of the multiple sclerosis-associated MHC class II Allele HLA-DRB1*1501 is regulated by vitamin D. PLoS Genet. 2009; 5:e1000369. doi: 10.1371/iournal.pgen.1000369 PMID: 19197344
- O’Connell K, Kelly S, Kinsella K, Jordan S, Kenny O, Murphy D, et al. Dose-related effects of vitamin D on immune responses in patients with clinically isolated syndrome and healthy control participants: study protocol for an exploratory randomized double- blind placebo-controlled trial. Trials. 2013;
- 14:272. doi: 10.1186/1745-6215-14-272 PMID: 23981773
- Thouvenot E, Suehs CM. Efficacy of cholecalciferol (vitamin D3) for delaying the diagnosis of MS after a clinically isolated syndrome (D-Lay-MS). ClinicalTrials.gov. Available: https://clinicaltrials.gov/ct2/ show/NCT01817166?term=vitamin+d+multiple+sclerosis&rank=13. Accessed 22 July 2015.
- Mokry LE, Ahmad O, Forgetta V, Thanassoulis G, Richards JB. Mendelian randomisation applied to drug development in cardiovascular disease: a review. J Med Genet. 2015; 52:71-79. doi: 10.1136/ imedgenet-2014-102438 PMID: 25515070
- Smith GD, Ebrahim S. “Mendelian randomization”: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003; 32:1-22. PMID: 12689998
- Leong A, Rehman W, Dastani Z, Greenwood C, Timpson N, Langsetmo L, et al. The causal effect of vitamin D binding protein (DBP) levels on calcemic and cardiometabolic diseases: a Mendelian randomization study. PLoS Med. 2014; 11:e1001751. doi: 10.1371/iournal.pmed.1001751 PMID: 25350643
- Ye Z, Sharp SJ, Burgess S, Scott RA, Imamura F, Consortium I, etal. Association between circulating 25-hydroxyvitamin D and incident type 2 diabetes: a Mendelian randomisation study. Lancet Diabetes Endocrinol. 2015; 3:35-42. doi: 10.1016/s2213-8587(14)70184-6 PMID: 25281353
- Cooper JD, Smyth DJ, Walker NM, Stevens H, Burren OS, Wallace C, etal. Inherited variation in vitamin D genes is associated with predisposition to autoimmune disease type 1 diabetes. Diabetes. 2011; 60:1624-1631. doi: 10.2337/db10-1656 PMID: 21441443
- Vimaleswaran KS, Cavadino A, Berry DJ, LifeLines Cohort Study Investigators, Jorde R, Dieffenbach AK, etal. Association of vitamin D status with arterial blood pressure and hypertension risk: a Mendelian randomisation study. Lancet Diabetes Endocrinol. 2014; 2:719-729. doi: 10.1016/s2213-8587(14) 70113-5 PMID: 24974252
- Vimaleswaran KS, Berry DJ, Lu C, Tikkanen E, Pilz S, Hiraki LT, etal. Causal relationship between
- obesity and vitamin D status: bi-directional Mendelian randomization analysis of multiple cohorts. PLoS Med. 2013; 10:e1001383. doi: 10.1371/iournal.pmed.1001383 PMID: 23393431
This list of resources contains links that can be accessed when viewing the PDF on a device or via the online version of the article at http://dx.doi.org/10.1371/iournal.pmed.1001866.
- The UK National Health Service Choices website provides information about multiple sclerosis, including personal stories, and information about vitamin D and how to get vitamin D from sunlight; A “Behind the Headlines” article examines health claims made for vitamin D
- The US National Institute of Neurological Disorders and Stroke provides information about multiple sclerosis and links to MS societies
- The National Multiple Sclerosis Society, a US not-for-profit organization, provides information about all aspects of MS, including information about vitamin D and MS and personal stories about MS
- The Multiple Sclerosis Society, a UK not-for-profit organization, also provides information about MS for patients, carers, and professionals in several languages, including news about the latest research, personal blogs, and forums
- The US Office of Dietary Supplements provides information about vitamin D (in English and Spanish)
- MedlinePlus has links to further information about multiple sclerosis and vitamin D (in English and Spanish)
- Wikipedia has a page on Mendelian randomization (note: Wikipedia is a free online encyclopedia that anyone can edit; available in several languages)
Dr Benjamin Jacobs, Consultant Paediatrician and Director of Children’s Service at the Royal National Orthopaedic Hospital, said:
“This study reveals important new evidence of a link between Vitamin D Deficiency and Multiple sclerosis (MS). Data from two massive genetic studies were combined: The European study of 34,000 people’s genes associated with Vitamin D Deficiency (the SUNLIGHT study); and the IMSGC Immunochip study, the largest international genetic study of MS, involving 14,498 MS cases and 24,091 healthy controls.
“The results show that if a baby is born with genes associated with Vitamin D Deficiency they are twice as likely as other babies to develop MS as an adult. This could be because Vitamin D deficiency causes MS or possibly because there are other complex genetic interactions.
We do not yet know if giving healthy children and adults Vitamin D will decrease their risk of developing MS, but clinical trials are being conducted now to study this.”
Prof. Danny Altmann, Professor of Immunology at Imperial College London, said:
“While MS has been understood for several decades to be a neuroinflammatory disease with a cocktail of contributory causes both from genes and the environment (‘nature and nurture’), we have tended to acquire more and more data without perhaps crystallising a clear notion of what this means about the origin of disease or how to prevent it.
“The ‘environmental’ component of susceptibility certainly encompasses many factors, however, considerable attention has focused on individual differences in exposure to sunlight and thus, levels of vitamin D. This vitamin is vital to many biochemical pathways in the body, notably including diverse aspects of efficient function of immunity. It has long been appreciated that there is a geographical correlation between living at latitudes with reduced sunlight and susceptibility to MS.
“What this keynote study by Richards and colleagues offers is an elegant utilisation of data from very large human cohorts to show that there are parts of the human genome that predispose to high or low levels of vitamin D in the blood, this is in turn being related to risk of developing MS. Until now many of the biggest hits from MS genetic studies have implicated genes controlling immunity, but this study draws our attention specifically to the relationship between exposure to sunlight, immune function and MS risk. This places MS alongside other diseases of immune susceptibility, including many infectious diseases, where risk is associated with low levels of serum vitamin D.
“In many of the other examples tried, the notion that all will be well if we simply give patients a large dose of vitamin D supplements has proved too simple. However, vitamin D is relatively cheap, safe and many of us would be all the healthier if we could achieve the serum levels that our ancient ancestors presumably acquired when roaming outdoors in temperate climates, unclothed and eating a diverse diet including oily fish. While it may be too much to expect therapeutic vitamin D to treat or reverse ongoing MS, this paper will add to the weight of argument for routine vitamin D supplementation of foodstuffs as a broad, preventative, public health measure.”