Comparative analysis of secondary structure of insect mitochondrial small subunit ribosomal RNA using maximum weighted matching

Comparative analysis is the preferred method of inferring RNA secondary structure, but its use requires considerable expertise and manual effort. As the importance of secondary structure for accurate sequence alignment and phylogenetic analysis becomes increasingly realised, the need for secondary structure models for diverse taxonomic groups becomes more pressing. The number of available structures bears little relation to the relative diversity or importance of the different taxonomic groups. Insects, for example, comprise the largest group of animals and yet are very poorly represented in secondary structure databases. This paper explores the utility of maximum weighted matching (MWM) to help automate the process of comparative analysis by inferring secondary structure for insect mitochondrial small subunit (12S) rRNA sequences. By combining information on correlated changes in substitutions and helix dot plots, MWM can rapidly generate plausible models of secondary structure. These models can be further refined using standard comparative techniques. This paper presents a secondary structure model for insect 12S rRNA based on an alignment of 225 insect sequences and an alignment for 16 exemplar insect sequences. This alignment is used as a template for a web server that automatically generates secondary structures for insect sequences.     

Mutational patterns in RNA secondary structure evolution examined in three RNA families

The goal of this work was to study mutational patterns in the evolution of RNA secondary structure. We analyzed bacterial tmRNA, RNaseP and eukaryotic telomerase RNA secondary structures, mapping structural variability onto phylogenetic trees constructed primarily from rRNA sequences. We found that secondary structures evolve both by whole stem insertion/deletion, and by mutations that create or disrupt stem base pairing. We analyzed the evolution of stem lengths and constructed substitution matrices describing the changes responsible for the variation in the RNA stem length. In addition, we used principal component analysis of the stem length data to determine the most variable stems in different families of RNA. This data provides new insights into the evolution of RNA secondary structures and patterns of variation in the lengths of double helical regions of RNA molecules. Our findings will facilitate design of improved mutational models for RNA structure evolution.     

Large Subunit Mitochondrial rRNA Secondary Structures and Site-Specific Rate Variation in Two Lizard Lineages

A phylogenetic-comparative approach was used to assess and refine existing secondary structure models for a frequently studied region of the mitochondrial encoded large subunit (16S) rRNA in two large lizard lineages within the Scincomorpha, namely the Scincidae and the Lacertidae. Potential pairings and mutual information were analyzed to identify site interactions present within each lineage and provide consensus secondary structures. Many of the interactions proposed by previous models were supported, but several refinements were possible. The consensus structures allowed a detailed analysis of rRNA sequence evolution. Phylogenetic trees were inferred from Bayesian analyses of all sites, and the topologies used for maximum likelihood estimation of sequence evolution parameters. Assigning gamma-distributed relative rate categories to all interacting sites that were homologous between lineages revealed substantial differences between helices. In both lineages, sites within helix G2 were mostly conserved, while those within helix E18 evolved rapidly. Clear evidence of substantial site-specific rate variation (covarion-like evolution) was also detected, although this was not strongly associated with specific helices. This study, in conjunction with comparable findings on different, higher-level taxa, supports the ubiquitous nature of site-specific rate variation in this gene and justifies the incorporation of covarion models in phylogenetic inference.Reviewing Editor: Dr. Yves Van de Peer     

Sequence and secondary structure variation in the Gyrodactylus (Platyhelminthes: Monogenea) ribosomal RNA gene array

Nucleotide sequences were determined for the rRNA internal transcribed spacers 1 and 2 (ITS1 and 2) and the 5' terminus of the large subunit rRNA in selected Gyrodactylus species. Examination of primary sequence variation and secondary structure models in ITS2 and variable region V4 of the small subunit rRNA revealed that structure was largely conserved despite significant variation in sequence. ITS1 sequences were highly variable, and models of structure were unreliable but, despite this, show some resemblance to structures predicted in Digenea. ITS2 models demonstrated binding of the 3' end of 5.8S rRNA to the 5' end of the large subunit rRNA and enabled the termini of these genes to be defined with greater confidence than previously. The structure model shown here may prove useful in future phylogenetic analyses.     

