Genetically-Proxied Levels of Vitamin D and Risk of Intracerebral Hemorrhage.

Background The evidence linking vitamin D (VitD) levels and spontaneous intracerebral hemorrhage (ICH) remains inconclusive. We tested the hypothesis that lower genetically determined VitD levels are associated with higher risk of ICH. Methods and Results We conducted a 2 sample Mendelian Randomization (MR) study using publicly available summary statistics from published genome-wide association studies of VitD levels (417 580 study participants) and ICH (1545 ICH cases and 1481 matched controls). We used the inverse-variance weighted approach to generate causal estimates and the MR Pleiotropy Residual Sum and Outlier and MR-Egger approaches to assess for horizontal pleiotropy. To account for known differences in their underlying mechanism, we implemented stratified analysis based on the location of the hemorrhage within the brain (lobar or nonlobar). Our primary analysis indicated that each SD decrease in genetically instrumented VitD levels was associated with a 60% increased risk of ICH (odds ratio [OR], 1.60; [95% CI, 1.05-2.43]; P=0.029). We found no evidence of horizontal pleiotropy (MR-Egger intercept and MR Pleiotropy Residual Sum and Outlier global test with P>0.05). Stratified analyses indicated that the association was stronger for nonlobar ICH (OR, 1.87; [95% CI, 1.18-2.97]; P=0.007) compared with lobar ICH (OR, 1.43; [95% CI, 0.86-2.38]; P=0.17). Conclusions Lower levels of genetically proxied VitD levels are associated with higher ICH risk. These results provide evidence for a causal role of VitD metabolism in ICH.

Vitamin D (VitD) metabolism is an appealing target for preventive interventions in cerebrovascular disease, as it can be easily and safely intervened upon. A number of studies have explored the role of VitD in ischemic stroke. , Two observational studies demonstrated that depleted levels of VitD were associated with increased mortality and worse outcome in this condition. , Additionally, Zhou et al found that lower VitD levels were related to a higher risk of ischemic stroke in a large meta‐analysis. Genetic studies have also examined the relationship between VitD levels and ischemic stroke. A Mendelian Randomization (MR) study evaluating the association between genetically determined VitD and ischemic stroke did not find an association between these 2 traits. Despite the significant amount of research outlined above, the role of VitD metabolism in intracerebral hemorrhage (ICH) remains relatively unexplored. While ICH shares some biological features with ischemic stroke, the overall pathophysiology of these 2 types of cerebrovascular diseases is significantly different. We therefore aimed to use MR analyses to evaluate the role of circulating levels of VitD in ICH, hypothesizing that genetically reduced levels of VitD lead to higher risk of ICH.

METHODS

Because of the sensitive nature of the data collected for this study, requests to access the data set from qualified researchers trained in human subject confidentiality protocols may be sent to dbGap. Institutional Review Board approval was not required for this study.

Study Design

We conducted a 2‐sample MR study, a study design where the genetic information for the exposure (VitD levels) and the outcome (risk of ICH) come from different studies, using summary statistics from large genetic studies of these traits. Importantly, there was no overlap between both study populations.

Instruments

We used independent (r 2>0.1) single nucleotide polymorphisms (SNPs) associated with VitD levels at genome‐wide levels (P<5×10−8) , from the largest GWAS (Genome‐Wide Association Study) on VitD available to date. Palindromic SNPs were excluded from all analyses. For the final list of VitD associated instruments, we abstracted effect estimates and standard errors from the largest GWAS of ICH conducted to date, using proxy variants (r 2>0.9) where available. The average F‐statistic for the genetic instrument was 216 and the median was 39. The I2 for the genetic instrument was 0.99. In combination, these metrics suggest that the genetic instrument had sufficient strength. We did not adjust for ancestry in the MR analysis because both GWAS (for the exposure and for the outcome) conducted principal component analyses, a robust method to identify genetic ancestry, to identify study participants from European ancestry and exclude population outliers.

Statistical Analysis

Our primary MR analysis used the inverse‐variance weighted method. We tested for horizontal pleiotropy (the possibility that the effect of the instrument on the outcome of interest is exerted through a pathway other than VitD levels) using the MR‐Egger intercept and Mendelian Randomization Pleiotropy Residual Sum and Outlier (MR‐PRESSO) global test. In sensitivity analyses, we implemented the weighted median method and calculated the MR‐PRESSO outlier corrected estimate.

We also examined whether any of the instruments were potentially associated with known confounders of the association between VitD and ICH. To this end, we performed sensitivity analyses using Cook's distance, defining outlier instruments as those with a Cook's distance >4/number of SNPs. We also looked for evidence of SNP outliers with RadialMR. In addition, we performed sensitivity analyses using multivariable MR to adjust for the genetic contribution of VitD going through other pathways associated with ICH. Finally, we completed leave‐one‐out plots (Figures S1 through S3). We recalibrated the effect estimates to express the change in ICH risk associated with an SD decrease in genetically instrumented VitD levels. We declared statistical significance at P<0.05 (2‐tailed) when testing the single primary hypothesis that lower genetically determined VitD levels are associated with higher risk of ICH. We used R (version 3.6.0) and its TwoSampleMR and MR‐PRESSO packages.

