Discovery of a minimal form of RNase P in Pyrobaculum

Coevolution of RNA and protein subunits in RNase P and RNase MRP, two RNA processing enzymes

RNase P and RNase mitochondrial RNA processing (MRP) are ribonucleoproteins (RNPs) that consist of a catalytic RNA and a varying number of protein cofactors. RNase P is responsible for precursor tRNA maturation in all three domains of life, while RNase MRP, exclusive to eukaryotes, primarily functions in rRNA biogenesis. While eukaryotic RNase P is associated with more protein cofactors and has an RNA subunit with fewer auxiliary structural elements compared to its bacterial cousin, the double-anchor precursor tRNA recognition mechanism has remarkably been preserved during evolution. RNase MRP shares evolutionary and structural similarities with RNase P, preserving the catalytic core within the RNA moiety inherited from their common ancestor. By incorporating new protein cofactors and RNA elements, RNase MRP has established itself as a distinct RNP capable of processing ssRNA substrates. The structural information on RNase P and MRP helps build an evolutionary trajectory, depicting how emerging protein cofactors harmonize with the evolution of RNA to shape different functions for RNase P and MRP. Here, we outline the structural and functional relationship between RNase P and MRP to illustrate the coevolution of RNA and protein cofactors, a key driver for the extant, diverse RNP world.

Multiple structural flavors of RNase P in precursor tRNA processing

The precursor transfer RNAs (pre-tRNAs) require extensive processing to generate mature tRNAs possessing proper fold, structural stability, and functionality required to sustain cellular viability. The road to tRNA maturation follows an ordered process: 5'-processing, 3'-processing, modifications at specific sites, if any, and 3'-CCA addition before aminoacylation and recruitment to the cellular protein synthesis machinery. Ribonuclease P (RNase P) is a universally conserved endonuclease in all domains of life, performing the hydrolysis of pre-tRNA sequences at the 5' end by the removal of phosphodiester linkages between nucleotides at position -1 and +1. Except for an archaeal species: Nanoarchaeum equitans where tRNAs are transcribed from leaderless-position +1, RNase P is indispensable for life and displays fundamental variations in terms of enzyme subunit composition, mechanism of substrate recognition and active site architecture, utilizing in all cases a two metal ion-mediated conserved catalytic reaction. While the canonical RNA-based ribonucleoprotein RNase P has been well-known to occur in bacteria, archaea, and eukaryotes, the occurrence of RNA-free protein-only RNase P in eukaryotes and RNA-free homologs of Aquifex RNase P in prokaryotes has been discovered more recently. This review aims to provide a comprehensive overview of structural diversity displayed by various RNA-based and RNA-free RNase P holoenzymes towards harnessing critical RNA-protein and protein-protein interactions in achieving conserved pre-tRNA processing functionality. Furthermore, alternate roles and functional interchangeability of RNase P are discussed in the context of its employability in several clinical and biotechnological applications. This article is categorized under: RNA Processing > tRNA Processing RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Elucidation of structure-function relationships in Methanocaldococcus jannaschii RNase P, a multi-subunit catalytic ribonucleoprotein

RNase P is a ribonucleoprotein (RNP) that catalyzes removal of the 5' leader from precursor tRNAs in all domains of life. A recent cryo-EM study of Methanocaldococcus jannaschii (Mja) RNase P produced a model at 4.6-Å resolution in a dimeric configuration, with each holoenzyme monomer containing one RNase P RNA (RPR) and one copy each of five RNase P proteins (RPPs; POP5, RPP30, RPP21, RPP29, L7Ae). Here, we used native mass spectrometry (MS), mass photometry (MP), and biochemical experiments that (i) validate the oligomeric state of the Mja RNase P holoenzyme in vitro, (ii) find a different stoichiometry for each holoenzyme monomer with up to two copies of L7Ae, and (iii) assess whether both L7Ae copies are necessary for optimal cleavage activity. By mutating all kink-turns in the RPR, we made the discovery that abolishing the canonical L7Ae-RPR interactions was not detrimental for RNase P assembly and function due to the redundancy provided by protein-protein interactions between L7Ae and other RPPs. Our results provide new insights into the architecture and evolution of RNase P, and highlight the utility of native MS and MP in integrated structural biology approaches that seek to augment the information obtained from low/medium-resolution cryo-EM models.

The Dynamic Network of RNP RNase P Subunits

