Multiple Sequence Alignment Using Partial Order Graphs





Christopher Lee*, Catherine Grasso, Mark Sharlow
Department of Chemistry and Biochemistry
University of California, Los Angeles
Los Angeles, CA 90095-1570

TEL: 310-825-7374
FAX: 310-267-0248
EMAIL: leec@mbi.ucla.edu

*Author to whom correspondence should be addressed.


Draft August 2, 2001
Printed August 2, 2001


ABSTRACT
Motivation: Progressive multiple sequence alignment (MSA) methods depend on reducing an MSA to a linear profile for each alignment step.  However, this leads to loss of information needed for accurate alignment, and gap scoring artifacts.  
Results: We present a graph representation of an MSA that can itself be aligned directly by pairwise dynamic programming, eliminating the need to reduce the MSA to a profile.  This enables our algorithm (POA: partial order alignment) to guarantee that the optimal alignment of each new sequence versus each sequence in the MSA will be considered.  Moreover, this algorithm introduces a new edit operator, homologous recombination, important for multidomain sequences.  The algorithm has improved speed (linear time complexity) over existing MSA algorithms, enabling construction of massive and complex alignments (e.g. an alignment of 5000 sequences in 4 hours on a Pentium II).  We demonstrate the utility of this algorithm on a family of multidomain SH2 proteins, and on EST assemblies containing alternative splicing and polymorphism.  
Availability: The partial order alignment program POA is available at http://www.bioinformatics.ucla.edu/poa.  
Contact: Christopher Lee, leec@mbi.ucla.edu.