Human ribosomal RNA variants from a single individual and their expression in different tissues.

We have investigated the extent of sequence variation in human ribosomal RNA (rRNA) genes and the expression of specific rRNA gene variants in different tissues of an individual. Focusing on the fifth variable region (V5; nt 2065-2244) of the 28S rRNA gene, we find that sequence differences between rRNA genes of a single individual are characterized by differences in number of repeats of simple sequences at four specific sites. These data support and extend previous findings which show similar V5 sequence variation in rRNA genes from a group of individuals. We performed experiments to determine if there is differential gene expression within the rRNA multigene family. From the analysis of data of six variant V5 probes protected from RNase digestion by rRNAs isolated from different tissues of the individual, we conclude that each variant rRNA is present in a similar proportion in these tissues, whereas the actual contributions of variants differ, their relative proportion is maintained from tissue to tissue in an individual. We favor the explanation of a gene dosage effect over that of a regulated gene effect to account for this pattern of rRNA gene expression. In addition, computer generated secondary structure models of each V5 clone structure predict the same three helix structure with the regions of sequence variation contained in one stem-loop structure.     

Evolution of large subunit rRNA structure. The 3' terminal domain contains elements of secondary structure specific to major phylogenetic groups

Refined secondary structure models supported by phylogenetic evidence have been derived for the 3' terminal domain of large subunit rRNA (the region that exists as a separate 4.5 S molecular entity in chloroplast ribosomes) through a comparative analysis of all the pro- and eukaryotic sequences at present available. While several universally conserved features of secondary structure are found, a few diversified structural elements are also detected which are specific to one of the primary kingdoms, eubacteria, archaebacteria, or eukaryotes. Remarkably, some appear to be selectively preserved during the evolution of the primary kindgom, suggesting they represent functionally important structures. Thus, although the role of this 3' terminal domain in ribosomal function still remains unknown, its mode of sequence variation clearly points to a significant diversification of its function among the primary kindgoms.     

Proposed secondary structure of eukaryote specific expansion segment 15 in 28S rRNA from mice, rats, and rabbits

The expansion segments in eukaryotic ribosomal RNAs are additional RNA sequences not found in the RNA core common to both prokaryotes and eukaryotes. These regions show large species-dependent variations in sequence and size. This makes it difficult to create secondary structure models for the expansion segments exclusively based on phylogenetic sequence comparison. Here we have used a combination of experimental data and computational methods to generate secondary structure models for expansion segment 15 in 28S rRNA in mice, rats, and rabbits. The experimental data were collected using the structure sensitive reagents DMS, CMCT, kethoxal, micrococcal nuclease, RNase T(1), RNase CL3, RNase V(1), and lead(II) acetate. ES15 was folded with the computer program RNAStructure 3.5 using modification data and phylogenetic similarities between different ES15 sequences. This program uses energy minimization to find the most stable secondary structure of an RNA sequence. The presented secondary structure models include several common structural motifs, but they also have characteristics unique to each organism. Overall, the secondary structure models showed indications of an energetically stable but dynamic structure, easily accessible from the solution by the modification reagents, suggesting that the expansion segment is located on the ribosomal surface.     

Natural selection is not required to explain universal compositional patterns in rRNA secondary structure categories

We have encountered an unexpected property of rRNA secondary structures that may generalize to all RNAs. Analysis of 8892 ribosomal RNA sequences and structures from a wide range of species revealed unexpected universal compositional trends. First, different categories of rRNA secondary structure (stems, loops, bulges, and junctions) have distinct, characteristic base compositions. Second, the observed patterns of variation are similar among sequences from large and small rRNA subunits and all domains of life, despite extensive evolutionary divergence. Surprisingly, these differences do not seem to be related to selection for different compositions in different structural categories, but rather relate to the overall composition of the molecule: Randomized RNAs with no evolutionary history show the same structure-dependent compositional biases as rRNAs. These compositional trends may improve the accuracy of RNA secondary structure prediction, because they allow us to compare predicted structures against known compositional preferences. They also suggest caution in interpreting differences in the rate of change of the GC content in different parts of the molecule as evidence of differential selection.     