Results

From a GWAS of VitD that analyzed data on 417 580 people of European ancestry, we identified 108 SNPs that met our criteria to select instruments for this trait (Table 1). Of these, 88 SNPs (Table S1) had available summary results in a GWAS of ICH that included 1545 ICH cases (664 lobar and 881 nonlobar ICH cases) and 1481 matched controls of European ancestry (Table 1). Our primary analysis using the inverse‐variance weighted MR method indicated that each SD decrease in genetically instrumented VitD levels was associated with an 60% increased risk of ICH (odds ratio [OR], 1.60; [95% CI, 1.05–2.43]; P=0.029) (Figure). To rule out heterogeneity in estimates across SNPs, we calculated Cochran Q statistics, which did not reveal heterogeneity. For all ICH, Q was 105.01, with 87 degrees of freedom and P value of 0.09. For lobar ICH, Q was 89.43 with 87 degrees of freedom and P value of 0.41. For nonlobar ICH, Q was 88.77, with 88 degrees of freedom and P value of 0.46.

Table 1

Characteristics of GWAS of Vitamin D and Intracerebral Hemorrhage

Characteristic	GWAS of vitamin D	GWAS of intracerebral hemorrhage
Sample size	417 580	3026
Mean age, y	47–66 [0.34–14.12]	67 [SD, 10]
Female sex	51%	45%
Genotyping platform	UK Biobank Axiom Array	Affymetrix 6.0 Illumina HumanHap610
No. of SNPs evaluated	8 806 780	5 258 103
No. of genome‐wide significant loci	143	2

GWAS indicates Genome‐Wide Association Study; and SNP, single nucleotide polymorphism.

Figure Figure. 

Mendelian randomization plot.

The plot presents the effect estimates of association tests between the single nucleotide polymorphisms and vitamin D levels (X axis) and risk of intracerebral hemorrhage (Y axis). The blue summary line corresponds to the slope of the inverse‐variance weighted method. ICH indicates intracerebral hemorrhage.

There was no evidence of significant horizontal pleiotropy, as evidenced by the nonsignificant results of the MR‐Egger intercept and MR‐PRESSO Global tests (both P>0.05, Table 2). The direction and size of the effect were consistent in secondary analyses aimed at reducing the impact of any pleiotropy (Table 2). Summary statistics for analyses stratified by location were based on the analysis of 664 lobar and 881 nonlobar ICH cases (Table 1). In these stratified analyses according to location (Tables S2 and S3), the results remained significant for nonlobar ICH (OR, 1.87; [95% CI, 1.18–2.97]; P=0.007) but not for lobar ICH (OR, 1.43; [95% CI, 0.86–2.38]; P=0.17) (Table 2). For nonlobar ICH, there was no evidence of significant horizontal pleiotropy (P>0.05 for both the MR‐Egger intercept and MR‐PRESSO Global tests) and the direction and size of the effect were consistent in secondary analyses aimed at reducing pleiotropy (Table 2). The direction of effect for the MR‐Egger estimate was consistent with other methods. For all ICH, the MR‐Egger estimate was 1.42 (95% CI, 0.82–2.46; P=0.21). For nonlobar ICH, the MR‐Egger estimate was 1.51 (95% CI, 0.83–2.73; P=0.19). For lobar ICH, the MR‐Egger estimate was 1.21 (95% CI, 0.62–2.36; P=0.57). After performing Cook’s distance, we found 2 outliers for all ICH, 1 outlier for lobar ICH and 3 outliers for nonlobar ICH. The Cook’s distance plots can be found in Figures S4 through S6. After removing these outliers, the association remained significant for both all ICH (OR, 1.95; [95% CI, 1.08–3.55]; P=0.027) and nonlobar ICH (OR, 2.37; [95% CI, 1.20–4.67]; P=0.013).

Table 2

Results of Different Mendelian Randomization Analyses

Mendelian randomization method	All intracerebral hemorrhage		Lobar intracerebral hemorrhage		Nonlobar intracerebral hemorrhage
	OR (95% CI)	P value	OR (95% CI)	P value	OR (95% CI)	P value
Association tests
Inverse‐variance weighted	1.60 (1.05–2.43)	0.029	1.43 (0.86–2.38)	0.17	1.87 (1.18–2.97)	0.007
Weighted median	1.42 (0.81–2.50)	0.22	1.36 (0.66–2.79)	0.40	1.46 (0.76–2.80)	0.26
MR‐Egger estimate	1.42 (0.82–2.46)	0.21	1.21 (0.62–2.36)	0.57	1.51 (0.83–2.73)	0.19
Weighted mode	1.37 (0.89–2.11)	0.15	1.11 (0.62–1.98)	0.73	1.45 (0.84–2.52)	0.19
Tests for horizontal pleiotropy*
MR‐Egger intercept	1.01 (0.99–1.02)	0.52	1.01 (0.99–1.03)	0.45	1.01 (0.99–1.03)	0.26
MR‐PRESSO global test	…	0.11	…	0.44	…	0.48

