HMGB1 prefers to interact with structural RNAs and regulates rRNA methylation modification and translation in HeLa cells

BACKGROUND

High-mobility group B1 (HMGB1) is both a DNA binding nuclear factor modulating transcription and a crucial cytokine that mediates the response to both infectious and noninfectious inflammation such as autoimmunity, cancer, trauma, and ischemia reperfusion injury. HMGB1 has been proposed to control ribosome biogenesis, similar as the other members of a class of HMGB proteins.

RESULTS

Here, we report that HMGB1 selectively promotes transcription of genes involved in the regulation of transcription, osteoclast differentiation and apoptotic process. Improved RNA immunoprecipitation by UV cross-linking and deep sequencing (iRIP-seq) experiment revealed that HMGB1 selectively bound to mRNAs functioning not only in signal transduction and gene expression, but also in axon guidance, focal adhesion, and extracellular matrix organization. Importantly, HMGB1-bound reads were strongly enriched in specific structured RNAs, including the domain II of 28S rRNA, H/ACA box snoRNAs including snoRNA63 and scaRNAs. RTL-P experiment showed that overexpression of HMGB1 led to a decreased methylation modification of 28S rRNA at position Am2388, Cm2409, and Gm2411. We further showed that HMGB1 overexpression increased ribosome RNA expression levels and enhanced protein synthesis.

CONCLUSION

Taken together, our results support a model in which HMGB1 binds to multiple RNA species in human cancer cells, which could at least partially contribute to HMGB1-modulated rRNA modification, protein synthesis function of ribosomes, and differential gene expression including rRNA genes. These findings provide additional mechanistic clues to HMGB1 functions in cancers and cell differentiation.

The Efg1-Bud22 dimer associates with the U14 snoRNP contacting the 5' rRNA domain of an early 90S pre-ribosomal particle

The DEAD-box helicase Dbp4 plays an essential role during the early assembly of the 40S ribosome, which is only poorly understood to date. By applying the yeast two-hybrid method and biochemical approaches, we discovered that Dbp4 interacts with the Efg1-Bud22 dimer. Both factors associate with early pre-90S particles and smaller complexes, each characterized by a high presence of the U14 snoRNA. A crosslink analysis of Bud22 revealed its contact to the U14 snoRNA and the 5' domain of the nascent 18S rRNA, close to its U14 snoRNA hybridization site. Moreover, depletion of Bud22 or Efg1 specifically affects U14 snoRNA association with pre-ribosomal complexes. Accordingly, we concluded that the role of the Efg1-Bud22 dimer is linked to the U14 snoRNA function on early 90S ribosome intermediates chaperoning the 5' domain of the nascent 18S rRNA. The successful rRNA folding of the 5' domain and the release of Efg1, Bud22, Dpb4, U14 snoRNA and associated snoRNP factors allows the subsequent recruitment of the Kre33-Bfr2-Enp2-Lcp5 module towards the 90S pre-ribosome.

Development of an RNA-protein crosslinker to capture protein interactions with diverse RNA structures in cells

Characterization of RNA-protein interaction is fundamental for understanding the metabolism and function of RNA. UV crosslinking has been widely used to map the targets of RNA-binding proteins, but is limited by low efficiency, requirement for zero-distance contact, and biases for single-stranded RNA structure and certain residues of RNA and protein. Here, we report the development of an RNA-protein crosslinker (AMT-NHS) composed of a psoralen derivative and an N-hydroxysuccinimide ester group, which react with RNA bases and primary amines of protein, respectively. We show that AMT-NHS can penetrate into living yeast cells and crosslink Cbf5 to H/ACA snoRNAs with high specificity. The crosslinker induced different crosslinking patterns than UV and targeted both single- and double-stranded regions of RNA. The crosslinker provides a new tool to capture diverse RNA-protein interactions in cells.

SnoRNA guide activities: real and ambiguous

In eukaryotes, rRNAs and spliceosomal snRNAs are heavily modified post-transcriptionally. Pseudouridylation and 2'-O-methylation are the most abundant types of RNA modifications. They are mediated by modification guide RNAs, also known as small nucleolar (sno)RNAs and small Cajal body-specific (sca)RNAs. We used yeast and vertebrate cells to test guide activities predicted for a number of snoRNAs, based on their regions of complementarity with rRNAs. We showed that human SNORA24 is a genuine guide RNA for 18S-Ψ609, despite some noncanonical base-pairing with its target. At the same time, we found quite a few snoRNAs that have the ability to base-pair with rRNAs and can induce predicted modifications in artificial substrate RNAs, but do not modify the same target sequence within endogenous rRNA molecules. Furthermore, certain fragments of rRNAs can be modified by the endogenous yeast modification machinery when inserted into an artificial backbone RNA, even though the same sequences are not modified in endogenous yeast rRNAs. In Xenopus cells, a guide RNA generated from scaRNA, but not from snoRNA, could induce an additional pseudouridylation of U2 snRNA at position 60; both guide RNAs were equally active on a U2 snRNA-specific substrate in yeast cells. Thus, post-transcriptional modification of functionally important RNAs, such as rRNAs and snRNAs, is highly regulated and more complex than simply strong base-pairing between a guide RNA and substrate RNA. We discuss possible regulatory roles for these unexpected modifications.

Discovery of 17 conserved structural RNAs in fungi

Many non-coding RNAs with known functions are structurally conserved: their intramolecular secondary and tertiary interactions are maintained across evolutionary time. Consequently, the presence of conserved structure in multiple sequence alignments can be used to identify candidate functional non-coding RNAs. Here, we present a bioinformatics method that couples iterative homology search with covariation analysis to assess whether a genomic region has evidence of conserved RNA structure. We used this method to examine all unannotated regions of five well-studied fungal genomes (Saccharomyces cerevisiae, Candida albicans, Neurospora crassa, Aspergillus fumigatus, and Schizosaccharomyces pombe). We identified 17 novel structurally conserved non-coding RNA candidates, which include four H/ACA box small nucleolar RNAs, four intergenic RNAs and nine RNA structures located within the introns and untranslated regions (UTRs) of mRNAs. For the two structures in the 3' UTRs of the metabolic genes GLY1 and MET13, we performed experiments that provide evidence against them being eukaryotic riboswitches.

HIV-1 genomic RNA U3 region forms a stable quadruplex-hairpin structure

The U3 promoter region of the HIV-1 long terminal repeat (LTR) has previously been shown to fold into a series of dynamic G-quadruplex structures. Among the G-quadruplexes identified in the LTR sequence, LTR-III was shown to be the most stable in vitro. NMR studies of this 28-nucleotide (nt) DNA revealed a unique quadruplex-hairpin structure. Whether the hairpin forms in RNA element is unknown and the role of the hairpin in the structure and stability of quadruplexes has not been characterized. Here, we used optical and thermodynamic studies to address these questions. The wild-type LTR-III RNA formed a monomolecular quadruplex with a parallel topology using only propeller loops, including the hairpin loop element. By comparison to the WT and variant RNAs, LTR-III DNA structures were more heterogeneous and less stable. Increased stability of the RNA suggests that the RNA quadruplex-hairpin structure may be a more attractive therapeutic target than the analogous DNA element.

Insertions of codons encoding basic amino acids in H7 hemagglutinins of influenza A viruses occur by recombination with RNA at hotspots near snoRNA binding sites

The presence of multiple basic amino acids in the protease cleavage site of the hemagglutinin (HA) protein is the main molecular determinant of virulence of highly pathogenic avian influenza (HPAI) viruses. Recombination of HA RNA with other RNA molecules of host or virus origin is a dominant mechanism of multibasic cleavage site (MBCS) acquisition for H7 subtype HA. Using alignments of HA RNA sequences from documented cases of MBCS insertion due to recombination, we show that such recombination with host RNAs is most likely to occur at particular hotspots in ribosomal RNAs (rRNAs), transfer RNAs (tRNAs), and viral RNAs. The locations of these hotspots in highly abundant RNAs indicate that RNA recombination is facilitated by the binding of small nucleolar RNA (snoRNA) near the recombination points.

scaRNAs and snoRNAs: Are they limited to specific classes of substrate RNAs?

