An atlas of single-cell eQTLs dissects autoimmune disease genes and identifies novel drug classes for treatment.

Most variants identified from genome-wide association studies (GWASs) are non-coding and regulate gene expression. However, many risk loci fail to colocalize with expression quantitative trait loci (eQTLs), potentially due to limited GWAS and eQTL analysis power or cellular heterogeneity. Population-scale single-cell RNA-sequencing (scRNA-seq) datasets are emerging, enabling mapping of eQTLs in different cell types (sc-eQTLs). Compared to eQTL data from bulk tissues (bk-eQTLs), sc-eQTL datasets are smaller. We propose a joint model of bk-eQTLs as a weighted sum of sc-eQTLs (JOBS) from constituent cell types to improve power. Applying JOBS to One1K1K and eQTLGen data, we identify 586% more eQTLs, matching the power of 4× the sample sizes of OneK1K. Integrating sc-eQTLs with GWAS data creates an atlas for 14 immune-mediated disorders, colocalizing 29.9% or 32.2% more loci than using sc-eQTL or bk-eQTL alone. Extending JOBS, we develop a drug-repurposing pipeline and identify novel drugs validated by real-world data.