Genome-wide analysis of 102,084 migraine cases identifies 123 risk loci and subtype-specific risk alleles.

Migraine affects over a billion individuals worldwide but its genetic underpinning remains largely unknown. Here, we performed a genome-wide association study of 102,084 migraine cases and 771,257 controls and identified 123 loci, of which 86 are previously unknown. These loci provide an opportunity to evaluate shared and distinct genetic components in the two main migraine subtypes: migraine with aura and migraine without aura. Stratification of the risk loci using 29,679 cases with subtype information indicated three risk variants that seem specific for migraine with aura (in HMOX2, CACNA1A and MPPED2), two that seem specific for migraine without aura (near SPINK2 and near FECH) and nine that increase susceptibility for migraine regardless of subtype. The new risk loci include genes encoding recent migraine-specific drug targets, namely calcitonin gene-related peptide (CALCA/CALCB) and serotonin 1F receptor (HTR1F). Overall, genomic annotations among migraine-associated variants were enriched in both vascular and central nervous system tissue/cell types, supporting unequivocally that neurovascular mechanisms underlie migraine pathophysiology.

Migraine, Stroke, and Cervical Arterial Dissection: Shared Genetics for a Triad of Brain Disorders With Vascular Involvement.

BACKGROUND AND OBJECTIVES: Migraine, stroke, and cervical artery dissection (CeAD) represent a triad of cerebrovascular disorders with pairwise comorbid relationships and vascular involvement. Larger samples and recent advances in methodology invite systematic exploration of their shared genetics. METHODS: Genetic analyses leveraged summary statistics from genome-wide association studies of the largest available samples of each disorder, including subtypes of stroke (ischemic stroke, large artery stroke, small vessel stroke, and cardioembolic stroke) and migraine (with aura and without aura). For each pair of disorders, genetic correlation was assessed both on a genome-wide basis and within independent segments across the genome including known specific loci for each disorder. A cross-trait meta-analysis was used to identify novel candidate loci. Finally, potential causality of migraine susceptibility on stroke and CeAD was assessed by Mendelian randomization. RESULTS: Among all pairs of disorders, genome-wide genetic correlation was observed only between CeAD and migraine, particularly MO. Local genetic correlations were more extensive between migraine and CeAD than those between migraine and stroke or CeAD and stroke and revealed evidence for novel CeAD associations at rs6693567 (ADAMTSL4/ECM1), rs11187838 (PLCE1), and rs7940646 (MRVI1) while strengthening prior subthreshold evidence at rs9486725 (FHL5) and rs650724 (LRP1). At known migraine loci, novel associations with stroke had concordant risk alleles for small vessel stroke at rs191602009 (CARF) and for cardioembolic stroke at rs55884259 (NKX2-5). Known migraine loci also revealed novel associations but with opposite risk alleles for all stroke, ischemic stroke, and small vessel stroke at rs55928386 (HTRA1), for large artery stroke at rs11172113 (LRP1), and for all stroke and ischemic stroke at rs1535791 and rs4942561 (both LRCH1), respectively. rs182923402 (near PTCH1) was a novel concordant locus for migraine and cardioembolic stroke. Mendelian randomization supported potential causal influences of migraine on CeAD (odds ratio [95% confidence interval] per doubling migraine prevalence = 1.69 [1.24-2.3], p = 0.0009) with concordant risk, but with opposite risk on large artery stroke (0.86 [0.76-0.96], p = 0.0067). DISCUSSION: The findings emphasize shared genetic risk between migraine and CeAD while identifying loci with likely vascular function in migraine and shared but opposite genetic risk between migraine and stroke subtypes, and a central role of LRP1 in all 3 cerebrovascular disorders.

Polygenic risk provides biological validity for the ICHD-3 criteria among Finnish migraine families.