IVW indicates inverse‐variance weighted; MR, Mendelian Randomization; MR‐PRESSO, Mendelian Randomization Pleiotropy Residual Sum and Outlier; and OR, odds ratio.

Horizontal pleiotropy is absent if P>0.05.

In additional sensitivity analyses using multivariable MR to adjust for the genetic contribution of VitD to other pathways associated with ICH, we found that the association between VitD and both all and nonlobar ICH remained significant (Table S4). Finally, we have conducted a sensitivity analysis using an R 2 threshold of 0.001 for clumping, that resulted in the selection of 64 SNPs. Results were similar and consistent with the original analyses for all ICH (OR, 1.55; [95% C,I 1.01–2.40]; P=0.047), lobar ICH (OR, 1.47; [95% CI, 0.86–2.54]; P=0.16), and nonlobar ICH (OR, 1.78; [95% CI, 1.08–2.92]; P=0.02). Moreover, when testing for potential outliers with RadialMR, we did not detect any significant outliers for all ICH, lobar ICH, or nonlobar ICH using a Bonferroni‐corrected threshold. Using a nominal threshold of 0.05, there were 9 outliers for all ICH, 6 outliers for lobar ICH, and 4 outliers for nonlobar ICH. The inverse‐variance weighted MR remained significant for both all ICH (P=0.0467) and nonlobar ICH (P=0.016) after removing these outliers.

Discussion

We report the results of a 2‐sample MR study that used summary level results to evaluate the relationship between genetically determined VitD levels and risk of ICH. We found that genetically reduced levels of VitD associated with a modest increase in the risk of ICH, without evidence of significant horizontal pleiotropy and consistent direction of effects in secondary, more conservative analyses. However, only substantial pleiotropic effects can be ruled out because of the limited statistical power of methods such as MR‐Egger. Stratification by location indicated that these associations were predominantly driven by nonlobar ICH, where associations were stronger despite a significant reduction in sample size.

Prior studies investigated the link between VitD deficit and ischemic stroke. A systematic review and meta‐analysis demonstrated a 62% higher risk of stroke among individuals in the lowest versus highest categories of VitD concentrations. Importantly, a randomized controlled trial evaluating cardiovascular disease broadly failed to show a protective effect of VitD supplementation on stroke risk. From a population genetics perspective, a prior MR study that evaluated the relationship between genetically instrumented VitD levels and risk of ischemic stroke failed to find a protective association.

Our study adds important new evidence to the field of VitD metabolism and ICH, an area of research that remains relatively unexplored. A prior MR study examined causal associations between VitD levels and several cardiovascular traits, including myocardial infarction, ischemic heart disease, ischemic stroke, subarachnoid hemorrhage, and ICH among Chinese. Importantly, the effect size for the association between genetically determined VitD levels and ICH risk observed in this study was similar (hazard ratio, 1.09 [95% CI, 1.01–1.18] per 25 nmol/L higher plasma 25(OH)D) to that obtained in our study. Another study implemented MR analyses using data from the China Kadoorie Biobank, the Copenhagen City Heart Study, and the Copenhagen General Population Study, and did not find any significant associations. However, the study only considered VitD‐related variants in the CYP2R1 and DCHR7 genes. In contrast, the present study used a more powerful approach using several dozen genetic risk variants known to modify VitD levels. In addition, by focusing on a single outcome, we increase our discovery power by reducing multiple testing.

There are several possible pathophysiological mechanisms that could mediate the inverse association between VitD and ICH. First, VitD plays an important role in the occurrence of various cardiometabolic traits that predispose to ICH, including diabetes, metabolic syndrome, and renin‐angiotensin system activation. Another potential mechanism involves an upregulation of inflammatory pathways, as VitD inhibits the production of various inflammation factors.

In conclusion, our results point to an inverse association between genetically determined VitD levels and risk of ICH. This association was driven by nonlobar (also known as deep) hemorrhages. Further research is needed to confirm these findings, extend these results to non‐Europeans, and identify the putative mechanisms behind this association.

Sources of Funding

Acosta is supported by the American Heart Association (AHA) Bugher Fellowship in Hemorrhagic Stroke Research. Leasure is supported by the AHA Medical Student Research Fellowship. Gill is supported by the Yale Claude D. Pepper Older Americans Independence Center and by the National Institute on Aging (P30AG021342). Sheth is supported by the National Institutes of Health and the AHA, and reports grants from Hyperfine, Biogen, and Bard unrelated to this work. Falcone is supported by the National Institutes of Health and the AHA.

Disclosures

Sheth received compensation from Sense, Cerevasc, Certus, Rhaeos, Zoll Medical Corporation, CSL Behring, Astrocyte, and a patent for Alva Health. The remaining authors have no disclosures to report.

Tables S1–S4

Figures S1–S6

Click here for additional data file.