title: Hap-seq: an optimal algorithm for haplotype phasing with imputation using sequencing data.
abstract: Inference of haplotypes, or the sequence of alleles along each chromosome, is a fundamental problem in genetics and is important for many analyses, including admixture mapping, identifying regions of identity by descent, and imputation. Traditionally, haplotypes are inferred from genotype data obtained from microarrays using information on population haplotype frequencies inferred from either a large sample of genotyped individuals or a reference dataset such as the HapMap. Since the availability of large reference datasets, modern approaches for haplotype phasing along these lines are closely related to imputation methods. When applied to data obtained from sequencing studies, a straightforward way to obtain haplotypes is to first infer genotypes from the sequence data and then apply an imputation method. However, this approach does not take into account that alleles on the same sequence read originate from the same chromosome. Haplotype assembly approaches take advantage of this insight and predict haplotypes by assigning the reads to chromosomes in such a way that minimizes the number of conflicts between the reads and the predicted haplotypes. Unfortunately, assembly approaches require very high sequencing coverage and are usually not able to fully reconstruct the haplotypes. In this work, we present a novel approach, Hap-seq, which is simultaneously an imputation and assembly method that combines information from a reference dataset with the information from the reads using a likelihood framework. Our method applies a dynamic programming algorithm to identify the predicted haplotype, which maximizes the joint likelihood of the haplotype with respect to the reference dataset and the haplotype with respect to the observed reads. We show that our method requires only low sequencing coverage and can reconstruct haplotypes containing both common and rare alleles with higher accuracy compared to the state-of-the-art imputation methods.
authors: IBM T.J. Watson Research, Yorktown Heights, NY 10598, USA. dhe@us.ibm.com
authors: Dan
authors: D
authors: He
authors: Buhm
authors: B
authors: Han
authors: Eleazar
authors: E
authors: Eskin
mesh_terms: D000465:Algorithms; D000483:Alleles; D002874:Chromosome Mapping / Q000379:methods / Q000706:statistics & numerical data*; D003198:Computer Simulation; D030541:Databases, Genetic; D005787:Gene Frequency; D014644:Genetic Variation; D060148:HapMap Project; D006239:Haplotypes; D006801:Humans; D016013:Likelihood Functions; D058977:Molecular Sequence Annotation / Q000379:methods / Q000706:statistics & numerical data*
references: Am J Hum Genet. 2004 May;74(5):979-1000
references: 15088269
references: Am J Hum Genet. 2001 Apr;68(4):978-89
references: 11254454
references: Bioinformatics. 2008 Aug 15;24(16):i153-9
references: 18689818
references: Bioinformatics. 2004 Aug 12;20(12):1842-9
references: 14988101
references: Science. 2001 Nov 23;294(5547):1719-23
references: 11721056
references: Nature. 2007 Oct 18;449(7164):851-61
references: 17943122
references: Am J Hum Genet. 2011 Feb 11;88(2):173-82
references: 21310274
references: Mol Biol Evol. 1990 Mar;7(2):111-22
references: 2108305
references: Nat Genet. 2007 Jul;39(7):906-13
references: 17572673
references: J Bioinform Comput Biol. 2003 Apr;1(1):1-20
references: 15290779
references: Nature. 2008 Apr 17;452(7189):872-6
references: 18421352
references: Am J Hum Genet. 2007 Nov;81(5):1084-97
references: 17924348
references: Genet Epidemiol. 2010 Dec;34(8):816-34
references: 21058334
references: PLoS Biol. 2007 Sep 4;5(10):e254
references: 17803354
references: Genome Res. 2008 Aug;18(8):1336-46
references: 18676820
references: Genome Res. 2009 Feb;19(2):318-26
references: 18971310
references: J Comput Biol. 2010 Mar;17(3):547-60
references: 20377463
references: Bioinformatics. 2010 Jun 15;26(12):i183-90
references: 20529904
references: Am J Hum Genet. 2010 Apr 9;86(4):526-39
references: 20303063