eQTLs were computed using a linear model on the imputed genotype dosages using Matrix EQTL69 (Shabalin, 2012). If the most significant variant was filtered out in this sample due to QC, the next available variant with the lowest p value and in LD (r2>0.9) with the most significant variant was tested instead. The gene expression data were adjusted for sex, age, post-mortem interval, RIN, batch, and decedent phenotype (AD or control)