Quartets and unrooted level-k networks

Similar documents
Paths. Connectivity. Euler and Hamilton Paths. Planar graphs.

LEO VAN IERSEL TU DELFT

b. How many ternary words of length 23 with eight 0 s, nine 1 s and six 2 s?

Constructive Geometric Constraint Solving

Math 61 : Discrete Structures Final Exam Instructor: Ciprian Manolescu. You have 180 minutes.

More Foundations. Undirected Graphs. Degree. A Theorem. Graphs, Products, & Relations

Computational Biology, Phylogenetic Trees. Consensus methods

12/3/12. Outline. Part 10. Graphs. Circuits. Euler paths/circuits. Euler s bridge problem (Bridges of Konigsberg Problem)

5/9/13. Part 10. Graphs. Outline. Circuits. Introduction Terminology Implementing Graphs

Algorithmic and NP-Completeness Aspects of a Total Lict Domination Number of a Graph

Graphs. Graphs. Graphs: Basic Terminology. Directed Graphs. Dr Papalaskari 1

An undirected graph G = (V, E) V a set of vertices E a set of unordered edges (v,w) where v, w in V

V={A,B,C,D,E} E={ (A,D),(A,E),(B,D), (B,E),(C,D),(C,E)}

Graphs. CSC 1300 Discrete Structures Villanova University. Villanova CSC Dr Papalaskari

Graph Isomorphism. Graphs - II. Cayley s Formula. Planar Graphs. Outline. Is K 5 planar? The number of labeled trees on n nodes is n n-2

arxiv: v1 [cs.ds] 20 Feb 2008

V={A,B,C,D,E} E={ (A,D),(A,E),(B,D), (B,E),(C,D),(C,E)}

CSE 373. Graphs 1: Concepts, Depth/Breadth-First Search reading: Weiss Ch. 9. slides created by Marty Stepp

Planar Upward Drawings

CS200: Graphs. Graphs. Directed Graphs. Graphs/Networks Around Us. What can this represent? Sometimes we want to represent directionality:

Outline. Circuits. Euler paths/circuits 4/25/12. Part 10. Graphs. Euler s bridge problem (Bridges of Konigsberg Problem)

0.1. Exercise 1: the distances between four points in a graph

(a) v 1. v a. v i. v s. (b)

Section 10.4 Connectivity (up to paths and isomorphism, not including)

Solutions for HW11. Exercise 34. (a) Use the recurrence relation t(g) = t(g e) + t(g/e) to count the number of spanning trees of v 1

Cycles and Simple Cycles. Paths and Simple Paths. Trees. Problem: There is No Completely Standard Terminology!

CSE 373: More on graphs; DFS and BFS. Michael Lee Wednesday, Feb 14, 2018

The University of Sydney MATH2969/2069. Graph Theory Tutorial 5 (Week 12) Solutions 2008

CSC Design and Analysis of Algorithms. Example: Change-Making Problem

learning objectives learn what graphs are in mathematical terms learn how to represent graphs in computers learn about typical graph algorithms

Problem solving by search

1. Determine whether or not the following binary relations are equivalence relations. Be sure to justify your answers.

12. Traffic engineering

Announcements. Not graphs. These are Graphs. Applications of Graphs. Graph Definitions. Graphs & Graph Algorithms. A6 released today: Risk

Module graph.py. 1 Introduction. 2 Graph basics. 3 Module graph.py. 3.1 Objects. CS 231 Naomi Nishimura

CSI35 Chapter 11 Review

N=4 L=4. Our first non-linear data structure! A graph G consists of two sets G = {V, E} A set of V vertices, or nodes f

COMP108 Algorithmic Foundations

CS 241 Analysis of Algorithms

CS61B Lecture #33. Administrivia: Autograder will run this evening. Today s Readings: Graph Structures: DSIJ, Chapter 12

Trees as operads. Lecture A formalism of trees

A Simple Code Generator. Code generation Algorithm. Register and Address Descriptors. Example 3/31/2008. Code Generation

Outline. 1 Introduction. 2 Min-Cost Spanning Trees. 4 Example

CS 461, Lecture 17. Today s Outline. Example Run

Similarity Search. The Binary Branch Distance. Nikolaus Augsten.

QUESTIONS BEGIN HERE!

