title: BAYSIC: a Bayesian method for combining sets of genome variants with improved specificity and sensitivity.
abstract: BACKGROUND
  Accurate genomic variant detection is an essential step in gleaning medically useful information from genome data. However, low concordance among variant-calling methods reduces confidence in the clinical validity of whole genome and exome sequence data, and confounds downstream analysis for applications in genome medicine.Here we describe BAYSIC (BAYeSian Integrated Caller), which combines SNP variant calls produced by different methods (e.g. GATK, FreeBayes, Atlas, SamTools, etc.) into a more accurate set of variant calls. BAYSIC differs from majority voting, consensus or other ad hoc intersection-based schemes for combining sets of genome variant calls. Unlike other classification methods, the underlying BAYSIC model does not require training using a "gold standard" of true positives. Rather, with each new dataset, BAYSIC performs an unsupervised, fully Bayesian latent class analysis to estimate false positive and false negative error rates for each input method. The user specifies a posterior probability threshold according to the user's tolerance for false positive and false negative errors; lowering the posterior probability threshold allows the user to trade specificity for sensitivity while raising the threshold increases specificity in exchange for sensitivity.
  
  
  RESULTS
  We assessed the performance of BAYSIC in comparison to other variant detection methods using ten low coverage (~5X) samples from The 1000 Genomes Project, a tumor/normal exome pair (40X), and exome sequences (40X) from positive control samples previously identified to contain clinically relevant SNPs. We demonstrated BAYSIC's superior variant-calling accuracy, both for somatic mutation detection and germline variant detection.
  
  
  CONCLUSIONS
  BAYSIC provides a method for combining sets of SNP variant calls produced by different variant calling programs. The integrated set of SNP variant calls produced by BAYSIC improves the sensitivity and specificity of the variant calls used as input. In addition to combining sets of germline variants, BAYSIC can also be used to combine sets of somatic mutations detected in the context of tumor/normal sequencing experiments.
authors.forename: Brandi L
authors.initials: BL
authors.lastname: Cantarel
authors.forename: Daniel
authors.initials: D
authors.lastname: Weaver
authors.forename: Nathan
authors.initials: N
authors.lastname: McNeill
authors.forename: Jianhua
authors.initials: J
authors.lastname: Zhang
authors.forename: Aaron J
authors.initials: AJ
authors.lastname: Mackey
authors.affiliation: Genformatic, LLC, Austin, TX 78731, USA. jreese@genformatic.com.
authors.forename: Justin
authors.initials: J
authors.lastname: Reese
mesh_terms: D000465:Algorithms; D001499:Bayes Theorem; D059472:Exome; D015894:Genome, Human; D006801:Humans; D009154:Mutation; D020641:Polymorphism, Single Nucleotide; D011336:Probability; D012985:Software Design
references.citation: BMC Genomics. 2012 May 20;13:194
references.pmid: 22607156
references.citation: PLoS One. 2007 Apr 18;2(4):e383
references.pmid: 17440619
references.citation: Nat Biotechnol. 2013 Mar;31(3):213-9
references.pmid: 23396013
references.citation: Genome Res. 2010 Sep;20(9):1297-303
references.pmid: 20644199
references.citation: Nat Biotechnol. 2014 Mar;32(3):246-51
references.pmid: 24531798
references.citation: Genome Med. 2013 Mar 27;5(3):28
references.pmid: 23537139
references.citation: Bioinformatics. 2011 Aug 1;27(15):2156-8
references.pmid: 21653522
references.citation: Bioinformatics. 2013 Sep 15;29(18):2223-30
references.pmid: 23842810
references.citation: Hum Mol Genet. 2009 Oct 15;18(R2):R163-8
references.pmid: 19808792
references.citation: Bioinformatics. 2012 Jul 15;28(14):1811-7
references.pmid: 22581179
references.citation: Genome Biol. 2007;8(1):R13
references.pmid: 17241472
references.citation: Bioinformatics. 2013 Jun 15;29(12):1498-503
references.pmid: 23620360
references.citation: Genome Res. 2003 Oct;13(10):2271-6
references.pmid: 14525928
references.citation: Nature. 2012 Nov 1;491(7422):56-65
references.pmid: 23128226
references.citation: Nucleic Acids Res. 2011 Jan;39(Database issue):D945-50
references.pmid: 20952405
references.citation: Nat Genet. 2011 May;43(5):491-8
references.pmid: 21478889
references.citation: Bioinformatics. 2014 May 1;30(9):1198-204
references.pmid: 24443148
references.citation: Genome Med. 2013 Jun 27;5(6):57
references.pmid: 23806086
references.citation: Genome Biol. 2011 Jul 25;12(7):R68
references.pmid: 21787409
references.citation: Am J Hum Genet. 2002 Jun;70(6):1490-7
references.pmid: 11992255
references.citation: Genome Res. 2012 Mar;22(3):568-76
references.pmid: 22300766
references.citation: Bioinformatics. 2009 Aug 15;25(16):2078-9
references.pmid: 19505943
references.citation: BMC Bioinformatics. 2012 Jan 12;13:8
references.pmid: 22239737
references.citation: Bioinformatics. 2013 Sep 1;29(17):2208-10
references.pmid: 23803469
references.citation: Neurology. 2004 May 11;62(9):1598-600
references.pmid: 15136689
references.citation: JAMA. 2014 Mar 12;311(10):1035-45
references.pmid: 24618965