Structural implications of novel diversity in eucaryal RNase P RNA

RNase P: Beyond Precursor tRNA Processing

Ribonuclease P (RNase P) was first described in the 1970's as an endoribonuclease acting in the maturation of precursor transfer RNAs (tRNAs). More recent studies, however, have uncovered non-canonical roles for RNase P and its components. Here, we review the recent progress of its involvement in chromatin assembly, DNA damage response, and maintenance of genome stability with implications in tumorigenesis. The possibility of RNase P as a therapeutic target in cancer is also discussed.

Piece by piece: Building a ribozyme

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

rPredictorDB: a predictive database of individual secondary structures of RNAs and their formatted plots

Secondary data structure of RNA molecules provides insights into the identity and function of RNAs. With RNAs readily sequenced, the question of their structural characterization is increasingly important. However, RNA structure is difficult to acquire. Its experimental identification is extremely technically demanding, while computational prediction is not accurate enough, especially for large structures of long sequences. We address this difficult situation with rPredictorDB, a predictive database of RNA secondary structures that aims to form a middle ground between experimentally identified structures in PDB and predicted consensus secondary structures in Rfam. The database contains individual secondary structures predicted using a tool for template-based prediction of RNA secondary structure for the homologs of the RNA families with at least one homolog with experimentally solved structure. Experimentally identified structures are used as the structural templates and thus the prediction has higher reliability than de novo predictions in Rfam. The sequences are downloaded from public resources. So far rPredictorDB covers 7365 RNAs with their secondary structures. Plots of the secondary structures use the Traveler package for readable display of RNAs with long sequences and complex structures, such as ribosomal RNAs. The RNAs in the output of rPredictorDB are extensively annotated and can be viewed, browsed, searched and downloaded according to taxonomic, sequence and structure data. Additionally, structure of user-provided sequences can be predicted using the templates stored in rPredictorDB.

The Diversity of Ribonuclease P: Protein and RNA Catalysts with Analogous Biological Functions

Ribonuclease P (RNase P) is an essential endonuclease responsible for catalyzing 5' end maturation in precursor transfer RNAs. Since its discovery in the 1970s, RNase P enzymes have been identified and studied throughout the three domains of life. Interestingly, RNase P is either RNA-based, with a catalytic RNA subunit, or a protein-only (PRORP) enzyme with differential evolutionary distribution. The available structural data, including the active site data, provides insight into catalysis and substrate recognition. The hydrolytic and kinetic mechanisms of the two forms of RNase P enzymes are similar, yet features unique to the RNA-based and PRORP enzymes are consistent with different evolutionary origins. The various RNase P enzymes, in addition to their primary role in tRNA 5' maturation, catalyze cleavage of a variety of alternative substrates, indicating a diversification of RNase P function in vivo. The review concludes with a discussion of recent advances and interesting research directions in the field.

Multiple Layers of Stress-Induced Regulation in tRNA Biology

tRNAs are the fundamental components of the translation machinery as they deliver amino acids to the ribosomes during protein synthesis. Beyond their essential function in translation, tRNAs also function in regulating gene expression, modulating apoptosis and several other biological processes. There are multiple layers of regulatory mechanisms in each step of tRNA biogenesis. For example, tRNA 3′ trailer processing is altered upon nutrient stress; tRNA modification is reprogrammed under various stresses; nuclear accumulation of tRNAs occurs upon nutrient deprivation; tRNA halves accumulate upon oxidative stress. Here we address how environmental stresses can affect nearly every step of tRNA biology and we describe the possible regulatory mechanisms that influence the function or expression of tRNAs under stress conditions.

Distribution of Ribonucleoprotein and Protein-Only RNase P in Eukarya

RNase P is the endonuclease that removes 5' leader sequences from tRNA precursors. In Eukarya, separate RNase P activities exist in the nucleus and mitochondria/plastids. Although all RNase P enzymes catalyze the same reaction, the different architectures found in Eukarya range from ribonucleoprotein (RNP) enzymes with a catalytic RNA and up to 10 protein subunits to single-subunit protein-only RNase P (PRORP) enzymes. Here, analysis of the phylogenetic distribution of RNP and PRORP enzymes in Eukarya revealed 1) a wealth of novel P RNAs in previously unexplored phylogenetic branches and 2) that PRORP enzymes are more widespread than previously appreciated, found in four of the five eukaryal supergroups, in the nuclei and/or organelles. Intriguingly, the occurrence of RNP RNase P and PRORP seems mutually exclusive in genetic compartments of modern Eukarya. Our comparative analysis provides a global picture of the evolution and diversification of RNase P throughout Eukarya.

RNase P RNA from the recently evolved plastid of Paulinella and from algae

The RNase P RNA catalytic subunit (RPR) encoded in some plastids has been found to be functionally defective. The amoeba Paulinella chromatophora contains an organelle (chromatophore) that is derived from the recent endosymbiotic acquisition of a cyanobacterium, and therefore represents a model of the early steps in the acquisition of plastids. In contrast with plastid RPRs the chromatophore RPR retains functionality similar to the cyanobacterial enzyme. The chromatophore RPR sequence deviates from consensus at some positions but those changes allow optimal activity compared with mutated chromatophore RPR with the consensus sequence. We have analyzed additional RPR sequences identifiable in plastids and have found that it is present in all red algae and in several prasinophyte green algae. We have assayed in vitro a subset of the plastid RPRs not previously analyzed and confirm that these organelle RPRs lack RNase P activity in vitro.

Vicinal: a method for the determination of ncRNA ends using chimeric reads from RNA-seq experiments

