We estimated the loss of function (LoF) burden of each gene on every endpoint.
Then we used the new summarised variables as input in the SAIGE GWAS, replacing the genotype matrix that was used in the regular GWAS.
First, we calculated per individual and gene whether any loss of function variant(s) was present, yielding a matrix with 0 and 1 values ( being the number of individuals and the number of genes).