Modifiable factors for migraine prophylaxis: A mendelian randomization analysis

Objective: To examine the causal effect of potentially modifiable risk factors contributing to migraine pathogenesis.

Methods: We performed Mendelian randomization analyses and acquired data from United Kingdom Biobank, FinnGen Biobank, and the MRC IEU OpenGWAS data infrastructure. An inverse-variance weighted (IVW) model was used to examine the relationship between 51 potentially modifiable risk factors and migraine in 3215 participants with migraine without aura (MwoA), 3541 participants with migraine with aura (MwA), and 176,107 controls. We adopted a Bonferroni-corrected threshold of p = 9.8 × 10–4 (.05 divided by 51 exposures) as a sign of significant effect, and a p < .05 was considered as the sign of a suggestive association.

Results: More years of schooling significantly correlated with lower odds of MwoA pathogenesis (OR .57 [95%CI .44 to .75], p < .0001). More vitamin B12 intake (OR .49 [95%CI .24 to .99], p = .046) and lower level of stress [OR 8.17 (95%CI 1.5 to 44.36), p = .015] or anxiety disorder (OR 1.92 × 109 [95%CI 8.76 to 4.23*1017], p = .029) were suggestive to be correlated lower odds of MwoA pathogenesis. More coffee intake (OR .39 [95%CI .22 to .7], p = .001), lower level of eicosapentaenoic acid status (OR 2.54 [95%CI 1.03 to 6.26], p = .043), and more light physical activity (OR .09 [95%CI .01 to .94], p = .046) were suggestive to be associated with lower odds of MwA.

Conclusion: The years of schooling, light physical activity, vitamin B12 intake, and coffee intake were the protective factors for migraine; stress, anxiety, and eicosapentaenoic acid status were harmful factors. Interventions could be developed based on modifying these factors for migraine prophylaxis.

Highlights

• Potentially modifiable risk factors for migraine pathogenesis are unclear. The study addressed which are the modifiable risk factors contributing to migraine attacks.

• Based on Mendelian randomization analyses, we found that more years of schooling, vitamin B12 intake, and light physical activity were correlated with lower odds of migraine, while lower levels of stress, anxiety disorder, and eicosapentaenoic acid status were associated with lower odds of migraine.

• Comprehensive interventions could be developed based on modifying the years of schooling, vitamin B12 intake, light physical activity, levels of stress, anxiety disorder, and eicosapentaenoic acid status for migraine prevention.

Introduction

Migraine is a highly prevalent neurological disorder, characterized by frequent headache attacks accompanying photophobia, phonophobia, nausea, or vomiting. It is reported that migraine has a global prevalence of 10%, which brings a heavy burden on healthcare expenditure and becomes a major public health problem.

The management of migraine includes lifestyle modifications, pharmacological treatments, non-pharmacological treatments, and complementary and alternative therapies. Pharmacological treatments are recommended by guidelines and commonly practiced in clinical settings. However, the effectiveness of the pharmacological treatments varies in differential drug classes, and many of the treatments have low response rates for migraine prophylaxis (Shamliyan et al., 2013). The reason for the variation in the treatment effects is that most of the pharmacological treatments are not specifically developed for the preventive treatment of migraine. For example, calcium channel blockers and βblockers are recommended as first-line treatments for migraine prophylaxis, and they are developed for hypertension. High-quality evidence is still lacking for non-pharmacological treatments and complementary and alternative therapies, hampering their generation in practice.

Regarding the unmet treatment needs for migraine management, it is proposed that clinicians should use evidence-based multidisciplinary approaches to optimize clinical practices (Ashina et al., 2021). Clarifying the modifiable risk factors that are correlated with the pathogenesis of migraine is essential for developing multidisciplinary strategies for migraine prophylaxis and individualized approaches for patients with migraine. Many risk factors have been reported in recent 20 years, and most of them were revealed through observational studies. Owing to the disadvantages of the observational design, these findings are susceptible to confounding bias and reverse causation (Yarmolinsky et al., 2018; Hammerton and Munafò, 2021). Mendelian randomization, using genetic instruments as probes to detect causal relationships between exposures and disease outcomes, is perceived as a new study design that can minimize the risk of confounding bias and reverse causation. The nature of the randomized assignment of genetic effect alleles makes a Mendelian randomization study analog to a randomized trial (Davies et al., 2018). No study has been conducted to systematically search for potential risk factors that have causal effects on the pathogenesis of migraine, although several MR studies have fragmentarily explored a few of them.

