Genetic variants in migraine: a field synopsis and systematic re-analysis of meta-analyses
The Journal of Headache and Pain volume 21, Article number: 13 (2020)
Numerous genetic variants from meta-analyses of observational studies and GWAS were reported to be associated with migraine susceptibility. However, due to the random errors in meta-analyses, the noteworthiness of the results showing statistically significant remains doubtful. Thus, we performed this field synopsis and re-analysis study to evaluate the noteworthiness using a Bayesian approach in hope of finding true associations.
Relevant meta-analyses from observational studies and GWAS examining correlation between all genetic variants and migraine risk were included in our study by a PubMed search. Identification of noteworthy associations were analyzed by false-positive rate probability (FPRP) and Bayesian false discovery probability (BFDP). Using noteworthy variants, GO enrichment analysis were conducted through DAVID online tool. Then, the PPI network and hub genes were performed using STRING database and CytoHubba software.
As for 8 significant genetic variants from observational studies, none of which showed noteworthy at prior probability of 0.001. Out of 47 significant genetic variants in GWAS, 36 were noteworthy at prior probability of 0.000001 via FPRP or BFDP. We further found the pathways “positive regulation of cytosolic calcium ion concentration” and “inositol phosphate-mediated signaling” and hub genes including MEF2D, TSPAN2, PHACTR1, TRPM8 and PRDM16 related to migraine susceptibility.
Herein, we have identified several noteworthy variants for migraine susceptibility in this field synopsis. We hope these data would help identify novel genetic biomarkers and potential therapeutic target for migraine.
Migraine is a complex and incapacitating neurologic condition with more than one billion individuals suffering from and imposes a huge socioeconomic burden worldwide [1, 2]. It is characterized by recurrent episodes of unilateral throbbing pain often accompanied with nausea, phonophobia and photophobia, leading to a decline in life quality or even disability . In the Global Burden of Disease (GBD) 2015, migraine ranked the seventh among the leading causes of years lived with disability for all ages and the third for ages 15 to 49 years . Recently, it has been reported that, in the United States, migraine affected almost 15% adults and the annual economic cost was over $ 2600 each person diagnosed with episodic migraine and $ 8000 for those with chronic migraine [5, 6]. Thus, a large amount of studies was performed to explore the risk factors and pathogenesis of migraine over the past decades.
Except for the common risk factors such as obesity, medication overuse, poor sleep, caffeine and stressful life events, the genetic factors for migraine susceptibility was drawing more and more attention [7, 8]. Numerous genetic polymorphisms from meta-analyses of observational studies and GWAS were reported to be associated with migraine susceptibility . However, due to the random errors in meta-analyses leading to false-positive results, the noteworthiness of the results with statistical significance remains doubtful. Therefore, in this field synopsis, we summarized and re-analyzed all significant genetic variants from meta-analyses of observational studies and GWAS, then assessed their noteworthiness using Bayesian procedures including false-positive rate probability (FPRP) and Bayesian false discovery probability (BFDP) and discussed possible molecular mechanisms for migraine occurrence.
Search strategy and data extraction
A comprehensive literature search was conducted in the PubMed database up to 31 July 2019, by using the following terms: “migraine/ headache” and “meta-analysis” and “polymorphism/genome-wide association study/SNP/GWAS/variant/allele/genotype”. Studies were all selected according to the following criteria: (1) meta-analysis design study; (2) evaluating the association between genetic polymorphisms and migraine risk; (3) raw data available including odds ratios (ORs), 95% confidence intervals (CIs) or other information necessary for FPRP and BFDP calculation; (4) studies published in English. Data including author, published year, genetic variant, genetic model, OR, 95% CI, P-value, ethnicity, type of migraine, number of cases and controls, heterogeneity and publication bias were extracted from the meta-analyses included according to above criteria.
Assessment methods for meta-analysis
In order to figure out the noteworthiness of meta-analysis on the association between genetic polymorphisms and migraine susceptibility, two novel statistic methods, FPRP and BFDP, were applied in our study. FPRP is the probability that no true association exists between genetic variant and disease drawing a statistically significant finding . The magnitude of the FPRP is determined by prior probability, statistical power and observed P-value. Owing to the highly subjective prior probability, we analyzed a wide range of prior probability. We calculated FPRP values using two levels of prior probabilities: at 0.05/0.001 (medium/low prior level) that would be expected for candidate SNPs, and at 0.001/0.000001 (medium/low prior level) for GWAS SNPs. Moreover, the lower the prior probability, the more reliable the result. Besides, we used statistical power to detect ORs of 1.2 and 1.5 for computing FPRP at each prior probability. For the statistically significant SNPs (95% CI that excluded 1 and P-value which was lower than 5 × 10− 8 for meta-analysis of GWAS or 0.05 for observational studies), we calculated FPRP by using the Excel spreadsheet offered by Wacholder (http://jncicancerspectum.oupjournals.org/jnci/content/vol96/issue6) . FPRP values lower than 0.2 were considered to be noteworthy.
BFDP is the fact that if an association is reported as noteworthy, BFDP is the probability of a false discovery . Relevant data and prior probability applied in the calculation of BFDP were same as it of FPRP. The BFDP was estimated by using the excel Calculation Spreadsheet (http://faculty.washington.edu/jonno/cv.html) . And BFDP with the values of less than 0.8 were considered to be a noteworthy significant association. Different from FPRP, BFDP is a new statistical method based on logistic regression model rather than standard normal distribution and doesn’t rely on the statistical power. Thus, BFDP has a sounder methodological basis. Nevertheless, Wakefield admits that there is no significant difference in the overall behavior between these two approaches . Therefore, we presented the both results of FPRP and BFDP which allowed readers to reach a more comprehensive judgment.
Besides, for the meta-analysis of observational studies, summary evidence was also evaluated using Venice criteria, which have been described in detail previously [12, 13]. Briefly, we classified the strength of credibility into A, B and C grades that were separately characterized as strong, moderated and weak in three parameters including amount of evidence, replication of association and the protection from bias. According to this criteria, high credibility was defined as including A grades only, intermediate credibility was composed of A and B grades and low credibility was one or more C grades.
Joint population attributable risk calculation
We evaluated the cumulative effect of all noteworthy SNPs on migraine susceptibility. We used the minor allelic frequency (MAF) to calculated the population attributable risk (PAR) and further estimated the Joint PAR% for the SNPs showing noteworthiness during the computing of BFDP (BFDP < 0.8) or FPRP (FPRP< 0.2) at a prior probability of 10− 6 assumed for GWAS SNPs and 10− 3 for candidate SNPs at a statistical power to detect the OR of 1.5.
GO and enrichment pathway analysis and PPI network construction
Gene ontology (GO) analysis is a useful bioinformatic method for annotating genes. In our study, GO enrichment analyses were carried out using a list of genes with noteworthy SNPs through DAVID online tool (http://david.abcc.ncifcrf.gov/). Then, we applied the STRING 11.0 network database to construct a protein-protein interaction (PPI) network. We set the minimum required interaction score at 0.15, no more than five interactors and four active interaction sources (experiments, text mining, co-expression and database) for PPI construction. Finally, we used the Cytoscape (version 3.4.2) and cytoHubba to detect hub genes meanwhile the cut-off criterion of hub genes was setting as degree ≥7.
First, a total of 89 articles were identified according to our search strategy. Second, 38 obviously irrelevant articles were excluded by screening the titles and abstracts. After reviewing the full-text articles, 16 were further excluded, among which 9 articles were not meta-analyses and 7 were not correlated with migraine susceptibility. At last, 35 articles were selected into our study [14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48]. The screening process of the articles was performed by two independent researchers and shown in Fig. 1. Overall data retrieved from meta-analyses of observational studies and GWAS on risk of migraine were summarized in Additional file 1: Table S1 and Additional file 2: Table S2. Third, when more than one meta-analysis was available for certain SNP, we only included the results from more recent meta-analysis with larger sample size. Finally, the statistically significant results with P < 0.05 for observational studies and P < 5 × 10− 8 for GWAS were summarized in Tables 1 and 2.
As shown in Table 1, 5 genes with 8 genetic variants from observational studies were found to be significant after excluding the overlapping data. At a prior probability of 0.05, we identified the genetic variants, MTHFR/rs1801133, noteworthy via FPRP estimation with a statistical power to detect OR of 1.2. Likewise, 4 genetic variants including MTHFR/rs1801133, BDNF/rs6265, ESR1/rs1801132 and 5-HTT/VNTR showed noteworthy FPRP values to detect OR of 1.5. As for the statistical method BFDP, the same 4 genetic variants were noteworthy. However, when it comes to the re-analysis at a prior probability of 0.001, no noteworthy relationship between genetic variants and migraine risk could be detected via FPRP and BFDP.
Moreover, we performed subgroup analysis of observational studies based on the migraine subtypes (Table 3). In the migraine with aura subgroup, only 1 (TNF-α/rs1800629) and 3 (MTHFR/rs1801133, ESR1/rs1801132 and TNF-α/rs1800629) genetic variants was noteworthy in FPRP at the prior probability of 0.05 with a statistical power to detect an OR of 1.2 and 1.5, respectively. And all the variants with noteworthy FPRP values were also noteworthy in BFDP estimation. Similarly, none was identified to be noteworthy in FPRP and BFDP when the prior probability was 0.001. In migraine without aura, we did not observe any noteworthy SNPs. Subsequent subgroup analyses based on ethnicity were performed and displayed in Table 4. Compared with 4 noteworthy SNPs identified at the prior probability of 0.05 containing MTHFR/rs1801133, BDNF/rs6265, ESR1/rs1801132 and ESR1/rs2228480 for Caucasian population, none of candidate polymorphisms was considered to be noteworthy in non-Caucasian population.
In addition to FPRP and BFDP, we also used the Venice criteria to evaluate the credibility of the meta-analyses characterized by low, intermediate and high level (Table 1). We observed a consistency between the noteworthiness measured by FPRP and BFDP and Venice criteria score for candidate SNPs, that is, most of the noteworthy SNPs observed in our re-analysis were with high or intermediate level of evidence, with the exception of MTHFR/rs1801133.
As shown in Table 2, within the data extracted from the meta-analyses of GWAS, 40 genes with 47 genetic variants were statistically significant (P < 5 × 10− 8). According to the results of the re-analyses, 32 and 26 SNPs were found to be noteworthy via FPRP estimation at the statistical power to detect the OR of 1.2 with the prior probability of 0.001 and 0.000001. 34 and 26 SNPs were identified as noteworthy to detect the OR of 1.5 via FPRP. As for BFDP, 40 and 35 SNPs were noteworthy under the prior probability of 0.001 and 0.000001. Almost all noteworthy variants calculated through FPRP also showed noteworthy BFDP values. Only three variants (YAP1/rs10895275, TGFBR2/rs6791480 and MPPED2/ rs11031122) were noteworthy in the calculation of FPRP, but not via BFDP. In addition, we detected 9 noteworthy SNPs (LRP1/rs11172113, PRDM16/rs10218452, FHL5/rs67338227, TSPAN2/rs2078371, PLCE1/ rs10786156, KCNK5/rs10456100, CFDP1/rs77505915, NRP1/rs2506142 and ADAMTSL4/rs6693567) in the computation of BFDP rather than FPRP, which could be explained by the fact that, in some cases, noteworthiness could not be assessed on account of a mathematical error during the calculation of the inverse of cumulative normal distribution in FPRP. Furthermore, in the subgroup analysis, we only performed the re-analysis of migraine without aura but not migraine with aura owing to the lack of raw data. We found that all ten statistically significant variants within 8 genes were noteworthy in migraine without aura subgroup (Table 5). In addition, due to the population included in the GWAS were almost Caucasians, subgroup re-analysis based on ethnicity cannot be performed.
In order to provide readers with the predictive potential of the migraine risk, we calculated the joint PAR which was recently argued to be more credible than P-values or ORs by combining 36 noteworthy variants via FPRP or BFDP. And the value of joint PAR in our study was 44.2%, indicating that these involved SNPs together may attribute to an obvious increment in the risk of migraine and this method may be a useful way for screening migraine susceptibility and identify at-risk populations (Table 6).
To explore the potential biology mechanism for migraine, GO enrichment analysis was performed. As shown in Fig. 2, twelve GO terms were identified including “cell-cell signaling” (GO: 0007267), “inositol phosphate-mediated signaling” (GO:0048016), “positive regulation of cytosolic calcium ion concentration” (GO:0007204), “integral component of plasma membrane” (GO:0005887), “adult heart development” (GO:0007512), “regulation of smooth muscle contraction” (GO:0006940), “activating transcription factor binding” (GO:0033613), “sprouting angiogenesis” (GO:0002040), “patterning of blood vessels” (GO:000156), “angiogenesis” (GO:0001525), “receptor activity” (GO:0004872) and “protein kinase C-activating G-protein coupled receptor signaling pathway” (GO:0007205). To detect the interactive relationships and potential hub nodes, PPI network of all genes with noteworthy SNPs were constructed and displayed in Additional file 3: Figure S1. Moreover, when degree ≥7 as the cut-off criterion, 5 genes (MEF2D, TSPAN2, PHACTR1, TRPM8 and PRDM16) were selected as hub genes by using CytoHubba software.
Despite the serious impact of migraine on human health and its burden on public, understanding of the pathogenesis, valid prevention and treatment remains limited [49, 50]. Recent years, the genetic component of migraine has received widespread attention. Moreover, a growing number of meta-analyses from observational studies and GWAS were performed and further identified numerous significant genetic variants, which displayed important insights into the mechanisms underlying migraine development. However, due to false-positive possibility in meta-analyses, we performed the first re-analysis of meta-analyses of genetic association studies in migraine in hope of finding noteworthy associations. We used Bayesian methods and Venice criteria to evaluate the credibility of genetic associations. In this work, we synthesized all relevant data from the meta-analyses which detected the association between genetic polymorphisms and migraine risk before 31 July 2019.
As for the candidate genes, we found 4 genetic variants noteworthy at prior probability of 0.05 which were consistent with the results of Venice criteria. However, when we raised the prior probability standard from 0.05 to 0.001, we did not find any gene noteworthy under FPRP or BFDP, which indicated that the results of observational studies should be interpreted with caution. In comparison to the observational studies, according to GWAS, out of the 47 significant genetic variants, 36 were considered to be noteworthy at prior probability of 0.000001 via FPRP or BFDP. Most significant variants identified in GWAS remained noteworthy at the prior probability of 0.000001 demonstrated that the results from GWAS were more reliable.
We also detected the pathways involved in migraine by conducting functional enrichment analysis and further explore the possible molecular mechanisms. We found several significant pathways (cell-cell signaling, inositol phosphate-mediated signaling, positive regulation of cytosolic calcium ion concentration, integral component of plasma membrane and adult heart development, etc.) and 5 hub genes (MEF2D, TSPAN2, PHACTR1, TRPM8 and PRDM16) which were considered to play vital roles in migraine occurrence.
To date, the pathophysiology of migraine is partially understood that the headache of migraine is associated with activation and sensitization of trigeminovascular system [51, 52]. This hypothesis is based on a fact that migraine attack was originated from the activation of nociceptors that are initialed from trigeminal ganglion and innervate arachnoid, pial, blood vessels of dural, large cerebral arteries and sinuses . These activated nociceptors released some inflammatory mediators and vasoactive neuropeptides such as calcitonin gene-related peptide (CGRP), neurokinin A and substance P causing vasodilation of dural and pial vessels and neurogenic inflammation which finally results in central sensitization causing headache of migraine to some extent [53,54,55].
PLCE1 and NMUR2, whose variants showed noteworthy association with migraine susceptibility, were found to be involved in two pathways called “positive regulation of cytosolic calcium ion concentration” and “inositol phosphate-mediated signaling” detected in our study. PLCE1 is a gene encodes an enzyme called phospholipase C that promotes the generation of inositol triphosphate (IP3) and further direct the calcium mobilization by initiate the release of calcium ion stored in the endoplasmic reticulum through IP3 receptor leading to the increment of cytosolic calcium concentration in neurons [56, 57]. As for NMUR2, an impaired pain response was observed in NMUR2-deficient mice, supporting the hypothesis that nociceptive effects may be partially mediated through NMUR2 . Similar with PLCE1, NMUR2 is abundantly expressed in the central nervous system and considered as a regulator for intracellular calcium mobilization via IP3 as well . Furthermore, the increasing of cytosolic calcium may cause the release of CGRP that finally results in migraine development [60, 61].
Within the hub genes detected in our study, TRPM8 gene was found to play an important role in the pathophysiology of migraine. As a cation channel, TRPM8 is characterized as a cold temperature detector firstly [62, 63]. However, recently, it is also discovered as a parameter for ongoing persistent pain with several evidence. TRPM8 null mice showed a significant decrease in the injury-induced response . On the opposite, the activation of TRPM8 can lead to the depolarization for nerve endings and afferent impulse into central nervous system. Moreover, these activated neurons by noxious cold temperature have properties similar with nociceptors . Mechanistically, a functional linkage was observed between TRPM8 and CGRP, that is, CGRP release was deficient in neurons without TRPM8 triggering and the release of CGRP was closely related to neurogenic inflammation and future migraine occurrence .
As for PHACTR1, it might be involved in migraine attack by regulating vasomotor tone. Concretely, PHACTR1 could bind with protein phosphatase 1 (PP1) and its gene silencing has been demonstrated to decrease the activity of PP1 . In arteries, PP1 may play an important role in vasomotor tone by mediating calcium cycling and contractility in smooth muscle cells . Also, PHACTR1 could regulate dendritic morphology and synaptic activity by interacting with PP1 in nervous system, and was further thought to be implicated with the pathophysiology of migraine [69, 70].
Different with the hub genes mentioned above, PRDM16, a zinc-finger nuclear protein, works as an important activator of brown adipogenesis, and depletion of PRDM16 may lead to a significant loss in brown adipocyte identity . Recently, brown adipose tissue has been in the focus of metabolism research, which can dissipate energy by the regulation of uncoupling protein-1, increase fatty acid oxidation and heat production and counteract obesity [72, 73]. Moreover, PRDM16 rs2651899 is an intron variant that may affect the splicing of PRDM16 and its downstream regulatory elements, reduce PRDM16 expression and thus increase body mass index . Large cohort studies suggest that obesity is a risk factor for migraine even after adjusting for comorbidities . In obese individuals, the expression of many inflammatory mediators was increased including CGRP and interleukins, which in turn could cause central sensitization in migraine pathophysiology [76, 77]. In addition, Obesity is a state of sympathetic activation, that may contribute to increase in migraine attack .
For other hub genes for migraine, TSPAN2 is highly expressed in oligodendrocyte cell lines and may regulate the differentiation process of oligodendrocytes to myelin-forming glia . As a transcription factor in neurons, MEF2D is concerned to be involved in neurogenesis, neuronal survival and differentiation by controlling MEF2D-dependent gene transcription . However, up till now, the potential role of TSPAN2 and MEF2D in migraine development still remains unknown. Further functional studies are required to explore the underlying biological mechanisms.
Although our study included a large number of relevant articles on migraine, some limitations should be acknowledged in our re-analysis. First, most of original studies did not consider the potential confounders such as gene-gene and gene-environment interactions. Second, there exist some inherent methodological flaws in observational studies, such as selection bias, publication bias and small sample size. Although Venice criteria were applied, these potential biases were difficult to measure. Additionally, owing to the inadequate raw data, Venice Score cannot be applied for GWAS. Third, the GWAS involved in our report were mostly performed in Caucasian population, thus, the results of which was restricted to be applied in diverse populations. Also, GWAS included in our study were lack of subgroup analysis on migraine subtype. Thus, we could only study the noteworthiness of GWAS SNPs in migraine without aura subgroup. Lastly, we constructed a PPI network to explore the underlying biological mechanism for migraine. However, the criteria for PPI construction was relatively subjective, so the result of PPI might need to be interpreted with caution.
The current findings identify several noteworthy variants for migraine susceptibility. We hope this field synopsis and systematic re-analysis would help identify novel genetic biomarkers and potential therapeutic target for migraine.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article.
Diener HC et al (2015) Integrated care for chronic migraine patients: epidemiology, burden, diagnosis and treatment options. Clin Med (Lond) 15(4):344–350
Collaborators GBDH (2018) 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 17(11):954–976
Younger DS (2016) Epidemiology of migraine. Neurol Clin 34(4):849–861
Disease, G.B.D., I. Injury, and C. Prevalence (2016) Global, regional, and national incidence, prevalence, and years lived with disability for 310 diseases and injuries, 1990–2015: a systematic analysis for the global burden of Disease study 2015. Lancet 388(10053):1545–1602
MacGregor EA (2017) Migraine. Ann Intern Med 166(7):ITC49–ITC64
Messali A et al (2016) Direct and indirect costs of chronic and episodic migraine in the United States: a web-based survey. Headache 56(2):306–322
Lipton RB, Bigal ME (2005) Migraine: epidemiology, impact, and risk factors for progression. Headache 45(Suppl 1):S3–S13
Bohm PE, Stancampiano FF, Rozen TD (2018) Migraine headache: updates and future developments. Mayo Clin Proc 93(11):1648–1653
Kondratieva N et al (2016) Biomarkers of migraine: part 1 - genetic markers. J Neurol Sci 369:63–76
Wacholder S et al (2004) Assessing the probability that a positive report is false: an approach for molecular epidemiology studies. J Natl Cancer Inst 96(6):434–442
Wakefield J (2007) A Bayesian measure of the probability of false discovery in genetic epidemiology studies. Am J Hum Genet 81(2):208–227
Khoury MJ et al (2009) Genome-wide association studies, field synopses, and the development of the knowledge base on genetic variation and human diseases. Am J Epidemiol 170(3):269–279
Langevin SM et al (2010) Assessment of cumulative evidence for the association between glutathione S-transferase polymorphisms and lung cancer: application of the Venice interim guidelines. Pharmacogenet Genomics 20(10):586–597
Oterino A et al (2007) Genetic association study and meta-analysis of the HTR2C Cys23Ser polymorphism and migraine. J Headache Pain 8(4):231–235
Rubino E et al (2009) Association of the C677T polymorphism in the MTHFR gene with migraine: a meta-analysis. Cephalalgia 29(8):818–825
Anttila V et al (2010) Genome-wide association study of migraine implicates a common susceptibility variant on 8q22.1. Nat Genet 42(10):869–873
Schurks M, Rist PM, Kurth T (2010) Sex hormone receptor gene polymorphisms and migraine: a systematic review and meta-analysis. Cephalalgia 30(11):1306–1328
Schurks M, Rist PM, Kurth T (2010) 5-HTTLPR polymorphism in the serotonin transporter gene and migraine: a systematic review and meta-analysis. Cephalalgia 30(11):1296–1305
Schurks M, Rist PM, Kurth T (2010) STin2 VNTR polymorphism in the serotonin transporter gene and migraine: pooled and meta-analyses. J Headache Pain 11(4):317–326
Schurks M, Rist PM, Kurth T (2010) MTHFR 677C>T and ACE D/I polymorphisms in migraine: a systematic review and meta-analysis. Headache 50(4):588–599
Chasman DI et al (2011) Genome-wide association study reveals three susceptibility loci for common migraine in the general population. Nat Genet 43(7):695–698
Ligthart L et al (2011) Meta-analysis of genome-wide association for migraine in six population-based European cohorts. Eur J Hum Genet 19(8):901–907
Liu H et al (2011) Association of 5-HTT gene polymorphisms with migraine: a systematic review and meta-analysis. J Neurol Sci 305(1–2):57–66
Samaan Z et al (2011) Methylenetetrahydrofolate reductase gene variant (MTHFR C677T) and migraine: a case control study and meta-analysis. BMC Neurol 11:66
Schurks M et al (2011) Tumour necrosis factor gene polymorphisms and migraine: a systematic review and meta-analysis. Cephalalgia 31(13):1381–1404
Gu L et al (2012) The TNF-alpha-308G/a polymorphism is associated with migraine risk: a meta-analysis. Exp Ther Med 3(6):1082–1086
Miao J, Wang F, Fang Y (2012) Association of 231G>a polymorphism of endothelin type a receptor gene with migraine: a meta-analysis. J Neurol Sci 323(1–2):232–235
Tammimaki A, Mannisto PT (2012) Catechol-O-methyltransferase gene polymorphism and chronic human pain: a systematic review and meta-analysis. Pharmacogenet Genomics 22(9):673–691
Anttila V et al (2013) Genome-wide meta-analysis identifies new susceptibility loci for migraine. Nat Genet 45(8):912–917
Esserlind AL et al (2013) Replication and meta-analysis of common variants identifies a genome-wide significant locus in migraine. Eur J Neurol 20(5):765–772
Liu R et al (2014) MTHFR C677T polymorphism and migraine risk: a meta-analysis. J Neurol Sci 336(1–2):68–73
Liu R et al (2014) Effects of tumor necrosis factor-beta (TNF-beta) 252A>G polymorphism on the development of migraine: a meta-analysis. PLoS One 9(6):e100189
Peng JM et al (2014) Meta-analysis of 5-hydroxytryptamine type 2A receptor polymorphisms and migraine susceptibility. Int J Neurosci 124(12):882–889
Chen H et al (2015) Association between polymorphisms of DRD2, COMT, DBH, and MAO-A genes and migraine susceptibility: a meta-analysis. Medicine (Baltimore) 94(47):e2012
Chen M et al (2015) Tumor necrosis factor (TNF) -308G>a, nitric oxide synthase 3 (NOS3) +894G>T polymorphisms and migraine risk: a meta-analysis. PLoS One 10(6):e0129372
Jacobsen KK et al (2015) Genome wide association study identifies variants in NBEA associated with migraine in bipolar disorder. J Affect Disord 172:453–461
Li L et al (2015) Impact of ESR1 gene polymorphisms on migraine susceptibility: a meta-analysis. Medicine (Baltimore) 94(35):e0976
Lippi G, Mattiuzzi C, Cervellin G (2015) Meta-analysis of factor V Leiden and prothrombin G20210A polymorphism in migraine. Blood Coagul Fibrinolysis 26(1):7–12
Miao J et al (2015) Association of the Apolipoprotein E polymorphism with migraine: a meta-analysis. BMC Neurol 15:138
Gormley P et al (2016) Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine. Nat Genet 48(8):856–866
Wan D et al (2016) Association between angiotensin-converting enzyme insertion/deletion polymorphism and migraine: a meta-analysis. Int J Neurosci 126(5):393–399
Cai X et al (2017) The association between brain-derived neurotrophic factor gene polymorphism and migraine: a meta-analysis. J Headache Pain 18(1):13
Liao YJ, Jiang JR, Jin SQ (2017) The association between COMT Val158Met polymorphism and migraine risk: a meta-analysis. Cephalalgia 37(6):592–598
Terrazzino S et al (2017) Brain-derived Neurotrophic factor Val66Met gene polymorphism impacts on migraine susceptibility: a meta-analysis of case-control studies. Front Neurol 8:159
Chang X et al (2018) Common variants at 5q33.1 predispose to migraine in African-American children. J Med Genet 55(12):831–836
Dong H et al (2018) Endothelial nitric oxide synthase (−786T>C) polymorphism and migraine susceptibility: a meta-analysis. Medicine (Baltimore) 97(36):e12241
Gao X, Wang J (2018) Quantitative assessment of the association between GRIA1 polymorphisms and migraine risk. Biosci Rep 38(6). https://doi.org/10.1042/BSR20181347
Liu L et al (2019) Effects of MTHFR C677T and A1298C polymorphisms on migraine susceptibility: a meta-analysis of 26 studies. Headache 59(6):891–905
Dodick DW, Gargus JJ (2008) Why migraines strike. Sci Am 299(2):56–63
Wessman M et al (2007) Migraine: a complex genetic disorder. Lancet Neurol 6(6):521–532
Olesen J et al (2009) Origin of pain in migraine: evidence for peripheral sensitisation. Lancet Neurol 8(7):679–690
Pietrobon D, Striessnig J (2003) Neurobiology of migraine. Nat Rev Neurosci 4(5):386–398
Goadsby PJ et al (2017) Pathophysiology of migraine: a disorder of sensory processing. Physiol Rev 97(2):553–622
Noseda R, Burstein R (2013) Migraine pathophysiology: anatomy of the trigeminovascular pathway and associated neurological symptoms, CSD, sensitization and modulation of pain. Pain 154(Suppl 1):S44–S53
Peck KR, Johnson YL, Smitherman TA (2016) Migraine. Handb Clin Neurol 138:283–293
Gargus JJ (2009) Genetic calcium signaling abnormalities in the central nervous system: seizures, migraine, and autism. Ann N Y Acad Sci 1151:133–156
Shi Y et al (2012) Genetic variation in the calcium/calmodulin-dependent protein kinase (CaMK) pathway is associated with antidepressant response in females. J Affect Disord 136(3):558–566
Torres R et al (2007) Mice genetically deficient in neuromedin U receptor 2, but not neuromedin U receptor 1, have impaired nociceptive responses. Pain 130(3):267–278
Liu JJ et al (2009) Discovery and pharmacological characterization of a small-molecule antagonist at neuromedin U receptor NMUR2. J Pharmacol Exp Ther 330(1):268–275
Han F et al (2016) Paeoniflorin protects diabetic mice against myocardial ischemic injury via the transient receptor potential vanilloid 1/calcitonin gene-related peptide pathway. Cell Biosci 6:37
Xiao Y, Richter JA, Hurley JH (2008) Release of glutamate and CGRP from trigeminal ganglion neurons: role of calcium channels and 5-HT1 receptor signaling. Mol Pain 4:12
Tsuzuki K et al (2004) Menthol-induced Ca2+ release from presynaptic Ca2+ stores potentiates sensory synaptic transmission. J Neurosci 24(3):762–771
Peier AM et al (2002) A TRP channel that senses cold stimuli and menthol. Cell 108(5):705–715
Colburn RW et al (2007) Attenuated cold sensitivity in TRPM8 null mice. Neuron 54(3):379–386
Dussor G, Cao YQ (2016) TRPM8 and migraine. Headache 56(9):1406–1417
de Jong PR et al (2015) TRPM8 on mucosal sensory nerves regulates colitogenic responses by innate immune cells via CGRP. Mucosal Immunol 8(3):491–504
Allain B et al (2012) Neuropilin-1 regulates a new VEGF-induced gene, Phactr-1, which controls tubulogenesis and modulates lamellipodial dynamics in human endothelial cells. Cell Signal 24(1):214–223
Patel RS et al (2012) A genetic risk variant for myocardial infarction on chromosome 6p24 is associated with impaired central hemodynamic indexes. Am J Hypertens 25(7):797–803
Allen PB et al (2004) Phactrs 1-4: a family of protein phosphatase 1 and actin regulatory proteins. Proc Natl Acad Sci U S A 101(18):7187–7192
Winsvold BS et al (2017) Shared genetic risk between migraine and coronary artery disease: a genome-wide analysis of common variants. PLoS One 12(9):e0185663
Kajimura S (2015) Promoting brown and beige adipocyte biogenesis through the PRDM16 pathway. Int J Obes Suppl 5(Suppl 1):S11–S14
Fenzl A, Kiefer FW (2014) Brown adipose tissue and thermogenesis. Horm Mol Biol Clin Invest 19(1):25–37
Seale P et al (2007) Transcriptional control of brown fat determination by PRDM16. Cell Metab 6(1):38–54
AlAmrani A, AbdelKarim M, AlZoghaibi M (2018) PRDM16 gene polymorphism is associated with obesity and blood lipids profiles in Saudi population. J Clin Med 7(6):141
Bigal ME, Lipton RB (2008) Putative mechanisms of the relationship between obesity and migraine progression. Curr Pain Headache Rep 12(3):207–212
Bigal ME et al (2007) Obesity, migraine, and chronic migraine: possible mechanisms of interaction. Neurology 68(21):1851–1861
Zelissen PM et al (1991) Calcitonin gene-related peptide in human obesity. Peptides 12(4):861–863
Peroutka SJ (2004) Migraine: a chronic sympathetic nervous system disorder. Headache 44(1):53–64
Fang J et al (2018) Involvement of the Tetraspanin 2 (TSPAN2) gene in migraine: a case-control study in Han Chinese. Front Neurol 9:714
Wang P et al (2017) Dual-specificity tyrosine-phosphorylation regulated kinase 1A gene transcription is regulated by Myocyte enhancer factor 2D. Sci Rep 7(1):7240
This work was supported by the National Natural Science Foundation of China (Grant No. 81400950, 81501006) and Natural Science Foundation of Liaoning Province (Grant No. 2019-MS-365, 2019-MS-364).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Overall summary of results from meta-analyses of observation studies on risk of migraine (including duplicates, statistically significant and non-significant results).
Overall summary of results from meta-analyses of GWAS on risk of migraine (including duplicates, statistically significant and non-significant results).
Protein-protein interaction network of noteworthy genes related with migraine. The active interaction sources included text mining, experiments, databases as well as co-expression. MCL clustering method was used with the inflation parameter setting at 3.
About this article
Cite this article
Zhao, Y., Zhu, R., Xiao, T. et al. Genetic variants in migraine: a field synopsis and systematic re-analysis of meta-analyses. J Headache Pain 21, 13 (2020). https://doi.org/10.1186/s10194-020-01087-5