Causal relationships between ischemic stroke and epilepsy: a Mendelian randomization study

DOI: https://doi.org/10.21203/rs.3.rs-2570721/v1

Abstract

Background

Epidemiological evidence suggests that ischemic stroke may increase the risk of epilepsy. However, the causal relationship between these two conditions is currently inconclusive. We aimed to elucidate the possible causal relationship between ischemic stroke and related subtypes of epilepsy by a two-sample Mendelian randomization (MR) analysis.

Methods

We extracted 11 single nucleotide polymorphisms (SNPs) significantly associated with ischemic stroke from a large genomic association study (GWAS) of European individuals (440328 individuals) as instrumental variables and used the GWAS of epilepsy-related subtypes provided by the International League Against Epilepsy as an outcome data source. Inverse variance weighting (IVW) was used as the primary method to assess effect estimates. Sensitivity analysis was performed using weighted median, MR-Egger and MR-PRESSO to assess heterogeneity and pleiotropy.

Results

In the IVW random-effects model, 11 SNPs were used as instrumental variables, and a significant causal relationship was found between genetic susceptibility to ischemic stroke and generalized epilepsy (OR, 1.24; 95% confidence interval (CI), 1.09, 1.42; P = 0.0009). The results of the weighted median and MR-Egger methods were in general agreement. Sensitivity analysis showed no heterogeneity or polymorphism. Reverse MR analysis failed to detect a causal relationship between generalized epilepsy and ischemic stroke. In addition, we found no causal relationship between ischemic stroke and focal epilepsy (OR, 1.06; 95% confidence interval (CI), 0.93, 1.20; P = 0.3957).

Conclusion

There is a causal relationship between ischemic stroke and the occurrence of generalized epilepsy, but not with focal epilepsy.

Introduction

Epilepsy, one of the most common neurological disorders, affects approximately 70 million people worldwide, and at a rate of 50.4 to 81.7 cases per 100,000 people per year[1; 2]. Idiopathic epilepsy is thought to have no obvious triggers, but triggered epilepsy seem to be associated with many conditions, such as cranial trauma[3], intracranial infections[4], intracranial tumors, and cerebrovascular diseases[5]. Although the pathogenesis of epilepsy is unclear, some studies suggest a preference for syndromes with multiple risk factors and strong genetic susceptibility to induce epilepsy[6], such as stroke.

Stroke is a common cardiovascular disease and is the second leading cause of death in the world [7], has been reported with approximately 10% of epilepsy and 55% of newly diagnosed epilepsy in older adults being caused by stroke[8]. What’s more, the distribution of seizures is bimodal: peak seizures are always seen in childhood and after age 60[6], which suggests that damage to the central system due to stroke may be a cause of epilepsy that cannot be ignored. However, the data on the occurrence of stroke-induced epilepsy vary considerably at this stage[9; 10] and are related to many confounding factors such as timing of diagnosis and treatment, rehabilitation care and economic conditions[11; 12; 13; 14]. Strokes are classified as ischemic strokes and hemorrhagic strokes. In a cohort study, children with ischemic stroke were found to have a 27.8-fold increased risk of future epilepsy compared to control children[15].Therefore, establishing a causal relationship between epilepsy and ischemic stroke can help improve the understanding of the pathogenesis of the disease and provide new evidence for clinical management and prognostic assessment.

Mendelian randomization (MR) as a statistical method is to use genetic variation associated with exposure factors as instrumental variables to assess the relationship between exposure factors and outcomes which relies on the natural, random classification of genetic variation during meiosis, resulting in a random distribution of genetic variation in the population.[16; 17; 18] Thus, MR is similar to randomized controlled studies in that it avoids the effects of confounding factors and reverse causality. In the present study, we performed a comprehensive two-sample MR analysis with the aim of exploring the causal relationship between ischemic stroke and epilepsy.

Materials And Methods

Data sources

For ischemic stroke, we obtained summary statistics on stroke from the publicly available database provided by the MEGASTROKE consortium, which includes 34217 cases and 406111 controls of European ancestry. The MEGASTROKE consortium defined stroke as a rapidly developing symptom of a local (or global) disturbance of brain function with no apparent etiology other than a vascular source that lasts more than 24 hours or results in death. The study received prior approval from the relevant institutional review board and informed consent was obtained from each study participant[19].

