Abstract
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.
Glossary
- ρHESS=
- Heritability Estimation from Summary Statistics;
- AS=
- all stroke;
- CE=
- cardioembolic stroke;
- CeAD=
- cervical artery dissection;
- CI=
- confidence interval;
- GNOVA=
- Genetic Covariation Analyzer;
- GWAS=
- genome-wide association study;
- GWAS-PW=
- GWAS-pairwise;
- IS=
- ischemic stroke;
- IVW=
- inverse variance–weighted;
- LAS=
- large artery stroke;
- LD=
- linkage disequilibrium;
- LDSc=
- linkage disequilibrium score regression;
- MA=
- migraine with aura;
- MR=
- Mendelian randomization;
- MO=
- migraine without aura;
- MTAG=
- multitrait analysis of GWAS;
- OR=
- odds ratio;
- PPA3=
- posterior probability of association in mode 3;
- SNV=
- single nucleotide variation;
- SVS=
- small vessel stroke;
- WGHS=
- Women's Genome Health Study
Migraine, stroke, and extracranial cervical artery dissection (CeAD) represent a triad of brain disorders with vascular involvement and pairwise comorbid relationships that are pertinent to risk assessment and clinical care.1 While the shared clinical features of all 3 disorders point to vasculature as the basis of the comorbidity, precise underlying mechanisms are not established. Understanding the shared and distinct biological mechanisms has the potential to clarify the basis of shared risk while also informing potential prophylactic and treatment strategies.
From this perspective, previous investigations have leveraged the unique properties of human genetics to reveal shared biology among the 3 disorders while limiting the potential influence of reverse causality and confounding that may arise in conventional observational epidemiology. One study found that genome-wide genetic overlap with migraine was most significant for large artery stroke and significant for cardioembolic stroke (CE), contrary to observational associations that had linked migraine to small vessel disease.2,3 Associations were stronger for migraine without aura (MO) than those for either overall migraine or migraine with aura (MA), though the latter is a stronger risk factor of ischemic stroke (IS).4,-,6 Similarly, genetic associations at specific loci diverged from conventional observational associations. At the 9p21 locus, associations with stroke and MO had concordant direction, but there was no association with MA, and there remained uncertainty about whether the causal variants for MO and stroke at the locus were the same.3 At the FHL5 locus on chromosome 6, the associations with stroke and migraine were in opposite directions, while the same locus has been noted for a concordant association between migraine and CeAD, the latter from a genome-wide association study (GWAS) at subgenome-wide significance.7 The GWAS of CeAD also noted concordant effects with migraine at loci implicating the PHACTR1/EDN1 and LRP1 genes,8 but only the former was replicated in an independent follow-up sample.
Recent GWASs of stroke and migraine incorporating much larger samples than previously available (therefore with much greater power), as well as novel genetic methods and the lack of a systematic comparison among all 3 cerebrovascular disorders, invite a new genetic analysis toward resolving several outstanding questions. First, what is the extent of shared genetics among the 3 disorders? Second, which specific susceptibility loci are shared on a genome-wide basis? Finally, does human genetics support causal relationships underlying the increased risk of stroke and CeAD among individuals susceptible to migraine?
Methods
Overview
Pairwise genetic relationships among migraine, CeAD, and subtypes of stroke were examined using 4 analytic strategies. Genome-wide genetic correlations were calculated to assess for overall genetic sharing. Local genetic correlations were calculated to assess shared genetics within disjoint segments across the genome and at specific candidate loci previously identified for association with at least one of the disorders. A genome-wide cross-trait association analysis was implemented to identify novel variant associations for each trait. Finally, the Mendelian randomization (MR) analysis was performed to assess potential causal influences of migraine on the other cerebrovascular disorders.
Summary Statistics
Analyses used discovery summary statistics from published, consortium-based GWASs of migraine,9 CeAD,7 and stroke.10 The total numbers of samples included in these summary statistics were as follows: any migraine (59,674 cases/316,078 controls), MA (6,332 cases/144,883 controls), MO (8,348 cases/139,622 controls), CeAD (carotoid and vertebral, 1,393 cases/14,416 controls), all stroke (AS) (40,584 cases/406,111 controls), IS (34,217 cases/406,111 controls), large artery stroke (LAS) (4,373 cases/297,290 controls), CE (7,193 cases/355,4468 controls), and small vessel stroke (SVS) (5,386 cases/343,560 controls). All summary statistics were derived from study populations exclusively with European ancestry. The migraine and stroke GWASs were based on 1000 Genomes Project imputed data (hg19) and included approximately 8 million single nucleotide variations ([SNVs], formerly SNPs), while the CeAD GWAS was based on HapMap (hg18) and 1000 Genomes Project (August 2010 release) imputed data and included approximately 6.6 million SNVs. The Women's Genome Health Study (WGHS)11 provided a substantial proportion of cases and controls to the GWAS for MA (1,177 WGHS cases/6,332 total cases) and MO (1,826 WGHS cases/8,348 total cases) while also contributing a smaller proportion of cases to overall migraine (N = 5,122 WGHS cases/59,674 total cases) and stroke (N = 422 cases of all stroke/22,795 controls) GWASs. Other cohorts may have contributed smaller numbers of samples to the GWAS for migraine and stroke. Potential bias due to this overlap was addressed in 2 ways. First, migraine summary statistics from a meta-analysis were as described but omitted the WGHS contribution.9 Second, the genetic correlation and cross-trait association methods intrinsically account for any minimal residual overlap.12,-,14
Summary statistics for the 23andMe cohort were obtained under an agreement with 23andMe that protects the privacy of the 23andMe participants. The participants of 23andMe provided informed consent and participated in the research under a protocol approved by the external AAHRPP-accredited IRB, Ethical & Independent Review Services. The full GWAS summary statistics for the 23andMe discovery data set will be made available through 23andMe to qualified researchers under an agreement with 23andMe that protects the privacy of the 23andMe participants. Information about access the data from 23andMe can be found at research.23andme.com/collaborate/#dataset-access/. Use of other summary statistics wass consistent with the local IRBs of each of the contributing cohorts or samples. All genomic coordinates refer to genome build hg19.
Genome-wide Genetic Correlation
Two established methods were used to estimate a genome-wide genetic correlation with the GWAS summary statistics: the conventional approach, linkage disequilibrium (LD) score regression (LDSc, version 1.0.0), and a similar but potentially more powerful approach, Genetic Covariation Analyzer (GNOVA) (downloaded in December 2017).12,13 An analysis with LDSc incorporated precomputed LD measures for approximately 1.3 million common SNVs based on the HapMap with minimum minor allele frequency of approximately 10%12 and shared across all of the summary statistics. An analysis with GNOVA included a step to calculate LD relationships among individuals with European ancestry in the 1000 Genomes reference panel including approximately 6.1 million SNVs with minimum minor allele frequency 5%. Genome-wide genetic correlation is the principal estimate in LDSc. By contrast, genetic covariance is the principal estimate in GNOVA, and genetic correlation is derived by scaling with the single-trait heritability estimates. As such, p values refer to the genetic correlation in LDSc and to the genetic covariance in GNOVA. Both LDSc and GNOVA have options to compute genetic correlation or covariance while adjusting for potential sample overlap or other potential causes of inflation. These options were invoked in all analyses. Although LDSc and GNOVA are substantially similar methods and the use of both provides cross-validation, differences in the minor allele frequency thresholds are expected to influence genetic correlation estimates and significance to some extent.
Locally Shared Genetic Effects
Local (as opposed to genome-wide) genetic correlation was estimated by 2 approaches, applied initially to approximately 1,704 prespecified disjoint segments across the genome with minimal intersegment LD within the GWAS summary statistics.15 One approach, ρ Heritability Estimation from Summary Statistics (ρHESS) (version 0.5),16 provided a frequentist estimate of the local covariance of genetic effects, while the other approach, GWAS-pairwise (GWAS-PW) (version 0.21),14 provided a posterior probability of a locally shared genetic association (posterior probability of association in mode 3 [PPA3]) within each segment in an empirical Bayes framework based on GWAS p values. Genome-wide significance across the prespecified segments was 2.9 × 10−5 (=0.05/1,704) in ρHESS and PPA3 >0.9 in GWAS-PW. Candidate genes were assigned for regions showing PPA3 >0.9 based on proximity to potentially shared causal variants. ρHESS was also adapted to examine candidate segments surrounding genome-wide significant loci for stroke or migraine, as previously performed.16 In this study, candidate segments were defined to include the SNVs neighboring each of GWAS index SNV, such that all SNVs outside of the segment had LD r2 < 0.1 to the index SNV. Testing 69 candidate loci for each pairwise comparison, p < 7 × 10−4 (=0.05/69) was considered significant. As needed, pairwise SNV LD was estimated using a European ancestry panel from the 1000 Genomes Project with PLINK or LDlink.17,-,19
A Cross-Trait Association Analysis With a Multitrait Analysis of GWAS
A cross-trait association analysis was performed using a multitrait analysis of GWAS (MTAG, version 1.0.7), which leverages the pairwise genome-wide trait genetic correlation to boost power in association testing.20 In MTAG, LDSc provided the estimates of pairwise genetic correlation. The significance threshold in the MTAG required pMTAG < 1.67 × 10−8 (=5.00 × 10−8/3 phenotypes) but was also restricted to SNVs that also had nominal significance (p < 0.05) for each phenotype separately in the preexisting univariate GWAS.
MR Instrumental Analysis
Genetic instrumental variable analysis with MR was performed using the 2-sample method, prioritizing the random-effects inverse variance–weighted (IVW) estimator within the package TwoSampleMR in the R computing environment.21,22 Sensitivity analysis included MR-weighted median, MR-Egger, MR-Egger bootstrap, MR-robust adjusted profile score, and MR-PRESSO.23,-,26 The latter detects and excludes instruments that are consistent with horizontal pleiotropy, a violation of assumptions underlying MR, and then evaluates an overall estimate with the remaining instruments using the IVW method. Pleiotropy was also assessed with Cochran Q for heterogeneity and the MR-Egger intercept as recommended.25 MR was limited to effects of migraine on the other disorders because (1) only the GWAS of any migraine, i.e., not the other disorders, had a sufficient number of genome-wide significant SNVs (N ≥ 10) for use as exposure instruments and (2) because migraine typically precedes stroke or CeAD, potential causal effects of migraine on the other disorders are most consistent with temporal plausibility. Because migraine is a binary exposure, MR effect estimates were scaled by 0.693 for reporting to represent the odds of the outcome for a doubling of the odds of the exposure.27
Standard Protocol Approvals, Registrations, and Patient Consents
The GWAS summary statistics were all derived by a meta-analysis. Participants who contributed to cohort-level summary statistics constituting the meta-analyses provided written informed consent, and each of the cohort protocols was approved by a local institutional review board.
Data Availability
This study used only GWAS summary statistics from published reports as described earlier. The availability of these data or procedures for accessing them is documented in the cited publications. Summary statistics for the GWAS of migraine lacking contribution from the WGHS will be made available by sending an application to the International Headache Genetics Consortium via the corresponding author using the same procedure that governs access to the summary statistics for the published migraine study.9
Results
Genome-wide Genetic Correlation
Genetic correlations involving CeAD and stroke or migraine were generally positive (i.e., concordant) and comparable between LDSc and GNOVA, but significant only with GNOVA for the combinations of migraine (rg [ρ (ρ SE)] = 0.22 [0.048 (0.012)], p = 4.9e-05) or MO (0.29 [0.051 (0.016)], 0.0017) with CeAD, after accounting for multiple testing (Table 1). The estimated genetic correlations were larger but only nominally significant with LDSc (rg = 0.45 and 0.41, respectively).
Genome-wide Genetic Correlations Between Pairs of Brain Disorders
Local Genetic Correlation
At the experiment-wide significance threshold (PPA3 >0.9), GWAS-PW implicated novel locally concordant associations of migraine and CE on chromosome 9q22.32 (top SNV rs113154802 near PTCH1) (Table 2). The remaining significant segments all include loci previously recognized by GWAS for 1 or more of the disorders, although many are newly implicated for an additional disorder (indicated in bold). These loci were ADAMTSL4/ECM1: CeAD-migraine, CARF: SVS-migraine, NKX2: migraine-CE, HDAC9: migraine-AS/IS, ARMS2: AS/IS/SVS-migraine, LRCH1: migraine-AS/IS, and COL4A1: migraine-AS. At COL4A1, neither migraine nor AS is genome-wide significant in the summary statistics used in this study, which are derived from population of European ancestry alone, but the locus has been recently identified for stroke by a trans-ancestry meta-analysis with index SNV rs9521634.28 However, this variant is not in LD (r2 = 0.02, D′ = 0.52) with the top SNV for the joint analysis of stroke and migraine, rs650724, which is also the top SNV for stroke alone in the current summary statistics derived from European ancestry and in high LD (r2 = 0.86, D′ = 0.94) with the top SNV for migraine alone (rs2000660). Identification of shared, concordant associations involving CeAD in segments encompassing LRP1 and FHL5 is consistent with strong, subgenome-wide significant associations in these regions previously noted.7 Although already known for any migraine, these 2 loci are newly implicated in both MA and MO. None of the segments was significant for combinations of CeAD and stroke. In contrast to the GWAS-PW method, local genetic covariance assessed with ρHESS did not meet experiment-wide significance for any pairwise combination of the cerebrovascular disorders (all p for local rg > 2.9 × 10−5 [=0.05/1704]).
Local Prespecified Segments With Significant Joint GWAS-PW Association Among Cerebrovascular Conditions
However, with ρHESS, local genetic covariance was also assessed at candidate regions defined by LD r2 >0.1 around the 69 known genome-wide significant loci for each of the disorders (Methods, Table 3).16 Local genetic covariance was concordant and met significance thresholds for candidate analysis (p < 0.0007 [=0.05/69]) between migraine and CeAD at PHACTR1/EDN1 and LRP1 as previously suggested but not formally demonstrated.7 At nominal significance (p < 0.05), concordant covariance was observed similarly not only at these candidate loci between CeAD and MO but also at FHL5 between CeAD and both any migraine and MO. The nominal associations also suggested shared signals between migraine and various stroke subtypes at known migraine loci mapping to PRDM16 (SVS, opposite directionality), ARMS2/HTRA1 (AS and IS, opposite), LRP1 (LAS, opposite), LRCH1 (AS and IS, opposite), and RNF213 (IS, opposite); and between migraine and CeAD at PLCE1 (concordant) and FGF6 (opposite). None of the nominally significant local correlations implicated MA.
Nominally Significant Local Genetic Covariance (ρg) at 69 Candidate Loci From Previous GWAS
Novel Genome-wide Significant SNVs in Cross-Trait Association Analysis
MTAG, which leverages pairwise genome-wide genetic correlations to boost univariate association signals, identified novel genome-wide signals among SNVs that were also nominally significant in single-trait analysis (Methods, Table 4). Combining migraine with CeAD, there were novel associations for CeAD at SNVs mapping to the PLCE1 (chr. 10, rs57866767) and MRVI1 (chr 11, rs7940646) genes. The former was also nominally significant in the candidate local genetic correlation analysis (previous section). MTAG associations also recapitulated the findings from GWAS-PW for CeAD at FHL5 (chr. 6, rs2971603 or rs9486725), LRP1 (chr. 12, rs11172113), and at ADAMTSL4/ECM1 (chr. 1, rs6693567) for SNVs that were previously genome-wide significant for migraine. Two SNVs in moderately high LD (R2 = 0.5, D′ = 1.0), rs2971603 and rs9486725, represent the top associations at FHL5, the former more significant with any migraine and the latter with CeAD. Combining migraine and stroke, MTAG identified 2 novel stroke loci: LRP1 (chr. 12, rs11172113) for LAS at a SNV shared with CeAD (and migraine, opposite effect), and CARF (chr. 2, rs191602009) for SVS also previously genome-wide significant for migraine (concordant effect). No locus was genome-wide significant for CeAD in combination with any of the stroke outcomes.
Pairwise MTAG Genome-wide Significant Associations That Are Also Nominally Significant in Original GWAS
Mendelian Randomization
In MR analysis (i.e., genetic instrumental analysis), liability to migraine was supported as causal for increased CeAD risk (odds ratio [OR] [95% confidence interval (CI)] = 1.69 [1.24–2.3], p = 0.0009) but protective for LAS (0.86 [0.76–0.96], p = 0.007) (Figure 1). There were no significant effects on either AS or other stroke subtypes, including all IS. There was significant heterogeneity detected for the migraine-CeAD effect (Cochran Q = 109, df = 39, p = 1.48 × 10−8), which was diminished but not eliminated by exclusion of 2 clearly pleiotropic SNVs, rs11172113 (LRP1) and rs9349379 (PHACTR1/EDN1) (Q = 65, df = 37, p = 0.003), leading to an attenuated but still nominally significant effect (OR [95% CI] = 1.33 [1.02–1.73], p = 0.04) (eTables 1 and 2, links.lww.com/NXG/A511). The MR-Egger intercept test for directional pleiotropy was null (p = 0.85), as suggested also by the largely consistent estimates of the effect in the sensitivity analyses (eTable 3, links.lww.com/NXG/A511). By contrast, there was no significant heterogeneity in the effect of migraine on LAS (Cochran Q = 46, df = 39, p = 0.21), and sensitivity analyses for pleiotropy yielded consistent and significant estimated protective effects of migraine, with the exception of MR-Egger (1.01 [0.73–1.39], p = 0.97). However, the effect obtained from the MR-Egger bootstrap test (0.85 [0.66–1.08], p = 0.091) was directionally consistent with the primary analysis, suggesting a potential undue influence of outliers on the estimate from MR-Egger when including all instruments without bootstrapping.
MR was performed using the random-effects inverse variance–weighted method. ORs and 95% CIs are scaled to reflect the effect of migraine liability on CeAD and stroke per doubling of migraine prevalence (see Methods). CeAD = cervical artery dissection; CI = confidence interval; MR = Mendelian randomization; OR = odds ratio.
Discussion
The preceding analysis was undertaken to investigate the etiologic basis of comorbidity among each pair of 3 brain disorders with known vascular involvement through the unique properties of genetics. Both genome-wide and at specific loci, the findings emphasized extensive sharing of biology between migraine and CeAD. Genetic sharing was less for migraine and stroke but still implicated a few loci, while still less sharing was detected for stroke and CeAD. Figure 2 summarizes the significant pairwise associations from all analyses.
Only pairs with significant associations are shown. For Mendelian randomization, the direction of each arrow indicates whether the direction of inferred causality is from the first disorder of each pair to the second (rightward arrow) or the opposite (leftward arrow) with concordant (blue) or opposite (red) functional relationship. Symbols following gene name labels for each locus refer to the disorder(s) for which association has been previously reported, if any, as follows: m = migraine, c = CeAD, s = stroke, and - = no known association. CeAD = cervical artery dissection; GWAS-PW = genome-wide association study-pairwise; MTAG = multitrait analysis of GWAS; ρHESS = ρ Heritability Estimation from Summary Statistics.
While the challenge of recruiting large samples of CeAD cases for genome-wide genetic analysis had limited power for previous genome-wide analysis, local genetic correlation with migraine boosted genetic signals to highlight novel candidate loci for CeAD and reinforced existing candidates, all with concordant effects on migraine. Novel genomic regions on at chr.1q21.3 (ADAMTSL4/ECM1 candidate genes), chr10q23.33 (PLCE1), and chr11p15.4 (MRVI1) were all previously associated with migraine9 and are now implicated also with CeAD. All reached genome-wide significance in the MTAG, and the chromosome 1 and 10 loci were further supported by GWAS-PW and ρHESS analyses, respectively. The extracellular matrix protein 1 (ECM1) gene at the chr1q21.3 locus has been suggested for involvement in vascular development.29 The association at the MRVI1 gene, encoding murine retrovirus integration site 1 homolog, a tumor suppressor, arose previously as a candidate influencing both migraine and IS.3 The lead SNV, rs7940646, and its LD proxies associated with migraine in the region are also associated with blood pressure, arterial stiffness, airway inflammatory diseases, platelet aggregation, brain region volume, and measures for white matter integrity.30,-,34 PLCE1 encodes the phospholipase C, epsilon 1 protein, at which GWAS has implicated other vascular conditions including CVD and blood pressure.33,35 The remaining CeAD loci shared with migraine at the PHACTR1/EDN1, LRP1, and FHL5 genes were all noted in the original CeAD GWAS (FHL5 at suggestive significance), but only PHACTR1/EDN1 was formally genome-wide significant in replication, again likely due to limited sample.7 All have been previously annotated as playing a role in vascular development or function.29 At the PHACTR1/EDN1 locus, recent functional work has not fully resolved which of 2 candidate genes may underlie the functional effects.8,36 Though there is no strong genetic overlap between stroke and CeAD, the genetic overlap involving LRP1 implicates both with migraine. Owing to the vascular etiology of CeAD and stroke, the preceding shared loci may be particularly relevant to vascular etiologies of migraine.
Despite the modest genome-wide correlations, local comparisons with GWAS-PW revealed potential new candidate loci with concordant effects for various stroke subtypes and any migraine. CARF (rs191602009), a known migraine locus encoding the calcium-response transcription factor, likely mediates calcium signaling in neurons, including regulation of the brain-derived neurotrophic factor,37 and was implicated in SVS; NKX2-5 (rs55884259), a known CE stroke locus encoding a homeobox protein for which variations cause congenital heart defects, was implicated in migraine; and PTCH1 (rs182923402), encoding patched 1, is a member of a protein family of receptors that are ligands for sonic hedgehog signaling peptides in development and was implicated in both migraine and CE. Other pairwise shared loci revealed novel associations for 1 disorder but with opposite risk alleles compared with the known association at the second disorder: (1) HDAC9/TWIST1 (rs2107595), a known LAS locus for which the candidate gene is not yet definitively identified,38,39 is implicated in migraine; (2) ARMS2/HTRA1 (3 SNVs [all D′ = 1, low R2 with each other and with rs10490924]), a known a migraine locus for which HTRA1 encoding a serine peptidase may be the best candidate gene, is implicated in AS, IS, and SVS, reinforcing previous subthreshold associations.10 Mendelian associations with small vessel disease also support this association because monogenic variations in HTRA1 lead to a rare autosomal dominant form of SVD, CARASIL40,41; (3) LRCH1 (2 SNVs, LD R2 = 1), a migraine locus, is now implicated in AS and IS; and (4) COL4A1 (rs650724), encoding collagen type 4 alpha 1, is now implicated in migraine having been previously identified for stroke by GWAS and Mendelian genetics of SVD.10,42 Signals at ARMS2/HTRA1 and LRCH1 were also supported at nominal significance by the ρHESS candidate local genetic correlations, while the signal at rs191602009 (CARF) was supported in the MTAG.
The less extensive genetic sharing of migraine (particularly MA) and IS across the genome is contrary to their strong comorbidity in epidemiologic studies.4 Genome-wide genetic correlations were not only modest but also emphasized an opposite relationship rather than concordance, particularly between migraine and LAS. Similarly, causality in an opposite relationship of migraine liability to LAS was supported by the instrumental analysis. This finding is consistent with prior MR analyses, which identified opposite instrumental relationships of migraine liability with coronary artery disease, a disorder that shares pathophysiology with LAS.43 Similarly, the findings are reminiscent of previous analysis of shared genetics between migraine, especially MO, and stroke that also found overlap, especially for large artery and CE types.3 Although the genetic correlation of migraine with SVS, which has been suspected in the mechanism of migraine comorbidity,2 was concordant with findings from conventional epidemiology and with migraine being an important feature of monogenic forms of small vessel disease, the estimate was only marginally significant. This observation may be qualified, however, by the low power of GWAS for IS subtypes, as well as likely imperfect ascertainment of SVS in many studies. The genetic relationship between migraine and small vessel disease deserves further investigation using more specific MRI markers of SVD, such as white matter hyperintensity burden.44
An MTAG-based genome-wide significant association at rs11172113 for LAS that was supported by local correlation at nominal significance in ρHESS with migraine implicates LRP1, which is the only locus influencing risk of all 3 cerebrovascular disorders, although opposite in its effect on stroke compared with that on migraine or CeAD. This same locus has recently also been implicated in aortic and coronary dissection and abdominal aortic aneurysm with the same directionality for CeAD and migraine, placing LRP1 at the center of shared biology and deserving further study.45,-,47 LRP1, a member of the LDL receptor family, has been implicated by GWAS also in pulmonary function and CHD, the latter likely related in pathophysiology to the association with LAS.33,48,49 LRP1 protein is involved in endocytosis of a wide variety of ligands, including lipoproteins, and understanding mechanism(s) of its contribution to the shared susceptibilities will require further research.50
The strengths of this study are the very large sample sizes and therefore power represented by the GWAS summary statistics for migraine and stroke. The study is limited in its restriction to populations of European ancestry, although multiancestry meta-analysis for stroke subtypes supports the top loci, implying that relevant biological functions are shared among European and other ancestries. However, it remains possible that genetic relationships in non-European ancestries among the 3 disorders would highlight additional relationships, including those that may contribute to health disparities. The study is also limited by the modest sample size underlying the summary statistics for CeAD, a consequence of the challenge in accumulating genome-wide genetic data for extremely low prevalence events. The incidence of CeAD is only on the order of approximately 2.6 per 100,000 per year.51 Similarly, despite the relatively large total sample for the stroke GWAS, heterogeneity in stroke mechanism and intrinsic difficulties in assigning stroke subclassifications may have limited the ability to detect genetic overlap with either migraine or CeAD. An additional consequence of the limitations in the CeAD and stroke GWASs was an insufficient number of qualifying instruments to perform MR for assessing potential causal effects of liability to these disorders on migraine.
Taken together, the results thus provide novel support for the contribution of vascular functions to migraine and enhance understanding of the comorbidity among migraine, CeAD, and stroke. Future functional studies prioritizing specific loci identified through this genetic analysis may reveal deeper insights into corresponding vascular mechanisms leading to susceptibility to the 3 brain disorders.
Study Funding
This work was supported by grants from the NINDS of the NIH (R21NS092963 and R21NS104398 to D.I.C.). The CADISP study has been supported by Inserm, Lille 2 University, Institut Pasteur de Lille, and Lille University Hospital and received funding from the ERDF (FEDER funds) and Région Nord-Pas de Calais in the frame of Contrat de Projets Etat-Region 2007–2013 Région Nord-Pas-de-Calais—Grant No. 09120030, Centre National de Genotypage, Emil Aaltonen Foundation, Paavo Ilmari Ahvenainen Foundation, Helsinki University Central Hospital Research Fund, Helsinki University Medical Foundation, Päivikki and Sakari Sohlberg Foundation, Aarne Koskelo Foundation, Maire Taponen Foundation, Aarne and Aili Turunen Foundation, Lilly Foundation, Alfred Kordelin Foundation, Finnish Medical Foundation, Orion Farmos Research Foundation, Maud Kuistila Foundation, the Finnish Brain Foundation, Biomedicum Helsinki Foundation, Projet Hospitalier de Recherche Clinique Régional, Fondation de France, Génopôle de Lille, Adrinord, Basel Stroke-Funds, Käthe-Zingg-Schwichtenberg-Fonds of the Swiss Academy of Medical Sciences, and Swiss Heart Foundation. Murali Sargurupremraj has been supported by a grant from the European Commission under H2020 Grant no. 640643.
Disclosure
P. Gormley is a current employee and/or stockholder of Glaxo SmithKline, but the work was conducted while employed by Massachusetts General Hospital, Boston, MA. T. Kurth reports to have received personal compensation for outside the submitted work from Eli Lilly & Company, Teva, Total S.E., and the BMJ. The other authors report no disclosures. Go to Neurology.org/NG for full disclosures.
Acknowledgment
The authors thank the IHGC researchers and member cohort participants, the staff, and participants of all CADISP centers and the research participants and employees of 23andMe for their important contributions that made this work possible. See supplementary material for listings of IHGC and CADISP researchers (links.lww.com/NXG/A511).
Appendix Authors
Footnotes
Go to Neurology.org/NG for full disclosures. Funding information is provided at the end of the article.
↵* These authors contributed equally to this work.
↵† These authors equally supervised this work.
The Article Processing Charge was funded by the authors.
- Received October 28, 2021.
- Accepted in final form December 7, 2021.
- Copyright © 2022 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Academy of Neurology.
This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND), which permits downloading and sharing the work provided it is properly cited. The work cannot be changed in any way or used commercially without permission from the journal.
References
Letters: Rapid online correspondence
REQUIREMENTS
You must ensure that your Disclosures have been updated within the previous six months. Please go to our Submission Site to add or update your Disclosure information.
Your co-authors must send a completed Publishing Agreement Form to Neurology Staff (not necessary for the lead/corresponding author as the form below will suffice) before you upload your comment.
If you are responding to a comment that was written about an article you originally authored:
You (and co-authors) do not need to fill out forms or check disclosures as author forms are still valid
and apply to letter.
Submission specifications:
- Submissions must be < 200 words with < 5 references. Reference 1 must be the article on which you are commenting.
- Submissions should not have more than 5 authors. (Exception: original author replies can include all original authors of the article)
- Submit only on articles published within 6 months of issue date.
- Do not be redundant. Read any comments already posted on the article prior to submission.
- Submitted comments are subject to editing and editor review prior to posting.
