Evolutionary conservation of structure and function of high molecular weight ribosomal RNA

Fossil-calibrated phylogenies of Southern cave wētā show dispersal and extinction confound biogeographic signal

The biota of continents and islands are commonly considered to have a source-sink relationship, but small islands can harbour distinctive taxa. The distribution of four monotypic genera of Orthoptera on young subantarctic islands indicates a role for long-distance dispersal and extinction. Phylogenetic relationships were inferred from whole mtDNA genomes and nuclear sequences (45S cassette; four histones). We used a fossil and one palaeogeographic event to calibrate molecular clock analysis. We confirm that neither the Australian nor Aotearoa-New Zealand Rhaphidophoridae faunas are monophyletic. The radiation of Macropathinae may have begun in the late Jurassic, but trans-oceanic dispersal is required to explain the current distribution of some lineages within this subfamily. Dating the most recent common ancestor of seven island endemic species with their nearest mainland relative suggests that each existed long before their island home was available. Time estimates from our fossil-calibrated molecular clock analysis suggest several lineages have not been detected on mainland New Zealand, Australia, or elsewhere most probably due to their extinction, providing evidence that patterns of extinction, which are not consistently linked to range size or lineage age, confound biogeographic signal.

Mwhitiwhiti Aotearoa: Phylogeny and synonymy of the silent alpine grasshopper radiation of New Zealand (Orthoptera: Acrididae)

Aotearoa New Zealand has a fauna of endemic alpine grasshoppers, consisting of thirteen species distributed among four genera. The many re-classifications of species within this group and the presence of species complexes highlight the uncertainty that surrounds relationships within and between these genera. High-throughput Next Generation Sequencing was used to assemble the complete mitochondrial genomes, 45S ribosomal cassettes and histone sequences of New Zealands four endemic alpine genera: Alpinacris, Brachaspis, Paprides and Sigaus. Phylogenetic analysis of these molecular datasets, as individual genes, partitions and combinations returned a consistent topology that is incompatible with the current classification. The genera Sigaus, Alpinacris, and Paprides all exhibit paraphyly. A consideration of the pronotum, epiphallus and terminalia of adult specimens reveals species-specific differences, but fails to provide compelling evidence for species groups justifying distinct genera. In combination with phylogenetic, morphological and spatial evidence we propose a simplified taxonomy consisting of a single genus for the mwhitiwhiti Aotearoa species radiation.

TOPAS: network-based structural alignment of RNA sequences

MOTIVATION

For many RNA families, the secondary structure is known to be better conserved among the member RNAs compared to the primary sequence. For this reason, it is important to consider the underlying folding structures when aligning RNA sequences, especially for those with relatively low sequence identity. Given a set of RNAs with unknown structures, simultaneous RNA alignment and folding algorithms aim to accurately align the RNAs by jointly predicting their consensus secondary structure and the optimal sequence alignment. Despite the improved accuracy of the resulting alignment, the computational complexity of simultaneous alignment and folding for a pair of RNAs is O(N6), which is too costly to be used for large-scale analysis.

RESULTS

In order to address this shortcoming, in this work, we propose a novel network-based scheme for pairwise structural alignment of RNAs. The proposed algorithm, TOPAS, builds on the concept of topological networks that provide structural maps of the RNAs to be aligned. For each RNA sequence, TOPAS first constructs a topological network based on the predicted folding structure, which consists of sequential edges and structural edges weighted by the base-pairing probabilities. The obtained networks can then be efficiently aligned by using probabilistic network alignment techniques, thereby yielding the structural alignment of the RNAs. The computational complexity of our proposed method is significantly lower than that of the Sankoff-style dynamic programming approach, while yielding favorable alignment results. Furthermore, another important advantage of the proposed algorithm is its capability of handling RNAs with pseudoknots while predicting the RNA structural alignment. We demonstrate that TOPAS generally outperforms previous RNA structural alignment methods on RNA benchmarks in terms of both speed and accuracy.

AVAILABILITY AND IMPLEMENTATION

Source code of TOPAS and the benchmark data used in this paper are available at https://github.com/bjyoontamu/TOPAS.

Post-transcriptional regulation by cytosine-5 methylation of RNA

The recent advent of high-throughput sequencing technologies coupled with RNA modifications detection methods has allowed the detection of RNA modifications at single nucleotide resolution giving a more comprehensive landscape of post-transcriptional gene regulation pathways. In this review, we focus on the occurrence of 5-methylcytosine (m⁵C) in the transcriptome. We summarise the main findings of the molecular role in post-transcriptional regulation that governs m⁵C deposition in RNAs. Functionally, m⁵C deposition can regulate several cellular and physiological processes including development, differentiation and survival to stress stimuli. Despite many aspects concerning m⁵C deposition in RNA, such as position or sequence context and the fact that many readers and erasers still remain elusive, the overall recent findings indicate that RNA cytosine methylation is a powerful mechanism to post-transcriptionally regulate physiological processes. In addition, mutations in RNA cytosine-5 methyltransferases are associated to pathological processes ranging from neurological syndromes to cancer.

Comparative RNA function analysis reveals high functional similarity between distantly related bacterial 16 S rRNAs

The 16 S rRNA sequence has long been used uncritically as a molecular clock to infer phylogenetic relationships among prokaryotes without fully elucidating the evolutionary changes that this molecule undergoes. In this study, we investigated the functional evolvability of 16 S rRNA, using comparative RNA function analyses between the 16 S rRNAs of Escherichia coli (Proteobacteria) and Acidobacteria (78% identity, 334 nucleotide differences) in the common genetic background of E. coli. While the growth phenotype of an E. coli mutant harboring the acidobacterial gene was disrupted significantly, it was restored almost completely following introduction of a 16 S rRNA sequence with a single base-pair variation in helix 44; the remaining 332 nucleotides were thus functionally similar to those of E. coli. Our results suggest that 16 S rRNAs share an inflexible cradle structure formed by ribosomal proteins and have evolved by accumulating species-specific yet functionally similar mutations. While this experimental evidence suggests the neutral evolvability of 16 S rRNA genes and hence satisfies the necessary requirements to use the sequence as a molecular clock, it also implies the promiscuous nature of the 16 S rRNA gene, i.e., the occurrence of horizontal gene transfer among bacteria.

Yessotoxin triggers ribotoxic stress

This work tests the hypothesis that the marine algal toxin yessotoxin (YTX) can trigger ribotoxic stress response in L6 and BC3H1 myoblast cells. YTX exposure at a concentration of 100 nM displays the characteristics of a ribotoxic stress response in such cells. The exposure leads to activation of the p38 mitogen-activated protein kinase, the stress-activated protein kinase c-jun, and the double-stranded RNA-activated protein kinase (PKR). YTX treatment also causes ribosomal RNA cleavage and inhibits protein synthesis. These observations support the idea that YTX can act as a ribotoxin.

Characterization of active ribosomal RNA harboring MITEs insertion in microsporidian Nosema bombycis genome

