Iterative sure independence screening EM-Bayesian LASSO algorithm for multi-locus genome-wide association studies

• Published in: PLoS computational biology 2017-01, Vol.13 (1), p.e1005357-e1005357
• Main Authors: Tamba, Cox Lwaka, Ni, Yuan-Li, Zhang, Yuan-Ming
• Format: Article
• Language: English
• Subjects: Algorithms; Arabidopsis thaliana; Bayes Theorem; Biology and Life Sciences; Chromosome Mapping - methods; Computer Simulation; Data analysis; Data Interpretation, Statistical; Datasets; DNA Mutational Analysis - methods; Experiments; Funding; Genetic Association Studies - methods; Genetic Markers - genetics; Genetics; Genome-wide association studies; Genomes; Genomics; Germplasm; Iteration (Mathematics); Laboratories; Likelihood Functions; Methods; Models, Genetic; Models, Statistical; Monte Carlo simulation; Optimization theory; Physical sciences; Polymorphism, Single Nucleotide - genetics; Quantitative Trait, Heritable; Research and analysis methods; Studies; Technology application; Variables
• ISSN: 1553-7358; 1553-734X; 1553-7358
• DOI: 10.1371/journal.pcbi.1005357

Genome-wide association study (GWAS) entails examining a large number of single nucleotide polymorphisms (SNPs) in a limited sample with hundreds of individuals, implying a variable selection problem in the high dimensional dataset. Although many single-locus GWAS approaches under polygenic background and population structure controls have been widely used, some significant loci fail to be detected. In this study, we used an iterative modified-sure independence screening (ISIS) approach in reducing the number of SNPs to a moderate size. Expectation-Maximization (EM)-Bayesian least absolute shrinkage and selection operator (BLASSO) was used to estimate all the selected SNP effects for true quantitative trait nucleotide (QTN) detection. This method is referred to as ISIS EM-BLASSO algorithm. Monte Carlo simulation studies validated the new method, which has the highest empirical power in QTN detection and the highest accuracy in QTN effect estimation, and it is the fastest, as compared with efficient mixed-model association (EMMA), smoothly clipped absolute deviation (SCAD), fixed and random model circulating probability unification (FarmCPU), and multi-locus random-SNP-effect mixed linear model (mrMLM). To further demonstrate the new method, six flowering time traits in Arabidopsis thaliana were re-analyzed by four methods (New method, EMMA, FarmCPU, and mrMLM). As a result, the new method identified most previously reported genes. Therefore, the new method is a good alternative for multi-locus GWAS.