using Illumina Sentrix HumanHap300 BeadChip (ILMN300K, including 307,981 autosomal SNPs) according to the manufacturers' instructions (Dixon et al. 2007; Moffatt et al. 2007). Before analysis, we excluded 4050 SNPs with a call rate <95%, 96 SNPs with Hardy-Weinberg equilibrium P-value <10−6, and 4310 SNPs with MAF < 2% from ILMN100K (a total of 8313 SNPs excluded), and 3921 SNPs with call rate <95%, 34 SNPs with Hardy-Weinberg equilibrium P-value <10−6, and 483 SNPs with MAF < 2% from ILMN300K (a total of 4420 SNPs excluded).