Posttranscriptional modifications of rRNA occur in the nucleolus where rRNA modification guide RNAs, or snoRNAs, concentrate. On the other hand, scaRNAs, the modification guide RNAs for spliceosomal snRNAs, concentrate in the Cajal body (CB). It is generally assumed, therefore, that snRNAs must accumulate in CBs to be modified by scaRNAs. Here we demonstrate that the evidence for the latter postulate is not consistent. In the nucleus, scaRNA localization is not limited to CBs. Furthermore, canonical scaRNAs can modify rRNAs. We suggest that the conventional view that scaRNAs function only in the CB needs revision.

The DEAD-box protein Dbp2p is linked to noncoding RNAs, the helicase Sen1p, and R-loops

The DEAD-box RNA helicase Dbp2p is highly conserved in eukaryotes and has been implicated in transcription, ribosome biogenesis, mRNP assembly, nuclear export, and long noncoding RNA (lncRNA) function. It is not understood how Dbp2p performs these seemingly unrelated biological roles. An important step toward addressing this question is the determination of cellular RNA binding sites of Dbp2p. Here, we identify transcriptome-wide RNA binding sites of Dbp2p from Saccharomyces cerevisiae using UV-crosslinking, denaturing tandem affinity purification, and next generation sequencing. We find that Dbp2p crosslinks to mRNAs and ribosomal RNAs, and markedly to noncoding RNAs, including snoRNA, snRNAs, and tRNAs. In snoRNAs, Dbp2p preferentially crosslinks at sites near the 3' ends. These sites coincide with regions where RNA-DNA hybrids (R-loops) form and with binding sites of Sen1p, another RNA helicase that functions in transcription termination and 3' processing of noncoding RNAs. We show that Dbp2p interacts in an RNA-independent manner with Sen1p in vivo. Dbp2p crosslinks to tRNAs and other RNAs also at sites where R-loops form. Collectively, our data link Dbp2p to noncoding RNAs, Sen1p, and R-loops. The transcriptome-wide connection to R-loops provides a unifying theme for diverse cellular roles of Dbp2p.

Knock-Down of a Novel snoRNA in Tetrahymena Reveals a Dual Role in 5.8S rRNA Processing and Generation of a 26S rRNA Fragment

In eukaryotes, 18S, 5.8S, and 28S rRNAs are transcribed as precursor molecules that undergo extensive modification and nucleolytic processing to form the mature rRNA species. Central in the process are the small nucleolar RNAs (snoRNAs). The majority of snoRNAs guide site specific chemical modifications but a few are involved in defining pre-rRNA cleavages. Here, we describe an unusual snoRNA (TtnuCD32) belonging to the box C/D subgroup from the ciliate Tetrahymena thermophila. We show that TtnuCD32 is unlikely to function as a modification guide snoRNA and that it is critical for cell viability. Cell lines with genetic knock-down of TtnuCD32 were impaired in growth and displayed two novel and apparently unrelated phenotypes. The most prominent phenotype is the accumulation of processing intermediates of 5.8S rRNA. The second phenotype is the decrease in abundance of a ~100 nt 26S rRNA fragment of unknown function. Sequence analysis demonstrated that TtnuCD32 share features with the essential snoRNA U14 but an alternative candidate (TtnuCD25) was more closely related to other U14 sequences. This, together with the fact that the observed rRNA processing phenotypes were not similar to what has been observed in U14 depleted cells, suggests that TtnuCD32 is a U14 homolog that has gained novel functions.

The 3' end formation in small RNAs

Small RNAs are a major class of RNAs along with transfer RNAs, ribosomal RNAs, and messenger RNAs. They vary in size from less than 100 nucleotides to several thousand nucleotides and have been identified and characterized both in prokaryotes and eukaryotes. Small RNAs participate in a variety of cellular functions including regulating RNA synthesis, RNA processing, guiding modifications in RNA, and in transport of proteins. Small RNAs are generated by a series of posttranscriptional processing steps following transcription. While RNA 5' end structure, 5' cap formation, and RNA processing mechanisms have been fairly well characterized, the 3' end processing is poorly understood. Recent data point to an emerging theme in small RNAs metabolism in which the 3' end processing is mediated by the exosome, a large multienzyme complex. In addition to removal of nucleotides by the exosome, there is simultaneous rebuilding of the 3' end of some small RNA by adenylation and/or uridylation. This review presents a picture of both degradative and rebuilding reactions operative on the 3' end of some small RNA molecules in prokaryotes and eukaryotes.

Small nucleolar RNAs: versatile trans-acting molecules of ancient evolutionary origin

The small nucleolar RNAs (snoRNAs) are an abundant class of trans-acting RNAs that function in ribosome biogenesis in the eukaryotic nucleolus. Elegant work has revealed that most known snoRNAs guide modification of pre-ribosomal RNA (pre-rRNA) by base pairing near target sites. Other snoRNAs are involved in cleavage of pre-rRNA by mechanisms that have not yet been detailed. Moreover, our appreciation of the cellular roles of the snoRNAs is expanding with new evidence that snoRNAs also target modification of small nuclear RNAs and messenger RNAs. Many snoRNAs are produced by unorthodox modes of biogenesis including salvage from introns of pre-mRNAs. The recent discovery that homologs of snoRNAs as well as associated proteins exist in the domain Archaea indicates that the RNA-guided RNA modification system is of ancient evolutionary origin. In addition, it has become clear that the RNA component of vertebrate telomerase (an enzyme implicated in cancer and cellular senescence) is related to snoRNAs. During its evolution, vertebrate telomerase RNA appears to have co-opted a snoRNA domain that is essential for the function of telomerase RNA in vivo. The unique properties of snoRNAs are now being harnessed for basic research and therapeutic applications.

Assembly and trafficking of box C/D and H/ACA snoRNPs

Box C/D and box H/ACA snoRNAs are abundant non-coding RNAs that localize in the nucleolus and mostly function as guides for nucleotide modifications. While a large pool of snoRNAs modifies rRNAs, an increasing number of snoRNAs could also potentially target mRNAs. ScaRNAs belong to a family of specific RNAs that localize in Cajal bodies and that are structurally similar to snoRNAs. Most scaRNAs are involved in snRNA modification, while telomerase RNA, which contains H/ACA motifs, functions in telomeric DNA synthesis. In this review, we describe how box C/D and H/ACA snoRNAs are processed and assembled with core proteins to form functional RNP particles. Their biogenesis involve several transport factors that first direct pre-snoRNPs to Cajal bodies, where some processing steps are believed to take place, and then to nucleoli. Assembly of core proteins involves the HSP90/R2TP chaperone-cochaperone system for both box C/D and H/ACA RNAs, but also several factors specific for each family. These assembly factors chaperone unassembled core proteins, regulate the formation and disassembly of pre-snoRNP intermediates, and control the activity of immature particles. The AAA+ ATPase RUVBL1 and RUVBL2 belong to the R2TP co-chaperones and play essential roles in snoRNP biogenesis, as well as in the formation of other macro-molecular complexes. Despite intensive research, their mechanisms of action are still incompletely understood.

Transcriptome-wide identification of NMD-targeted human mRNAs reveals extensive redundancy between SMG6- and SMG7-mediated degradation pathways

Besides degrading aberrant mRNAs that harbor a premature translation termination codon (PTC), nonsense-mediated mRNA decay (NMD) also targets many seemingly "normal" mRNAs that encode for full-length proteins. To identify a bona fide set of such endogenous NMD targets in human cells, we applied a meta-analysis approach in which we combined transcriptome profiling of knockdowns and rescues of the three NMD factors UPF1, SMG6, and SMG7. We provide evidence that this combinatorial approach identifies NMD-targeted transcripts more reliably than previous attempts that focused on inactivation of single NMD factors. Our data revealed that SMG6 and SMG7 act on essentially the same transcripts, indicating extensive redundancy between the endo- and exonucleolytic decay routes. Besides mRNAs, we also identified as NMD targets many long noncoding RNAs as well as miRNA and snoRNA host genes. The NMD target feature with the most predictive value is an intron in the 3' UTR, followed by the presence of upstream open reading frames (uORFs) and long 3' UTRs. Furthermore, the 3' UTRs of NMD-targeted transcripts tend to have an increased GC content and to be phylogenetically less conserved when compared to 3' UTRs of NMD insensitive transcripts.

The DEAD-box Protein Rok1 Orchestrates 40S and 60S Ribosome Assembly by Promoting the Release of Rrp5 from Pre-40S Ribosomes to Allow for 60S Maturation