In addition, we obtained genome-wide association study (GWAS) statistics from the International League Against Epilepsy (ILAE) that included epilepsy-related subtypes. The focal epilepsy dataset contained 9,671 cases and 29,677 controls, and the generalized epilepsy dataset contained 3,769 cases and 29,677 controls, with specific diagnostic criteria, demographic characteristics. More details about epilepsy cases available in the original study[20].

All GWAS summary statistics used in the MR analysis are shown in Table 1. All data for this study were obtained from publicly available summary statistics from GWAS and ethical approval can be found in the original study.

Data processing

First, 18 single nucleotide polymorphisms (SNPs) with significant (P < 5×10− 8) association with ischemic stroke were screened from the GWAS dataset of ischemic stroke and then the parameter r2 threshold was set to 0.1 and kilobase pair (kb) to 10000 to exclude linkage disequilibrium. The remaining 17 SNPs were then matched with SNPs in the epilepsy database, and finally 11 SNPs significantly associated with stroke were obtained as instrumental variables for subsequent statistical processing. SNPs with a minor allele frequency of < 0.01 were also excluded to avoid potential the statistical bias from the original GWAS since they usually carry with low confidence. The R2 and F values were also calculated to determine the degree of explanation of exposure by instrumental variables and the presence of weak instrumental variables, respectively. R2 was calculated as R2 = 2 × EAF × (1-EAF) × Beta2, and F was calculated as F = R2 × (N-2)/(1- R2)[21].

Statistical analysis

This MR study is based on three hypotheses, as described in Fig. 1.

We performed a series of Mendelian randomization analyses to investigate the causal relationship between ischemic stroke and epilepsy.

First, we regressed exposure (genetic variance in ischemic stroke) and outcome (genetic variance in epilepsy) and analyzed the causal relationship between ischemic stroke and epilepsy using the inverse variance weighted random effects (IVW) method, which is mainly used for causal inference of genes and will provide the most precise results, and therefore we used it as the primary method for MR [22]. In addition, we supplemented it with MR Egger and weighted median methods to identify any bias in assessing the MR hypothesis. The major difference between the MR-Egger method and IVW is that MR-Egger considers the presence of the intercept term in the regression, and in addition it uses the inverse of the ending variance as a weight for the fit[23]. Also we applied the MR Egger intercept test to detect unbalanced horizontal pleiotropy [23]. The weighted median method provides robust estimates even if up to 50% of the genetic variance violates the assumptions [24]. In addition, we performed MR multiplicity residuals and outliers (MR − PRESSO) tests to assess the presence of multiplicity [25], which mainly detects outliers and removes them if any, and re-calculate the MR estimates. In addition, we used leave-one-out analysis to check the robustness and consistency of the results[26]. Finally, Cochran's Q statistic was used in this study to detect heterogeneity.

All MR analyses in this study were performed using R software (version 4.2.1) and the "TwoSampleMR" software package (version 0.5.6).

Results

Identification of instrumental variables

Eleven SNPs significantly associated with ischemic stroke matched in the epilepsy database were used as instrumental variables, and these instrumental variables were calculated to explain approximately 1.87% of the variance in exposure, with F-values greater than 10, indicating that the presence of weak instrumental variables is less likely to bias the results. Details of all instrumental variables can be found in Table 2.

Mendelian randomization results

The results of the IVW model based on the ILAE dataset suggest that ischemic stroke may increase the relative risk of generalized epilepsy occurrence (OR, 1.24; 95% CI, 1.09, 1.42; P = 0.0009), but not the relative risk of focal epilepsy (OR, 1.06; 95% CI, 0.93, 1.20; P = 0.3957). The specific results are shown in Fig. 2, and the related scatter plots are shown in Fig. 3.

The Cochran’Q test showed no heterogeneity between ischemic stroke and generalized epilepsy (p = 0.4506), and the results of the polymorphism test also showed no cross-sectional polymorphism in this study (p = 0.2618). In addition, no polymorphic variants or outliers were detected using the MR-PRESSO method. The results of the sensitivity test analysis of the exclusion showed that the inclusion of SNPs had no significant effect on the results, as shown in Fig. 4.

