Functional equivalence of structurally distinct ribosomes in the malaria parasite, Plasmodium berghei

Unlike most eukaryotes, many apicomplexan parasites contain only a few unlinked copies of ribosomal RNA (rRNA) genes. Based on stage-specific expression of these genes and structural differences among the rRNA molecules it has been suggested that Plasmodium spp. produce functionally different ribosomes in different developmental stages. This hypothesis was investigated through comparison of the structure of the large subunit rRNA molecules of the rodent malaria parasite, Plasmodium berghei, and by disruption of both of the rRNA gene units that are transcribed exclusively during development of this parasite in the mosquito (S-type rRNA gene units). In contrast to the human parasite, Plasmodium falciparum, we did not find evidence of structural differences in core regions of the distinct large subunit rRNAs which are known to be associated with catalytic activity including the GTPase site that varies in P. falciparum. Knockout P. berghei parasites lacking either of the S-type gene units were able to complete development in both the vertebrate and mosquito hosts. These results formally exclude the hypothesis that two functionally different ribosome types distinct from the predominantly blood stage-expressed A-type ribosomes, are required for development of all Plasmodium species in the mosquito. The maintenance of two functionally equivalent rRNA genes might now be explained as a gene dosage phenomenon.     

Evolution and phylogeny of the heterogeneous cytosolic SSU rRNA genes in the genus Plasmodium

Unlike other eukaryotes, malaria parasites in the genus Plasmodium have structurally and functionally different paralogous copies of the cytosolic (cyto-) SSU rRNA (18S rRNA) gene that are expressed at different developmental stages. In P. falciparum, P. vivax, and P. berghei, A-type cyto-SSU rRNA is expressed in asexual stage, while S-type in sporozoite stage. A third type (O-type) has been described in P. vivax. It is expressed only in oocyst stage in the mosquito. Recently, it has been shown that the maintenance of heterogeneous cyto-SSU rRNAs in Plasmodium can be modeled as a birth-and-death process under strong purifying selection [Rooney, A.P., 2004. Mechanisms underlying the evolution and maintenance of functionally heterogeneous 18S rRNA genes in Apicomplexans. Mol. Biol. Evol. 21, 1704-1711]. In this study, we performed detailed phylogenetic analyses of Plasmodium cyto-SSU rRNAs with special emphasis on the evolution of multi-copy genes in simian Plasmodium species. We sequenced paralogous copies of the cyto-SSU rRNA genes from an African simian Plasmodium species, P. gonderi, and Asian simian Plasmodium species, P. fragile, P. coatneyi, P. inui, P. hylobati, P. fieldi, P. simiovale, and P. cynomolgi. Interestingly, all Asian simian Plasmodium species have a single S-type-like gene and several A-type-like genes. Alignment analysis demonstrated for the first time that an approximately 50-residue insertion in the V7 variable region near the stem 43 is shared exclusively by the S-type-like sequences of the Asian simian Plasmodium species and the S- and O-type sequences of P. vivax. We comprehensively analyzed all cyto-SSU rRNA sequences of the genus Plasmodium currently available in the database. Phylogenetic analyses of all publicly available cyto-SSU rRNA sequences for the genus Plasmodium clearly demonstrated that gene duplication events giving rise to A- and S-type-like sequences took place independently at least three times in the Plasmodium evolution, supporting the hypothesis that these genes evolve according to a birth-and-death model.     

Replacement of the L11 binding region within E.coli 23S ribosomal RNA with its homologue from yeast: in vivo and in vitro analysis of hybrid ribosomes altered in the GTPase centre.