Sequence and secondary structure of mouse 28S rRNA 5''terminal domain. Organisation of the 5.8S-28S rRNA complex.

We present the sequence of the 5' terminal 585 nucleotides of mouse 28S rRNA as inferred from the DNA sequence of a cloned gene fragment. The comparison of mouse 28S rRNA sequence with its yeast homolog, the only known complete sequence of eukaryotic nucleus-encoded large rRNA (see ref. 1, 2) reveals the strong conservation of two large stretches which are interspersed with completely divergent sequences. These two blocks of homology span the two segments which have been recently proposed to participate directly in the 5.8S-large rRNA complex in yeast (see ref. 1) through base-pairing with both termini of 5.8S rRNA. The validity of the proposed structural model for 5.8S-28S rRNA complex in eukaryotes is strongly supported by comparative analysis of mouse and yeast sequences: despite a number of mutations in 28S and 5.8S rRNA sequences in interacting regions, the secondary structure that can be proposed for mouse complex is perfectly identical with yeast's, with all the 41 base-pairings between the two molecules maintained through 11 pairs of compensatory base changes. The other regions of the mouse 28S rRNA 5'terminal domain, which have extensively diverged in primary sequence, can nevertheless be folded in a secondary structure pattern highly reminiscent of their yeast' homolog. A minor revision is proposed for mouse 5.8S rRNA sequence.     

A new method to predict the consensus secondary structure of a set of unaligned RNA sequences

MOTIVATION: To predict the consensus secondary structure, possibly including pseudoknots, of a set of RNA unaligned sequences. RESULTS: We have designed a method based on a new representation of any RNA secondary structure as a set of structural relationships between the helices of the structure. We refer to this representation as a structural pattern. In a first step, we use thermodynamic parameters to select, for each sequence, the best secondary structures according to energy minimization and we represent each of them using its corresponding structural pattern. In a second step, we search for the repeated structural patterns, i.e. the largest structural patterns that occur in at least one sequence, i.e. included in at least one of the structural patterns associated to each sequence. Thanks to an efficient encoding of structural patterns, this search comes down to identifying the largest repeated word suffixes in a dictionary. In a third step, we compute the plausibility of each repeated structural pattern by checking if it occurs more frequently in the studied sequences than in random RNA sequences. We then suppose that the consensus secondary structure corresponds to the repeated structural pattern that displays the highest plausibility. We present several experiments concerning tRNA, fragments of 16S rRNA and 10Sa RNA (including pseudoknots); in each of them, we found the putative consensus secondary structure.     

Consideration of RNA secondary structure significantly improves likelihood-based estimates of phylogeny: examples from the bilateria

Sequences from ribosomal RNA (rRNA) genes have made a huge contribution to our current understanding of metazoan phylogeny and indeed the phylogeny of all of life. That said, some parts of this rRNA-based phylogeny remain unresolved. One approach to increase the resolution of these trees would be to use more appropriate models of sequence evolution in phylogenetic analysis. RNAs transcribed from rRNA genes have a complex secondary structure mediated by base pairing between sometimes distant regions of the rRNA molecule. The pairing between the stem nucleotides has important consequences for their evolution which differs from that of unpaired loop nucleotides. These differences in evolution should ideally be accounted for when using rRNA sequences for phylogeny estimation. We use a novel permutation approach to demonstrate the significant superiority of models of sequence evolution that allow stem and loop regions to evolve according to separate models and, in common with previous studies, we show that 16-state models that take base pairing of stems into account are significantly better than simpler, 4-state, single-nucleotide models. One of these 16-state models has been applied to the phylogeny of the Bilateria using small subunit rRNA (SSU) sequences. Our optimal tree largely echoes previous results based on SSU in particular supporting the tripartite Bilaterian tree of deuterostomes, lophotrochozoans, and ecdysozoans. There are also a number of differences, however, perhaps most important of which is the observation of a clade consisting of the gastrotrichs plus platyheminthes that is basal to all other lophotrochozoan taxa. Use of 16-state models also appears to reduce the Bayesian support given to certain biologically improbable groups found using standard 4-state models.     

