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

Efficient test for deviation from Hardy Weinberg Equilibrium with known or ambiguous typing in highly polymorphic loci

Efficient test for deviation from Hardy–Weinberg equilibrium with known or ambiguous typing in highly polymorphic loci

A Shrinkage Method for Testing the Hardy–Weinberg Equilibrium in Case‐Control Studies

A straightforward multiallelic significance test for the Hardy-Weinberg equilibrium law

An unconditional exact test for the Hardy-Weinberg equilibrium law: sample-space ordering using the Bayes factor.

Population Admixture: Detection by Hardy-Weinberg Test and Its Quantitative Effects on Linkage-Disequilibrium Methods for Localizing Genes Underlying Complex Traits.

Regression-based Approach for Testing the Association Between Multi-Region Haplotype Configuration and Complex Trait

Genuine Bayesian multiallelic significance test for the Hardy-Weinberg equilibrium law

A Fast and Flexible Statistical Model for Large-Scale Population Genotype Data: Applications to Inferring Missing Genotypes and Haplotypic Phase

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

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

Association Test Between Haplotypes and Longitudinal Traits in Complex Pedigrees.

Robust identification of local adaptation from allele frequencies

Demography-adjusted tests of neutrality based on genome-wide SNP data

Statistical Tests for Detecting Positive Selection by Utilizing High-Frequency Variants

Testing Linkage Disequilibrium from Pooled DNA: A Contingency Table Perspective

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

A Probabilistic Method for Testing and Estimating Selection Differences Between Populations.

Efficient variant set mixed model association tests for continuous and binary traits in large-scale whole genome sequencing studies

Compound Tests for the Detection of Hitchhiking Under Positive Selection

A Model-Embedded Trend Test with Incorporating Hardy-Weinberg Equilibrium Information