Melvin's digital garden

Ancestral reconstruction over continuous characters and distributions

CREATED: 200810140800 TITLE: Ancestral reconstruction by asymmetric Wagner parsimony over continuous characters and squared parsimony over distributions SPEAKER: Miklos Csuros ** Squared parsimony

  • Evolutionary character may be over real values, for example probability values
  • Parsimony labeling (steiner tree labeling), can be solved by DP in the discrete case (Sankoff and Rousseau 1975)
  • Squared parsimony, $\Delta(x,y) = (x - y)^2$
  • SquaP for a distribution can be computed by Maddison’s algorithm in each position independently ** Gene family evolution using Wagner parsimony
  • events: duplication, speciation, gene gain/loss
  • character is the size of the gene family, eg COG0247 (Kinesin-like protein)
  • Wagner parsimony: $\Delta(x,y) = |x - y|$
  • Asymmetric Wagner parsimony with gain loss penalty: $\Delta(x,y) = (y > x) ? a(y - x) : b(x - y)$
  • For a tree of n nodes and height h, can be computed in $O(nh)$
  • symmetric Wagner is like a cup, asymmetric Wagner parsimony is a piece-wise linear, continuous convex function
  • can be computed in postorder traversal ** Gene content evolution in Archaea
  • arCOG data set: 7538 gene families (Makarova et al. 2007)
  • gains followed by specialization

Links to this note