, each of which is a tree, and whose roots r 1. , respectively, are children of r. Data Structures & File Management

CS September 2018

Why the Junction Tree Algorithm? The Junction Tree Algorithm. Clique Potential Representation. Overview. Chris Williams 1.

COMPLEXITY OF COUNTING PLANAR TILINGS BY TWO BARS

# 1 ' 10 ' 100. Decimal point = 4 hundred. = 6 tens (or sixty) = 5 ones (or five) = 2 tenths. = 7 hundredths.

5/7/13. Part 10. Graphs. Theorem Theorem Graphs Describing Precedence. Outline. Theorem 10-1: The Handshaking Theorem

QUESTIONS BEGIN HERE!

Numbering Boundary Nodes

Present state Next state Q + M N

Exam 1 Solution. CS 542 Advanced Data Structures and Algorithms 2/14/2013

Outline. Computer Science 331. Computation of Min-Cost Spanning Trees. Costs of Spanning Trees in Weighted Graphs

Announcements. These are Graphs. This is not a Graph. Graph Definitions. Applications of Graphs. Graphs & Graph Algorithms

From graph classes to phylogenetic networks Philippe Gambette

MULTIPLE-LEVEL LOGIC OPTIMIZATION II

MAT3707. Tutorial letter 201/1/2017 DISCRETE MATHEMATICS: COMBINATORICS. Semester 1. Department of Mathematical Sciences MAT3707/201/1/2017

The University of Sydney MATH 2009

NP-Completeness. CS3230 (Algorithm) Traveling Salesperson Problem. What s the Big Deal? Given a Problem. What s the Big Deal? What s the Big Deal?

ECE 407 Computer Aided Design for Electronic Systems. Circuit Modeling and Basic Graph Concepts/Algorithms. Instructor: Maria K. Michael.

Counting Paths Between Vertices. Isomorphism of Graphs. Isomorphism of Graphs. Isomorphism of Graphs. Isomorphism of Graphs. Isomorphism of Graphs

Chapter 9. Graphs. 9.1 Graphs

On Contract-and-Refine Transformations Between Phylogenetic Trees

Lecture 20: Minimum Spanning Trees (CLRS 23)

Graph Theory. Vertices. Vertices are also known as nodes, points and (in social networks) as actors, agents or players.

d e c b a d c b a d e c b a a c a d c c e b

Graph width-parameters and algorithms

1 Introduction to Modulo 7 Arithmetic

A 43k Kernel for Planar Dominating Set using Computer-Aided Reduction Rule Discovery

ECE COMBINATIONAL BUILDING BLOCKS - INVEST 13 DECODERS AND ENCODERS

Discovering Pairwise Compatibility Graphs

Weighted Graphs. Weighted graphs may be either directed or undirected.

Outline. Binary Tree

In which direction do compass needles always align? Why?

IEEE TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS, VOL. TK, NO. TK, MONTHTK YEARTK 1. Hamiltonian Walks of Phylogenetic Treespaces

CSE303 - Introduction to the Theory of Computing Sample Solutions for Exercises on Finite Automata

Can transitive orientation make sandwich problems easier?

Prefix-Free Regular-Expression Matching

Motivation. Problems. Puzzle 1

Weighted graphs -- reminder. Data Structures LECTURE 15. Shortest paths algorithms. Example: weighted graph. Two basic properties of shortest paths

Solutions for HW9. Bipartite: put the red vertices in V 1 and the black in V 2. Not bipartite!

CIT 596 Theory of Computation 1. Graphs and Digraphs

10/30/12. Today. CS/ENGRD 2110 Object- Oriented Programming and Data Structures Fall 2012 Doug James. DFS algorithm. Reachability Algorithms

Nefertiti. Echoes of. Regal components evoke visions of the past MULTIPLE STITCHES. designed by Helena Tang-Lim

Phylogenetic networks: overview, subclasses and counting problems Philippe Gambette

Multi-Section Coupled Line Couplers

Register Allocation. How to assign variables to finitely many registers? What to do when it can t be done? How to do so efficiently?

Greedy Algorithms, Activity Selection, Minimum Spanning Trees Scribes: Logan Short (2015), Virginia Date: May 18, 2016

OpenMx Matrices and Operators

Seven-Segment Display Driver

