Haplotype reconstruction from SNP fragments by minimum error correction.
BACKGROUND
Haplotype reconstruction based on aligned single nucleotide polymorphism (SNP) fragments is to infer a pair of haplotypes from localized polymorphism data gathered through short genome fragment assembly. An important computational model of this problem is the minimum error correction (MEC) model, which has been mentioned in several literatures. The model retrieves a pair of haplotypes by correcting minimum number of SNPs in given genome fragments coming from an individual's DNA.


RESULTS
In the first part of this paper, an exact algorithm for the MEC model is presented. Owing to the NP-hardness of the MEC model, we also design a genetic algorithm (GA). The designed GA is intended to solve large size problems and has very good performance. The strength and weakness of the MEC model are shown using experimental results on real data and simulation data. In the second part of this paper, to improve the MEC model for haplotype reconstruction, a new computational model is proposed, which simultaneously employs genotype information of an individual in the process of SNP correction, and is called MEC with genotype information (shortly, MEC/GI). Computational results on extensive datasets show that the new model has much higher accuracy in haplotype reconstruction than the pure MEC model.
Wang
Rui-Sheng
RS
Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100080, China. wangrsh@amss.ac.cn
Wu
Ling-Yun
LY
Li
Zhen-Ping
ZP
Zhang
Xiang-Sun
XS
D000465:Algorithms; D002874:Chromosome Mapping / Q000379:methods*; D003198:Computer Simulation; D004252:DNA Mutational Analysis / Q000379:methods*; D015894:Genome, Human; D006239:Haplotypes / Q000235:genetics*; D006801:Humans; D008957:Models, Genetic; D015233:Models, Statistical; D010802:Phylogeny; D020641:Polymorphism, Single Nucleotide / Q000235:genetics*; D016415:Sequence Alignment / Q000379:methods*; D017422:Sequence Analysis, DNA / Q000379:methods*