Ribonuclease P (RNase P) is an important ribonucleoprotein (RNP), responsible for the maturation of the 5′ end of precursor tRNAs (pre-tRNAs). In all organisms, the cleavage activity of a single phosphodiester bond adjacent to the first nucleotide of the acceptor stem is indispensable for cell viability and lies within an essential catalytic RNA subunit. Although RNase P is a ribozyme, its kinetic efficiency in vivo, as well as its structural variability and complexity throughout evolution, requires the presence of one protein subunit in bacteria to several protein partners in archaea and eukaryotes. Moreover, the existence of protein-only RNase P (PRORP) enzymes in several organisms and organelles suggests a more complex evolutionary timeline than previously thought. Recent detailed structures of bacterial, archaeal, human and mitochondrial RNase P complexes suggest that, although apparently dissimilar enzymes, they all recognize pre-tRNAs through conserved interactions. Interestingly, individual protein subunits of the human nuclear and mitochondrial holoenzymes have additional functions and contribute to a dynamic network of elaborate interactions and cellular processes. Herein, we summarize the role of each RNase P subunit with a focus on the human nuclear RNP and its putative role in flawless gene expression in light of recent structural studies.

The many faces of RNA-based RNase P, an RNA-world relic

RNase P is an essential enzyme that catalyzes removal of the 5' leader from precursor transfer RNAs. The ribonucleoprotein (RNP) form of RNase P is present in all domains of life and comprises a single catalytic RNA (ribozyme) and a variable number of protein cofactors. Recent cryo-electron microscopy structures of representative archaeal and eukaryotic (nuclear) RNase P holoenzymes bound to tRNA substrate/product provide high-resolution detail on subunit organization, topology, and substrate recognition in these large, multisubunit catalytic RNPs. These structures point to the challenges in understanding how proteins modulate the RNA functional repertoire and how the structure of an ancient RNA-based catalyst was reshaped during evolution by new macromolecular associations that were likely necessitated by functional/regulatory coupling.

Protein cofactors and substrate influence Mg2+-dependent structural changes in the catalytic RNA of archaeal RNase P

The ribonucleoprotein (RNP) form of archaeal RNase P comprises one catalytic RNA and five protein cofactors. To catalyze Mg²⁺-dependent cleavage of the 5′ leader from pre-tRNAs, the catalytic (C) and specificity (S) domains of the RNase P RNA (RPR) cooperate to recognize different parts of the pre-tRNA. While ∼250–500 mM Mg²⁺ renders the archaeal RPR active without RNase P proteins (RPPs), addition of all RPPs lowers the Mg²⁺ requirement to ∼10–20 mM and improves the rate and fidelity of cleavage. To understand the Mg²⁺- and RPP-dependent structural changes that increase activity, we used pre-tRNA cleavage and ensemble FRET assays to characterize inter-domain interactions in Pyrococcus furiosus (Pfu) RPR, either alone or with RPPs ± pre-tRNA. Following splint ligation to doubly label the RPR (Cy3-RPR_{C domain} and Cy5-RPR_{S domain}), we used native mass spectrometry to verify the final product. We found that FRET correlates closely with activity, the Pfu RPR and RNase P holoenzyme (RPR + 5 RPPs) traverse different Mg²⁺-dependent paths to converge on similar functional states, and binding of the pre-tRNA by the holoenzyme influences Mg²⁺ cooperativity. Our findings highlight how Mg²⁺ and proteins in multi-subunit RNPs together favor RNA conformations in a dynamic ensemble for functional gains.

A method for achieving complete microbial genomes and improving bins from metagenomics data

Metagenomics facilitates the study of the genetic information from uncultured microbes and complex microbial communities. Assembling complete genomes from metagenomics data is difficult because most samples have high organismal complexity and strain diversity. Some studies have attempted to extract complete bacterial, archaeal, and viral genomes and often focus on species with circular genomes so they can help confirm completeness with circularity. However, less than 100 circularized bacterial and archaeal genomes have been assembled and published from metagenomics data despite the thousands of datasets that are available. Circularized genomes are important for (1) building a reference collection as scaffolds for future assemblies, (2) providing complete gene content of a genome, (3) confirming little or no contamination of a genome, (4) studying the genomic context and synteny of genes, and (5) linking protein coding genes to ribosomal RNA genes to aid metabolic inference in 16S rRNA gene sequencing studies. We developed a semi-automated method called Jorg to help circularize small bacterial, archaeal, and viral genomes using iterative assembly, binning, and read mapping. In addition, this method exposes potential misassemblies from k-mer based assemblies. We chose species of the Candidate Phyla Radiation (CPR) to focus our initial efforts because they have small genomes and are only known to have one ribosomal RNA operon. In addition to 34 circular CPR genomes, we present one circular Margulisbacteria genome, one circular Chloroflexi genome, and two circular megaphage genomes from 19 public and published datasets. We demonstrate findings that would likely be difficult without circularizing genomes, including that ribosomal genes are likely not operonic in the majority of CPR, and that some CPR harbor diverged forms of RNase P RNA. Code and a tutorial for this method is available at https://github.com/lmlui/Jorg and is available on the DOE Systems Biology KnowledgeBase as a beta app.