DUET WITH DIAMONDS COLOR SHIFTING BRACELET By Leslie Rogalski

A 4-state solution to the Firing Squad Synchronization Problem based on hybrid rule 60 and 102 cellular automata

Math 166 Week in Review 2 Sections 1.1b, 1.2, 1.3, & 1.4

Solutions to Homework 5

Multipoint Alternate Marking method for passive and hybrid performance monitoring

XML and Databases. Outline. Recall: Top-Down Evaluation of Simple Paths. Recall: Top-Down Evaluation of Simple Paths. Sebastian Maneth NICTA and UNSW

Transcription:

Phylogntis Workshop, Is Nwton Institut or Mthmtil Sins Cmrig 21/06/2011 Qurtts n unroot lvl-k ntworks Philipp Gmtt

Outlin Astrt n xpliit phylognti ntworks Lvl-k ntworks Unroot lvl-1 ntworks n irulr split systms Ronstrution rom triplts n qurtts

Outlin Astrt n xpliit phylognti ntworks Lvl-k ntworks Unroot lvl-1 ntworks n irulr split systms Ronstrution rom triplts n qurtts

Astrt n xpliit phylognti ntworks Astrt phylognti ntworks: Visuliz rtiult volution t Expliit phylognti ntworks: Eh no n intrprt s urrnt or nstrl spis split ntwork SplitsTr lvl-2 ntwork Simplisti min ntwork gll ntwork Ntwork Dnrosop

Split ntworks Astrt phylognti ntworks: Visuliz rtiult volution t Split ntwork: Visuliz st o splits Split 1 Björn M. Hllström, Axl Jnk - Mmmlin volution my not stritly iurting, MBE, 2010

Split ntworks Astrt phylognti ntworks: Visuliz rtiult volution t Split ntwork: Visuliz st o splits Split 2 Björn M. Hllström, Axl Jnk - Mmmlin volution my not stritly iurting, MBE, 2010

Split ntworks Astrt phylognti ntworks: Visuliz rtiult volution t Split ntwork: Visuliz st o splits Split 3 Björn M. Hllström, Axl Jnk - Mmmlin volution my not stritly iurting, MBE, 2010

Ntwork sulss hirrhy split systm xpliit root wkly omptil k-omptil rgulr tr-siling tr-hil lvl k lvl-2 gll ntwork nst 2-omptil norml lvl-1 irulr uniyli omptil tr simpl lvl-1 split systms ontins xpliit root ntworks http://www.lirmm.r/~gmtt/rphylogntintworks.php

Ntwork sulss hirrhy split systm xpliit root wkly omptil k-omptil rgulr tr-siling tr-hil lvl k lvl-2 gll ntwork nst 2-omptil norml lvl-1 irulr uniyli omptil tr simpl lvl-1 split systms ontins xpliit root ntworks http://www.lirmm.r/~gmtt/rphylogntintworks.php

Outlin Astrt n xpliit phylognti ntworks Lvl-k ntworks Unroot lvl-1 ntworks n irulr split systms Ronstrution rom triplts n qurtts

Lvl-k ntworks lvl: how r is th ntwork rom tr? smll lvl tr strutur st lgorithms h 1 h 3 h 2 lvl = mximum numr o rtiultions y riglss omponnt (lo) o th unrlying unirt grph. g h i lvl-2 ntwork j k Choy, Jnsson, Skn & Sung, TCS, 2005

Lvl-k ntworks lvl: how r is th ntwork rom tr? smll lvl tr strutur st lgorithms lvl = mximum numr o rtiultions y lo. g h i lvl-2 ntwork j k lvl-1 ntwork ( gll tr ) g h i j k

Unroot lvl-k ntworks lvl: how r is th ntwork rom n unroot tr? smll lvl tr strutur st lgorithms i h g lvl = mximum numr o gs to rmov, y lo, to otin tr. unroot lvl-2 ntwork Gmtt, Brry & Pul, mnusript, 2011

Unroot lvl-k ntworks lvl: how r is th ntwork rom n unroot tr? smll lvl tr strutur st lgorithms i h g lvl = mximum numr o gs to rmov, y lo, to otin tr. = mximum ylomti numr o th los unroot lvl-2 ntwork