Microsporidia are a group of obligate intracellular parasites of medical and agricultural importance, which can infect almost all animals, including human beings. Using the genome data of Nosema bombycis, four families of miniature inverted-repeat transposable elements (MITEs) in ribosomal DNA (rDNA) were characterized in the microsporidian N. bombycis and were named LSUME1, ITSME1, SSUME1, and SSUME2, respectively. The genome-wide investigation of these MITEs shows that these MITEs families distribute randomly in N. bombycis genome. All insertion sequences have conserved characteristics of MITEs, the direct repeat sequence and terminal inverted-repeat sequence at both ends of each MITEs sequence. Additionally, using the CLC RNA Workbench Software, secondary structures of rRNA containing MITEs sequence have been predicted and were located in variable region or expansion segment. Furthermore, using two different probes, one is prepared by MITE sequence only (short probe) and the other is prepared by MITE sequence flanking partial rDNA sequence (long probe); northern blotting and dot blotting have been performed to detect the transcriptional and functional activity of the rDNA containing MITEs insertion. Fortunately, we found that the rDNA, which harbors the MITE, not only can be transcripted but also can form a complete ribosome. This is an interesting thing that one gene can keep active even when it has been inserted with another sequence. But the biological and structural significance of this observation is not readily apparent.

Identification and characterization of the Thermus thermophilus 5-methylcytidine (m5C) methyltransferase modifying 23 S ribosomal RNA (rRNA) base C1942

Methylation of cytidines at carbon-5 is a common posttranscriptional RNA modification encountered across all domains of life. Here, we characterize the modifications of C1942 and C1962 in Thermus thermophilus 23 S rRNA as 5-methylcytidines (m(5)C) and identify the two associated methyltransferases. The methyltransferase modifying C1942, named RlmO, has not been characterized previously. RlmO modifies naked 23 S rRNA, but not the assembled 50 S subunit or 70 S ribosomes. The x-ray crystal structure of this enzyme in complex with the S-adenosyl-l-methionine cofactor at 1.7 Å resolution confirms that RlmO is structurally related to other m(5)C rRNA methyltransferases. Key residues in the active site are located similar to the further distant 5-methyluridine methyltransferase RlmD, suggestive of a similar enzymatic mechanism. RlmO homologues are primarily found in mesophilic bacteria related to T. thermophilus. In accordance, we find that growth of the T. thermophilus strain with an inactivated C1942 methyltransferase gene is not compromised at non-optimal temperatures.

The Era GTPase recognizes the GAUCACCUCC sequence and binds helix 45 near the 3' end of 16S rRNA

Era, composed of a GTPase domain and a K homology domain, is essential for bacterial cell viability. It is required for the maturation of 16S rRNA and assembly of the 30S ribosomal subunit. We showed previously that the protein recognizes nine nucleotides ( ) near the 3' end of 16S rRNA, and that this recognition stimulates GTP-hydrolyzing activity of Era. In all three kingdoms of life, the sequence and helix 45 (h45) (nucleotides 1506-1529) are highly conserved. It has been shown that the to double mutation severely affects the viability of bacteria. However, whether Era interacts with G1530 and/or h45 and whether such interactions (if any) contribute to the stimulation of Era's GTPase activity were not known. Here, we report two RNA structures that contain nucleotides 1506-1542 (RNA301), one in complex with Era and GDPNP (GNP), a nonhydrolysable GTP-analogue, and the other in complex with Era, GNP, and the KsgA methyltransferase. The structures show that Era recognizes 10 nucleotides, including G1530, and that Era also binds h45. Moreover, GTPase assay experiments show that G1530 does not stimulate Era's GTPase activity. Rather, A1531 and A1534 are most important for stimulation and h45 further contributes to the stimulation. Although G1530 does not contribute to the intrinsic GTPase activity of Era, its interaction with Era is important for binding and is essential for the protein to function, leading to the discovery of a new cold-sensitive phenotype of Era.

Novel microbial diversity retrieved by autonomous robotic exploration of the world's deepest vertical phreatic sinkhole

The deep phreatic thermal explorer (DEPTHX) is an autonomous underwater vehicle designed to navigate an unexplored environment, generate high-resolution three-dimensional (3-D) maps, collect biological samples based on an autonomous sampling decision, and return to its origin. In the spring of 2007, DEPTHX was deployed in Zacatón, a deep (approximately 318 m), limestone, phreatic sinkhole (cenote) in northeastern Mexico. As DEPTHX descended, it generated a 3-D map based on the processing of range data from 54 onboard sonars. The vehicle collected water column samples and wall biomat samples throughout the depth profile of the cenote. Post-expedition sample analysis via comparative analysis of 16S rRNA gene sequences revealed a wealth of microbial diversity. Traditional Sanger gene sequencing combined with a barcoded-amplicon pyrosequencing approach revealed novel, phylum-level lineages from the domains Bacteria and Archaea; in addition, several novel subphylum lineages were also identified. Overall, DEPTHX successfully navigated and mapped Zacatón, and collected biological samples based on an autonomous decision, which revealed novel microbial diversity in a previously unexplored environment.

Characterization of internal transcribed spacer (ITS1)-ITS2 region of ribosomal RNA gene from 25 species of Culicoides biting midges (Diptera: Ceratopogonidae) in Japan

We determined nucleotide sequences of the nuclear rDNA internal transcribed spacer (ITS)1-5.8S-ITS2a-2S-ITS2 region in 103 individuals of 25 Culicoides species (Diptera: Ceratopogonidae) from 11 locations in Japan. Ribosomal RNA genes, 5.8S and 2S rDNA, were highly conserved among the species with few variations. The ITS2a region showed length variation among species. Both ITS1 and ITS2 showed highly varied sequences among species. The noticeable indel regions among ITS1 sequences are present in some Culicoides species, separating species into two types having long or short ITS1 region. However, Culicoides cylindratus Kitaoka possesses both types of ITS1 in each individual; these results seem to indicate that the ITS1-long type was the prototype and the short type was produced through deletion in many Culicoides species. One species, belonging to subgenus Avaritia, possessed an Avaritia-specific sequence in ITS1 and phylogenetically formed a monophyletic group. Geographical genotypes in a species were not clear. Species-specific sequence features were observed, enabling molecular identification of Culicoides species.

Structure of ERA in complex with the 3' end of 16S rRNA: implications for ribosome biogenesis

ERA, composed of an N-terminal GTPase domain followed by an RNA-binding KH domain, is essential for bacterial cell viability. It binds to 16S rRNA and the 30S ribosomal subunit. However, its RNA-binding site, the functional relationship between the two domains, and its role in ribosome biogenesis remain unclear. We have determined two crystal structures of ERA, a binary complex with GDP and a ternary complex with a GTP-analog and the 1531AUCACCUCCUUA1542 sequence at the 3' end of 16S rRNA. In the ternary complex, the first nine of the 12 nucleotides are recognized by the protein. We show that GTP binding is a prerequisite for RNA recognition by ERA and that RNA recognition stimulates its GTP-hydrolyzing activity. Based on these and other data, we propose a functional cycle of ERA, suggesting that the protein serves as a chaperone for processing and maturation of 16S rRNA and a checkpoint for assembly of the 30S ribosomal subunit. The AUCA sequence is highly conserved among bacteria, archaea, and eukaryotes, whereas the CCUCC, known as the anti-Shine-Dalgarno sequence, is conserved in noneukaryotes only. Therefore, these data suggest a common mechanism for a highly conserved ERA function in all three kingdoms of life by recognizing the AUCA, with a "twist" for noneukaryotic ERA proteins by also recognizing the CCUCC.

