For the allelic test, expression values were binned by rs1344706 allele and the means of the bins were then compared52. Each individual’s expression value is used twice: for homozygotes, twice in the same bin, for heterozygotes, once in each. As a result, the variation in the homozygotes may be less than that in heterozygotes, and we therefore compare the two group means using the unpaired t-test with Welch’s correction to account for possible heteroscedastic variances.