Unroot lvl-k ntworks lvl: how r is th ntwork rom n unroot tr? smll lvl tr strutur st lgorithms i h g lvl = mximum numr o gs to rmov, y lo, to otin tr. unroot lvl-1 ntwork (unroot gll tr) tr o yls

Unroot lvl-k ntworks i h g lvl = mximum numr o gs to rmov, y lo, to otin tr. unroot lvl-k ntwork tr o los tr o gnrtors o lvl k Unroot lvl-k gnrtors: riglss looplss 3-rgulr multigrphs with 2k- 2 vrtis lvl-2 gnrtor lvl-3 gnrtors Gmtt, Brry & Pul, mnusript, 2011

Counting ll lvl-k ntworks Unroot lvl-1 ntworks: xpliit ormul or n lvs, yls, m gs ross yls Smpl & Stl, TCBB, 2006

Counting ll lvl-k ntworks Unroot lvl-1 ntworks: xpliit ormul or n lvs, yls, m gs ross yls + symptoti vlution or n lvs: 0.207 n n-1 1.890 n Root lvl-1 ntworks : Expliit ormul or n lvs, yls, m gs ross yls + symptoti vlution or n lvs: 0.134 2.943 n n n-1 Smpl & Stl, TCBB, 2006 Unroot lvl-2 ntworks : Expliit ormul or n lvs : n+i-1 4i+j-1 i k p q -4 s 9 i 23 k -21 q (n-1)! Σ i j k p q s 21 2 9 (-1) p 46 0 s q p k i n-1 j=n-1-i-k-p-q-s 0 i 0 numr o lvs 2 3 4 5 6 7 unroot lvl-1-2 15 192 3 450 79 740 root lvl-1 3 36 723 20 280 730 755 32 171 580 unroot lvl-2-9 282 14 697 1 071 720 100 461 195 Bouvl, Gmtt, Brry & Pul, mnusript, 2011

Equivln twn root n unroot lvl g i h rooting Rooting: - hoosing root - hoosing n orinttion or th gs g h i Gmtt, Brry & Pul, mnusript, 2011

Equivln twn root n unroot lvl g i h rooting Rooting: - hoosing root - hoosing n orinttion or th gs g h i Gmtt, Brry & Pul, mnusript, 2011

Equivln twn root n unroot lvl g i h rooting Rooting: - hoosing root - hoosing n orinttion or th gs g h i - mny possil rootings (possily xponntil in th lvl) - sm lvl (invrint) Gmtt, Brry & Pul, mnusript, 2011

Outlin Astrt n xpliit phylognti ntworks Lvl-k ntworks Unroot lvl-1 ntworks n irulr split systms Ronstrution rom triplts n qurtts

Splits in unroot lvl-k ntworks Split: Split o tr ontin in th ntwork? Lvs sprt y miniml ut in th ntwork? Woolly, Pos & Crnll, PLoS On, 2008 Brns & Cornlsn, DAM, 2009 i h g

Splits in unroot lvl-k ntworks Biprtition : Split o tr ontin in th ntwork Lvs sprt y miniml ut in th ntwork i h g quivlnt! i h g Gmtt, Brry & Pul, mnusript, 2011

Rprsnting splits y unroot ntworks rl ntwork g h i trs g h i g h i

Rprsnting splits y unroot ntworks rl ntwork g h i trs g h i g h i splits SplitsTr i h g split ntwork

Rprsnting splits y unroot ntworks rl ntwork unroot lvl-1 ntwork i g h i h trs splits g h i splits g h i SplitsTr ghi i g ghi h g split ntwork

Splits in unroot lvl-1 ntworks Splits Σ(N) o lvl-1 ntwork i h g N i h g Cirulr split systm Σ

Splits in unroot lvl-1 ntworks Splits Σ(N) o lvl-1 ntwork i h g Σ(N) irulr N i h g Cirulr split systm Σ Gmtt, Brry & Pul, mnusript, 2011

Splits in unroot lvl-1 ntworks Splits Σ(N) o lvl-1 ntwork i h g Σ(N) irulr N i h g Cirulr split systm Σ Thr xists lvl-1 ntwork N suh tht Σ Σ(N) Gmtt, Brry & Pul, mnusript, 2011

Outlin Astrt n xpliit phylognti ntworks Lvl-k ntworks Unroot lvl-1 ntworks n irulr split systms Ronstrution rom triplts n qurtts

