An empirical analysis of mt 16S rRNA covarion-like evolution in insects: site-specific rate variation is clustered and frequently detected.

The structural and functional analysis of rRNA molecules has attracted considerable scientific interest. Empirical studies have demonstrated that sequence variation is not directly translated into modifications of rRNA secondary structure. Obviously, the maintenance of secondary structure and sequence variation are in part governed by different selection regimes. The nature of those selection regimes still remains quite elusive. The analysis of individual bacterial models cannot adequately explore this topic. Therefore, we used primary sequence data and secondary structures of a mitochondrial 16S rRNA fragment of 558 insect species from 15 monophyletic groups to study patterns of sequence variation, and variation of secondary structure. Using simulation studies to establish significance levels of change, we found that despite conservation of secondary structure, the location of sequence variation within the conserved rRNA structure changes significantly between groups of insects. Despite our conservative estimation procedure we found significant site-specific rate changes at 56 sites out of 184. Additionally, site-specific rate variation is somewhat clustered in certain helices. Both results confirm what has been predicted from an application of non-stationary maximum likelihood models to rRNA sequences. Clearly, constraints on sequence variation evolve and leave footprints in the form of evolutionary plasticity in rRNA sequences. Here, we show that a better understanding of the evolution of rRNA sequences can be obtained by integrating both phylogenetic and structural information.