The population analysis software Structure 2.3.3 (developed by J. Pritchard from University of Chicago) was utilized to check for population genetic substructure in our sample using the 23 markers genotyped as described above. The first two ancestry proportions were adjusted for in the association analysis.