DEAD-box proteins are ubiquitous regulators of RNA biology. While commonly dubbed “helicases,” their activities also include duplex annealing, adenosine triphosphate (ATP)-dependent RNA binding, and RNA-protein complex remodeling. Rok1, an essential DEAD-box protein, and its cofactor Rrp5 are required for ribosome assembly. Here, we use in vivo and in vitro biochemical analyses to demonstrate that ATP-bound Rok1, but not adenosine diphosphate (ADP)-bound Rok1, stabilizes Rrp5 binding to 40S ribosomes. Interconversion between these two forms by ATP hydrolysis is required for release of Rrp5 from pre-40S ribosomes in vivo, thereby allowing Rrp5 to carry out its role in 60S subunit assembly. Furthermore, our data also strongly suggest that the previously described accumulation of snR30 upon Rok1 inactivation arises because Rrp5 release is blocked and implicate a previously undescribed interaction between Rrp5 and the DEAD-box protein Has1 in mediating snR30 accumulation when Rrp5 release from pre-40S subunits is blocked.

During ribosomal biogenesis, Rrp5 is unusual in being required for assembly of both small and large subunits. This study demonstrates a role for ATP hydrolysis by the DEAD-box protein Rok1 in releasing Rrp5 from pre-40S subunits.

Assembly of the small and large ribosomal subunits requires two separate machineries. The assembly factor Rrp5 is unusual in being one of only three proteins required for assembly of both subunits. While it binds cotranscriptionally during early stages of small subunit assembly, it departs with large subunit intermediates after the separation of these precursors. How Rrp5 switches from interacting with small subunit precursors to binding large subunit precursors remains unknown but is potentially important, as it could regulate the interplay between small and large subunit assembly. Here, we show that the DEAD-box protein Rok1, a member of a ubiquitous class of RNA-dependent ATPases, releases Rrp5 from assembling small subunits to allow for its function in large subunit assembly. We show that a complex of Rrp5, Rok1, and adenosine triphosphate (ATP) binds small subunits or mimics of ribosomal RNA more tightly than does a complex of Rrp5, Rok1, and adenosine diphosphate (ADP). In cells, interconversion between the ATP and the ADP-form of Rok1 is required for release of Rrp5 from nascent small subunits and for binding to assembling large subunits. Furthermore, we show that the release of snR30, which leads to formation of a large substructure on small subunits, also requires Rok1-mediated release of Rrp5.

Network of epistatic interactions within a yeast snoRNA

Epistatic interactions play a fundamental role in molecular evolution, but little is known about the spatial distribution of these interactions within genes. To systematically survey a model landscape of intragenic epistasis, we quantified the fitness of ~60,000 Saccharomyces cerevisiae strains expressing randomly mutated variants of the 333-nucleotide-long U3 small nucleolar RNA (snoRNA). The fitness effects of individual mutations were correlated with evolutionary conservation and structural stability. Many mutations had small individual effects but had large effects in the context of additional mutations, which indicated negative epistasis. Clusters of negative interactions were explained by local thermodynamic threshold effects, whereas positive interactions were enriched among large-effect sites and between base-paired nucleotides. We conclude that high-throughput mapping of intragenic epistasis can identify key structural and functional features of macromolecules.

Enhanced snoMEN Vectors Facilitate Establishment of GFP-HIF-1α Protein Replacement Human Cell Lines

The snoMEN (snoRNA Modulator of gene ExpressioN) vector technology was developed from a human box C/D snoRNA, HBII-180C, which contains an internal sequence that can be manipulated to make it complementary to RNA targets, allowing knock-down of targeted genes. Here we have screened additional human nucleolar snoRNAs and assessed their application for gene specific knock-downs to improve the efficiency of snoMEN vectors. We identify and characterise a new snoMEN vector, termed 47snoMEN, that is derived from box C/D snoRNA U47, demonstrating its use for knock-down of both endogenous cellular proteins and G/YFP-fusion proteins. Using multiplex 47snoMEM vectors that co-express multiple 47snoMEN in a single transcript, each of which can target different sites in the same mRNA, we document >3-fold increase in knock-down efficiency when compared with the original HBII-180C based snoMEN. The multiplex 47snoMEM vector allowed the construction of human protein replacement cell lines with improved efficiency, including the establishment of novel GFP–HIF-1α replacement cells. Quantitative mass spectrometry analysis confirmed the enhanced efficiency and specificity of protein replacement using the 47snoMEN-PR vectors. The 47snoMEN vectors expand the potential applications for snoMEN technology in gene expression studies, target validation and gene therapy.

Revisiting the structure/function relationships of H/ACA(-like) RNAs: a unified model for Euryarchaea and Crenarchaea

A structural and functional classification of H/ACA and H/ACA-like motifs is obtained from the analysis of the H/ACA guide RNAs which have been identified previously in the genomes of Euryarchaea (Pyrococcus) and Crenarchaea (Pyrobaculum). A unified structure/function model is proposed based on the common structural determinants shared by H/ACA and H/ACA-like motifs in both Euryarchaea and Crenarchaea. Using a computational approach, structural and energetic rules for the guide:target RNA-RNA interactions are derived from structural and functional data on the H/ACA RNP particles. H/ACA(-like) motifs found in Pyrococcus are evaluated through the classification and their biological relevance is discussed. Extra-ribosomal targets found in both Pyrococcus and Pyrobaculum might support the hypothesis of a gene regulation mediated by H/ACA(-like) guide RNAs in archaea.

Identification of a new ribose methylation in the 18S rRNA of S. cerevisiae

Methylation of ribose sugars at the 2'-OH group is one of the major chemical modifications in rRNA, and is catalyzed by snoRNA directed C/D box snoRNPs. Previous biochemical and computational analyses of the C/D box snoRNAs have identified and mapped a large number of 2'-OH ribose methylations in rRNAs. In the present study, we systematically analyzed ribose methylations of 18S rRNA in Saccharomyces cerevisiae, using mung bean nuclease protection assay and RP-HPLC. Unexpectedly, we identified a hitherto unknown ribose methylation at position G562 in the helix 18 of 5' central domain of yeast 18S rRNA. Furthermore, we identified snR40 as being responsible to guide snoRNP complex to catalyze G562 ribose methylation, which makes it only second snoRNA known so far to target three ribose methylation sites: Gm562, Gm1271 in 18S rRNA, and Um898 in 25S rRNA. Our sequence and mutational analysis of snR40 revealed that snR40 uses the same D' box and methylation guide sequence for both Gm562 and Gm1271 methylation. With the identification of Gm562 and its corresponding snoRNA, complete set of ribose methylations of 18S rRNA and their corresponding snoRNAs have finally been established opening great prospects to understand the physiological function of these modifications.

The human nuclear exosome targeting complex is loaded onto newly synthesized RNA to direct early ribonucleolysis

The RNA exosome complex constitutes the major nuclear eukaryotic 3'-5' exonuclease. Outside of nucleoli, the human nucleoplasmic exosome is directed to some of its substrates by the nuclear exosome targeting (NEXT) complex. How NEXT targets RNA has remained elusive. Using an in vivo crosslinking approach, we report global RNA binding sites of RBM7, a key component of NEXT. RBM7 associates broadly with RNA polymerase II-derived RNA, including pre-mRNA and short-lived exosome substrates such as promoter upstream transcripts (PROMPTs), enhancer RNAs (eRNAs), and 3'-extended products from snRNA and replication-dependent histone genes. Within pre-mRNA, RBM7 accumulates at the 3' ends of introns, and pulse-labeling experiments demonstrate that RBM7/NEXT defines an early exosome-targeting pathway for 3'-extended snoRNAs derived from such introns. We propose that RBM7 is generally loaded onto newly synthesized RNA to accommodate exosome action in case of available unprotected RNA 3' ends.

A pre-ribosomal RNA interaction network involving snoRNAs and the Rok1 helicase

