I had a look at the region using ENCODE data on chromatin openness. I used the SCREEN API which organizes the evidence into candidate regulatory elements (CCREs). I looked at 11 core SNPs for the Decode hit (which have more than 80% of the probability of the signal) and found that 7 were in a CCRE.
Next, I could look up scores for these 7 CCRE for various experiments and biosamples. It includes z-scores for open DNA (DNase or ATAC), the active-enhancer mark H3K27ac, the promoter mark H3K4me3, and binding of CTCF (a protein that anchors DNA loops). Using AI I also tried to group the biosamples in meaningful categories (broad_group and fine_group) such as immune, CNS, etc. This help with screening and interpretability but be aware that these are my own categories, not official labels from ENCODE.
Overall, this region doesn’t stand out in chromatin openness and other markers, not more than what would be expected by chance. The most notable results were for immune cells for the SNPs 20:48911205: A: C and 20:48909956:T:G. HL-60 refers to a leukemia cell line.
The graph below plots the CRE position and distance with the location of our core SNPs added. The ARFGEF2 transcription start site is also shown as a comparison of how an active region looks like. The average across all biosamples didn’t show anything interesting in this region but when I restrict to the immune system data, there seems to be a peak at the location where 20:48911205: A: C and 20:48909956:T:G lie.
The average for central nervous system samples didn’t show peaks in our SNPs, so I had a look at a more refined dataset: the Li et al. 2023 brain atlas that used sn-ATAC-seq normalized per 107 cell types. Here I found a minor peak at 20:48917118:T:C (closer to ARFGEF2) that seems to point to GABAergic neurons from the basal ganglia. The FOXP2 finding is interesting because it
might be point to eccentric medium spiny neurons found in the RNA analyses of DecodeME (but not so much in the ATAC-analysis).
So, what does all of this mean? I’m not sure. None of hypothesis’s currently seems more convincing that the others. I plan to look closer the eQTL, coloc and AphaGenome predictions but I suspect we’ll need some experimental data from this region to know what it actually points to. It’s the biggest hit from DecodeME so I think it deserves a follow-up!