Non-coding (nc)RNAs are important structural and regulatory molecules. Accurate determination of the primary sequence and secondary structure of ncRNAs is important for understanding their functions. During cDNA synthesis, RNA 3' end stem-loops can self-prime reverse transcription, creating RNA-cDNA chimeras. We found that chimeric RNA-cDNA fragments can also be detected at 5' end stem-loops, although at much lower frequency. Using the Gubler-Hoffman method, both types of chimeric fragments can be converted to cDNA during library construction, and they are readily detectable in high-throughput RNA sequencing (RNA-seq) experiments. Here, we show that these chimeric reads contain valuable information about the boundaries of ncRNAs. We developed a bioinformatic method, called Vicinal, to precisely map the ends of numerous fruitfly, mouse and human ncRNAs. Using this method, we analyzed chimeric reads from over 100 RNA-seq datasets, the results of which we make available for users to find RNAs of interest. In summary, we show that Vicinal is a useful tool for determination of the precise boundaries of uncharacterized ncRNAs, facilitating further structure/function studies.

Structural studies of RNase P

Ribonuclease P (RNase P) is one of the first ribozymes discovered and it is found in all phylogenetic groups. It is responsible for processing the 5' end of pre-tRNAs as well as other RNA molecules. RNase P is formed by an RNA molecule responsible for catalysis and one or more proteins. Structural studies of the proteins from different organisms, the bacterial RNA component, and a bacterial RNase P holoenzyme/tRNA complex provide insights into the mechanism of this universal ribozyme. Together with the existing wealth of biochemical information, these studies provide atomic-level information on the mechanism of RNase P and continue to expand our understanding of the structure and architecture of large RNA molecules and ribonucleoprotein complexes, the nature of catalysis by ribozymes, the structural basis of recognition of RNA by RNA molecules, and the evolution of enzymes from the prebiotic, RNA-based world to the modern world.

RNase P branches out from RNP to protein: organelle-triggered diversification?

RNase P is the enzyme that removes 5' leader sequences from precursor tRNAs. Remarkably, in most organisms, RNase P is a ribonucleoprotein particle where the RNA component is responsible for catalysis. In this issue of Genes & Development, Gutmann and colleagues (pp. 1022-1027) report the first organism, Arabidopsis thaliana, to employ protein-only RNase P in both its nucleus and organelles. An intriguing possibility is that replacement of RNase P ribonucleoprotein particles (RNPs) by proteins may have been triggered by the acquisition of organelles.

Prostaglandin A2 influences gene expression in an established insect cell line (BCIRL-HzAM1) cells

Prostaglandins (PGs) and other eicosanoids are oxygenated metabolites of arachidonic acid and two other C(20) polyunsaturated fatty acids. While most well studied in mammals, PGs exert important actions in insects and virtually all other invertebrates. We have been researching the mechanisms of PG actions in established insect cell lines and reported earlier that two PGs, PGA(1) and PGE(1), influence gene and protein expression in HzAM1 cells. Here we report on further experiments with three 2-series PGs, PGA(2), PGE(2) and PGF(2α). In separate experiments we treated cells with each of the three PGs for 12 and 24h and then analyzed cell lysates by 2-D electrophoresis. Analysis of the gels by Delta2D software showed that PGA(2) influenced expression of 60 proteins while PGE(2) and PGF(2α) treatments led to expression changes for only a few proteins. All spots representing changes in protein expression were processed for analysis by MALDI TOF/TOF mass spectrometry. Bioinformatic analysis of the resulting sequences yielded in silico identifications of all proteins. The apparent changes in some proteins were confirmed by quantitative PCR, which demonstrated that changes in protein expression were parallel to changes in mRNA expression. We assorted the proteins into functional categories, including 1/cell structure and function; 2/cell protection and immunity; 3/energetics and metabolism; 4/nucleotide processing; 5/protein action and processing and 6/signal transduction. These findings substantially extend our idea that one mechanism of PG actions in insect cells is the modulation of gene and protein expression.

RNase MRP RNA and RNase P activity in plants are associated with a Pop1p containing complex

RNase P processes the 5'-end of tRNAs. An essential catalytic RNA has been demonstrated in Bacteria, Archaea and the nuclei of most eukaryotes; an organism-specific number of proteins complement the holoenzyme. Nuclear RNase P from yeast and humans is well understood and contains an RNA, similar to the sister enzyme RNase MRP. In contrast, no protein subunits have yet been identified in the plant enzymes, and the presence of a nucleic acid in RNase P is still enigmatic. We have thus set out to identify and characterize the subunits of these enzymes in two plant model systems. Expression of the two known Arabidopsis MRP RNA genes in vivo was verified. The first wheat MRP RNA sequences are presented, leading to improved structure models for plant MRP RNAs. A novel mRNA encoding the central RNase P/MRP protein Pop1p was identified in Arabidopsis, suggesting the expression of distinct protein variants from this gene in vivo. Pop1p-specific antibodies precipitate RNase P activity and MRP RNAs from wheat extracts. Our results provide evidence that in plants, Pop1p is associated with MRP RNAs and with the catalytic subunit of RNase P, either separately or in a single large complex.

A functional RNase P protein subunit of bacterial origin in some eukaryotes

RNase P catalyzes 5'-maturation of tRNAs. While bacterial RNase P comprises an RNA catalyst and a protein cofactor, the eukaryotic (nuclear) variant contains an RNA and up to ten proteins, all unrelated to the bacterial protein. Unexpectedly, a nuclear-encoded bacterial RNase P protein (RPP) homolog is found in several prasinophyte algae including Ostreococcus tauri. We demonstrate that recombinant O. tauri RPP can functionally reconstitute with bacterial RNase P RNAs (RPRs) but not with O. tauri organellar RPRs, despite the latter's presumed bacterial origins. We also show that O. tauri PRORP, a homolog of Arabidopsis PRORP-1, displays tRNA 5'-processing activity in vitro. We discuss the implications of the striking diversity of RNase P in O. tauri, the smallest known free-living eukaryote.

Discovery of a minimal form of RNase P in Pyrobaculum

