How slippage-derived sequences are incorporated into rRNA variable-region secondary structure: implications for phylogeny reconstruction.

We analyzed the type and frequency of mutational changes in hypervariable rRNA regions, using the highly length-variable region V4 of the small subunit rRNA locus of tiger beetles (Cicindelidae) as an example. Phylogenetic analysis of indels in closely related species showed that (1) most indels are single nucleotides (usually A or T and sometimes G) or di-nucleotides of A and T. These occur at numerous foci, and they exhibit a strong bias for duplication of 5' single and di-nucleotide motifs but not 3' motifs. (2) Insertions/deletions in stem-forming regions affected paired and unpaired bases with about equal frequency but they did not disrupt the secondary structure. (3) Recurring mutations involving short repeats of the same bases caused parallel evolution of similar sequence motifs in the rRNA of different lineages. The observed types of change are consistent with the propostion that slippage is the main mutational mechanism. Slippage-derived sequences tend to be self-complementary, and therefore the stem-loop structure could be self-organizing as a consequence of the underlying mutational mechanism. Thus, the secondary structure in the cicindelid V4 region may be conserved due to the dynamics of the mutational mechanism rather than to functional constraints. These processes may also have a tendency to produce similar primary sequences irrespective of phylogenetic associations. The findings have implications for sequence alignment in phylogenetic analysis and should caution against the use of secondary structure to improve the determination of positional homology in hypervariable regions.

[1]  E. Herniou,et al.  Acoel flatworms: earliest extant bilaterian Metazoans, not members of Platyhelminthes. , 1999, Science.

[2]  John M. Hancock,et al.  High sequence turnover in the regulatory regions of the developmental gene hunchback in insects. , 1999, Molecular biology and evolution.

[3]  M. Zuker,et al.  Using reliability information to annotate RNA secondary structures. , 1998, RNA.

[4]  John M. Hancock,et al.  Modelling the secondary structures of slippage-prone hypervariable RNA regions: the example of the tiger beetle 18S rRNA variable region V4. , 1998, Nucleic acids research.

[5]  M. Whiting.  Long-branch distraction and the Strepsiptera. , 1998, Systematic biology.

[6]  R. Raff,et al.  Evidence for a clade of nematodes, arthropods and other moulting animals , 1997, Nature.

[7]  D A Morrison,et al.  Effects of nucleotide sequence alignment on phylogeny estimation: a case study of 18S rDNAs of apicomplexa. , 1997, Molecular biology and evolution.

[8]  W C Wheeler,et al.  The Strepsiptera problem: phylogeny of the holometabolous insect orders inferred from 18S and 28S ribosomal DNA sequences and morphology. , 1997, Systematic biology.

[9]  John M. Hancock,et al.  Phylogenetic analysis of slippage-like sequence variation in the V4 rRNA expansion segment in tiger beetles (Cicindelidae). , 1997, Molecular biology and evolution.

[10]  A. Vogler,et al.  A molecular phylogeny of the tiger beetles (Cicindelidae): congruence of mitochondrial and nuclear rDNA data sets. , 1996, Molecular phylogenetics and evolution.

[11]  V. Schawaroch,et al.  THREE STEPS OF HOMOLOGY ASSESSMENT , 1996, Cladistics : the international journal of the Willi Hennig Society.

[12]  D. Frost,et al.  Molecular homology assessment and phylogeny in the lizard family opluridae (Squamata: Iguania). , 1996, Molecular phylogenetics and evolution.

[13]  D. Penny,et al.  Conserved sequence motifs, alignment, and secondary structure for the third domain of animal 12S rRNA. , 1996, Molecular biology and evolution.

[14]  K. Kjer,et al.  Use of rRNA secondary structure in phylogenetic studies to identify homologous positions: an example of alignment and data presentation from the frogs. , 1995, Molecular phylogenetics and evolution.

[15]  M. Zuker,et al.  "Well-determined" regions in RNA secondary structure prediction: analysis of small subunit ribosomal RNA. , 1995, Nucleic acids research.

[16]  R. Gutell,et al.  Lessons from an evolving rRNA: 16S and 23S rRNA structures from a comparative perspective. , 1994, Microbiological reviews.

[17]  M. Zuker.  Prediction of RNA secondary structure by energy minimization. , 1994, Methods in molecular biology.

[18]  W. Brown,et al.  Rates and patterns of base change in the small subunit ribosomal RNA gene. , 1993, Genetics.

[19]  D. Hillis,et al.  Ribosomal DNA: Molecular Evolution and Phylogenetic Inference , 1991, The Quarterly Review of Biology.

[20]  Mário C. C. Pinna.  CONCEPTS AND TESTS OF HOMOLOGY IN THE CLADISTIC PARADIGM , 1991 .

[21]  D. Turner,et al.  A comparison of optimal and suboptimal RNA secondary structures predicted by free energy minimization with structures determined by phylogenetic comparison. , 1991, Nucleic acids research.

[22]  John M. Hancock,et al.  'Compensatory slippage' in the evolution of ribosomal RNA genes. , 1990, Nucleic acids research.

[23]  D. Turner,et al.  Predicting optimal and suboptimal secondary structure for RNA. , 1990, Methods in enzymology.

[24]  M. Zuker.  On finding all suboptimal foldings of an RNA molecule. , 1989, Science.

[25]  John M. Hancock,et al.  Molecular coevolution among cryptically simple expansion segments of eukaryotic 26S/28S rRNAs. , 1988, Molecular biology and evolution.

[26]  W. Wheeler,et al.  Paired sequence difference in ribosomal RNAs: evolutionary and phylogenetic implications. , 1988, Molecular biology and evolution.

[27]  G. Gutman,et al.  Slipped-strand mispairing: a major mechanism for DNA sequence evolution. , 1987, Molecular biology and evolution.

[28]  S. Gerbi.  The evolution of eukaryotic ribosomal DNA. , 1986, Bio Systems.

[29]  G. Dover,et al.  Molecular drive: a cohesive mode of species evolution , 1982, Nature.