BACKGROUND: Migraine is diagnosed using the extensively field-tested International Classification of Headache Disorders (ICHD-3) consensus criteria derived by the International Headache Society. To evaluate the criteria in respect to a measurable biomarker, we studied the relationship between the main ICHD-3 criteria and the polygenic risk score, a measure of common variant burden in migraine. METHODS: We used linear mixed models to study the correlation of ICHD-3 diagnostic criteria, underlying symptoms, and main diagnoses with the polygenic risk score of migraine in a cohort of 8602 individuals from the Finnish Migraine Genome Project. RESULTS: Main diagnostic categories and all underlying diagnostic criteria formed a consistent continuum along the increasing polygenic burden. Polygenic risk was associated with the heterogeneous clinical picture starting from the non-migraine headache (mean 0.07; 95% CI 0.02-0.12; p = 0.008 compared to the non-headache group), to probable migraine (mean 0.13; 95% CI 0.08-0.18; p < 0.001), migraine headache (mean 0.17; 95% CI 0.14-0.21; p < 0.001) and migraine with typical visual aura (mean 0.29; 95% CI 0.26-0.33; p < 0.001), all the way to the hemiplegic aura (mean 0.37; 95% CI 0.31-0.43; p < 0.001). All individual ICHD-3 symptoms and the total number of reported symptoms, a surrogate of migraine complexity, demonstrated a clear inclination with an increasing polygenic risk. CONCLUSIONS: The complex migraine phenotype progressively follows the polygenic burden from individuals with no headache to non-migrainous headache and up to patients with attacks manifesting all the features of the ICHD-3 headache and aura. Results provide further biological support for the ICHD-3 diagnostic criteria.

Phenotypic and Genotypic Associations Between Migraine and Lipoprotein Subfractions.

BACKGROUND AND OBJECTIVE: To evaluate phenotypic and genetic relationships between migraine and lipoprotein subfractions. METHODS: We evaluated phenotypic associations between migraine and 19 lipoprotein subfraction measures in the Women's Genome Health Study (n = 22,788). We then investigated genetic relationships between these traits using summary statistics from the International Headache Genetics Consortium for migraine (ncase = 54,552, ncontrol = 297,970) and combined summary data for lipoprotein subfractions (n up to 47,713). RESULTS: There was a significant phenotypic association (odds ratio 1.27 [95% confidence interval 1.12-1.44]) and a significant genetic correlation at 0.18 (p = 0.001) between migraine and triglyceride-rich lipoproteins (TRLPs) concentration but not for low-density lipoprotein or high-density lipoprotein subfractions. Mendelian randomization (MR) estimates were largely null, implying that pleiotropy rather than causality underlies the genetic correlation between migraine and lipoprotein subfractions. Pleiotropy was further supported in cross-trait meta-analysis, revealing significant shared signals at 4 loci (chr2p21 harboring THADA, chr5q13.3 harboring HMGCR, chr6q22.31 harboring HEY2, and chr7q11.23 harboring MLXIPL) between migraine and lipoprotein subfractions. Three of these loci were replicated for migraine (p < 0.05) in a smaller sample from the UK Biobank. The shared signal at chr5q13.3 colocalized with expression of HMGCR, ANKDD1B, and COL4A3BP in multiple tissues. CONCLUSIONS: The study supports the association between certain lipoprotein subfractions, especially for TRLP, and migraine in populations of European ancestry. The corresponding shared genetic components may help identify potential targets for future migraine therapeutics. CLASSIFICATION OF EVIDENCE: This study provides Class I evidence that migraine is significantly associated with some lipoprotein subfractions.

Gene family information facilitates variant interpretation and identification of disease-associated genes in neurodevelopmental disorders.