Comintoril phylognti ntwork ronstrution spis 1 : AATTGCAG TAGCCCAAAAT spis 2 : ACCTGCAG TAGACCAAT spis 3 : GCTTGCCG TAGACAAGAAT spis 4 : ATTTGCAG AAGACCAAAT spis 5 : TAGACAAGAAT spis 6 : ACTTGCAG TAGCACAAAAT spis 7 : ACCTGGTG TAAAAT G1 G2 {gn squns} {trs} T1 T2 HOGENOM ts Duyr, Durt, Pnl, Gouy, Rhnmnn & Prrièr, BioIn, 2005 ntwork ontins th trs + optiml

Comintoril phylognti ntwork ronstrution spis 1 : AATTGCAG TAGCCCAAAAT spis 2 : ACCTGCAG TAGACCAAT spis 3 : GCTTGCCG TAGACAAGAAT spis 4 : ATTTGCAG AAGACCAAAT spis 5 : TAGACAAGAAT spis 6 : ACTTGCAG TAGCACAAAAT spis 7 : ACCTGGTG TAAAAT G1 G2 {gn squns} {trs} T1 T2 HOGENOM ts Duyr, Durt, Pnl, Gouy, Rhnmnn & Prrièr, BioIn, 2005 > 500 spis, >70 000 trs ntwork ontins th trs + optiml NP-omplt or 2 root trs Borwih & Smpl, DAM, 2007

Ronstrution rom triplts / qurtts {gn squns} {trs} {triplts} {qurtts} ntwork ontins th qurtts/triplts + optiml

Ronstrution rom triplts / qurtts {gn squns} {trs} {triplts} {qurtts} triplt ntwork ontins th qurtts/triplts + optiml

Ronstrution rom triplts / qurtts {gn squns} {trs} triplt {triplts} {qurtts} qurtt ntwork ontins th qurtts/triplts + optiml

Ronstrution rom triplts / qurtts {trs} {triplts} {qurtts} ntwork N ntwork N' N'=N?

Ronstrution rom triplts / qurtts A ntwork ontining ll qurtts o tr T os not lwys ontin T. T N N ontins ll qurtts o T ut not T

