We removed samples with genotype call rates below 97%, excess heterozygosity, duplicates, samples genetically related to other individuals in the cohort or sample mix-up (PIHAT > 0.1875). If a sex discrepancy was detected, the sample was removed unless the discrepancy was safely resolved. To detect population outliers of non-European ancestry (>6 SD from European population mean), principal component analysis (PCA) was conducted using SMARTPCA from EIGENSOFT 6.1.4.
We removed variants with a call rate < 95% or that grossly deviated from Hardy–Weinberg equilibrium in controls (P-value ≤ 1 × 10–4), markers with a different missing rate between case and control (P-value < 5 × 10–4 for the difference) or minor allele frequency (MAF) below 0.01. Imputation was carried out using Haplotype reference consortium (HRC) panel in Michigan Imputation servers