Sequence determinants of quaternary structure in lumazine synthase

Riboflavin, an essential cofactor for all organisms, is biosynthesized in plants, fungi and microorganisms. The penultimate step in the pathway is catalyzed by the enzyme lumazine synthase. One of the most distinctive characteristics of this enzyme is that it is found in different species in two different quaternary structures, pentameric and icosahedral, built from practically the same structural monomeric unit. In fact, the icosahedral structure is best described as a capsid of twelve pentamers. Despite this noticeable difference, the active sites are virtually identical in all structurally studied members. Furthermore, the main regions involved in the catalysis are located at the interface between adjacent subunits in the pentamer. Thus, the two quaternary forms of the enzyme must meet similar structural requirements to achieve their function, but, at the same time, they should differ in the sequence traits responsible for the different quaternary structures observed. Here, we present a combined analysis that includes sequence-structure and evolutionary studies to find the sequence determinants of the different quaternary assemblies of this enzyme. A data set containing 86 sequences of the lumazine synthase family was recovered by sequence similarity searches. Seven of them had resolved three-dimensional structures. A subsequent phylogenetic reconstruction by maximum parsimony (MP) allowed division of the total set into two clusters in accord with their quaternary structure. The comparison between the patterns of three-dimensional contacts derived from the known three-dimensional structures and variation in sequence conservation revealed a significant shift in structural constraints of certain positions. Also, to explore the changes in functional constraints between the two groups, site-specific evolutionary rate shifts were analyzed. We found that the positions involved in icosahedral contacts suffer a larger increase in constraints than the rest. We found eight sequence sites that would be the most important icosahedral sequence determinants. We discuss our results and compare them with previous work. These findings should contribute to refinement of the current structural data, to the design of assays that explore the role of these positions, to the structural characterization of new sequences, and to initiation of a study of the underlying evolutionary mechanisms.     

Selection on the structural stability of a ribosomal RNA expansion segment in Daphnia obtusa

The high rate of sequence divergence in nuclear ribosomal RNA (rRNA) expansion segments offers a unique opportunity to study the importance of natural selection in their evolution. To this end, we polymerase chain reaction amplified and cloned a 589-nt fragment of the 18S rRNA gene containing expansion segments 43/e1 and 43/e4 from six individual Daphnia obtusa from four populations. We screened 2,588 clones using single-stranded conformation polymorphism analysis and identified 103 unique haplotype sequences. We detected two pairs of indel sites in segment 43/e4 that complement each other when the secondary structure of the linear sequence is formed. Seven of the 12 observed combinations of length variants at these four sites (haplotypes) are shared between individuals from different populations, which may suggest that some of the length variation was present in their common ancestor. Haplotypes with uncompensated indels were only observed at low frequencies, while compensated indel haplotypes were found at a wide range of frequencies, supporting the hypothesis that the energetic stability of expansion segments is a trait under natural selection. In addition, there was strong linkage disequilibrium between the four complementary indel sites, particularly those that pair with one another in the secondary structure. Despite selection against unpaired bulges at these four indel sites, some nucleotides that form unpaired bulges are highly conserved in segment 43/e4, indicating that they are under a different selective constraint, possibly due to their role in higher level structural interactions.     

Probing the secondary structure of expansion segment ES6 in 18S ribosomal RNA