Ronstrution rom triplts / qurtts {trs} {triplts} {qurtts} ntwork N ntwork N' N'=N? Not lwys, ut {N} {N'}

Ronstrution rom triplts / qurtts Fining ll triplts o root ntwork: O(n 3 ) Byrk, Gwryhowski, Hur & Klk, JDA, 2010

Ronstrution rom triplts / qurtts Fining ll triplts o root ntwork: O(n 3 ) Byrk, Gwryhowski, Hur & Klk, JDA, 2010 Fining ll qurtts o n unroot ntwork? i h g qurtt

Ronstrution rom triplts / qurtts Fining ll triplts o root ntwork: O(n 3 ) Byrk, Gwryhowski, Hur & Klk, JDA, 2010 Fining ll qurtts o n unroot ntwork? i h g qurtt 2-isjoint pths -,-

Ronstrution rom triplts / qurtts Fining ll triplts o root ntwork: O(n 3 ) Byrk, Gwryhowski, Hur & Klk, JDA, 2010 Fining ll qurtts o n unroot ntwork: O(n 6 ) 2-Disjoint Pths in grph o gr 3: O(n(1+α(n,n))) h g i Tholy, SOFSEM'09, 2009 qurtt 2-isjoint pths -,-

Tr ronstrution unroot rom qurtts root rom triplts gnrl NP-omplt Stl, JOC, 1992 polynomil Aho, Sgiv, Szymnski & Ullmn, SJOC, 1981 Hnzingr, King & Wrnow, ALG, 1999 Jnsson, Ng, Skn & Sung, ALG, 2005 ns t lst on qurtt or h st o 4 lvs O(n 4 ) Brry & Gsul, TCS, 2000 O(n 3 ) Aho t l., SJOC, 1981

Ronstrution o lvl-k ntworks unroot rom qurtts root rom triplts lvl 1 lvl k>1 lvl 1 lvl k>1 gnrl NP-omplt Grünwl, Moulton & Spillnr, DAM, 2009? NP-omplt Jnsson, Nguyn & Sung, SJOC, 2006 NP-omplt Vn Irsl, Klk & Mnih, JBCB, 2009 ns t lst on qurtt or h st o 4 lvs? (omposition in polynomil tim)? O(n 3 ) Jnsson, Nguyn & Sung, SJOC, 2006 O(n 5k+4 ) To & Hi, CPM'09 omplt ll qurtts o th ntwork O(n 4 )? (omposition in polynomil tim) O(n 3 ) Jnsson, Nguyn & Sung, SJOC, 2006 O(n 3k+3 ) Vn Irsl & Klk, ALG, 2010 Gmtt, Brry & Pul, mnusript, 2011

Qurtt st omposition A A SN-split o th qurtt st Q: For ll lvs x,y A, z,t A, Th only qurtt o Q on {x,y,z,t} is xy zt Root ontxt: SN-st Jnsson & Sung, TCS, 2006 To & Hi, CPM'09 Q is ns qurtt st SN-splits o Q r omptil (n rprsnt y n unroot tr) Computing th SN-splits: O(n 4 ) vrint o th Q* lgorithm Brry & Gsul, TCS, 2001 Root ontxt: Computing th SN-sts: O(n 3 ) Jnsson, Nguyn & Sung, TCS, 2006 Gmtt, Brry & Pul, mnusript, 2011

Qurtt st omposition Q(N), th st o ll qurtts ontin in n unroot lvl-k ntwork N SN-splits o Q(N) r ijtivly ssoit with ut-gs o N. N h g gh gh gh gh h g gh gh gh g h gh gh Gmtt, Brry & Pul, mnusript, 2011

Blo ronstrution rom qurtts W n sprt th los o N rom Q(N). How to ronstrut h lo o N? Lvl-1 ntwork ronstrution rom th omplt qurtt st: N, lvl-1 ntwork Fining n orring o th lvs roun th yl: O(n 2 ) Qurtt irulr puzzling Fix our lvs,,,:

Blo ronstrution rom qurtts W n sprt th los o N rom Q(N). How to ronstrut h lo o N? Lvl-1 ntwork ronstrution rom th omplt qurtt st: N, lvl-1 ntwork Fining n orring o th lvs roun th yl: O(n 2 ) Qurtt irulr puzzling Fix our lvs,,, (, ) For h othr l: For h possil position on th yl: Tst in O(1) i it is th orrt on.

Blo ronstrution rom qurtts W n sprt th los o N rom Q(N). How to ronstrut h lo o N? Lvl-1 ntwork ronstrution rom th omplt qurtt st: N, lvl-1 ntwork Fining n orring o th lvs roun th yl: O(n 2 ) x Qurtt irulr puzzling Fix our lvs,,, (, ) For h othr l: For h possil position on th yl: Tst in O(1) i it is th orrt on. x + x?

Link with non-btwnnss Btwnnss: Input: - st X o lmnts - st C o twnnss onstrints: twn n Output: - orring σ whih rspts th onstrints o C

Link with non-btwnnss Non-Btwnnss: Input: - st X o lmnts - st C o non-twnnss onstrints: not twn n Output: - orring σ whih rspts th onstrints o C

Link with non-btwnnss Cirulr Non-Btwnnss: Input: - st X o lmnts - st C o irulr non-twnnss onstrints: sn rom, not twn n Output: - irulr orring σ whih rspts th onstrints o C

Link with non-btwnnss Cirulr Non-Btwnnss: Input: - st X o lmnts - st C o irulr non-twnnss onstrints: sn rom, not twn n Output: - irulr orring σ whih rspts th onstrints o C quivlnt to unroot lvl-1 lo ronstrution rom qurtts

Link with non-btwnnss Gnrl s Dns s Btwnnss Non-Btwnnss Cirulr Non-Btwnnss NP-omplt Optrny, JOC, 1979 NP-omplt Guttmn & Muhr, IFIP-TCS'06 NP-omplt Grünwl, Moulton & Spillnr, DAM, 2009 polynomil??

Thnk you! Couthors o ths rsults: Vinnt Brry & Christoph Pul (Montpllir), Mthil Bouvl (Borux) Unroot lvl-2 ntwork uilt rom th mnusript with TrClou SplitsTr T-Rx http://www.trlou.org