Replacement of the protein L11 binding domain within Escherichia coli 23S ribosomal RNA (rRNA) by the equivalent region from yeast 26S rRNA appeared to have no effect on the growth rate of E.coli cells harbouring a plasmid carrying the mutated rrnB operon. The hybrid rRNA was correctly processed and assembled into ribosomes, which accumulated normally in polyribosomes. Of the total ribosomal population, < 25% contained wild-type, chromosomally encoded rRNA; the remainder were mutant. The hybrid ribosomes supported GTP hydrolysis dependent upon E.coli elongation factor G, although at a somewhat reduced rate compared with wild-type particles, and were sensitive to the antibiotic, thiostrepton, a potent inhibitor of ribosomal GTPase activity that binds to 23S rRNA within the L11 binding domain. That thiostrepton could indeed bind to the mutant ribosomes, although at a reduced level relative to that seen with wild-type ribosomes, was confirmed in a non-equilibrium assay. The rationale for the ability of the hybrid ribosomes to bind the antibiotic, given that yeast ribosomes do not, was provided when yeast rRNA was shown by equilibrium dialysis to bind thiostrepton only 10-fold less tightly than did E.coli rRNA. The extreme conservation of secondary, but not primary, structure in this region between E.coli and yeast rRNAs allows the hybrid ribosomes to function competently in protein synthesis and also preserves the interaction with thiostrepton.     

Poliovirus genome RNA hybridizes specifically to higher eukaryotic rRNAs.

The RNA genome of poliovirus hybridizes to 28S and 18S rRNAs of higher eukaryotes under stringent conditions. The hybridization detected by Northern blot analyses is specific since little or no signal was detected for yeast or prokaryotic rRNAs or other major cellular RNAs. Southern blot analysis of DNA clones of mouse rRNA genes leads us to conclude that several regions of 28S rRNA, and at least one region in 18S rRNA, are involved in the hybridization to polio RNA, and that G/C regions are not responsible for this phenomenon. We have precisely mapped one of these hybridizing regions in both molecules. Computer analysis confirms that extensive intermolecular base-pairing (81 out of 104 contiguous bases in the rRNA strand) could be responsible for this one particular site of interaction (polio genome, bases 5075-5250; 28S rRNA, bases 1097-1200). We discuss the possible functional and/or evolutionary significance of this novel type of interaction.     

Differentiation of oocyte- and somatic-type 5S rRNAs in animals

In some amphibians and bony fishes, oocyte- and somatic-type 5S rRNA genes are expressed differently in oocytes and somatic cells. In order to determine at what stage of animal evolution this differential expression system appeared and how it is regulated, the sequences of oocyte and somatic 5S rRNAs from three invertebrates (sea urchin, sea hare, and silkworm) and two vertebrates (lamprey and chick) were analyzed. It was found that the oocyte 5S rRNA from lamprey consists of two components, while its somatic 5S rRNA consists of only one. In other animals, such differential expression of 5S rRNA in oocytes and somatic cells was not seen. A phylogenetic tree of 63 animal 5S rRNAs was constructed by means of the parsimony method, and the evolution of oocyte and somatic-type 5S rRNAs was discussed.     

A comparison of the solution structures and conformational properties of the somatic and oocyte 5S rRNAs of Xenopus laevis

The secondary and tertiary structures of Xenopus oocyte and somatic 5S rRNAs were investigated using chemical and enzymatic probes. The accessibility of both RNAs towards single-strand specific nucleases (T1, T2, A and S1) and a helix-specific ribonuclease from cobra venom (RNase V1) was determined. The reactivity of nucleobase N7, N3 and N1 positions towards chemical probes was investigated under native (5 mM MgCl2, 100 mM KCl, 20 degrees C) and semi-denaturing (1 mM EDTA, 20 degrees C) conditions. Ethylnitrosourea was used to identify phosphates not reactive towards alkylation under native conditions. The results obtained confirm the presence of the five helical stems predicted by the consensus secondary structure model of 5S rRNA. The chemical reactivity data indicate that loops C and D are involved in a number of tertiary interactions, and loop E folds into an unusual secondary structure. A comparison of the data obtained for the two types of Xenopus 5S rRNA indicates that the conformations of the oocyte and somatic 5S rRNAs are very similar. However, the data obtained with nucleases under native conditions, and chemical probes under semi-denaturing conditions, reveal that helices III and IV in the somatic 5S rRNA are less stable than the same structures in oocyte 5S rRNA. Using chimeric 5S rRNAs, it was possible to demonstrate that the relative resistance of oocyte 5S rRNA to partial denaturation in 4 M urea is conferred by the five oocyte-specific nucleotide substitutions in loop B/helix III. In contrast, the superior stability of oocyte 5S rRNA in the presence of EDTA is related to a single C substitution at position 79.     

