370 likes | 473 Views
A Fundamental Decomposition Theory for Phylogenetic Networks and Incompatible Characters. D. Gusfield, V. Bansal (Recomb 2005). Alternative Title. The Continuing Role of Incompatibility Graphs in the Study of Phylogenetic Networks. Geneological or Phylogenetic Networks.
E N D
A Fundamental Decomposition Theory for Phylogenetic Networks and Incompatible Characters D. Gusfield, V. Bansal (Recomb 2005)
Alternative Title The Continuing Roleof Incompatibility Graphs in the Study of Phylogenetic Networks
Geneological or Phylogenetic Networks • The major biological motivation comes from genetics and attempts to reconstruct the history of recombination in populations. • The results also have phylogenetic applications, for example in hybrid speciation, lateral gene transfer.
Reconstructing the Evolution of Binary Bio-Sequences (SNPs) • Perfect Phylogeny (tree) model • Phylogenetic Networks (DAG) with recombination (ARG) • Blobbed Trees • Incompatibility Graph and Connected its Components • Prior uses of Connected Components • Decomposition Theorem and Proof Sketch • Optimality Conjecture and Progress
The Perfect Phylogeny Model for binary sequences sites 12345 Ancestral sequence 00000 one mutation per site 1 4 Site mutations on edges 3 00010 The tree derives the set M: 10100 10000 01011 01010 00010 2 10100 5 10000 01010 01011 Extant sequences at the leaves
When can a set of sequences be derived on a perfect phylogeny? Classic NASC: Arrange the sequences in a matrix. Then (with no duplicate columns), the sequences can be generated on a unique perfect phylogeny if and only if no two columns (sites) contain all four pairs: 0,0 and 0,1 and 1,0 and 1,1 This is the 4-Gamete Test
A richer model 10100 10000 01011 01010 00010 10101 added 12345 00000 1 4 3 00010 2 10100 5 Pair 4, 5 fails the four gamete-test. The sites 4, 5 are ``incompatible” 10000 01010 01011 Real sequence histories often involve recombination.
Sequence Recombination 01011 10100 S P 5 Single crossover recombination 10101 A recombination of P and S at recombination point 5. The first 4 sites come from P (Prefix) and the sites from 5 onward come from S (Suffix). Called ``crossing over” in genetics
Network with Recombination 10100 10000 01011 01010 00010 10101 new 12345 00000 1 4 3 00010 2 10100 5 P 10000 01010 The previous tree with one recombination event now derives all the sequences. 01011 5 S 10101
Multiple Crossover Recombination 4-crossovers 2-crossovers = ``gene conversion”
A Phylogenetic Network 00000 4 00010 a:00010 3 1 10010 00100 5 00101 2 01100 S b:10010 S P 3 4 01101 p g:00101 c:00100 10100 d:10100 f:01101 e:01100
Minimizing Recombinations • Any set M of sequences can be generated by a phylogenetic network with enough recombinations, and one mutation per site. This is not interesting or useful. • However, the number of (observable) recombinations is small in realistic sets of sequences. ``Observable” depends on n and m relative to the number of recombinations. • Problem: Given a set of sequences M, find a phylogenetic network generating M, minimizing the number of recombinations. NP-hard (Wang et al 2000, Semple et al 2004)
Decomposition can help First we introduce the viewpoint needed.
Blobs in Networks • In a Phylogenetic Network, with a recombination node x, if we trace two paths backwards from x, then the paths will eventually meet. • The cycle specified by those two paths is called a ``recombination cycle”. • In a phylogenetic Network a maximal set of (edge) intersecting cycles is called a blob.
A Phylogenetic Network with one Blob 00000 4 00010 a:00010 3 1 10010 00100 5 00101 2 01100 S b:10010 P S 4 01101 c:00100 p g:00101 3 d:10100 f:01101 e:01100
Blobbed-trees • Contracting each blob to a single node results in a directed, rooted tree, otherwise one of the “blobs” was not maximal. • So every phylogenetic network can be viewed as a directed tree of blobs - a blobbed-tree. The blobs are the non-tree-like parts of the network.
Every network is a tree of blobs. How do the tree parts and the blobs relate? How can we exploit this relationship? Ugly tangled network inside the blob.
1 2 3 4 5 Incompatibility Graph G(M) a b c d e f g 0 0 0 1 0 1 0 0 1 0 0 0 1 0 0 1 0 1 0 0 0 1 1 0 0 0 1 1 0 1 0 0 1 0 1 4 M 1 3 2 5 Two nodes are connected iff the pair of sites are incompatible, i.e, fail the 4-gamete test. G(M) has two connected components.
The connected components of G(M) are very informative • The number of non-trivial connected components is a lower-bound on the number of recombinations needed in any network (Bafna, Bansal; Gusfield, Hickerson). • When each blob is a single-cycle (galled-tree case) all the incompatible sites in a blob must come from a single connected component C, and that blob must contain all the sites from C. Compatible sites need not be inside any blob. (Gusfield et al 2003-5)
Galled-Tree Structure So when each blob contains only a single cycle, there is a one-one correspondence between the blobs and the non-trivial connected components of the incompatibility graph. This is the central fact used in polynomial-time solutions to the (NP-hard) recombination minimization problem, when a galled-tree for M exists. Motivating Question: To what extent does this clean one-one structure carry over to general phylogenetic networks? How do we exploit the general structure?
The Decomposition Theorem (Recomb 2005) For any set of sequences M, there is a blobbed-tree T(M) that derives M, where each blob contains all and only the sites in one non-trivial connected component of G(M). The compatible sites can always be put on edges outside of any blob. A blobbed-tree with this structure is called fully-decomposed.
General Structure So, for any set of sequences M, there is a phylogenetic network where there is a one-one correspondence between the blobs and the non-trivial connected components of G(M). Moreover, the tree part of T(M) is unique. And it is easy to find the tree part.
A fully-decomposed network for the sequences generated by the prior network. Incompatibility Graph 4 4 3 1 3 2 5 1 s p a: 00010 2 c: 00100 b: 10010 d: 10100 2 5 s p 4 g: 00101 e: 01100 f: 01101
Proof Ideas Let C be a connected component of G(M). Define M[C] as the sequences in M restricted to the sites in C.
1 2 3 4 5 a b c d e f g 0 0 0 1 0 1 0 0 1 0 0 0 1 0 0 1 0 1 0 0 0 1 1 0 0 0 1 1 0 1 0 0 1 0 1 4 C2 C1 M 1 3 2 5 B1 B2 1 3 4 2 5 a b c d e f g 0 0 0 0 0 0 0 0 1 0 1 1 0 1 a 0 0 1 b 1 0 1 c 0 1 0 d 1 1 0 e 0 1 0 f 0 1 0 g 0 1 0 M[C1] M[C2]
Faux Proof Pick one site from each connected component C in G(M) to ``represent” C. No pair of those sites are incompatible, so by the NASC for a perfect phylogeny, there will be a perfect phylogeny T for the sites. Expand each node to a network generating the sequences in M[C]. Incorrect, because the structure of T can be wrong. We need to use information about all the sites in each C.
Now for each connected component C in G(M), call each distinct sequence in M[C] a supercharacter, and let W be the indicator matrix for the supercharacters. So W indicates which rows of M contain which particular supercharacters. 1 2 3 4 5 6 7 8 1 3 4 2 5 a b c d e f g 1 0 0 0 1 0 0 0 0 1 0 0 1 0 0 0 0 0 1 0 1 0 0 0 0 0 0 1 1 0 0 0 0 0 1 0 0 1 0 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 0 1 5 5 5 5 6 7 8 a b c d e f g 0 0 0 0 0 0 0 0 1 0 1 1 0 1 1 2 3 4 3 3 3 a 0 0 1 b 1 0 1 c 0 1 0 d 1 1 0 e 0 1 0 f 0 1 0 g 01 0 W M[C1] M[C2]
Proof Ideas Lemma: No pair of supercharacters are incompatible. So by the NASC for a Perfect Phylogeny, there is a unique perfect phylogeny T for W.
Proof Ideas For each connected component C of G(M), all supercharacters that originate from C label edges in T that are incident with one single node v[C] in T. So, if we expand each node v[C] to be a network that generates the supercharacters from C (the sequences in M[C]), and connect each network correctly to the edges in T, the resulting network is a fully-decomposed blobbed-tree that generates M.
Algorithmically, T is easy to find and is the tree resulting from contracting each blob in the fully-decomposed blobbed-tree T(M) for M. T can be constructed from M in O(nm^2) time.
Broader Biological Applications Our major interest is in recombination, but the proof of the decomposition theorem does not explicitly use recombination. So it holds for whatever biological phenomena caused the incompatibility of sites. For example, back or recurrent mutation, gene-conversion, lateral gene transfer etc.
What is the most tree-like network? • Simple definition: The ``treeness’’ of a network is the number of edges in the tree after contracting each blob to a single node. • Simple fact: In any phylogenetic network N for M, all sites from a single non-trivial connected component must be together in a single blob of N. • Hence, under this simple definition, a fully-decomposed blobbed tree is the most tree-like network for M.
The supercharacters from M play the role in phylogenetic networks that normal binary characters play in perfect phylogeny trees. So supercharacters are the fundamental characters of phylogenetic networks.
The main open question The Decomposition Theorem says there is always a fully-decomposed blobbed-tree for any M, but Is there always a fully-decomposed blobbed-tree that minimizes the number of recombinations over all possible phylogenetic networks for M?
We conjecture the answer is yes. If true, then we can decompose the problem of minimizing the total number of recombinations into separate problems on each connected component, and also find lower bounds on the needed number of recombinations, in each component separately, adding those bounds to get a valid overall lower bound for M. This computation of lower bounds is known to be correct for certain lower bounds (Bafna, Bansal 2004).
Progress on Proving the Conjecture Definition: If N is a phylogenetic network for M, and a node v in N is labeled with a sequence in M, then v is said to be visible in N. Theorem: If every node in N is visible, then there is a fully-decomposed network for M where the number of recombinations is at most the number in N. Corollary: The conjecture is true for any M where the Haplotype or History lower bounds (S. Myers) on the number of recombinations needed to generate M, is tight.
Papers and Software wwwcsif.cs.ucdavis.edu/~gusfield/