P finder: genomic and metagenomic annotation of RNase P RNA gene (rnpB)

Background

The rnpB gene encodes for an essential catalytic RNA (RNase P). Like other essential RNAs, RNase P’s sequence is highly variable. However, unlike other essential RNAs (i.e. tRNA, 16 S, 6 S,...) its structure is also variable with at least 5 distinct structure types observed in prokaryotes. This structural variability makes it labor intensive and challenging to create and maintain covariance models for the detection of RNase P RNA in genomic and metagenomic sequences. The lack of a facile and rapid annotation algorithm has led to the rnpB gene being the most grossly under annotated essential gene in completed prokaryotic genomes with only a 24% annotation rate. Here we describe the coupling of the largest RNase P RNA database with the local alignment scoring algorithm to create the most sensitive and rapid prokaryote rnpB gene identification and annotation algorithm to date.

Results

Of the 2772 completed microbial genomes downloaded from GenBank only 665 genomes had an annotated rnpB gene. We applied P Finder to these genomes and were able to identify 2733 or nearly 99% of the 2772 microbial genomes examined. From these results four new rnpB genes that encode the minimal T-type P RNase P RNAs were identified computationally for the first time. In addition, only the second C-type RNase P RNA was identified in Sphaerobacter thermophilus. Of special note, no RNase P RNAs were detected in several obligate endosymbionts of sap sucking insects suggesting a novel evolutionary adaptation.

Conclusions

The coupling of the largest RNase P RNA database and associated structure class identification with the P Finder algorithm is both sensitive and rapid, yielding high quality results to aid researchers annotating either genomic or metagenomic data. It is the only algorithm to date that can identify challenging RNAse P classes such as C-type and the minimal T-type RNase P RNAs. P Finder is written in C# and has a user-friendly GUI that can run on multiple 64-bit windows platforms (Windows Vista/7/8/10). P Finder is free available for download at https://github.com/JChristopherEllis/P-Finder as well as a small sample RNase P RNA file for testing.

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.

Cryo-electron microscopy structure of an archaeal ribonuclease P holoenzyme

Ribonuclease P (RNase P) is an essential ribozyme responsible for tRNA 5′ maturation. Here we report the cryo-EM structures of Methanocaldococcus jannaschii (Mja) RNase P holoenzyme alone and in complex with a tRNA substrate at resolutions of 4.6 Å and 4.3 Å, respectively. The structures reveal that the subunits of MjaRNase P are strung together to organize the holoenzyme in a dimeric conformation required for efficient catalysis. The structures also show that archaeal RNase P is a functional chimera of bacterial and eukaryal RNase Ps that possesses bacterial-like two RNA-based anchors and a eukaryal-like protein-aided stabilization mechanism. The 3′-RCCA sequence of tRNA, which is a key recognition element for bacterial RNase P, is dispensable for tRNA recognition by MjaRNase P. The overall organization of MjaRNase P, particularly within the active site, is similar to those of bacterial and eukaryal RNase Ps, suggesting a universal catalytic mechanism for all RNase Ps.

Ribonulease P is a conserved ribozyme present in all kingdoms of life that is involved in the 5′ maturation step of tRNAs. Here the authors determine the structure of an archaeal RNase P holoenzyme that reveals how archaeal RNase P recognizes its tRNA substrate and suggest a conserved catalytic mechanism amongst RNase Ps despite structural variability.

Insights into RNA-processing pathways and associated RNA-degrading enzymes in Archaea

