Secondary structure of two regions in expansion segments ES3 and ES6 with the potential of forming a tertiary interaction in eukaryotic 40S ribosomal subunits

The 18S rRNA of the small eukaryotic ribosomal subunit contains several expansion segments. Electron microscopy data indicate that two of the largest expansion segments are juxtaposed in intact 40S subunits, and data from phylogenetic sequence comparisons indicate that these two expansion segments contain complementary sequences that could form a direct tertiary interaction on the ribosome. We have investigated the secondary structure of the two expansion segments in the region around the putative tertiary interaction. Ribosomes from yeast, wheat, and mouse-three organisms representing separate eukaryotic kingdoms-were isolated, and the structure of ES3 and part of the ES6 region were analyzed using the single-strand-specific chemical reagents CMCT and DMS and the double-strand-specific ribonuclease V1. The modification patterns were analyzed by primer extension and gel electrophoresis on an ABI 377 automated DNA sequencer. The investigated sequences were relatively exposed to chemical and enzymatic modification. This is in line with their indicated location on the surface at the solvent side of the subunit. The complementary ES3 and ES6 sequences were clearly inaccessible to single-strand modification, but available for cleavage by double-strand-specific RNase V1. The results are compatible with a direct helical interaction between bases in ES3 and ES6. Almost identical results were obtained with ribosomes from the three organisms investigated.     

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.     

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.     

Comparison of fungal 80 S ribosomes by cryo-EM reveals diversity in structure and conformation of rRNA expansion segments

Compared to the prokaryotic 70 S ribosome, the eukaryotic 80 S ribosome contains additional ribosomal proteins and extra segments of rRNA, referred to as rRNA expansion segments (ES). These eukaryotic-specific rRNA ES are mainly on the periphery of the 80 S ribosome, as revealed by cryo-electron microscopy (cryo-EM) studies, but their precise function is not known. To address the question of whether the rRNA ES are structurally conserved among 80 S ribosomes of different fungi we performed cryo-electron microscopy on 80 S ribosomes from the thermophilic fungus Thermomyces lanuginosus and compared it to the Saccharomyces cerevisiae 80 S ribosome. Our analysis reveals general structural conservation of the rRNA expansion segments but also changes in ES27 and ES7/39, as well as the absence of a tertiary interaction between ES3 and ES6 in T. lanuginosus. The differences provide a hint on the role of rRNA ES in regulating translation. Furthermore, we show that the stalk region and interactions with elongation factor 2 (eEF2) are different in T. lanuginosus, exhibiting a more extensive contact with domain I of eEF2.     

Molecular coevolution among cryptically simple expansion segments of eukaryotic 26S/28S rRNAs

The set of "expansion segments" of any eukaryotic 26S/28S ribosomal RNA (rRNA) gene is responsible for the bulk of the difference in length between the prokaryotic 23S rRNA gene and the eukaryotic 26S/28S rRNA gene. The expansion segments are also responsible for interspecific fluctuations in length during eukaryotic evolution. They show a consistent bias in base composition in any species; for example, they are AT rich in Drosophila melanogaster and GC rich in vertebrate species. Dot-matrix comparisons of sets of expansion segments reveal high similarities between members of a set within any 28S rRNA gene of a species, in contrast to the little or spurious similarity that exists between sets of expansion segments from distantly related species. Similarities among members of a set of expansion segments within any 28S rRNA gene cannot be accounted for by their base-compositional bias alone. In contrast, no significant similarity exists within a set of "core" segments (regions between expansion segments) of any 28S rRNA gene, although core segments are conserved between species. The set of expansion segments of a 26S/28S gene is coevolving as a unit in each species, at the same time as the family of 28S rRNA genes, as a whole, is undergoing continual homogenization, making all sets of expansion segments from all ribosomal DNA (rDNA) arrays in a species similar in sequence. Analysis of DNA simplicity of 26S/28S rRNA genes shows a direct correlation between significantly high relative simplicity factors (RSFs) and sequence similarity among a set of expansion segments. A similar correlation exists between RSF values, overall rDNA lengths, and the lengths of individual expansion segments. Such correlations suggest that most length fluctuations reflect the gain and loss of simple sequence motifs by slippage-like mechanisms. We discuss the molecular coevolution of expansion segments, which takes place against a background of slippage-like and unequal crossing-over mechanisms of turnover that are responsible for the accumulation of interspecific differences in rDNA sequences.      

Analysis of the secondary structure of expansion segment 39 in ribosomes from fungi, plants and mammals