Expansion segment ES6 in 18S ribosomal RNA is, unlike many other expansion segments, present in all eukaryotes. The available data suggest that ES6 is located on the surface of the small ribosomal subunit. Here we have analyzed the secondary structure of the complete ES6 sequence in intact ribosomes from three eukaryotes, wheat, yeast, and mouse, representing different eukaryotic kingdoms. The availability of the ES6 sequence for modification and cleavage by structure sensitive chemicals and enzymatic reagents was analyzed by primer extension and gel electrophoresis on an ABI 377 automated DNA sequencer. The experimental results were used to restrict the number of possible secondary structure models of ES6 generated by the folding software MFOLD. The modification data obtained from the three experimental organisms were very similar despite the sequence variation. Consequently, similar secondary structure models were obtained for the ES6 sequence in wheat, yeast, and mouse ribosomes. A comparison of sequence data from more than 6000 eukaryotes showed that similar structural elements could also be formed in other organisms. The comparative analysis also showed that the extent of compensatory base changes in the suggested helices was low. The in situ structure analysis was complemented by a secondary structure analysis of wheat ES6 transcribed and folded in vitro. The obtained modification data indicate that the secondary structure of the in vitro transcribed sequence differs from that observed in the intact ribosome. These results suggest that chaperones, ribosomal proteins, and/or tertiary rRNA interactions could be involved in the in vivo folding of ES6.     

Secondary structure of mouse 28S rRNA and general model for the folding of the large rRNA in eukaryotes.

We present a secondary structure model for the entire sequence of mouse 28S rRNA (1) which is based on an extensive comparative analysis of the available eukaryotic sequences, i.e. yeast (2, 3), Physarum polycephalum (4), Xenopus laevis (5) and rat (6). It has been derived with close reference to the models previously proposed for yeast 26S rRNA (2) and for prokaryotic 23S rRNA (7-9). Examination of the recently published eukaryotic sequences confirms that all pro- and eukaryotic large rRNAs share a largely conserved secondary structure core, as already apparent from the previous analysis of yeast 26S rRNA (2). These new comparative data confirm most features of the yeast model (2). They also provide the basis for a few modifications and for new proposals which extend the boundaries of the common structural core (now representing about 85% of E. coli 23S rRNA length) and bring new insights for tracing the structural evolution, in higher eukaryotes, of the domains which have no prokaryotic equivalent and are inserted at specific locations within the common structural core of the large subunit rRNA.     

Phylogenetic analysis of molluscan mitochondrial LSU rDNA sequences and secondary structures

Mollusks are an extraordinarily diverse group of animals with an estimated 200,000 species, second only to the phylum Arthropoda. We conducted a comparative analysis of complete mitochondrial ribosomal large subunit sequences (LSU) of a chiton, two bivalves, six gastropods, and a cephalopod. In addition, we determined secondary structure models for each of them. Comparative analyses of nucleotide variation revealed substantial length variation among the taxa, with stylommatophoran gastropods possessing the shortest lengths. Phylogenetic analyses of the nucleotide sequence data supported the monophyly of Albinaria, Euhadra herklotsi + Cepaea nemoralis, Stylommatophora, Cerithioidea, and when only transversions are included, the Bivalvia. The phylogenetic limits of the mitochondrial LSU rRNA gene within mollusks appear to be up to 400 million years, although this estimate will have to be tested further with additional taxa. Our most novel finding was the discovery of phylogenetic signal in the secondary structure of rRNA of mollusks. The absence of entire stem/loop structures in Domains II, III, and V can be viewed as three shared derived characters uniting the stylommatophoran gastropods. The absence of the aforementioned stem/loop structure explains much of the observed length variation of the mitochondrial LSU rRNA found within mollusks. The distribution of these unique secondary structure characters within mollusks should be examined.     

rRNA二级结构序列用于真菌系统学研究的方法初探

本文首次利用核酸二级结构特征代替核酸碱基作为探讨类群之间亲缘关系的信号,构建了基于结构特征的子囊菌部分类群的系统进化树。该方法以S(规范的碱基对),Q(不规范的碱基对),I(单链),B(侧环),M(多分枝环)和H(发卡结构)为代码将二级结构特征区分为6种不同的亚结构类型,然后将二级结构特征转换为结构序列,并进行结构序列分析。该方法使rRNA不只局限于碱基比较,拓展了其应用范围,为揭示分子的功能与进化的关系提供了线索。结果表明,结构序列分析可用于子囊菌的系统学研究;相对于核酸序列分析,结构分析的结果似乎更加清晰地体现子囊果的演化过程。     

