Molecular Biology and Evolution 17:530-539 (2000)
© 2000 Society for Molecular Biology and Evolution
Regular Articles |
The Performance of Several Multiple-Sequence Alignment Programs in Relation to Secondary-Structure Features for an rRNA Sequence
*Department of Genetics and Molecular Biology, University of Hawaii at Manoa;
Department of Ecology and Evolutionary Biology, University of Connecticut;
and
Department of Mathematics, Massey University, Palmerston North, New Zealand
Abstract
The performances of five global multiple-sequence alignment programs (CLUSTAL W, Divide and Conquer, Malign, PileUp, and TreeAlign) were evaluated using part of the animal mitochondrial small subunit (12S) rRNA molecule. Conserved sequence motifs derived from an alignment based on secondary structural information were used to score how well each program aligned a data set of five vertebrate and five invertebrate taxa over a range of parameter values. All of the programs could align the motifs with reasonable accuracy for at least one set of parameter conditions, although if the whole sequence was considered, similarity to the structural alignment was only 25%–34%. Use of small gap costs generally gave more accurate results, although Malign and TreeAlign generated longer alignments when gap costs were low. The programs differed in the consistency of the alignments when gap cost was varied; CLUSTAL W, Divide and Conquer, and TreeAlign were the most accurate and robust, while PileUp performed poorly as gap cost values increased, and the accuracy of Malign fluctuated. Default settings for the programs did not give the best results, and attempting to select similar parameter values in different programs did not always result in more similar alignments. Poor alignment of even well-conserved motifs can occur if these are near sites with insertions or deletions. Since there is no a priori way to determine gap costs and because such costs can vary over the gene, alignment of rRNA sequences, particularly the less well conserved regions, should be treated carefully and aided by secondary structure and conserved motifs. Some motifs are single bases and so are often invisible to alignment programs. Our tests involved the most conserved regions of the 12S rRNA gene, and alignment of less well conserved regions will be more problematical. None of the alignments we examined produced a fully resolved phylogeny for the data set, indicating that this portion of 12S rRNA is insufficient for resolution of distant evolutionary relationships.
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  G. Landan and D. Graur Heads or Tails: A Simple Reliability Check for Multiple Sequence Alignments Mol. Biol. Evol., June 1, 2007; 24(6): 1380 - 1383. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 T. Lassmann and E. L. L. Sonnhammer Automatic assessment of alignment quality Nucleic Acids Res., December 16, 2005; 33(22): 7120 - 7128. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 S. A. Kelchner Group II introns as phylogenetic tools: structure, function, and evolutionary constraints Am. J. Botany, October 1, 2002; 89(10): 1651 - 1669. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 E. Carosa, Z. Kozmik, J. E. Rall, and J. Piatigorsky Structure and Expression of the Scallop Omega -Crystallin Gene. EVIDENCE FOR CONVERGENT EVOLUTION OF PROMOTER SEQUENCES J. Biol. Chem., January 4, 2002; 277(1): 656 - 664. [Abstract] [Full Text]
|
|
|
|
|
|
R. D. M. Page Comparative analysis of secondary structure of insect mitochondrial small subunit ribosomal RNA using maximum weighted matching Nucleic Acids Res., October 15, 2000; 28(20): 3839 - 3845. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. Piatigorsky, Z. Kozmik, J. Horwitz, L. Ding, E. Carosa, W. G. Robison Jr., P. J. Steinbach, and E. R. Tamm Omega -Crystallin of the Scallop Lens. A DIMERIC ALDEHYDE DEHYDROGENASE CLASS 1/2 ENZYME-CRYSTALLIN J. Biol. Chem., December 22, 2000; 275(52): 41064 - 41073. [Abstract] [Full Text] [PDF]
|
|