We performed a comprehensive search for potentially modifiable risk factors that are reported in previous reviews or studies and implemented MR analyses to clarify whether the factors have causal effects on migraine pathogenesis.

Materials and methods

We acquired deidentified summary-level data from large-scale genome-wide association studies (GWASs). The data were openly available, and the GWAS studies had ethical approval from the participating centers. The design and conduct of this Mendelian randomization study conformed to the STROBE-MR (Strengthening the Reporting of Observational Studies in Epidemiology Using Mendelian Randomization) statement (Skrivankova et al., 2021).

Exposure data

Before the conduct of our MR study, we performed a literature search in Medline, looking for systematic reviews, meta-analyses, reviews with evidence summary, or studies in recent 3 years of modifiable factors for the prevention of migraine, using the keyword: migraine, risk factors, association factors, lifestyle, diet, physical activity, sleep, supplements, vitamin, and metabolism. We read the retrieved articles and screened for potentially modifiable risk factors that were associated with migraine. The reference articles were provided in Supplementary eTable S1.

Fifty-one risk factors were finally identified after the literature review. Nine factors were diet-related; eighteen were dietary supplements or serum indexes-related factors; thirteen were lifestyle-related; eleven were cardio-metabolic factors. The factors were listed in Table 1. Correlated SNPs of each factor were identified from the MRC IEU OpenGWAS data infrastructure (Hemani et al., 2018; Elsworth et al., 2020), which included a database of 214,846,534,918 genetic associations from 39,994 GWAS summary datasets. Each exposure factor and its related SNPs were further confirmed by reading the corresponding article IDs provided in the database.

TABLE 1

TABLE 1. Summary of modifiable risk factors for migraine in the Mendelian randomization analysis.

Most of the exposure factors were observed in the United Kingdom Biobank, which is a prospective cohort study with deep genetic, physical, and health data collected on ∼500,000 individuals across the United Kingdom from 2006 to 2010 (Bycroft et al., 2018). The other exposures were extracted from 12 GWAS studies, of which the PMIDs were listed in Table 1. We described the number of SNPs for each exposure trait and the pooled effect size of the SNPs on the trait.

Selection of instrumental variables

The instrumental variables (IVs), normally refer as SNPs in the Mendelian randomization studies, were selected according to predefined criteria. First, we selected SNPs with a correlation p-value <5×10⁸ in the reported GWAS studies. In the traits that had no SNP surviving the p-value limit, we relaxed the p-value to <1×10⁶ according to previous MR studies. Second, we clumped the data with an r² threshold of .001 and a distance window of 10,000 kB for minimizing bias from linkage disequilibrium (LD). Third, we harmonized the exposure data and the outcome data by removing SNPs for incompatible alleles or SNPs for being palindromic with intermediate allele frequencies, and we finally obtained a dataset containing data for both exposure traits and migraine outcomes which was used for Mendelian analysis.

Migraine outcomes