RNA-processing pathways are at the centre of regulation of gene expression. All RNA transcripts undergo multiple maturation steps in addition to covalent chemical modifications to become functional in the cell. This includes destroying unnecessary or defective cellular RNAs. In Archaea, information on mechanisms by which RNA species reach their mature forms and associated RNA-modifying enzymes are still fragmentary. To date, most archaeal actors and pathways have been proposed in light of information gathered from Bacteria and Eukarya. In this context, this review provides a state of the art overview of archaeal endoribonucleases and exoribonucleases that cleave and trim RNA species and also of the key small archaeal proteins that bind RNAs. Furthermore, synthetic up-to-date views of processing and biogenesis pathways of archaeal transfer and ribosomal RNAs as well as of maturation of stable small non-coding RNAs such as CRISPR RNAs, small C/D and H/ACA box guide RNAs, and other emerging classes of small RNAs are described. Finally, prospective post-transcriptional mechanisms to control archaeal messenger RNA quality and quantity are discussed.

Insight into the functional role of unique determinants in RNA component of RNase P of Mycobacterium tuberculosis

RNase P, an essential ribonucleoprotein enzyme is involved in processing 5' end of pre-tRNA molecules. All bacterial RNase P holoenzymes, including that of Mycobacterim tuberculosis, an important human pathogen contain a catalytically active RNA subunit and a protein subunit. However, the mycobacterial RNA is larger than typical bacterial RNase P RNAs. It contains the essential core structure and many unique features in the peripheral elements. In the current study, an extensive mutational analysis was performed to analyze the function of the unique features in P12, P15.A, P18 and P19 helices in the mycobacterial RNase P RNA. The study demonstrates that P12 interacts with monovalent and divalent ions and is important for the function of mycobacterial holoenzyme. The helices, P15.A and P18 appear to interact with ammonium and magnesium ions, respectively. P19 is involved in the thermostability of the RNA component as well as interaction with ammonium ions. A homology model of M. tuberculosis RNase P RNA indicates many new inter- and intra-helical interactions. The significance of the unique interactions paves way towards understanding the differential functioning of M. tuberculosis RNase P RNA, for exploring specific inhibition of the same in the pathogen to contain infection.

Critical domain interactions for type A RNase P RNA catalysis with and without the specificity domain

The natural trans-acting ribozyme RNase P RNA (RPR) is composed of two domains in which the catalytic (C-) domain mediates cleavage of various substrates. The C-domain alone, after removal of the second specificity (S-) domain, catalyzes this reaction as well, albeit with reduced efficiency. Here we provide experimental evidence indicating that efficient cleavage mediated by the Escherichia coli C-domain (Eco CP RPR) with and without the C5 protein likely depends on an interaction referred to as the "P6-mimic". Moreover, the P18 helix connects the C- and S-domains between its loop and the P8 helix in the S-domain (the P8/ P18-interaction). In contrast to the "P6-mimic", the presence of P18 does not contribute to the catalytic performance by the C-domain lacking the S-domain in cleavage of an all ribo model hairpin loop substrate while deletion or disruption of the P8/ P18-interaction in full-size RPR lowers the catalytic efficiency in cleavage of the same model hairpin loop substrate in keeping with previously reported data using precursor tRNAs. Consistent with that P18 is not required for cleavage mediated by the C-domain we show that the archaeal Pyrococcus furiosus RPR C-domain, which lacks the P18 helix, is catalytically active in trans without the S-domain and any protein. Our data also suggest that the S-domain has a larger impact on catalysis for E. coli RPR compared to P. furiosus RPR. Finally, we provide data indicating that the absence of the S-domain and P18, or the P8/ P18-interaction in full-length RPR influences the charge distribution near the cleavage site in the RPR-substrate complex to a small but reproducible extent.

A novel double kink-turn module in euryarchaeal RNase P RNAs

RNase P is primarily responsible for the 5΄ maturation of transfer RNAs (tRNAs) in all domains of life. Archaeal RNase P is a ribonucleoprotein made up of one catalytic RNA and five protein cofactors including L7Ae, which is known to bind the kink-turn (K-turn), an RNA structural element that causes axial bending. However, the number and location of K-turns in archaeal RNase P RNAs (RPRs) are unclear. As part of an integrated approach, we used native mass spectrometry to assess the number of L7Ae copies that bound the RPR and site-specific hydroxyl radical-mediated footprinting to localize the K-turns. Mutagenesis of each of the putative K-turns singly or in combination decreased the number of bound L7Ae copies, and either eliminated or changed the L7Ae footprint on the mutant RPRs. In addition, our results support an unprecedented ‘double K-turn’ module in type A and type M archaeal RPR variants.

Minimal and RNA-free RNase P in Aquifex aeolicus

