Fine-mapping

We fine-mapped each region from the GWASs where a variant reached p < 1-6. Each region was fine-mapped with SuSiE 0.8.1.0545 and FINEMAP v1.4_0510.

We used a 3-megabase window (+- 1.5M) around each lead variant and merged overlapping regions into one. After merging, a handful of regions became too large to computationally handle with SuSiE. For such regions, we only merged two overlapping pieces when the LD between the two lead variants was r2 > 0.2. When LD was less than that, we made each of the two overlapping regions non-overlapping by splitting the overlap in half.

The codebase, workflow, and inputs we used for R5 fine-mapping is here: https://github.com/FINNGEN/finemapping-pipeline/releases/tag/r5

Integration to PheWeb

The "Credible Sets"-table on a phenotype page in the PheWeb browser shows the SuSiE-fine-mapped credible sets of that phenotype. The variant shown per credible set is the maximum PIP (posterior inclusion probability) variant of that credible set. In addition to the causal variants, variants that were in sufficient LD (pearsonr^2 > 0.05), had a small enough p-value (pval < 0.01), and were close enough to the lead variant (distance to lead variant < 1.5 megabases) were clumped together with the credible set. Variants have been compared against GWAS Catalog and annotated. The LD grouping, annotation and GWAS Catalog comparison were done using the autoreporting pipeline.

The columns of the table are explained below:

Last updated