An empirical Bayes method for genotyping and SNP detection using multi-sample next-generation sequencing data

MOTIVATION The development of next generation sequencing technology provides an efficient and powerful approach to rare variant detection. To identify genetic variations, the essential question is how to quantity the sequencing error rate in the data. Because of the advantage of easy implementation and the ability to integrate data from different sources, the empirical Bayes method is popularly employed to estimate the sequencing error rate for SNP detection. RESULTS We propose a novel statistical model to fit the observed non-reference allele frequency data, and utilize the empirical Bayes method for both genotyping and SNP detection, where an ECM algorithm is implemented to estimate the model parameters. The performance of our proposed method is investigated via simulations and real data analysis. It is shown that our method makes less genotype-call errors, and with the parameter estimates from the ECM algorithm, it attains high detection power with FDR being well controlled. AVAILABILITY AND IMPLEMENTATION The proposed algorithm is wrapped in the R package ebGenotyping, which can be downloaded from http://cran.r-project.org/web/packages/ebGenotyping/ CONTACT: youn@mail.sysu.edu.cnSupplementary information: Supplementary data are available at Bioinformatics online.

[1]  Si Quang Le,et al.  SNP detection and genotyping from low-coverage sequencing data on multiple diploid samples. , 2011, Genome research.

[2]  Xiao-Li Meng,et al.  Maximum likelihood estimation via the ECM algorithm: A general framework , 1993 .

[3]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[4]  Mark Gerstein,et al.  Personal genome sequencing: current approaches and challenges. , 2010, Genes & development.

[5]  Michael A. Schmidt,et al.  SeqEM: an adaptive genotype-calling approach for next-generation sequencing studies , 2010, Bioinform..

[6]  S. Schuster.  Next-generation sequencing transforms today's biology , 2008, Nature Methods.

[7]  John G. Cleary,et al.  Joint Variant and De Novo Mutation Identification on Pedigrees from High-Throughput Sequencing Data , 2014, bioRxiv.

[8]  Zhi Wei,et al.  An empirical Bayes testing procedure for detecting variants in analysis of next generation sequencing data , 2013, 1401.2278.

[9]  Suzanne M. Leal,et al.  Discovery of Rare Variants via Sequencing: Implications for the Design of Complex Trait Association Studies , 2009, PLoS genetics.

[10]  Bradley Efron,et al.  Large-scale inference , 2010 .

[11]  Kevin P. Murphy,et al.  SNVMix: predicting single nucleotide variants from next-generation sequencing of tumors , 2010, Bioinform..

[12]  B. Efron.  Size, power and false discovery rates , 2007, 0710.2245.

[13]  Nancy R. Zhang,et al.  A cross-sample statistical model for SNP detection in short-read sequencing data , 2011, Nucleic acids research.

[14]  M. DePristo,et al.  A framework for variation discovery and genotyping using next-generation DNA sequencing data , 2011, Nature Genetics.

[15]  Huanming Yang,et al.  SNP detection for massively parallel whole-genome resequencing. , 2009, Genome research.

[16]  Nancy R. Zhang,et al.  Detecting Mutations in Mixed Sample Sequencing Data Using Empirical Bayes , 2012, 1209.6453.

[17]  Wei Cui,et al.  MultiGeMS: detection of SNVs from multiple samples using model selection on high-throughput sequencing data , 2016, Bioinform..

[18]  B. Efron.  Large-Scale Simultaneous Hypothesis Testing , 2004 .

[19]  Jian Xu,et al.  SNP calling using genotype model selection on high-throughput sequencing data , 2012, Bioinform..

[20]  Baiyu Zhou,et al.  An empirical Bayes mixture model for SNP detection in pooled sequencing data , 2012, Bioinform..

[21]  B. Efron.  Correlation and Large-Scale Simultaneous Significance Testing , 2007 .

[22]  G. Weinstock,et al.  A SNP discovery method to assess variant allele probability from next-generation resequencing data. , 2010, Genome research.