The structure of expansion segment 39, ES39, in eukaryotic 23 S-like ribosomal RNA was analysed using a combination of chemical and enzymic reagents. Ribosomes were isolated from yeast, wheat, mouse, rat and rabbit, five organisms representing three different eukaryotic kingdoms. The isolated ribosomes were treated with structure-sensitive chemical and enzymic reagents and the modification patterns analysed by primer extension and gel electrophoresis on an ABI 377 automated DNA sequencer. The expansion segment was relatively accessible to modification by both enzymic and chemical probes, suggesting that ES39 was exposed on the surface of the ribosomes. The collected modification data were used in secondary structure modelling of the expansion segment. Despite considerable variation in both sequence and length between organisms from different kingdoms, the structure analysis of the expansion segment gave rise to structural fingerprints that allowed identification of homologous structures in ES39 from fungi, plants and mammals. The homologous structures formed an initial helix and an invariant hairpin connected to the initial helix via a long single-stranded loop. The remaining part of the ES39 sequences accounted for most of the length variation seen between the analysed species. This part could form additional, albeit less similar, hairpins. A comparison of ES39 sequences from other fungi, plants and mammals showed that identical structures could be formed in these organisms.     

The evolution of eukaryotic ribosomal DNA

Mutations occur randomly throughout the ribosomal DNA (rDNA) sequence. Molecular drive (unequal crossing-over, gene conversion, and transposition) spreads these variations through the multiple copies of rDNA. Forces of selection act upon the variants to favor and fix them or disfavor and eliminate them. Selection has not permitted changes in regions within rRNA vital for its function; these sequences are evolutionarily conserved between diverse species. Possible functions for some of these conserved sequences are discussed. The secondary structure of rRNA is also highly conserved during evolution. However, eukaryotic rRNA is larger than prokaryotic rRNA due to blocks of "expansion segments". Arguments are put forward that expansion segments might not play any functional role. Other examples are reviewed of rDNA sequence insertion or deletion, including introns and the internal transcribed spacer 2.     

G-Quadruplexes in Human Ribosomal RNA

rRNA is the single most abundant polymer in most cells. Mammalian rRNAs are nearly twice as large as those of prokaryotes. Differences in rRNA size are due to expansion segments, which contain extended tentacles in metazoans. Here we show that the terminus of an rRNA tentacle of Homo sapiens contains 10 tandem G-tracts that form highly stable G-quadruplexes in vitro. We characterized rRNA of the H. sapiens large ribosomal subunit by computation, circular dichroism, UV melting, fluorescent probes, nuclease accessibility, electrophoretic mobility shifts, and blotting. We investigated Expansion Segment 7 (ES7), oligomers derived from ES7, intact 28S rRNA, 80S ribosomes, and polysomes. We used mass spectrometry to identify proteins that bind to rRNA G-quadruplexes in cell lysates. These proteins include helicases (DDX3, CNBP, DDX21, DDX17) and heterogeneous nuclear ribonucleoproteins. Finally, by multiple sequence alignments, we observe that G-quadruplex-forming sequences are a general feature of LSU rRNA of Chordata but not, as far as we can tell, of other species. Chordata ribosomes present polymorphic tentacles with the potential to switch between inter- and intramolecular G-quadruplexes. To our knowledge, G-quadruplexes have not been reported previously in ribosomes.     

The role of expansion segment of human ribosomal protein L35 in nuclear entry, translation activity, and endoplasmic reticulum docking.

The phylogenic alignment of homologous L35 protein suggests that human large subunit ribosomal protein L35 carries a 54 aa eukaryotic expansion segment (ES) at the C-terminal end. Within this ES, the first 25 amino acid residues were found to be essential for the nuclear import of the protein. The last 29 residues of the ES were shown to be uninvolved in the ribosome's structural and translational functions, although this region proved to be one of the contact sites for ribosomal docking to endoplasmic reticulum, as evident from the results of an in vivo recombinant ribosome analysis.     

Examination of protein sequence homologies: III. Ribosomal protein YS25 fromSaccharomyces cerevisiae and its counterparts fromSchizosaccharomyces pombe,rat liver,andEscherichia coli

Summary The sequences of the ribosomal proteins YS25, SP-S28, RL-S21, and Ec-S6, fromSaccharomyces cerevisiae, Schizosaccharomyces pombe, rat liver, andEscherichia coli, respectively, have been examined using a computer program that searches for homologous tertiary structures. Matrices of comparisons among the eukaryotic sequences show that they match each other sequentially without any internal gaps. The average values of the correlation coefficients obtained from the comparison matrices are higher for the first halves of the sequences than for the latter halves. This result suggests that the first halves of the sequences may represent a more important domain than the latter halves. The comparison matrices between the eukaryotic and bacterial sequences of ribosomal proteins, however, do not show sequentially arranged homology, though there are six well-matching segments arranged in different orders in the two types of sequences. This implies that the eukaryotic sequences of the ribosomal protein were reconstituted by two internal transpositions and six deletions of 4–12 residues each from the ancestral sequence during the divergence between bacterial and eukaryotic genes. These findings may give insight into structural and quantitative studies of evolutionary divergence between eukaryotes and prokaryotes.     

