To address potential population stratification, we used the PLINK nearest neighbor method and filtered out individuals outside three standard deviations of the first four multidimensional scaling (MDS) factors. The first four MDS factors were also adjusted in the GWAS analyses to control for population substructure, as when the first 10 MDS factors were adjusted, the genomic inflation factor λs were not improved.