Different mechanisms for pseudouridine formation in yeast 5S and 5.8S rRNAs

The selection of sites for pseudouridylation in eukaryotic cytoplasmic rRNA occurs by the base pairing of the rRNA with specific guide sequences within the RNA components of box H/ACA small nucleolar ribonucleoproteins (snoRNPs). Forty-four of the 46 pseudouridines (Psis) in the cytoplasmic rRNA of Saccharomyces cerevisiae have been assigned to guide snoRNAs. Here, we examine the mechanism of Psi formation in 5S and 5.8S rRNA in which the unassigned Psis occur. We show that while the formation of the Psi in 5.8S rRNA is associated with snoRNP activity, the pseudouridylation of 5S rRNA is not. The position of the Psi in 5.8S rRNA is guided by snoRNA snR43 by using conserved sequence elements that also function to guide pseudouridylation elsewhere in the large-subunit rRNA; an internal stem-loop that is not part of typical yeast snoRNAs also is conserved in snR43. The multisubstrate synthase Pus7 catalyzes the formation of the Psi in 5S rRNA at a site that conforms to the 7-nucleotide consensus sequence present in other substrates of Pus7. The different mechanisms involved in 5S and 5.8S rRNA pseudouridylation, as well as the multiple specificities of the individual trans factors concerned, suggest possible roles in linking ribosome production to other processes, such as splicing and tRNA synthesis.

Functional conservation between structurally diverse ribosomal proteins from Drosophila melanogaster and Saccharomyces cerevisiae: fly L23a can substitute for yeast L25 in ribosome assembly and function

The proposed Drosophila melanogaster L23a ribosomal protein features a conserved C-terminal amino acid signature characteristic of other L23a family members and a unique N-terminal extension [Koyama et al. (Poly(ADP-ribose) polymerase interacts with novel Drosophila ribosomal proteins, L22 and l23a, with unique histone-like amino-terminal extensions. Gene 1999; 226: 339-345)], absent from Saccharomyces cerevisiae L25 that nearly doubles the size of fly L23a. The ability of fly L23a to replace the role of yeast L25 in ribosome biogenesis was determined by creating a yeast strain carrying an L25 chromosomal gene disruption and a plasmid-encoded FLAG-tagged L23a gene. Though affected by a reduced growth rate, the strain is dependent on fly L23a-FLAG function for survival and growth, demonstrating functional compatibility between the fly and yeast proteins. Pulse-chase experiments reveal a delay in rRNA processing kinetics, most notably at a late cleavage step that converts precursor 27S rRNA into mature 25S rRNA, likely contributing to the strain's slower growth pattern. Yet, given the essential requirement for L23(a)/L25 in ribosome biogenesis, there is a remarkable tolerance for accommodating the fly L23a N-terminal extension within the structure of the yeast ribosome. A search of available databases shows that the unique N-terminal extension is shared by multiple insect lineages. An evolutionary perspective on L23a structure and function within insect lineages is discussed.

Characterization of the 16S-23S internal transcribed spacer among 34 higher plants: suitability for interspecific plastid transformation

Biomanufacturing by chloroplast transgene expression has the potential to produce significant amounts of biopharmaceuticals, endow plants with novel commercial or humanitarian capabilities, enhance phytoremediation methods and harden plants against adverse environments. Plastid bioengineering exploits the phenomenon of homologous recombination to specifically integrate heterologous sequences into the plastid genome. Previous research suggests the plastid genome 16S-23S internal transcribed spacer provides an advantageous integration site for transgene expression. To characterize the suitability of the 16S-23S region for interspecific recombination, we developed primers against conserved plastid sequences and amplified approximately 2.6 kb from 25 plant species. We analyzed the amplicons with nine species from Genbank for homeology, phylogenetic relationships, potential to form chimeric rDNA elements disruptive to translational/replication systems, and the potential number of recombination events for various minimal essential processing segments (MEPS) lengths. Multiple sequence alignment of the 34 species revealed considerable conservation, with identities exceeding 95% among the angiosperms. Substitutions were statistically clustered, generally in noncoding sites, although proposed functional elements such as the OriA region and 3' terminus of the 16S rRNA exhibited unexpected variation. The nonrandom distribution of substitutions undermines the established, statistical method of estimating the number of recombination initiation sites. This finding is further substantiated by comparing statistical estimates of the number of MEPS sites to a direct count at three different MEPS lengths. We frame this in silico analysis in terms of the potential of the 16S-23S region as a target for interspecific transformation, and describe a 'primer-to-plastid' system to rapidly generate species-specific flanking regions for transformation vectors.

Polarity and hydrophobicity interactions in protein synthesis process

About 30 years ago, experiments found that there are polarity and hydrophobicity (P and H) correlations and affinity between amino acids and their anticodons. Although it is shown that these experimental findings are important for explaining the origins of the genetic code, the great potential of P and H interactions in investigating other bio-functions have not been fully explored. Here, through raising, discussing and answering seven relevant questions hidden in tRNA aminoacylation, the formation of peptide bonds, and the ending of translations, etc., we show our theoretical findings that the P and H correlations and affinity take vital roles in the protein synthesis process. We found the relationship between the 3' end ACCN sequences of tRNA molecules and the activated amino acids and its biological significance, the rRNAs' consensus sequences 5'NCC...TGG3' or 5'TGG...NCC3' which may perform as functional segments of rRNAs to help triggering the reaction of peptide formation, and common nature of releasing factors that the first amino acid residue of releasing factors ERF, RF1 and RF2 are all Methionine, except a few Alanine, which may be necessary for releasing the translated polypeptide and stopping the translating process. In the terms of P and H correlations and affinity, we provide explanations of why only using the poly (G) as mRNA template cannot get the poly (Gly) in experiments deciphering the genetic code, why Gly often appears in beta turns and why translational bypassing might occur when translating 5'GGAUGA on mRNA. Since amino acids and nucleotides are the subunits, respectively, for composing proteins and nucleic acids, these findings will help in further understanding interactions among the bio-macromolecules. These findings are also helpful for investigating rRNAs, further understanding the protein synthesis process and analysing similar bio-problems, and should be proved useful for experimental biologists.

Interaction of Era with the 30S ribosomal subunit implications for 30S subunit assembly

Era (E. coliRas-like protein) is a highly conserved and essential GTPase in bacteria. It binds to the 16S ribosomal RNA (rRNA) of the small (30S) ribosomal subunit, and its depletion leads to accumulation of an unprocessed precursor of the 16S rRNA. We have obtained a three-dimensional cryo-electron microscopic map of the Thermus thermophilus 30S-Era complex. Era binds in the cleft between the head and platform of the 30S subunit and locks the subunit in a conformation that is not favorable for association with the large (50S) ribosomal subunit. The RNA binding KH motif present within the C-terminal domain of Era interacts with the conserved nucleotides in the 3' region of the 16S rRNA. Furthermore, Era makes contact with several assembly elements of the 30S subunit. These observations suggest a direct involvement of Era in the assembly and maturation of the 30S subunit.

The complete set of H/ACA snoRNAs that guide rRNA pseudouridylations in Saccharomyces cerevisiae