Ribosome biogenesis in yeast requires 75 small nucleolar RNAs (snoRNAs) and a myriad of cofactors for processing, modification, and folding of the ribosomal RNAs (rRNAs). For the 19 RNA helicases implicated in ribosome synthesis, their sites of action and molecular functions have largely remained unknown. Here, we have used UV cross-linking and analysis of cDNA (CRAC) to reveal the pre-rRNA binding sites of the RNA helicase Rok1, which is involved in early small subunit biogenesis. Several contact sites were identified in the 18S rRNA sequence, which interestingly all cluster in the "foot" region of the small ribosomal subunit. These include a major binding site in the eukaryotic expansion segment ES6, where Rok1 is required for release of the snR30 snoRNA. Rok1 directly contacts snR30 and other snoRNAs required for pre-rRNA processing. Using cross-linking, ligation and sequencing of hybrids (CLASH) we identified several novel pre-rRNA base-pairing sites for the snoRNAs snR30, snR10, U3, and U14, which cluster in the expansion segments of the 18S rRNA. Our data suggest that these snoRNAs bridge interactions between the expansion segments, thereby forming an extensive interaction network that likely promotes pre-rRNA maturation and folding in early pre-ribosomal complexes and establishes long-range rRNA interactions during ribosome synthesis.

Targeting vertebrate intron-encoded box C/D 2'-O-methylation guide RNAs into the Cajal body

Post-transcriptional pseudouridylation and 2'-O-methylation of splicesomal small nuclear ribonucleic acids (snRNAs) is mediated by box H/ACA and box C/D small Cajal body (CB)-specific ribonucleoproteins (scaRNPs), respectively. The WD-repeat protein 79 (WDR79) has been proposed to interact with both classes of modification scaRNPs and target them into the CB. The box H/ACA scaRNAs carry the common CAB box motif (consensus, ugAG) that is required for both WDR79 binding and CB-specific accumulation. Thus far, no cis-acting CB-localization element has been reported for vertebrate box C/D scaRNAs. In this study, systematic mutational analysis of the human U90 and another newly identified box C/D scaRNA, mgU2-47, demonstrated that the CB-specific accumulation of vertebrate intron-encoded box C/D scaRNAs relies on GU- or UG-dominated dinucleotide repeat sequences which are predicted to form the terminal stem-loop of the RNA apical hairpin. While the loop nucleotides are unimportant, the adjacent terminal helix that is composed mostly of consecutive G.U and U.G wobble base-pairs is essential for CB-specific localization of box C/D scaRNAs. Co-immunoprecipitation experiments confirmed that the newly identified CB localization element, called the G.U/U.G wobble stem, is crucial for in vivo association of box C/D scaRNPs with WDR79.

Imp3 unfolds stem structures in pre-rRNA and U3 snoRNA to form a duplex essential for small subunit processing

Eukaryotic ribosome biogenesis requires rapid hybridization between the U3 snoRNA and the pre-rRNA to direct cleavages at the A0, A1, and A2 sites in pre-rRNA that liberate the small subunit precursor. The bases involved in hybridization of one of the three duplexes that U3 makes with pre-rRNA, designated the U3-18S duplex, are buried in conserved structures: box A/A' stem-loop in U3 snoRNA and helix 1 (H1) in the 18S region of the pre-rRNA. These conserved structures must be unfolded to permit the necessary hybridization. Previously, we reported that Imp3 and Imp4 promote U3-18S hybridization in vitro, but the mechanism by which these proteins facilitate U3-18S duplex formation remained unclear. Here, we directly addressed this question by probing base accessibility with chemical modification and backbone accessibility with ribonuclease activity of U3 and pre-rRNA fragments that mimic the secondary structure observed in vivo. Our results demonstrate that U3-18S hybridization requires only Imp3. Binding to each RNA by Imp3 provides sufficient energy to unfold both the 18S H1 and the U3 box A/A' stem structures. The Imp3 unfolding activity also increases accessibility at the U3-dependent A0 and A1 sites, perhaps signaling cleavage at these sites to generate the 5' mature end of 18S. Imp4 destabilizes the U3-18S duplex to aid U3 release, thus differentiating the roles of these proteins. Protein-dependent unfolding of these structures may serve as a switch to block U3-pre-rRNA interactions until recruitment of Imp3, thereby preventing premature and inaccurate U3-dependent pre-rRNA cleavage and folding events in eukaryotic ribosome biogenesis.

Generation of a mouse model with down-regulated U50 snoRNA (SNORD50) expression and its organ-specific phenotypic modulation

Box C/D-type small nucleolar RNAs (snoRNAs) are functional RNAs responsible for mediating 2'-O-ribose methylation of ribosomal RNAs (rRNAs) within the nucleolus. In the past years, evidence for the involvement of human U50 snoRNA in tumorigenesis has been accumulating. We previously identified U50HG, a non-protein-coding gene that hosted a box C/D-type U50 snoRNA, in a chromosomal breakpoint in a human B-cell lymphoma. Mouse genome analysis revealed four mouse U50 (mU50) host-genes: three mU50HG-a gene variants that were clustered in the genome and an mU50HG-b gene that we supposed to be the U50HG ortholog. In this study, to investigate the physiological importance of mU50 snoRNA and its involvement in tumorigenesis, we eliminated mU50 snoRNA sequences from the mU50HG-b gene. The established mouse line (ΔmU50(HG-b)) showed a significant reduction of mU50 snoRNA expression without alteration of the host-gene length and exon-intron structure, and the corresponding target rRNA methylation in various organs was reduced. Lifelong phenotypic monitoring showed that the ΔmU50(HG-b) mice looked almost normal without accelerated tumorigenicity; however, a notable difference was the propensity for anomalies in the lymphoid organs. Transcriptome analysis showed that dozens of genes, including heat shock proteins, were differentially expressed in ΔmU50(HG-b) mouse lymphocytes. This unique model of a single snoRNA knockdown with intact host-gene expression revealed further new insights into the discrete transcriptional regulation of multiple mU50 host-genes and the complicated dynamics involved in organ-specific processing and maintenance of snoRNAs.

Comparative study of two box H/ACA ribonucleoprotein pseudouridine-synthases: relation between conformational dynamics of the guide RNA, enzyme assembly and activity

Multiple RNA-guided pseudouridine synthases, H/ACA ribonucleoprotein particles (RNPs) which contain a guide RNA and four proteins, catalyze site-specific post-transcriptional isomerization of uridines into pseudouridines in substrate RNAs. In archaeal particles, the guide small RNA (sRNA) is anchored by the pseudouridine synthase aCBF5 and the ribosomal protein L7Ae. Protein aNOP10 interacts with both aCBF5 and L7Ae. The fourth protein, aGAR1, interacts with aCBF5 and enhances catalytic efficiency. Here, we compared the features of two H/ACA sRNAs, Pab21 and Pab91, from Pyrococcus abyssi. We found that aCBF5 binds much more weakly to Pab91 than to Pab21. Surprisingly, the Pab91 sRNP exhibits a higher catalytic efficiency than the Pab21 sRNP. We thus investigated the molecular basis of the differential efficiencies observed for the assembly and catalytic activity of the two enzymes. For this, we compared profiles of the extent of lead-induced cleavages in these sRNAs during a stepwise reconstitution of the sRNPs, and analyzed the impact of the absence of the aNOP10–L7Ae interaction. Such probing experiments indicated that the sRNAs undergo a series of conformational changes upon RNP assembly. These changes were also evaluated directly by circular dichroism (CD) spectroscopy, a tool highly adapted to analyzing RNA conformational dynamics. In addition, our results reveal that the conformation of helix P1 formed at the base of the H/ACA sRNAs is optimized in Pab21 for efficient aCBF5 binding and RNP assembly. Moreover, P1 swapping improved the assembly of the Pab91 sRNP. Nonetheless, efficient aCBF5 binding probably also relies on the pseudouridylation pocket which is not optimized for high activity in the case of Pab21.

Structural and functional analysis of the U3 snoRNA binding protein Rrp9

The U3 snoRNA is required for 18S rRNA processing and small subunit ribosome formation in eukaryotes. Different from other box C/D snoRNAs, U3 contains an extra 5' domain that pairs with pre-rRNA and a unique B/C motif essential for recruitment of the U3-specific Rrp9 protein. Here, we analyze the structure and function of Rrp9 with crystallographic, biochemical, and cellular approaches. Rrp9 is composed of a WD repeat domain and an N-terminal region. The crystal structures of the WD domain of yeast Rrp9 and its human ortholog U3-55K were determined, revealing a typical seven-bladed propeller fold. Several conserved surface patches on the WD domain were identified, and their function in RNP assembly and yeast growth were analyzed by mutagenesis. Prior association of Snu13 with the B/C motif was found to enhance the specific binding of the WD domain. We show that a conserved 7bc loop is crucial for specific recognition of U3, nucleolar localization of Rrp9, and yeast growth. In addition, we show that the N-terminal region of Rrp9 contains a bipartite nuclear localization signal that is dispensable for nucleolar localization. Our results provide insight into the functional sites of Rrp9.