Evolution of large-subunit rRNA structure. The diversification of divergent D3 domain among major phylogenetic groups

During evolution, the potential for sequence (and length) variation of large-subunit rRNA has been mostly restricted over 12 divergent domains (termed D1-D12) interspersed along the molecule. Here, we have focused our attention onto the D3 divergent domain, through a detailed analysis of its pattern of variation in the phylogeny, both in terms of primary and secondary structures. We have systematically compared all the procaryotic and eucaryotic sequences published so far (i.e. 36 species), together with a series of 10 additional eucaryotic specimens, which were determined by direct RNA sequencing. Secondary structures supported by comparative evidence have been derived for archaebacteria, eubacteria and eucaryotes respectively, which shows that the D3 domain contains a subset of universally conserved structural features interspersed with four variable subdomains. Within the four portions where a structural diversification has taken place, elementary structures specific of large phylogenetic groups can be identified. Remarkably such diversified structures appear to be preserved despite sequence divergence, suggesting they correspond to functionally important structures. Accordingly, the mode of sequence variation of the D3 domain suggests this region of the molecule may encode elementary functions of rRNA which could have significantly diversified during the evolution of the major groups of organisms.     

Intra-Genomic Variation in the Ribosomal Repeats of Nematodes

Ribosomal loci represent a major tool for investigating environmental diversity and community structure via high-throughput marker gene studies of eukaryotes (e.g. 18S rRNA). Since the estimation of species’ abundance is a major goal of environmental studies (by counting numbers of sequences), understanding the patterns of rRNA copy number across species will be critical for informing such high-throughput approaches. Such knowledge is critical, given that ribosomal RNA genes exist within multi-copy repeated arrays in a genome. Here we measured the repeat copy number for six nematode species by mapping the sequences from whole genome shotgun libraries against reference sequences for their rRNA repeat. This revealed a 6-fold variation in repeat copy number amongst taxa investigated, with levels of intragenomic variation ranging from 56 to 323 copies of the rRNA array. By applying the same approach to four C. elegans mutation accumulation lines propagated by repeated bottlenecking for an average of ~400 generations, we find on average a 2-fold increase in repeat copy number (rate of increase in rRNA estimated at 0.0285-0.3414 copies per generation), suggesting that rRNA repeat copy number is subject to selection. Within each Caenorhabditis species, the majority of intragenomic variation found across the rRNA repeat was observed within gene regions (18S, 28S, 5.8S), suggesting that such intragenomic variation is not a product of selection for rRNA coding function. We find that the dramatic variation in repeat copy number among these six nematode genomes would limit the use of rRNA in estimates of organismal abundance. In addition, the unique pattern of variation within a single genome was uncorrelated with patterns of divergence between species, reflecting a strong signature of natural selection for rRNA function. A better understanding of the factors that control or affect copy number in these arrays, as well as their rates and patterns of evolution, will be critical for informing estimates of global biodiversity.     

Nucleotide sequences of Acanthamoeba castellanii 5S and 5.8S ribosomal ribonucleic acids: phylogenetic and comparative structural analyses.

Sequences of 5S and 5.8S rRNAs of the amoeboid protist Acanthamoeba castellanii have been determined by gel sequencing of terminally-labeled RNAs which were partially degraded with chemical reagents or ribonucleases. The sequence of the 5S rRNA is (formula, see text). This sequence is compared to eukaryotic 5S rRNA sequences previously published and fitted to a secondary structure model which incorporates features of several previously proposed models. All reported eukaryotic 5S rRNAs fit this model. The sequence of the 5.8S rRNA is (formula, see text). This sequence does not fit parts of existing secondary structure models for 5.8S rRNA, and we question the significance of such models.