RNase P RNA is an ancient, nearly universal feature of life. As part of the ribonucleoprotein RNase P complex, the RNA component catalyzes essential removal of 5' leaders in pre-tRNAs. In 2004, Li and Altman computationally identified the RNase P RNA gene in all but three sequenced microbes: Nanoarchaeum equitans, Pyrobaculum aerophilum, and Aquifex aeolicus (all hyperthermophiles) [Li Y, Altman S (2004) RNA 10:1533-1540]. A recent study concluded that N. equitans does not have or require RNase P activity because it lacks 5' tRNA leaders. The "missing" RNase P RNAs in the other two species is perplexing given evidence or predictions that tRNAs are trimmed in both, prompting speculation that they may have developed novel alternatives to 5' pre-tRNA processing. Using comparative genomics and improved computational methods, we have now identified a radically minimized form of the RNase P RNA in five Pyrobaculum species and the related crenarchaea Caldivirga maquilingensis and Vulcanisaeta distributa, all retaining a conventional catalytic domain, but lacking a recognizable specificity domain. We confirmed 5' tRNA processing activity by high-throughput RNA sequencing and in vitro biochemical assays. The Pyrobaculum and Caldivirga RNase P RNAs are the smallest naturally occurring form yet discovered to function as trans-acting precursor tRNA-processing ribozymes. Loss of the specificity domain in these RNAs suggests altered substrate specificity and could be a useful model for finding other potential roles of RNase P. This study illustrates an effective combination of next-generation RNA sequencing, computational genomics, and biochemistry to identify a divergent, formerly undetectable variant of an essential noncoding RNA gene.

tRNA biology charges to the front

tRNA biology has come of age, revealing an unprecedented level of understanding and many unexpected discoveries along the way. This review highlights new findings on the diverse pathways of tRNA maturation, and on the formation and function of a number of modifications. Topics of special focus include the regulation of tRNA biosynthesis, quality control tRNA turnover mechanisms, widespread tRNA cleavage pathways activated in response to stress and other growth conditions, emerging evidence of signaling pathways involving tRNA and cleavage fragments, and the sophisticated intracellular tRNA trafficking that occurs during and after biosynthesis.

Archaeal/eukaryal RNase P: subunits, functions and RNA diversification

RNase P, a catalytic ribonucleoprotein (RNP), is best known for its role in precursor tRNA processing. Recent discoveries have revealed that eukaryal RNase P is also required for transcription and processing of select non-coding RNAs, thus enmeshing RNase P in an intricate network of machineries required for gene expression. Moreover, the RNase P RNA seems to have been subject to gene duplication, selection and divergence to generate two new catalytic RNPs, RNase MRP and MRP-TERT, which perform novel functions encompassing cell cycle control and stem cell biology. We present new evidence and perspectives on the functional diversification of the RNase P RNA to highlight it as a paradigm for the evolutionary plasticity that underlies the extant broad repertoire of catalytic and unexpected regulatory roles played by RNA-driven RNPs.

Dissecting functional cooperation among protein subunits in archaeal RNase P, a catalytic ribonucleoprotein complex

RNase P catalyzes the Mg²⁺-dependent 5′-maturation of precursor tRNAs. Biochemical studies on the bacterial holoenzyme, composed of one catalytic RNase P RNA (RPR) and one RNase P protein (RPP), have helped understand the pleiotropic roles (including substrate/Mg²⁺ binding) by which a protein could facilitate RNA catalysis. As a model for uncovering the functional coordination among multiple proteins that aid an RNA catalyst, we use archaeal RNase P, which comprises one catalytic RPR and at least four RPPs. Exploiting our previous finding that these archaeal RPPs function as two binary RPP complexes (POP5•RPP30 and RPP21•RPP29), we prepared recombinant RPP pairs from three archaea and established interchangeability of subunits through homologous/heterologous assemblies. Our finding that archaeal POP5•RPP30 reconstituted with bacterial and organellar RPRs suggests functional overlap of this binary complex with the bacterial RPP and highlights their shared recognition of a phylogenetically-conserved RPR catalytic core, whose minimal attributes we further defined through deletion mutagenesis. Moreover, single-turnover kinetic studies revealed that while POP5•RPP30 is solely responsible for enhancing the RPR’s rate of precursor tRNA cleavage (by 60-fold), RPP21•RPP29 contributes to increased substrate affinity (by 16-fold). Collectively, these studies provide new perspectives on the functioning and evolution of an ancient, catalytic ribonucleoprotein.

Bcheck: a wrapper tool for detecting RNase P RNA genes

BACKGROUND

Effective bioinformatics solutions are needed to tackle challenges posed by industrial-scale genome annotation. We present Bcheck, a wrapper tool which predicts RNase P RNA genes by combining the speed of pattern matching and sensitivity of covariance models. The core of Bcheck is a library of subfamily specific descriptor models and covariance models.

RESULTS

Scanning all microbial genomes in GenBank identifies RNase P RNA genes in 98% of 1024 microbial chromosomal sequences within just 4 hours on single CPU. Comparing to existing annotations found in 387 of the GenBank files, Bcheck predictions have more intact structure and are automatically classified by subfamily membership. For eukaryotic chromosomes Bcheck could identify the known RNase P RNA genes in 84 out of 85 metazoan genomes and 19 out of 21 fungi genomes. Bcheck predicted 37 novel eukaryotic RNase P RNA genes, 32 of which are from fungi. Gene duplication events are observed in at least 20 metazoan organisms. Scanning of meta-genomic data from the Global Ocean Sampling Expedition, comprising over 10 million sample sequences (18 Gigabases), predicted 2909 unique genes, 98% of which fall into ancestral bacteria A type of RNase P RNA and 66% of which have no close homolog to known prokaryotic RNase P RNA.

CONCLUSIONS

The combination of efficient filtering by means of a descriptor-based search and subsequent construction of a high-quality gene model by means of a covariance model provides an efficient method for the detection of RNase P RNA genes in large-scale sequencing data. Bcheck is implemented as webserver and can also be downloaded for local use from http://rna.tbi.univie.ac.at/bcheck.

The ancient history of the structure of ribonuclease P and the early origins of Archaea

BACKGROUND

Ribonuclease P is an ancient endonuclease that cleaves precursor tRNA and generally consists of a catalytic RNA subunit (RPR) and one or more proteins (RPPs). It represents an important macromolecular complex and model system that is universally distributed in life. Its putative origins have inspired fundamental hypotheses, including the proposal of an ancient RNA world.