In the reverse MR analysis from epilepsy to ischemic stroke, 11 SNPs were selected as IVs. (Supplementary Table S1). In the IVW method, we found no relevant evidence of a causal relationship between ischemic stroke and the risk of generalized epilepsy (OR, 1.04; 95% CI, 0.97, 1.12; p = 0.2571) (Supplementary Table S2). This was also confirmed by other MR approaches as shown in Supplementary Table S2. Sensitivity analysis verifies that epilepsy is not heterogeneous for ischemic stroke (p = 0.4354). The pleiotropy test indicated no pleiotropy in the analysis (p = 0.6877), and the results were stable when excluding any single SNP. Collectively, there is currently no evidence for an inverse causal relationship between ischemic stroke and generalized epilepsy.

Discussion

In this study, we used the two − sample MR framework by applying public data from different GWAS to decipher the causality between ischemic stroke and the risk of epilepsy. We provide genetic evidence that ischemic stroke leads to an elevated risk of developing epilepsy. In addition, no heterogeneity or horizontal pleiotropy was detected. Our findings suggest the need for timely intervention after ischemic stroke to prevent epilepsy.

At one time the relationship between stroke and epilepsy was thought to be bidirectional. One study showed that a quarter of all epilepsy patients over the age of 65 experience a stroke seizure[27]. A retrospective cohort study published in 2015 included a total of 21,035 epilepsy cases and 16,638 controls, and the results indicated that the incidence of stroke after seizures was 2.5 times higher in adults than in controls. After adjusting for comorbidities and other factors, the risk of stroke was as high as 60% in epilepsy cases compared to controls[28]. As a subtype of stroke, ischemic stroke is one of the leading causes of death and a major cause of disability worldwide. This study is the first to use Mendelian randomization to discuss the potential causal relationship between ischemic stroke and epilepsy.

The occurrence of epilepsy after ischemic stroke involves a variety of mechanisms, firstly, the ischemia and hypoxia caused by stroke leads to sodium pump failure, Na+ inward flow causes cell depolarization, when the depolarization reaches a certain level, calcium channels are activated and Ca2+ enters intracellularly in large quantities thus leading to abnormal neuronal discharge[29], meanwhile, neurotransmitter dysregulation is also an important cause of seizures, and studies have shown that with ischemia and hypoxia After glutamate and other excitatory neurotransmitters increase, which increases membrane excitability[30], and in addition, GABA inhibitory neurotransmitters undergo degeneration, weakening the inhibitory effect and lowering the membrane excitation threshold [30], which together cause abnormal neuronal firing leading to epilepsy. Stroke-induced disruption of the blood-brain barrier causes alterations in astrocyte gene expression, resulting in a decrease in the ability of astrocytes to take up K+. Elevated extracellular K+ concentrations can cause excessive neuronal firing, leading to seizures [31].

Several advantages exist in the present study. First, this study primarily used a large sample of genome-wide association studies (GWAS), which allows for a better overall analysis of stroke events. Second, in the experimental design of MR, because SNPs are randomly assigned, bias due to reverse causality or confounding factors is greatly reduced, which is not available in observational studies. Finally, we also performed several sensitivity analyses to test the consistency of the results and to correct for possible pleiotropy.

However, there are some limitations to this study. First, the epilepsy dataset provided by the ILAE used in this study had approximately 14% of non-European descent, Therefore, population stratification is likely to introduce biases here. Second, the results of the study were limited to Europeans and could not be generalized to other populations. Finally, because ischemic stroke and epilepsy are heterogeneous disorders with different genetic factors and may be differentially associated with each other, larger subgroup MR studies may be needed in the future to further explore causal relationships.

Conclusions

In conclusion, this MR study provides genetic evidence that ischemic stroke is associated with an increased risk of generalized epilepsy, but not focal epilepsy, in a European population. These findings may provide strong evidence for the future prevention and treatment of post-stroke epilepsy.

Declarations

Acknowledgments

The authors sincerely thank the MEGASTROKE consortium and the International League Against Epilepsy Consortium for providing summary-level data.

Authors’ contributions

SC designed the study. WH and XJ analyzed the data and prepared the original draft. TH and MZ revised and edited the paper. HX supervised this work. All authors have read and agreed to the published version of the manuscript.

Funding

Not applicable.

Availability of data and materials

