# A double classification tree search algorithm for index SNP selection

- Peisen Zhang
^{1}Email author, - Huitao Sheng
^{2}and - Ryuhei Uehara
^{3}

**5**:89

**DOI: **10.1186/1471-2105-5-89

© Zhang et al; licensee BioMed Central Ltd. 2004

**Received: **22 January 2004

**Accepted: **06 July 2004

**Published: **06 July 2004

## Abstract

### Background

In population-based studies, it is generally recognized that single nucleotide polymorphism (SNP) markers are not independent. Rather, they are carried by haplotypes, groups of SNPs that tend to be coinherited. It is thus possible to choose a much smaller number of SNPs to use as indices for identifying haplotypes or haplotype blocks in genetic association studies. We refer to these characteristic SNPs as index SNPs. In order to reduce costs and work, a minimum number of index SNPs that can distinguish all SNP and haplotype patterns should be chosen. Unfortunately, this is an NP-complete problem, requiring brute force algorithms that are not feasible for large data sets.

### Results

We have developed a double classification tree search algorithm to generate index SNPs that can distinguish all SNP and haplotype patterns. This algorithm runs very rapidly and generates very good, though not necessarily minimum, sets of index SNPs, as is to be expected for such NP-complete problems.

### Conclusions

A new algorithm for index SNP selection has been developed. A webserver for index SNP selection is available at

## Background

Because SNPs are often coinherited as components of a haplotype, they can be highly correlated. Because of this, it is theoretically possible to choose a much smaller number of SNPs to be used as an index set in identifying haplotype or SNP patterns. Johnson and his collaborators [1] have referred to such characteristic SNPs as haplotype tagging SNPs (htSNPs). Bafna et al. [2] refer to them as informative SNPs, using the language of probability theory. We prefer the use of the more general "index SNPs" to indicate not only haplotype but any SNP patterns. The use of index SNPs can reduce the work in SNP-based genotyping research. Clayton [3] provides computer software for htSNP selection. In his program, he uses five as the default maximum htSNP number and implements a brute force search algorithm to browse over subsets of SNP numbers up to a given maximum, choosing the subset according to predetermined criteria. However, if a large number of index SNPs is required, this algorithm fails. Similarly, Sebastiani and his collaborators [4] have developed a program called BEST (Best Enumeration of SNP Tags); again, use of this program is not feasible with very large sets of SNPs. In the HapScope project, Zhang et al. [5] have developed two programs for selection of index SNPs: BFA, a brute force algorithm and GPA, a greedy partition algorithm. We have re-formulated the index SNP selection problem and developed a new greedy algorithm for index SNP selection based on a double classification tree search algorithm similar to the double search algorithm we previously developed for physical mapping [6]. This is not an enumeration algorithm. It runs rapidly and generates very reasonable results, though not guaranteeing generation of a minimum set, as is expected by the NP-complete nature of problem. The NP-complete property has been proved by Bafna et al [2]. For the reader's convenience, we have attached a brief proof as an appendix.

### Algorithm

#### Classification tree search algorithm for SNP generation

A data sample to show the algorithm.

Haplotype1 | ACAGATG |
---|---|

Haplotype2 | ACGAATG |

Haplotype3 | ATGGGTG |

Haplotype4 | GTAAGTG |

Haplotype5 | GTGGGCA |

Haplotype6 | GTAGACA |

Haplotype7 | ATAAGCA |

Haplotype8 | GTGGACA |

*i.e.*, the tree has a minimum height of three. It is easy to appreciate that there is no classification tree for the above haplotype set with a height less then three. We propose here a greedy algorithm to generate a classification tree with a "good" height, but no guarantee that it is the minimum height. Our algorithm can be divided into two phases: a greedy phase to choose the classifiers and a tree-building phase to divide the haplotype patterns into the subtrees. A classification tree will be built by recurrently switching from greedy phase to tree building phase until all leaves of the tree have only one haplotype pattern. It is the purpose of our greedy method to choose a classifier from among the SNPs based on its possessing the smallest maximum sized subtree compared to those of the SNPs that have not yet been used as classifiers. If more than one SNP generates smallest maximum subtrees of the same height, we then examine the second maximum subtrees. If they are also the same size, we check the third, and so on. If all classifiers have smallest maximum subtrees of the same size, we can choose any one of them. In the above example, the first SNP has 4 as the maximum size of its subtrees. In contrast, the second SNP has 6 as the maximum size of its subtrees, so it would be rejected. The algorithm is described in Figure 2.

