Three new small nucleolar RNAs that are psoralen cross-linked in vivo to unique regions of pre-rRNA

The RNA Base-Pairing Problem and Base-Pairing Solutions

RNA molecules are folded into structures and complexes to perform a wide variety of functions. Determination of RNA structures and their interactions is a fundamental problem in RNA biology. Most RNA molecules in living cells are large and dynamic, posing unique challenges to structure analysis. Here we review progress in RNA structure analysis, focusing on methods that use the "cross-link, proximally ligate, and sequence" principle for high-throughput detection of base-pairing interactions in living cells. Beginning with a comparison of commonly used methods in structure determination and a brief historical account of psoralen cross-linking studies, we highlight the important features of cross-linking methods and new biological insights into RNA structures and interactions from recent studies. Further improvement of these cross-linking methods and application to previously intractable problems will shed new light on the mechanisms of the "modern RNA world."

RNA interactomics: recent advances and remaining challenges

Tight regulation of cellular processes is key to the development of complex organisms but also vital for simpler ones. During evolution, different regulatory systems have emerged, among them RNA-based regulation that is carried out mainly by intramolecular and intermolecular RNA-RNA interactions. However, methods for the transcriptome-wide detection of these interactions were long unavailable. Recently, three publications described high-throughput methods to directly detect RNA duplexes in living cells. This promises to enable in-depth studies of RNA-based regulation and will narrow the gaps in our understanding of RNA structure and function. In this review, we highlight the benefits of these methods and their commonalities and differences and, in particular, point to methodological shortcomings that hamper their wider application. We conclude by presenting ideas for how to overcome these problems and commenting on the prospects we see in this area of research.

Matrin 3 binds and stabilizes mRNA

Matrin 3 (MATR3) is a highly conserved, inner nuclear matrix protein with two zinc finger domains and two RNA recognition motifs (RRM), whose function is largely unknown. Recently we found MATR3 to be phosphorylated by the protein kinase ATM, which activates the cellular response to double strand breaks in the DNA. Here, we show that MATR3 interacts in an RNA-dependent manner with several proteins with established roles in RNA processing, and maintains its interaction with RNA via its RRM2 domain. Deep sequencing of the bound RNA (RIP-seq) identified several small noncoding RNA species. Using microarray analysis to explore MATR3′s role in transcription, we identified 77 transcripts whose amounts depended on the presence of MATR3. We validated this finding with nine transcripts which were also bound to the MATR3 complex. Finally, we demonstrated the importance of MATR3 for maintaining the stability of several of these mRNA species and conclude that it has a role in mRNA stabilization. The data suggest that the cellular level of MATR3, known to be highly regulated, modulates the stability of a group of gene transcripts.

Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements

As an accompanying manuscript to the release of the honey bee genome, we report the entire sequence of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) ribosomal RNA (rRNA)-encoding gene sequences (rDNA) and related internally and externally transcribed spacer regions of Apis mellifera (Insecta: Hymenoptera: Apocrita). Additionally, we predict secondary structures for the mature rRNA molecules based on comparative sequence analyses with other arthropod taxa and reference to recently published crystal structures of the ribosome. In general, the structures of honey bee rRNAs are in agreement with previously predicted rRNA models from other arthropods in core regions of the rRNA, with little additional expansion in non-conserved regions. Our multiple sequence alignments are made available on several public databases and provide a preliminary establishment of a global structural model of all rRNAs from the insects. Additionally, we provide conserved stretches of sequences flanking the rDNA cistrons that comprise the externally transcribed spacer regions (ETS) and part of the intergenic spacer region (IGS), including several repetitive motifs. Finally, we report the occurrence of retrotransposition in the nuclear large subunit rDNA, as R2 elements are present in the usual insertion points found in other arthropods. Interestingly, functional R1 elements usually present in the genomes of insects were not detected in the honey bee rRNA genes. The reverse transcriptase products of the R2 elements are deduced from their putative open reading frames and structurally aligned with those from another hymenopteran insect, the jewel wasp Nasonia (Pteromalidae). Stretches of conserved amino acids shared between Apis and Nasonia are illustrated and serve as potential sites for primer design, as target amplicons within these R2 elements may serve as novel phylogenetic markers for Hymenoptera. Given the impending completion of the sequencing of the Nasonia genome, we expect our report eventually to shed light on the evolution of the hymenopteran genome within higher insects, particularly regarding the relative maintenance of conserved rDNA genes, related variable spacer regions and retrotransposable elements.

