Migraine and Early-onset Ischemic Stroke: A Mendelian Randomization Study
=========================================================================

* Peng-Peng Niu
* Rui Zhang
* Shuo Li
* Yu-Sheng Li

## Abstract

**Objective** Through the utilization of the data specifically related to early-onset ischemic stroke, we aimed to investigate the causal effects of migraine and its subtypes on the risk of early-onset ischemic stroke using the two-sample Mendelian randomization method.

**Methods** Genetic instrumental variables were acquired from two sources with the largest sample sizes available. Summary data for early-onset ischemic stroke was acquired from a study encompassing individuals aged 18 to 59 years, comprising 16,730 cases and 599,237 non-stroke controls. The random-effects inverse variance weighted method was used as the primary analysis approach.

**Results** The Mendelian randomization analysis revealed no association between overall migraine and migraine without aura with the risk of early-onset ischemic stroke. However, migraine with aura showed a suggestive association with an elevated risk of early-onset ischemic stroke, with odds ratios of 1.114 (95% confidence interval = 1.005 to 1.236, p-value = 0.040) and 1.062 (95% confidence interval = 1.002 to 1.126, p-value = 0.042) based on instruments from two independent sources. The odds ratio was 1.074 (95% confidence interval = 1.022 to 1.130, p-value = 0.005) based on instruments from both two sources. No evidence of heterogeneity or horizontal pleiotropy was found. Furthermore, a positive genetic correlation was found between migraine with aura and early-onset ischemic stroke (genetic correlation = 0.208, 95% confidence interval = 0.038 to 0.377, p-value = 0.016). By contrast, migraine with aura was not related to ischemic stroke of all adults.

**Conclusion** This study provides evidence of a causal relationship between migraine with aura and the risk of early-onset ischemic stroke.

Keywords
*   Causality
*   genetic correlation
*   ischemic stroke
*   Mendelian randomization
*   migraine

## Introduction

Migraine, a prevalent neurological disorder, is characterized by recurrent and often severe headaches, frequently accompanied by a range of associated symptoms. It has a global age-standardized prevalence of approximately 14%.1 Observational studies have shown that there was an increased risk of ischemic stroke for individuals who are afflicted with migraine, especially for those with migraine with aura. However, the preceding Mendelian randomization (MR) studies have yet to succeed in confirming the causal connection between them.2, 3 The lack of conclusive MR findings may be due to migraine’s potential association with early-onset rather than late-onset ischemic stroke.4

We employed the two-sample MR method to assess the causal effects of migraine and its subtypes on early-onset ischemic stroke (Figure 1).

![Figure 1](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/02/2024.08.01.24311334/F1.medium.gif)