Conversion of uridines into pseudouridines (Psis) is the most frequent base modification in ribosomal RNAs (rRNAs). In eukaryotes, the pseudouridylation sites are specified by base-pairing with specific target sequences within H/ACA small nucleolar RNAs (snoRNAs). The yeast rRNAs harbor 44 Psis, but, when this work began, 15 Psis had completely unknown guide snoRNAs. This suggested that many snoRNAs remained to be discovered. To address this problem and further complete the snoRNA assignment to Psi sites, we identified the complete set of RNAs associated with the H/ACA snoRNP specific proteins Gar1p and Nhp2p by coupling TAP-tag purifications with genomic DNA microarrays experiments. Surprisingly, while we identified all the previously known H/ACA snoRNAs, we selected only three new snoRNAs. This suggested that most of the missing Psi guides were present in previously known snoRNAs but had been overlooked. We confirmed this hypothesis by systematically investigating the role of previously known, as well as of the newly identified snoRNAs, in specifying rRNA Psi sites and found all but one missing guide RNAs. During the completion of this work, another study, based on bioinformatic predictions, also reported the identification of most missing guide RNAs. Altogether, all Psi guides are now identified and we can tell that, in budding yeast, the 44 Psis are guided by 28 snoRNAs. Finally, aside from snR30, an atypical small RNA of heterogeneous length and at least one mRNA, all Gar1p and Nhp2p associated RNAs characterized by our work turned out to be snoRNAs involved in rRNA Psi specification.

The first complete monogenean ribosomal RNA gene operon: sequence and secondary structure of the Gyrodactylus salaris Malmberg, 1957, large subunit ribosomal RNA gene

The sequence of the Gyrodactylus salaris Malmberg, 1957, large subunit, or 28S, ribosomal RNA (rRNA) gene has been determined. This gene is the final portion of the Gyrodactylus rRNA gene operon to be sequenced and results in the first complete sequence of all rRNA genes and spacers from a monogenean. The nucleotide sequence was used to predict the secondary structure of the large subunit rRNA, and regions of conserved and variable sequence and structure were identified. The site where the 5' terminus of the 5.8S rRNA binds to a region within the large subunit rRNA was predicted and complements the anticipated interaction of the 3' terminus of the 5.8S with the 5' terminus of the large subunit rRNA. The large subunit gene may be useful in phylogenetic analysis of the Monogenea or Platyhelminthes and comparisons with other eukaryotes. The variable domains C and H may be most suitable for this purpose.

Expression and purification of human ribosomal proteins S3, S5, S10, S19, and S26

The cDNAs for the human ribosomal proteins S3, S5, S10, S19, and S26 were introduced into a pET-15b vector and recombinant proteins containing an N-(His)(6)-fusion tag were expressed in high yields. To resolve the problem of frameshift during expression of S26 caused by the presence of tandem arginine codons in its mRNA that are rare in Escherichia coli, we substituted the rare AGA codon with the more frequent arginine codon (CGC) using a primer with this mutation for PCR amplification of S26 cDNA. All proteins were expressed mainly in the form of inclusion bodies and purified to homogeneity by metal affinity chromatography in one step (except for S3). Expression of the full-length S3 was accompanied by the formation of a low molecular weight polypeptide that was co-purified with S3 by metal affinity chromatography. Complete purification of S3 required an additional gel-filtration step. The proteins were refolded by stepwise dialysis. Both identity and purity of the proteins were confirmed by 2D PAGE. The proteins obtained could be used in a wide range of applications in biophysics, biochemistry, and molecular biology.

Evolution and phylogenetic information content of the ribosomal DNA repeat unit in the Blattodea (Insecta)

The organization, structure, and nucleotide variability of the ribosomal repeat unit was compared among families, genera, and species of cockroaches (Insecta:Blattodea). Sequence comparisons and molecular phylogenetic analyses were used to describe rDNA repeat unit variation at differing taxonomic levels. A reverse similar 1200 bp fragment of the 28S rDNA sequence was assessed for its potential utility in reconstructing higher-level phylogenetic relationships in cockroaches. Parsimony and maximum likelihood analyses of these data strongly support the expected pattern of relationships among cockroach groups. The examined 5' end of the 28S rDNA is shown to be an informative marker for larger studies of cockroach phylogeny. Comparative analysis of the nucleotide sequences of the rDNA internal transcribed spacers (ITS1 and ITS2) among closely related species of Blattella and Periplaneta reveals that ITS sequences can vary widely in primary sequence, length, and folding pattern. Secondary structure estimates for the ITS region of Blattella species indicate that variation in this spacer region can also influence the folding pattern of the 5.8S subunit. These results support the idea that ITS sequences play an important role in the stability and function of the rRNA cluster.

RNA:(guanine-N2) methyltransferases RsmC/RsmD and their homologs revisited--bioinformatic analysis and prediction of the active site based on the uncharacterized Mj0882 protein structure

BACKGROUND

Escherichia coli guanine-N2 (m2G) methyltransferases (MTases) RsmC and RsmD modify nucleosides G1207 and G966 of 16S rRNA. They possess a common MTase domain in the C-terminus and a variable region in the N-terminus. Their C-terminal domain is related to the YbiN family of hypothetical MTases, but nothing is known about the structure or function of the N-terminal domain.

RESULTS

Using a combination of sequence database searches and fold recognition methods it has been demonstrated that the N-termini of RsmC and RsmD are related to each other and that they represent a "degenerated" version of the C-terminal MTase domain. Novel members of the YbiN family from Archaea and Eukaryota were also indentified. It is inferred that YbiN and both domains of RsmC and RsmD are closely related to a family of putative MTases from Gram-positive bacteria and Archaea, typified by the Mj0882 protein from M. jannaschii (1dus in PDB). Based on the results of sequence analysis and structure prediction, the residues involved in cofactor binding, target recognition and catalysis were identified, and the mechanism of the guanine-N2 methyltransfer reaction was proposed.

CONCLUSIONS

Using the known Mj0882 structure, a comprehensive analysis of sequence-structure-function relationships in the family of genuine and putative m2G MTases was performed. The results provide novel insight into the mechanism of m2G methylation and will serve as a platform for experimental analysis of numerous uncharacterized N-MTases.

Trichinella spiralis mtDNA: a nematode mitochondrial genome that encodes a putative ATP8 and normally structured tRNAS and has a gene arrangement relatable to those of coelomate metazoans

The complete mitochondrial DNA (mtDNA) of the nematode Trichinella spiralis has been amplified in four overlapping fragments and 16,656 bp of its sequence has been determined. This sequence contains the 37 genes typical of metazoan mtDNAs, including a putative atp8, which is absent from all other nematode mtDNAs examined. The genes are transcribed from both mtDNA strands and have an arrangement relatable to those of coelomate metazoans, but not to those of secernentean nematodes. All protein genes appear to initiate with ATN codons, typical for metazoans. Neither TTG nor GTT start codons, inferred for several genes of other nematodes, were found. The 22 T. spiralis tRNA genes fall into three categories: (i) those with the potential to form conventional "cloverleaf" secondary structures, (ii) those with TPsiC arm + variable arm replacement loops, and (iii) those with DHU-arm replacement loops. Mt-tRNA(R) has a 5'-UCG-3' anticodon, as in most other metazoans, instead of the very unusual 5'-ACG-3' present in the secernentean nematodes. The sequence also contains a large repeat region that is polymorphic in size at the population and/or individual level.

