RNA Sequence Evolution With Secondary Structure Constraints: Comparison of Substitution Rate Models Using Maximum-Likelihood Methods

Corresponding author: Paul G. Higgs, School of Biological Sciences, University of Manchester, Oxford Rd., Manchester M13 9PT, United Kingdom., paul.higgs{at}man.ac.uk (E-mail)

Communicating editor: S. YOKOYAMA

We test models for the evolution of helical regions of RNA sequences, where the base pairing constraint leads to correlated compensatory substitutions occurring on either side of the pair. These models are of three types: 6-state models include only the four Watson-Crick pairs plus GU and UG; 7-state models include a single mismatch state that combines all of the 10 possible mismatches; 16-state models treat all mismatch states separately. We analyzed a set of eubacterial ribosomal RNA sequences with a well-established phylogenetic tree structure. For each model, the maximum-likelihood values of the parameters were obtained. The models were compared using the Akaike information criterion, the likelihood-ratio test, and Cox's test. With a high significance level, models that permit a nonzero rate of double substitutions performed better than those that assume zero double substitution rate. Some models assume symmetry between GC and CG, between AU and UA, and between GU and UG. Models that relaxed this symmetry assumption performed slightly better, but the tests did not all agree on the significance level. The most general time-reversible model significantly outperformed any of the simplifications. We consider the relative merits of all these models for molecular phylogenetics.

This article has been cited by other articles:

				A. Varriale, G. Torelli, and G. Bernardi Compositional properties and thermal adaptation of 18S rRNA in vertebrates RNA, August 1, 2008; 14(8): 1492 - 1500. [Abstract] [Full Text] [PDF]

				J. L. Thorne, S. C. Choi, J. Yu, P. G. Higgs, and H. Kishino Population Genetics Without Intraspecific Data Mol. Biol. Evol., August 1, 2007; 24(8): 1667 - 1677. [Abstract] [Full Text] [PDF]

				C. Kosiol, I. Holmes, and N. Goldman An Empirical Codon Model for Protein Sequence Evolution Mol. Biol. Evol., July 1, 2007; 24(7): 1464 - 1479. [Abstract] [Full Text] [PDF]

				L. M. F. Merlo, M. Lunzer, and A. M. Dean An empirical test of the concomitantly variable codon hypothesis PNAS, June 26, 2007; 104(26): 10938 - 10943. [Abstract] [Full Text] [PDF]

				V. Gowri-Shankar and M. Rattray A Reversible Jump Method for Bayesian Phylogenetic Inference with a Nonhomogeneous Substitution Model Mol. Biol. Evol., June 1, 2007; 24(6): 1286 - 1299. [Abstract] [Full Text] [PDF]

				S. L. Kosakovsky Pond, F. V. Mannino, M. B. Gravenor, S. V. Muse, and S. D. W. Frost Evolutionary Model Selection with a Genetic Algorithm: A Case Study Using Stem RNA Mol. Biol. Evol., January 1, 2007; 24(1): 159 - 170. [Abstract] [Full Text] [PDF]

				S. L. Kosakovsky Pond, D. Posada, M. B. Gravenor, C. H. Woelk, and S. D. W. Frost Automated Phylogenetic Detection of Recombination Using a Genetic Algorithm Mol. Biol. Evol., October 1, 2006; 23(10): 1891 - 1901. [Abstract] [Full Text] [PDF]

				G. Baele, J. Raes, Y. Van de Peer, and S. Vansteelandt An Improved Statistical Method for Detecting Heterotachy in Nucleotide Sequences Mol. Biol. Evol., July 1, 2006; 23(7): 1397 - 1405. [Abstract] [Full Text] [PDF]

				V. Gowri-Shankar and M. Rattray On the Correlation Between Composition and Site-Specific Evolutionary Rate: Implications for Phylogenetic Inference Mol. Biol. Evol., February 1, 2006; 23(2): 352 - 364. [Abstract] [Full Text] [PDF]

				J. Dutheil, T. Pupko, A. Jean-Marie, and N. Galtier A Model-Based Approach for Detecting Coevolving Positions in a Molecule Mol. Biol. Evol., September 1, 2005; 22(9): 1919 - 1928. [Abstract] [Full Text] [PDF]

				S. Sengupta and P. G. Higgs A Unified Model of Codon Reassignment in Alternative Genetic Codes Genetics, June 1, 2005; 170(2): 831 - 840. [Abstract] [Full Text] [PDF]

				M. J Telford, M. J Wise, and V. Gowri-Shankar Consideration of RNA Secondary Structure Significantly Improves Likelihood-Based Estimates of Phylogeny: Examples from the Bilateria Mol. Biol. Evol., April 1, 2005; 22(4): 1129 - 1136. [Abstract] [Full Text] [PDF]

				J. S. Pedersen, R. Forsberg, I. M. Meyer, and J. Hein An Evolutionary Model for Protein-Coding Regions with Conserved RNA Structure Mol. Biol. Evol., October 1, 2004; 21(10): 1913 - 1922. [Abstract] [Full Text] [PDF]

				S. Whelan and N. Goldman Estimating the Frequency of Events That Cause Multiple-Nucleotide Changes Genetics, August 1, 2004; 167(4): 2027 - 2043. [Abstract] [Full Text] [PDF]

				Y. Iwasa, F. Michor, and M. A. Nowak Stochastic Tunnels in Evolutionary Dynamics Genetics, March 1, 2004; 166(3): 1571 - 1579. [Abstract] [Full Text] [PDF]

				A. D. Smith, T. W. H. Lui, and E. R. M. Tillier Empirical Models for Substitution in Ribosomal RNA Mol. Biol. Evol., March 1, 2004; 21(3): 419 - 427. [Abstract] [Full Text] [PDF]

				D. Jameson, A. P. Gibson, C. Hudelot, and P. G. Higgs OGRe: a relational database for comparative analysis of mitochondrial genomes Nucleic Acids Res., January 1, 2003; 31(1): 202 - 206. [Abstract] [Full Text] [PDF]

				D. C. Hoyle and P. G. Higgs Factors Affecting the Errors in the Estimation of Evolutionary Distances Between Sequences Mol. Biol. Evol., January 1, 2003; 20(1): 1 - 9. [Abstract] [Full Text] [PDF]

				H. Jow, C. Hudelot, M. Rattray, and P. G. Higgs Bayesian Phylogenetics Using an RNA Substitution Model Applied to Early Mammalian Evolution Mol. Biol. Evol., September 1, 2002; 19(9): 1591 - 1601. [Abstract] [Full Text] [PDF]

				G. Caetano-Anolles Tracing the evolution of RNA structure in ribosomes Nucleic Acids Res., June 1, 2002; 30(11): 2575 - 2587. [Abstract] [Full Text] [PDF]

				D. B. Carlini, Y. Chen, and W. Stephan The Relationship Between Third-Codon Position Nucleotide Content, Codon Bias, mRNA Secondary Structure and Gene Expression in the Drosophilid Alcohol Dehydrogenase Genes Adh and Adhr Genetics, October 1, 2001; 159(2): 623 - 633. [Abstract] [Full Text] [PDF]

				H. Innan and W. Stephan Selection Intensity Against Deleterious Mutations in RNA Secondary Structures and Rate of Compensatory Nucleotide Substitutions Genetics, September 1, 2001; 159(1): 389 - 399. [Abstract] [Full Text] [PDF]