Kinetic analysis of pre-ribosome structure in vivo

Pre-ribosomal particles undergo numerous structural changes during maturation, but their high complexity and short lifetimes make these changes very difficult to follow in vivo. In consequence, pre-ribosome structure and composition have largely been inferred from purified particles and analyzed in vitro. Here we describe techniques for kinetic analyses of the changes in pre-ribosome structure in living cells of Saccharomyces cerevisiae. To allow this, in vivo structure probing by DMS modification was combined with affinity purification of newly synthesized 20S pre-rRNA over a time course of metabolic labeling with 4-thiouracil. To demonstrate that this approach is generally applicable, we initially analyzed the accessibility of the region surrounding cleavage site D site at the 3' end of the mature 18S rRNA region of the pre-rRNA. This revealed a remarkably flexible structure throughout 40S subunit biogenesis, with little stable RNA-protein interaction apparent. Analysis of folding in the region of the 18S central pseudoknot was consistent with previous data showing U3 snoRNA-18S rRNA interactions. Dynamic changes in the structure of the hinge between helix 28 (H28) and H44 of pre-18S rRNA were consistent with recently reported interactions with the 3' guide region of U3 snoRNA. Finally, analysis of the H18 region indicates that the RNA structure matures early, but additional protection appears subsequently, presumably reflecting protein binding. The structural analyses described here were performed on total, affinity-purified, newly synthesized RNA, so many classes of RNA and RNA-protein complex are potentially amenable to this approach.

Kinetic and thermodynamic characterization of the reaction pathway of box H/ACA RNA-guided pseudouridine formation

The box H/ACA RNA-guided pseudouridine synthase is a complicated ribonucleoprotein enzyme that recruits substrate via both the guide RNA and the catalytic subunit Cbf5. Structural studies have revealed multiple conformations of the enzyme, but a quantitative description of the reaction pathway is still lacking. Using fluorescence correlation spectroscopy, we here measured the equilibrium dissociation constants and kinetic association and dissociation rates of substrate and product complexes mimicking various reaction intermediate states. These data support a sequential model for substrate loading and product release regulated by the thumb loop of Cbf5. The uridine substrate is first bound primarily through interaction with the guide RNA and then loaded into the active site while progressively interacted with the thumb. After modification, the subtle chemical structure change from uridine to pseudouridine at the target site triggers the release of the thumb, resulting in an intermediate complex with the product bound mainly by the guide RNA. By dissecting the role of Gar1 in individual steps of substrate turnover, we show that Gar1 plays a major role in catalysis and also accelerates product release about 2-fold. Our biophysical results integrate with previous structural knowledge into a coherent reaction pathway of H/ACA RNA-guided pseudouridylation.

Box H/ACA snoRNAs are preferred substrates for the trimethylguanosine synthase in the divergent unicellular eukaryote Trichomonas vaginalis

The 2,2,7-trimethylguanosine caps of eukaryal snRNAs and snoRNA are formed by the enzyme Tgs1, which catalyzes sequential guanine-N2 methylations of m(7)G caps. Atypically, in the divergent unicellular eukaryote Trichomonas vaginalis, spliceosomal snRNAs lack a guanosine cap and the recombinant T. vaginalis trimethylguanosine synthase (TvTgs) produces only m(2,7)G in vitro. Here, we show by direct metabolic labeling that endogenous T. vaginalis RNAs contain m(7)G, m(2,7)G, and m(2,2,7)G caps. Immunodepletion of TvTgs from cell extracts and TvTgs add-back experiments demonstrate that TvTgs produces m(2,7)G and m(2,2,7)G caps. Expression of TvTgs in yeast tgs1Δ cells leads to the formation of m(2,7)G and m(2,2,7)G caps and complementation of the lethality of a tgs1Δ mud2Δ strain. Whereas TvTgs is present in the nucleus and cytosol of T. vaginalis cells, TMG-containing RNAs are localized primarily in the nucleolus. Molecular cloning of anti-TMG affinity-purified T. vaginalis RNAs identified 16 box H/ACA snoRNAs, which are implicated in guiding RNA pseudouridylation. The ensemble of new T. vaginalis H/ACA snoRNAs allowed us to predict and partially validate an extensive map of pseudouridines in T. vaginalis rRNA.

Structural and functional characterization of ribosomal protein gene introns in sponges

Ribosomal protein genes (RPGs) are a powerful tool for studying intron evolution. They exist in all three domains of life and are much conserved. Accumulating genomic data suggest that RPG introns in many organisms abound with non-protein-coding-RNAs (ncRNAs). These ancient ncRNAs are small nucleolar RNAs (snoRNAs) essential for ribosome assembly. They are also mobile genetic elements and therefore probably important in diversification and enrichment of transcriptomes through various mechanisms such as intron/exon gain/loss. snoRNAs in basal metazoans are poorly characterized. We examined 449 RPG introns, in total, from four demosponges: Amphimedon queenslandica, Suberites domuncula, Suberites ficus and Suberites pagurorum and showed that RPG introns from A. queenslandica share position conservancy and some structural similarity with "higher" metazoans. Moreover, our study indicates that mobile element insertions play an important role in the evolution of their size. In four sponges 51 snoRNAs were identified. The analysis showed discrepancies between the snoRNA pools of orthologous RPG introns between S. domuncula and A. queenslandica. Furthermore, these two sponges show as much conservancy of RPG intron positions between each other as between themselves and human. Sponges from the Suberites genus show consistency in RPG intron position conservation. However, significant differences in some of the orthologous RPG introns of closely related sponges were observed. This indicates that RPG introns are dynamic even on these shorter evolutionary time scales.

Extensive terminal and asymmetric processing of small RNAs from rRNAs, snoRNAs, snRNAs, and tRNAs

Deep sequencing studies frequently identify small RNA fragments of abundant RNAs. These fragments are thought to represent degradation products of their precursors. Using sequencing, computational analysis, and sensitive northern blot assays, we show that constitutively expressed non-coding RNAs such as tRNAs, snoRNAs, rRNAs and snRNAs preferentially produce small 5' and 3' end fragments. Similar to that of microRNA processing, these terminal fragments are generated in an asymmetric manner that predominantly favors either the 5' or 3' end. Terminal-specific and asymmetric processing of these small RNAs occurs in both mouse and human cells. In addition to the known processing of some 3' terminal tRNA-derived fragments (tRFs) by the RNase III endonuclease Dicer, we show that several RNase family members can produce tRFs, including Angiogenin that cleaves the TψC loop to generate 3' tRFs. The 3' terminal tRFs but not the 5' tRFs are highly complementary to human endogenous retroviral sequences in the genome. Despite their independence from Dicer processing, these tRFs associate with Ago2 and are capable of down regulating target genes by transcript cleavage in vitro. We suggest that endogenous 3' tRFs have a role in regulating the unwarranted expression of endogenous viruses through the RNA interference pathway.

Maturation of mammalian H/ACA box snoRNAs: PAPD5-dependent adenylation and PARN-dependent trimming

Small nucleolar and small Cajal body RNAs (snoRNAs and scaRNAs) of the H/ACA box and C/D box type are generated by exonucleolytic shortening of longer precursors. Removal of the last few nucleotides at the 3' end is known to be a distinct step. We report that, in human cells, knock-down of the poly(A) specific ribonuclease (PARN), previously implicated only in mRNA metabolism, causes the accumulation of oligoadenylated processing intermediates of H/ACA box but not C/D box RNAs. In agreement with a role of PARN in snoRNA and scaRNA processing, the enzyme is concentrated in nucleoli and Cajal bodies. Oligo(A) tails are attached to a short stub of intron sequence remaining beyond the mature 3' end of the snoRNAs. The noncanonical poly(A) polymerase PAPD5 is responsible for addition of the oligo(A) tails. We suggest that deadenylation is coupled to clean 3' end trimming, which might serve to enhance snoRNA stability.

Features of a unique intronless cluster of class I small heat shock protein genes in tandem with box C/D snoRNA genes on chromosome 6 in tomato (Solanum lycopersicum)