Conformational changes induced in the Saccharomyces cerevisiae GTPase-associated rRNA by ribosomal stalk components and a translocation inhibitor

The yeast ribosomal GTPase associated center is made of parts of the 26S rRNA domains II and VI, and a number of proteins including P0, P1alpha, P1beta, P2alpha, P2beta and L12. Mapping of the rRNA neighborhood of the proteins was performed by footprinting in ribosomes from yeast strains lacking different GTPase components. The absence of protein P0 dramatically increases the sensitivity of the defective ribosome to degradation hampering the RNA footprinting. In ribosomes lacking the P1/P2 complex, protection of a number of nucleotides is detected around positions 840, 880, 1100, 1220-1280 and 1350 in domain II as well as in several positions in the domain VI alpha-sarcin region. The protection pattern resembles the one reported for the interaction of elongation factors in bacterial systems. The results exclude a direct interaction of these proteins with the rRNA and are compatible with an increase in the ribosome affinity for EF-2 in the absence of the acidic P proteins. Interestingly, a sordarin derivative inhibitor of EF-2 causes an opposite effect, increasing the reactivity in positions protected by the absence of P1/P2. Similarly, a deficiency in protein L12 exposes nucleotides G1235, G1242, A1262, A1269, A1270 and A1272 to chemical modification, thus situating the protein binding site in the most conserved part of the 26S rRNA, equivalent to the bacterial protein L11 binding site.

The plastid ribosomal proteins. Identification of all the proteins in the 50 S subunit of an organelle ribosome (chloroplast)

We have completed identification of all the ribosomal proteins (RPs) in spinach plastid (chloroplast) ribosomal 50 S subunit via a proteomic approach using two-dimensional electrophoresis, electroblotting/protein sequencing, high performance liquid chromatography purification, polymerase chain reaction-based screening of cDNA library/nucleotide sequencing, and mass spectrometry (reversed-phase HPLC coupled to electrospray ionization mass spectrometry and electrospray ionization mass spectrometry). Spinach plastid 50 S subunit comprises 33 proteins, of which 31 are orthologues of Escherichia coli RPs and two are plastid-specific RPs (PSRP-5 and PSRP-6) having no homologues in other types of ribosomes. Orthologues of E. coli L25 and L30 are absent in spinach plastid ribosome. 25 of the plastid 50 S RPs are encoded in the nuclear genome and synthesized on cytosolic ribosomes, whereas eight of the plastid RPs are encoded in the plastid organelle genome and synthesized on plastid ribosomes. Sites for transit peptide cleavages in the cytosolic RP precursors and formyl Met processing in the plastid-synthesized RPs were established. Post-translational modifications were observed in several mature plastid RPs, including multiple forms of L10, L18, L31, and PSRP-5 and N-terminal/internal modifications in L2, L11 and L16. Comparison of the RPs in gradient-purified 70 S ribosome with those in the 30 and 50 S subunits revealed an additional protein, in approximately stoichiometric amount, specific to the 70 S ribosome. It was identified to be plastid ribosome recycling factor. Combining with our recent study of the proteins in plastid 30 S subunit (Yamaguchi, K., von Knoblauch, K., and Subramanian, A. R. (2000) J. Biol. Chem. 275, 28455-28465), we show that spinach plastid ribosome comprises 59 proteins (33 in 50 S subunit and 25 in 30 S subunit and ribosome recycling factor in 70 S), of which 53 are E. coli orthologues and 6 are plastid-specific proteins (PSRP-1 to PSRP-6). We propose the hypothesis that PSRPs were evolved to perform functions unique to plastid translation and its regulation, including protein targeting/translocation to thylakoid membrane via plastid 50 S subunit.

The interrelationships of proseriata (Platyhelminthes: seriata) tested with molecules and morphology

Proseriate flatworms are common members of the interstitial benthic fauna worldwide, predominantly occupying marine environments. As minute animals, having relatively few characters useful for cladistic analysis, they have been difficult to present in a phylogenetic framework using morphology alone. Here we present a new morphological matrix consisting of 16 putatively homologous characters and two molecular data sets to investigate further this major group of free-living members of the Platyhelminthes. Complete 18S rDNA (representing 277 parsimony-informative characters) from 17 ingroup taxa and partial 28S rDNA spanning variable expansion regions D1 to D3 and D1 to D6 (representing 219 and 361 parsimony-informative characters, respectively) from 27 and 14 ingroup taxa, respectively, were determined and aligned as complementary data sets. Morphological and molecular data sets were analyzed separately and together to determine underlying phylogenetic patterns and to resolve conflict between published scenarios based on morphology alone. The monophyly of the Proseriata cannot be confirmed categorically with any of these data sets. However, the constituent taxa are confirmed as basal members of the Neoophora, and a sister group relationship with Tricladida is rejected. Similarly, the monophyly of one of the two subtaxa of the Proseriata, the Lithophora, could not be confirmed with molecules. Concerning intragroup relationships, we could reject one of the two phylogenetic trees formerly proposed, as well as the clade Otoplanidae + Coelogynoporidae. However, a clade Otoplanidae + Archimonocelididae + Monocelididae (to which the Monotoplanidae belong) was supported, and the position of the genus Calviria shifted from the Archimonocelididae to the Coelogynoporidae.

Identities and phylogenetic comparisons of posttranscriptional modifications in 16 S ribosomal RNA from Haloferax volcanii

Small subunit (16 S) rRNA from the archaeon Haloferax volcanii, for which sites of modification were previously reported, was examined using mass spectrometry. A census of all modified residues was taken by liquid chromatography/electrospray ionization-mass spectrometry analysis of a total nucleoside digest of the rRNA. Following rRNA hydrolysis by RNase T(1), accurate molecular mass values of oligonucleotide products were measured using liquid chromatography/electrospray ionization-mass spectrometry and compared with values predicted from the corresponding gene sequence. Three modified nucleosides, distributed over four conserved sites in the decoding region of the molecule, were characterized: 3-(3-amino-3-carboxypropyl)uridine-966, N(6)-methyladenosine-1501, and N(6),N(6)-dimethyladenosine-1518 and -1519 (all Escherichia coli numbering). Nucleoside 3-(3-amino-3-carboxypropyl)uridine, previously unknown in rRNA, occurs at a highly conserved site of modification in all three evolutionary domains but for which no structural assignment in archaea has been previously reported. Nucleoside N(6)-methyladenosine, not previously placed in archaeal rRNAs, frequently occurs at the analogous location in eukaryotic small subunit rRNA but not in bacteria. H. volcanii small subunit rRNA appears to reflect the phenotypically low modification level in the Crenarchaeota kingdom and is the only cytoplasmic small subunit rRNA shown to lack pseudouridine.

Caspase-dependent cleavage of nucleic acids

Autoimmune diseases are frequently characterized by the presence of autoantibodies directed against nucleic acid-protein complexes present in the nucleus of the cell. The mechanisms by which these autoantigenic molecules escape immunological tolerance are largely unknown, although a number of recent observations suggest that modified self-proteins generated during apoptosis may play an important role in the development of autoimmunity. It has been hypothesized that the recognition of these modified self-proteins by the immune system may promote autoantibody production. While apoptosis is specifically characterized by posttranslational modification of proteins, recent findings also show that nucleic acids are modified. This review summarizes the specific cleavages of some of these key nucleic acids, i.e. chromosomal DNA, ribosomal RNA and small structural RNAs (U1 snRNA, Y RNA), in apoptotic cells.

