Expression quantitative trait loci (eQTL) analysis was performed for all SNPs in each cell type within 100,000 bases of a differentially expressed gene overlapping a GWAS loci and within an ATAC-seq peak from that cell type. tensorQTL106 was used to perform the analysis using the map_cis function with default parameters. Variants with minor allele frequency < 0.05 were excluded from the analysis. Sex, age, and ethnic origin were used as covariates.