Bayesian colocalization of GWAS and eQTL signals reveals cell type-specific genes and regulatory variants for susceptibility to subtypes of ischemic stroke

Abstract

A colocalization analysis of genome-wide association study (GWAS) signals and expression quantitative trait loci (eQTL) was conducted to pinpoint target genes and their regulatory nucleotide variants for subtypes of ischemic stroke. We utilized GWAS data from prominent meta-analysis consortia (MEGASTROKE and GIGASTROKE) and single-cell eQTL data in brain and blood tissues to enhance accuracy and minimize noise inherent in bulk RNA-seq. Employing Bayesian colocalization methods, we identified ten shared loci between GWAS and eQTL signals, targeting five eGenes. Specifically, RAPH1 and ICA1L were discovered for small vessel stroke (SVS), whereas SCYL3, CAV1, and CAV2 were for cardioembolic stroke (CS). However, no findings have been made for large artery stroke. The exploration and subsequent functional analysis of causal variants within the colocalized regions revealed their regulatory roles, particularly as enhancer variants (e.g., rs144505847 and rs72932755 targeting ICA1L; rs629234 targeting SCYL3; rs3807989 targeting CAV1 and CAV2). Notably, our study unveiled that all eQTL for CS were identified in oligodendrocytes, while those for SVS were across excitatory neurons, astrocytes, and oligodendrocyte precursor cells. This underscores the heterogeneous tissue-specific genetic factors by subtypes of ischemic stroke. The study emphasizes the need for intensive research efforts to discover causative genes and variants, unravelling the cell type-specific genetic architecture of ischemic stroke subtypes. This knowledge is crucial for advancing our understanding of the underlying pathophysiology and paving the way for precision neurology applications.