Bioinformatics Advance Access originally published online on June 23, 2008
Bioinformatics 2008 24(17):1942-1948; doi:10.1093/bioinformatics/btn324
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
PoooL: an efficient method for estimating haplotype frequencies from large DNA pools
1Department of Statistics and Finance, University of Science and Technology of China, Anhui 230026 and 2Institute of Statistical Science, Academia Sinica, Taipei 11529, Taiwan
*To whom correspondence should be addressed.
 |
Abstract |
|---|
Motivation: Pooling DNA is a cost-effective alternative to individual genotyping method. It is often used for initial screening in genome-wide association analysis. In some studies, large pools with sizes up to several hundreds were applied in order to significantly reduce genotyping cost. However, method for estimating haplotype frequencies from large DNA pools has not been available due to computational complexity involved.
Methods: We propose a novel constrained EM algorithm, PoooL, to estimate frequencies of single-nucleotide polymorphism (SNP) haplotypes from DNA pools. A quantity called importance factor is introduced to measure the contribution of a haplotype to the likelihood. Under the assumption of asymptotic normality of the estimated allele frequencies and a system of linear constraints on haplotype frequencies the importance factor remains a constant in the iterative maximization process. The maximization problem in the EM algorithm is then formulated into a constrained maximum entropy model and solved by the improved iterative scaling method.
Results: Simulation study shows that our algorithm can efficiently estimate haplotype frequencies from DNA pools with arbitrarily large sizes. The algorithm works equally well for large pools with sizes up to hundreds or thousands and for pools with sizes as small as one or two individuals. The computational complexity of the PoooL algorithm is independent of pool sizes, and the computational efficiency for large pools is thus substantially improved over existing estimating methods. Simulation results also show that the proposed method is robust to genotype errors and population admixture.
Availability: http://staff.ustc.edu.cn/~ynyang/poool
Contact: zhanghan{at}mail.ustc.edu.cn; ynyang{at}ustc.edu.cn
Associate Editor: Limsoon Wong
Received on January 19, 2008; revised on June 20, 2008; accepted on June 20, 2008
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  J. E. Craig, A. W. Hewitt, A. E. McMellon, A. K. Henders, L. Ma, L. Wallace, S. Sharma, K. P. Burdon, P. M. Visscher, G. W. Montgomery, et al. Rapid inexpensive genome-wide association using pooled whole blood Genome Res., November 1, 2009; 19(11): 2075 - 2080. [Abstract] [Full Text] [PDF]
|
|