BACKGROUND: Classifying pathogenicity of missense variants represents a major challenge in clinical practice during the diagnoses of rare and genetic heterogeneous neurodevelopmental disorders (NDDs). While orthologous gene conservation is commonly employed in variant annotation, approximately 80% of known disease-associated genes belong to gene families. The use of gene family information for disease gene discovery and variant interpretation has not yet been investigated on a genome-wide scale. We empirically evaluate whether paralog-conserved or non-conserved sites in human gene families are important in NDDs. METHODS: Gene family information was collected from Ensembl. Paralog-conserved sites were defined based on paralog sequence alignments; 10,068 NDD patients and 2078 controls were statistically evaluated for de novo variant burden in gene families. RESULTS: We demonstrate that disease-associated missense variants are enriched at paralog-conserved sites across all disease groups and inheritance models tested. We developed a gene family de novo enrichment framework that identified 43 exome-wide enriched gene families including 98 de novo variant carrying genes in NDD patients of which 28 represent novel candidate genes for NDD which are brain expressed and under evolutionary constraint. CONCLUSION: This study represents the first method to incorporate gene family information into a statistical framework to interpret variant data for NDDs and to discover new NDD-associated genes.

Contribution of rare and common variants to intellectual disability in a sub-isolate of Northern Finland.

The contribution of de novo variants in severe intellectual disability (ID) has been extensively studied whereas the genetics of mild ID has been less characterized. To elucidate the genetics of milder ID we studied 442 ID patients enriched for mild ID (>50%) from a population isolate of Finland. Using exome sequencing, we show that rare damaging variants in known ID genes are observed significantly more often in severe (27%) than in mild ID (13%) patients. We further observe a significant enrichment of functional variants in genes not yet associated with ID (OR: 2.1). We show that a common variant polygenic risk significantly contributes to ID. The heritability explained by polygenic risk score is the highest for educational attainment (EDU) in mild ID (2.2%) but lower for more severe ID (0.6%). Finally, we identify a Finland enriched homozygote variant in the CRADD ID associated gene.

Analysis of shared heritability in common disorders of the brain.

Disorders of the brain can exhibit considerable epidemiological comorbidity and often share symptoms, provoking debate about their etiologic overlap. We quantified the genetic sharing of 25 brain disorders from genome-wide association studies of 265,218 patients and 784,643 control participants and assessed their relationship to 17 phenotypes from 1,191,588 individuals. Psychiatric disorders share common variant risk, whereas neurological disorders appear more distinct from one another and from the psychiatric disorders. We also identified significant sharing between disorders and a number of brain phenotypes, including cognitive measures. Further, we conducted simulations to explore how statistical power, diagnostic misclassification, and phenotypic heterogeneity affect genetic correlations. These results highlight the importance of common genetic variation as a risk factor for brain disorders and the value of heritability-based methods in understanding their etiology.

Common Variant Burden Contributes to the Familial Aggregation of Migraine in 1,589 Families.

Complex traits, including migraine, often aggregate in families, but the underlying genetic architecture behind this is not well understood. The aggregation could be explained by rare, penetrant variants that segregate according to Mendelian inheritance or by the sufficient polygenic accumulation of common variants, each with an individually small effect, or a combination of the two hypotheses. In 8,319 individuals across 1,589 migraine families, we calculated migraine polygenic risk scores (PRS) and found a significantly higher common variant burden in familial cases (n = 5,317, OR = 1.76, 95% CI = 1.71-1.81, p = 1.7 × 10-109) compared to population cases from the FINRISK cohort (n = 1,101, OR = 1.32, 95% CI = 1.25-1.38, p = 7.2 × 10-17). The PRS explained 1.6% of the phenotypic variance in the population cases and 3.5% in the familial cases (including 2.9% for migraine without aura, 5.5% for migraine with typical aura, and 8.2% for hemiplegic migraine). The results demonstrate a significant contribution of common polygenic variation to the familial aggregation of migraine.

The contribution of CACNA1A, ATP1A2 and SCN1A mutations in hemiplegic migraine: A clinical and genetic study in Finnish migraine families.

