Haplotype estimation explained

In genetics, haplotype estimation (also known as "phasing") refers to the process of statistical estimation of haplotypes from genotype data. The most common situation arises when genotypes are collected at a set of polymorphic sites from a group of individuals. For example in human genetics, genome-wide association studies collect genotypes in thousands of individuals at between 200,000-5,000,000 SNPs using microarrays. Haplotype estimation methods are used in the analysis of these datasets and allow genotype imputation^{[1] [2]} of alleles from reference databases such as the HapMap Project and the 1000 Genomes Project.

Genotypes and haplotypes

Genotypes measure the unordered combination of alleles at each locus, whereas haplotypes represent the genetic information on multiple loci that have been inherited together from an individual's parents. Theoretically the number of possible haplotypes equals to the product of allele numbers of each locus in consideration. Specially, most of the SNPs are bi-allelic; Therefore when considering

N heterozygous bi-allelic loci, there will be 2 ^{N} possible pairs of haplotypes that could underlie the genotypes. For example, when considering two bi-allelic loci A and B (N =2 ), of which the genotypes are a₁ and a₂, b₁ and b₂, respectively, we will have the following haplotypes: a₁_b₁, a₁_b₂, a₂_b₁, and a₂_b₂ ("_" indicates that the alleles are on the same chromosome).

Haplotype estimation methods

Many statistical methods have been proposed for estimation of haplotypes. Some of the earliest approaches used a simple multinomial model in which each possible haplotype consistent with the sample was given an unknown frequency parameter and these parameters were estimated with an Expectation–maximization algorithm. These approaches were only able to handle small numbers of sites at once, although sequential versions were later developed, specifically the SNPHAP method.

The most accurate and widely used methods for haplotype estimation utilize some form of hidden Markov model (HMM) to carry out inference. For a long time PHASE^[3] was the most accurate method. PHASE was the first method to utilize ideas from coalescent theory concerning the joint distribution of haplotypes. This method used a Gibbs sampling approach in which each individuals haplotypes were updated conditional upon the current estimates of haplotypes from all other samples. Approximations to the distribution of a haplotype conditional upon a set of other haplotypes were used for the conditional distributions of the Gibbs sampler. PHASE was used to estimate the haplotypes from the HapMap Project. PHASE was limited by its speed and was not applicable to datasets from genome-wide association studies.

The fastPHASE ^[4] and BEAGLE methods ^[5] introduced haplotype cluster models applicable to GWAS-sized datasets. Subsequently the IMPUTE2^[6] and MaCH^[7] methods were introduced that were similar to the PHASE approach but much faster. These methods iteratively update the haplotype estimates of each sample conditional upon a subset of K haplotype estimates of other samples. IMPUTE2 introduced the idea of carefully choosing which subset of haplotypes to condition on to improve accuracy. Accuracy increases with K but with quadratic

O(K ^{2)} computational complexity.

The SHAPEIT1 method made a major advance by introducing a linear

O(K) complexity method that operates only on the space of haplotypes consistent with an individual’s genotypes.^[8] The HAPI-UR method subsequently proposed a very similar method.^[9] SHAPEIT2 ^[10] combines the best features of SHAPEIT1 and IMPUTE2 to improve efficiency and accuracy.

See also

• List of haplotype estimation and genotype imputation software
• imputation: predict missing genotypes using known haplotypes

Notes and References

1. Marchini . J. . Howie . B. . 10.1038/nrg2796 . Genotype imputation for genome-wide association studies . Nature Reviews Genetics . 11 . 7 . 499–511 . 2010 . 20517342 . 1465707 .
2. Howie . B. . Fuchsberger . C. . Stephens . M. . Marchini . J. . Abecasis . G. A. R. . Fast and accurate genotype imputation in genome-wide association studies through pre-phasing . 10.1038/ng.2354 . Nature Genetics . 44 . 8 . 955–959 . 2012 . 22820512 . 3696580.
3. Stephens . M. . Smith . N. J. . Donnelly . P. . 10.1086/319501 . A New Statistical Method for Haplotype Reconstruction from Population Data . The American Journal of Human Genetics . 68 . 4 . 978–989 . 2001 . 11254454 . 1275651 .
4. Scheet . P. . Stephens . M. . 10.1086/502802 . A Fast and Flexible Statistical Model for Large-Scale Population Genotype Data: Applications to Inferring Missing Genotypes and Haplotypic Phase . The American Journal of Human Genetics . 78 . 4 . 629–644 . 2006 . 16532393 . 1424677 .
5. Browning . S. R. . Browning . B. L. . 10.1086/521987 . Rapid and Accurate Haplotype Phasing and Missing-Data Inference for Whole-Genome Association Studies by Use of Localized Haplotype Clustering . The American Journal of Human Genetics . 81 . 5 . 1084–1097 . 2007 . 17924348 . 2265661 .
6. Howie . B. N. . Donnelly . P. . Marchini . J. . Schork . Nicholas J . A Flexible and Accurate Genotype Imputation Method for the Next Generation of Genome-Wide Association Studies . 10.1371/journal.pgen.1000529 . PLOS Genetics . 5 . 6 . e1000529 . 2009 . 19543373 . 2689936 . free .
7. Li . Y. . Willer . C. J. . Ding . J. . Scheet . P. . Abecasis . G. A. R. . MaCH: Using sequence and genotype data to estimate haplotypes and unobserved genotypes . 10.1002/gepi.20533 . Genetic Epidemiology . 34 . 8 . 816–834 . 2010 . 21058334 . 3175618 .
8. Delaneau . O. . Marchini . J. . Zagury . J. F. O. . 10.1038/nmeth.1785 . A linear complexity phasing method for thousands of genomes . Nature Methods . 9 . 2 . 179–181 . 2011 . 22138821 . 13765612 .
9. Williams . A. L. . Patterson . N. . Glessner . J. . Hakonarson . H. . Reich . D. . Phasing of Many Thousands of Genotyped Samples . 10.1016/j.ajhg.2012.06.013 . The American Journal of Human Genetics . 91 . 2 . 238–251 . 2012 . 22883141 . 3415548 .
10. Delaneau . O. . Zagury . J. F. . Marchini . J. . 10.1038/nmeth.2307 . Improved whole-chromosome phasing for disease and population genetic studies . Nature Methods . 10 . 1 . 5–6 . 2012 . 23269371 . 205421216 .