37/hg19 (http://genome.ucsc.edu/). For the 100 most significantly associated SNPs from logistic regression, flanking regions of 50kb on either side were also searched. Finally, in order to identify independent association signals in our data, we applied the clumping procedure in PLINK2 [21], taking all SNPs associated with our phenotype with p≤0.001, termed index SNPs, and identifying SNPs that are in LD with the index SNP (r2 > 0.5) within a sliding window of 250kb.