title: An EM algorithm based on an internal list for estimating haplotype distributions of rare variants from pooled genotype data.
abstract: BACKGROUND
  Pooling is a cost effective way to collect data for genetic association studies, particularly for rare genetic variants. It is of interest to estimate the haplotype frequencies, which contain more information than single locus statistics. By viewing the pooled genotype data as incomplete data, the expectation-maximization (EM) algorithm is the natural algorithm to use, but it is computationally intensive. A recent proposal to reduce the computational burden is to make use of database information to form a list of frequently occurring haplotypes, and to restrict the haplotypes to come from this list only in implementing the EM algorithm. There is, however, the danger of using an incorrect list, and there may not be enough database information to form a list externally in some applications.
  
  
  RESULTS
  We investigate the possibility of creating an internal list from the data at hand. One way to form such a list is to collapse the observed total minor allele frequencies to "zero" or "at least one", which is shown to have the desirable effect of amplifying the haplotype frequencies. To improve coverage, we propose ways to add and remove haplotypes from the list, and a benchmarking method to determine the frequency threshold for removing haplotypes. Simulation results show that the EM estimates based on a suitably augmented and trimmed collapsed data list (ATCDL) perform satisfactorily. In two scenarios involving 25 and 32 loci respectively, the EM-ATCDL estimates outperform the EM estimates based on other lists as well as the collapsed data maximum likelihood estimates.
  
  
  CONCLUSIONS
  The proposed augmented and trimmed CD list is a useful list for the EM algorithm to base upon in estimating the haplotype distributions of rare variants. It can handle more markers and larger pool size than existing methods, and the resulting EM-ATCDL estimates are more efficient than the EM estimates based on other lists.
authors.affiliation: Department of Statistics and Applied Probability, National University of Singapore, Singapore, 117546, Singapore. lx@nus.edu.sg.
authors.forename: Anthony Y C
authors.initials: AY
authors.lastname: Kuk
authors.forename: Xiang
authors.initials: X
authors.lastname: Li
authors.forename: Jinfeng
authors.initials: J
authors.lastname: Xu
mesh_terms: D000465:Algorithms; D000581:Amidohydrolases / Q000235:genetics; D005787:Gene Frequency; D056426:Genetic Loci; D014644:Genetic Variation; D005838:Genotype; D006239:Haplotypes; D006801:Humans; D008994:Monoacylglycerol Lipases / Q000235:genetics; D009154:Mutation; D009765:Obesity / Q000235:genetics; D000099252:Fatty Acid Amide Hydrolases
chemical_list: D000581:Amidohydrolases; D008994:Monoacylglycerol Lipases; D000099252:Fatty Acid Amide Hydrolases
references.citation: Genet Epidemiol. 2004 Dec;27(4):321-33
references.pmid: 15368617
references.citation: Bioinformatics. 2008 Sep 1;24(17):1896-902
references.pmid: 18617537
references.citation: Bioinformatics. 2010 Oct 15;26(20):2556-63
references.pmid: 20801910
references.citation: Nucleic Acids Res. 2008 Apr;36(6):e35
references.pmid: 18276640
references.citation: Mol Biol Evol. 1995 Sep;12(5):921-7
references.pmid: 7476138
references.citation: Genet Epidemiol. 2002 Oct;23(3):221-33
references.pmid: 12384975
references.citation: Am J Hum Genet. 2006 Mar;78(3):437-50
references.pmid: 16465620
references.citation: Stat Med. 2013 Apr 15;32(8):1343-60
references.pmid: 22855289
references.citation: PLoS Comput Biol. 2010 Oct 14;6(10):e1000954
references.pmid: 20976246
references.citation: Bioinformatics. 2007 Nov 15;23(22):3048-55
references.pmid: 17895275
references.citation: Genet Res (Camb). 2008 Dec;90(6):509-24
references.pmid: 19123969
references.citation: Am J Hum Genet. 2005 Mar;76(3):449-62
references.pmid: 15700229
references.citation: Nat Rev Genet. 2010 Jun;11(6):446-50
references.pmid: 20479774
references.citation: Bioinformatics. 2008 Sep 1;24(17):1942-8
references.pmid: 18573795
references.citation: Ann Med. 2004;36(2):146-52
references.pmid: 15119834
references.citation: Annu Rev Genomics Hum Genet. 2008;9:387-402
references.pmid: 18576944
references.citation: BMC Genet. 2012 Oct 30;13:94
references.pmid: 23110720
references.citation: Bioinformatics. 2009 Dec 15;25(24):3296-302
references.pmid: 19861356
references.citation: Genet Epidemiol. 2012 Dec;36(8):870-81
references.pmid: 22972696
references.citation: Nat Rev Genet. 2002 Nov;3(11):862-71
references.pmid: 12415316
references.citation: Nat Rev Genet. 2011 Feb;12(2):77
references.pmid: 21191422
references.citation: Am J Hum Genet. 2011 Sep 9;89(3):382-97
references.pmid: 21855840
references.citation: Genet Epidemiol. 2004 Dec;27(4):334-47
references.pmid: 15368348
references.citation: Am J Hum Genet. 2003 Feb;72(2):384-98
references.pmid: 12533787
references.citation: Genet Epidemiol. 2004 Dec;27(4):348-64
references.pmid: 15543638
references.citation: Bioinformatics. 2009 Feb 1;25(3):379-86
references.pmid: 19050036
references.citation: Genet Epidemiol. 2010 Jul;34(5):479-91
references.pmid: 20552648
references.citation: IEEE/ACM Trans Comput Biol Bioinform. 2011 Jan-Mar;8(1):36-44
references.pmid: 21071795
references.citation: Nat Rev Genet. 2011 Mar;12(3):215-23
references.pmid: 21301473
references.citation: Nucleic Acids Res. 2006 Feb 14;34(4):e27
references.pmid: 16478714