RESULTS

To study the evolution of this complex, we constructed rooted phylogenetic trees of RPR molecules and substructures and estimated RPP age using a cladistic method that embeds structure directly into phylogenetic analysis. The general approach was used previously to study the evolution of tRNA, SINE RNA and 5S rRNA, the origins of metabolism, and the evolution and complexity of the protein world, and revealed here remarkable evolutionary patterns. Trees of molecules uncovered the tripartite nature of life and the early origin of archaeal RPRs. Trees of substructures showed molecules originated in stem P12 and were accessorized with a catalytic P1-P4 core structure before the first substructure was lost in Archaea. This core currently interacts with RPPs and ancient segments of the tRNA molecule. Finally, a census of protein domain structure in hundreds of genomes established RPPs appeared after the rise of metabolic enzymes at the onset of the protein world.

CONCLUSIONS

The study provides a detailed account of the history and early diversification of a fundamental ribonucleoprotein and offers further evidence in support of the existence of a tripartite organismal world that originated by the segregation of archaeal lineages from an ancient community of primordial organisms.

De novo prediction of structured RNAs from genomic sequences

Growing recognition of the numerous, diverse and important roles played by non-coding RNA in all organisms motivates better elucidation of these cellular components. Comparative genomics is a powerful tool for this task and is arguably preferable to any high-throughput experimental technology currently available, because evolutionary conservation highlights functionally important regions. Conserved secondary structure, rather than primary sequence, is the hallmark of many functionally important RNAs, because compensatory substitutions in base-paired regions preserve structure. Unfortunately, such substitutions also obscure sequence identity and confound alignment algorithms, which complicates analysis greatly. This paper surveys recent computational advances in this difficult arena, which have enabled genome-scale prediction of cross-species conserved RNA elements. These predictions suggest that a wealth of these elements indeed exist.

High throughput genome-wide survey of small RNAs from the parasitic protists Giardia intestinalis and Trichomonas vaginalis

RNA interference (RNAi) is a set of mechanisms which regulate gene expression in eukaryotes. Key elements of RNAi are small sense and antisense RNAs from 19 to 26 nt generated from double-stranded RNAs. MicroRNAs (miRNAs) are a major type of RNAi-associated small RNAs and are found in most eukaryotes studied to date. To investigate whether small RNAs associated with RNAi appear to be present in all eukaryotic lineages, and therefore present in the ancestral eukaryote, we studied two deep-branching protozoan parasites, Giardia intestinalis and Trichomonas vaginalis. Little is known about endogenous small RNAs involved in RNAi of these organisms. Using Illumina Solexa sequencing and genome-wide analysis of small RNAs from these distantly related deep-branching eukaryotes, we identified 10 strong miRNA candidates from Giardia and 11 from Trichomonas. We also found evidence of Giardia short-interfering RNAs potentially involved in the expression of variant-specific surface proteins. In addition, eight new small nucleolar RNAs from Trichomonas are identified. Our results indicate that miRNAs are likely to be general in ancestral eukaryotes and therefore are likely to be a universal feature of eukaryotes.

The making of tRNAs and more - RNase P and tRNase Z

Transfer-RNA (tRNA) molecules are essential players in protein biosynthesis. They are transcribed as precursors, which have to be extensively processed at both ends to become functional adaptors in protein synthesis. Two endonucleases that directly interact with the tRNA moiety, RNase P and tRNase Z, remove extraneous nucleotides on the molecule's 5'- and 3'-side, respectively. The ribonucleoprotein enzyme RNase P was identified almost 40 years ago and is considered a vestige from the "RNA world". Here, we present the state of affairs on prokaryotic RNase P, with a focus on recent findings on its role in RNA metabolism. tRNase Z was only identified 6 years ago, and we do not yet have a comprehensive understanding of its function. The current knowledge on prokaryotic tRNase Z in tRNA 3'-processing is reviewed here. A second, tRNase Z-independent pathway of tRNA 3'-end maturation involving 3'-exonucleases will also be discussed.

Detection of Balamuthia mandrillaris DNA by real-time PCR targeting the RNase P gene

Background

The free-living amoeba Balamuthia mandrillaris may cause fatal encephalitis both in immunocompromised and in – apparently – immunocompetent humans and other mammalian species. Rapid, specific, sensitive, and reliable detection requiring little pathogen-specific expertise is an absolute prerequisite for a successful therapy and a welcome tool for both experimental and epidemiological research.

Results

A real-time polymerase chain reaction assay using TaqMan^® probes (real-time PCR) was established specifically targeting the RNase P gene of B. mandrillaris amoebae. The assay detected at least 2 (down to 0.5) genomes of B. mandrillaris grown in axenic culture. It did not react with DNA from closely related Acanthamoeba (3 species), nor with DNA from Toxoplasma gondii, Leishmania major, Pneumocystis murina, Mycobacterium bovis (BCG), human brain, various mouse organs, or from human and murine cell lines. The assay efficiently detected B. mandrillaris DNA in spiked cell cultures, spiked murine organ homogenates, B. mandrillaris-infected mice, and CNS tissue-DNA preparations from 2 patients with proven cerebral balamuthiasis. This novel primer set was successfully combined with a published set that targets the B. mandrillaris 18S rRNA gene in a duplex real-time PCR assay to ensure maximum specificity and as a precaution against false negative results.

Conclusion

A real-time PCR assay for B. mandrillaris amoebae is presented, that is highly specific, sensitive, and reliable and thus suited both for diagnosis and for research.

Computational identification of four spliceosomal snRNAs from the deep-branching eukaryote Giardia intestinalis

