A Fast Method for Estimating Statistical Power of Multivariate GWAS in Real Case Scenarios: Examples from the Field of Imaging Genetics

In GWAS of imaging phenotypes (e.g., by the ENIGMA and CHARGE consortia), the growing number of phenotypes considered presents a statistical challenge that other fields are not experiencing (e.g. psychiatry and the Psychiatric Genetics Consortium). However, the multivariate nature of MRI measurements may also be an advantage as many of the MRI phenotypes are correlated and multivariate methods could be considered. Here, we compared the statistical power of a multivariate GWAS versus the current univariate approach, which consists of multiple univariate analyses. To do so, we used results from twin models to estimate pertinent vectors of SNP effect sizes on brain imaging phenotypes, as well as the residual correlation matrices, necessary to estimate analytically the statistical power. We showed that for subcortical structure volumes and hippocampal subfields, a multivariate GWAS yields similar statistical power to the current univariate approach. Our analytical approach is as accurate but ~ 1000 times faster than simulations and we have released the code to facilitate the investigation of other scenarios, may they be outside the field of imaging genetics.