The vertebrate E1/U17 small nucleolar ribonucleoprotein particle

Each of the many different box H/ACA ribonucleoprotein particles (RNPs) present in eukaryotes and archaea consists of four common core proteins and one specific H/ACA small RNA, which bears the sequence elements H (ANANNA) and ACA. Most of the H/ACA RNPs are small nucleolar RNPs (snoRNPs), which are localized in nucleoli, and are one of the two major classes of snoRNPs. Most H/ACA RNPs direct pseudouridine synthesis in pre-rRNA and other RNAs. One H/ACA small nucleolar RNA (snoRNA), vertebrate E1/U17 (snR30 in yeast), is required for pre-rRNA cleavage processing that generates mature 18S rRNA. E1 snoRNA is encoded in introns of protein-coding genes, and the evidence suggests that human E1 RNA undergoes uridine insertional RNA editing. The vertebrate E1 RNA consensus secondary structure shows several features that are absent in other box H/ACA snoRNAs. The available UV-induced RNA-protein crosslinking results suggest that the E1 snoRNP is asymmetrical in vertebrate cells, in contrast to other H/ACA snoRNPs. The vertebrate E1 snoRNP in cells is surprisingly complex: (i) E1 RNA contacts directly and specifically several proteins which do not appear to be any of the H/ACA RNP four core proteins; and (ii) multiple E1 RNA sites are needed for E1 snoRNP formation, E1 RNA stability, and E1 RNA-protein direct interactions.

Predicted Secondary Structure for 28S and 18S rRNA from Ichneumonoidea (Insecta: Hymenoptera: Apocrita): Impact on Sequence Alignment and Phylogeny Estimation

