The SNP-based association is performed by the standard allelic test and logistic regression to obtain the P values for individual SNPs, and PLINK software (version 1.07) was used for analysis [58]. Meanwhile, a list of SNP pairs in linkage disequilibrium (LD) (r 2 > 0.2) is calculated for the gene-based association test.