[Figure 1](http://medrxiv.org/content/early/2024/08/02/2024.08.01.24311334/F1)

Figure 1 Design of the present study
In MR study, three key assumptions are crucial: the relevance assumption, stating that the genetic variant strongly associates with the exposure; the independence assumption, ensuring the genetic variant’s independence from confounders; and the exclusion restriction assumption, which dictates that the genetic variant solely affects the outcome through the exposure factor.

For instrumental SNPs used in both data sources, we first combined un-clumped SNPs from the two sources and then conducted the clumping procedure.

**Abbreviations**: Any, any migraine; MA, migraine with aura; MO, migraine without aura; SNP, single nucleotide polymorphism; IHGC, International Headache Genetics Consortium.

## Methods

### Summary data for migraine

Three sets of instrumental SNPs were used. The first set of instrumental SNPs were chosen from the genome-wide association study (GWAS) performed by Bjornsdottir et al recently (Supplementary Table 1).5 The second set of instrumental SNPs were chosen from two GWASs performed by the International Headache Genetics Consortium (IHGC) (Supplementary Table 2).6, 7 The third set of instrumental SNPs were obtained by combining the SNPs from both the two sources together to carry out the clumping procedure. Significant variants at the genome-wide level (p-value < 5 × 10-8) were utilized as instrumental variables for any migraine. As there is a limited number of variants with a p-value < 5 × 10-8 for migraine subtypes, similar to the previous MR studies associated with migraine subtypes, the threshold for instrumental variables was relaxed to < 1 × 10-5.2, 3, 8 The clumping procedure was carried out by employing an r2 value of 0.01 and a physical window of 10 MB. Supplementary Figure 1 shows the flow diagram of instrumental variables selection. SNPs associated with cardiometabolic risk factors were excluded (Supplementary Table 3 and Supplementary Table 4).9

### Summary data for early-onset ischemic stroke

Jaworek et al conducted a GWAS meta-analysis on early-onset ischemic stroke, with individuals ranging in age from 18 to 59 years.10 We employed summary data specific to Europeans, which encompassed 11,114 cases and 435,540 controls, to carry out the MR analysis. To make a comparison, we also carried out MR analyses by utilizing the summary data of ischemic stroke which is not specifically for early onset.11 If instrumental SNPs are not found in the outcome data, proxy SNPs (with an r2 > 0.8) will be used if available.

### Mendelian Randomization Analysis

We employed the random-effects inverse variance weighted method as the primary analysis approach. We also utilized MR-Egger, weighted median, and mode-based MR methods to account for the impact of horizontal pleiotropy in MR studies. Additionally, we employed the MR Pleiotropy RESidual Sum and Outlier approach to identify potential outliers.

We carried out MR analysis by utilizing the TwoSampleMR R package within R (version 3.6.1).12 We estimated the genetic correlation between migraine and ischemic stroke using the LD score software (LDSC v1.0.1).13

## Results

Supplementary Table 5 show the characteristics of SNP included. The MR analysis demonstrated that overall migraine and migraine without aura were not related to the risk of early-onset ischemic stroke (Figure 2 and Supplementary Table 6).

![Figure 2](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/02/2024.08.01.24311334/F2.medium.gif)

[Figure 2](http://medrxiv.org/content/early/2024/08/02/2024.08.01.24311334/F2)

Figure 2 Inverse variance weighted estimates of migraine on risk of ischemic stroke
Het_p indicates p-value of heterogeneity test.

Ple_p indicates p-value of pleiotropy test.

**Abbreviations**: SNP, single nucleotide polymorphism; IHGC, International Headache Genetics Consortium; OR, odds ratio; CI, confidence interval.

For migraine with aura, the odds ratios (OR) based on instrumental SNPs from Bjornsdottir et al., the International Headache Genetics Consortium, and both sources were 1.114 (95% confidence interval [CI] = 1.005 to 1.236, p-value = 0.040), 1.062 (95% CI = 1.002 to 1.126, p-value = 0.042), and 1.074 (95% CI = 1.022 to 1.130, p-value = 0.005), respectively (Figure 2, Supplementary Figure 2, Supplementary Figure 3, and Supplementary Table 6). There was no evidence of heterogeneity or horizontal pleiotropy. The MR Pleiotropy RESidual Sum and Outlier method did not detect any outliers. Leave-one-out analysis demonstrated that the results were not driven by any single SNP (Supplementary Figure 4). The Steiger test did not find reverse causal SNPs. Sensitivity analyses showed similar results (Supplementary Figure 5, Supplementary Table 7, and Supplementary Table 8). Lastly, migraine with aura was not related to the risk of all adult ischemic strokes (Figure 2).

Due to the limited sample size, MR analysis did not find consistent associations for early-onset ischemic stroke subtypes (Supplementary Figure 6). Reverse MR analysis showed that early-onset ischemic stroke was not associated with the risk of any migraine and its subtypes (Supplementary Figure 7).

Linkage disequilibrium score regression analysis showed a positive genetic correlation between migraine with aura and early-onset ischemic stroke (rg = 0.208, se = 0.086, p-value = 0.016) and a negative genetic correlation between migraine without aura and early-onset ischemic stroke (rg = -0.197, se = 0.099, p-value = 0.047). The genetic correlation between migraine with aura and ischemic stroke of all adults attenuated to 0.154 (se = 0.064, p-value = 0.015).

## Discussion

Previous studies suggested that migraine, especially migraine with aura, might be associated with ischemic stroke in young adults, but not in older.4 The age-related differential risk might be explained by the fact that, migraine and age-related risk factors such as hypertension might have shared underlying mechanisms. In older individuals, age-related risk factors condition may saturate the biological mechanisms that ultimately lead to ischemic stroke, so that there are few open pathways for migraine to affect ischemic stroke.14

There were several possible underlying mechanisms for the association between migraine and ischemic stroke. First, previous studies have revealed that migraine shares genetic basis with cardiometabolic risk factors and thrombosis.10 On the other hand, migraine with aura shows a higher prevalence of patent foramen ovale compared to migraine without aura, potentially elevating the risk of paradoxical embolism. Third, spreading depolarizations likely play a pivotal role in the increased risk of ischemic stroke in individuals with migraine, especially those with migraine accompanied by aura. Last, additional pathophysiological mechanisms might involve inflammatory markers, endothelial and vascular indicators, as well as sex hormones.15

The present study has several limitations. First, due to the unavailability of sex-specific summary data, we were unable to carry out sex-stratified analysis. Second, since only a few SNPs associated with migraine with aura reach genome-wide significance, we employed a more lenient threshold for instrumental variables. Third, due to the limited sample size of stroke subtypes, MR analysis did not identify consistent associations among the three migraine data sources. Lastly, the conclusions drawn from this study rely on summary data from European populations, potentially limiting their generalizability to other ethnic groups.

## Conclusion

In conclusion, this study confirmed the causal effect of migraine with aura on early-onset ischemic stroke risk. Further studies are needed to confirm the findings and to address the limitations of this study.

## Availability of data

The supporting data for this study can be found in the supplementary files of the article as well as in the referenced studies within this article.

Summary data for migraine by Bjornsdottir et al. were obtained from [https://www.decode.com/summarydata/](https://www.decode.com/summarydata/). Summary data for migraine by International Headache Genetics Consortium were obtained from [https://www.headachegenetics.org/datasets-cohorts](https://www.headachegenetics.org/datasets-cohorts). Summary data for early-onset ischemic stroke were obtained from the Cerebrovascular Disease Knowledge Portal ([https://cd.hugeamp.org/](https://cd.hugeamp.org/)).

## Declaration of interest

We declare no competing interests.

## Supporting information

Supplementary Table [[supplements/311334_file02.xlsx]](pending:yes)

## Data Availability

The supporting data for this study can be found in the supplementary files of the article as well as in the referenced studies within this article. Summary data for migraine by Bjornsdottir et al. were obtained from https://www.decode.com/summarydata/. Summary data for migraine by International Headache Genetics Consortium were obtained from https://www.headachegenetics.org/datasets-cohorts. Summary data for early-onset ischemic stroke were obtained from the Cerebrovascular Disease Knowledge Portal (https://cd.hugeamp.org/).

[https://www.decode.com/summarydata/](https://www.decode.com/summarydata/) 

[https://www.headachegenetics.org/datasets-cohorts](https://www.headachegenetics.org/datasets-cohorts) 

[https://cd.hugeamp.org/](https://cd.hugeamp.org/) 

## Funding

None.

![Supplementary Figure 1](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/02/2024.08.01.24311334/F3.medium.gif)

[Supplementary Figure 1](http://medrxiv.org/content/early/2024/08/02/2024.08.01.24311334/F3)

Supplementary Figure 1 Flow diagram of instrumental variables selection
Significant single nucleotide polymorphisms (SNPs) at the genome-wide level (with a p-value < 5 × 10-8) were utilized as instrumental variables for any migraine. However, due to a limited number of SNPs with a p-value < 5 × 10-8 for migraine with aura and migraine without aura, the threshold for instrumental variables was relaxed to < 1 × 10-5.1, 2

SNPs that are absent in the 1000G v3 reference panel of the European population were excluded before being screened.

Heterogeneity was determined to exist when the I2 exceeded 50% or the p-value of the Q statistic was less than 0.05.

Palindromic SNPs were excluded since the effect allele frequency was missing in the outcome data.

For the instrumental SNPs from both two sources, they were obtained through clumping on all the SNPs detailed in the box which is next to the last one.

All instrumental SNPs exhibit an F statistics value greater than 10, and this was arrived at by computing beta2 divided by se2.

References

1. Lee KJ, Lee SJ, Bae HJ, Sung J. Exploring the causal inference of migraine on stroke: A Mendelian randomization study. *Eur J Neurol* 2022;29(1):335-338. doi: 10.1111/ene.15101

2. Wu XP, Niu PP, Liu H. Association between migraine and venous thromboembolism: a Mendelian randomization and genetic correlation study. *Front Genet* 2024;15:1272599. doi: 10.3389/fgene.2024.1272599

![Supplementary Figure 2](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/02/2024.08.01.24311334/F4.medium.gif)

[Supplementary Figure 2](http://medrxiv.org/content/early/2024/08/02/2024.08.01.24311334/F4)

Supplementary Figure 2 Forest plots of the effect of migraine with aura on early-onset ischemic stroke
The horizontal axis represents the estimate of migraine with aura on early-onset ischemic stroke. A: instrumental SNPs based on study by Bjornsdottir et al. 2023

B: instrumental SNPs based on study by IHGC

C: instrumental SNPs based on both two sources.

![Supplementary Figure 3](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/02/2024.08.01.24311334/F5.medium.gif)

[Supplementary Figure 3](http://medrxiv.org/content/early/2024/08/02/2024.08.01.24311334/F5)

Supplementary Figure 3 Scatter plots of SNP effect on migraine with aura and SNP effect on early-onset ischemic stroke
A: instrumental SNPs based on study by Bjornsdottir et al. 2023;

B: instrumental SNPs based on study by IHGC;

C: instrumental SNPs based on both two sources.

Abbreviations: SNP, single nucleotide polymorphism; IHGC, International Headache Genetics Consortium.

![Supplementary Figure 4](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/02/2024.08.01.24311334/F6.medium.gif)

[Supplementary Figure 4](http://medrxiv.org/content/early/2024/08/02/2024.08.01.24311334/F6)

Supplementary Figure 4 Leave-one-out sensitivity analysis of the effect of migraine with aura on early-onset ischemic stroke
A: instrumental SNPs based on study by Bjornsdottir et al. 2023

B: instrumental SNPs based on study by IHGC

C: instrumental SNPs based on both two sources.

![Supplementary Figure 5](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/02/2024.08.01.24311334/F7.medium.gif)

[Supplementary Figure 5](http://medrxiv.org/content/early/2024/08/02/2024.08.01.24311334/F7)

Supplementary Figure 5 Sensitivity analyses of migraine with aura on early-onset ischemic stroke
A: Sensitivity analysis was conducted by exclusively relying on SNPs that exhibited low evidence of heterogeneity across different stroke cohorts (with I2 less than 50% and a p-value of the Q statistic greater than 0.05)

B: Sensitivity analysis was conducted using SNPs, without consideration of whether they were associated with potential confounders.

Het\_p indicates p-value of heterogeneity test. Ple\_p indicates p-value of pleiotropy test.

Abbreviations: SNP, single nucleotide polymorphism; IHGC, International Headache Genetics Consortium; OR, odds ratio; CI, confidence interval; MA, migraine with aura.

![Supplementary Figure 6](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/02/2024.08.01.24311334/F8.medium.gif)

[Supplementary Figure 6](http://medrxiv.org/content/early/2024/08/02/2024.08.01.24311334/F8)

Supplementary Figure 6 Inverse variance weighted estimates of migraine on risk of early-onset ischemic stroke subtypes
Abbreviations: MR, IHGC, International Headache Genetics Consortium; OR, odds ratio; CI, confidence interval.

![Supplementary Figure 7](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/02/2024.08.01.24311334/F9.medium.gif)

[Supplementary Figure 7](http://medrxiv.org/content/early/2024/08/02/2024.08.01.24311334/F9)

Supplementary Figure 7 Reverse Mendelian randomization analysis of early-onset ischemic stroke on risk of migraine
Due to a limited number of SNPs with a p-value < 5 × 10-8 for early-onset ischemic stroke, the threshold for instrumental variables was relaxed to < 1 × 10-5

Abbreviations: MR, Mendelian randomization; SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval.

## Acknowledgments

We would like to acknowledge the participants and investigators of the included genome wide association studies.

*   Received August 1, 2024.
*   Revision received August 1, 2024.
*   Accepted August 2, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.GBD 2016 Headache Collaborators. Global, regional, and national burden of migraine and tension-type headache, 1990-2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet Neurol 2018;17(11):954–976. doi: 10.1016/s1474-4422(18)30322-3
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1474-4422(18)30322-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30353868&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.08.01.24311334.atom) 

2.  2.Lee KJ, Lee SJ, Bae HJ, Sung J. Exploring the causal inference of migraine on stroke: A Mendelian randomization study. Eur J Neurol 2022;29(1):335–338. doi: 10.1111/ene.15101
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/ene.15101&link_type=DOI) 

3.  3.Shu MJ, Li JR, Zhu YC, Shen H. Migraine and Ischemic Stroke: A Mendelian Randomization Study. Neurol Ther 2022;11(1):237–246. doi: 10.1007/s40120-021-00310-y
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s40120-021-00310-y&link_type=DOI) 

4.  4.Hvitfeldt Fuglsang C, Pedersen L, Schmidt M, Vandenbroucke JP, Bøtker HE, Sørensen HT. Migraine and risk of premature myocardial infarction and stroke among men and women: A Danish population-based cohort study. PLoS Med 2023;20(6):e1004238. doi: 10.1371/journal.pmed.1004238
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1004238&link_type=DOI) 

5.  5.Bjornsdottir G, Chalmer MA, Stefansdottir L, Skuladottir AT, Einarsson G, Andresdottir M, et al. Rare variants with large effects provide functional insights into the pathology of migraine subtypes, with and without aura. Nat Genet 2023;55(11):1843–1853. doi: 10.1038/s41588-023-01538-0
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-023-01538-0&link_type=DOI) 

6.  6.Hautakangas H, Winsvold BS, Ruotsalainen SE, Bjornsdottir G, Harder AVE, Kogelman LJA, et al. Genome-wide analysis of 102,084 migraine cases identifies 123 risk loci and subtype-specific risk alleles. Nat Genet 2022;54(2):152–160. doi: 10.1038/s41588-021-00990-0
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00990-0&link_type=DOI) 

7.  7.Gormley P, Anttila V, Winsvold BS, Palta P, Esko T, Pers TH, et al. Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine. Nat Genet 2016;48(8):856–866. doi: 10.1038/ng.3598
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3598&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27322543&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.08.01.24311334.atom) 

8.  8.Wu X-P, Niu P-P, Liu H. Association between migraine and venous thromboembolism: a Mendelian randomization and genetic correlation study. Frontiers in Genetics 2024;15. doi: 10.3389/fgene.2024.1272599
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fgene.2024.1272599&link_type=DOI) 

9.  9.He C, Wu X, Lin L, Liu C, Li M, Jiang C, et al. Causal relationship between atrial fibrillation and stroke risk: a Mendelian randomization. J Stroke Cerebrovasc Dis 2023;32(12):107446. doi: 10.1016/j.jstrokecerebrovasdis.2023.107446
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jstrokecerebrovasdis.2023.107446&link_type=DOI) 

10. 10.Jaworek T, Xu H, Gaynor BJ, Cole JW, Rannikmae K, Stanne TM, et al. Contribution of Common Genetic Variants to Risk of Early-Onset Ischemic Stroke. Neurology 2022;99(16):e1738–e1754. doi: 10.1212/wnl.0000000000201006
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0000000000201006&link_type=DOI) 

11. 11.Mishra A, Malik R, Hachiya T, Jürgenson T, Namba S, Posner DC, et al. Stroke genetics informs drug discovery and risk prediction across ancestries. Nature 2022;611(7934):115–123. doi: 10.1038/s41586-022-05165-3
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-022-05165-3&link_type=DOI) 

12. 12.Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife 2018;7. doi: 10.7554/eLife.34408
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.34408&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29846171&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.08.01.24311334.atom) 

13. 13.Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Loh P-R, et al. An atlas of genetic correlations across human diseases and traits. Nature Genetics 2015;47(11):1236–1241. doi: 10.1038/ng.3406
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3406&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26414676&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.08.01.24311334.atom) 

14. 14.Hvitfeldt Fuglsang C, Pedersen L, Schmidt M, Vandenbroucke JP, Bøtker HE, Toft Sørensen H. Combined Impact of Migraine and Pregnancy-Induced Hypertension on Long-term Risk of Premature Myocardial Infarction and Stroke. Neurology 2024;102(1):e207813. doi: 10.1212/wnl.0000000000207813
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/wnl.0000000000207813&link_type=DOI) 

15. 15.van Welie FC, Kreft LA, Huisman JMA, Terwindt GM. Sex-specific metabolic profiling to explain the increased CVD risk in women with migraine: a narrative review. J Headache Pain 2023;24(1):64. doi: 10.1186/s10194-023-01601-5
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s10194-023-01601-5&link_type=DOI)