The datasets generated and/or analyzed during the current study are publicly available in the IEU Open GWAS Project repository(http://gwas.mrcieu.au.uk).

Ethics approval and consent to participate

Not applicable. Ethical approval and informed consent for studies included in the analyses was provided in the original publications.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

References

  1. R.S. Fisher, C. Acevedo, A. Arzimanoglou, A. Bogacz, J.H. Cross, C.E. Elger, et al., ILAE official report: a practical clinical definition of epilepsy. Epilepsia 55 (2014) 475-82.
  2. G.B.D.N. Collaborators, Global, regional, and national burden of neurological disorders, 1990-2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet Neurol 18 (2019) 459-480.
  3. S. Lapalme-Remis, and G.D. Cascino, Imaging for Adults With Seizures and Epilepsy. Continuum (Minneap Minn) 22 (2016) 1451-1479.
  4. A. Vezzani, R.S. Fujinami, H.S. White, P.M. Preux, I. Blumcke, J.W. Sander, et al., Infections, inflammation and epilepsy. Acta Neuropathol 131 (2016) 211-234.
  5. J.W. Sander, Y.M. Hart, A.L. Johnson, and S.D. Shorvon, National General Practice Study of Epilepsy: newly diagnosed epileptic seizures in a general population. Lancet 336 (1990) 1267-71.
  6. R.D. Thijs, R. Surges, T.J. O'Brien, and J.W. Sander, Epilepsy in adults. Lancet 393 (2019) 689-701.
  7. M. Endres, M.A. Moro, C.H. Nolte, C. Dames, M.S. Buckwalter, and A. Meisel, Immune Pathways in Etiology, Acute Phase, and Chronic Sequelae of Ischemic Stroke. Circ Res 130 (2022) 1167-1186.
  8. A.M. Feyissa, T.F. Hasan, and J.F. Meschia, Stroke-related epilepsy. Eur J Neurol 26 (2019) 18-e3.
  9. S. Zou, X. Wu, B. Zhu, J. Yu, B. Yang, and J. Shi, The pooled incidence of post-stroke seizure in 102 008 patients. Top Stroke Rehabil 22 (2015) 460-7.
  10. J.Z. Wang, M.V. Vyas, G. Saposnik, and J.G. Burneo, Incidence and management of seizures after ischemic stroke: Systematic review and meta-analysis. Neurology 89 (2017) 1220-1228.
  11. J.G. Burneo, J. Fang, G. Saposnik, and N. Investigators of the Registry of the Canadian Stroke, Impact of seizures on morbidity and mortality after stroke: a Canadian multi-centre cohort study. Eur J Neurol 17 (2010) 52-8.
  12. C.W. Huang, G. Saposnik, J. Fang, D.A. Steven, and J.G. Burneo, Influence of seizures on stroke outcomes: a large multicenter study. Neurology 82 (2014) 768-76.
  13. M.R. Keezer, G.S. Bell, A. Neligan, J. Novy, and J.W. Sander, Cause of death and predictors of mortality in a community-based cohort of people with epilepsy. Neurology 86 (2016) 704-12.
  14. S. Fazel, A. Wolf, N. Langstrom, C.R. Newton, and P. Lichtenstein, Premature mortality in epilepsy and the role of psychiatric comorbidity: a total population study. Lancet 382 (2013) 1646-54.
  15. H.E.K. Sundelin, T. Tomson, J. Zelano, J. Soderling, P. Bang, and J.F. Ludvigsson, Pediatric Ischemic Stroke and Epilepsy: A Nationwide Cohort Study. Stroke 52 (2021) 3532-3540.
  16. D.A. Lawlor, R.M. Harbord, J.A. Sterne, N. Timpson, and G. Davey Smith, Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med 27 (2008) 1133-63.
  17. G. Davey Smith, and G. Hemani, Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet 23 (2014) R89-98.
  18. C.A. Emdin, A.V. Khera, and S. Kathiresan, Mendelian Randomization. JAMA 318 (2017) 1925-1926.
  19. R. Malik, G. Chauhan, M. Traylor, M. Sargurupremraj, Y. Okada, A. Mishra, et al., Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes. Nat Genet 50 (2018) 524-537.
  20. E. International League Against Epilepsy Consortium on Complex, Genome-wide mega-analysis identifies 16 loci and highlights diverse biological mechanisms in the common epilepsies. Nat Commun 9 (2018) 5269.
  21. S. Burgess, S.G. Thompson, and C.C.G. Collaboration, Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol 40 (2011) 755-64.
  22. J. Bowden, M.F. Del Greco, C. Minelli, G. Davey Smith, N. Sheehan, and J. Thompson, A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med 36 (2017) 1783-1802.
  23. J. Bowden, G. Davey Smith, and S. Burgess, Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol 44 (2015) 512-25.
  24. J. Bowden, G. Davey Smith, P.C. Haycock, and S. Burgess, Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol 40 (2016) 304-14.
  25. M. Verbanck, C.Y. Chen, B. Neale, and R. Do, Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet 50 (2018) 693-698.
  26. Q.F. Gronau, and E.J. Wagenmakers, Limitations of Bayesian Leave-One-Out Cross-Validation for Model Selection. Comput Brain Behav 2 (2019) 1-11.
  27. M. Holtkamp, E. Beghi, F. Benninger, R. Kalviainen, R. Rocamora, H. Christensen, et al., European Stroke Organisation guidelines for the management of post-stroke seizures and epilepsy. Eur Stroke J 2 (2017) 103-115.
  28. B.B. Wannamaker, D.A. Wilson, A.M. Malek, and A.W. Selassie, Stroke after adult-onset epilepsy: a population-based retrospective cohort study. Epilepsy Behav 43 (2015) 93-9.
  29. M. Gooshe, A.H. Abdolghaffari, A.R. Aleyasin, L. Chabouk, S. Tofigh, G.R. Hassanzadeh, et al., Hypoxia/ischemia a key player in early post stroke seizures: modulation by opioidergic and nitrergic systems. Eur J Pharmacol 746 (2015) 6-13.
  30. S.M. Mizielinska, S.M. Greenwood, H. Tummala, and C.N. Connolly, Rapid dendritic and axonal responses to neuronal insults. Biochem Soc Trans 37 (2009) 1389-93.
  31. U. Heinemann, S. Gabriel, R. Jauch, K. Schulze, A. Kivi, A. Eilers, et al., Alterations of glial cell function in temporal lobe epilepsy. Epilepsia 41 Suppl 6 (2000) S185-9.

tables

Table 1 | Details of studies and datasets used for analyses.

Data Source

Phenotype

Sample size

Cases

Population

References

 


The MEGASTROKE consortium

ischemic stroke

440328

34217

European

Malik R et al.

 




The International League Against Epilepsy Consortium

focal epilepsy

39348

9671

86% European

Abou-Khalil B et al.



generalized epilepsy

33446

3769



 

Table 2 | Association strength of SNP allele frequencies and effect alleles with ischemic stroke and epilepsy

SNP

EA

OA

EAF

R2(%)

F

Beta

SE

P-value

rs12445022

A

G

0.31 

0.16 

694.33 

0.061 

0.010 

1.28 ×10-10

rs17035646

A

G

0.41 

0.14 

610.53 

0.054 

0.009 

1.34 ×10-09

rs2107595

A

G

0.23 

0.20 

890.35 

0.076 

0.010 

9.25 ×10-14

rs3184504

C

T

0.55 

0.28 

1233.78 

-0.075 

0.010 

2.17 ×10-14

rs35436

T

C

0.38 

0.12 

509.69 

-0.050 

0.009 

3.21 ×10-08

rs42039

T

C

0.23 

0.15 

667.45 

-0.066 

0.011 

6.55 ×10-09

rs4932370

A

G

0.33 

0.12 

527.51 

0.052 

0.009 

2.88 ×10-08

rs4959130

A

G

0.14 

0.16 

722.82 

0.083 

0.014 

2.83 ×10-09

rs6825454

C

T

0.31 

0.14 

597.66 

0.056 

0.009 

7.43 ×10-10

rs6847935

T

A

0.33 

0.27 

1192.02 

0.078 

0.010 

3.50 ×10-16

rs9526212

G

A

0.76 

0.14 

605.95 

0.062 

0.010 

9.19 ×10-10

SNP , single-nucleotide polymorphism; EA, effect allele; EAF , effect allele frequency; R2, percentage of the variation explained by the SNP; F , F statistic; Beta, the per-allele effect on cannabis use; SE, standard error of Beta; P-value is for the genetic association