@kaurao found this issue: https://www.biostars.org/p/9490208/ So the idea is to use cat-bgen and then prsice to concatenate al bgen files and then compute the PRS score. This will require R with prsice