On the minimum number of topologies explaining a sample of DNA sequences.

In this article I derive an alternative algorithm to Hudson and Kaplan's (Genetics 111, 147-165) algorithm that gives a lower bound to the number of recombination events in a sample's history. It is shown that the number, T(M), found by the algorithm is the least number of topologies required to explain a set of DNA sequences sampled under the infinite-site assumption. Let Tao = (T(1),...,T(r)) be a list of topologies compatible with the sequences, i.e., T(k) is compatible with an interval, I(k), of sites in the alignment. A characterization of all lists having T(M) topologies is given and it is shown that T(M) relates to specific patterns in the alignment, here called chain series. Further, a number of theorems relating general lists of topologies to the number T(M) is presented. The results are discussed in relation to the true minimum number of recombination events required to explain an alignment. Copyright 2002 Elsevier Science (USA)