Haplotype associations were tested by the program HBAT in the FBAT package (Laird et al 2000). We adopted the global haplotype association test with Monte Carlo simulation, i.e., repeated random sampling, to compute an exact global p-value for multi-haplotype based association (http://www.biostat.harvard.edu/~fbat/default.html).