Physical clustering of genes has been shown in plants; however, little is known about gene clusters that have different functions, particularly those expressed in the tomato fruit. A class I 17.6 small heat shock protein (Sl17.6 shsp) gene was cloned and used as a probe to screen a tomato (Solanum lycopersicum) genomic library. An 8.3-kb genomic fragment was isolated and its DNA sequence determined. Analysis of the genomic fragment identified intronless open reading frames of three class I shsp genes (Sl17.6, Sl20.0, and Sl20.1), the Sl17.6 gene flanked by Sl20.1 and Sl20.0, with complete 5' and 3' UTRs. Upstream of the Sl20.0 shsp, and within the shsp gene cluster, resides a box C/D snoRNA cluster made of SlsnoR12.1 and SlU24a. Characteristic C and D, and C' and D', boxes are conserved in SlsnoR12.1 and SlU24a while the upstream flanking region of SlsnoR12.1 carries TATA box 1, homol-E and homol-D box-like cis sequences, TM6 promoter, and an uncharacterized tomato EST. Molecular phylogenetic analysis revealed that this particular arrangement of shsps is conserved in tomato genome but is distinct from other species. The intronless genomic sequence is decorated with cis elements previously shown to be responsive to cues from plant hormones, dehydration, cold, heat, and MYC/MYB and WRKY71 transcription factors. Chromosomal mapping localized the tomato genomic sequence on the short arm of chromosome 6 in the introgression line (IL) 6-3. Quantitative polymerase chain reaction analysis of gene cluster members revealed differential expression during ripening of tomato fruit, and relatively different abundances in other plant parts.

hNUDT16: a universal decapping enzyme for small nucleolar RNA and cytoplasmic mRNA

Human NUDT16 (hNUDT16) is a decapping enzyme initially identified as the human homolog to the Xenopus laevis X29. As a metalloenzyme, hNUDT16 relies on divalent cations for its cap-hydrolysis activity to remove m⁷GDP and m²²⁷GDP from RNAs. Metal also determines substrate specificity of the enzyme. So far, only U8 small nucleolar RNA (snoRNA) has been identified as the substrate of hNUDT16 in the presence of Mg²(+). Here we demonstrate that besides U8, hNUDT16 can also actively cleave the m⁷GDP cap from mRNAs in the presence of Mg²(+) or Mn²(+). We further show that hNUDT16 does not preferentially recognize U8 or mRNA substrates by our cross-inhibition and quantitative decapping assays. In addition, our mutagenesis analysis identifies several key residues involved in hydrolysis and confirms the key role of the REXXEE motif in catalysis. Finally an investigation into the subcellular localization of hNUDT16 revealed its abundance in both cytoplasm and nucleus. These findings extend the substrate spectrum of hNUDT16 beyond snoRNAs to also include mRNA, demonstrating the pleiotropic decapping activity of hNUDT16.

Structural and functional evidence of high specificity of Cbf5 for ACA trinucleotide

Cbf5 is the catalytic subunit of the H/ACA small nucleolar/Cajal body ribonucleoprotein particles (RNPs) responsible for site specific isomerization of uridine in ribosomal and small nuclear RNA. Recent evidence from studies on archaeal Cbf5 suggests its second functional role in modifying tRNA U55 independent of guide RNA. In order to act both as a stand-alone and a RNP pseudouridine synthase, Cbf5 must differentiate features in H/ACA RNA from those in tRNA or rRNA. Most H/ACA RNAs contain a hallmark ACA trinucleotide downstream of the H/ACA motif. Here we challenged an archaeal Cbf5 (in the form of a ternary complex with its accessory proteins Nop10 and Gar1) with T-stem-loop RNAs with or without ACA trinucleotide in the stem. Although these substrates were previously shown to be substrates for the bacterial stand-alone pseudouridine synthase TruB, the Cbf5-Nop10-Gar1 complex was only able to modify those without ACA trinucleotide. A crystal structure of Cbf5-Nop10-Gar1 trimer bound with an ACA-containing T-stem-loop revealed that the ACA trinucleotide detracted Cbf5 from the stand-alone binding mode, thereby suggesting that the H/ACA RNP-associated function of Cbf5 likely supersedes its stand-alone function.

Efficient and specific knockdown of small non-coding RNAs in mammalian cells and in mice

Hundreds of small nuclear non-coding RNAs, including small nucleolar RNAs (snoRNAs), have been identified in different organisms, with important implications in regulating gene expression and in human diseases. However, functionalizing these nuclear RNAs in mammalian cells remains challenging, due to methodological difficulties in depleting these RNAs, especially snoRNAs. Here we report a convenient and efficient approach to deplete snoRNA, small Cajal body RNA (scaRNA) and small nuclear RNA in human and mouse cells by conventional transfection of chemically modified antisense oligonucleotides (ASOs) that promote RNaseH-mediated cleavage of target RNAs. The levels of all seven tested snoRNA/scaRNAs and four snRNAs were reduced by 80-95%, accompanied by impaired endogenous functions of the target RNAs. ASO-targeting is highly specific, without affecting expression of the host genes where snoRNAs are embedded in the introns, nor affecting the levels of snoRNA isoforms with high sequence similarities. At least five snoRNAs could be depleted simultaneously. Importantly, snoRNAs could be dramatically depleted in mice by systematic administration of the ASOs. Together, our findings provide a convenient and efficient approach to characterize nuclear non-coding RNAs in mammalian cells, and to develop antisense drugs against disease-causing non-coding RNAs.

A novel Nop5-sRNA interaction that is required for efficient archaeal box C/D sRNP formation

Archaeal and eukaryotic box C/D RNPs catalyze the 2'-O-methylation of ribosomal RNA, a modification that is essential for the correct folding and function of the ribosome. Each archaeal RNP contains three core proteins--L7Ae, Nop5, and fibrillarin (methyltransferase)--and a box C/D sRNA. Base-pairing between the sRNA guide region and the rRNA directs target site selection with the C/D and related C'/D' motifs functioning as protein binding sites. Recent structural analysis of in vitro assembled archaeal complexes has produced two divergent models of box C/D sRNP structure. In one model, the complex is proposed to be monomeric, while the other suggests a dimeric sRNP. The position of the RNA in the RNP is significantly different in each model. We have used UV-cross-linking to characterize protein-RNA contacts in the in vitro assembled Pyrococcus furiosus box C/D sRNP. The P. furiosus sRNP components assemble into complexes that are the expected size of di-sRNPs. Analysis of UV-induced protein-RNA cross-links revealed a novel interaction between the ALFR motif, in the Nop domain of Nop5, and the guide/spacer regions of the sRNA. We show that the ALFR motif and the spacer sequence adjacent to box C or C' are important for box C/D sRNP assembly in vitro. These data therefore reveal new RNA-protein contacts in the box C/D sRNP and suggest a role for Nop5 in substrate binding and/or release.

Structure and RNA recognition by the snRNA and snoRNA transport factor PHAX

Small nuclear and small nucleolar RNAs (snRNAs and snoRNAs) are critical components of snRNPs and snoRNPs and play an essential role in the maturation of, respectively, mRNAs and rRNAs within the nucleus of eukaryotic cells. Complex and specific pathways exist for the assembly of snRNPs and snoRNPs, involving, for instance, nucleocytoplasmic transport of snRNAs and intranuclear transport between compartments of snoRNAs. The phosphorylated adaptor for nuclear export (PHAX) is required for nuclear export of snRNAs in metazoans and also involved in the intranuclear transport of snoRNAs to Cajal bodies. PHAX contains a conserved single-stranded nucleic acid binding domain (RNA_GG_bind domain) with no sequence homology with any other known RNA-binding module. Here, we report NMR and X-ray crystallography studies that elucidate the structural basis for RNA recognition by the PHAX RNA-binding domain (PHAX-RBD). The crystal structure of the RNA_GG_bind domain from the parasite Cryptosporidium parvum (Cp RBD) forms well-folded dimers in solution in the absence of any ligand. The human PHAX-RBD is monomeric and only adopts a tertiary fold upon RNA binding. The PHAX-RBD represents a novel helical fold and binds single-stranded RNA with micromolar affinity without sequence specificity. RNA recognition by human PHAX-RBD is consistent with mutational analysis that affects RNA binding and PHAX-mediated nuclear export. Our data suggest that the PHAX-RBD mediates auxiliary RNA contacts with the snRNA and snoRNA substrates that are required for transport and/or substrate release.