A census of protein repeats.

In this study, we analyzed all known protein sequences for repeating amino acid segments. Although duplicated sequence segments occur in 14 % of all proteins, eukaryotic proteins are three times more likely to have internal repeats than prokaryotic proteins. After clustering the repetitive sequence segments into families, we find repeats from eukaryotic proteins have little similarity with prokaryotic repeats, suggesting most repeats arose after the prokaryotic and eukaryotic lineages diverged. Consequently, protein classes with the highest incidence of repetitive sequences perform functions unique to eukaryotes. The frequency distribution of the repeating units shows only weak length dependence, implicating recombination rather than duplex melting or DNA hairpin formation as the limiting mechanism underlying repeat formation. The mechanism favors additional repeats once an initial duplication has been incorporated. Finally, we show that repetitive sequences are favored that contain small and relatively water-soluble residues. We propose that error-prone repeat expansion allows repetitive proteins to evolve more quickly than non-repeat-containing proteins.     

The sequence of 28S ribosomal RNA varies within and between human cell lines.

The primary structure of 28S ribosomal RNA constitutes a conserved core which is similar among most 23S-like rRNAs and expansion segments which occur at specific positions in the sequence. The expansion segments account for most of the size difference between prokaryotic (archaeal and eubacterial) and eukaryotic rRNAs and they exhibit a sequence variation which is unique among rRNAs. We have investigated the sequence variation of one of the expansion segments, V8, by sequencing a total of 111 V8 segments from 9 different human cell lines and tissues and have found 35 different variants. The variation occur mainly at two 'hot spots' which are separated by 170 nucleotides in the primary sequence but are neighbours in the secondary structure. The sequence of V8 segments varies both within and between human cell lines and tissues. The implications for the evolution of the eukaryotic 28S rRNA are discussed together with possible functions of the expansion segments. We also present a secondary structure model for the V8 segment based on comparative sequence analysis and chemical and enzymatic foot printing.     

Trapped in freshwater: the internal anatomy of the entoproct Loxosomatoides sirindhornae

Background

The identification of vast numbers of unknown organisms using DNA sequences becomes more and more important in ecological and biodiversity studies. In this context, a fragment of the mitochondrial cytochrome c oxidase I (COI) gene has been proposed as standard DNA barcoding marker for the identification of organisms. Limitations of the COI barcoding approach can arise from its single-locus identification system, the effect of introgression events, incomplete lineage sorting, numts, heteroplasmy and maternal inheritance of intracellular endosymbionts. Consequently, the analysis of a supplementary nuclear marker system could be advantageous.

Results

We tested the effectiveness of the COI barcoding region and of three nuclear ribosomal expansion segments in discriminating ground beetles of Central Europe, a diverse and well-studied invertebrate taxon. As nuclear markers we determined the 18S rDNA: V4, 18S rDNA: V7 and 28S rDNA: D3 expansion segments for 344 specimens of 75 species. Seventy-three species (97%) of the analysed species could be accurately identified using COI, while the combined approach of all three nuclear markers provided resolution among 71 (95%) of the studied Carabidae.

Conclusion

Our results confirm that the analysed nuclear ribosomal expansion segments in combination constitute a valuable and efficient supplement for classical DNA barcoding to avoid potential pitfalls when only mitochondrial data are being used. We also demonstrate the high potential of COI barcodes for the identification of even closely related carabid species.     

Intra-genomic ribosomal RNA polymorphism and morphological variation in Elphidium macellum suggests inter-specific hybridization in foraminifera

Elphidium macellum is a benthic foraminifer commonly found in the Patagonian fjords. To test whether its highly variable morphotypes are ecophenotypes or different genotypes, we analysed 70 sequences of the SSU rRNA gene from 25 specimens. Unexpectedly, we identified 11 distinct ribotypes, with up to 5 ribotypes co-occurring within the same specimen. The ribotypes differ by varying blocks of sequence located at the end of stem-loop motifs in the three expansion segments specific to foraminifera. These changes, distinct from typical SNPs and indels, directly affect the structure of the expansion segments. Their mosaic distribution suggests that ribotypes originated by recombination of two or more clusters of ribosomal genes. We propose that this expansion segment polymorphism (ESP) could originate from hybridization of morphologically different populations of Patagonian Elphidium. We speculate that the complex geological history of Patagonia enhanced divergence of coastal foraminiferal species and contributed to increasing genetic and morphological variation.     

