Study of association of migraine susceptibility genes with common migraine in 200,000 exome-sequenced UK Biobank participants ============================================================================================================================= * Katherine Alexis Markel * David Curtis ## Abstract **Background** A number of genes have been implicated in rare familial syndromes which have migraine as part of their phenotype but these genes have not previously been implicated in the common form of migraine. **Methods** Among exome-sequenced participants in the UK Biobank we identified 7,194 migraine cases with the remaining 193,433 participants classified as controls. We investigated ten genes previously reported to be implicated in conditions with migraine as a prominent part of the phenotype and carried out gene and variant based tests for association. **Results** We found no evidence for association of these genes or variants with the common form of migraine seen in our subjects. In particular, a frameshift variant in *KCNK18*, F139Wfs*24, which had been shown to segregate with migraine with aura in a multiply affected pedigree was found in 196 (0.10%) controls as well as in 10 (0.14%) cases (χ2 = 0.96, 1 df, p = 0.33). **Conclusions** Since there is no other reported evidence to implicate *KCNK18*, we conclude that this gene and its product, TRESK, should no longer be regarded as being involved in migraine aetiology. Overall, we do not find that rare, functional variants in genes previously implicated to be involved in familial syndromes including migraine as part of the phenotype make a contribution to the commoner forms of migraine observed in this population. ## Introduction Migraine is a disabling primary headache disorder which has a 1 year prevalence of about 12%, peaking in middle age and then reducing thereafter (Lipton et al., 2007). It is characterized by severe pain and is sometimes preceded by neurosensory aberrations known as aura, while associated symptoms include nausea, vomiting, photophobia, phonophobia, tiredness, irritability, reduced concentration, and problems with cognition (Goadsby et al., 2017; Sutherland, Albury and Griffiths, 2019). Migraine may have major socio-economic impacts on individuals’ quality of life, their families, and society due to their greater healthcare needs and reduced work productivity (Becker et al., 2007; Ferrari et al., 2015). Migraine has been shown to have a strong genetic component, with heritability estimates ranging from 34-64%, with common genetic variants with small effect sizes contributing to risk of the common form of the disease while some rare monogenic conditions can include migraine as part of the phenotype (Honkasalo et al., 1995; Russell and Olesen, 1995; Mulder et al., 2003; Bron, Sutherland and Griffiths, 2021). A recent meta-analysis of genome-wide association studies (GWASs) of migraine identified 38 implicated loci (Gormley et al., 2016). Monogenic conditions with migraine as part of the phenotype include familial hemiplegic migraine (FHM), sporadic hemiplegic migraine (SHM), retinal vasculopathy with cerebral leukoencephalopathy and systemic manifestations (RVCL-S), familial advanced sleep-phase syndrome (FASPS), and cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL). In FHM and SHM migraine is a core feature of the phenotype whereas in the other conditions migraine can appear before other symptoms emerge. A number of ion-regulating genes have been implicated in FHM, which is a rare monogenic form of migraine with aura that also includes motor disturbances and hemiparalysis (Goadsby et al., 2017). These include *CACNA1A, ATP1A2, SCN1A*, and possibly *PRRT2*, though there is debate as to whether the latter is disease-causing or modifying (Ophoff et al., 1996; De Fusco et al., 2003; Dichgans et al., 2005; Riant et al., 2012; Pelzer et al., 2014). Variants which have been implicated in FHM differ in their predicted effects, with *ATP1A2* and *PRRT2* having loss-of-function (LOF) variants, *CACNA1A* having gain-of-function (GOF) variants, and *SCN1A* variants having a complex impact on protein function (Pietrobon and Moskowitz, 2013; Sutherland, Albury and Griffiths, 2019). Paroxysmal non-kinesigenic dyskinesia (PNKD) is an autosomal dominant episodic movement disorder caused by variants in the *PNKD* gene (Shen et al., 2011, 2015). There is a report of a family with PNKD accompanied by FHM associated with a frameshift variant in PNKD (Gardiner et al., 2015). SHM is similar to FHM, except that it develops in the absence of a family history of hemiplegic migraine (Kovermann et al., 2017). Variants in *SLC1A3* have been identified in individuals with SHM and are thought to induce the migraine phenotype via a LOF mechanism that disrupts glutamate reuptake and leads to neuronal hyperexcitability (Jen et al., 2005; Kovermann et al., 2017; Sutherland, Albury and Griffiths, 2019; Paucar et al., 2020). With respect to the other conditions, variants in *TREX1* have been implicated in RVCL-S (Richards et al., 2007; Stam et al., 2009), variants in *CSNK1D* in FASPS (Brennan et al., 2013) and variants in *NOTCH3* in CADASIL (Joutel et al., 1996). Two previous studies have specifically looked at *NOTCH3* variants associated with CADASIL in UK Biobank participants. Masoli and colleagues studied the imputed data of 451,424 UKBB participants of European descent for two missense variants predicted to be pathogenic by VEP (Variant Effect Predictor), p.R1231C and p.A1020P, and found that individuals with the former variant had higher diastolic and systolic blood pressure and an increased rate of incident stroke, while individuals with the latter variant only had higher diastolic blood pressure (Masoli et al., 2019). Since almost all pathologic variants in *NOTCH3* involve alterations in the number of cysteines in the protein’s epidermal growth factor-like repeat domains, Rutten and colleagues examined these variants in 50,000 UKBB participants and found that these factors were associated with a very broad range of phenotypes from normal to CADASIL, with migraines often presenting decades before other symptoms (Rutten et al., 2020). A single pedigree has been reported in which migraine with aura segregated with a frameshift variant (F139Wfs*24) in the *KCNK18* gene which encodes the TRESK two-pore potassium channel protein (Lafrenière et al., 2010). Subsequent studies suggested another frameshift variant (Y121Lfs*44) may impair protein function in a similar manner (Royal et al., 2019). These variants produce an additional start codon (ATG) within the reading frame, which may lead to the production of TRESK protein fragments that are thought to downregulate *KCNK18* expression via a feedback loop and also interfere with normal channel activity in a dominant negative fashion (Lafrenière et al., 2010; Lafrenière and Rouleau, 2011, 2012; Andres-Enguix et al., 2012; Royal et al., 2019). We sought to investigate whether variants impacting the function of these ten genes might contribute to risk of developing commoner, less familial forms of migraine in the general population. In order to do this we investigated whether functional variants were associated with a clinically defined migraine phenotype in 200,000 exome-sequenced UK Biobank participants. ## Methods The approach used was similar to that previously described for other phenotypes (Curtis, 2021). The UK Biobank consists of 500,000 volunteers who have undergone extensive phenotyping and who have provided biological samples. They are on average somewhat older and more healthy than the British population as a whole. Exome sequence data has been released for 200,627 participants and these variant call files were downloaded after genotype-calling by the UK Biobank Exome Sequencing Consortium using the GRCh38 assembly with coverage 20X at 95.6% of sites on average (Szustakowski et al., 2020). UK Biobank had obtained ethics approval from the North West Multi-Centre Research Ethics Committee which covers the UK (approval number: 11/NW/0382) and had obtained written informed consent from all participants. The UK Biobank approved an application for use of the data (ID 51119) and ethics approval for the analyses was obtained from the UCL Research Ethics Committee (11527/001). All variants were annotated using the standard software packages VEP, PolyPhen and SIFT (Kumar, Henikoff and Ng, 2009; Adzhubei, Jordan and Sunyaev, 2013; McLaren et al., 2016). As described previously, population principal components were obtained using version 2.0 of *plink* ([https://www.cog-gemonics.org/plink/2.0/](https://www.cog-gemonics.org/plink/2.0/)) with the *commands –maf 0*.*1 –pca 20 approx* (Chang et al., 2015; Galinsky et al., 2016; Curtis, 2021). In order to define participants as cases we used a similar approach as that described previously and combined information about recorded diagnoses and about medication (Curtis, 2020a). Participants who had an ICD-10 diagnosis code of G43.* and/or who reported taking one of a number of medications which are relatively specifically indicated as treatments for migraine were classified as cases while the remainder were taken to be controls. The list of the medications used to define cases is provided in Supplementary Table S1. Tests for association were carried out at the level of the gene, of categories of variant within each gene and at the level of individual variants. In order to test for association at the level of the gene, weighted burden analysis was carried out using the SCOREASSOC program (Curtis, 2012, 2019). For each variant a weight was assigned according to its predicted effect on gene function, with variants with a more severe impact being allocated a higher weight so that, for example, stop gained variants were assigned a weight of 100 while missense variants were assigned a weight of 5. The full set of variant types and weights is presented in Supplementary Table S2. Attention was restricted to rare variants with minor allele frequency (MAF) <= 0.01 in both cases and controls. As previously described, variants were also weighted by MAF so that variants with MAF=0.01 were given a weight of 1 while very rare variants with MAF close to zero were given a weight of 10 (Curtis, 2020b). For each variant the functional weight was multiplied by the frequency weight to produce an overall weight and then for each subject the weights of the variants carried by that subject were summed to produce a weighted burden score which was included in a logistic regression analysis with sex and 20 principal components as covariates. We have previously shown that this process adequately controls test statistic inflation in this population (Curtis, 2021). The statistical significance for association between migraine and the weighted burden score for each gene was summarised as a signed log p value (SLP), which is the log base 10 of the p value given a positive sign if the score is higher in cases and negative if it is higher in controls. The variant types were also grouped into broader categories such as intronic, splice site, protein altering, as shown in Supplemental Table S2, in order to test whether any particular category of variant within a gene was associated with migraine risk. As described previously, logistic regression analyses were performed using the counts of the separate categories of variant as predictor variables, again including principal components and sex as covariates, to estimate the effect size for each category (Curtis, 2021). The odds ratios associated with each category were estimated along with their standard errors and the Wald statistic was used to obtain a p value, except for categories in which variants occurred fewer than 50 times in which case Fisher’s exact test was applied to the variant counts. The associated p value was converted to an SLP, again with the sign being positive if the mean count was higher in cases than controls. We also examined specific missense variants which had previously been reported as implicated in the migraine-associated syndromes listed above. After a thorough literature review, we compiled a list of 222 such variants in the ten genes included for this study, as listed in Supplementary Table S3. For variants observed in 20 or more subjects, counts were compared between cases and controls. For rarer variants, the pooled counts for each gene were compared. Additionally, because of the report implicating a *KCNK18* frameshift variant we compared the counts of all frameshift variants in this gene individually and collectively. Data manipulation and statistical analyses were performed using GENEVARASSOC, SCOREASSOC and R (R Core Team, 2014). ## Results Using the process described to define cases based on assigned ICD10 diagnosis and/or migraine specific medication, there were 7,194 cases, of whom 78.2% were female with mean age 55.2 (SD=7.9), and 193,433 controls, of whom 54.2% were female with mean age 56.5 (SD=8.1). Gene-wise weighted burden analysis did not find evidence for association of rare, functional variants with migraine for any of the ten genes tested. Table 1 shows the SLP obtained for each gene of which the highest is 1.58 for *ATP1A2*, equivalent to *p*=0.026. This result is not statistically significant after correction for multiple testing and none of the SLPs for the remaining genes approaches significance. View this table: [Table 1.](http://medrxiv.org/content/early/2021/12/03/2021.12.02.21267184/T1) Table 1. Results of weighted burden analysis testing for association with migraine. Statistical significance is expressed as the signed log *p* value (SLP), with positive values indicating that rate, functional variants in the gene are associated with increased risk of migraine. Likewise, there was no category of variants within any gene that produced significant evidence of association with the migraine phenotype. Detailed results are provided in Supplementary Table S4. Of the individual missense variants identified from the literature, 13 were observed in 20 or more participants. Results for these variants are shown in Table 2 and it can be seen that none of them is markedly more common in cases than controls. E492K in *ATP1A2* occurs in 18 cases versus 298 controls, a results which is nominally significant at p=0.04 but again this is not statistically significant after correction for multiple testing. The cumulative totals for rarer variants in each gene are shown in Table 3 and again it can be seen that these are not commoner among cases. View this table: [Table 2.](http://medrxiv.org/content/early/2021/12/03/2021.12.02.21267184/T2) Table 2. Counts of missense variants previously reported to be implicated in migraine phenotypes for variants observed in 20 or more participants. View this table: [Table 3.](http://medrxiv.org/content/early/2021/12/03/2021.12.02.21267184/T3) Table 3. Cumulative counts of missense variants previously reported to be implicated in migraine phenotypes for variants observed in fewer than 20 participants. The frameshift variant in *KCNK18* which had been reported to segregate with migraine with aura, F139Wfs*24, was found in 196 (0.10%) controls and 10 (0.14%) cases (χ2 = 0.96, p = 0.33). The Y121Lfs*44 frameshift variant was found in 205 (0.11%) controls and 6 (0.08%) cases (χ2 = 0.34, p = 0.56). We identified an additional 20 controls and 0 additional cases with other frameshift variants in *KCNK18*, bringing the overall numbers to 421 (0.22%) controls and 16 (0.22%) cases (χ2 = 0.01, p = 0.93). ## Discussion Overall, we do not observe an association between migraine risk and an excess of rare variants predicted to impact function in any of the ten genes previously implicated in rare familial disorders with migraine as part of the phenotype. This applies to the three canonical FHM genes (*CACNA1A, ATP1A2*, and *SCN1A*) and a potential fourth FHM gene *(PRRT2*), as well as the other genes reported to cause syndromes in which migraine may be a feature. These findings contrast with those we obtained from a similar study of hyperlipidaemia in the same UK Biobank sample, in which we found that for genes known to cause rare, severe forms of familial hypercholesterolaemia there were very rare variants with large effects on risk but also less rare variants which made a wider contribution to hyperlipidaemia in the general population (Curtis, 2021). Even when we focus on specific variants previously reported to be associated with migraine syndromes we do not observe a significant excess among our cases. Of course, the method we have used to define caseness is not as rigorous as could be applied in a specially designed study and the control samples will include some participants with migraine who have not been assigned a formal diagnosis and who are not taking a migraine-specific medication. Nevertheless, it is also possible that not all of the previously reported variants do in fact raise migraine risk. Our findings for *KCNK18* are particularly noteworthy. The only evidence to support the claim that this gene has an aetiological role in migraine is the observation that a frameshift variant, F139Wfs*24, perfectly segregated with the phenotype of migraine with aura in a single large pedigree (Lafrenière et al., 2010). It has been the subject of several functional studies, most of which were conducted on cultured cells or in animal models (Lafrenière et al., 2010; Lafrenière and Rouleau, 2011; Andres-Enguix et al., 2012; Liu et al., 2013; Guo et al., 2014, 2019; Royal et al., 2019). However we cannot find any report that this variant or other functional variants in *KCNK18* have been observed in other migraine patients. We observe this specific variant in 196 controls as well as in 10 cases and overall frameshift variants in this gene occur at equal frequency in controls and cases. These findings are not compatible with the notion that F139Wfs*24 causes a very rare quasi-mendelian form of migraine and hence we conclude that *KCNK18* and its product TRESK should no longer be regarded as influencing migraine susceptibility. We note that the haplotype segregating with migraine in the original pedigree extends over 10 megabases and harbours dozens of genes. While we can appreciate that a frameshift variant would have appeared to be a plausible candidate, our new results rule out a role for this variant and we speculate that perhaps another gene within this haplotype is actually responsible for the migraine cases observed in the pedigree. Using an exome-sequenced sample of 200,000 participants we have been unable to identify a role for previously identified genes implicated in rare familial disorders as contributing to the risk of the more common migraine phenotype. Sequence data for the remaining 300,000 participants will be made available in due course and with the increased sample size it will become feasible to systematically study all genes in order to attempt to identify some which do influence migraine risk. ## Data Availability Raw data is available on application from UK Biobank. Relevant derived variables will be deposited in UK Biobank. [https://www.ukbiobank.ac.uk/](https://www.ukbiobank.ac.uk/) ## Declaration of conflicting interests The authors declare no conflicting interests. ## Contributions DC conceived the methods, developed the programs, and contributed to the analysis and writing of this paper. KM helped adapt the methods for the migraine phenotype, selected the genes and variants included in the study, contributed to the analysis of the data, and wrote the paper. ## Data Availability Statment Raw data is available on application from UK Biobank. Relevant derived variables will be deposited in UK Biobank. The scripts and software used to carry out the analyses are available at: [https://github.com/davenomiddlenamecurtis](https://github.com/davenomiddlenamecurtis). ## Acknowledgements The authors wish to acknowledge the staff of University College London’s Computer Science Department and their High Performance Computing Cluster. We would also like to thank the participants of the UK Biobank study who provided the data for this study. This work was carried out in part using resources provided by BBSRC equipment grant BB/R01356X/1. KM would like to thank The Migraine Trust for awarding her a Susan Haydon Bursary. ## Footnotes * There was a typo in the name of a key variant in the abstract. I had F138Wfs*24 instead of F139Wfs*24. As this variant has previously been reported in the literature and is key to our findings, I think it's important that we name it correctly in the abstract. Apologies for this error. * Received December 2, 2021. * Revision received December 3, 2021. * Accepted December 3, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. Adzhubei, I., Jordan, D. M. and Sunyaev, S. R. (2013) ‘Predicting functional effect of human missense mutations using PolyPhen-2.’, Current protocols in human genetics. NIH Public Access, 7 Unit7.20. doi: 10.1002/0471142905.hg0720s76. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/0471142905.hg0720s76&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23315928&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 2. Andres-Enguix, I. et al. (2012) ‘Functional analysis of missense variants in the TRESK (KCNK18) K + channel’, Scientific Reports. Nature Publishing Group, 2(1), pp. 1–7. doi: 10.1038/srep00237. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/srep00237&link_type=DOI) 3. Becker, C. et al. (2007) ‘Migraine Incidence, Comorbidity and Health Resource Utilization in the UK’:, [https://doi.org/10.1111/j.1468-2982.2007.01469.x](https://doi.org/10.1111/j.1468-2982.2007.01469.x). SAGE PublicationsSage UK: London, England, 28(1), xpp. 57–64. doi: 10.1111/J.1468-2982.2007.01469.X. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/J.1468-2982.2007.01469.X&link_type=DOI) 4. Brennan, K. C. et al. (2013) ‘Casein kinase Iδ mutations in familial migraine and advanced sleep phase’, Science Translational Medicine. NIH Public Access, 5(183), p. 183ra56. doi: 10.1126/scitranslmed.3005784. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6InNjaXRyYW5zbWVkIjtzOjU6InJlc2lkIjtzOjEzOiI1LzE4My8xODNyYTU2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMDMvMjAyMS4xMi4wMi4yMTI2NzE4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 5. Bron, C., Sutherland, H. G. and Griffiths, L. R. (2021) ‘Exploring the Hereditary Nature of Migraine’, Neuropsychiatric Disease and Treatment, Volume 17, pp. 1183–1194. doi: 10.2147/NDT.S282562. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2147/NDT.S282562&link_type=DOI) 6. Chang, C. C. et al. (2015) ‘Second-generation PLINK: rising to the challenge of larger and richer datasets’, GigaScience. BioMed Central, 4(1), p. 7. doi: 10.1186/s13742-015-0047-8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13742-015-0047-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25722852&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 7. Curtis, D. (2012) ‘A rapid method for combined analysis of common and rare variants at the level of a region, gene, or pathway’, Adv Appl Bioinform Chem, 5, pp. 1–9. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22888262&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 8. Curtis, D. (2019) ‘A weighted burden test using logistic regression for integrated analysis of sequence variants, copy number variants and polygenic risk score’, European Journal of Human Genetics, 27(1), pp. 114–124. doi: 10.1038/s41431-018-0272-6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41431-018-0272-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 9. Curtis, D. (2020a) ‘Analysis of exome-sequenced UK Biobank subjects implicates genes affecting risk of hyperlipidaemia’, Molecular Genetics and Metabolism. Academic Press, 131, pp. 277–283. doi: 10.1016/j.ymgme.2020.07.009. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ymgme.2020.07.009&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 10. Curtis, D. (2020b) ‘Multiple Linear Regression Allows Weighted Burden Analysis of Rare Coding Variants in an Ethnically Heterogeneous Population’, Human Heredity. Karger Publishers, pp. 1–10. doi: 10.1159/000512576. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000512576&link_type=DOI) 11. Curtis, D. (2021) ‘Analysis of 200 000 exome-sequenced UK Biobank subjects illustrates the contribution of rare genetic variants to hyperlipidaemia’, Journal of Medical Genetics. BMJ Publishing Group Ltd, p. jmedgenet-2021-107752. doi: 10.1136/jmedgenet-2021-107752. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToiam1lZGdlbmV0IjtzOjU6InJlc2lkIjtzOjIzOiJqbWVkZ2VuZXQtMjAyMS0xMDc3NTJ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzAzLzIwMjEuMTIuMDIuMjEyNjcxODQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 12. Dichgans, M.⍰; et al. (2005) ‘Mutation in the neuronal voltage-gated sodium channel SCN1A in familial hemiplegic migraine’, The Lancet, 366(9483), p. 371. 13. Ferrari, M. D. et al. (2015) ‘Migraine pathophysiology: lessons from mouse models and human genetics’, The Lancet Neurology. Elsevier, 14(1), pp. 65–80. doi: 10.1016/S1474-4422(14)70220-0. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1474-4422(14)70220-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25496898&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 14. De Fusco, M. et al. (2003) ‘Haploinsufficiency of ATP1A2 encoding the Na+/K+ pump α2 subunit associated with familial hemiplegic migraine type 2’, Nature Genetics. Nature Publishing Group, 33(2), pp. 192–196. doi: 10.1038/ng1081. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng1081&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12539047&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000180773700019&link_type=ISI) 15. Galinsky, K. J. et al. (2016) ‘Fast Principal-Component Analysis Reveals Convergent Evolution of ADH1B in Europe and East Asia’, American Journal of Human Genetics. Cell Press, 98(3), pp. 456–472. doi: 10.1016/j.ajhg.2015.12.022. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2015.12.022&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26924531&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 16. Gardiner, A. R. et al. (2015) ‘The clinical and genetic heterogeneity of paroxysmal dyskinesias’, Brain. Oxford University Press, 138(12), pp. 3567–3580. doi: 10.1093/brain/awv310. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awv310&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26598494&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 17. Goadsby, P. J. et al. (2017) ‘Pathophysiology of migraine: A disorder of sensory processing’, Physiological Reviews. American Physiological Society, 97(2), pp. 553–622. doi: 10.1152/physrev.00034.2015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/physrev.00034.2015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28179394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 18. Gormley, P. et al. (2016) ‘Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine’, Nature Genetics. Nature Research, 48(8), pp. 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%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 19. Guo, Z. et al. (2014) ‘Nonmigraine-associated TRESK K+ channel variant C110R does not increase the excitability of trigeminal ganglion neurons’, [https://doi.org/10.1152/jn.00267.2014](https://doi.org/10.1152/jn.00267.2014). American Physiological Society Bethesda, MD, 112(3), xpp. 568–579. doi: 10.1152/JN.00267.2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/JN.00267.2014&link_type=DOI) 20. Guo, Z. et al. (2019) ‘TRESK K+ Channel Activity Regulates Trigeminal Nociception and Headache’, eNeuro. Society for Neuroscience, 6(4). doi: 10.1523/ENEURO.0236-19.2019. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZW5ldXJvIjtzOjU6InJlc2lkIjtzOjIzOiI2LzQvRU5FVVJPLjAyMzYtMTkuMjAxOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzAzLzIwMjEuMTIuMDIuMjEyNjcxODQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 21. Honkasalo, M. -L et al. (1995) ‘Migraine and Concomitant Symptoms Among 8167 Adult Twin Pairs’, Headache: The Journal of Head and Face Pain, 35(2), pp. 70–78. doi: 10.1111/J.1526-4610.1995.HED3502070.X. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/J.1526-4610.1995.HED3502070.X&link_type=DOI) 22. Jen, J. C. et al. (2005) ‘Mutation in the glutamate transporter EAAT1 causes episodic ataxia, hemiplegia, and seizures’, Neurology. Lippincott Williams and Wilkins, 65(4), pp. 529–534. doi: 10.1212/01.WNL.0000172638.58172.5A. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/01.WNL.0000172638.58172.5a&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16116111&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 23. Joutel, A. et al. (1996) ‘Notch3 mutations in CADASIL, a hereditary adult-onset condition causing stroke and dementia’, Nature 1996 383:6602. Nature Publishing Group, 383(6602), pp. 707–710. doi: 10.1038/383707a0. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/383707a0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8878478&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VN91800050&link_type=ISI) 24. Kovermann, P. et al. (2017) ‘Impaired K+ binding to glial glutamate transporter EAAT1 in migraine’, Scientific Reports. Nature Publishing Group, 7(1), pp. 1–10. doi: 10.1038/s41598-017-14176-4. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-017-14176-4&link_type=DOI) 25. Kumar, P., Henikoff, S. and Ng, P. C. (2009) ‘Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm’, Nature Protocols, 4(8), pp. 1073–1081. doi: 10.1038/nprot.2009.86. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nprot.2009.86&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19561590&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268858700007&link_type=ISI) 26. Lafrenière, R. G. et al. (2010) ‘A dominant-negative mutation in the TRESK potassium channel is linked to familial migraine with aura’, Nature Medicine. Nature Publishing Group, 16(10), pp. 1157–1160. doi: 10.1038/nm.2216. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm.2216&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20871611&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000282644800048&link_type=ISI) 27. Lafrenière, R. G. and Rouleau, G. A. (2011) ‘Migraine: Role of the TRESK two-pore potassium channel’, The International Journal of Biochemistry & Cell Biology. Pergamon, 43(11), pp. 1533–1536. doi: 10.1016/J.BIOCEL.2011.08.002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.BIOCEL.2011.08.002&link_type=DOI) 28. Lafrenière, R. G. and Rouleau, G. A. (2012) ‘Identification of novel genes involved in migraine.’, Headache. John Wiley & Sons, Ltd, 52 Suppl 2, pp. 107–110. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1526-4610.2012.02237.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23030542&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000309447000012&link_type=ISI) 29. Lipton, R. B. et al. (2007) ‘Migraine prevalence, disease burden, and the need for preventive therapy’, Neurology. Wolters Kluwer Health, Inc. on behalf of the American Academy of Neurology, 68(5), pp. 343–349. doi: 10.1212/01.WNL.0000252808.97649.21. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/01.wnl.0000252808.97649.21&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17261680&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 30. Liu, P. et al. (2013) ‘Functional Analysis of a Migraine-Associated TRESK K+ Channel Mutation’, Journal of Neuroscience. Society for Neuroscience, 33(31), pp. 12810–12824. doi: 10.1523/JNEUROSCI.1237-13.2013. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjExOiIzMy8zMS8xMjgxMCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzAzLzIwMjEuMTIuMDIuMjEyNjcxODQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 31. Masoli, J. A. H. et al. (2019) ‘Clinical Outcomes of CADASIL-Associated NOTCH3 Mutations in 451,424 European Ancestry Community Volunteers’, Translational Stroke Research. Springer, 10(4), p. 339. doi: 10.1007/S12975-018-0671-6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/S12975-018-0671-6&link_type=DOI) 32. McLaren, W. et al. (2016) ‘The Ensembl Variant Effect Predictor.’, Genome biology, 17(1), p. 122. doi: 10.1186/s13059-016-0974-4. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-016-0974-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27268795&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 33. Mulder, E. J. et al. (2003) ‘Genetic and Environmental Influences on Migraine: A Twin Study Across Six Countries’, Twin Research and Human Genetics. Cambridge University Press, 6(5), pp. 422–431. doi: 10.1375/TWIN.6.5.422. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1375/TWIN.6.5.422&link_type=DOI) 34. Ophoff, R. A. et al. (1996) ‘Familial Hemiplegic Migraine and Episodic Ataxia Type-2 Are Caused by Mutations in the Ca2+ Channel Gene CACNL1A4’, Cell. Elsevier, 87(3), pp. 543–552. doi: 10.1016/S0092-8674(00)81373-2. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0092-8674(00)81373-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8898206&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VQ46600019&link_type=ISI) 35. Paucar, M. et al. (2020) ‘SLC1A3 variant associated with hemiplegic migraine and acetazolamide-responsive MRS changes’, Neurology: Genetics. Lippincott Williams and Wilkins, 6(4). doi: 10.1212/NXG.0000000000000474. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoibm5nIjtzOjU6InJlc2lkIjtzOjg6IjYvNC9lNDc0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMDMvMjAyMS4xMi4wMi4yMTI2NzE4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 36. Pelzer, N. et al. (2014) ‘PRRT2 and hemiplegic migraine: A complex association’, Neurology. Lippincott Williams and Wilkins, 83(3), pp. 288–290. doi: 10.1212/WNL.0000000000000590. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0000000000000590&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24928127&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 37. Pietrobon, D. and Moskowitz, M. A. (2013) ‘Pathophysiology of migraine’, Annual Review of Physiology, pp. 365–391. doi: 10.1146/annurev-physiol-030212-183717. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-physiol-030212-183717&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23190076&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000316381400017&link_type=ISI) 38. R Core Team (2014) R: A language and environment for statistical computing. Vienna, Austria.: R Foundation for Statistical Computing. 39. Riant, F. et al. (2012) ‘PRRT2 mutations cause hemiplegic migraine’, Neurology. Lippincott Williams and Wilkins, 79(21), pp. 2122–2124. doi: 10.1212/WNL.0b013e3182752cb8. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibmV1cm9sb2d5IjtzOjU6InJlc2lkIjtzOjEwOiI3OS8yMS8yMTIyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMDMvMjAyMS4xMi4wMi4yMTI2NzE4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 40. Richards, A. et al. (2007) ‘C-terminal truncations in human 31-51 DNA exonuclease TREX1 cause autosomal dominant retinal vasculopathy with cerebral leukodystrophy’, Nature Genetics. Nature Publishing Group, 39(9), pp. 1068–1070. doi: 10.1038/ng2082. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng2082&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17660820&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249122400014&link_type=ISI) 41. Royal, P. et al. (2019) ‘Migraine-Associated TRESK Mutations Increase Neuronal Excitability through Alternative Translation Initiation and Inhibition of TREK’, Neuron, 101(2), pp. 232–245.e6. doi: 10.1016/j.neuron.2018.11.039. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2018.11.039&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30573346&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) 42. Russell, M. B. and Olesen, J. (1995) ‘Increased familial risk and evidence of genetic factor in migraine’, BMJ. British Medical Journal Publishing Group, 311(7004), pp. 541–544. doi: 10.1136/BMJ.311.7004.541. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjEyOiIzMTEvNzAwNC81NDEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMi8wMy8yMDIxLjEyLjAyLjIxMjY3MTg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 43. Rutten, J. W. et al. (2020) ‘Broad phenotype of cysteine-altering NOTCH3 variants in UK Biobank: CADASIL to nonpenetrance’, Neurology, 95(13), pp. e1835–e1843. doi: 10.1212/WNL.0000000000010525. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibmV1cm9sb2d5IjtzOjU6InJlc2lkIjtzOjExOiI5NS8xMy9lMTgzNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzAzLzIwMjEuMTIuMDIuMjEyNjcxODQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 44. Shen, Y. et al. (2011) ‘Mutations in PNKD causing paroxysmal dyskinesia alters protein cleavage and stability’, Human Molecular Genetics. Oxford Academic, 20(12), pp. 2322–2332. doi: 10.1093/HMG/DDR125. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddr125&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21487022&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000290849200003&link_type=ISI) 45. Shen, Y. et al. (2015) ‘Protein mutated in paroxysmal dyskinesia interacts with the active zone protein RIM and suppresses synaptic vesicle exocytosis’, Proceedings of the National Academy of Sciences. National Academy of Sciences, 112(10), pp. 2935–2941. doi: 10.1073/PNAS.1501364112. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTEyLzEwLzI5MzUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMi8wMy8yMDIxLjEyLjAyLjIxMjY3MTg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 46. Stam, A. H. et al. (2009) ‘Migraine and genetic and acquired vasculopathies’, Cephalalgia. SAGE PublicationsSage UK: London, England, pp. 1006–1017. doi: 10.1111/j.1468-2982.2009.01940.x. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1468-2982.2009.01940.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19689610&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F03%2F2021.12.02.21267184.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268683800014&link_type=ISI) 47. Sutherland, H. G., Albury, C. L. and Griffiths, L. R. (2019) ‘Advances in genetics of migraine’, Journal of Headache and Pain. BioMed Central Ltd., pp. 1–20. doi: 10.1186/s10194-019-1017-9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s10194-019-1017-9&link_type=DOI) 48. Szustakowski, J. D. et al. (2020) ‘Advancing Human Genetics Research and Drug Discovery through Exome Sequencing of the UK Biobank’, medRxiv, p. 2020.11.02.20222232. doi: 10.1101/2020.11.02.20222232. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4xMS4wMi4yMDIyMjIzMnYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMDMvMjAyMS4xMi4wMi4yMTI2NzE4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=)