eQTL Catalogue analysis
Had a closer look at eQTL that match the DecodeME signal in this region. This uses data from experiments that measured RNA (gene expression) and linked it to the DNA of the participants. eQTL are the DNA letters that change a gene's expression. Often a GWAS hit cannot be matched with an eQTL but if it does, it helps to show which gene it changes and in what cell type or tissue.
I used the eQTLCatalogue which offers multiple of these datasources in the same format through their API. Their latest version also includes the full data for the Genotype-Tissue Expression (GTEx) Project version 10, a popular database that tested multiple tissue that was also used in the DecodeME preprint. I've included all other datasets related to brain and the immune system using the overview they have available here:
eQTL-Catalogue-resources/data_tables/dataset_metadata_r8.tsv at master · eQTL-Catalogue/eQTL-Catalogue-resources
I've only looked at gene expression (quant_method ="ge" or "microarray"). The data is large but it's in parquet and tabix formats that allows you to download only the locus that you need rather than the full file. So that's how I ended up with 259 datasets for 36 genes, giving 3,705 gene × dataset tests. I've tested colocalization with the DecodeME data in this region, which gives the probability (PP.H4) that both results are driven by the same (single) causal SNP. I've pasted the top findings below, the full results are in the excel file attached.
There are multiple datasets x genes that reach a high PP.H4 but given the large number of tests we should be wary of false positives. Most of the eQTL datasets also have a small sample size (often below 300) and findings below the significance threshold of 5*10^-8.
Interestingly, the best match was for LPS-triggered monocytes in the
Kim-Hellmuth 2017 study. I've plotted the data below
This is, however, only one of the 21 tests they did using different microarray probes for the gene and e.g. a test 6 hours instead of 90 minutes after LPS administration. These others didn't match that well. There's also another LPS-triggered monocyte dataset (
Fairfax 2014) that didn't gave a good match.
The GTEx data for CSE1L in the nucleus accumbens (part of the striatum and basal ganglia) also forms a good match. Nonetheless, plotting the data I often have the impression that all eQTLs are driven by the broad stroke to the right of the core SNPS from the DecodeME finemapping. There's some noise in that so if test multiple of these datasets, sometimes it shifts a bit to the left and seems to match the DecodeME signal. So despite the high PP.H4 I remain skeptical that both point to the same causal SNP and that we really have a match here.
ARFGEF2, the closest gene, is notably absent from the top eQTL results here but based on its position and uncertainty about whether the eQTL really match, I wouldn't write it off.
The core SNPs of the DecodeME signal lie in the introns of a lncRNA (it's called ENSG00000294533) which isn't included in these datasets. But I figured it wouldn't mess up the results (it would simply be an intermediary) unless it indirectly influences a gene that is really far from the locus.
Hope this is useful to someone. It didn't really give me much clarity, unfortunately.