The enrichment analysis was performed by assessing the association between AD and whole gene eQTLs SNP-set in the R package gskat. Gskat performs a set-based test for the effects of a SNP set in association studies for both quantitative and discrete phenotypes using the generalized estimation equation approach. Specifically, a Kernel Machine (KM) estimating based statistic was constructed to test for the association between an AD and other alcohol related phenotype and a SNP set.