To minimize the possibility of false positives, the SNP-expression probe pairs with corresponding to a Bonferroni-corrected P value of <0.05 (i.e., uncorrected P<0.05/[30,465×534,404] = 3.1×10−12) were also examined using Spearman’s rank correlation in a similar method as described above. The SNP-probe pairs with Bonferroni-corrected P value of the average Spearman’s rank correlation <0.05 (i.e., uncorrected P<3.1×10−12) were considered significant.