Inter- and intra-strain variation in the 5.8S ribosomal RNA and internal transcribed spacer sequences of Entamoeba histolytica and comparison with Entamoeba dispar, Entamoeba moshkovskii and Entamoeba invadens

The ribosomal RNA genes in Entamoeba histolytica are located on circular DNA molecules in about 200 copies per genome equivalent. Nucleotide sequence analysis of the 5.8S rRNA gene and the flanking internal transcribed spacers was carried out to determine the degree of sequence divergence in the multiple rRNA gene copies of a given strain; amongst three different E. histolytica strains (HM-1:IMSS, Rahman and HK-9); and amongst four species of Entamoeba (Entamoeba histolytica, Entamoeba dispar, Entamoeba moshkovskii and Entamoeba invadens). The results show that all rRNA gene copies of a given strain are identical. Few nucleotide positions varied between strains of a species but the differences were very pronounced amongst species. In general, the internal transcribed spacer 2 sequence was more variable and may be useful for strain- and species-identification. The 5.8S rRNA gene and the internal transcribed spacer 2 of E. invadens were unusually small in size.

Ribosomal DNA, species structure, and biogeography of the cactophilic yeast Clavispora opuntiae

The ribosomal DNA of the cactophilic yeast species Clavispora opuntiae was studied in order to clarify the global distribution of the yeast. Over 500 strains, including isolates from several new localities worldwide, were characterized by rDNA restriction mapping. An unusual restriction pattern previously encountered only in one strain, from Conception Island in the Bahamas, was found in several Brazilian isolates. Sequences of the D1/D2 and D7/D8 divergent domains of the large subunit (LSU) and of the intergenic spacers (IGS) confirmed that these strains represent a genetically distinct variety of Clavispora opuntiae. This divergence had previously been hypothesized on the basis of reduced genetic recombination in inter-varietal crosses and the presence of a polymorphic ApaI restriction site located in the LSU. The exact position of the ApaI site in the D8 divergent domain and the nature of the variation that it reveals were determined. The complete sequences of 12 intergenic spacers clarified the significance of the species-wide variation uncovered by restriction mapping. Most of the polymorphic sites occur in the IGS1 and IGS2 regions, on either side of the 5S gene, and the variation is largely due to differences in the numbers and the sequences of internal repeats. Two other polymorphic sites are located in the external transcribed spacer (ETS) region. The reliability of various sites as indicators of overall spacer sequence divergence differed from one case to another. Variety-specific probes were devised and used to screen 120 strains for the presence of recombinant rDNA spacers. Three strains gave ambiguous results, but these did not constitute evidence that inter-varietal recombination has taken place in nature. The hypothesis that the global movement of Clavispora opuntiae has been influenced by the worldwide biological control of prickly pear with Cactoblastis cactorum, a moth of Argentinian origin, has received additional support from the demonstration that Argentinian strains have rDNAs similar to those found where the moth has been introduced. A dramatic founder effect was identified in a yeast population collected in cacti (Maui, Hawaii) in a site where the moth had been recently introduced.

Domain III of Saccharomyces cerevisiae 25 S ribosomal RNA: its role in binding of ribosomal protein L25 and 60 S subunit formation

Domain III of Saccharomyces cerevisiae 25 S rRNA contains the recognition site for the primary rRNA-binding ribosomal protein L25, which belongs to the functionally conserved EL23/L25 family of ribosomal proteins. The EL23/L25 binding region is very complex, consisting of several irregular helices held together by long-distance secondary and tertiary interactions. Moreover, it contains the eukaryote-specific V9 (D7a) expansion segment. Functional characterisation of the structural elements of this site by a detailed in vitro and in vivo mutational analysis indicates the presence of two separate regions that are directly involved in L25 binding. In particular, mutation of either of two conserved nucleotides in the loop of helix 49 significantly reduces in vitro L25 binding, thus strongly supporting their role as attachment sites for the r-protein. Two other helices appear to be primarily required for the correct folding of the binding site. Mutations that abolish in vitro binding of L25 block accumulation of 25 S rRNA in vivo because they stall pre-rRNA processing at the level of its immediate precursor, the 27 S(B) pre-rRNA. Surprisingly, several mutations that do not significantly affect L25 binding in vitro cause the same lethal defect in 27 S(B) pre-rRNA processing. Deletion of the V9 expansion segment also leads to under-accumulation of mature 25 S rRNA and a twofold reduction in growth rate. We conclude that an intact domain III, including the V9 expansion segment, is essential for normal processing and assembly of 25 S rRNA.

The RNA interacting domain but not the protein interacting domain is highly conserved in ribosomal protein P0

Protein P0 interacts with proteins P1alpha, P1beta, P2alpha, and P2beta, and forms the Saccharomyces cerevisiae ribosomal stalk. The capacity of RPP0 genes from Aspergillus fumigatus, Dictyostelium discoideum, Rattus norvegicus, Homo sapiens, and Leishmania infantum to complement the absence of the homologous gene has been tested. In S. cerevisiae W303dGP0, a strain containing standard amounts of the four P1/P2 protein types, all heterologous genes were functional except the one from L. infantum, some of them inducing an osmosensitive phenotype at 37 degrees C. The polymerizing activity and the elongation factor-dependent functions but not the peptide bond formation capacity is affected in the heterologous P0 containing ribosomes. The heterologous P0 proteins bind to the yeast ribosomes but the composition of the ribosomal stalk is altered. Only proteins P1alpha and P2beta are found in ribosomes carrying the A. fumigatus, R. norvegicus, and H. sapiens proteins. When the heterologous genes are expressed in a conditional null-P0 mutant whose ribosomes are totally deprived of P1/P2 proteins, none of the heterologous P0 proteins complemented the conditional phenotype. In contrast, chimeric P0 proteins made of different amino-terminal fragments from mammalian origin and the complementary carboxyl-terminal fragments from yeast allow W303dGP0 and D67dGP0 growth at restrictive conditions. These results indicate that while the P0 protein RNA-binding domain is functionally conserved in eukaryotes, the regions involved in protein-protein interactions with either the other stalk proteins or the elongation factors have notably evolved.

Sequence differences in the internal transcribed spacers of DNA among four species of hookworm (Ancylostomatoidea: Ancylostoma)

The two ribosomal DNA internal transcribed spacers (1 and 2) of the hookworms Ancylostoma caninum, A. tubaeforme, A. ceylanicum and A. duodenale were sequenced. The sequence lengths were similar among the four species, except that A. ceylanicum had slightly longer (by 5-7 bp) internal transcribed spacer 1 and 2 sequences. The predicted secondary structure of the internal transcribed spacer 2 precursor rRNA was similar for all species, despite interspecific differences in primary sequence ranging from 0.9% to 13.2%. Interspecific differences in internal transcribed spacer 1 sequence ranged from 0.9% to 7.5%. A cladistic analysis of the sequence data, using the human hookworm Necator americanus as the outgroup, provided little resolution of the phylogenetic relationships, except that A. ceylanicum occurred on a branch external to the other three species. Nonetheless, internal transcribed spacers 1 and 2 may provide useful phylogenetic information at higher taxonomic levels within the superfamily Ancylostomatoidea.

