Selecting additional tag SNPs for tolerating missing data in genotyping
1 Department of Computer Science and Information Engineering, National Taiwan University, Taipei, Taiwan
2 Graduate Institute of Networking and Multimedia, National Taiwan University, Taipei, Taiwan
3 Section on Statistical Genetics, Department of Biostatistics, University of Alabama at Birmingham, USA
4 Department of Biological Sciences, University of Southern California, Los Angeles, CA 90089, USA
BMC Bioinformatics 2005, 6:263 doi:10.1186/1471-2105-6-263Published: 1 November 2005
Recent studies have shown that the patterns of linkage disequilibrium observed in human populations have a block-like structure, and a small subset of SNPs (called tag SNPs) is sufficient to distinguish each pair of haplotype patterns in the block. In reality, some tag SNPs may be missing, and we may fail to distinguish two distinct haplotypes due to the ambiguity caused by missing data.
We show there exists a subset of SNPs (referred to as robust tag SNPs) which can still distinguish all distinct haplotypes even when some SNPs are missing. The problem of finding minimum robust tag SNPs is shown to be NP-hard. To find robust tag SNPs efficiently, we propose two greedy algorithms and one linear programming relaxation algorithm. The experimental results indicate that (1) the solutions found by these algorithms are quite close to the optimal solution; (2) the genotyping cost saved by using tag SNPs can be as high as 80%; and (3) genotyping additional tag SNPs for tolerating missing data is still cost-effective.
Genotyping robust tag SNPs is more practical than just genotyping the minimum tag SNPs if we can not avoid the occurrence of missing data. Our theoretical analysis and experimental results show that the performance of our algorithms is not only efficient but the solution found is also close to the optimal solution.