Secondary structure constraints on the evolution of Drosophila 28 S ribosomal RNA expansion segments

Eukaryotic ribosomal RNA genes contain rapidly evolving regions of unknown function termed expansion segments. We present the comparative analysis of the primary and secondary structure of two expansion segments from the large subunit rRNA gene of ten species of Drosophila and the tsetse fly species Glossina morsitans morsitans. At the primary sequence level, most of the differences observed in the sequences obtained are single base substitutions. This is in marked contrast with observations in vertebrate species in which the insertion or deletion of repetitive motifs, probably generated by a DNA-slippage mechanism, is a major factor in the evolution of these regions. The secondary structure of the two regions, supported by multiple compensatory base changes, is highly conserved between the species examined and supports the existence of a general folding pattern for all eukaryotes. Intriguingly, the evolutionary rate of expansion segments is very slow relative to other genic and non-genic regions of the Drosophila genome. These results suggest that the evolution of expansion segments in the rDNA multigene family is a balance between the homogenization of new mutations by unequal crossing over and a combination of selection against some such mutations per se and selection for subsequent compensatory mutations, in order to maintain a particular RNA secondary structure.     

Coevolution of ribosomal RNA expansion segment 7L and assembly factor Noc2p specializes the ribosome biogenesis pathway between Saccharomyces cerevisiae and Candida albicans

Ribosomes of different species share an evolutionarily conserved core, exhibiting flexible shells formed partially by the addition of species-specific ribosomal RNAs (rRNAs) with largely unexplored functions. In this study, we showed that by swapping the Saccharomyces cerevisiae 25S rRNA genes with non-S. cerevisiae homologs, species-specific rRNA variations caused moderate to severe pre-rRNA processing defects. Specifically, rRNA substitution by the Candida albicans caused severe growth defects and deficient pre-rRNA processing. We observed that such defects could be attributed primarily to variations in expansion segment 7L (ES7^L) and could be restored by an assembly factor Noc2p mutant (Noc2p-K384R). We showed that swapping ES7^L attenuated the incorporation of Noc2p and other proteins (Erb1p, Rrp1p, Rpl6p and Rpl7p) into pre-ribosomes, and this effect could be compensated for by Noc2p-K384R. Furthermore, replacement of Noc2p with ortholog from C. albicans could also enhance the incorporation of Noc2p and the above proteins into pre-ribosomes and consequently restore normal growth. Taken together, our findings help to elucidate the roles played by the species-specific rRNA variations in ribosomal biogenesis and further provide evidence that coevolution of rRNA expansion segments and cognate assembly factors specialized the ribosome biogenesis pathway, providing further insights into the function and evolution of ribosome.     

Complementarity between ferritin H mRNA and 28 S ribosomal RNA

We have found an interesting complementarity in sequences of human ferritin H mRNA and 28 S ribosomal RNA. Immediately upstream of the initiating AUG in the ferritin mRNA is a stretch of 67 nucleotides which contains sequences complementary to several regions in 28 S RNA. One such region can form 55 base pairings with the 5' noncoding region of the ferritin H mRNA. Most of the complementarity is due to repeats of CCG in the ferritin mRNA and GGC in the ribosomal RNA. The regions of complementarity in the 28 S RNA appear to be expansion sequences that have arisen in the evolution of eukaryotic ribosomal RNA. We suggest that interaction of ferritin mRNA and 28 S RNA may function to regulate the stability and/or translatability of ferritin mRNA.     

Ribosomal RNAs are tolerant toward genetic insertions: evolutionary origin of the expansion segments

Ribosomal RNAs (rRNAs), assisted by ribosomal proteins, form the basic structure of the ribosome, and play critical roles in protein synthesis. Compared to prokaryotic ribosomes, eukaryotic ribosomes contain elongated rRNAs with several expansion segments and larger numbers of ribosomal proteins. To investigate architectural evolution and functional capability of rRNAs, we employed a Tn5 transposon system to develop a systematic genetic insertion of an RNA segment 31 nt in length into Escherichia coli rRNAs. From the plasmid library harboring a single rRNA operon containing random insertions, we isolated surviving clones bearing rRNAs with functional insertions that enabled rescue of the E. coli strain (Δ7rrn) in which all chromosomal rRNA operons were depleted. We identified 51 sites with functional insertions, 16 sites in 16S rRNA and 35 sites in 23S rRNA, revealing the architecture of E. coli rRNAs to be substantially flexible. Most of the insertion sites show clear tendency to coincide with the regions of the expansion segments found in eukaryotic rRNAs, implying that eukaryotic rRNAs evolved from prokaryotic rRNAs suffering genetic insertions and selections.