In the SAGE dataset, a total of 646 SNPs showed association with maxdrinks at p < 1 × 10−4; however, none of these variants met criteria for genome-wide significance (Supplementary table 2, supplementary figure 3). The observed number of SNPs is almost 1.5 times more than the number of SNPs expected to show suggestive association at this level by chance (p < 1 × 10−4). rs67666182 near SNX16 (sorting nexin 16) on chromosome 8 showed suggestive association (9.9 × 10−5 <p<7.1 × 10−7). Other SNPs in high linkage disequilibrium with rs67666182 (r2=0.8 D′=0.95) also provided evidence of association (Supplementary table 2). Several SNPs in cytochrome P450, family 4, subfamily F, polypeptide 8 (CYP4F8), muskelin 1, intracellular mediator containing kelch motifs (MKLN1), LIM domain only 1 (LMO1), phenylalanyl-tRNA synthetase 2, mitochondria (FARS2) and SH3-domain GRB2-like 3 (SH3GL3) were also associated with maxdrinks (p < 1.0 × 10−4).