The outcome data, the summary-level genomic data of migraine, was acquired from the FinnGen study, which was launched in Finland in 2017. The FinnGen study, also known as Finn biobank, provides genome information and digital healthcare data of 506,800 participants (https://www.finngen.fi/en). We included 3215 participants with migraine without aura (MwoA), 3541 participants with migraine with aura (MwA), and 176,107 participants as controls. The participants with migraine were diagnosed with the International Classification of Diseases criteria (eighth version, code 346; ninth version, code 346; tenth version, code G43). These participants were of European ancestry, and both males and females were included. The data were released in the FinnGen biobank analysis round 5, and we imported the data from the MRC IEU OpenGWAS data infrastructure.

Statistical analysis

Using a two-sample MR analysis, we estimated the causal effects of each exposure (risk factor) on migraine, utilizing the inverse variance weighted (IVW) method. We extracted the SNP-exposure association (X) and SNP-migraine association estimates from the respective GWAS study, calculated the ratio of Y versus X for each SNP, and combined the ratios of all SNPs extracted for a specific exposure using the IVW method. For exposures reported with continuous data (eg, alcoholic intake frequency, years of schooling, and BMI), the combined estimates were explained as the log-odds of developing migraine with per genetically predicted standard deviation (SD) difference in the exposure. For exposures reported with binary data (eg, snoring, insomnia, and depression), the combined estimates were explained as the odds ratios (ORs) per genetically predicted unit difference in log-odds of having the relevant exposures. We drew scatter plots to show the direction of the causal effect and forest plots to show the effect size and whether the effect was statistically significant. To adjust for multiple comparison bias, we adopted a Bonferroni-corrected threshold of p = 9.8 × 10⁻⁴ (.05 divided by 51 exposures) as a sign of significant effect, while a p < .05 was considered as the sign of a suggestive association.

We performed heterogeneity tests using Cochran’s Q-statistic to downweigh or exclude SNPs with heterogeneous causal estimates. The Q-statistic is based on the first-order weights of the IVW model which assumes that there is no measurement error (NOME) in the genetic associations with the risk factor (Rees et al., 2019). We then performed a leave-one-out analysis to detect whether the causal effect of exposure was led by one dominant SNP.

We performed several sensitivity analyses to minimize confounding issues. First, we performed MR-Egger analysis, a meta-regression of the estimates of SNP-outcome association on the estimates of SNP-exposure association, to examine whether there was directional pleiotropy. We would assume no directional pleiotropy in the MR analysis when the intercept of the MR-Egger regression approached zero and the results were consistent between the IVW and the MR-Egger analysis. Second, we implemented the weighted median analysis, which had greater robustness than the IVW and the MR-Egger methods in handling individual genetic variants with strongly outlying causal estimates (Bowden and Holmes, 2019). Third, we carried out IVW radial analysis, an analysis that adopts a range of weighting specifications and determines outliers concerning their contribution to global heterogeneity. In the IVW radial model, a significance threshold of .05 and a tolerance threshold of .0001 were selected. Fourth, we implemented the Mendelian Randomization Pleiotropy RESidual Sum and Outlier (MR-PRESSO) test, which evaluates overall horizontal pleiotropy by comparing the observed distance of all the SNPs to the regression line (residual sum of squares) with the expected distance under the null hypothesis of no horizontal pleiotropy (Verbanck et al., 2018). We provided the results of the MR-PRESSO global test, and outlier-corrected estimates would be provided when there were indications of horizontal pleiotropy.

To assess weak instrument bias, we estimated the F-statistics for each exposure by using the method developed for two-sample MR (13). Larger F-statistic values indicate a stronger instrument variable, and we used a threshold of 10 to determine whether there was weak instrument bias (Burgess and Thompson, 2011). It is reported that F-statistic is not sufficient for evaluating weak instrument bias in the MR-Egger regression analysis, so we also reported I2 GX for each exposure. A high value of I2 GX, larger than 90%, indicates that weak instrument bias is unlikely to affect the MR-Egger analysis (Bowden et al., 2016). All the statistical analyses were performed in R 4.1.2 (www.r-project.org) with the TwoSampleMR .56 package.

Results

We included 51 modifiable risk factors (exposure traits) in this Mendelian randomization. Eleven traits were classified into diets, sixteen were correlated to diet supplements or diet metabolites, thirteen were correlated to lifestyle factors, and eleven were cardiometabolic parameters. The number of SNPs correlated to each trait ranged from 7 to 437. Forty-five traits had a proportion of variance explained by SNPs larger than 1%. All the included traits and their SNPs passed the test of weak instrument test, with the F-statistic larger than 10 and the I2 GX over 90%. Table 1 shows the summary of the included traits and corresponding SNPs.

Migraine without aura

Vitamin B12 intake was suggestive of a correlation with the pathogenesis of MwoA (OR .49 [95%CI .24 to .99], p = .046; Figure 1), and the weighted median analysis confirmed the finding (OR .35 [95%CI .15 to .85], p = .021; Supplementary eTable S2). As shown in Figure 2, the amount of vitamin B12 intake was inversely correlated to the incidence of migraine.

FIGURE 1

FIGURE 1. Odds ratios for associations between genetically predicted risk factors and migraine without aura. OR, odds ratio; SNP, single nucleotide polymorphisms. Annotation: We used the inverse-variance weighted method to calculate the odds ratios for individual SNPs. The odds ratio represents the odds ratio per genetically predicted SD unit increase in the risk factor. Take alcohol intake frequency for example, one SD unit increase in alcohol intake frequency was associated with 1.21 higher odds of migraine pathogenesis, which indicates that a higher frequency of alcohol intake was associated with a higher odds of migraine pathogenesis.

FIGURE 2

FIGURE 2. Scatter plots of the association between years of schooling, vitamin B12, and stress and migraine without aura. IVW, inverse variance weighted; MR, Mendelian randomization; SNP, single nucleotide polymorphisms. Annotation: Scatter plots of the migraine-SNP associations (y-axis) versus the exposures-SNP associations (x-axis) were shown, with vertical and horizontal lines showing 95% confidence intervals. The fitted curves among the dots refer to the differential statistical analysis methods. The MR analyses were performed with the inverse variance weighted method, the MR-Egger method, the weighted median method, and the IVW radial method. (A) shows the scatter plot for years of schooling, (B) shows the level of vitamin B12, and (C) shows the stress owing to financial difficulties.

The factor years of schooling was found significantly correlated to the pathogenesis of MwoA (OR .57 [95%CI .44 to .75], p < .0001; Figure 1), and the weighted median analysis confirmed the finding (OR .56 [95%CI .39 to .81], p = .002; Supplementary eTable S2). An increase in the years of schooling was correlated with a lower incidence of migraine (Figure 2).

The stress owing to financial difficulties was suggestive of a correlation with the pathogenesis of migraine (OR 8.17 [95%CI 1.5 to 44.36], p = .015; Figure 1), and the weighted median analysis confirmed the finding (OR 34.85 [95%CI 3.31 to 366.42], p = .003; Supplementary eTable S2). Figure 2 shows that higher odds of stress were correlated to a higher incidence of migraine.

The anxiety disorder was also suggestive of a correlation with migraine (OR 1.92*10⁹ [95%CI 8.76 to 4.23*10¹⁷], p = .029), but the finding was not confirmed in the sensitivity analyses.

Other factors, like body mass index (OR 1.17 [95%CI one to 1.38], p = .05), light physical activity (OR .11 [95%CI .01 to 1.25], p = .075; Figure 1), and diastolic blood pressure (OR 1.02 [95%CI one to 1.03], p = .068), were close to the suggestive boundary of a correlation with migraine pathogenesis. In addition, an MR-Egger analysis suggested a causal effect of total cholesterol level on migraine pathogenesis (OR .77 [95%CI .63 to .96], p = .068; Supplementary eTable S3), although it was not confirmed in the IVW analysis.

The other traits had no significant or causal effect on migraine without aura in all the analyses (Figures 1, 2, and Supplementary eTable S2–S4). Heterogeneity was not found in the analyses (Supplementary eTable S5).

Migraine with aura

Coffee intake was found suggestively correlated with the pathogenesis of MwA (OR .39 [95%CI .22 to .7], p = .001; Figure 3), but this finding was not supported by other sensitivity analyses (Supplementary eTable S6–S8). Coffee intake was inversely correlated with MwA, higher odds of genetically predicted coffee intake correlate with lower odds of MwA pathogenesis (Figure 4).

FIGURE 3

FIGURE 3. Odds ratios for associations between genetically predicted risk factors and migraine with aura. OR, odds ratio; SNP, single nucleotide polymorphisms. Annotation: We used the inverse-variance weighted method to calculate the odds ratios for individual SNPs. The odds ratio represents the odds ratio per genetically predicted SD unit increase in the risk factor.

FIGURE 4

FIGURE 4. Scatter plots of the association between coffee intake, EPA status, light physical activity and migraine with aura. EPA, Eicosapentaenoic acid; IVW, inverse variance weighted; MR, Mendelian randomization; SNP, single nucleotide polymorphisms. Annotation: Scatter plots of the migraine-SNP associations (y-axis) versus the exposures-SNP associations (x-axis) were shown, with vertical and horizontal lines showing 95% confidence intervals. The fitted curves among the dots refer to the differential statistical analysis methods. The MR analyses were performed with the inverse variance weighted method, the MR-Egger method, the weighted median method, and the IVW radial method. (A) shows the scatter plot of coffee intake, (B) shows EPA status, and (C) shows light physical activity.

Eicosapentaenoic acid status was suggestive of a correlation with MwA (OR 2.54 [95%CI 1.03 to 6.26], p = .043; Figure 3), but this finding was not confirmed in subsequent sensitivity analyses (Supplementary eTable S6–S8). A higher level of Eicosapentaenoic acid status was indicative of higher odds of MwA pathogenesis (Figure 4).

Light physical activity, instead of heavy physical activity, was suggested to correlate with MwA pathogenesis (OR .09 [95%CI .01 to .94], p = .046; Figure 3), although this finding was not confirmed in other sensitivity analyses. A higher frequency of light physical activity was correlated with lower odds of MwA pathogenesis (Figure 4).

Fasting glucose and HbA1c, as suggested by the weighted median analysis, were correlated with MwA pathogenesis (fasting glucose, OR 2.08 [95%CI 1.14 to 3.79], p = .016; HbA1c, OR 2.53 [95%CI 1.11 to 5.74], p = .02; Supplementary eTable S6). Higher levels of fasting glucose and HbA1c caused higher odds of MwA. The MR-Egger analysis showed that sleep duration correlated with MwA pathogenesis (OR 10.49 [95%CI 1.09 to 101.26], p = .046; Supplementary eTable S7).

Polyunsaturated fat, zinc status, vitamin D, years of schooling, depression, anxiety disorders, snoring, body mass index, and DBP were possibly correlated with MwA pathogenesis, which had p-values between .05 and .1 (Figure 3). The other risk factors had no significant or suggestive correlation with MwA pathogenesis (Figure 3, Supplementary eTable S6–S8). Heterogeneity was not found in the analyses (Supplementary eTable S9).

Discussion

We found that, in the Mendelian randomization analysis, genetically predicted years of schooling, intake of vitamin B12, and stress were causally correlated with MwoA. More years of schooling, a larger amount of intake of vitamin B12, and lower odds of stress were correlated with lower incidence of MwoA. We also found that genetically predicted coffee intake, EPA status, and light physical activity were causally correlated with MwA. A larger amount of coffee intake and higher odds of light physical activity were correlated with a lower risk of MwA, and a higher level of EPA was correlated with a higher risk of MwA. In the aforementioned correlations, only the years of schooling were statistically significant, the others were suggestive.

Our study adopted an MR design, which minimized confounding bias and reverse causation. Before performing the MR analysis, we implemented a comprehensive search for previously reported modifiable risk factors for migraine, so our study covered a wide range of risk factors for migraine, which is meaningful for developing multidisciplinary strategies for migraine prevention. In addition, we selected exposures mostly from United Kingdom Biobank and outcome data from the FinnGen study, which, to a certain extent, avoided the sample overlap and reduced bias from the overlap.

Years of schooling were found to be significantly associated with MwoA, and more years of schooling correlated with lower rates of MwoA. Previous studies reported that higher prevalence rates of migraine were correlated with lower education levels (Wang et al., 2016; Brusa et al., 2017; Brusa et al., 2019). We assumed that individuals with more years of education might be more voluntarily looking for access to migraine prophylaxis (Spadea et al., 2019) and acquiring a guideline-recommended strategy for migraine to avoid medication overuse headaches (Radtke and Neuhauser, 2012). This finding informed that headache education for a broader population is essential, and a recent randomized controlled trial showed that headache education could have a 2-day decrease in monthly headache days (Wells et al., 2021).

We found that a higher frequency of light physical activity was suggestively associated with lower odds of MwA and possibly correlative with MwoA (p = .075). These findings indicated a protective effect of light physical activity on migraine. Previous reviews concluded a beneficial effect of regular exercise on migraine prevention, but whether a higher exercise intensity leads to better prevention is unknown (Krøll et al., 2017; Amin et al., 2018; Lemmens et al., 2019; Barber and Pace, 2020). Our study demonstrated that light physical activity was more likely to be beneficial than heavy activity. The mechanism of exercise for migraine prophylaxis might include increased levels of beta-endorphin, endocannabinoid, and brain-derived neurotrophic factor in plasma after exercise (Amin et al., 2018).

Vitamin B12 and EPA were associated with the pathogenesis of migraine. More vitamin B12 intake caused lower odds of migraine, while higher EPA levels induced higher odds of migraine. Previous studies indicated that the high level of homocysteine was the cause of migraine, and B vitamins, especially vitamin B12, catalyzed homocysteine and therefore were beneficial for migraine prevention (Shaik et al., 2014; Urits et al., 2020). In addition, a recent systematic review showed that patients with migraine presented with a higher level of homocysteine and a lower level of folate, which also confirmed a relationship between B vitamins and migraine pathogenesis (Liampas et al., 2020). The mechanism of how EPA affects migraine is still unknown, which might warrant further investigation.

We found that more coffee intake was associated with a lower rate of MwoA, which might be the effect of caffeine. Caffeine, as a cure for acute migraine attacks, is assumed to antagonize adenosine and therefore decrease the production of NO, which is believed to be a major trigger of migraine attacks. Although migraine headaches might become chronic owing to caffeine withdrawal, it is reported that only 5%–30% of the patients with migraine had higher attack frequencies after caffeine withdrawal, and the migraine attack frequency might gradually decrease when patients keep a stable use of caffeine (<200 mg/day) (Nowaczewska et al., 2020).

Migraine was comorbid with several psychiatric disorders. In our study, we found that stress had a suggestive causal effect on migraine, and anxiety and depression were possibly associated with migraine pathogenesis. These findings were correlated with previous studies, and the effectiveness of amitriptyline for migraine confirmed the causal effect of the psychiatric disorders (Couch, 2011). The mechanism of stress, anxiety, and depression causing migraine mainly focused on the abnormal function, structure, and connectivity of brain regions—anterior cingulate cortex, anterior insula, prefrontal cortex, hippocampus, and amygdala (Minen et al., 2016). One recent study demonstrated that occipital bending was associated with migraine with visual aura, which in part explains psychiatric disorders in patients with MwA (Özkan and Gürsoy-Özdemir, 2021). Another study that found frontotemporal hypometabolism and increased frontal cortical thickness in patients with chronic migraine might explain some cognitive and behavioral pain-processing and sensory integration alterations in patients with migraine (Torres-Ferrus et al., 2021).

Taking prophylactic medications for migraine and the difference between episodic migraine and chronic migraine in disease nature might affect the results of the study. Owing to the lack of information on these points, we were unable to clarify whether these two factors are important confounding factors.

Our study still had several limitations. First, we adopted summary-level data in the analysis, which might ignore important covariates that influence the study result. Second, we screened potentially modifiable risk factors reported in previously published studies and examined these risk factors. Those that were unreported might also affect migraine pathogenesis. Third, genetically predicted exposures have a lifetime effect on migraine, so the effect size of these exposures might be exaggerated. Fourth, the GWAS study data were from the European population, and the results of this Mendelian randomization might only be generalized to the European population. Fifth, the risk factors that had an insignificant correlation with migraine in this study might be the consequence of low statistical power.

In conclusion, we found that years of schooling were significantly associated with migraine pathogenesis, and stress, light physical activity, vitamin B12 intake, coffee intake, and EPA status were suggestive of correlations. Comprehensive interventions might be developed based on modifying these factors for migraine prophylaxis.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: pan.baidu.com.

Ethics statement

The studies involving human participants were reviewed and approved by UK biobank and FinnGen biobank. The patients/participants provided their written informed consent to participate in this study.

Author contributions

All authors had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. HZ and MC designed the study. HZ, Y-ZS, and J-TL acquired the study data. HZ and L-LL analyzed and interpreted the data. MC provided advice and material support for conducting the study. HZ and Y-ZS wrote the first draft of the manuscript. All authors revised the manuscript and approved it for publication.

Funding

HZ received a grant from the Sichuan Youth Science and Technology Innovation Research Team (no. 2021JDTD0007). MC received a grant from the Hospital of Chengdu University of Traditional Chinese Medicine (Hundred Talents Program for Improving Scientific Research Capacity, no.20-B05).

Acknowledgments

We want to acknowledge the participants and investigators of the United Kingdom Biobank study and the FinnGen study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2023.1010996/full#supplementary-material

Supplementary eTable 1 | Reference articles reporting risk factors of migraine pathogenesis

Supplementary eTable 2 | Weighted median analysis for migraine without aura (MwoA).

Supplementary eTable 3 | MR-Egger analysis for migraine without aura (MwoA).

Supplementary eTable 4 | IVW radial analysis for migraine without aura (MwoA).

Supplementary eTable 5 | Heterogeneity test (MwoA).

Supplementary eTable 6 | Weighted median analysis of risk factors of MwA.

Supplementary eTable 7 | MR-Egger analysis of risk factors of MwA.

Supplementary eTable 8 | IVW Radial analysis (MwA).

Supplementary eTable 9 | Heterogeneity test (MwA).

References

Amin, F. M., Aristeidou, S., Baraldi, C., Czapinska-Ciepiela, E. K., Ariadni, D. D., and Di Lenola, D., (2018). The association between migraine and physical exercise. J. Headache Pain 19, 83. doi:10.1186/s10194-018-0902-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashina, M., Buse, D. C., Ashina, H., Pozo-Rosich, P., Peres, M. F. P., and Lee, M. J.,(2021). Migraine: Integrated approaches to clinical management and emerging treatments. Lancet 397, 1505–1518. doi:10.1016/S0140-6736(20)32342-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Barber, M., and Pace, A. (2020). Exercise and migraine prevention: A review of the literature. Curr. Pain Headache Rep. 24, 39. doi:10.1007/s11916-020-00868-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., Del Greco, M. F., Minelli, C., Davey Smith, G., Sheehan, N. A., and Thompson, J. R. (2016). Assessing the suitability of summary data for two-sample mendelian randomization analyses using MR-egger regression: The role of the I2 statistic. Int. J. Epidemiol. 45, 1961–1974. doi:10.1093/ije/dyw220

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., and Holmes, M. V. (2019). Meta-analysis and mendelian randomization: A review. Res. Synthesis Methods 10, 486–496. doi:10.1002/jrsm.1346

PubMed Abstract | CrossRef Full Text | Google Scholar

Brusa, P., Allais, G., Scarinzi, C., Baratta, F., Parente, M., and Rolando, S., (2019). Self-medication for migraine: A nationwide cross-sectional study in Italy. PLoS One 14, e0211191. doi:10.1371/journal.pone.0211191

PubMed Abstract | CrossRef Full Text | Google Scholar

Brusa, P., Parente, M., Allais, G., Rolando, S., Costa, G., and Gnavi, R., (2017). Community pharmacies as epidemiological sentinels of headache: First experience in Italy. Neurol. Sci. 38, 15–20. doi:10.1007/s10072-017-2908-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgess, S., and Thompson, S. G. (2011). Avoiding bias from weak instruments in Mendelian randomization studies. Int. J. Epidemiol. 40, 755–764. doi:10.1093/ije/dyr036

PubMed Abstract | CrossRef Full Text | Google Scholar

Bycroft, C., Freeman, C., Petkova, D., Band, G., Elliott, L. T., and Sharp, K., (2018). The UK Biobank resource with deep phenotyping and genomic data. Nature 562, 203–209. doi:10.1038/s41586-018-0579-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Couch, J. R. (2011). Amitriptyline in the prophylactic treatment of migraine and chronic daily headache. Headache 51, 33–51. doi:10.1111/j.1526-4610.2010.01800.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Davies, N. M., Holmes, M. V., and Smith, G. D. (2018). Reading mendelian randomisation studies: A guide, glossary, and checklist for clinicians. BMJ 362, k601. doi:10.1136/bmj.k601

PubMed Abstract | CrossRef Full Text | Google Scholar

Elsworth, B., Lyon, M., Alexander, T., Liu, Y., Matthews, P., and Hallett, J., (2020). The MRC IEU OpenGWAS data infrastructure. Genetics. doi:10.1101/2020.08.10.244293

CrossRef Full Text | Google Scholar

Hammerton, G., and Munafò, M. R. (2021). Causal inference with observational data: The need for triangulation of evidence. Psychol. Med. 51, 563–578. doi:10.1017/S0033291720005127

PubMed Abstract | CrossRef Full Text | Google Scholar

Hemani, G., Zheng, J., Elsworth, B., Wade, K. H., Haberland, V., and Baird, D., (2018). The MR-Base platform supports systematic causal inference across the human phenome. eLife 7, e34408. doi:10.7554/eLife.34408

PubMed Abstract | CrossRef Full Text | Google Scholar

Krøll, L. S., Hammarlund, C. S., Westergaard, M. L., Nielsen, T., Sloth, L. B., and Jensen, R. H., (2017). Level of physical activity, well-being, stress and self-rated health in persons with migraine and co-existing tension-type headache and neck pain. J. Headache Pain 18, 46. doi:10.1186/s10194-017-0753-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Lemmens, J., De Pauw, J., Van Soom, T., Michiels, S., Versijpt, J., and van Breda, E., (2019). The effect of aerobic exercise on the number of migraine days, duration and pain intensity in migraine: A systematic literature review and meta-analysis. J. Headache Pain 20, 16. doi:10.1186/s10194-019-0961-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Liampas, I., Siokas, V., Mentis, A-F. A., Aloizou, A-M., Dastamani, M., and Tsouris, Z., (2020). Serum homocysteine, pyridoxine, folate, and vitamin B12 levels in migraine: Systematic review and meta-analysis. Headache 60, 1508–1534. doi:10.1111/head.13892

PubMed Abstract | CrossRef Full Text | Google Scholar

Minen, M. T., Begasse De Dhaem, O., Kroon Van Diest, A., Powers, S., Schwedt, T. J., and Lipton, R., (2016). Migraine and its psychiatric comorbidities. J. Neurol. Neurosurg. Psychiatry 87, 741–749. doi:10.1136/jnnp-2015-312233

PubMed Abstract | CrossRef Full Text | Google Scholar

Nowaczewska, M., Wiciński, M., and Kaźmierczak, W. (2020). The ambiguous role of caffeine in migraine headache: From trigger to treatment. Nutrients 12, E2259. doi:10.3390/nu12082259

PubMed Abstract | CrossRef Full Text | Google Scholar

Özkan, E., and Gürsoy-Özdemir, Y. (2021). Occipital bending in migraine with visual aura. Headache 61, 1562–1567. doi:10.1111/head.14240

PubMed Abstract | CrossRef Full Text | Google Scholar

Radtke, A., and Neuhauser, H. (2012). Low rate of self-awareness and medical recognition of migraine in Germany. Cephalalgia 32, 1023–1030. doi:10.1177/0333102412454945

PubMed Abstract | CrossRef Full Text | Google Scholar

Rees, J. M. B., Wood, A. M., Dudbridge, F., and Burgess, S. (2019). Robust methods in Mendelian randomization via penalization of heterogeneous causal estimates. PLoS One 14, e0222362. doi:10.1371/journal.pone.0222362

PubMed Abstract | CrossRef Full Text | Google Scholar

Shaik, M. M., Tan, H. L., Kamal, M. A., and Gan, S. H. (2014). Do folate, vitamins B and B₁₂ play a role in the pathogenesis of migraine? The role of pharmacoepigenomics. CNS Neurol. Disord. Drug Targets 13, 828–835. doi:10.2174/18715273113129990112

PubMed Abstract | CrossRef Full Text | Google Scholar

Shamliyan, T. A., Choi, J-Y., Ramakrishnan, R., Miller, J. B., Wang, S-Y., and Taylor, F. R., (2013). Preventive pharmacologic treatments for episodic migraine in adults. J. Gen. Intern Med. 28, 1225–1237. doi:10.1007/s11606-013-2433-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Skrivankova, V. W., Richmond, R. C., Woolf, B. A. R., Yarmolinsky, J., Davies, N. M., and Swanson, S. A., (2021). Strengthening the reporting of observational studies in Epidemiology using mendelian randomization: The STROBE-MR statement. JAMA 326, 1614–1621. doi:10.1001/jama.2021.18236

PubMed Abstract | CrossRef Full Text | Google Scholar

Spadea, T., Scarinzi, C., Baratta, F., Allais, G., Rolando, S., and Manzoni, G. C., (2019). Access to migraine centres by educational level of patients and awareness of the disease. Neurol. Sci. 40, 207–209. doi:10.1007/s10072-019-03819-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Torres-Ferrus, M., Pareto, D., Gallardo, V. J., Cuberas-Borrós, G., Alpuente, A., and Caronna, E., (2021). Cortical metabolic and structural differences in patients with chronic migraine. An exploratory 18FDG-PET and MRI study. J. Headache Pain 22, 75. doi:10.1186/s10194-021-01289-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Urits, I., Yilmaz, M., Bahrun, E., Merley, C., Scoon, L., and Lassiter, G., (2020). Utilization of B12 for the treatment of chronic migraine. Best. Pract. Res. Clin. Anaesthesiol. 34, 479–491. doi:10.1016/j.bpa.2020.07.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Verbanck, M., Chen, C-Y., Neale, B., and Do, R. (2018). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 50, 693–698. doi:10.1038/s41588-018-0099-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Xing, Y., Sun, J., Zhou, H., Yu, H., and Zhao, Y., (2016). Prevalence, associated factors, and impact on quality of life of migraine in a community in northeast China. J. Oral Facial Pain Headache 30, 139–149. doi:10.11607/ofph.1584

PubMed Abstract | CrossRef Full Text | Google Scholar

Wells, R. E., O’Connell, N., Pierce, C. R., Estave, P., Penzien, D. B., and Loder, E., (2021). Effectiveness of mindfulness meditation vs headache education for adults with migraine: A randomized clinical trial. JAMA Intern Med. 181, 317–328. doi:10.1001/jamainternmed.2020.7090

PubMed Abstract | CrossRef Full Text | Google Scholar

Yarmolinsky, J., Wade, K. H., Richmond, R. C., Langdon, R. J., Bull, C. J., and Tilling, K. M., (2018). Causal inference in cancer Epidemiology: What is the role of mendelian randomization? Cancer Epidemiol. Biomarkers Prev. 27, 995–1010. doi:10.1158/1055-9965.EPI-17-1177

PubMed Abstract | CrossRef Full Text | Google Scholar