Abstract
Candidate gene studies have identified genetic variants associated with clinical outcomes following aneurysmal subarachnoid haemorrhage (aSAH), but no genome-wide association studies have been performed to date. Here we report the results of the discovery phase of a two-stage genome-wide meta-analysis of outcome after aSAH. We identified 157 independent loci harbouring 756 genetic variants associated with outcome after aSAH (p < 1 × 10−4), which require validation. A single variant (rs12949158), in SPNS2, achieved genome-wide significance (p = 4.29 × 10−8) implicating sphingosine-1-phosphate signalling in outcome after aSAH. A large multicentre international effort to recruit samples for validation is required and ongoing. Validation of these findings will provide significant insight into the pathophysiology of outcomes after aSAH with potential implications for treatment.
Similar content being viewed by others
Introduction
Aneurysmal subarachnoid haemorrhage (aSAH) is a devastating form of stroke with the worst outcomes and highest socioeconomic burden of any stroke type [1]. The pathophysiology of neurological injury following aSAH is incompletely understood. The mechanism is thought to be multifactorial with the initial surge in intracranial pressure following haemorrhage combined with the presence of blood breakdown products in the cerebrospinal fluid leading to a pattern of injury characterised by inflammation, cerebral vasospasm, microthrombosis, oxidative injury and cortical spreading depression [2,3,4]. Despite multiple clinical trials, nimodipine is the only therapeutic agent to improve outcomes [5]. It is our incomplete understanding of the mechanisms underlying neurological injury that is, at least in part, responsible for the lack of therapeutic innovation to improve outcomes.
The best outcome prediction model after aSAH, utilising clinical, demographic and imaging characteristics only explains up to 31% of the variation in outcome following aSAH [6]. Consequently, a large proportion of variation in outcome following aSAH is unexplained. There is a growing body of evidence from candidate gene studies that genetic background accounts for a proportion of this unexplained variation [7, 8]. However, no genome-wide analysis has been performed. Such a study would have the potential to provide valuable insights into the mechanisms underlying neurological injury following aSAH by identifying, as yet unstudied, genes associated with outcome and thus novel targets for therapeutic intervention.
In 2018, the HATCH consortium highlighted the need to better understand the pathophysiological mechanisms underlying outcome and proposed a large multicentre genetic analysis of outcome following aSAH [2]. The HATCH consortium has developed this proposal into an international collaboration to undertake a two-stage (discovery and validation) genome-wide association (GWA) study of outcome following aSAH, the protocol for which was published in this journal [9].
The aim of this manuscript is to (1) report the completion of the discovery stage of the study including preliminary results and (2) raise awareness of the study to recruit further samples for the validation stage.
Methods
This is a two-stage (discovery and validation) GWA meta-analysis of outcome following aSAH. The results of the discovery analysis are reported in this manuscript. All analyses were performed according to the published protocol [9]. The study has both national ethical (REC 19 SC 0485) and institutional (ERGO 49253) approval.
For the discovery analysis, individuals were identified from (1) six studies from the HATCH consortium network and (2) the UK Biobank, a major biomedical database with extensive genetic and clinical data, previously described in detail[10] (application number 49305).
In the HATCH dataset, the primary outcome was the modified Rankin Scale (mRS) [11, 12] or Glasgow Outcome Scale (GOS) [13, 14] dichotomised into good (mRS 0–2, GOS 4–5) and poor (mRS 3–6, GOS 1–3) outcomes in the first two years following aSAH. The mRS and/or GOS are not available in the UK Biobank and, therefore, a measure of cognitive performance, psychomotor reaction time, was used since cognition is highly correlated with mRS/GOS following aSAH [15] and reaction time is significantly slower in aSAH cases compared to controls in the UK Biobank [16]. Reaction times were ranked from fastest to slowest and then the UK Biobank was dichotomised into good (faster) and poor (slower) outcomes, generating an equivalent proportion of good outcome individuals to the HATCH dataset.
Genotype information from eligible patients underwent quality control and imputation as required (see protocol for details [9]). Within individual cohorts, genetic variants were tested for association with dichotomised outcome using multivariable logistic regression under an additive model, controlling for confounding variables (age and genetic ancestry). A fixed effects meta-analysis was performed to determine each genetic variant’s overall effect size and significance. The meta-analysis was performed on all datasets and repeated in the HATCH dataset alone (i.e. excluding the UK Biobank given its alternative outcome metric). Independent loci were identified for validation using a clumping procedure to group single nucleotide polymorphisms (SNPs) in linkage disequilibrium (LD) (R2 > 0.2) and within 250 kb of each other. Index SNPs with suggestive significance (p < 1 × 10−4) were selected for validation. The threshold for genome-wide significance was p ≤ 5 × 10−8 and all analyses were performed using PLINK, STATA (StataCorp. 2011. Stata Statistical Software: Release 16. College Station, TX: StataCorp LP), wANNOVAR [17] and FUMA [18].
Results
A total of 2489 samples were used for the discovery analysis following quality control. These samples were drawn from six datasets from the HATCH consortium [19,20,21,22] (n = 1685 patients) and 804 individuals from the UK Biobank. After dichotomisation of mRS/GOS within the HATCH consortium data, 1382 (82.0%) patients were classified as good outcome and 303 (18.0%) as poor outcome. Based on reaction times in the UK Biobank, 653 (81.2%) individuals were classified as good outcome and 151 (18.8%) as poor outcome. Table 1 details the demographics and other characteristics of the included datasets.
Analysis of samples from the HATCH consortium (n = 1685) identified 403 SNPs associated with clinical outcome (p < 1 × 10−4) within 97 independent loci after LD-based SNP clumping (Fig. 1A and Supplementary Table 1A and B). No genetic variants reached genome-wide significance.
Including all seven datasets (n = 2489) 85 independent loci were identified from 406 SNPs, associated with clinical outcome (p < 1 × 10−4) (Fig. 1B and Supplementary Table 2A and B). A single variant, rs12949158, reached genome-wide significance (p = 4.29 × 10−8). rs12949158 is located on chromosome 17 in an intronic region of the sphingolipid transporter 2 (SPNS2) gene, which codes for the major transporter of sphingosine-1-phosphate (S1P) (Supplementary Table 2A). The rs12949158 SNP was only genotyped in two datasets [UK Biobank (n = 744) and Korean datasets (n = 89)] (Fig. 1C). Arrays used in the other datasets did not include rs12949158 or any other SNPs with sufficient LD to allow reliable imputation. The rs12949158 alternate A allele was associated with an increased risk of poor outcome with an odds ratio of 2.15 (95% confidence interval 1.63–2.82). The association of rs12949158 with psychomotor reaction time was specific to aSAH, since in a previously published UK Biobank control cohort matched to the same aSAH population[16], this relationship was absent (p = 0.55).
Including both analyses, a total of 157 independent loci were identified from 756 unique SNPs and these will be taken forward for validation.
Discussion
In this discovery genome-wide meta-analysis, we identified 157 independent loci from 756 unique SNPs associated with outcome following aSAH (p < 1 × 10−4) for validation. We also report that the rs12949158 alternate A allele, located within the SPNS2 gene, was associated with an increased risk of poor outcome after aSAH (OR 2.15 95% CI 1.63–2.82) with genome-wide significance (p = 4.29 × 10−8). Although one possible alternative explanation is that rs12949158 associates with psychomotor reaction time independent of aSAH, this association was not observed in control individuals in the UK Biobank.
The genome-wide significant rs12949158 finding is not conclusive and requires validation. Firstly, the rs12949158 genotype was only typed in a subset of the discovery cohort (n = 833). Secondly, the finding is primarily driven by the UK Biobank, which uses an outcome measure that is different from the other datasets (psychomotor reaction time). In the second stage, a customised genotyping array will be used to directly capture all variants targeted for validation including rs12949158.
The rs12949158 variant is intronic, located within the gene SPNS2, a member of the S1P signalling pathway. While variation in a gene intron does not guarantee that the same gene is involved, intron-mediated enhancement of gene expression is increasingly recognised and the intronic variation is most likely to regulate the closest gene [23]. Moreover, the S1P signalling pathway is a biologically plausible candidate to influence outcomes after aSAH since S1P has been implicated in neurological injury following stroke via activation of S1P receptors (S1PR) leading to microglial activation, neuronal death, inflammation and blood–brain barrier disruption [24, 25]. Specifically, after human aSAH, S1P was found to be elevated in the cerebrospinal fluid where its concentration correlated with haemorrhage volume [26] and worse neurological outcome [26]. A possible mechanism linking S1P to clinical outcome is provided by studies showing that S1P induces cerebral vasospasm in canine basilar artery in vitro and in vivo [27] and in murine basilar artery in vitro via S1PR3 [28]. This pathway is of particular interest since S1PR-modulating drugs, already licensed in other neurological conditions (e.g. fingolimod) have been shown to be neuroprotective in ischemic stroke [29] and intracerebral haemorrhage [30] in humans and could be re-purposed for aSAH if this finding is validated.
The study’s population was biased towards participants with a good outcome almost certainly because individuals with a poor prognosis were less likely to be recruited. However, this study’s aim was to better understand the pathophysiological mechanisms underlying outcomes in survivors with a view to developing treatments to improve outcomes. Individuals dying in the acute phase or early after admission will unfortunately be unlikely to benefit from such interventions. Hence, while this study was biased towards participants with a good outcome, these are the individuals most likely to benefit from the findings.
This discovery study achieved its prespecified sample size in a timely manner and represents a highly successful international collaboration generating interesting results to take forward to validation. The total target sample size including the validation cohort is 5000 [9] which would be powered to detect common variation (minor allele frequency (MAF = 0.4) with an effect size of 1.39 and rare variation (MAF = 0.1) with an effect size of 1.66 at genome-wide significance. Recruitment is ongoing for validation and investigators wishing to collaborate using either retrospective or prospective data can find further information on the study website. The study has been designed to maximise inclusivity. It is open to any investigator able to provide the following de-identified biosamples or data from patients with aSAH: genome-wide genotype information (or DNA/cellular sample for genotyping), mRS/GOS within two years of haemorrhage, age, sex and evidence of institutional review board approval. In addition, funding will be provided to facilitate genotyping where local funding is not already in place. Where data availability allows we will also explore whether significant genetic variants mediate the effect of other factors, known to influence outcome after aSAH such as clinical and radiological features.
Data Availability
Study data will be available from the authors subject to institutional agreements and ethical approvals.
References
Taylor TN, Davis PH, Torner JC, Holmes J, Meyer JW, Jacobson MF. Lifetime cost of stroke in the United States, (in eng). Stroke. 1996;27(9):1459–66.
Bulters D, et al. Haemoglobin scavenging in intracranial bleeding: biology and clinical implications, (in eng). Nat Rev Neurol. 2018;14(7):416–32. https://doi.org/10.1038/s41582-018-0020-0.
Zolnourian A, Galea I, Bulters D. Neuroprotective Role of the Nrf2 Pathway in Subarachnoid Haemorrhage and Its Therapeutic Potential, (in eng). Oxid Med Cell Longev. 2019;2019:6218239. https://doi.org/10.1155/2019/6218239.
Macdonald RL, Pluta RM, Zhang JH. Cerebral vasospasm after subarachnoid hemorrhage: the emerging revolution, (in eng). Nat Clin Pract Neurol. 2007;3(5):256–63. https://doi.org/10.1038/ncpneuro0490.
Pickard JD, et al. Effect of oral nimodipine on cerebral infarction and outcome after subarachnoid haemorrhage: British aneurysm nimodipine trial, (in eng). BMJ. 1989;298(6674):636–42.
Jaja BNR, et al. Development and validation of outcome prediction models for aneurysmal subarachnoid haemorrhage: the SAHIT multinational cohort study. BMJ. 2018;360:j5745. https://doi.org/10.1136/bmj.j5745.
Morton MJ, et al. Haptoglobin genotype and outcome after aneurysmal subarachnoid haemorrhage," (in eng). J Neurol Neurosurg Psychiatry. 2020;91(3):305–13. https://doi.org/10.1136/jnnp-2019-321697.
Gallek MJ, Conley YP, Sherwood PR, Horowitz MB, Kassam A, Alexander SA. APOE genotype and functional outcome following aneurysmal subarachnoid hemorrhage, (in eng). Biol Res Nurs. 2009;10(3):205–12. https://doi.org/10.1177/1099800408323221.
Gaastra B, et al. Genome-Wide Association Study of Clinical Outcome After Aneurysmal Subarachnoid Haemorrhage: Protocol, (in eng). Transl Stroke Res. 2022. https://doi.org/10.1007/s12975-021-00978-2.
Sudlow C, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age, (in eng). PLoS Med. 2015;12(3):e1001779. https://doi.org/10.1371/journal.pmed.1001779.
Farrell B, Godwin J, Richards S, Warlow C. The United Kingdom transient ischaemic attack (UK-TIA) aspirin trial: final results. J Neurol Neurosurg Psychiatry. 1991;54(12):1044–54.
Rankin J. Cerebral vascular accidents in patients over the age of 60. II. Prognosis. Scott Med J. 1957;2(5):200–15. https://doi.org/10.1177/003693305700200504.
Jennett B, Bond M. Assessment of outcome after severe brain damage. Lancet. 1975;1(7905):480–4.
Jennett B, Snoek J, Bond MR, Brooks N. Disability after severe head injury: observations on the use of the Glasgow Outcome Scale, (in eng). J Neurol Neurosurg Psychiatry. 1981;44(4):285–93. https://doi.org/10.1136/jnnp.44.4.285.
Mayer SA, et al. Global and domain-specific cognitive impairment and outcome after subarachnoid hemorrhage, (in eng). Neurology. 2002;59(11):1750–8. https://doi.org/10.1212/01.wnl.0000035748.91128.c2.
Gaastra B, Ewbank F, Tapper W, Bulters D, Galea I. Long-Term Cognitive Outcome following Aneurysmal Subarachnoid Haemorrhage. J Stroke Cerebrovasc Dis. 2022;31(1):106184. https://doi.org/10.1016/j.jstrokecerebrovasdis.2021.106184.
Yang H, Wang K. Genomic variant annotation and prioritization with ANNOVAR and wANNOVAR, (in eng). Nat Protoc. 2015;10(10):1556–66. https://doi.org/10.1038/nprot.2015.105.
Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA, (in eng). Nat Commun. 2017;8(1):1826. https://doi.org/10.1038/s41467-017-01261-5.
Bakker MK, et al. Genome-wide association study of intracranial aneurysms identifies 17 risk loci and genetic overlap with clinical risk factors, (in eng). Nat Genet. 2020;52(12):1303–13. https://doi.org/10.1038/s41588-020-00725-7.
Kim H, et al. Cerebral vasospasm after sub-arachnoid hemorrhage as a clinical predictor and phenotype for genetic association study, (in eng). Int J Stroke. 2013;8(8):620–5. https://doi.org/10.1111/j.1747-4949.2012.00823.x.
Li J, et al. Polygenic Risk Scores Augment Stroke Subtyping, (in eng). Neurol Genet. 2021;7(2):e560. https://doi.org/10.1212/NXG.0000000000000560.
Hong EP et al. Genomic Variations in Susceptibility to Intracranial Aneurysm in the Korean Population," (in eng), J Clin Med, 2019;8(2). https://doi.org/10.3390/jcm8020275.
Rose AB. Introns as Gene Regulators: A Brick on the Accelerator, (in eng). Front Genet. 2018;9:672. https://doi.org/10.3389/fgene.2018.00672.
Gaire BP, Choi JW. Sphingosine 1-Phosphate Receptors in Cerebral Ischemia, (in eng). Neuromolecular Med. 2021;2(1):211–23. https://doi.org/10.1007/s12017-020-08614-2.
Sun N, Keep RF, Hua Y, Xi G. Critical Role of the Sphingolipid Pathway in Stroke: a Review of Current Utility and Potential Therapeutic Targets, (in eng). Transl Stroke Res. 2016;7(5):420–38. https://doi.org/10.1007/s12975-016-0477-3.
Männer A, et al. Sphingosine 1-phosphate levels in cerebrospinal fluid after subarachnoid hemorrhage, (in eng). Neurol Res Pract. 2020;2:49. https://doi.org/10.1186/s42466-020-00093-x.
Tosaka M, et al. Sphingosine 1-phosphate contracts canine basilar arteries in vitro and in vivo: possible role in pathogenesis of cerebral vasospasm, (in eng). Stroke. 2001;32(12):2913–9. https://doi.org/10.1161/hs1201.099525.
Salomone S, et al. Analysis of sphingosine 1-phosphate receptors involved in constriction of isolated cerebral arteries with receptor null mice and pharmacological tools, (in eng). Br J Pharmacol. 2008;153(1):140–7. https://doi.org/10.1038/sj.bjp.0707581.
Fu Y, et al. Impact of an immune modulator fingolimod on acute ischemic stroke, (in eng). Proc Natl Acad Sci U S A. 2014;111(51):18315–20. https://doi.org/10.1073/pnas.1416166111.
Fu Y, et al. Fingolimod for the treatment of intracerebral hemorrhage: a 2-arm proof-of-concept study, (in eng). JAMA Neurol. 2014;71(9):1092–101. https://doi.org/10.1001/jamaneurol.2014.1065.
Acknowledgements
This research has been conducted using the UK Biobank Resource under Application Number 49305. The authors acknowledge the use of the IRIDIS High Performance Computing Facility, and associated support services at the University of Southampton, in the completion of this work. We also acknowledge the contributions of the Clinical Research Center, University Hospital and Faculty of Medicine, Geneva.
Funding
BG is funded by is funded by the Royal College of Surgeons, Society of British Neurological Surgeons, Barrow Foundation and Guarantors of Brain in addition to the Institute for Life Sciences, University of Southampton. DJW, HH and IH received funding for recruitment to the GOSH study and genotyping from the Stroke Association and the UCLH NIHR Biomedical Research Centre. SA received funding for data and sample collection from the National Institute of Nursing Research (R01NR004339). YR received funding from the Netherlands Cardiovascular Research Initiative: An initiative with the support of the Dutch Heart Foundation, CVON2015-08 ERASE and from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (PRYSM, grant agreement No. 852173). PB and SM received funding for the AneuX project from SystemsX.ch; it was evaluated by the Swiss National Science Foundation. SB is supported by NIH Grant (K23NS106054).
Author information
Authors and Affiliations
Contributions
All authors contributed to the study conception and design. BG and WT performed the analysis. The first draft of the manuscript was written by BG, WT, IG and DB, and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics Approvals
National (REC 19 SC 0485) and local (ERGO 49253) ethical approval in place.
Consent to Participate
Informed consent was obtained from all individual participants included in the study.
Consent for Publication
Not applicable.
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Will Tapper, Ian Galea and Diederik Bulters are joint senior authors.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Gaastra, B., Alexander, S., Bakker, M.K. et al. A Genome-Wide Association Study of Outcome After Aneurysmal Subarachnoid Haemorrhage: Discovery Analysis. Transl. Stroke Res. 14, 681–687 (2023). https://doi.org/10.1007/s12975-022-01095-4
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s12975-022-01095-4