We included 2,202 endpoints in the analysis. Endpoints with less than 80 cases among the 342,499 samples were excluded, as well as endpoints labeled with an OMIT tag in the endpoint definition file.
For regenie step 1 LOCO prediction computation for each endpoint, we used age, sex, 10 PCs, Finngen 1 or 2 chip or legacy genotyping batch as covariates. For sex-specific phenotypes, sample sex was left out from the covariates.
The null model didn’t initially converge for two phenotypes: PD_DEMENTIA_EXMORE and AD_U_EXMORE. For those two, we excluded legacy batches with less than 10 cases from covariates, after which the nulls converged successfully.
For calculating genetic relatedness in regenie step 1, we included variants 1) imputed with an INFO score > 0.95 in all batches and 2) > 97 % non-missing genotypes and 3) MAF > 1 %. The remaining variants were LD pruned with a 1Mb window and r2 threshold of 0.1. This resulted in a set of 61,289 well-imputed not rare variants for relatedness calculation.
We used a genotype block size of 1,000 in regenie step 1.
We ran association tests with regenie for each of the 2,202 endpoints for each variant with a minimum allele count of 5 among each phenotype’s cases and controls. We used the approximate Firth test for variants with an initial p-value of less than 0.01 and computed the standard error based on effect size and likelihood ratio test p-value (regenie options --firth --approx --pThresh 0.01 --firth-se).