Eigenstrat, implemented in Eigensoft, was used to identify outliers in the population who were subsequently removed from the analysis. The first two eigenvectors were used as covariates to correct for any remaining population structure. As a result, genetic data from 404 individuals (304 MSM and 100 DU) were used for further analysis.