Objective To study the position of hemiplegic migraine in the clinical spectrum of migraine with aura and to reveal the importance of CACNA1A, ATP1A2 and SCN1A in the development of hemiplegic migraine in Finnish migraine families. Methods The International Classification of Headache Disorders 3rd edition criteria were used to determine clinical characteristics and occurrence of hemiplegic migraine, based on detailed questionnaires, in a Finnish migraine family collection consisting of 9087 subjects. Involvement of CACNA1A, ATP1A2 and SCN1A was studied using whole exome sequencing data from 293 patients with hemiplegic migraine. Results Overall, hemiplegic migraine patients reported clinically more severe headache and aura episodes than non-hemiplegic migraine with aura patients. We identified two mutations, c.1816G>A (p.Ala606Thr) and c.1148G>A (p.Arg383His), in ATP1A2 and one mutation, c.1994C>T (p.Thr665Met) in CACNA1A. Conclusions The results highlight hemiplegic migraine as a clinically and genetically heterogeneous disease. Hemiplegic migraine patients do not form a clearly separate group with distinct symptoms, but rather have an extreme phenotype in the migraine with aura continuum. We have shown that mutations in CACNA1A, ATP1A2 and SCN1A are not the major cause of the disease in Finnish hemiplegic migraine patients, suggesting that there are additional genetic factors contributing to the phenotype.

Phenotypic Consequences of a Genetic Predisposition to Enhanced Nitric Oxide Signaling.

BACKGROUND: Nitric oxide signaling plays a key role in the regulation of vascular tone and platelet activation. Here, we seek to understand the impact of a genetic predisposition to enhanced nitric oxide signaling on risk for cardiovascular diseases, thus informing the potential utility of pharmacological stimulation of the nitric oxide pathway as a therapeutic strategy. METHODS: We analyzed the association of common and rare genetic variants in 2 genes that mediate nitric oxide signaling (Nitric Oxide Synthase 3 [NOS3] and Guanylate Cyclase 1, Soluble, Alpha 3 [GUCY1A3]) with a range of human phenotypes. We selected 2 common variants (rs3918226 in NOS3 and rs7692387 in GUCY1A3) known to associate with increased NOS3 and GUCY1A3 expression and reduced mean arterial pressure, combined them into a genetic score, and standardized this exposure to a 5 mm Hg reduction in mean arterial pressure. Using individual-level data from 335 464 participants in the UK Biobank and summary association results from 7 large-scale genome-wide association studies, we examined the effect of this nitric oxide signaling score on cardiometabolic and other diseases. We also examined whether rare loss-of-function mutations in NOS3 and GUCY1A3 were associated with coronary heart disease using gene sequencing data from the Myocardial Infarction Genetics Consortium (n=27 815). RESULTS: A genetic predisposition to enhanced nitric oxide signaling was associated with reduced risks of coronary heart disease (odds ratio, 0.37; 95% confidence interval [CI], 0.31-0.45; P=5.5*10-26], peripheral arterial disease (odds ratio 0.42; 95% CI, 0.26-0.68; P=0.0005), and stroke (odds ratio, 0.53; 95% CI, 0.37-0.76; P=0.0006). In a mediation analysis, the effect of the genetic score on decreased coronary heart disease risk extended beyond its effect on blood pressure. Conversely, rare variants that inactivate the NOS3 or GUCY1A3 genes were associated with a 23 mm Hg higher systolic blood pressure (95% CI, 12-34; P=5.6*10-5) and a 3-fold higher risk of coronary heart disease (odds ratio, 3.03; 95% CI, 1.29-7.12; P=0.01). CONCLUSIONS: A genetic predisposition to enhanced nitric oxide signaling is associated with reduced risks of coronary heart disease, peripheral arterial disease, and stroke. Pharmacological stimulation of nitric oxide signaling may prove useful in the prevention or treatment of cardiovascular disease.

Shared genetic risk between migraine and coronary artery disease: A genome-wide analysis of common variants.

