Comparative Sequence Analysis and Patterns of Covariation in RNA Secondary Structures

Corresponding author: Wolfgang Stephan, Department of Biology, University of Rochester, Rochester, NY 14627-0211., stephan{at}troi.cc.rochester.edu (E-mail)

Communicating editor: G. B. GOLDING

A novel method of RNA secondary structure prediction based on a comparison of nucleotide sequences is described. This method correctly predicts nearly all evolutionarily conserved secondary structures of five different RNAs: tRNA, 5S rRNA, bacterial ribonuclease P (RNase P) RNA, eukaryotic small subunit rRNA, and the 3' untranslated region (UTR) of the Drosophila bicoid (bcd) mRNA. Furthermore, covariations occurring in the helices of these conserved RNA structures are analyzed. Two physical parameters are found to be important determinants of the evolution of compensatory mutations: the length of a helix and the distance between base-pairing nucleotides. For the helices of bcd 3' UTR mRNA and RNase P RNA, a positive correlation between the rate of compensatory evolution and helix length is found. The analysis of Drosophila bcd 3' UTR mRNA further revealed that the rate of compensatory evolution decreases with the physical distance between base-pairing residues. This result is in qualitative agreement with Kimura's model of compensatory fitness interactions, which assumes that mutations occurring in RNA helices are individually deleterious but become neutral in appropriate combinations.

This article has been cited by other articles:

				E. BINDEWALD and B. A. SHAPIRO RNA secondary structure prediction from sequence alignments using a network of k-nearest neighbor classifiers. RNA, March 1, 2006; 12(3): 342 - 352. [Abstract] [Full Text] [PDF]

				M. C. Cowperthwaite, J. J. Bull, and L. A. Meyers Distributions of Beneficial Fitness Effects in RNA Genetics, August 1, 2005; 170(4): 1449 - 1457. [Abstract] [Full Text] [PDF]

				T. Castrignano, A. Canali, G. Grillo, S. Liuni, F. Mignone, and G. Pesole CSTminer: a web tool for the identification of coding and noncoding conserved sequence tags through cross-species genome comparison Nucleic Acids Res., July 1, 2004; 32(suppl_2): W624 - W627. [Abstract] [Full Text] [PDF]

				Y. Chen and W. Stephan Compensatory evolution of a precursor messenger RNA secondary structure in the Drosophila melanogaster Adh gene PNAS, September 30, 2003; 100(20): 11499 - 11504. [Abstract] [Full Text] [PDF]

				D. B. Carlini and W. Stephan In Vivo Introduction of Unpreferred Synonymous Codons Into the Drosophila Adh Gene Results in Reduced Levels of ADH Protein Genetics, January 1, 2003; 163(1): 239 - 243. [Abstract] [Full Text] [PDF]

				J. F. Baines, Y. Chen, A. Das, and W. Stephan DNA Sequence Variation at a Duplicated Gene: Excess of Replacement Polymorphism and Extensive Haplotype Structure in the Drosophila melanogasterbicoid Region Mol. Biol. Evol., July 1, 2002; 19(7): 989 - 998. [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]

				J. Parsch, C. D. Meiklejohn, and D. L. Hartl Patterns of DNA Sequence Variation Suggest the Recent Action of Positive Selection in the janus-ocnus Region of Drosophila simulans Genetics, October 1, 2001; 159(2): 647 - 657. [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]

				J. Gorodkin, S. L. Stricklin, and G. D. Stormo Discovering common stem-loop motifs in unaligned RNA sequences Nucleic Acids Res., May 15, 2001; 29(10): 2135 - 2144. [Abstract] [Full Text] [PDF]