RNAs processing other RNAs is very general in eukaryotes, but is not clear to what extent it is ancestral to eukaryotes. Here we focus on pre-mRNA splicing, one of the most important RNA-processing mechanisms in eukaryotes. In most eukaryotes splicing is predominantly catalysed by the major spliceosome complex, which consists of five uridine-rich small nuclear RNAs (U-snRNAs) and over 200 proteins in humans. Three major spliceosomal introns have been found experimentally in Giardia; one Giardia U-snRNA (U5) and a number of spliceosomal proteins have also been identified. However, because of the low sequence similarity between the Giardia ncRNAs and those of other eukaryotes, the other U-snRNAs of Giardia had not been found. Using two computational methods, candidates for Giardia U1, U2, U4 and U6 snRNAs were identified in this study and shown by RT-PCR to be expressed. We found that identifying a U2 candidate helped identify U6 and U4 based on interactions between them. Secondary structural modelling of the Giardia U-snRNA candidates revealed typical features of eukaryotic U-snRNAs. We demonstrate a successful approach to combine computational and experimental methods to identify expected ncRNAs in a highly divergent protist genome. Our findings reinforce the conclusion that spliceosomal small-nuclear RNAs existed in the last common ancestor of eukaryotes.

De novo search for non-coding RNA genes in the AT-rich genome of Dictyostelium discoideum: performance of Markov-dependent genome feature scoring

Genome data are increasingly important in the computational identification of novel regulatory non-coding RNAs (ncRNAs). However, most ncRNA gene-finders are either specialized to well-characterized ncRNA gene families or require comparisons of closely related genomes. We developed a method for de novo screening for ncRNA genes with a nucleotide composition that stands out against the background genome based on a partial sum process. We compared the performance when assuming independent and first-order Markov-dependent nucleotides, respectively, and used Karlin-Altschul and Karlin-Dembo statistics to evaluate the significance of hits. We hypothesized that a first-order Markov-dependent process might have better power to detect ncRNA genes since nearest-neighbor models have been shown to be successful in predicting RNA structures. A model based on a first-order partial sum process (analyzing overlapping dinucleotides) had better sensitivity and specificity than a zeroth-order model when applied to the AT-rich genome of the amoeba Dictyostelium discoideum. In this genome, we detected 94% of previously known ncRNA genes (at this sensitivity, the false positive rate was estimated to be 25% in a simulated background). The predictions were further refined by clustering candidate genes according to sequence similarity and/or searching for an ncRNA-associated upstream element. We experimentally verified six out of 10 tested ncRNA gene predictions. We conclude that higher-order models, in combination with other information, are useful for identification of novel ncRNA gene families in single-genome analysis of D. discoideum. Our generalizable approach extends the range of genomic data that can be searched for novel ncRNA genes using well-grounded statistical methods.

Activation of bacterial ribonuclease P by macrolides

The effect of macrolide antibiotic spiramycin on RNase P holoenzyme and M1 RNA from Escherichia coli was investigated. Ribonuclease P (RNase P) is a ribozyme that is responsible for the maturation of 5' termini of tRNA molecules. Spiramycin revealed a dose-dependent activation on pre-tRNA cleavage by E. coli RNase P holoenzyme and M1 RNA. The K s and V max, as well as the K s(app) and V max(app) values of RNase P holoenzyme and M1 RNA in the presence or absence of spiramycin, were calculated from primary and secondary kinetic plots. It was found that the activity status of RNase P holoenzyme and M1 RNA is improved by the presence of spiramycin 18- and 12-fold, respectively. Primer extension analysis revealed that spiramycin induces a conformational change of the P10/11 structural element of M1 RNA, which is involved in substrate recognition.

RNase P of the Cyanophora paradoxa cyanelle: A plastid ribozyme

Ribonuclease P (RNase P) is a ribonucleoprotein enzyme that generates the mature 5' ends of tRNAs. Ubiquitous across all three kingdoms of life, the composition and functional contributions of the RNA and protein components of RNase P differ between the kingdoms. RNA-alone catalytic activity has been reported throughout bacteria, but only for some archaea, and only as trace activity for eukarya. Available information for RNase P from photosynthetic organelles points to large differences to bacterial as well as to eukaryotic RNase P: for spinach chloroplasts, protein-alone activity has been discussed; for RNase P from the cyanelle of the glaucophyte Cyanophora paradoxa, a type of organelle sharing properties of both cyanobacteria and chloroplasts, the proportion of protein was found to be around 80% rather than the usual 10% in bacteria. Furthermore, the latter RNase P was previously found catalytically inactive in the absence of protein under a variety of conditions; however, the RNA could be activated by a cyanobacterial protein, but not by the bacterial RNase P protein from Escherichia coli. Here we demonstrate that, under very high enzyme concentrations, the RNase P RNA from the cyanelle of C. paradoxa displays RNA-alone activity well above the detection level. Moreover, the RNA can be complemented to a functional holoenzyme by the E. coli RNase P protein, further supporting its overall bacterial-like architecture. Mutational analysis and domain swaps revealed that this A,U-rich cyanelle RNase P RNA is globally optimized but conformationally unstable, since changes as little as a single point mutation or a base pair identity switch at positions that are not part of the universally conserved catalytic core led to a complete loss of RNA-alone activity. Likely related to this low robustness, extensive structural changes towards an E. coli-type P5-7/P15-17 subdomain as a canonical interaction site for tRNA 3'-CCA termini could not be coaxed into increased ribozyme activity.

Function of heterologous and truncated RNase P proteins in Bacillus subtilis

Bacterial RNase P is composed of an RNA subunit and a single protein (encoded by the rnpB and rnpA genes respectively). The Bacillus subtilis rnpA knockdown strain d7 was used to screen for functional conservation among bacterial RNase P proteins from a representative spectrum of bacterial subphyla. We demonstrate conserved function of bacterial RNase P (RnpA) proteins despite low sequence conservation. Even rnpA genes from psychrophilic and thermophilic bacteria rescued growth of B. subtilis d7 bacteria; likewise, terminal extensions and insertions between beta strands 2 and 3, in the so-called metal binding loop, were compatible with RnpA function in B. subtilis. A deletion analysis of B. subtilis RnpA defined the structural elements essential for bacterial RNase P function in vivo. We further extended our complementation analysis in B. subtilis strain d7 to the four individual RNase P protein subunits from three different Archaea, as well as to human Rpp21 and Rpp29 as representatives of eukaryal RNase P. None of these non-bacterial RNase P proteins showed any evidence of being able to replace the B. subtilis RNase P protein in vivo, supporting the notion that archaeal/eukaryal RNase P proteins are evolutionary unrelated to the bacterial RnpA protein.

