Title
Leveraging reads that span multiple single nucleotide polymorphisms for haplotype inference from sequencing data.
Abstract
Motivation: Haplotypes, defined as the sequence of alleles on one chromosome, are crucial for many genetic analyses. As experimental determination of haplotypes is extremely expensive, haplotypes are traditionally inferred using computational approaches from genotype data, i.e. the mixture of the genetic information from both haplotypes. Best performing approaches for haplotype inference rely on Hidden Markov Models, with the underlying assumption that the haplotypes of a given individual can be represented as a mosaic of segments from other haplotypes in the same population. Such algorithms use this model to predict the most likely haplotypes that explain the observed genotype data conditional on reference panel of haplotypes. With rapid advances in short read sequencing technologies, sequencing is quickly establishing as a powerful approach for collecting genetic variation information. As opposed to traditional genotyping-array technologies that independently call genotypes at polymorphic sites, short read sequencing often collects haplotypic information; a read spanning more than one polymorphic locus (multi-single nucleotide polymorphic read) contains information on the haplotype from which the read originates. However, this information is generally ignored in existing approaches for haplotype phasing and genotype-calling from short read data. Results: In this article, we propose a novel framework for haplotype inference from short read sequencing that leverages multi-single nucleotide polymorphic reads together with a reference panel of haplotypes. The basis of our approach is a new probabilistic model that finds the most likely haplotype segments from the reference panel to explain the short read sequencing data for a given individual. We devised an efficient sampling method within a probabilistic model to achieve superior performance than existing methods. Using simulated sequencing reads from real individual genotypes in the HapMap data and the 1000 Genomes projects, we show that our method is highly accurate and computationally efficient. Our haplotype predictions improve accuracy over the basic haplotype copying model by similar to 20% with comparable computational time, and over another recently proposed approach Hap-SeqX by similar to 10% with significantly reduced computational time and memory usage.
Year
DOI
Venue
2013
10.1093/bioinformatics/btt386
BIOINFORMATICS
Field
DocType
Volume
Population,International HapMap Project,Computer science,Haplotype,1000 Genomes Project,Statistical model,Bioinformatics,Locus (genetics),Hidden Markov model,Haplotype estimation
Journal
29
Issue
ISSN
Citations 
18
1367-4803
6
PageRank 
References 
Authors
0.59
8
6
Name
Order
Citations
PageRank
Wen-Yun Yang1417.56
Farhad Hormozdiari211611.73
Zhanyong Wang3507.04
Dan He413312.54
Bogdan Paşaniuc59515.06
Eleazar Eskin61790170.53