In the GWAS (including genome-wide imputation data) and replication study, the chi-square test was applied to a 2-by-2 contingency table in the allele frequency model. The odds ratio (OR) and the con dence interval (CI) were calculated using the major alleles as references. We considered P < 5 x 10 - 8 as the threshold for genome-wide signi cance in the combined analysis.
In the replication stage, we selected 134 SNPs with P values < 10 - 5 and linkage disequilibrium (LD) < 0.9 from the results of the chi-square test in the GWAS using genome-wide imputed data. We additionally selected 49 SNPs located on the functionally interested genes. DigiTag2 assay and TaqMan SNP genotyping assay (Applied Biosystems, Foster City, CA, USA) were used to con rm the genotypes at each SNP. We genotyped 80 cases and 333 controls to validate the GWAS results and for the replication study.