Migraine is a recurrent pain condition traditionally viewed as a neurovascular disorder, but little is known of its vascular basis. In epidemiological studies migraine is associated with an increased risk of cardiovascular disease, including coronary artery disease (CAD), suggesting shared pathogenic mechanisms. This study aimed to determine the genetic overlap between migraine and CAD, and to identify shared genetic risk loci, utilizing a conditional false discovery rate approach and data from two large-scale genome-wide association studies (GWAS) of CAD (C4D, 15,420 cases, 15,062 controls; CARDIoGRAM, 22,233 cases, 64,762 controls) and one of migraine (22,120 cases, 91,284 controls). We found significant enrichment of genetic variants associated with CAD as a function of their association with migraine, which was replicated across two independent CAD GWAS studies. One shared risk locus in the PHACTR1 gene (conjunctional false discovery rate for index SNP rs9349379 < 3.90 x 10-5), which was also identified in previous studies, explained much of the enrichment. Two further loci (in KCNK5 and AS3MT) showed evidence for shared risk (conjunctional false discovery rate < 0.05). The index SNPs at two of the three loci had opposite effect directions in migraine and CAD. Our results confirm previous reports that migraine and CAD share genetic risk loci in excess of what would be expected by chance, and highlight one shared risk locus in PHACTR1. Understanding the biological mechanisms underpinning this shared risk is likely to improve our understanding of both disorders.

Heterogeneous contribution of microdeletions in the development of common generalised and focal epilepsies.

BACKGROUND: Microdeletions are known to confer risk to epilepsy, particularly at genomic rearrangement 'hotspot' loci. However, microdeletion burden not overlapping these regions or within different epilepsy subtypes has not been ascertained. OBJECTIVE: To decipher the role of microdeletions outside hotspots loci and risk assessment by epilepsy subtype. METHODS: We assessed the burden, frequency and genomic content of rare, large microdeletions found in a previously published cohort of 1366 patients with genetic generalised epilepsy (GGE) in addition to two sets of additional unpublished genome-wide microdeletions found in 281 patients with rolandic epilepsy (RE) and 807 patients with adult focal epilepsy (AFE), totalling 2454 cases. Microdeletions were assessed in a combined and subtype-specific approaches against 6746 controls. RESULTS: When hotspots are considered, we detected an enrichment of microdeletions in the combined epilepsy analysis (adjusted p=1.06×10-6,OR 1.89, 95% CI 1.51 to 2.35). Epilepsy subtype-specific analyses showed that hotspot microdeletions in the GGE subgroup contribute most of the overall signal (adjusted p=9.79×10-12, OR 7.45, 95% CI 4.20-13.5). Outside hotspots , microdeletions were enriched in the GGE cohort for neurodevelopmental genes (adjusted p=9.13×10-3,OR 2.85, 95% CI 1.62-4.94). No additional signal was observed for RE and AFE. Still, gene-content analysis identified known (NRXN1, RBFOX1 and PCDH7) and novel (LOC102723362) candidate genes across epilepsy subtypes that were not deleted in controls. CONCLUSIONS: Our results show a heterogeneous effect of recurrent and non-recurrent microdeletions as part of the genetic architecture of GGE and a minor contribution in the aetiology of RE and AFE.

Pitfalls in genetic testing: the story of missed SCN1A mutations.

BACKGROUND: Sanger sequencing, still the standard technique for genetic testing in most diagnostic laboratories and until recently widely used in research, is gradually being complemented by next-generation sequencing (NGS). No single mutation detection technique is however perfect in identifying all mutations. Therefore, we wondered to what extent inconsistencies between Sanger sequencing and NGS affect the molecular diagnosis of patients. Since mutations in SCN1A, the major gene implicated in epilepsy, are found in the majority of Dravet syndrome (DS) patients, we focused on missed SCN1A mutations. METHODS: We sent out a survey to 16 genetic centers performing SCN1A testing. RESULTS: We collected data on 28 mutations initially missed using Sanger sequencing. All patients were falsely reported as SCN1A mutation-negative, both due to technical limitations and human errors. CONCLUSION: We illustrate the pitfalls of Sanger sequencing and most importantly provide evidence that SCN1A mutations are an even more frequent cause of DS than already anticipated.

Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine.