A computational screen for methylation guide snoRNAs in yeast

Small nucleolar RNAs (snoRNAs) are required for ribose 2'-O-methylation of eukaryotic ribosomal RNA. Many of the genes for this snoRNA family have remained unidentified in Saccharomyces cerevisiae, despite the availability of a complete genome sequence. Probabilistic modeling methods akin to those used in speech recognition and computational linguistics were used to computationally screen the yeast genome and identify 22 methylation guide snoRNAs, snR50 to snR71. Gene disruptions and other experimental characterization confirmed their methylation guide function. In total, 51 of the 55 ribose methylated sites in yeast ribosomal RNA were assigned to 41 different guide snoRNAs.

rRNA-complementarity in the 5' untranslated region of mRNA specifying the Gtx homeodomain protein: evidence that base- pairing to 18S rRNA affects translational efficiency

Numerous eukaryotic mRNAs contain sequences complementary to segments of the 18S and 28S rRNAs. Previous studies raised the possibility that these complementarities might allow mRNA-rRNA interactions and affect rates of translation. In the present study, we investigated the mRNA encoding the mouse Gtx homeodomain protein. This mRNA contains within its 5' untranslated region (UTR) a segment that is complementary to two regions of the 18S rRNA, located at nucleotides 701-741 and 1104-1136. A Gtx RNA probe containing this complementarity could be photochemically cross-linked to ribosomal subunits through a linkage to 18S rRNA but not to 28S rRNA. Oligonucleotide-directed RNase H digestion of the rRNA and a reverse transcription analysis localized the cross-linked probe to the complementary segment of 18S rRNA at nucleotides 1104-1136 but not at nucleotides 701-741. To determine whether complementarity in the Gtx mRNA affected translation, a mutational analysis was performed with a Gtx-luciferase fusion construct and four related constructs with altered complementarity to the 18S rRNA. These constructs were examined for their ability to be translated in cell-free lysates prepared from P19 embryonal carcinoma and C6 glioma cell lines and after cellular transfection into these same cell lines. In both cell-free translation and transfection studies, the rate of translation decreased more than 9-fold as the degree of complementarity to nucleotides 1104-1136 of the 18S rRNA increased. We hypothesize that segments complementary to rRNA, such as those contained within the Gtx mRNA, form a category of cis-acting regulatory elements in mRNAs that affect translation by base pairing to rRNA within ribosomes.

Purification, cloning, and characterization of the 16 S RNA m2G1207 methyltransferase from Escherichia coli

The methyltransferase that forms m2G1207 in Escherichia coli small subunit rRNA has been purified, cloned, and characterized. The gene was identified from the N-terminal sequence of the purified enzyme as the open reading frame yjjT (SWISS-PROT accession number ). The gene, here renamed rsmC in view of its newly established function, codes for a 343-amino acid protein that has homologs in prokaryotes, Archaea, and possibly also in lower eukaryotes. The enzyme reacted well with 30 S subunits reconstituted from 16 S RNA transcripts and 30 S proteins but was almost inactive with the corresponding free RNA. By hybridization and protection of appropriate segments of 16 S RNA that had been extracted from 30 S subunits methylated by the enzyme, it was shown that of the three naturally occurring m2G residues, only m2G1207 was formed. Whereas close to unit stoichiometry of methylation could be achieved at 0.9 mM Mg2+, both 2 mM EDTA and 6 mM Mg2+ markedly reduced the level of methylation, suggesting that the optimal substrate may be a ribonucleoprotein particle less structured than a 30 S ribosome but more so than free RNA.

Common secondary structures for the second internal transcribed spacer pre-rRNA of two subfamilies of trichostrongylid nematodes

Sequences of the second internal transcribed spacer ribosomal DNA for the parasitic trichostrongylid nematodes Trichostrongylus probolurus, Trichostrongylus rugatus and Camelostrongylus mentulatus were compared with previously published sequences for five other species within the genus Trichostrongylus. The secondary structures of the second internal transcribed spacer pre-rRNA for these nematodes were predicted using an energy minimisation method. The results indicate that a common secondary structure of the second internal transcribed spacer of these nematodes is maintained despite distinct differences in primary sequence between species. Sequence differences among Trichostrongylus species ranged from 1.3 to 7.6%, but each species differed by 22-26% in sequence when compared with C. mentulatus which belongs to a different subfamily.

Reconstitution of peptide bond formation with Escherichia coli 23S ribosomal RNA domains

It was recently demonstrated that peptide bond formation can occur using an Escherichia coli naked 23S ribosomal RNA without any of the ribosomal proteins. Here, the six domains of the 23S ribosomal RNA were individually synthesized and shown to be capable, when complexed together, of stimulating the reaction. Omission and addition experiments indicated that the activity could be reconstituted solely by domain V at a concentration 10 times higher than that of the intact 23S ribosomal RNA, whereas domain VI could enhance the activity in trans. These findings suggest that fragments of an RNA molecule have the ability to associate into a functional whole.

Correlation of the expansion segments in mammalian rRNA with the fine structure of the 80 S ribosome; a cryoelectron microscopic reconstruction of the rabbit reticulocyte ribosome at 21 A resolution

Samples of 80 S ribosomes from rabbit reticulocytes were subjected to electron cryomicroscopy combined with angular reconstitution. A three-dimensional reconstruction at 21 A resolution was obtained, which was compared with the corresponding (previously published) reconstruction of Escherichia coli 70 S ribosomes carrying tRNAs at the A and P sites. In the region of the intersubunit cavity, the principal features observed in the 70 S ribosome (such as the L1 protuberance, the central protuberance and A site finger in the large subunit) could all be clearly identified in the 80 S particle. On the other hand, significant additional features were observed in the 80 S ribosomes on the solvent sides and lower regions of both subunits. In the case of the small (40 S) subunit, the most prominent additions are two extensions at the base of the particle. By comparing the secondary structure of the rabbit 18 S rRNA with our model for the three-dimensional arrangement of E. coli 16 S rRNA, these two extensions could be correlated with the rabbit expansion segments (each totalling ca 170 bases) in the regions of helix 21, and of helices 8, 9 and 44, respectively. A similar comparison of the secondary structures of mammalian 28 S rRNA and E. coli 23 S rRNA, combined with preliminary modelling studies on the 23 S rRNA within the 50 S subunit, enabled the additional features in the 60 S subunit to be sub-divided into five groups. The first (corresponding to a total of ca 335 extra bases in helices 45, 98 and 101) is located on the solvent side of the 60 S subunit, close to the L7/L12 area. The second (820 bases in helices 25 and 38) is centrally placed on the solvent side of the subunit, whereas the third group (totaling 225 bases in helices 18/19, 27/29, 52 and 54) lies towards the L1 side of the subunit. The fourth feature (80 bases in helices 78 and 79) lies within or close to the L1 protuberance itself, and the fifth (560 bases in helix 63) is located underneath the L1 protuberance on the interface side of the 60 S subunit.

The rRNA-processing function of the yeast U14 small nucleolar RNA can be rescued by a conserved RNA helicase-like protein

