Interaction association analysis of imputed SNPs in case-control and follow-up studies

A new method is described to assess the interactions of imputed SNPs (single nucleotide polymorphisms) in case-control and follow-up studies, properly incorporating SNP imputation uncertainty in the likelihood model. Using simulation studies and analysis of real data obtained from the Framingham study cohort, we compare the performance of this new method to DOSAGE and NAIVE (also known as Best-Guess) methods, developed and commonly used in the context of single SNP and extended to SNP-by-SNP interaction. The results show that only our new method is unbiased under all examined scenarios regarding allele frequencies, imputation uncertainty degree, and interaction effect size. In addition, our method achieves at least as much power as the other two, and exceeds their statistical power in certain follow-up analysis situations. This method is fast enough to perform Genome Wide Interaction Studies (GWIS) with hundreds of thousands of interactions. By performing an exhaustive simulation study let us to provide recommendations for selecting the most appropriated method depending on MAF, interaction effect size, and uncertainty degree. In general, DOSAGE and our proposed method are recommended in most situations being our method more powerful and accurate when uncertainty and effect increase.