This algorithm runs very fast. Let the number of SNPs be N and the number of haplotype patterns be M. The major calculation is on the loop of step 2 and step 3. Since the loop can run no more than the number of SNPs: N, step 2 needs less than **O**(NM) operations. Step 3 needs less than **O**(NM) operations also. The total complexity of this algorithm is below the order of **O**(N^{2}M).

#### Properties of classifiers

**complete set**if and only if the set of classifiers can distinguish haplotypes. If no proper subset of a complete set is a complete set, we will call it

**minimal complete set**. The smallest minimal complete set will be called the

**minimum complete set**.

- (1)
The whole SNP set for the group of haplotypes is a complete set.

- (2)
For SNPs with only two variations (the major and the minor), the size of a complete set of classifiers cannot be less than

**log**_{ 2 }**N**where**N**is the number of haplotypes. - (3)
Any complete set of classifiers can be used to build up a classification tree. If the complete set is a minimal set, the height of the tree is equal to the number of classifiers in the set.

- (4)
The classification tree algorithm generates a complete set.

#### A double classification tree search algorithm

A data sample to show the second round search is needed.

1 | 2 | 3 | 4 | 5 | |
---|---|---|---|---|---|

A | 1 | 1 | 1 | 1 | 1 |

B | 1 | 1 | 1 | 1 | 0 |

C | 1 | 1 | 1 | 0 | 1 |

D | 0 | 1 | 1 | 0 | 0 |

E | 1 | 1 | 0 | 1 | 1 |

F | 1 | 1 | 0 | 1 | 0 |

G | 0 | 1 | 0 | 0 | 1 |

H | 0 | 1 | 0 | 0 | 0 |

I | 1 | 0 | 1 | 1 | 1 |

J | 0 | 0 | 1 | 1 | 0 |

K | 0 | 0 | 1 | 0 | 1 |

L | 0 | 0 | 0 | 1 | 1 |

Using the previously described classification tree search algorithm:

SNP 1 splits the 12 patterns into groups of 6 and 6;

SNP 2 splits the 12 patterns into groups of 8 and 4;

SNP 3 splits the 12 patterns into groups of 7 and 5;

SNP 4 splits the 12 patterns into groups of 7 and 5;

SNP 5 splits the 12 patterns into groups of 7 and 5.

Based on the algorithm, we choose SNP 1 first as classifier. But no set of four SNPs containing SNP 1 suffice to distinguish all 12 patterns:

SNPs 1, 2, 3, 4 cannot distinguish A from B;

SNPs 1, 2, 3, 5 cannot distinguish A from C;

SNPs 1, 2, 4, 5 cannot distinguish A from E;

SNPs 1, 3, 4, 5 cannot distinguish A from I.

Hence, the algorithm will have to choose all five SNPs to distinguish all the patterns. But SNPs 2, 3, 4, 5 will distinguish these patterns, and clearly that is a minimal set. We have been trapped by SNP 1. In order to avoid such a trap, a second round tree search is needed. For the second round search, we force the last classifier of the first round to be used as the first classifier in the second round. The same rule is followed for choosing the second classifier, and so on. By the double search algorithm, in the first search we may generate classifiers in the order: SNP1, SNP5, SNP3, SNP4, and SNP2; in the second search we will generate in order: SNP2, SNP3, SNP4, and SNP5.

#### Index SNP selection with constraints

Sometimes it is necessary to select some important and interesting SNPs as the index SNPs. In that case, we can use those SNPs as classifiers first in building up the trees. Then the greedy algorithm is used to choose additional classifiers. On our webserver, the user can provide a list for those SNPs that definitely should be included.

## Discussion