Regulation of pre-mRNA splicing in Xenopus oocytes by targeted 2'-O-methylation

The 2'-OH group of the branch point adenosine is a key moiety to initiate pre-mRNA splicing. We use RNA-guided RNA modification to target the pre-mRNA branch point adenosine for 2'-O-methylation, with the aim of blocking pre-mRNA splicing in vertebrate cells. We show that, under certain conditions, injection of a branch point-specific artificial box C/D RNA into Xenopus oocytes effectively 2'-O-methylates adenovirus pre-mRNA at the target nucleotide. However, 2'-O-methylation at the authentic branch point activates a host of cryptic branch points, thus allowing splicing to continue. These cryptic sites are mapped, and mutated. Upon injection, pre-mRNA free of cryptic branch points fails to splice when the branch point-specific box C/D RNA is present. However, 2'-O-methylation at the branch point does not prevent pre-mRNA from being assembled into pre-catalytic spliceosome-like complexes prior to the first chemical step of splicing. Our results demonstrate that RNA-guided pre-mRNA modification can occur in the nucleoplasm of vertebrate cells, thus offering a powerful tool for molecular biology research.

Plant U13 orthologues and orphan snoRNAs identified by RNomics of RNA from Arabidopsis nucleoli

Small nucleolar RNAs (snoRNAs) and small Cajal body-specific RNAs (scaRNAs) are non-coding RNAs whose main function in eukaryotes is to guide the modification of nucleotides in ribosomal and spliceosomal small nuclear RNAs, respectively. Full-length sequences of Arabidopsis snoRNAs and scaRNAs have been obtained from cDNA libraries of capped and uncapped small RNAs using RNA from isolated nucleoli from Arabidopsis cell cultures. We have identified 31 novel snoRNA genes (9 box C/D and 22 box H/ACA) and 15 new variants of previously described snoRNAs. Three related capped snoRNAs with a distinct gene organization and structure were identified as orthologues of animal U13snoRNAs. In addition, eight of the novel genes had no complementarity to rRNAs or snRNAs and are therefore putative orphan snoRNAs potentially reflecting wider functions for these RNAs. The nucleolar localization of a number of the snoRNAs and the localization to nuclear bodies of two putative scaRNAs was confirmed by in situ hybridization. The majority of the novel snoRNA genes were found in new gene clusters or as part of previously described clusters. These results expand the repertoire of Arabidopsis snoRNAs to 188 snoRNA genes with 294 gene variants.

The small nucleolar ribonucleoprotein (snoRNP) database

Small nucleolar ribonucleoproteins (snoRNPs) are widely studied and characterized as guide RNAs for sequence-specific 2'-O-ribose methylation and psuedouridylation of ribosomal RNAs. In addition, snoRNAs have also been shown to interact with some tRNAs and direct alternative splicing in mRNA biogenesis. Recent advances in bioinformatics have resulted in new algorithms able to rapidly identify noncoding RNAs generally and snoRNAs specifically in genomic and metagenomic sequences, resulting in a rapid increase in the number and diversity of identified snoRNA sequences. The snoRNP database is a web-based collection of snoRNA and snoRNA-associated protein sequences from a wide range of species. The database currently contains 8994 snoRNA sequences from Bacteria, Archaea, and Eukaryotes and 589 snoRNA-associated protein sequences. The snoRNP database can be found at: http://evolveathome.com/snoRNA/snoRNA.php.

Identification of protein factors and U3 snoRNAs from a Brassica oleracea RNP complex involved in the processing of pre-rRNA

We report on the structural characterization of a functional U3 snoRNA ribonucleoprotein complex isolated from Brassica oleracea. The BoU3 snoRNP complex (formerly NF D) binds ribosomal DNA (rDNA), specifically cleaves pre-rRNA at the primary cleavage site in vitro and probably links transcription to early pre-rRNA processing in vivo. Using a proteomic approach we have identified 62 proteins in the purified BoU3 snoRNP fraction, including small RNA associated proteins (Fibrillarin, NOP5/Nop58p, Diskerin/Cbf5p, SUS2/PRP8 and CLO/GFA1/sn114p) and 40S ribosomal associated proteins (22 RPS and four ARCA-like proteins). Another major protein group is composed of chaperones/chaperonins (HSP81/TCP-1) and at least one proteasome subunit (RPN1a). Remarkably, RNA-dependent RNA polymerase (RdRP) and Tudor staphylococcal nuclease (TSN) proteins, which have RNA- and/or DNA-associated activities, were also revealed in the complex. Furthermore, three U3 snoRNA variants were identified in the BoU3 snoRNP fraction, notably an evolutionarily conserved and variable stem loop structure located just downstream from the C-box domain of the U3 sequence structures. We conclude that the BoU3 snoRNP complex is mainly required for 40S pre-ribosome synthesis. It is also expected that U3 snoRNA variants and interacting proteins might play a major role in BoU3 snoRNP complex assembly and/or function. This study provides a basis for further investigation of these novel ribonucleoprotein factors and their role in plant ribosome biogenesis.

Identification and classification of small RNAs in transcriptome sequence data

Current methods for high throughput sequencing (HTS) for the first time offer the opportunity to investigate the entire transcriptome in an essentially unbiased way. In many species, small non-coding RNAs with specific secondary structures constitute a significant part of the transcriptome. Some of these RNA classes, in particular microRNAs and snoRNAs, undergo maturation processes that lead to the production of shorter RNAs. After mapping the sequences to the reference genome specific patterns of short reads can be observed. These read patterns seem to reflect the processing and thus are specific for the RNA transcripts of which they are derived from. We explore here the potential of short read sequence data in the classification and identification of non-coding RNAs.

Human miRNA precursors with box H/ACA snoRNA features

MicroRNAs (miRNAs) and small nucleolar RNAs (snoRNAs) are two classes of small non-coding regulatory RNAs, which have been much investigated in recent years. While their respective functions in the cell are distinct, they share interesting genomic similarities, and recent sequencing projects have identified processed forms of snoRNAs that resemble miRNAs. Here, we investigate a possible evolutionary relationship between miRNAs and box H/ACA snoRNAs. A comparison of the genomic locations of reported miRNAs and snoRNAs reveals an overlap of specific members of these classes. To test the hypothesis that some miRNAs might have evolved from snoRNA encoding genomic regions, reported miRNA-encoding regions were scanned for the presence of box H/ACA snoRNA features. Twenty miRNA precursors show significant similarity to H/ACA snoRNAs as predicted by snoGPS. These include molecules predicted to target known ribosomal RNA pseudouridylation sites in vivo for which no guide snoRNA has yet been reported. The predicted folded structures of these twenty H/ACA snoRNA-like miRNA precursors reveal molecules which resemble the structures of known box H/ACA snoRNAs. The genomic regions surrounding these predicted snoRNA-like miRNAs are often similar to regions around snoRNA retroposons, including the presence of transposable elements, target site duplications and poly (A) tails. We further show that the precursors of five H/ACA snoRNA-like miRNAs (miR-151, miR-605, mir-664, miR-215 and miR-140) bind to dyskerin, a specific protein component of functional box H/ACA small nucleolar ribonucleoprotein complexes suggesting that these molecules have retained some H/ACA snoRNA functionality. The detection of small RNA molecules that share features of miRNAs and snoRNAs suggest that these classes of RNA may have an evolutionary relationship.

The major functions known for RNA were long believed to be either messenger RNAs, which function as intermediates between genes and proteins, or ribosomal RNAs and transfer RNAs which carry out the translation process. In recent years, however, newly discovered classes of small RNAs have been shown to play important cellular roles. These include microRNAs (miRNAs), which can regulate the production of specific proteins, and small nucleolar RNAs (snoRNAs), which recognise and chemically modify specific sequences in ribosomal RNA. Although miRNAs and snoRNAs are currently believed to be generated by different cellular pathways and to function in different cellular compartments, members of these two types of small RNAs display numerous genomic similarities, and a small number of snoRNAs have been shown to encode miRNAs in several organisms. Here we systematically investigate a possible evolutionary relationship between snoRNAs and miRNAs. Using computational analysis, we identify twenty genomic regions encoding miRNAs with highly significant similarity to snoRNAs, both on the level of their surrounding genomic context as well as their predicted folded structure. A subset of these miRNAs display functional snoRNA characteristics, strengthening the possibility that these miRNA molecules might have evolved from snoRNAs.