Comparison of the nucleotide sequence of soybean 18S rRNA with the sequences of other small-subunit rRNAs

We present the sequence of the nuclear-encoded ribosomal small-subunit RNA from soybean. The soybean 18S rRNA sequence of 1807 nucleotides (nt) is contained in a gene family of approximately 800 closely related members per haploid genome. This sequence is compared with the ribosomal small-subunit RNAs of maize (1805 nt), yeast (1789 nt), Xenopus (1825 nt), rat (1869 nt), and Escherichia coli (1541 nt). Significant sequence homology is observed among the eukaryotic small-subunit rRNAs examined, and some sequence homology is observed between eukaryotic and prokaryotic small-subunit rRNAs. Conserved regions are found to be interspersed among highly diverged sequences. The significance of these comparisons is evaluated using computer simulation of a random sequence model. A tentative model of the secondary structure of soybean 18S rRNA is presented and discussed in the context of the functions of the various conserved regions within the sequence. On the basis of this model, the short base-paired sequences defining the four structural and functional domains of all 18S rRNAs are seen to be well conserved. The potential roles of other conserved soybean 18S rRNA sequences in protein synthesis are discussed.     

Comparison of the nucleotide sequence of soybean 18S rRNA with the sequences of other small-subunit rRNAs

Summary We present the sequence of the nuclearencoded ribosomal small-subunit RNA from soybean. The soybean 18S rRNA sequence of 1807 nucleotides (nt) is contained in a gene family of approximately 800 closely related members per haploid genome. This sequence is compared with the ribosomal small-subunit RNAs of maize (1805 nt), yeast (1789 nt),Xenopus (1825 nt), rat (1869 nt), andEscherichia coli (1541 nt). Significant sequence homology is observed among the eukaryotic small-subunit rRNAs examined, and some sequence homology is observed between eukaryotic and prokaryotic small-subunit rRNAs. Conserved regions are found to be interspersed among highly diverged sequences. The significance of these comparisons is evaluated using computer simulation of a random sequence model. A tentative model of the secondary structure of soybean 18S rRNA is presented and discussed in the context of the functions of the various conserved regions within the sequence. On the basis of this model, the short basepaired sequences defining the four structural and functional domains of all 18S rRNAs are seen to be well conserved. The potential roles of other conserved soybean 18S rRNA sequences in protein synthesis are discussed.     

Functional role of rRNAs in plant translation system tested with antisense strategy.

There are regions in rRNA which are evolutionary conserved and exposed on ribosomal surface. We selected in plant material (Lupinus luteus) two of them: the alpha-sarcin domain of 26S rRNA (L-rRNA) and C loop of 5S rRNA, to be further investigated using antisense oligomers as research tools. We found inhibition of the model polypeptide biosynthesis (up to 80%) due to specific hybridization of oligomers addressed to alpha-sarcin domain and loop C. Based on our results we present the evidence for the key role played by these regions of rRNAs during protein biosynthesis in plant system. According to our hypothesis, conformational changes of these two regions are synchronised and cooperative during transition of pre- to post-translocation state of the ribosome. The correlation of structure and activity of rRNA domains in translation is shown.     

Secondary structure and domain architecture of the 23S and 5S rRNAs