The index SNP selection problem is a very important and practical problem. Since it is an NP-complete problem, there is no polynomial algorithm so far for an exact solution. Brute force algorithms have been developed that are useful for small sets of data. In contrast, the double search algorithm is good for both small and large data sets. This algorithm gives a quite reasonable solution but is not guaranteed to generate the minimum index set. Given the NP-complete nature of the problem, it may be possible to develop different approximation algorithms in the future.

This program is designed for haplotype data. It can be extended for genotype data. It is our strategy to select a minimum set of index SNPs after a small set of data has been genotyped and haplotypes have been generated. Then the selected minimum index SNPs will be used to genotype the whole sample set.

This program is limited to deal with biallelic SNP. The non-biallelic case and the missing data case can be developed using a SNP pattern extension.

## Appendix: A brief proof of the N-P completeness

We reduce the following NP-complete problem known as the minimum test set problem [7] to the minimum index SNP set problem:

### Input

Collection C of subsets of a finite set S, positive integer k ≤ C.

### Question

Is there a subcollection C'⊆ C with C'≤ k such that for each pair of distinct elements u, v ∈ S, there is some set **c** ∈ C' that contains exactly one of u and v?

**c**

_{ 1 }, ...,

**c**

_{ n }} and S= {

**s**

_{ 1 }, ...,

**s**

_{ m }}. We then construct a set of SNPs as follows;

- (1)
the number of SNPs is

*n*(the number of the size of C), - (2)
the number of SNP patterns is

*m*(the number of the size of S), - (3)
the

*i*th letter of the*j*th SNP pattern is '1' if s_{ j }∈ c_{ i }, otherwise '0'.

Intuitively, the *j* th SNP pattern describes if the element s_{
j
}∈ S is in each subset c_{
i
}or not.

The reduction can be done in linear time, and the solution of the minimum index SNP set problem directly gives the solution of the minimum test set problem.

## Declarations

## Authors’ Affiliations

## References

- Johnson GCL, Esposito L, Barratt BJ, Smith AN, Heward J, Genova GD, Ueda H, Cordell HJ, Eaves IA, Dudbridge F, Twells RCJ, Payne F, Hughes W, Nutland S, Stevens H, Carr P, Tuomilehto-Wolf E, Tuomilehto J, Gough SCL, Clayton DG, Todd JA:
**Haplotype tagging for the identification of common disease genes.***Nature Genet*2001,**29:**233–237. 10.1038/ng1001-233View ArticlePubMedGoogle Scholar - Bafna V, Halldorsson BV, Schwartz R, Clark AG, Istrail S:
**Haplotypes and informative SNP selection algorithms:don't block out information.***RECOMB Berlin, Germany*2003.Google Scholar - Clayton DG:
**Choosing a set of haplotype tagging SNPs from a larger set of diallelic loci.**[http://www.nature.com/ng/journal/v29/n2/extref/ng1001–233-S10.pdf] - Sebastiani P, Lazarus R, Weiss ST, Kunkel LM, Kohane IS, Ramoni MF:
**Minimal haplotype tagging.***Proc Natl Acad Sci USA*2003,**100:**9900–9905. 10.1073/pnas.1633613100PubMed CentralView ArticlePubMedGoogle Scholar - Zhang J, Rowe WL, Struewing JP, Buetow KH:
**HapScope: A Software System for Automated and Visual Analysis of Functionally Annotated Haplotypes.***Nucleic Acids Research*2002,**30:**5213–5221. 10.1093/nar/gkf654PubMed CentralView ArticlePubMedGoogle Scholar - Zhang P, Schon EA, Fischer SG, Cayanis E, Weiss J, Kistler S, Bourne PE:
**An algorithm based on graph theory for the assembly of contigs in physical mapping of DNA.***CABIOS*1994,**10:**309–317.PubMedGoogle Scholar - Garey MR, Johnson DS:
**Computers and Intractability – A Guide to the Theory of NP-Completeness.***Freeman*1979.Google Scholar **UW-FHCRC Variation Discovery Resource (SeattleSNPs)**[http://pga.gs.washington.edu]

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.