AP-SKAT: highly-efficient genome-wide rare variant association test

BMC Genomics

Software

Takahiro Mimori¹  Kaname Kojima¹  Yosuke Kawai¹  Masao Nagasaki¹  Takanori Hasegawa¹  Kazuharu Misawa¹ 

[1] Department of Integrative Genomics, Tohoku Medical Megabank Organization, Tohoku University, 2-1 Seiryo-machi, Aoba-ku, Sendai, Miyagi, Japan;

关键词: Genome wide association study;    Multiple test;    Rare variants;   

DOI  :  10.1186/s12864-016-3094-3

 received in 2016-01-21, accepted in 2016-09-15,  发布年份 2016

来源: Springer

【 摘 要 】

BackgroundGenome-wide association studies have revealed associations between single-nucleotide polymorphisms (SNPs) and phenotypes such as disease symptoms and drug tolerance. To address the small sample size for rare variants, association studies tend to group gene or pathway level variants and evaluate the effect on the set of variants. One of such strategies, known as the sequential kernel association test (SKAT), is a widely used collapsing method. However, the reported p-values from SKAT tend to be biased because the asymptotic property of the statistic is used to calculate the p-value. Although this bias can be corrected by applying permutation procedures for the test statistics, the computational cost of obtaining p-values with high resolution is prohibitive.ResultsTo address this problem, we devise an adaptive SKAT procedure termed AP-SKAT that efficiently classifies significant SNP sets and ranks them according to the permuted p-values. Our procedure adaptively stops the permutation test when the significance level is outside some confidence interval of the estimated p-value for a binomial distribution. To evaluate the performance, we first compare the power and sample size calculation and the type I error rates estimate of SKAT, SKAT-O, and the proposed procedure using genotype data in the SKAT R package and from 1000 Genome Project. Through computational experiments using whole genome sequencing and SNP array data, we show that our proposed procedure is highly efficient and has comparable accuracy to the standard procedure.ConclusionsFor several types of genetic data, the developed procedure could achieve competitive power and sample size under small and large sample size conditions with controlling considerable type I error rates, and estimate p-values of significant SNP sets that are consistent with those estimated by the standard permutation test within a realistic time. This demonstrates that the procedure is sufficiently powerful for recent whole genome sequencing and SNP array data with increasing numbers of phenotypes. Additionally, this procedure can be used in other association tests by employing alternative methods to calculate the statistics.

【 授权许可】

CC BY   
© The Author(s) 2016

【 图 表 】

12864_2016_2756_Article_IEq2.gif

12864_2017_4186_Article_IEq32.gif

12864_2017_4132_Article_IEq5.gif

12864_2015_1945_Article_IEq4.gif

12864_2017_4186_Article_IEq29.gif

12864_2017_4132_Article_IEq2.gif

12864_2015_2055_Article_IEq8.gif

12888_2016_811_Article_IEq1.gif

12864_2017_4004_Article_IEq7.gif

12864_2016_3105_Article_IEq9.gif

12864_2017_3777_Article_IEq18.gif

12864_2017_3777_Article_IEq17.gif

12864_2016_3098_Article_IEq87.gif