RNase P is an essential tRNA-processing enzyme in all domains of life. We identified an unknown type of protein-only RNase P in the hyperthermophilic bacterium Aquifex aeolicus: Without an RNA subunit and the smallest of its kind, the 23-kDa polypeptide comprises a metallonuclease domain only. The protein has RNase P activity in vitro and rescued the growth of Escherichia coli and Saccharomyces cerevisiae strains with inactivations of their more complex and larger endogenous ribonucleoprotein RNase P. Homologs of Aquifex RNase P (HARP) were identified in many Archaea and some Bacteria, of which all Archaea and most Bacteria also encode an RNA-based RNase P; activity of both RNase P forms from the same bacterium or archaeon could be verified in two selected cases. Bioinformatic analyses suggest that A. aeolicus and related Aquificaceae likely acquired HARP by horizontal gene transfer from an archaeon.

Structural Roles of Noncoding RNAs in the Heart of Enzymatic Complexes

Over billions of years of evolution, nature has embraced proteins as the major workhorse molecules of the cell. However, nearly every aspect of metabolism is dependent upon how structured RNAs interact with proteins, ligands, and other nucleic acids. Key processes, including telomere maintenance, RNA processing, and protein synthesis, require large RNAs that assemble into elaborate three-dimensional shapes. These RNAs can (i) act as flexible scaffolds for protein subunits, (ii) participate directly in substrate recognition, and (iii) serve as catalytic components. Here, we juxtapose the near atomic level interactions of three ribonucleoprotein complexes: ribonuclease P (involved in 5' pre-tRNA processing), the spliceosome (responsible for pre-mRNA splicing), and telomerase (an RNA-directed DNA polymerase that extends the ends of chromosomes). The focus of this perspective is profiling the structural and dynamic roles of RNAs at the core of these enzymes, highlighting how large RNAs contribute to molecular recognition and catalysis.

Cleavage of Model Substrates by Arabidopsis thaliana PRORP1 Reveals New Insights into Its Substrate Requirements

Two broad classes of RNase P trim the 5' leader of precursor tRNAs (pre-tRNAs): ribonucleoprotein (RNP)- and proteinaceous (PRORP)-variants. These two RNase P types, which use different scaffolds for catalysis, reflect independent evolutionary paths. While the catalytic RNA-based RNP form is present in all three domains of life, the PRORP family is restricted to eukaryotes. To obtain insights on substrate recognition by PRORPs, we examined the 5' processing ability of recombinant Arabidopsis thaliana PRORP1 (AtPRORP1) using a panel of pre-tRNASer variants and model hairpin-loop derivatives (pATSer type) that consist of the acceptor-T-stem stack and the T-/D-loop. Our data indicate the importance of the identity of N-1 (the residue immediately 5' to the cleavage site) and the N-1:N+73 base pair for cleavage rate and site selection of pre-tRNASer and pATSer. The nucleobase preferences that we observed mirror the frequency of occurrence in the complete suite of organellar pre-tRNAs in eight algae/plants that we analyzed. The importance of the T-/D-loop in pre-tRNASer for tight binding to AtPRORP1 is indicated by the 200-fold weaker binding of pATSer compared to pre-tRNASer, while the essentiality of the T-loop for cleavage is reflected by the near-complete loss of activity when a GAAA-tetraloop replaced the T-loop in pATSer. Substituting the 2'-OH at N-1 with 2'-H also resulted in no detectable cleavage, hinting at the possible role of this 2'-OH in coordinating Mg2+ ions critical for catalysis. Collectively, our results indicate similarities but also key differences in substrate recognition by the bacterial RNase P RNP and AtPRORP1: while both forms exploit the acceptor-T-stem stack and the elbow region in the pre-tRNA, the RNP form appears to require more recognition determinants for cleavage-site selection.

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.

Sequence Analysis and Comparative Study of the Protein Subunits of Archaeal RNase P

RNase P, a ribozyme-based ribonucleoprotein (RNP) complex that catalyzes tRNA 5′-maturation, is ubiquitous in all domains of life, but the evolution of its protein components (RNase P proteins, RPPs) is not well understood. Archaeal RPPs may provide clues on how the complex evolved from an ancient ribozyme to an RNP with multiple archaeal and eukaryotic (homologous) RPPs, which are unrelated to the single bacterial RPP. Here, we analyzed the sequence and structure of archaeal RPPs from over 600 available genomes. All five RPPs are found in eight archaeal phyla, suggesting that these RPPs arose early in archaeal evolutionary history. The putative ancestral genomic loci of archaeal RPPs include genes encoding several members of ribosome, exosome, and proteasome complexes, which may indicate coevolution/coordinate regulation of RNase P with other core cellular machineries. Despite being ancient, RPPs generally lack sequence conservation compared to other universal proteins. By analyzing the relative frequency of residues at every position in the context of the high-resolution structures of each of the RPPs (either alone or as functional binary complexes), we suggest residues for mutational analysis that may help uncover structure-function relationships in RPPs.

Hammerhead Ribozymes in Archaeal Genomes: A Computational Hunt

Hammerhead ribozymes (HHRs) are small self-cleaving RNAs, first discovered in viroids and satellite RNAs of plant viruses. They are composed of a catalytic core of conserved nucleotides flanked by three helices. More recently, hammerhead-encoding sequences have been identified in the genomes of many eukaryotes, prokaryotes and other non-viral species regulating various functions. In this study we have explored the Archaeal domain to identify HHRs using three different bioinformatics approach. Our study reveals four putative hits of HHRs type I and type II in the group Thaumarchaeota and Euryarchaeota in the Archaeal domain, one of which is the instance of HHR 1 in C. symbiosum A, already identified in a previous study. These HHRs are very similar to those previously described in terms of the conservation of their catalytic core. Based on 3-D structure analysis and free energy, these instances were concluded as putative HHRs. Our findings reveal that the catalytic core contains the conserved motifs that are essential for cleavage activity, but there are some instances in which compensatory core variations are present. However, no instances of HHRs have been found in Crenarchaeota. This study reveals a very scarce presence of HHRs in Archaea which suggests the involvement of other ncRNA elements in gene regulatory system like RNase P which are abundantly found in the Archaeal domain.

Conservation and losses of non-coding RNAs in avian genomes

Here we present the results of a large-scale bioinformatics annotation of non-coding RNA loci in 48 avian genomes. Our approach uses probabilistic models of hand-curated families from the Rfam database to infer conserved RNA families within each avian genome. We supplement these annotations with predictions from the tRNA annotation tool, tRNAscan-SE and microRNAs from miRBase. We identify 34 lncRNA-associated loci that are conserved between birds and mammals and validate 12 of these in chicken. We report several intriguing cases where a reported mammalian lncRNA, but not its function, is conserved. We also demonstrate extensive conservation of classical ncRNAs (e.g., tRNAs) and more recently discovered ncRNAs (e.g., snoRNAs and miRNAs) in birds. Furthermore, we describe numerous “losses” of several RNA families, and attribute these to either genuine loss, divergence or missing data. In particular, we show that many of these losses are due to the challenges associated with assembling avian microchromosomes. These combined results illustrate the utility of applying homology-based methods for annotating novel vertebrate genomes.

Transcriptional control of an essential ribozyme in Drosophila reveals an ancient evolutionary divide in animals

Ribonuclease P (RNase P) is an essential enzyme required for 5'-maturation of tRNA. While an RNA-free, protein-based form of RNase P exists in eukaryotes, the ribonucleoprotein (RNP) form is found in all domains of life. The catalytic component of the RNP is an RNA known as RNase P RNA (RPR). Eukaryotic RPR genes are typically transcribed by RNA polymerase III (pol III). Here we showed that the RPR gene in Drosophila, which is annotated in the intron of a pol II-transcribed protein-coding gene, lacks signals for transcription by pol III. Using reporter gene constructs that include the RPR-coding intron from Drosophila, we found that the intron contains all the sequences necessary for production of mature RPR but is dependent on the promoter of the recipient gene for expression. We also demonstrated that the intron-coded RPR copurifies with RNase P and is required for its activity. Analysis of RPR genes in various animal genomes revealed a striking divide in the animal kingdom that separates insects and crustaceans into a single group in which RPR genes lack signals for independent transcription and are embedded in different protein-coding genes. Our findings provide evidence for a genetic event that occurred approximately 500 million years ago in the arthropod lineage, which switched the control of the transcription of RPR from pol III to pol II.

Class-specific prediction of ncRNAs

Many RNA families, i.e., groups of homologous RNA genes, belong to RNA classes, such as tRNAs, snoRNAs, or microRNAs, that are characterized by common sequence motifs and/or common secondary structure features. The detection of new members of RNA classes, as well as the comprehensive annotation of genomes with members of RNA classes is a challenging task that goes beyond simple homology search. Computational methods addressing this problem typically use a three-tiered approach: In the first step an efficient and sensitive filter is employed. In the second step the candidate set is narrowed down using computationally expensive methods geared towards specificity. In the final step the hits are annotated with class-specific features and scored. Here we review the tools that are currently available for a diverse set of RNA classes.

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.

Ribonucleases P/MRP and the expanding ribonucleoprotein world

One of the hallmarks of life is the widespread use of certain essential ribozymes. The ubiquitous ribonuclease P (RNase P) and eukaryotic RNase MRP are essential complexes where a structured, noncoding RNA acts in catalysis. Recent discoveries have elucidated the three-dimensional structure of the ancestral ribonucleoprotein complex, suggested the possibility of a protein-only composition in organelles, and even noted the absence of RNase P in a non-free-living organism. With respect to these last two findings, import mechanisms for RNases P/MRP into mitochondria have been demonstrated, and RNase P is present in organisms with some of the smallest known genomes. Together, these results have led to an ongoing debate regarding the precise definition of how "essential" these ribozymes truly are.

Modeling the Thermoproteaceae RNase P RNA

The RNA component of the RNase P complex is found throughout most branches of the tree of life and is principally responsible for removing the 5' leader sequence from pre-tRNA transcripts during tRNA maturation. RNase P RNA has a number of universal core features, however variations in sequence and structure found in homologs across the tree of life require multiple Rfam covariance search models to detect accurately. We describe a new Rfam search model to enable efficient detection of the diminutive archaeal Type T RNase P RNAs, which are missed by existing Rfam models. Using the new model, we establish effective score detection thresholds, and detect four new RNase P RNA genes in recently completed genomes from the crenarchaeal family Thermoproteaceae.

Complete genome sequence of Pyrobaculum oguniense

Pyrobaculum oguniense TE7 is an aerobic hyperthermophilic crenarchaeon isolated from a hot spring in Japan. Here we describe its main chromosome of 2,436,033 bp, with three large-scale inversions and an extra-chromosomal element of 16,887 bp. We have annotated 2,800 protein-coding genes and 145 RNA genes in this genome, including nine H/ACA-like small RNA, 83 predicted C/D box small RNA, and 47 transfer RNA genes. Comparative analyses with the closest known relative, the anaerobe Pyrobaculum arsenaticum from Italy, reveals unexpectedly high synteny and nucleotide identity between these two geographically distant species. Deep sequencing of a mixture of genomic DNA from multiple cells has illuminated some of the genome dynamics potentially shared with other species in this genus.

RNA processing in the minimal organism Nanoarchaeum equitans

BACKGROUND

The minimal genome of the tiny, hyperthermophilic archaeon Nanoarchaeum equitans contains several fragmented genes and revealed unusual RNA processing pathways. These include the maturation of tRNA molecules via the trans-splicing of tRNA halves and genomic rearrangements to compensate for the absence of RNase P.

RESULTS

Here, the RNA processing events in the N. equitans cell are analyzed using RNA-Seq deep sequencing methodology. All tRNA half precursor and tRNA termini were determined and support the tRNA trans-splicing model. The processing of CRISPR RNAs from two CRISPR clusters was verified. Twenty-seven C/D box small RNAs (sRNAs) and a H/ACA box sRNA were identified. The C/D box sRNAs were found to flank split genes, to form dicistronic tRNA-sRNA precursors and to be encoded within the tRNAMet intron.

CONCLUSIONS

The presented data provide an overview of the production and usage of small RNAs in a cell that has to survive with a highly reduced genome. N. equitans lost many essential metabolic pathways but maintains highly active CRISPR/Cas and rRNA modification systems that appear to play an important role in genome fragmentation.

Nuclear RNase P of Trypanosoma brucei: a single protein in place of the multicomponent RNA-protein complex

RNase P is the endonuclease that removes 5′ extensions from tRNA precursors. In its best-known form, the enzyme is composed of a catalytic RNA and a protein moiety variable in number and mass. This ribonucleoprotein enzyme is widely considered ubiquitous and apparently reached its highest complexity in the eukaryal nucleus, where it is typically composed of at least ten subunits. Here, we show that in the protist Trypanosoma brucei, two proteins are the sole forms of RNase P. They localize to the nucleus and the mitochondrion, respectively, and have RNase P activity each on their own. The protein-RNase P is, moreover, capable of replacing nuclear RNase P in yeast cells. This shows that complex ribonucleoprotein structures and RNA catalysis are not necessarily required to support tRNA 5′ end formation in eukaryal cells.

Cleavage mediated by the catalytic domain of bacterial RNase P RNA

Like other RNA molecules, RNase P RNA (RPR) is composed of domains, and these have different functions. Here, we provide data demonstrating that the catalytic (C) domain of Escherichia coli (Eco) RPR when separated from the specificity (S) domain mediates cleavage using various model RNA hairpin loop substrates. Compared to full-length Eco RPR, the rate constant, k(obs), of cleavage for the truncated RPR (CP RPR) was reduced 30- to 13,000-fold depending on substrate. Specifically, the structural architecture of the -1/+73 played a significant role where a C(-1)/G(+73) pair had the most dramatic effect on k(obs). Substitution of A(248) (E. coli numbering), positioned near the cleavage site in the RNase P-substrate complex, with G in the CP RPR resulted in 30-fold improvement in rate. In contrast, strengthening the interaction between the RPR and the 3' end of the substrate only had a modest effect. Interestingly, although deleting the S-domain gave a reduction in the rate, it resulted in a less erroneous RPR with respect to cleavage site selection. These data support and extend our understanding of the coupling between the distal interaction between the S-domain and events at the active site. Our findings will also be discussed with respect to the structure of RPR derived from different organisms.

Of P and Z: mitochondrial tRNA processing enzymes

Mitochondrial tRNAs are generally synthesized as part of polycistronic transcripts. Release of tRNAs from these precursors is thus not only required to produce functional adaptors for translation, but also responsible for the maturation of other mitochondrial RNA species. Cleavage of mitochondrial tRNAs appears to be exclusively accomplished by endonucleases. 5'-end maturation in the mitochondria of different Eukarya is achieved by various kinds of RNase P, representing the full range of diversity found in this enzyme family. While ribonucleoprotein enzymes with RNA components of bacterial-like appearance are found in a few unrelated protists, algae, and fungi, highly degenerate RNAs of dramatic size variability are found in the mitochondria of many fungi. The majority of mitochondrial RNase P enzymes, however, appear to be pure protein enzymes. Human mitochondrial RNase P, the first to be identified and possibly the prototype of all animal mitochondrial RNases P, is composed of three proteins. Homologs of its nuclease subunit MRPP3/PRORP, are also found in plants, algae and several protists, where they are apparently responsible for RNase P activity in mitochondria (and beyond) without the help of extra subunits. The diversity of RNase P enzymes is contrasted by the uniformity of mitochondrial RNases Z, which are responsible for 3'-end processing. Only the long form of RNase Z, which is restricted to eukarya, is found in mitochondria, even when an additional short form is present in the same organism. Mitochondrial tRNA processing thus appears dominated by new, eukaryal inventions rather than bacterial heritage. This article is part of a Special Issue entitled: Mitochondrial Gene Expression.

Cleavage mediated by the P15 domain of bacterial RNase P RNA

Independently folded domains in RNAs frequently adopt identical tertiary structures regardless of whether they are in isolation or are part of larger RNA molecules. This is exemplified by the P15 domain in the RNA subunit (RPR) of the universally conserved endoribonuclease P, which is involved in the processing of tRNA precursors. One of its domains, encompassing the P15 loop, binds to the 3'-end of tRNA precursors resulting in the formation of the RCCA-RNase P RNA interaction (interacting residues underlined) in the bacterial RPR-substrate complex. The function of this interaction was hypothesized to anchor the substrate, expose the cleavage site and result in re-coordination of Mg(2+) at the cleavage site. Here we show that small model-RNA molecules (~30 nt) carrying the P15-loop mediated cleavage at the canonical RNase P cleavage site with significantly reduced rates compared to cleavage with full-size RPR. These data provide further experimental evidence for our model that the P15 domain contributes to both substrate binding and catalysis. Our data raises intriguing evolutionary possibilities for 'RNA-mediated' cleavage of RNA.

RNA binding properties of conserved protein subunits of human RNase P

Human nuclear RNase P is required for transcription and processing of tRNA. This catalytic RNP has an H1 RNA moiety associated with ten distinct protein subunits. Five (Rpp20, Rpp21, Rpp25, Rpp29 and Pop5) out of eight of these protein subunits, prepared in refolded recombinant forms, bind to H1 RNA in vitro. Rpp20 and Rpp25 bind jointly to H1 RNA, even though each protein can interact independently with this transcript. Nuclease footprinting analysis reveals that Rpp20 and Rpp25 recognize overlapping regions in the P2 and P3 domains of H1 RNA. Rpp21 and Rpp29, which are sufficient for reconstitution of the endonucleolytic activity, bind to separate regions in the catalytic domain of H1 RNA. Common themes and discrepancies in the RNA-protein interactions between human nuclear RNase P and its related yeast and archaeal counterparts provide a rationale for the assembly of the fully active form of this enzyme.