Estimating population haplotype frequencies from pooled SNP data using incomplete database information.

MOTIVATION: Information about haplotype structures gives a more detailed picture of genetic variation between individuals than single-locus analyses. Databases that contain the most frequent haplotypes of certain populations are developing rapidly (e.g. the HapMap database for single-nucleotide polymorphisms in humans). Utilization of such prior information about the prevailing haplotype structures makes it possible to estimate the haplotype frequencies also from large DNA pools. When genetic material from dozens of individuals is pooled together and analysed in a single genotyping, the overall number of genotypings and the costs of the genetic studies are reduced.
RESULTS: A Bayesian model for estimating the haplotypes and their frequencies from pooled allelic observations is introduced. The model combines an idea of using database information for haplotype estimation with a computationally efficient multinormal approximation. In addition, the model treats the number and structures of the unknown haplotypes as random variables whose joint posterior distribution is estimated. The results on real human data from the HapMap database show that the proposed method provides significant improvements over the existing methods. AVAILABILITY: A reversible-jump Markov chain Monte Carlo algorithm for analysing the model is implemented in a program called Hippo (Haplotype estimation under incomplete prior information using pooled observations). For comparisons, an approximate expectation-maximization algorithm (EM-algorithm) that utilizes database information about the existing haplotypes is implemented in a program called AEML. The source codes written in C (using GNU Scientific Library) are available at www.iki.fi/ approximately mpirinen. CONTACT: matti.pirinen@iki.fi