On the evolution of ribosomal RNA

The Peptidyl Transferase Center: a Window to the Past

In his 2001 article, "Translation: in retrospect and prospect," the late Carl Woese made a prescient observation that there was a need for the then-current view of translation to be "reformulated to become an all-embracing perspective about which 21st century Biology can develop" (RNA 7:1055-1067, 2001, https://doi.org/10.1017/s1355838201010615). The quest to decipher the origins of life and the road to the genetic code are both inextricably linked with the history of the ribosome. After over 60 years of research, significant progress in our understanding of how ribosomes work has been made. Particularly attractive is a model in which the ribosome may facilitate an ∼180° rotation of the CCA end of the tRNA from the A-site to the P-site while the acceptor stem of the tRNA would then undergo a translation from the A-site to the P-site. However, the central question of how the ribosome originated remains unresolved. Along the path from a primitive RNA world or an RNA-peptide world to a proto-ribosome world, the advent of the peptidyl transferase activity would have been a seminal event. This functionality is now housed within a local region of the large-subunit (LSU) rRNA, namely, the peptidyl transferase center (PTC). The PTC is responsible for peptide bond formation during protein synthesis and is usually considered to be the oldest part of the modern ribosome. What is frequently overlooked is that by examining the origins of the PTC itself, one is likely going back even further in time. In this regard, it has been proposed that the modern PTC originated from the association of two smaller RNAs that were once independent and now comprise a pseudosymmetric region in the modern PTC. Could such an association have survived? Recent studies have shown that the extant PTC is largely depleted of ribosomal protein interactions. It is other elements like metallic ion coordination and nonstandard base/base interactions that would have had to stabilize the association of RNAs. Here, we present a detailed review of the literature focused on the nature of the extant PTC and its proposed ancestor, the proto-ribosome.

Piece by piece: Building a ribozyme

The ribosome and RNase P are cellular ribonucleoprotein complexes that perform peptide bond synthesis and phosphodiester bond cleavage, respectively. Both are ancient biological assemblies that were already present in the last universal common ancestor of all life. The large subunit rRNA in the ribosome and the RNA subunit of RNase P are the ribozyme components required for catalysis. Here, we explore the idea that these two large ribozymes may have begun their evolutionary odyssey as an assemblage of RNA "fragments" smaller than the contemporary full-length versions and that they transitioned through distinct stages along a pathway that may also be relevant for the evolution of other non-coding RNAs.

Origin and evolution of the ribosome

The modern ribosome was largely formed at the time of the last common ancestor, LUCA. Hence its earliest origins likely lie in the RNA world. Central to its development were RNAs that spawned the modern tRNAs and a symmetrical region deep within the large ribosomal RNA, (rRNA), where the peptidyl transferase reaction occurs. To understand pre-LUCA developments, it is argued that events that are coupled in time are especially useful if one can infer a likely order in which they occurred. Using such timing events, the relative age of various proteins and individual regions within the large rRNA are inferred. An examination of the properties of modern ribosomes strongly suggests that the initial peptides made by the primitive ribosomes were likely enriched for l-amino acids, but did not completely exclude d-amino acids. This has implications for the nature of peptides made by the first ribosomes. From the perspective of ribosome origins, the immediate question regarding coding is when did it arise rather than how did the assignments evolve. The modern ribosome is very dynamic with tRNAs moving in and out and the mRNA moving relative to the ribosome. These movements may have become possible as a result of the addition of a template to hold the tRNAs. That template would subsequently become the mRNA, thereby allowing the evolution of the code and making an RNA genome useful. Finally, a highly speculative timeline of major events in ribosome history is presented and possible future directions discussed.

Evolution of rDNA ITS1 and ITS2 sequences and RNA secondary structures within members of the fungal genera Grosmannia and Leptographium

The two internal transcribed spacers (ITS) of the nuclear ribosomal (r) DNA tandem repeat were examined in ophiostomatoid fungi belonging to the genera Grosmannia and Leptographium and closely-related taxa. Although the DNA sequence of the ITS region evolves rapidly, core features of the RNA secondary structure of the ITS1 and ITS2 segments are conserved. The results demonstrate that structural conservation of GC-rich helical regions is facilitated primarily through compensatory base changes (CBCs), hemi-CBCs, and compensating insertions/deletions (indels), although slippage of the RNA strand is potentially an additional mechanism for maintaining basepairing interactions. The major conclusion of the structural analysis of both ITS segments is that two factors appear to be involved in limiting the type of changes observed: a high GC bias for both ITS1 and ITS2 and structural constraints at the RNA level.

Ribosome origins: the relative age of 23S rRNA Domains

The modern ribosome and its component RNAs are quite large and it is likely that at an earlier time they were much smaller. Hence, not all regions of the modern ribosomal RNAs (rRNA) are likely to be equally old. In the work described here, it is hypothesized that the oldest regions of the RNAs will usually be highly integrated into the machinery. When this is the case, an examination of the interconnectivity between local RNA regions can provide insight to the relative age of the various regions. Herein, we describe an analysis of all known long-range RNA/RNA interactions within the 23S rRNA and between the 23S rRNA and the 16S rRNA in order to assess the interconnectivity between the usual Domains as defined by secondary structure. Domain V, which contains the peptidyl transferase center is centrally located, extensively connected, and therefore likely to be the oldest region. Domain IV and Domain II are extensively interconnected with both themselves and Domain V. A portion of Domain IV is also extensively connected with the 30S subunit and hence Domain IV may be older than Domain II. These results are consistent with other evidence relating to the relative age of RNA regions. Although the relative time of addition of the GTPase center can not be reliably deduced it is pointed out that the development of this may have dramatically affected the progenotes that preceded the last common ancestor.

Comparative analysis of the intergenic spacer regions and population structure of the species complex of the pathogenic yeast

Cryptococcus neoformans is an opportunistic basidiomycete responsible for the high incidence of cryptococcosis in patients with AIDS and in other immune-compromised individuals. This study, which focused on the molecular structure and genetic variability of the two varieties in the C. neoformans and Cryptococcus gattii species complex, employed sequence analysis of the intergenic spacer regions, IGSI and IGSII. The IGS region is the most rapidly evolving region of the rDNA families. The IGSI displayed the most genetic variability represented by nucleotide base substitutions and the presence of long insertions/deletions (indels). In contrast, the IGSII region exhibited less heterogeneity and the indels were not as extensive as those displayed in the IGSI region. Both intergenic spacers contained short, interspersed repeat motifs, which can be related to length polymorphisms observed between sequences. Phylogenetic analysis undertaken in the IGSI, IGSII and IGSI +5S rRNA + IGSII regions revealed the presence of six major phylogenetic lineages, some of which segregated into subgroups. The major lineages are represented by genotypes 1 (C. neoformans var. grubii), genotype 2 (C. neoformans var. neoformans), and genotypes 3, 4, 5 and 6 represented by C. gattii. Genotype 6 is a newly described IGS genotypic group within the C. neoformans species complex. With the inclusion of IGS subgenotypic groups, our sequence analysis distinguished 12 different lineages. Sequencing of clones, which was performed to determine the presence of multiple alleles at the IGS locus in several hybrid strains, yielded a single IGS sequence type per isolate, thus suggesting that the selected group of cloned strains was mono-allelic at this locus. IGS sequence analyses proved to be a powerful technique for the delineation of the varieties of C. neoformans and C. gattii at genotypic and subgenotypic levels.

A secondary structural model of the 28S rRNA expansion segments D2 and D3 from rootworms and related leaf beetles (Coleoptera: Chrysomelidae; Galerucinae)

We analysed the secondary structure of two expansion segments (D2, D3) of the 28S rRNA gene from 229 leaf beetles (Coleoptera: Chrysomelidae), the majority of which are in the subfamily Galerucinae. The sequences were compared in a multiple sequence alignment, with secondary structure inferred primarily from the compensatory base changes in the conserved helices of the rRNA molecules. This comparative approach yielded thirty helices comprised of base pairs with positional covariation. Based on these leaf beetle sequences, we report an annotated secondary structural model for the D2 and D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related beetle taxa. This predicted structure, consisting of seven major compound helices, is mostly consistent with previously proposed models for the D2 and D3 expansion segments in insects. Despite a lack of conservation in the primary structure of these regions of insect 28S rRNA, the evolution of the secondary structure of these seven major motifs may be informative above the nucleotide level for higher-order phylogeny reconstruction of major insect lineages.

A possible tertiary rRNA interaction between expansion segments ES3 and ES6 in eukaryotic 40S ribosomal subunits

Eukaryotic 16S-like ribosomal RNAs contain 12 so-called expansion segments, i.e., sequences not included in the RNA secondary structure core common to eubacteria, archaea, and eukarya. Two of these expansion segments, ES3 and ES6, are juxtaposed in the recent three-dimensional model of the eukaryotic 40S ribosomal subunit. We have analyzed ES3 and ES6 sequences from more than 2900 discrete eukaryotic species, for possible sequence complementarity between the two expansion segments. The data show that ES3 and ES6 could interact by forming a helix consisting of seven to nine contiguous base pairs in almost all analyzed species. We, therefore, suggest that ES3 and ES6 form a direct RNA-RNA contact in the ribosome.

Completion of the sequence of the nuclear ribosomal DNA subunit of Simulium sanctipauli, with descriptions of the 18S, 28S genes and the IGS

We describe the IGS-ETS, 18S and 28S ribosomal gene sequences of Simulium sanctipauli Vajime & Dunbar, a member of the S. damnosum Theobald (Diptera: Simuliidae) complex of blackflies (Diptera: Simuliidae). These regions, together with the ITS-1, ITS-2 and 5.8S rDNA presented elsewhere (accession number U36206), constitute the composite sequence of the entire rDNA unit, making S. sanctipauli the second dipteran species of medical importance for which the entire rDNA has been sequenced. Despite the lack of sequence identity, the IGS of S. sanctipauli showed some structural similarities to other Diptera, i.e. the mosquito Aedes albopictus Skuse (Culicidae), the fruitfly Drosophila melanogaster Meigen (Drosophilidae) and the tsetse Glossina (Glossinidae). Two blocks of tandemly repeated subunits were present in the IGS of S. sanctipauli and, unlike other species of Diptera, they contained no duplications of promoter-like sequences. However, two promoter-like sequences were identified in the unique DNA stretches of the IGS by their sequence similarity to the promoter of Aedes aegypti L. (Diptera: Culicidae). The observed sequence variation can be explained, as in the case of Drosophila spp., by the occurrence of slippage-like and point mutation processes, with unequal crossing-over homogenizing (to a certain extent) the region throughout the gene family and blackfly population. The 18S and 28S rDNA genes show more intraspecific variability within the expansion segments than in the core regions. This is also the case in the interspecific comparison of these genes from S. sanctipauli with those of Simulium vittatum, Ae. albopictus and D. melanogaster. This pattern is typical of many eukaryotes and likely to be the result of a more relaxed functional selection in the expansion segments than on the core regions. The A + T content of the S. sanctipauli genes is high and similar to those of other Diptera. This could be the result of a change in the mutation pressure towards AT in the Diptera lineage.

Tracing the evolution of RNA structure in ribosomes

The elucidation of ribosomal structure has shown that the function of ribosomes is fundamentally confined to dynamic interactions established between the RNA components of the ribosomal ensemble. These findings now enable a detailed analysis of the evolution of ribosomal RNA (rRNA) structure. The origin and diversification of rRNA was studied here using phylogenetic tools directly at the structural level. A rooted universal tree was reconstructed from the combined secondary structures of large (LSU) and small (SSU) subunit rRNA using cladistic methods and considerations in statistical mechanics. The evolution of the complete repertoire of structural ribosomal characters was formally traced lineage-by-lineage in the tree, showing a tendency towards molecular simplification and a homogeneous reduction of ribosomal structural change with time. Character tracing revealed patterns of evolution in inter-subunit bridge contacts and tRNA-binding sites that were consistent with the proposed coupling of tRNA translocation and subunit movement. These patterns support the concerted evolution of tRNA-binding sites in the two subunits and the ancestral nature and common origin of certain structural ribosomal features, such as the peptidyl (P) site, the functional relay of the penultimate stem helix of SSU rRNA, and other structures participating in ribosomal dynamics. Overall results provide a rare insight into the evolution of ribosomal structure.

Distribution of substitution rates and location of insertion sites in the tertiary structure of ribosomal RNA

The relative substitution rate of each nucleotide site in bacterial small subunit rRNA, large subunit rRNA and 5S rRNA was calculated from sequence alignments for each molecule. Two-dimensional and three-dimensional variability maps of the rRNAs were obtained by plotting the substitution rates on secondary structure models and on the tertiary structure of the rRNAs available from X-ray diffraction results. This showed that the substitution rates are generally low near the centre of the ribosome, where the nucleotides essential for its function are situated, and that they increase towards the surface. An inventory was made of insertions characteristic of the Archaea, Bacteria and Eucarya domains, and for additional insertions present in specific eukaryotic taxa. All these insertions occur at the ribosome surface. The taxon-specific insertions seem to arise randomly in the eukaryotic evolutionary tree, without any phylogenetic relatedness between the taxa possessing them.

Ribin, a protein encoded by a message complementary to rRNA, modulates ribosomal transcription and cell proliferation

The control of rRNA transcription, tightly coupled to the cell cycle and growth state of the cell, is a key process for understanding the mechanisms that drive cell proliferation. Here we describe a novel protein, ribin, found in rodents, that binds to the rRNA promoter and stimulates its activity. The protein also interacts with the basal rRNA transcription factor UBF. The open reading frame encoding ribin is 96% complementary to a central region of the large rRNA. This demonstrates that ribosomal DNA-related sequences in higher eukaryotes can be expressed as protein-coding messages. Ribin contains two predicted nuclear localization sequence elements, and green fluorescent protein-ribin fusion proteins localize in the nucleus. Cell lines overexpressing ribin exhibit enhanced rRNA transcription and faster growth. Furthermore, these cells significantly overcome the suppression of rRNA synthesis caused by serum deprivation. On the other hand, the endogenous ribin level correlates positively with the amount of serum in the medium. The data show that ribin is a limiting stimulatory factor for rRNA synthesis in vivo and suggest its involvement in the pathway that adapts ribosomal transcription and cell proliferation to physiological changes.

Molecular evolution of the small subunit ribosomal DNA in woodlice (Crustacea, Isopoda, Oniscidea) and implications for Oniscidean phylogeny

The small subunit ribosomal DNA (ssu rDNA) of 13 isopods was sequenced. The entire length of the ribosomal gene is unusually long, resulting from the presence of five expansion elements accounting for more than 40% of the gene. We found that in terrestrial isopods the length of the ssu rDNA ranges from 2414 bp (Ligidium hypnorum) to 3537 bp (Cubaris murina). This is the longest metazoan ssu rDNA reported to date. The conserved regions are highly informative for analysis of the early nodes of the tree, whereas the variable expansion elements are better suited to reconstruction of the branching pattern between closely related taxa. The suggested relationship among Synochaeta, Crinochaeta, and Diplochaeta based on the conserved regions confirms that based on previous morphological analyses. In contrast, the phylogeny within the Crinochaeta based on the entire ssu rDNA including the variable domains is in conflict with that based on most of the morphological analyses. The phylogenetic analyses of the ssu rDNA support a repeated independent evolution of the three different types of pleopodal lungs in the Crinochaeta.

Cirripede phylogeny using a novel approach: molecular morphometrics

We present a new method using nucleic acid secondary structure to assess phylogenetic relationships among species. In this method, which we term "molecular morphometrics," the measurable structural parameters of the molecules (geometrical features, bond energies, base composition, etc.) are used as specific characters to construct a phylogenetic tree. This method relies both on traditional morphological comparison and on molecular sequence comparison. Applied to the phylogenetic analysis of Cirripedia, molecular morphometrics supports the most recent morphological analyses arguing for the monophyly of Cirripedia sensu stricto (Thoracica + Rhizocephala + Acrothoracica). As a proof, a classical multiple alignment was also performed, either using or not using the structural information to realign the sequence segments considered in the molecular morphometrics analysis. These methods yielded the same tree topology as the direct use of structural characters as a phylogenetic signal. By taking into account the secondary structure of nucleic acids, the new method allows investigators to use the regions in which multiple alignments are barely reliable because of a large number of insertions and deletions. It thus appears to be complementary to classical primary sequence analysis in phylogenetic studies.

Fibrillarin-associated box C/D small nucleolar RNAs in Trypanosoma brucei. Sequence conservation and implications for 2'-O-ribose methylation of rRNA

We report the identification of 17 box C/D fibrillarin-associated small nucleolar RNAs (snoRNAs) from the ancient eukaryote, Trypanosoma brucei. To systematically isolate and characterize these snoRNAs, the T. brucei cDNA for the box C/D snoRNA common protein, fibrillarin, was cloned and polyclonal antibodies to the recombinant fibrillarin protein were generated in rabbits. Immunoprecipitations from T. brucei extracts with the anti-fibrillarin antibodies indicated that this trypanosomatid has at least 30 fibrillarin-associated snoRNAs. We have sequenced seventeen of them and designated them TBR for T. brucei RNA 1-17. All of them bear conserved box C, D, C', and D' elements, a hallmark of fibrillarin-associated snoRNAs in eukaryotes. Fourteen of them are novel T. brucei snoRNAs. Fifteen bear potential guide regions to mature rRNAs suggesting that they are involved in 2'-O-ribose methylation. Indeed, eight ribose methylations have been mapped in the rRNA at sites predicted by the snoRNA sequences. Comparative genomics indicates that six of the seventeen are the first trypanosome homologs of known yeast and vertebrate methylation guide snoRNAs. Our results indicate that T. brucei has many fibrillarin-associated box C/D snoRNAs with roles in 2'-O-ribose methylation of rRNA and that the mechanism for targeting the nucleotide to be methylated at the fifth nucleotide upstream of box D or D' originated in early eukaryotes.

The role of recombination and mutation in 16S-23S rDNA spacer rearrangements

The intragenomic heterogeneity of the bacterial intergenic (16S-23S rDNA) spacer region (ISR) was analysed from the following species in which sequences for the complete rRNA operon (rrn) set have been determined (rrn number): Enterococcus faecalis (6) and E. faecium (6), Bacillus subtilis (10), Staphylococcus aureus (9), Vibrio cholerae (4), Haemophilus influenzae (6) and Escherichia coli (7). It was found that some spacer sequence blocks were highly conserved between operons of a genome, whereas the presence of others was variable. When these variations were analysed using the program PLATO and partial likelihood phylogenies determined by DNAml for each operon set, three regions showed significant (Z>3.3) spatial variation [Region I was 78-184 nt long (2.1<Z<49.4), Region II was 10-60 nt long (3.7<Z<23)] and Region III was 6 nt long (3.4<Z>4.4) possibly due to recombination or selection. Within Region I, there was sequence block variation in all operon sets [some operons contained tRNA genes (tRNAala, tRNAile or tRNAglu), whereas others had sequence blocks such as VS2 (S. aureus) or rsl (E. coli)]. Q Analysis of the ISR sequence from E. faecalis and E. faecium showed that there was more interspecies than intraspecies variation (both in DNA sequence and in the presence or absence of blocks). Dot matrix analysis of the sequence blocks in the nine rrn ISRs from S. aureus showed that there was significant homology between VS2 and VS5/VS6. Furthermore, repeat motifs with only A or T were present in higher copy numbers in VS5/VS6 than in VS2. Since these sequence blocks (VS2 and VS5-VS6) are related, intragenic evolution resulting in AT expansion may have occurred between these two regions. A model is proposed that postulates a role for recombination and AT-expansion in intra-genomic ISR variations. This process may represent a general mechanism of concerted evolution for bacterial ISR rearrangements.

Novel processing in a mammalian nuclear 28S pre-rRNA: tissue-specific elimination of an 'intron' bearing a hidden break site

Splitting and apparent splicing of ribosomal RNA, both previously unknown in vertebrates, were found in rodents of the genus Ctenomys. Instead of being formed by a single molecule of 4.4 kb, 28S rRNA is split in two molecules of 2.6 and 1.8 kb. A hidden break, mapping within a 106 bp 'intron' located in the D6 divergent region, is expressed in mature ribosomes of liver, lung, heart and spleen, as well as in primary fibroblast cultures. Testis-specific processing eliminates the intron and concomitantly the break site, producing non-split 28S rRNA molecules exclusively in this organ. The intron is flanked by two 9 bp direct repeats, revealing the acquisition by insertion of a novel rRNA processing strategy in the evolution of higher organisms.

The secondary structure of Nosema apis large subunit ribosomal RNA

The microsporidia are a group of obligate intracellular eukaryotic parasites, that lack mitochondria. Their ribosomes show several prokaryote-like features. This paper presents the secondary structure of the large subunit ribosomal RNA (LSU rRNA) of the microsporidium Nosema apis. With its 2481 bases, it is the shortest known non-mitochondrial LSU rRNA. The seemingly prokaryote-like features of the molecule cannot be used as evidence for the ancient origin of the microsporidia. The reduction in size can be attributed to changes in the regions of the LSU rRNA that are known to show great variability in length and sequence within the eukaryotes. The lack of fragmentation commonly seen in other eukaryotes may also be a derived feature.

Divergence towards a dead end? Cleavage of the divergent domains of ribosomal RNA in apoptosis

In several cases of apoptotic death the large ribosomal subunit 28S rRNA is specifically cleaved. The cleavages appear at specific sites within those domains of the rRNA molecule that have shown exceptional high divergence in evolution (D domains). The cleavages accompany rather than precede apoptosis, and there is a positive, but not complete, correlation between rRNA cleavage and internucleosomal DNA fragmentation. Most cell types studied so far show two alternative cleavage pathways that are mutually exclusive. Cleavage can either start in the D8 domain with secondary cuts within a subdomain of D2 (D2c), or in the D2 domain with subsequent excision of the D2c subdomain. The latter pathway is of particular interest since D2 (unlike D8) is normally inaccessible for RNase attack. That apoptosis specifically affects the ribosomal divergent domains suggests that these domains, which make up roughly 25% of total cellular RNA, might have evolved to serve functions related to apoptosis. Future studies will be directed to test the hypothesis that rRNA fragmentation may be part of an apoptotic program directed against the elimination of illegitimate (viral?) polynucleotides.

Secondary structures and features of the 18S, 5.8S and 26S ribosomal RNAs from the Apicomplexan parasite Toxoplasma gondii

The two major subunits of the ribosomal RNA (rRNA) of Toxoplasma gondii, 18S and 26S, as well as 5.8S, have been sequenced and folded according to known consensus and established secondary structures. Conserved and variable nucleotide (nt) regions were identified using multiple alignments with rRNA sequences of selected organisms. The 18S rRNA showed a well conserved core structure of 48 stems and a hypervariable V4 region identified four additional stems including a pseudoknot. The 18S rRNA contained an additional helix in the V2 region located between nt 204 to 258. We noted that T. gondii 18S does not have a true V6 region, but was organized as a motif of a simple stem. T. gondii 26S had a conserved core structure of 83 stems and its expansion segments, so-called divergent domains, demonstrated a high degree of similarity with secondary structures from rRNA of dinoflagellates and ciliates. For the T. gondii 26S sequence, we found two additional stems, D3d and D3e, composed of 140 nt having a higher deltaG value. These segments are absent from the prokaryotic rRNA structures, whereas the hypervariable V4 region of the small subunit is not as variable. The well preserved structures could indicate an additional function for the eukaryotic ribosome.

A quantitative map of nucleotide substitution rates in bacterial rRNA

A recently developed method for estimating the variability of nucleotide sites in a sequence alignment [Van de Peer, Y., Van der Auwera, G. and De Wachter, R. (1996) J. Mol. Evol. 42, 201-210] was applied to bacterial 16S, 5S and 23S rRNAs. In this method, the variability of each nucleotide site is defined as its evolutionary rate relative to the average evolutionary rate of all the nucleotide sites of the molecule. Spectra of evolutionary rates were calculated for each rRNA and show the fastest evolving sites substituting at rates more than 1000 times that of the slowest ones. Variability maps are presented for each rRNA, consisting of secondary structure models where the variability of each nucleotide site is indicated by means of a colored dot. The maps can be interpreted in terms of higher order structure, function and evolution of the molecules and facilitate the selection of areas suitable for the design of PCR primers and hybridization probes. Variability measurement is also important for the precise estimation of evolutionary distances and the inference of phylogenetic trees.

Sequence variations in the large-subunit ribosomal RNA gene of Ammonia (Foraminifera, Protozoa) and their evolutionary implications

An unusually high divergence was observed in the ribosomal RNA genes of a free-living population of foraminifera belonging to the genus Ammonia. The sequences of a large-subunit (LSU) rDNA expansion segment D1 and flanking regions were obtained from 20 specimens named Ammonia sp. 1 and Ammonia sp. 2. The sequence divergence between the two species averages 14%. Within each species it ranges from 0.2% to 7.1% in Ammonia sp. 1 and from 0.7% to 2.3% in Ammonia sp. 2. We did not find two specimens having identical sequences. Moreover, in opposition to the generally accepted view, rDNA sequence variations were also found within a single individual. The variations among several rDNA copies in a single specimen of Ammonia may reach up to 4.9%. Most of the observed variations result from multiplication of CA or TA serial repeats occurring in two particularly variable regions. For single base changes, C-T transitions are most frequently observed. We discuss the evolution of expansion segments and their use for phylogenetic studies.

Discrimination of S. aureus strains by PCR for r-RNA gene spacer size polymorphism and comparison to SmaI macrorestriction patterns

The size polymorphism of the internal spacer between the 16 S and the 23 S r-RNA genes was studied in S. aureus with the aid of PCR. The patterns of corresponding PCR products were compared with SmaI-generated macrorestriction patterns for definite propagating strains of S. aureus typing phages, for strains with phage pattern 29, phage-group II patterns, phage pattern 94, 96, phage pattern 95 and epidemic methicillin-resistant strains (MRSA). The spacer length polymorphism did not prove to be as discriminative as genomic DNA fragment patterns. However, as shown for S. aureus with phage patterns 29; group II; 94, 96; 95 and also for 4 out of 6 epidemic MRSA, unique patterns of r-RNA gene spacers probably indicate a relatedness among strains which is also suggested by SmaI macrorestriction patterns.

Simplicity-correlated size growth of the nuclear 28S ribosomal RNA D3 expansion segment in the crustacean order Isopoda

The expansion segments within the eukaryote nuclear 23S-like ribosomal RNA molecule are now well characterized in many diverse organisms. A different base compositional bias, a higher propensity for size variability, and an increased evolutionary rate distinguish these regions from the universally conserved "core" regions of the molecule. In addition, some expansion segments of higher eukaryotes exhibit significant sequence simplicity which is hypothesized to occur by slippage-mediated mutational processes. We describe the discovery of extreme size variation of the D3 expansion segment in the crustacean order Isopoda. Among 11 species D3 varies in size from 180 to 518 nucleotides but maintains a homologous secondary structure. The D3 size is significantly positively correlated to relative simplicity factor (RSF), indicating that growth is most likely by insertion of simple sequences. D3 size and RSF correlate approximately with a morphology-based phylogeny, and within oniscideans RSF increases as more recent divergences occur. The D3 of Armadillidium vulgare, with an RSF of 1.87, is the highest value recorded for any known expansion segment. Regions of high sequence simplicity in nuclear ribosomal RNA were previously only known from the higher vertebrate lineage. Here we demonstrate that this phenomenon occurs in a more extreme condition within a monophyletic invertebrate lineage. The extreme size changes identified could indicate that expansion segments are an extraneous element in the functioning ribosome.

The 5S rRNA and the rRNA intergenic spacer of the two varieties of Cryptococcus neoformans

The intergenic spacers (IGS) separating the 23S-like and 16S-like rDNAs of the two varieties of the human pathogenic fungus Cryptococcus neoformans were amplified, cloned and sequenced. The C. neoformans var. neoformans IGS was 2421 nt with 5S rRNA at positions 1228-1345 3' of the 23S-like rRNA. The C. neoformans var. gattii IGS was 2480 nt with 5S rRNA at positions 1268-1385 3' of the 23S-like rRNA. For both varieties the 5S rDNA genes were in the same orientation as the 16S-5.8-23S genes and encode a 118 nt molecule of identical sequence. Phylogenetic comparison of C. neoformans 5S rDNA with that of other fungi placed this fungus in close relationship with other basidiomycetes including Tremella mesenterica, Bullera alba, and Cryptococcus laurentii. A secondary structure model for the deduced 5S rRNA was constructed by comparative sequence analysis. Polymerase chain reaction-amplified IGS of 12 C. neoformans var. neoformans strains revealed extensive size variation ranging from 100 to 300 nt. Size variation between strains in the length of the IGS may be useful for distinguishing strains. Structurally, the IGS were characterized by the presence of occasional short direct GC-rich 19-nt repeats. Overall IGS sequence identity between the C. neoformans varieties was only 78.5%, in sharp contrast to the identical or nearly identical sequences for the rDNA genes, and suggests rapid evolution for IGS sequences.

Partial nucleotide sequence of a single ribosomal RNA coding region and secondary structure of the large subunit 25 s rRNA of Candida albicans

A rDNA cistron of Candida albicans strain WO-1 was cloned and the ITS1, ITS2, 5.8 s rDNA and 25 s rDNA coding regions sequenced in their entirety. These sequences were compared to those of three related yeast species (Saccharomyces cerevisiae, Saccharomyces carlsbergensis, and Thermomyces lanuginosus), and the 5.8 s rDNA was compared to seven additional 5.8 s rDNAs from organisms ranging in complexity from D. discoideum to H. sapiens. The C. albicans ITS regions are shorter than those of most other eukaryotes. The 25 s and 5.8 s rDNA sequences were folded into a secondary structure model based on comparative methods. In a comparison of regional similarities between the large subunit rDNAs of C. albicans, the three related yeasts and other eukaryotes, it is demonstrated that the additional sequences not present in the E. coli 23 s rDNA are more variable than the regions present in both prokaryotes and eukaryotes.

Processing of eukaryotic ribosomal RNA

In summary, it can be argued that the understanding of eukaryotic rRNA processing is no less important than the understanding of mRNA maturation, since the capacity of a cell to carry out protein synthesis is controlled, in part, by the abundance of ribosomes. Processing of pre-rRNA is highly regulated, involving many cellular components acting either alone or as part of a complex. Some of these components are directly involved in the modification and cleavage of the precursor rRNA, while others direct the packaging of the rRNA into ribosome subunits. As is the case for pre-mRNA processing, snoRNPs are clearly involved in eukaryotic rRNA processing, and have been proposed to assemble with other proteins into at least one complex called a "processosome" (17), which carries out the ordered processing of the pre-rRNA and its assembly into ribosomes. The formation of a processing complex clearly makes possible the regulation required to coordinate the abundance of ribosomes with the physiological and developmental changes of a cell. It may be that eukaryotic rRNA processing is even more complex than pre-mRNA maturation, since pre-rRNA undergoes extensive nucleotide modification and is assembled into a complex structure called the ribosome. Undoubtedly, features of the eukaryotic rRNA-processing pathway have been conserved evolutionarily, and the genetic approach available in yeast research (6) should provide considerable knowledge that will be useful for other investigators working with higher eukaryotic systems. Interestingly, it was originally hoped that the extensive work and understanding of bacterial ribosome formation would provide a useful paradigm for the process in eukaryotes. However, although general features of ribosome structure and function are highly conserved between bacterial and eukaryotic systems, the basic strategy in ribosome biogenesis seems to be, for the most part, distinctly different. Thus, the detailed molecular mechanisms for rRNA processing in each kingdom will have to be independently deciphered in order to elucidate the features and regulation of this important process for cell survival.

rRNA gene restriction patterns as an epidemiological marker in nosocomial outbreaks of Staphylococcus aureus infections

rRNA gene restriction patterns (ribotyping) were compared with phage typing, serotyping, enterotoxins and exfoliatin production in the analysis of 26 Staphylococcus aureus strains isolated from two different nosocomial outbreaks. Total DNA was cleaved by EcoRI restriction endonuclease. After agarose gel electrophoresis and Southern transfer, the hybridization of the membranes was done with radiolabelled 16S rRNA gene from Bacillus subtilis inserted into a plasmid vector. Six to 13 fragments were visualized. A core of common fragments was discerned for all strains tested. A full correlation between ribotyping and conventional markers was observed in only one of the outbreaks studied. In both outbreaks, ribotyping proved helpful in characterizing otherwise untypable strains.

Is higher-order structure conserved in eukaryotic ribosomal DNA intergenic spacers?

Computer-based structural analysis of the ribosomal DNA intergenic spacer (IGS) from the mosquito Aedes albopictus revealed a potential to form strong and extensive secondary structures throughout a 4.7-kilobase (kb) region. The predicted stability of secondary structures was particularly high within a 3.15-kb region containing 17 tandem 201 base-pair subrepeats. Similarly strong secondary structure potential was also found when IGS subrepeats were analyzed from 17 phylogenetically diverse eukaryotes, including vertebrates, invertebrates, and plants. Conservation of higher-order structure potential in the IGS region of ribosomal DNA may reflect evolutionary and functional constraints on chromatin organization, transcriptional regulation of the ribosomal RNA genes, and/or transcript processing and stability.

Extrachromosomal ribosomal DNA of Didymium iridis: sequence analysis of the large subunit ribosomal RNA gene and sub-telomeric region

The ribosomal DNA of the myxomycete Didymium iridis is organized as extrachromosomal linear molecules of about 20 kb, containing only one transcription unit of the ribosomal RNA genes. We have determined the sequence of the large subunit ribosomal RNA (LSU rRNA) gene as well as the sub-telomeric and telomeric regions. The LSU rRNA gene was found to encode a 3857 nucleotide-long LSU rRNA, interrupted by a transcribed spacer and two group I introns. A complete secondary structure model of D. iridis LSU rRNA has been constructed. The compact sub-telomeric region of D. iridis rDNA was found to contain several directly repeated sequence elements that include the simple telomere motif TTAGGG. Based on pairwise comparisons of LSU rRNA sequences, the time of divergence between the two myxomycete genera Didymium and Physarum was estimated.

A hidden break in the 28.0S rRNA from Diphyllobothrium dendriticum

Nondenatured and denatured total RNA from the tapeworm Diphyllobothrium dendriticum (Cestoda) was analysed by agarose gel electrophoresis. It was found that the large subunit ribosomal RNA (lrRNA) is 28.0S and the small subunit ribosomal RNA (srRNA) is 19.5S. Following denaturation the 28.0S rRNA was disrupted into a 19.5S subfragment and a 20.7S subfragment due to the presence of a centrally located hidden break. By hybridization of Northern blot membranes with oligonucleotide probes specific for the 5'- and 3'-ends of the lrRNA respectively, we have shown that the 19.5S subfragment is from the 5'-end (the alpha-subfragment) and the 20.7S subfragment from the 3'-end (the beta-subfragment) of the 28.0S rRNA of D. dendriticum.

Secondary structure as a constraint on the evolution of a plant viral satellite RNA

The genetic variability and evolution of the satellite RNA (satRNA) of cucumber mosaic virus (CMV) was analyzed. Twenty-five CMV-satRNAs compared clustered into three main groups, and no correlation was found between genetic proximity and other characteristics (pathogenicity, geographical origin) of the satRNAs. Values for the number of nucleotide substitutions per site between any two satRNAs suggest that divergence is checked by functional constraints. The analysis of mutations relative to an ancestral sequence, and the number of substitutions per site at first, second and third positions of codons in putative open reading frames, show that the variation of CMV-satRNAs does not follow a pattern typical of coding sequences, and indicates that preservation of the sequence of encoded products is not a constraint to evolution. On the other hand, when the observed variation was analyzed relative to a secondary structure model proposed for CMV-satRNAs, several lines of evidence indicated that the maintenance of the secondary structure is a constraint to evolution: the number of substitutions per site, the number of point insertions and deletions and the number of base substitutions that would disrupt base-pairing were significantly higher for unpaired than for base-paired positions. Also, compensatory mutations at base-paired positions occurred more frequently than expected from random. The results suggest that CMV-satRNAs are non-coding, functional RNAs whose biology would be determined by their direct interaction with components of the host and/or the helper virus.

Nucleotide sequence of a mosquito 18S ribosomal RNA gene

We have sequenced an 18S ribosomal RNA gene from the mosquito, Aedes albopictus. Computer alignment of the 1950 nucleotide coding region (56% A + T) with 18S rRNA sequences from two insect and three vertebrate species revealed greater sequence divergence among the insects than among the vertebrates. Sequence alignments showed that variable region V4, which has been considered to be the most poorly conserved domain in the 18S rRNA gene, was better conserved among insects and vertebrates than was the V6 domain.

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.

Six group I introns and three internal transcribed spacers in the chloroplast large subunit ribosomal RNA gene of the green alga Chlamydomonas eugametos

The chloroplast large subunit rRNA gene of Chlamydomonas eugametos and its 5' flanking region encoding tRNA(Ile) (GAU) and tRNA(Ala) (UGC) have been sequenced. The DNA sequence data along with the results of a detailed RNA analysis disclosed two unusual features of this green algal large subunit rRNA gene: (1) the presence of six group I introns (CeLSU.1-CeLSU.6) whose insertion positions have not been described previously, and (2) the presence of three short internal transcribed spacers that are post-transcriptionally excised to yield four rRNA species of 280, 52, 810 and 1720 nucleotides, positioned in this order (5' to 3') in the primary transcript. Together, these RNA species can assume a secondary structure that is almost identical to that proposed for the 23 S rRNA of Escherichia coli. All three internal transcribed spacers map to variable regions of primary sequence and/or potential secondary structure, whereas all six introns lie within highly conserved regions. The first three introns are inserted within the sequence encoding the 810 nucleotide rRNA species and map within domain II of the large subunit rRNA structure; the remaining introns, found in the sequence encoding the 1720 nucleotide rRNA species, lie within either domain IV or V, as is the case for all other large subunit rDNA introns that have been documented to date. CeLSU.5 and CeLSU.6 each contain a long open reading frame (ORF) of more than 200 codons. While the CeLSU.6 ORF is not related to any known ORFs, the CeLSU.5 ORF belongs to a family of ORFs that have been identified in Podospora and Neurospora mitochondrial group I introns. The finding that a polymorphic marker showing unidirectional gene conversion during crosses between C. eugametos and Chlamydomonas moewusii is located within the CeLSU.5 ORF makes it likely that this intron is a mobile element and that its ORF encodes a site-specific endonuclease promoting the transfer of the intron DNA sequence.

'Compensatory slippage' in the evolution of ribosomal RNA genes

The distribution patterns of shared short repetitive motifs in the expansion segments of the large subunit rRNA genes of different species show that these segments are coevolving as a set and that in two examined vertebrate species the RNA secondary structures are conserved as a consequence of runs of motifs in one region being compensated by complementary motifs in another. These unusual processes, involving replication-slippage, have implications for the evolution of ribosomal RNA and for the use of the rDNA multigene family as a 'molecular clock' for assessing relationships between species.

Evolutionary origin of the U6 small nuclear RNA intron

U6 is the most conserved of the five small nuclear RNAs known to participate in pre-mRNA splicing. In the fission yeast Schizosaccharomyces pombe, the single-copy gene encoding this RNA is itself interrupted by an intron (T. Tani and Y. Ohshima, Nature (London) 337:87-90, 1989). Here we report analysis of the U6 genes from all four Schizosaccharomyces species, revealing that each is interrupted at an identical position by a homologous intron; in other groups, including ascomycete and basidiomycete fungi, as well as more distantly related organisms, the U6 gene is colinear with the RNA. The most parsimonious interpretation of our data is that the ancestral U6 gene did not contain an intron, but rather, it was acquired via a single relatively recent insertional event.

Evolutionary origin of the U6 small nuclear RNA intron

U6 is the most conserved of the five small nuclear RNAs known to participate in pre-mRNA splicing. In the fission yeast Schizosaccharomyces pombe, the single-copy gene encoding this RNA is itself interrupted by an intron (T. Tani and Y. Ohshima, Nature (London) 337:87-90, 1989). Here we report analysis of the U6 genes from all four Schizosaccharomyces species, revealing that each is interrupted at an identical position by a homologous intron; in other groups, including ascomycete and basidiomycete fungi, as well as more distantly related organisms, the U6 gene is colinear with the RNA. The most parsimonious interpretation of our data is that the ancestral U6 gene did not contain an intron, but rather, it was acquired via a single relatively recent insertional event.

Fourteen internal transcribed spacers in the circular ribosomal DNA of Euglena gracilis

Cytoplasmic ribosomes from Euglena gracilis contain 16 rRNA components. These include the typical 5 S, 5.8 S and 19 S rRNAs that are found in other eukaryotes as well as 13 discrete small RNAs that interact to form the equivalent of eukaryotic 25-28 S rRNA (accompanying paper). We have utilized DNA sequencing techniques to establish that genes for all of these RNAs, with the exception of 5 S rRNA, are encoded by the 11,500 base-pair circular rDNA of E. gracilis. We have determined the relative positions of the coding regions for the 19 S rRNA and the 14 components (including 5.8 S rRNA) of the large subunit rRNA, thereby establishing that the genes for each of these rRNAs are separated by internal transcribed spacers. We conclude that sequences corresponding to these spacers are removed post-transcriptionally from a high molecular weight pre-rRNA, resulting in a multiply fragmented large subunit rRNA. Internal transcribed spacers, in positions analogous to some of these additional Euglena rDNA spacers, have been found in the rDNA of other organisms and organelles. This finding supports the view that at least some internal transcribed spacers may have been present at an early stage in the evolution of rRNA genes.

Sixteen discrete RNA components in the cytoplasmic ribosome of Euglena gracilis

We have isolated cytoplasmic ribosomes from Euglena gracilis and characterized the RNA components of these particles. We show here that instead of the four rRNAs (17-19 S, 25-28 S, 5.8 S and 5 S) found in typical eukaryotic ribosomes, Euglena cytoplasmic ribosomes contain 16 RNA components. Three of these Euglena rRNAs are the structural equivalents of the 17-19 S, 5.8 S and 5 S rRNAs of other eukaryotes. However, the equivalent of 25-28 S rRNA is found in Euglena as 13 separate RNA species. We demonstrate that together with 5 S and 5.8 S rRNA, these 13 RNAs are all components of the large ribosomal subunit, while a 19 S RNA is the sole RNA component of the small ribosomal subunit. Two of the 13 pieces of 25-28 S rRNA are not tightly bound to the large ribosomal subunit and are released at low (0 to 0.1 mM) magnesium ion concentrations. We present here the complete primary sequences of each of the 14 RNA components (including 5.8 S rRNA) of Euglena large subunit rRNA. Sequence comparisons and secondary structure modeling indicate that these 14 RNAs exist as a non-covalent network that together must perform the functions attributed to the covalently continuous, high molecular weight, large subunit rRNA from other systems.

Genome structure and evolution of Naegleria and its relatives

In the 4 yr since the molecular biology of DNA in Naegleria was last reviewed several major advances have been made, and these are reviewed here: isolation and characterization of mitochondrial and ribosomal DNAs; enumeration of chromosomal DNAs by pulsed field gel electrophoresis; sequence analysis of differentially expressed genes; phylogenetic placement of the genus Naegleria among the eukaryotes and Naegleria species within the genus.