Grouping preprocess to accurately extend application of EM algorithm to haplotype inference.
Haplotype inference is an indispensable technique in medical science, especially in genome-wide association studies. Although the conventional method of inference using the expectation-maximization (EM) algorithm by Excoffier and Slatkin is one standard approach, as its calculation cost is an exponential function of the maximum number of heterozygous loci, it has not been widely applied. We propose a method of haplotype inference that can empirically accommodate up to several tens of single nucleotide polymorphism loci in a single haplotype block while maintaining criteria that are exactly equivalent to those of the EM algorithm. The idea is to reduce the cost of calculating the EM algorithm by using a haplotype-grouping preprocess exploiting the symmetrical and inclusive relationships of haplotypes based on the Hardy-Weinberg equilibrium. Testing of the proposed method using real data sets revealed that it has a wider range of applications than the EM algorithm.
Department of Electrical Engineering and Bioscience, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo, 1698555, Japan.
Hiroyuki
H
Shindo
Department of Electrical Engineering and Bioscience, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo, 1698555, Japan.
Hiroshi
H
Chigira
Digital Information Technologies Corporation, daVinci Sakurabashi Bldg., 4-5-4 Hatchobori, Chuo-ku, Tokyo, 1040032, Japan.
Junji
J
Tanaka
Institute of Rheumatology, Tokyo Women's Medical University, 10-22, Kawada-cho, Shinjuku-ku, Tokyo, 1620054, Japan.
Naoyuki
N
Kamatani
Department of Electrical Engineering and Bioscience, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo, 1698555, Japan. masato.inoue@eb.waseda.ac.jp.
Masato
M
Inoue
D000465:Algorithms; D005826:Genetics, Medical / Q000379:methods* / Q000706:statistics & numerical data*; D006239:Haplotypes / Q000235:genetics*; D008957:Models, Genetic; D020641:Polymorphism, Single Nucleotide / Q000235:genetics; D012984:Software