Treasure hunt in an amoeba: non-coding RNAs in Dictyostelium discoideum

The traditional view of RNA being merely an intermediate in the transfer of genetic information, as mRNA, spliceosomal RNA, tRNA, and rRNA, has become outdated. The recent discovery of numerous regulatory RNAs with a plethora of functions in biological processes has truly revolutionized our understanding of gene regulation. Tiny RNAs such as microRNAs and small interfering RNAs play vital roles at different levels of gene control. Small nucleolar RNAs are much more abundant than previously recognized, and new functions beyond processing and modification of rRNA have recently emerged. Longer non-coding RNAs (ncRNAs) can also have important regulatory roles in the cell, e.g., antisense RNAs that control their target mRNAs. The majority of these important findings arose from analyses in various model organisms. In this review, we focus on ncRNAs in the social amoeba Dictyostelium discoideum. This important genetically tractable model organism has recently received renewed attention in terms of discovery, regulation and functional studies of ncRNAs. Old and recent findings are discussed and put in context of what we today know about ncRNAs in other organisms.

Combined experimental and computational approach to identify non-protein-coding RNAs in the deep-branching eukaryote Giardia intestinalis

Non-protein-coding RNAs represent a large proportion of transcribed sequences in eukaryotes. These RNAs often function in large RNA-protein complexes, which are catalysts in various RNA-processing pathways. As RNA processing has become an increasingly important area of research, numerous non-messenger RNAs have been uncovered in all the model eukaryotic organisms. However, knowledge on RNA processing in deep-branching eukaryotes is still limited. This study focuses on the identification of non-protein-coding RNAs from the diplomonad parasite Giardia intestinalis, showing that a combined experimental and computational search strategy is a fast method of screening reduced or compact genomes. The analysis of our Giardia cDNA library has uncovered 31 novel candidates, including C/D-box and H/ACA box snoRNAs, as well as an unusual transcript of RNase P, and double-stranded RNAs. Subsequent computational analysis has revealed additional putative C/D-box snoRNAs. Our results will lead towards a future understanding of RNA metabolism in the deep-branching eukaryote Giardia, as more ncRNAs are characterized.

DRpp20 and DRpp40: Two protein subunits involved in Dictyostelium discoideum ribonuclease P holoenzyme assembly

Ribonuclease P is an essential enzyme that matures the 5' ends of all primary tRNA transcripts. RNase P enzymes contain a similar in size RNA subunit which is absolutely required for catalysis. The holoenzyme from Dictyostelium discoideum possesses an essential for activity RNA subunit but the exact protein composition is still under investigation. Bioinformatic analysis of D. discoideum sequencing data returned seven ORFs homologous to previously characterized RNase P protein subunits from human. In the present study, DRpp20 and DRpp40 were cloned and characterized. These proteins apart from the noted similarity possess idiosyncratic regions. Immunobiochemical analysis presented herein indicates their direct involvement in the formation of the ribonucleoprotein complex of D. discoideum RNase P holoenzyme.

An important piece of the RNase P jigsaw solved

RNase P is the ribonucleoprotein enzyme that generates the mature 5' ends of tRNAs throughout all three kingdoms of life. Long known to function as a ribozyme in bacteria and several archaea, it has remained unclear if eukaryal RNase P has entirely lost this RNA-alone catalytic capacity (i.e. the ability to perform catalysis even if the protein part of the enzyme is removed). This controversial debate has now ended after the recent demonstration that eukaryal RNase P also exhibits ribozyme activity.

Identification of conserved secondary structures and expansion segments in enod40 RNAs reveals new enod40 homologues in plants

enod40 is a plant gene that participates in the regulation of symbiotic interaction between leguminous plants and bacteria or fungi. Furthermore, it has been suggested to play a general role in non-symbiotic plant development. Although enod40 seems to have multiple functions, being present in many land plants, the molecular mechanisms of its activity are unclear; they may be determined though, by short peptides and/or RNA structures encoded in the enod40 genes. We utilized conserved RNA structures in enod40 sequences to search nucleotide sequence databases and identified a number of new enod40 homologues in plant species that belong to known, but also, to yet unknown enod40-containing plant families. RNA secondary structure predictions and comparative sequence analysis of enod40 RNAs allowed us to determine the most conserved structural features, present in all known enod40 genes. Remarkably, the topology and evolution of one of the conserved structural domains are similar to those of the expansion segments found in structural RNAs such as rRNAs, RNase P and SRP RNAs. Surprisingly, the enod40 RNA structural elements are much more stronger conserved than the encoded peptides. This finding suggests that some general functions of enod40 gene could be determined by the encoded RNA structure, whereas short peptides may be responsible for more diverse functions found only in certain plant families.

A 40.7 kDa Rpp30/Rpp1 homologue is a protein subunit of Dictyostelium discoideum RNase P holoenzyme

RNase P is an essential and ubiquitous endonuclease that mediates the maturation of the 5' ends of all precursor tRNA molecules. The holoenzyme from Dictyostelium discoideum possesses RNA and protein subunits essential for activity, but the exact composition of the ribonucleoprotein complex is still under investigation. Bioinformatic analysis of D. discoideum genome identified seven open reading frames encoding candidate RNase P protein subunits. The gene named drpp30 encodes a protein with a predicted molecular mass of 40.7 kDa that clusters with Rpp1 and Rpp30 RNase P protein subunits from Saccharomyces cerevisiae and human respectively, which have significantly lower molecular masses. Cloning and heterologous expression of DRpp30 followed by immunochemical analysis of RNase P active fractions demonstrates its association with RNase P holoenzyme. Furthermore, we show that DRpp30 can bind D. discoideum RNase P RNA and tRNA transcripts in vitro, giving a first insight of its possible role in D. discoideum RNase P function. Homology modeling using as a template the archaeal Ph1887p, and molecular dynamics simulations of the modeled structure suggest that DRpp30 adopts a TIM-barrel fold.