We utilize the secondary structural properties of the 28S rRNA D2-D10 expansion segments to hypothesize a multiple sequence alignment for major lineages of the hymenopteran superfamily Ichneumonoidea (Braconidae, Ichneumonidae). The alignment consists of 290 sequences (originally analyzed in Belshaw and Quicke, Syst Biol 51:450-477, 2002) and provides the first global alignment template for this diverse group of insects. Predicted structures for these expansion segments as well as for over half of the 18S rRNA are given, with highly variable regions characterized and isolated within conserved structures. We demonstrate several pitfalls of optimization alignment and illustrate how these are potentially addressed with structure-based alignments. Our global alignment is presented online at (http://hymenoptera.tamu.edu/rna) with summary statistics, such as basepair frequency tables, along with novel tools for parsing structure-based alignments into input files for most commonly used phylogenetic software. These resources will be valuable for hymenopteran systematists, as well as researchers utilizing rRNA sequences for phylogeny estimation in any taxon. We explore the phylogenetic utility of our structure-based alignment by examining a subset of the data under a variety of optimality criteria using results from Belshaw and Quicke (2002) as a benchmark.

A Secondary Structural Model of the 28S rRNA Expansion Segments D2 and D3 for Chalcidoid Wasps (Hymenoptera: Chalcidoidea)

We analyze the secondary structure of two expansion segments (D2, D3) of the 28S ribosomal (rRNA)-encoding gene region from 527 chalcidoid wasp taxa (Hymenoptera: Chalcidoidea) representing 18 of the 19 extant families. The sequences are compared in a multiple sequence alignment, with secondary structure inferred primarily from the evidence of compensatory base changes in conserved helices of the rRNA molecules. This covariation analysis yielded 36 helices that are composed of base pairs exhibiting positional covariation. Several additional regions are also involved in hydrogen bonding, and they form highly variable base-pairing patterns across the alignment. These are identified as regions of expansion and contraction or regions of slipped-strand compensation. Additionally, 31 single-stranded locales are characterized as regions of ambiguous alignment based on the difficulty in assigning positional homology in the presence of multiple adjacent indels. Based on comparative analysis of these sequences, the largest genetic study on any hymenopteran group to date, we report an annotated secondary structural model for the D2, D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related apocritan taxa.

UV-crosslinking of E1 small nucleolar RNA to proteins in frog oocytes

E1/U17 small nucleolar RNA (snoRNA) is a box H/ACA snoRNA. To detect protein bands that UV-crosslink to E1 RNA primarily at uridines, frog oocytes were injected with [alpha-32P]UTP-labeled E1 RNA and incubated, isolated nuclei were UV irradiated, and nuclear contents were digested with RNase A. Wild-type E1 RNA specifically UV-crosslinked to several protein bands. To identify E1 RNA sites involved in these interactions, we tested 21 E1 RNA mutants, each consisting of substitutions in a conserved sequence or structure. UV-crosslinking of different protein bands to E1 RNA depended on one of the following sets of conserved E1 RNA segments: two 5' end RNA sites; five 5' half RNA sites; two 3' half RNA sites; or 14 sites located throughout E1 RNA. Of these conserved E1 RNA sites, UV-crosslinking apparently depended on sequences at 11 sites, and structures at 2 sites. Gel electrophoresis with and without RNA competition detected protein bands that are not common to all of the box H/ACA snoRNAs.

A secondary structural model of the 28S rRNA expansion segments D2 and D3 from rootworms and related leaf beetles (Coleoptera: Chrysomelidae; Galerucinae)

We analysed the secondary structure of two expansion segments (D2, D3) of the 28S rRNA gene from 229 leaf beetles (Coleoptera: Chrysomelidae), the majority of which are in the subfamily Galerucinae. The sequences were compared in a multiple sequence alignment, with secondary structure inferred primarily from the compensatory base changes in the conserved helices of the rRNA molecules. This comparative approach yielded thirty helices comprised of base pairs with positional covariation. Based on these leaf beetle sequences, we report an annotated secondary structural model for the D2 and D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related beetle taxa. This predicted structure, consisting of seven major compound helices, is mostly consistent with previously proposed models for the D2 and D3 expansion segments in insects. Despite a lack of conservation in the primary structure of these regions of insect 28S rRNA, the evolution of the secondary structure of these seven major motifs may be informative above the nucleotide level for higher-order phylogeny reconstruction of major insect lineages.

Multiple conserved segments of E1 small nucleolar RNA are involved in the formation of a ribonucleoprotein particle in frog oocytes

E1/U17 small nucleolar RNA (snoRNA) is a box H/ACA snoRNA. To identify E1 RNA elements required for its assembly into a ribonucleoprotein (RNP) particle, we have made substitution mutations in evolutionarily conserved sequences and structures of frog E1 RNA. After E1 RNA was injected into the nucleus of frog oocytes, assembly of this exogenous RNA into an RNP was monitored by non-denaturing gel electrophoresis. Unexpectedly, nucleotide substitutions in many phylogenetically conserved segments of E1 RNA produced RNPs with abnormal gel-electrophoresis patterns. These RNA segments were at least nine conserved sequences and an apparently conserved structure. In another region needed for RNP formation, the requirement may be sequence(s) and/or structure. Base substitutions in each of these and in one additional conserved E1 RNA segment reduced the stability of this snoRNA in frog oocytes. Nucleolar localization was assayed by fluorescence microscopy after injection of fluorescein-labelled RNA. The H box (ANANNA) and the ACA box are both needed for efficient nucleolar localization of frog E1 RNA.

Box H and box ACA are nucleolar localization elements of U17 small nucleolar RNA

The nucleolar localization elements (NoLEs) of U17 small nucleolar RNA (snoRNA), which is essential for rRNA processing and belongs to the box H/ACA snoRNA family, were analyzed by fluorescence microscopy. Injection of mutant U17 transcripts into Xenopus laevis oocyte nuclei revealed that deletion of stems 1, 2, and 4 of U17 snoRNA reduced but did not prevent nucleolar localization. The deletion of stem 3 had no adverse effect. Therefore, the hairpins of the hairpin-hinge-hairpin-tail structure formed by these stems are not absolutely critical for nucleolar localization of U17, nor are sequences within stems 1, 3, and 4, which may tether U17 to the rRNA precursor by base pairing. In contrast, box H and box ACA are major NoLEs; their combined substitution or deletion abolished nucleolar localization of U17 snoRNA. Mutation of just box H or just the box ACA region alone did not fully abolish the nucleolar localization of U17. This indicates that the NoLEs of the box H/ACA snoRNA family function differently from the bipartite NoLEs (conserved boxes C and D) of box C/D snoRNAs, where mutation of either box alone prevents nucleolar localization.

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

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

Site-specific pseudouridine formation in preribosomal RNA is guided by small nucleolar RNAs

During the nucleolar maturation of eukaryotic ribosomal RNAs, many selected uridines are converted into pseudouridine by a thus far undefined mechanism. The nucleolus contains a large number of small RNAs (snoRNAs) that share two conserved sequence elements, box H and ACA. In this study, we demonstrate that site-specific pseudouridylation of rRNAs relies on short ribosomal signal sequences that are complementary to sequences in box H/ACA snoRNAs. Genetic depletion and reconstitution studies on yeast snR5 and snR36 snoRNAs demonstrate that box H/ACA snoRNAs function as guide RNAs in rRNA pseudouridylation. These results define a novel function for snoRNAs and further reinforce the idea that base pairing is the most common way to obtain specific substrate-"enzyme" interactions during rRNA maturation.

Intracellular localization and unique conserved sequences of three small nucleolar RNAs

Three human small nucleolar RNAs (snoRNAs), E1, E2 and E3, were reported earlier that have unique sequences, interact directly with unique segments of pre-rRNA in vivo and are encoded in introns of protein genes. In the present report, human and frog E1, E2 and E3 RNAs injected into the cytoplasm of frog oocytes migrated to the nucleus and specifically to the nucleolus. This indicates that the nucleolar and nuclear localization signals of these snoRNAs reside within their evolutionarily conserved segments. Homologs of these snoRNAs from several vertebrates were sequenced and this information was used to develop RNA secondary structure models. These snoRNAs have unique phylogenetically conserved sequences.

Elements essential for processing intronic U14 snoRNA are located at the termini of the mature snoRNA sequence and include conserved nucleotide boxes C and D

Essential elements for intronic U14 processing have been analyzed by microinjecting various mutant hsc70/Ul4 pre-mRNA precursors into Xenopus oocyte nuclei. Initial truncation experiments revealed that elements sufficient for U14 processing are located within the mature snoRNA sequence itself. Subsequent deletions within the U14 coding region demonstrated that only the terminal regions of the folded U14 molecule containing con- served nucleotide boxes C and D are required for processing. Mutagenesis of either box C or box D completely blocked U14 processing. The importance of boxes C and D was confirmed with the excision of appropriately sized U3 and U8 fragments containing boxes C and D from an hsc7O pre-mRNA intron. Competition studies indicate that a trans-acting factor (protein?) is binding this terminal motif and is essential for U14 processing. Competition studies also revealed that this factor is common to both intronic and non-intronic snoRNAs possessing nucleotide boxes C and D. Immunoprecipitation of full-length and internally deleted U14 snoRNA molecules demonstrated that the terminal region containing boxes C and D does not bind fibrillarin. Collectively, our results indicate that a trans-acting factor (different from fibrillarin) binds to the box C- and D-containing terminal motif of U14 snoRNA, thereby stabilizing the intronic snoRNA sequence in an RNP complex during processing.

A mammalian gene with introns instead of exons generating stable RNA products

The nucleoli of eukaryotic cells are the sites of ribosomal RNA transcription and processing and of ribosomal subunit assembly. They contain multiple small nucleolar RNAs (snoRNAs), several of which are essential for rRNA maturation. The U3, U8 and U13 snoRNA genes are transcribed independently, whereas U14-U24, as well as E3, are located within introns of protein-coding genes, most of whose functions are linked to translation. These snoRNAs are co-transcribed with their host pre-mRNAs and released by processing from excised introns. Here we show that, in addition to U22, seven novel fibrillarin-associated snoRNAs, named U25-U31, are encoded within different introns of the unusually compact mammalian U22 host gene (UHG). All seven RNAs exhibit extensive (12-15 nucleotides) complementarity to different segments of the mature rRNAs, followed by a C/AUGA ('U-turn') sequence. The spliced UHG RNA, although it is associated with polysomes, has little potential for protein coding, is short-lived, and is poorly conserved between human and mouse. Thus, the introns rather than the exons specify the functional products of UHG.

Small nucleolar RNA

A growing list of small nucleolar RNAs (snoRNAs) has been characterized in eukaryotes. They are transcribed by RNA polymerase II or III; some snoRNAs are encoded in the introns of other genes. The nonintronic polymerase II transcribed snoRNAs receive a trimethylguanosine cap, probably in the nucleus, and move to the nucleolus. snoRNAs are complexed with proteins, sometimes including fibrillarin. Localization and maintenance in the nucleolus of some snoRNAs requires the presence of initial precursor rRNA (pre-rRNA). Many snoRNAs have conserved sequence boxes C and D and a 3' terminal stem; the role of these features are discussed. Functional assays done for a few snoRNAs indicate their roles in rRNA processing for cleavage of the external and internal transcribed spacers (ETS and ITS). U3 is the most abundant snoRNA and is needed for cleavage of ETS1 and ITS1; experimental results on U3 binding sites in pre-rRNA are reviewed. 18S rRNA production also needs U14, U22, and snR30 snoRNAs, whereas U8 snoRNA is needed for 5.8S and 28S rRNA production. Other snoRNAs that are complementary to 18S or 28S rRNA might act as chaperones to mediate RNA folding. Whether snoRNAs join together in a large rRNA processing complex (the "processome") is not yet clear. It has been hypothesized that such complexes could anchor the ends of loops in pre-rRNA containing 18S or 28S rRNA, thereby replacing base-paired stems found in pre-rRNA of prokaryotes.

Intron-encoded small nucleolar RNAs: new RNA sequence variants and genomic loci

Three small nucleolar RNAs (snoRNAs) whose 5' termini are monophosphorylated, termed E1, E2 and E3, were reported earlier, and E1 and E3 are encoded in pre-mRNA gene introns. In the present work, the ends of these snoRNAs were identified by analysis of terminal mononucleotides, and heterogeneity was observed at the 3' ends of E1 and E2 RNAs. Two new E1 RNA species were detected in HeLa cells by cDNA cloning. Four novel human genomic loci were identified that have E1 or E3 sequence homology. The sequence CTAGAGCACYSAATCTGGAT (where S = C or G and Y = C or T), that is present three nucleotides downstream from the coding region of an E1 RNA-encoding gene, lies in the same location in a different human genomic locus (which has a E1-homology sequence whose expression has not been detected yet), suggesting that this sequence may be functional.

Processing of eukaryotic ribosomal RNA

In summary, it can be argued that the understanding of eukaryotic rRNA processing is no less important than the understanding of mRNA maturation, since the capacity of a cell to carry out protein synthesis is controlled, in part, by the abundance of ribosomes. Processing of pre-rRNA is highly regulated, involving many cellular components acting either alone or as part of a complex. Some of these components are directly involved in the modification and cleavage of the precursor rRNA, while others direct the packaging of the rRNA into ribosome subunits. As is the case for pre-mRNA processing, snoRNPs are clearly involved in eukaryotic rRNA processing, and have been proposed to assemble with other proteins into at least one complex called a "processosome" (17), which carries out the ordered processing of the pre-rRNA and its assembly into ribosomes. The formation of a processing complex clearly makes possible the regulation required to coordinate the abundance of ribosomes with the physiological and developmental changes of a cell. It may be that eukaryotic rRNA processing is even more complex than pre-mRNA maturation, since pre-rRNA undergoes extensive nucleotide modification and is assembled into a complex structure called the ribosome. Undoubtedly, features of the eukaryotic rRNA-processing pathway have been conserved evolutionarily, and the genetic approach available in yeast research (6) should provide considerable knowledge that will be useful for other investigators working with higher eukaryotic systems. Interestingly, it was originally hoped that the extensive work and understanding of bacterial ribosome formation would provide a useful paradigm for the process in eukaryotes. However, although general features of ribosome structure and function are highly conserved between bacterial and eukaryotic systems, the basic strategy in ribosome biogenesis seems to be, for the most part, distinctly different. Thus, the detailed molecular mechanisms for rRNA processing in each kingdom will have to be independently deciphered in order to elucidate the features and regulation of this important process for cell survival.