We present a de novo re-determination of the secondary (2°) structure and domain architecture of the 23S and 5S rRNAs, using 3D structures, determined by X-ray diffraction, as input. In the traditional 2° structure, the center of the 23S rRNA is an extended single strand, which in 3D is seen to be compact and double helical. Accurately assigning nucleotides to helices compels a revision of the 23S rRNA 2° structure. Unlike the traditional 2° structure, the revised 2° structure of the 23S rRNA shows architectural similarity with the 16S rRNA. The revised 2° structure also reveals a clear relationship with the 3D structure and is generalizable to rRNAs of other species from all three domains of life. The 2° structure revision required us to reconsider the domain architecture. We partitioned the 23S rRNA into domains through analysis of molecular interactions, calculations of 2D folding propensities and compactness. The best domain model for the 23S rRNA contains seven domains, not six as previously ascribed. Domain 0 forms the core of the 23S rRNA, to which the other six domains are rooted. Editable 2° structures mapped with various data are provided (http://apollo.chemistry.gatech.edu/RibosomeGallery).     

Repeated parallel evolution of minimal rRNAs revealed from detailed comparative analysis

The concept of a minimal ribosomal RNA-containing ribosome, a structure with a minimal set of elements capable of providing protein biosynthesis, is essential for understanding this fundamental cellular process. Nematodes and trypanosomes have minimal mitochondrial rRNAs and detailed reconstructions of their secondary structures indicate that certain conserved helices have been lost in these taxa. In contrast, several recent studies on acariform mites have argued that minimal rRNAs may evolve via shortening of secondary structure elements but not the loss of these elements as shown for trypanosomes and nematodes. Based on extensive structural analysis of chelicerate arthropods, we demonstrate that extremely short rRNAs of acariform mites share certain structural modifications with nematodes and trypanosomes: loss of helices of the GTPase region and divergence in the evolutionarily conserved connecting loop between helices H1648 and H1764 of the large subunit rRNA. These highly concerted parallel modifications indicate that minimal rRNAs were generated under the strong selection that favored or tolerated reductions of helices in particular locations while maintaining the functionality of the rRNA molecules throughout evolution. We also discuss potential evolution of minimal rRNAs and atypical transfer RNAs.     

In vitro incorporation of eubacterial, archaebacterial and eukaryotic 5S rRNAs into large ribosomal subunits of Bacillus stearothermophilus.

Bacillus stearothermophilus large ribosomal subunits were reconstituted in the presence of 5S rRNAs from different origins and tested for their biological activities. The results obtained have shown that eubacterial and archaebacterial 5S rRNAs can easily substitute for B. stearothermophilus 5S rRNA in the reconstitution, while eukaryotic 5S rRNAs yield ribosomal subunits with reduced biological activities. From our results we propose an interaction between nucleotides 42-47 of 5S rRNA and nucleotides 2603-2608 of 23S rRNA during the assembly of the 50S ribosomal subunit. Other experiments with eukaryotic 5.8S rRNAs reveal, if at all, a very low incorporation of these RNA species into the reconstituted ribosomes.     

Sequence and secondary structure of Drosophila melanogaster 5.8S and 2S rRNAs and of the processing site between them.

Drosophila melanogaster 5.8S and 2S rRNAs were end-labeled with 32p at either the 5' or 3' end and were sequenced. 5.8S rRNA is 123 nucleotides long and homologous to the 5' part of sequenced 5.8S molecules from other species. 2S rRNA is 30 nucleotides long and homologous to the 3' part of other 5.8S molecules. The 3' end of the 5.8S molecule is able to base-pair with the 5' end of the 2S rRNA to generate a helical region equivalent in position to the "GC-rich hairpin" found in all previously sequenced 5.8S molecules. Probing the structure of the labeled Drosophila 5.8S molecule with S1 nuclease in solution verifies its similarity to other 5.8S rRNAs. The 2S rRNA is shown to form a stable complex with both 5.8S and 26S rRNAs separately and together. 5.8S rRNA can also form either binary or ternary complexes with 2S and 26S rRNA. It is concluded that the 5.8S rRNA in Drosophila melanogaster is very similar both in sequence and structure to other 5.8 rRNAs but is split into two pieces, the 2S rRNA being the 3' part. 2S anchors the 5.8S and 26S rRNA. The order of the rRNA coding regions in the ribosomal DNA repeating unit is shown to be 18S - 5.8S - 2S - 26S. Direct sequencing of ribosomal DNA shows that the 5.8S and 2S regions are separated by a 28 nucleotide spacer which is A-T rich and is presumably removed by a specific processing event. A secondary structure model is proposed for the 26S-5.8S ternary complex and for the presumptive precursor molecule.     

Deletion of EFL1 results in heterogeneity of the 60 S GTPase-associated rRNA conformation

Previous work suggested that the release of the nucleolar Tif6 from nascent 60 S subunits occurs in the cytoplasm and requires the cytoplasmic EF-2-like GTPase, Efl1. To check whether this release involves an rRNA structural rearrangement mediated by Efl1, we analyzed the rRNA conformation of the GTPase center of 80 S ribosomes in three contexts: wild-type, Deltaefl1 and a dominant suppressor R1 of Deltaefl1. This analysis was restricted to domain II and VI of 25 S rRNA. The rRNA analysis of R1 ribosomes allows us to distinguish the effects due to depletion of Efl1 from the resulting nucleolar deficit of Tif6. Efl1 inhibits the EF-2 GTPase activity, suggesting that the two proteins share a similar ribosome-binding site. The 80 S ribosomes from either type failed to show any difference of conformation in the two rRNA domains analyzed. However, the same analysis performed on the pool of free 60 S subunits reveals several rRNA conformational differences between wild-type and Deltaefl1 subunits, whereas that from the suppressor strain is similar to wild-type. This suggests that the nucleolar deficit of Tif6 during assembly of the 60 S preribosomes is responsible for the changes in rRNA conformation observed in Deltaefl1 60 S subunits. We also purified 60 S preribosomes from the three genetic contexts by TAP-tagging Tif6. The protein content of 60 S preribosomes associated with Tif6p in a Deltaefl1 strain are obtained at a lower yield but have, surprisingly, a protein composition that is a priori similar to that of wild-type and the suppressor strain.     

Structural analysis of prokaryotic and eukaryotic 5S rRNAs by RNase H

The availabilities of single-stranded 5S rRNA regions c, d and d' for base pairing interactions were analyzed by using synthetic DNA oligomers. Hybrid formation was detected by the endonucleolytical mode of the RNA-DNA specific action of RNase H. Provided that the hybrid interaction involved 6 successive base pairs, 5S rRNA loop c nucleotides 42-47 displayed accessibility in Escherichia coli, Bacillus stearothermophilus and Thermus thermophilus 5S rRNAs as well as in eukaryotic 5S rRNAs from Saccharomyces carlsbergensis, Rattus rattus and Equisetum arvense. Investigating eubacterial 5S rRNA regions d and d' (nucleotides 71-76 and 99-105, respectively), susceptibility was observed in E. coli 5S rRNA which, however, decreases in B. stearothermophilus and even more so in T. thermophilus 5S rRNA. For additional evaluation of the data obtained by RNase H cleavage, association constants of the hexanucleotides were determined by equilibrium dialysis at 4 degrees C for B. stearothermophilus 5S rRNA. The results obtained reveal that nucleotides 36-41 of B. stearothermophilus 5S rRNA are inaccessible for Watson-Crick interaction, which suggests that this part of loop c is in a structurally constrained configuration, or buried in the tertiary structure or involved in tertiary interactions.     

Secondary structure and phylogeny of Staphylococcus and Micrococcus 5S rRNAs

Nucleotide sequences of 5S rRNAs from four bacteria, Staphylococcus aureus Smith (diffuse), Staphylococcus epidermidis ATCC 14990, Micrococcus luteus ATCC 9341 and Micrococcus luteus ATCC 4698, were determined. The secondary structural models of S. aureus and S. epidermidis sequences showed characteristics of the gram-positive bacterial 5S rRNA (116-N type [H. Hori and S. Osawa, Proc. Natl. Acad. Sci. U.S.A. 76:381-385, 1979]). Those of M. luteus ATCC 9341 and M. luteus ATCC 4698 together with that of Streptomyces griseus (A. Simoncsits, Nucleic Acids Res. 8:4111-4124, 1980) showed intermediary characteristics between the gram-positive and gram-negative (120-N type [H. Hori and S. Osawa, 1979]) 5S rRNAs. This and previous studies revealed that there exist at least three major groups of eubacteria having distinct 5S rRNA and belonging to different stems in the 5S rRNA phylogenic tree.