Abstract:MOTIVATION:Pooling large number of DNA samples is a common practice in association study, especially for initial screening. However, the use of expectation-maximization (EM)-type algorithms in estimating haplotype distributions for even moderate pool sizes is hampered by the computational complexity involved. A novel constrained EM algorithm called PoooL has been proposed recently to bypass the difficulty via the use of asymptotic normality of the pooled allele frequencies. The resulting estimates are, however, not maximum likelihood estimates and hence not optimal. Furthermore, the assumption of Hardy-Weinberg equilibrium (HWE) made may not be realistic in practice.METHODS:Rather than carrying out constrained maximization as in PoooL, we revert to the usual EM algorithm but make it computationally feasible by using normal approximations. The resulting algorithm is much simpler to implement than PoooL because there is no need to invoke sophisticated iterative scaling methods as in PoooL. We also develop an estimating equation analogue of the EM algorithm for the case of Hardy-Weinberg disequilibrium (HWD) by conditioning on the haplotypes of both chromosomes of the same individual. Incorporated into the method is a way of estimating the inbreeding coefficient by relating it to overdispersion.RESULTS:Simulation study assuming HWE shows that our simplified implementation of the EM algorithm leads to estimates with substantially smaller SDs than PoooL estimates. Further simulations show that ignoring HWD will induce biases in the estimates. Our extended method with estimation of inbreeding coefficient incorporated is able to reduce the bias leading to estimates with substantially smaller mean square errors. We also present results to suggest that our method can cope with a certain degree of locus-specific inbreeding as well as additional overdispersion not caused by inbreeding.AVAILABILITY:http://staff.ustc.edu.cn/ approximately ynyang/aem-aes

An EM method based on entropy LD block partition for haplotype inference

Accurate Haplotype Inference for Multiple Linked Single-Nucleotide Polymorphisms Using Sibship Data

Haplotype Block Partitioning and Tag SNP Selection Using Genotype Data and Their Applications to Association Studies

A faster haplotyping algorithm based on block partition, and greedy ligation strategy

A Better Block Partition and Ligation Strategy for Individual Haplotyping

Nonparametric Disequilibrium Mapping of Functional Sites Using Haplotypes of Multiple Tightly Linked Single-Nucleotide Polymorphism Markers

The Effect of Haplotype-Block Definitions on Inference of Haplotype-Block Structure and Htsnps Selection

Computationally Feasible Estimation of Haplotype Frequencies from Pooled Dna with and Without Hardy-Weinberg Equilibrium

A Multilocus Linkage Disequilibrium Measure Based on Mutual Information Theory and Its Applications

Fast and accurate haplotype inference with hidden markov model

Inference of missing SNPs and information quantity measurements for haplotype blocks.

A New Method for Haplotype Inference Including Full-Sib Information

Patterns of Linkage Disequilibrium and Haplotype Distribution in Disease Candidate Genes

Dynamic Programming Algorithms for Haplotype Block Partitioning and Tag SNP Selection Using Haplotype Data or Genotype Data

A Dynamic Programming Algorithm for Haplotype Block Partitioning and Its Application in Association Studies.

A Fast Haplotype Inference Method for Large Population Genotype Data

Haplotypes Vs Single Marker Linkage Disequilibrium Tests: What Do We Gain?

Comparison of haplotype inference methods using genotypic data from unrelated individuals.

Haplotype Reconstruction In Half-Sib Pedigree Using A Combination Of EM Algorithm And Rule-Based Methods

Accurate and Efficient Estimation of Local Heritability using Summary Statistics and LD Matrix

THE HAPLOTYPE LINKAGE DISEQUILIBRIUM TEST FOR GENOME-WIDE SCREENS: ITS POWER AND STUDY DESIGN