Processing of a dicistronic tRNA-snoRNA precursor: combined analysis in vitro and in vivo reveals alternate pathways and coupling to assembly of snoRNP

The C/D box small nucleolar RNAs (snoRNAs) represent an essential class of small nucleolar RNAs that guide 2'-O-Rib methylation of ribosomal RNAs and other RNAs in eukaryotes. In Arabidopsis (Arabidopsis thaliana), >100 C/D snoRNAs have been identified, most of them encoded by polycistronic gene clusters, but little is known on the factors controlling their biogenesis. Here, we focus on the identification of factors controlling the processing of tRNA-snoRNA dicistronic precursors (pre-tsnoRNA) synthesized by RNA polymerase III and producing tRNA(Gly) and C/D snoR43. We produced radiolabeled RNA probes corresponding to different pre-tsnoRNA mutants to test their impact on processing in vitro by a recombinant tRNAse Z, the Arabidopsis endonuclease that processes the 3'end of tRNAs, and by nuclear extracts from cauliflower (Brassica oleracea) inflorescences that accurately process the pre-tsnoRNA. This was coupled to an in vivo analysis of the processing of tagged pre-tsnoRNA mutants expressed in Arabidopsis. Our results strongly implicate tRNase Z in endonucleolytic cleavage of the pre-tsnoRNA. In addition, they reveal an alternate pathway that could depend on a tRNA decay surveillance mechanism. Finally, we provide arguments showing that processing of pre-tsnoRNA, both in planta and by nuclear extracts, is coupled to the assembly of snoRNA with core proteins forming the functional snoRNP (for small nucleolar ribonucleoprotein complex).

RNA chaperones stimulate formation and yield of the U3 snoRNA-Pre-rRNA duplexes needed for eukaryotic ribosome biogenesis

Short duplexes between the U3 small nucleolar RNA and the precursor ribosomal RNA must form quickly and with high yield to satisfy the high demand for ribosome synthesis in rapidly growing eukaryotic cells. These interactions, designated the U3-ETS (external transcribed spacer) and U3-18S duplexes, are essential to initiate the processing of small subunit ribosomal RNA. Previously, we showed that duplexes corresponding to those in Saccharomyces cerevisiae are only observed in vitro after addition of one of two proteins: Imp3p or Imp4p. Here, we used fluorescence-based and other in vitro assays to determine whether these proteins possess RNA chaperone activities and to assess whether these activities are sufficient to satisfy the duplex yield and rate requirements expected in vivo. Assembly of both proteins with the U3 small nucleolar RNA into a chaperone complex destabilizes a U3 stem structure, apparently to expose its 18S base-pairing site. As a result, the chaperone complex accelerates formation of the U3-18S duplex from an undetectable rate to one comparable with the intrinsic rate observed for hybridizing short duplexes. The chaperone complex also stabilizes the U3-ETS duplex by 2.7 kcal/mol. These chaperone activities provide high U3-ETS duplex yield and rapid U3-18S duplex formation over a broad concentration range to help ensure that the U3-precursor ribosomal RNA interactions limit neither ribosome biogenesis nor rapid cell growth. The thermodynamic and kinetic framework used is general and thus suitable for investigating the mechanism of action of other RNA chaperones.

Non-coding RNA annotation of the genome of Trichoplax adhaerens

A detailed annotation of non-protein coding RNAs is typically missing in initial releases of newly sequenced genomes. Here we report on a comprehensive ncRNA annotation of the genome of Trichoplax adhaerens, the presumably most basal metazoan whose genome has been published to-date. Since blast identified only a small fraction of the best-conserved ncRNAs--in particular rRNAs, tRNAs and some snRNAs--we developed a semi-global dynamic programming tool, GotohScan, to increase the sensitivity of the homology search. It successfully identified the full complement of major and minor spliceosomal snRNAs, the genes for RNase P and MRP RNAs, the SRP RNA, as well as several small nucleolar RNAs. We did not find any microRNA candidates homologous to known eumetazoan sequences. Interestingly, most ncRNAs, including the pol-III transcripts, appear as single-copy genes or with very small copy numbers in the Trichoplax genome.

Deficiency of the tRNATyr:Psi 35-synthase aPus7 in Archaea of the Sulfolobales order might be rescued by the H/ACA sRNA-guided machinery

Up to now, Psi formation in tRNAs was found to be catalysed by stand-alone enzymes. By computational analysis of archaeal genomes we detected putative H/ACA sRNAs, in four Sulfolobales species and in Aeropyrum pernix, that might guide Psi 35 formation in pre-tRNA(Tyr)(GUA). This modification is achieved by Pus7p in eukarya. The validity of the computational predictions was verified by in vitro reconstitution of H/ACA sRNPs using the identified Sulfolobus solfataricus H/ACA sRNA. Comparison of Pus7-like enzymes encoded by archaeal genomes revealed amino acid substitutions in motifs IIIa and II in Sulfolobales and A. pernix Pus7-like enzymes. These conserved RNA:Psi-synthase- motifs are essential for catalysis. As expected, the recombinant Pyrococcus abyssi aPus7 was fully active and acted at positions 35 and 13 and other positions in tRNAs, while the recombinant S. solfataricus aPus7 was only found to have a poor activity at position 13. We showed that the presence of an A residue 3' to the target U residue is required for P. abyssi aPus7 activity, and that this is not the case for the reconstituted S. solfataricus H/ACA sRNP. In agreement with the possible formation of Psi 35 in tRNA(Tyr)(GUA) by aPus7 in P. abyssi and by an H/ACA sRNP in S. solfataricus, the A36G mutation in the P. abyssi tRNA(Tyr)(GUA) abolished Psi 35 formation when using P. abyssi extract, whereas the A36G substitution in the S. solfataricus pre-tRNA(Tyr) did not affect Psi 35 formation in this RNA when using an S. solfataricus extract.

Functionality and substrate specificity of human box H/ACA guide RNAs

A large number of box H/ACA RNAs have been identified in human cells, and have been predicted to account for nearly all pseudouridylation sites in human rRNAs. However, the function of these mammalian H/ACA RNAs in directing pseudouridylation has been verified experimentally in only two cases. In this study, we used three in vitro reconstitution systems, including yeast and mammalian systems, to test the function of seven H/ACA RNAs guiding16 pseudouridylation sites. Our results verified 12 of these sites; four predictions were incorrect. Further analyses indicated that three components, including the stability of the hairpin structure harboring the pseudouridylation pocket, the stability of guide sequence-target RNA base-pairing interaction, and the distance between the target uridine and the box H or ACA, were critical for the guide function, as changes in these components were sufficient to alter the functionality and specificity of the pseudouridylation pocket. The dynamic functional changes in response to changes in these three important components were further tested in vivo, and the results were completely consistent with the in vitro results. Finally, we compared our results with predictions made by two computer programs, as well as predictions made by human experts using visual inspection. We found that the predictions of one program (snoGPS) agreed with our experimental results with 100% sensitivity (12/12) and 75% specificity (3/4).

RNase P: increased versatility through protein complexity?

Ribonuclease P (RNase P) is an essential enzyme that catalyzes the 5' endonucleolytic cleavage of precursor transfer RNAs (pretRNAs). It is found in all phylogenetic domains: bacteria, archaea and eukaryotes. The bacterial enzyme consists of a single, catalytic RNA subunit and one small protein, while the archaeal and eukaryotic enzymes have 4-10 proteins in addition to a similar RNA subunit. The bacterial RNA acts as a ribozyme at high salt in vitro; however the added protein optimizes kinetics and makes specific contacts with the pre-tRNA substrate. The bacterial protein subunit also appears to be required for the processing of non-tRNA substrates by broadening recognition tolerance. In addition, the immense increase in protein content in the eukaryotic enzymes suggests substantially enlarged capacity for recognition of additional substrates. Recently intron-encoded box C/D snoRNAs were shown to be likely substrates for RNase P, with several lines of evidence suggesting that the nuclear holoenzyme binds tightly to, and can cleave single-stranded RNA in a sequence dependent fashion. The possible involvement of RNase P in additional RNA processing or turnover pathways would be consistent with previous findings that RNase MRP, a variant of RNase P that has evolved to participate in ribosomal RNA processing, is also involved in turnover of specific messenger RNAs. Here, involvement of RNase P in multiple RNA processing pathways is discussed.