RNase MRP and the RNA processing cascade in the eukaryotic ancestor

BACKGROUND

Within eukaryotes there is a complex cascade of RNA-based macromolecules that process other RNA molecules, especially mRNA, tRNA and rRNA. An example is RNase MRP processing ribosomal RNA (rRNA) in ribosome biogenesis. One hypothesis is that this complexity was present early in eukaryotic evolution; an alternative is that an initial simpler network later gained complexity by gene duplication in lineages that led to animals, fungi and plants. Recently there has been a rapid increase in support for the complexity-early theory because the vast majority of these RNA-processing reactions are found throughout eukaryotes, and thus were likely to be present in the last common ancestor of living eukaryotes, herein called the Eukaryotic Ancestor.

RESULTS

We present an overview of the RNA processing cascade in the Eukaryotic Ancestor and investigate in particular, RNase MRP which was previously thought to have evolved later in eukaryotes due to its apparent limited distribution in fungi and animals and plants. Recent publications, as well as our own genomic searches, find previously unknown RNase MRP RNAs, indicating that RNase MRP has a wide distribution in eukaryotes. Combining secondary structure and promoter region analysis of RNAs for RNase MRP, along with analysis of the target substrate (rRNA), allows us to discuss this distribution in the light of eukaryotic evolution.

CONCLUSION

We conclude that RNase MRP can now be placed in the RNA-processing cascade of the Eukaryotic Ancestor, highlighting the complexity of RNA-processing in early eukaryotes. Promoter analyses of MRP-RNA suggest that regulation of the critical processes of rRNA cleavage can vary, showing that even these key cellular processes (for which we expect high conservation) show some species-specific variability. We present our consensus MRP-RNA secondary structure as a useful model for further searches.

Eukaryotic RNase P RNA mediates cleavage in the absence of protein

The universally conserved ribonucleoprotein RNase P is involved in the processing of tRNA precursor transcripts. RNase P consists of one RNA and, depending on its origin, a variable number of protein subunits. Catalytic activity of the RNA moiety so far has been demonstrated only for bacterial and some archaeal RNase P RNAs but not for their eukaryotic counterparts. Here, we show that RNase P RNAs from humans and the lower eukaryote Giardia lamblia mediate cleavage of four tRNA precursors and a model RNA hairpin loop substrate in the absence of protein. Compared with bacterial RNase P RNA, the rate of cleavage (k(obs)) was five to six orders of magnitude lower, whereas the affinity for the substrate (appK(d)) was reduced approximately 20- to 50-fold. We conclude that the RNA-based catalytic activity of RNase P has been preserved during evolution. This finding opens previously undescribed ways to study the role of the different proteins subunits of eukaryotic RNase P.

Structure and function of eukaryotic Ribonuclease P RNA

Ribonuclease P (RNase P) is the ribonucleoprotein endonuclease that processes the 5' ends of precursor tRNAs. Bacterial and eukaryal RNase P RNAs had the same primordial ancestor; however, they were molded differently by evolution. RNase P RNAs of eukaryotes, in contrast to bacterial RNAs, are not catalytically active in vitro without proteins. By comparing the bacterial and eukaryal RNAs, we can begin to understand the transitions made between the RNA and protein-dominated worlds. We report, based on crosslinking studies, that eukaryal RNAs, although catalytically inactive alone, fold into functional forms and specifically bind tRNA even in the absence of proteins. Based on the crosslinking results and crystal structures of bacterial RNAs, we develop a tertiary structure model of the eukaryal RNase P RNA. The eukaryal RNA contains a core structure similar to the bacterial RNA but lacks specific features that in bacterial RNAs contribute to catalysis and global stability of tertiary structure.

Bacterial RNase P: a new view of an ancient enzyme

Ribonuclease P (RNase P) is a ubiquitous endonuclease that catalyses the maturation of the 5' end of transfer RNA (tRNA). Although it carries out a biochemically simple reaction, RNase P is a complex ribonucleoprotein particle composed of a single large RNA and at least one protein component. In bacteria and some archaea, the RNA component of RNase P can catalyse tRNA maturation in vitro in the absence of proteins. The discovery of the catalytic activity of the bacterial RNase P RNA triggered numerous mechanistic and biochemical studies of the reactions catalysed by the RNA alone and by the holoenzyme and, in recent years, structures of individual components of the RNase P holoenzyme have been determined. The goal of the present review is to summarize what is known about the bacterial RNase P, and to bring together the recent structural results with extensive earlier biochemical and phylogenetic findings.

Inventory and analysis of the protein subunits of the ribonucleases P and MRP provides further evidence of homology between the yeast and human enzymes

The RNases P and MRP are involved in tRNA and rRNA processing, respectively. Both enzymes in eukaryotes are composed of an RNA molecule and 9–12 protein subunits. Most of the protein subunits are shared between RNases P and MRP. We have here performed a computational analysis of the protein subunits in a broad range of eukaryotic organisms using profile-based searches and phylogenetic methods. A number of novel homologues were identified, giving rise to a more complete inventory of RNase P/MRP proteins. We present evidence of a relationship between fungal Pop8 and the protein subunit families Rpp14/Pop5 as well as between fungal Pop6 and metazoan Rpp25. These relationships further emphasize a structural and functional similarity between the yeast and human P/MRP complexes. We have also identified novel P and MRP RNAs and analysis of all available sequences revealed a K-turn motif in a large number of these RNAs. We suggest that this motif is a binding site for the Pop3/Rpp38 proteins and we discuss other structural features of the RNA subunit and possible relationships to the protein subunit repertoire.

RNase P: interface of the RNA and protein worlds