The phylogenetically conserved U14 small nucleolar RNA is required for processing of rRNA, and this function involves base pairing with conserved complementary sequences in 18S RNA. With a view to identifying other important U14 interactions, a stem-loop domain required for activity of Saccharomyces cerevisiae U14 RNAs (the Y domain) was first subjected to detailed mutational analysis. The mapping results showed that most nucleotides of the Y domain can be replaced without affecting function, except for loop nucleotides conserved among five different yeast species. Defective variants were then used to identify both intragenic and extragenic suppressor mutations. All of the intragenic mutations mapped within six nucleotides of the primary mutation, suggesting that suppression involves a change in conformation and that the loop element is involved in an essential intermolecular interaction rather than intramolecular base pairing. A high-copy extragenic suppressor gene, designated DBP4 (DEAD box protein 4), encodes an essential, putative RNA helicase of the DEAD-DEXH box family. Suppression by DBP4 (initially CA4 [T.-H. Chang, J. Arenas, and J. Abelson, Proc. Natl. Acad. Sci. USA 87:1571-1575, 1990]) restores the level of 18S rRNA and is specific for the Y domain but is not allele specific. DBP4 is predicted to function either in assembly of the U14 small nucleolar RNP or, more likely, in its interaction with other components of the rRNA processing apparatus. Mediating the interaction of U14 with precursor 18S RNA is an especially attractive possibility.

Eukaryotic ribosomal RNA: the recent excitement in the nucleotide modification problem

Eukaryotic ribosomal RNA (rRNA) contains numerous modified nucleotides: about 115 methyl groups and some 95 pseudouridines in vertebrates; about 65 methyl groups and some 45 pseudouridines in Saccharomyces cerevisiae. All but about ten of the methyl groups are ribose methylations. The remaining ten are on heterocyclic bases. The ribose methylations occur very rapidly upon the primary rRNA transcript in the nucleolus, probably on nascent chains, and they appear to play an important role in ribosome maturation, at least in vertebrates. All of the methyl groups occur in the conserved core of rRNA. However, there is no consensus feature of sequence or secondary structure for the methylation sites; thus the nature of the signal(s) for site-specific methylations had until recently remained a mystery. The situation changed dramatically with the discovery that many of the ribose methylation sites are in regions that are precisely complementary to small nucleolar RNA (snoRNA) species. Experimental evidence indicates that structural motifs within the snoRNA species do indeed pinpoint the precise nucleotides to be methylated by the putative 2'-O-methyl transferase(s). Regarding base methylations, the gene DIM1, responsible for modification of the conserved dimethyladenosines near the 3' end of 18S rRNA, has been shown to be essential for viability in S. cerevisiae and is suggested to play a role in the nucleocytoplasmic transport of the small ribosomal subunit. Recently nearly all of the pseudouridines have also been mapped in the rRNA of several eukaryotic species. As is the case for ribose methylations, most pseudouridine modifications occur rapidly upon precursor rRNA, within core sequences, and in a variety of local primary and secondary structure environments. In contrast to ribose methylation, no potentially unifying process has yet been identified for the enzymic recognition of the many pseudouridine modification sites. However, the new data afford the basis for a search for any potential involvement of snoRNAs in the recognition process.

Variable region V1 of Saccharomyces cerevisiae 18S rRNA participates in biogenesis and function of the small ribosomal subunit

The role of helix 6, which forms the major portion of the most 5'-located expansion segment of Saccharomyces cerevisiae 18S rRNA, was studied by in vivo mutational analysis. Mutations that increased the size of the helical part and/or the loop, even to a relatively small extent, abolished 18S rRNA formation almost completely. Concomitantly, 35S pre-rRNA and an abnormal 23S precursor species accumulated. rDNA units containing these mutations did not support cell growth. A deletion removing helix 6 almost completely, on the other hand, had a much less severe effect on the formation of 18S rRNA, and cells expressing only the mutant rRNA remained able to grow, albeit at a much reduced rate. Disruption of the apical A.U base pair by a single point mutation did not cause a noticeable reduction in the level of 18S rRNA but did result in a twofold lower growth rate of the cells. This effect could not be reversed by introduction of a second point mutation that restores base pairing. We conclude that both the primary and the secondary structure of helix 6 play an important role in the formation and the biological function of the 40S subunit.

Identification of additional rRNA fragments encoded by the Plasmodium falciparum 6 kb element

Sequences similar to mitochondrial large and small subunit rRNAs are found as small scattered fragments on a tandemly reiterated 6 kb element in the human malaria parasite Plasmodium falciparum. The rDNA sequences previously identified include strongly conserved portions of rRNA, suggesting that fragmented rRNAs derived from them are able to associate into functional ribosomes. However, sequences corresponding to other expected rRNA regions were not found. We here report that 10 of the 13 previously described rDNA regions have abundant small transcripts. An additional 10 transcripts were found from regions not previously known to contain genes. Five of the latter have been identified as rRNA fragments, including those corresponding to the 5'end and 790 loop sequences of small subunit rRNA and the sarcin/ ricin loop of large subunit rRNA. Demonstration that most of the previously described rDNA regions have abundant transcripts and the identification of new transcripts with other portions of conventional rRNAs provide support for the hypothesis that these small transcripts comprise functional rRNAs.

Processing of mammalian rRNA precursors at the 3' end of 18S rRNA. Identification of cis-acting signals suggests the involvement of U13 small nucleolar RNA

Molecular mechanisms involved in the nucleolytic cleavage at the 18S rRNA/internal transcribed spacer 1 (ITS 1) junction, a late step of small-subunit pre-rRNA processing in vertebrates, remain largely unknown, mostly due to the lack of faithful in vitro assays. To identify the minimal cis-acting signals required for this reaction, we studied the processing of truncated human rRNA gene transcripts transiently expressed upon transfection of rRNA minigenes into cultured mouse cells. We observed that processing at this site was faithfully reproduced with transcripts containing only 60 nucleotides of 18S rRNA and the adjacent 103 nucleotides of ITS 1, but was abolished or severely altered by further shortening of either sequence. Remarkably, this minimal transcript contains, within its 18S rRNA part, long sequences complementary to both U20 and U13 small nucleolar RNAs (snoRNAs). The cis-acting elements essential for the reaction were studied further by site-directed mutagenesis. The U20 snoRNA complementary region in 18S rRNA was not required for faithful processing at the 18S rRNA/ITS 1 junction. Also, processing at this site was not appreciably altered by random substitution of proximal ITS 1 sequences (including the 5' terminal nucleotide) or of the terminal nucleotide of mature 18S rRNA. Substitutions in the four-nucleotide loop of the 18S rRNA 3'-terminal stem-loop, including the two adenosine residues substrates of dimethylation, did not alter appreciably the formation of the 18S rRNA 3' end, showing that the (methyl)2A1850.(methyl)2A1851 doublet was not required for processing at this site. Two highly conserved 18S rRNA elements acted as major cis-acting signals for processing at the 3' end, the CAUU sequence immediately preceding the 3'-terminal nucleotide and the 3' strand of the 3'-terminal 18S rRNA helix, complementary to U13 snoRNA. Compensatory mutations, restoring the potential for helix formation, but not U13 snoRNA complementarity, did not restitute the cleavage at the 3' end of 18S rRNA. This suggests that U13 snoRNA may be a trans-acting factor in the nucleolytic cleavage at the 3' end of 18S rRNA.

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

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

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.