Migraine is a debilitating neurological disorder affecting around one in seven people worldwide, but its molecular mechanisms remain poorly understood. There is some debate about whether migraine is a disease of vascular dysfunction or a result of neuronal dysfunction with secondary vascular changes. Genome-wide association (GWA) studies have thus far identified 13 independent loci associated with migraine. To identify new susceptibility loci, we carried out a genetic study of migraine on 59,674 affected subjects and 316,078 controls from 22 GWA studies. We identified 44 independent single-nucleotide polymorphisms (SNPs) significantly associated with migraine risk (P < 5 × 10(-8)) that mapped to 38 distinct genomic loci, including 28 loci not previously reported and a locus that to our knowledge is the first to be identified on chromosome X. In subsequent computational analyses, the identified loci showed enrichment for genes expressed in vascular and smooth muscle tissues, consistent with a predominant theory of migraine that highlights vascular etiologies.

Evaluation of Presumably Disease Causing SCN1A Variants in a Cohort of Common Epilepsy Syndromes.

OBJECTIVE: The SCN1A gene, coding for the voltage-gated Na+ channel alpha subunit NaV1.1, is the clinically most relevant epilepsy gene. With the advent of high-throughput next-generation sequencing, clinical laboratories are generating an ever-increasing catalogue of SCN1A variants. Variants are more likely to be classified as pathogenic if they have already been identified previously in a patient with epilepsy. Here, we critically re-evaluate the pathogenicity of this class of variants in a cohort of patients with common epilepsy syndromes and subsequently ask whether a significant fraction of benign variants have been misclassified as pathogenic. METHODS: We screened a discovery cohort of 448 patients with a broad range of common genetic epilepsies and 734 controls for previously reported SCN1A mutations that were assumed to be disease causing. We re-evaluated the evidence for pathogenicity of the identified variants using in silico predictions, segregation, original reports, available functional data and assessment of allele frequencies in healthy individuals as well as in a follow up cohort of 777 patients. RESULTS AND INTERPRETATION: We identified 8 known missense mutations, previously reported as pathogenic, in a total of 17 unrelated epilepsy patients (17/448; 3.80%). Our re-evaluation indicates that 7 out of these 8 variants (p.R27T; p.R28C; p.R542Q; p.R604H; p.T1250M; p.E1308D; p.R1928G; NP_001159435.1) are not pathogenic. Only the p.T1174S mutation may be considered as a genetic risk factor for epilepsy of small effect size based on the enrichment in patients (P = 6.60 x 10-4; OR = 0.32, fishers exact test), previous functional studies but incomplete penetrance. Thus, incorporation of previous studies in genetic counseling of SCN1A sequencing results is challenging and may produce incorrect conclusions.

Simultaneous impairment of neuronal and metabolic function of mutated gephyrin in a patient with epileptic encephalopathy.

Synaptic inhibition is essential for shaping the dynamics of neuronal networks, and aberrant inhibition plays an important role in neurological disorders. Gephyrin is a central player at inhibitory postsynapses, directly binds and organizes GABAA and glycine receptors (GABAARs and GlyRs), and is thereby indispensable for normal inhibitory neurotransmission. Additionally, gephyrin catalyzes the synthesis of the molybdenum cofactor (MoCo) in peripheral tissue. We identified a de novo missense mutation (G375D) in the gephyrin gene (GPHN) in a patient with epileptic encephalopathy resembling Dravet syndrome. Although stably expressed and correctly folded, gephyrin-G375D was non-synaptically localized in neurons and acted dominant-negatively on the clustering of wild-type gephyrin leading to a marked decrease in GABAAR surface expression and GABAergic signaling. We identified a decreased binding affinity between gephyrin-G375D and the receptors, suggesting that Gly375 is essential for gephyrin-receptor complex formation. Surprisingly, gephyrin-G375D was also unable to synthesize MoCo and activate MoCo-dependent enzymes. Thus, we describe a missense mutation that affects both functions of gephyrin and suggest that the identified defect at GABAergic synapses is the mechanism underlying the patient's severe phenotype.