Ribonuclease P (RNase P) is an endonuclease involved in processing tRNA. It contains both RNA and protein subunits and occurs in all three domains of life: namely, Archaea, Bacteria and Eukarya. The RNase P RNA subunits from bacteria and some archaea are catalytically active in vitro, whereas those from eukaryotes and most archaea require protein subunits for activity. RNase P has been characterized biochemically and genetically in several systems, and detailed structural information is emerging for both RNA and protein subunits from phylogenetically diverse organisms. In vitro reconstitution of activity is providing insight into the role of proteins in the RNase P holoenzyme. Together, these findings are beginning to impart an understanding of the coevolution of the RNA and protein worlds.

Sequence analysis of RNase MRP RNA reveals its origination from eukaryotic RNase P RNA

RNase MRP is a eukaryote-specific endoribonuclease that generates RNA primers for mitochondrial DNA replication and processes precursor rRNA. RNase P is a ubiquitous endoribonuclease that cleaves precursor tRNA transcripts to produce their mature 5' termini. We found extensive sequence homology of catalytic domains and specificity domains between their RNA subunits in many organisms. In Candida glabrata, the internal loop of helix P3 is 100% conserved between MRP and P RNAs. The helix P8 of MRP RNA from microsporidia Encephalitozoon cuniculi is identical to that of P RNA. Sequence homology can be widely spread over the whole molecule of MRP RNA and P RNA, such as those from Dictyostelium discoideum. These conserved nucleotides between the MRP and P RNAs strongly support the hypothesis that the MRP RNA is derived from the P RNA molecule in early eukaryote evolution.

Ribonuclease P: the evolution of an ancient RNA enzyme

Ribonuclease P (RNase P) is an ancient and essential endonuclease that catalyses the cleavage of the 5' leader sequence from precursor tRNAs (pre-tRNAs). The enzyme is one of only two ribozymes which can be found in all kingdoms of life (Bacteria, Archaea, and Eukarya). Most forms of RNase P are ribonucleoproteins; the bacterial enzyme possesses a single catalytic RNA and one small protein. However, in archaea and eukarya the enzyme has evolved an increasingly more complex protein composition, whilst retaining a structurally related RNA subunit. The reasons for this additional complexity are not currently understood. Furthermore, the eukaryotic RNase P has evolved into several different enzymes including a nuclear activity, organellar activities, and the evolution of a distinct but closely related enzyme, RNase MRP, which has different substrate specificities, primarily involved in ribosomal RNA biogenesis. Here we examine the relationship between the bacterial and archaeal RNase P with the eukaryotic enzyme, and summarize recent progress in characterizing the archaeal enzyme. We review current information regarding the nuclear RNase P and RNase MRP enzymes in the eukaryotes, focusing on the relationship between these enzymes by examining their composition, structure and functions.

Type A and B RNase P RNAs are interchangeable in vivo despite substantial biophysical differences

We show that structural type A and B bacterial ribonuclease P (RNase P) RNAs can fully replace each other in vivo despite the many reported differences in their biogenesis, biochemical/biophysical properties and enzyme function in vitro. Our findings suggest that many of the reported idiosyncrasies of type A and B enzymes either do not reflect the in vivo situation or are not crucial for RNase P function in vivo, at least under standard growth conditions. The discrimination of mature tRNA by RNase P, so far thought to prevent product inhibition of the enzyme in the presence of a large cellular excess of mature tRNA relative to the precursor form, is apparently not crucial for RNase P function in vivo.

Computational analysis of RNAs

Genome sequence analysis of RNAs presents special challenges to computational biology, because conserved RNA secondary structure plays a large part in RNA analysis. Algorithms well suited for RNA secondary structure and sequence analysis have been borrowed from computational linguistics. These "stochastic context-free grammar" (SCFG) algorithms have enabled the development of new RNA gene-finding and RNA homology search software. The aim of this paper is to provide an accessible introduction to the strengths and weaknesses of SCFG methods and to describe the state of the art in one particular kind of application: SCFG-based RNA similarity searching. The INFERNAL and RSEARCH programs are capable of identifying distant RNA homologs in a database search by looking for both sequence and secondary structure conservation.

Secondary structure probing of the human RNase MRP RNA reveals the potential for MRP RNA subsets

RNase MRP is a ribonucleoprotein endoribonuclease involved in eukaryotic pre-rRNA processing. The enzyme possesses an RNA subunit, structurally related to that of RNase P RNA, that is thought to be catalytic. RNase MRP RNA sequences from Saccharomycetaceae species are structurally well defined through detailed phylogenetic and structural analysis. In contrast, higher eukaryote MRP RNA structure models are based on comparative sequence analysis of only five sequences and limited probing data. Detailed structural analysis of the Homo sapiens MRP RNA, entailing enzymatic and chemical probing, is reported. The data are consistent with the phylogenetic secondary structure model and demonstrate unequivocally that higher eukaryote MRP RNA structure differs significantly from that reported for Saccharomycetaceae species. Neither model can account for all of the known MRP RNAs and we thus propose the evolution of at least two subsets of RNase MRP secondary structure, differing predominantly in the predicted specificity domain.

Identification and analysis of ribonuclease P and MRP RNA in a broad range of eukaryotes

RNases P and MRP are ribonucleoprotein complexes involved in tRNA and rRNA processing, respectively. The RNA subunits of these two enzymes are structurally related to each other and play an essential role in the enzymatic reaction. Both of the RNAs have a highly conserved helical region, P4, which is important in the catalytic reaction. We have used a bioinformatics approach based on conserved elements to computationally analyze available genomic sequences of eukaryotic organisms and have identified a large number of novel nuclear RNase P and MRP RNA genes. For MRP RNA for instance, this investigation increases the number of known sequences by a factor of three. We present secondary structure models of many of the predicted RNAs. Although all sequences are able to fold into the consensus secondary structure of P and MRP RNAs, a striking variation in size is observed, ranging from a Nosema locustae MRP RNA of 160 nt to much larger RNAs, e.g. a Plasmodium knowlesi P RNA of 696 nt. The P and MRP RNA genes appear in tandem in some protists, further emphasizing the close evolutionary relationship of these RNAs.