Identification of novel non-coding RNAs as potential antisense regulators in the archaeon Sulfolobus solfataricus.

By generating a specialized cDNA library from the archaeon Sulfolobus solfataricus, we have identified 57 novel small non-coding RNA (ncRNA) candidates and confirmed their expression by Northern blot analysis. The majority was found to belong to one of two classes, either antisense or antisense-box RNAs, where the latter only exhibit partial complementarity to RNA targets. The most prominent group of antisense RNAs is transcribed in the opposite orientation to the transposase genes, encoded by insertion elements (transposons). Thus, these antisense RNAs may regulate transposition of insertion elements by inhibiting expression of the transposase mRNA. Surprisingly, the class of antisense RNAs also contained RNAs complementary to tRNAs or sRNAs (small-nucleolar-like RNAs). For the antisense-box ncRNAs, the majority could be assigned to the class of C/D sRNAs, which specify 2'-O-methylation sites on rRNAs or tRNAs. Five C/D sRNAs of this group are predicted to target methylation at six sites in 13 different tRNAs, thus pointing to the widespread role of these sRNA species in tRNA modification in Archaea. Another group of antisense-box RNAs, lacking typical C/D sRNA motifs, was predicted to target the 3'-untranslated regions of certain mRNAs. Furthermore, one of the ncRNAs that does not show antisense elements is transcribed from a repeat unit of a cluster of small regularly spaced repeats in S. solfataricus which is potentially involved in replicon partitioning. In conclusion, this is the first report of stably expressed antisense RNAs in an archaeal species and it raises the prospect that antisense-based mechanisms are also used widely in Archaea to regulate gene expression.

Experimental RNomics: a global approach to identifying small nuclear RNAs and their targets in different model organisms.

Non-messenger RNAs (nmRNAs) play a wide and essential role in cellular functions. Computational identification of novel nmRNAs in genomes of model organisms is severely restricted owing to their lack of an open reading frame. Hence, we describe experimental approaches for their identification by generating cDNA libraries derived from nmRNAs for which we coined the term experimental RNomics. Two different procedures are introduced for cDNA library construction. First, we describe the construction of a general purpose cDNA library from sized RNA fractions. Second, we introduce a more specialized RNomics strategy employing this approach to generate a cDNA library from a specific abundant class of nmRNAs. This is illustrated using as a paradigm the two families of small nucleolar RNAs that guide modification of nucleotides in rRNAs or spliceosomal RNAs small nuclear RNAs (snRNAs) by short antisense elements complementary to the modification site. Following the identification of novel members from the class of small nuclear RNAs by experimental RNomics, we demonstrate how their target sequences in rRNAs or snRNAs can be identified.

In vitro RNP assembly and methylation guide activity of an unusual box C/D RNA, cis-acting archaeal pre-tRNA(Trp).

Among the large family of C/D methylation guide RNAs, the intron of euryarchaeal pre-tRNA(Trp) represents an outstanding specimen able to guide in cis, instead of in trans, two 2'-O-methylations in the pre-tRNA exons. Remarkably, both sites of methylation involve nucleotides within the bulge-helix-bulge (BHB) splicing motif, while the RNA-guided methylation and pre-tRNA splicing events depend on mutually exclusive RNA folding patterns. Using the three recombinant core proteins of archaeal C/D RNPs, we have analyzed in vitro RNP assembly of the pre-tRNA and tested its site-specific methylation activity. Recognition by L7Ae of hallmark K-turns at the C/D and C'/D' motifs appears as a crucial assembly step required for subsequent binding of a Nop5p-aFib heterodimer at each site. Unexpectedly, however, even without L7Ae but at a higher concentration of Nop5p-aFib, a substantially active RNP complex can still form, possibly reflecting the higher propensity of the cis-acting system to form guide RNA duplex(es) relative to classical trans- acting C/D RNA guides. Moreover, footprinting data of RNPs, consistent with Nop5p interacting with the non-canonical stem of the K-turn, suggest that binding of Nop5p-aFib to the pre-tRNA-L7Ae complex might direct transition from a splicing-competent structure to an RNA conformer displaying the guide RNA duplexes required for site-specific methylation.

Identification of 13 novel human modification guide RNAs.

Members of the two expanding RNA subclasses termed C/D and H/ACA RNAs guide the 2'-O-methylations and pseudouridylations, respectively, of rRNA and spliceosomal RNAs (snRNAs). Here, we report on the identification of 13 novel human intron-encoded small RNAs (U94-U106) belonging to the two subclasses of modification guides. Seven of them are predicted to direct 2'-O-methylations in rRNA or snRNAs, while the remainder represent novel orphan RNA modification guides. From these, U100, which is exclusively detected in Cajal bodies (CBs), is predicted to direct modification of a U6 snRNA uridine, U(9), which to date has not been found to be pseudouridylated. Hence, within CBs, U100 might function in the folding pathway or other aspects of U6 snRNA metabolism rather than acting as a pseudouridylation guide. U106 C/D snoRNA might also possess an RNA chaperone activity only since its two conserved antisense elements match two rRNA sequences devoid of methylated nucleotides and located remarkably close to each other within the 18S rRNA secondary structure. Finally, we have identified a retrogene for U99 snoRNA located within an intron of the Siat5 gene, supporting the notion that retro-transposition events might have played a substantial role in the mobility and diversification of snoRNA genes during evolution.

Imprinted microRNA genes transcribed antisense to a reciprocally imprinted retrotransposon-like gene.

MicroRNAs (miRNAs) are an abundant class of RNAs that are approximately 21-25 nucleotides (nt) long, interact with mRNAs and trigger either translation repression or RNA cleavage (RNA interference, RNAi) depending on the degree of complementarity with their targets. Here we show that the imprinted mouse distal chromosome 12 locus encodes two miRNA genes expressed from the maternally inherited chromosome and antisense to a retrotransposon-like gene (Rtl1) expressed only from the paternal allele.

RNomics in Drosophila melanogaster: identification of 66 candidates for novel non-messenger RNAs.

By generating a specialised cDNA library from four different developmental stages of Drosophila melanogaster, we have identified 66 candidates for small non-messenger RNAs (snmRNAs) and have confirmed their expression by northern blot analysis. Thirteen of them were expressed at certain stages of D.melanogaster development, only. Thirty-five species belong to the class of small nucleolar RNAs (snoRNAs), divided into 15 members from the C/D subclass and 20 members from the H/ACA subclass, which mostly guide 2'-O-methylation and pseudouridylation, respectively, of rRNA and snRNAs. These also include two outstanding C/D snoRNAs, U3 and U14, both functioning as pre-rRNA chaperones. Surprisingly, the sequence of the Drosophila U14 snoRNA reflects a major change of function of this snoRNA in Diptera relative to yeast and vertebrates. Among the 22 snmRNAs lacking known sequence and structure motifs, five were located in intergenic regions, two in introns, five in untranslated regions of mRNAs, eight were derived from open reading frames, and two were transcribed opposite to an intron. Interestingly, detection of two RNA species from this group implies that certain snmRNA species are processed from alternatively spliced pre-mRNAs. Surprisingly, a few snmRNA sequences could not be found on the published D.melanogaster genome, which might suggest that more snmRNA genes (as well as mRNAs) are hidden in unsequenced regions of the genome.

Binding of L7Ae protein to the K-turn of archaeal snoRNAs: a shared RNA binding motif for C/D and H/ACA box snoRNAs in Archaea.

Small nucleolar RNAs (designated as snoRNAs in Eukarya or sRNAs in Archaea) can be grouped into H/ACA or C/D box snoRNA (sRNA) subclasses. In Eukarya, H/ACA snoRNAs assemble into a ribonucleoprotein (RNP) complex comprising four proteins: Cbf5p, Gar1p, Nop10p and Nhp2p. A homolog for the Nhp2p protein has not been identified within archaeal H/ACA RNPs thus far, while potential orthologs have been identified for the other three proteins. Nhp2p is related, particularly in the middle portion of the protein sequence, to the archaeal ribosomal protein and C/D box protein L7Ae. This finding suggests that L7Ae may be able to substitute for the Nhp2p protein in archaeal H/ACA sRNAs. By band shift assays, we have analyzed in vitro the interaction between H/ACA box sRNAs and protein L7Ae from the archaeon Archaeoglobus fulgidus. We present evidence that L7Ae forms specific complexes with three different H/ACA sRNAs, designated as Afu-4, Afu-46 and Afu-190 with an apparent K(d) ranging from 28 to 100 nM. By chemical and enzymatic probing we show that distinct bases located within bulges or loops of H/ACA sRNAs interact with the L7Ae protein. These findings are corroborated by mutational analysis of the L7Ae binding site. Thereby, the RNA motif required for L7Ae binding exhibits a structure, designated as the K-turn, which is present in all C/D box sRNAs. We also identified four H/ACA RNAs from the archaeal species Pyrococcus which exhibit the K-turn motif at a similar position in their structure. These findings suggest a triple role for L7Ae protein in Archaea, e.g. in ribosomes as well as H/ACA and C/D box sRNP biogenesis and function by binding to the K-turn motif.

Plant snoRNA database.

The Plant snoRNA database (http://www.scri.sari.ac.uk/plant_snoRNA/) provides information on small nucleolar RNAs from Arabidopsis and eighteen other plant species. Information includes sequences, expression data, methylation and pseudouridylation target modification sites, initial gene organization (polycistronic, single gene and intronic) and the number of gene variants. The Arabidopsis information is divided into box C/D and box H/ACA snoRNAs, and within each of these groups, by target sites in rRNA, snRNA or unknown. Alignments of orthologous genes and gene variants from different plant species are available for many snoRNA genes. Plant snoRNA genes have been given a standard nomenclature, designed wherever possible, to provide a consistent identity with yeast and human orthologues.

Experimental RNomics: identification of 140 candidates for small non-messenger RNAs in the plant Arabidopsis thaliana.

BACKGROUND: Genomes from all organisms known to date express two types of RNA molecules: messenger RNAs (mRNAs), which are translated into proteins, and non-messenger RNAs, which function at the RNA level and do not serve as templates for translation. RESULTS: We have generated a specialized cDNA library from Arabidopsis thaliana to investigate the population of small non-messenger RNAs (snmRNAs) sized 50-500 nt in a plant. From this library, we identified 140 candidates for novel snmRNAs and investigated their expression, abundance, and developmental regulation. Based on conserved sequence and structure motifs, 104 snmRNA species can be assigned to novel members of known classes of RNAs (designated Class I snmRNAs), namely, small nucleolar RNAs (snoRNAs), 7SL RNA, U snRNAs, as well as a tRNA-like RNA. For the first time, 39 novel members of H/ACA box snoRNAs could be identified in a plant species. Of the remaining 36 snmRNA candidates (designated Class II snmRNAs), no sequence or structure motifs were present that would enable an assignment to a known class of RNAs. These RNAs were classified based on their location on the A. thaliana genome. From these, 29 snmRNA species located to intergenic regions, 3 located to intronic sequences of protein coding genes, and 4 snmRNA candidates were derived from annotated open reading frames. Surprisingly, 15 of the Class II snmRNA candidates were shown to be tissue-specifically expressed, while 12 are encoded by the mitochondrial or chloroplast genome. CONCLUSIONS: Our study has identified 140 novel candidates for small non-messenger RNA species in the plant A. thaliana and thereby sets the stage for their functional analysis.

RNomics: identification and function of small, non-messenger RNAs.

In the past few years, our knowledge about small non-mRNAs (snmRNAs) has grown exponentially. Approaches including computational and experimental RNomics have led to a plethora of novel snmRNAs, especially small nucleolar RNAs (snoRNAs). Members of this RNA class guide modification of ribosomal and spliceosomal RNAs. Novel targets for snoRNAs were identified such as tRNAs and potentially mRNAs, and several snoRNAs were shown to be tissue-specifically expressed. In addition, previously unknown classes of snmRNAs have been discovered. MicroRNAs and small interfering RNAs of about 21-23 nt, were shown to regulate gene expression by binding to mRNAs via antisense elements. Regulation of gene expression is exerted by degradation of mRNAs or translational regulation. snmRNAs play a variety of roles during regulation of gene expression. Moreover, the function of some snmRNAs known for decades, has been finally elucidated. Many other RNAs were identified by RNomics studies lacking known sequence and structure motifs. Future challenges in the field of RNomics include identification of the novel snmRNA's biological roles in the cell.

The expanding snoRNA world.

In eukaryotes, the site-specific formation of the two prevalent types of rRNA modified nucleotides, 2'-O-methylated nucleotides and pseudouridines, is directed by two large families of snoRNAs. These are termed box C/D and H/ACA snoRNAs, respectively, and exert their function through the formation of a canonical guide RNA duplex at the modification site. In each family, one snoRNA acts as a guide for one, or at most two modifications, through a single, or a pair of appropriate antisense elements. The two guide families now appear much larger than anticipated and their role not restricted to ribosome synthesis only. This is reflected by the recent detection of guides that can target other cellular RNAs, including snRNAs, tRNAs and possibly even mRNAs, and by the identification of scores of tissue-specific specimens in mammals. Recent characterization of homologs of eukaryotic modification guide snoRNAs in Archaea reveals the ancient origin of these non-coding RNA families and offers new perspectives as to their range of function.

Identification of tandemly-repeated C/D snoRNA genes at the imprinted human 14q32 domain reminiscent of those at the Prader-Willi/Angelman syndrome region.

A human imprinted domain at 14q32 contains two co-expressed and reciprocally imprinted genes, DLK1 and GTL2, which are expressed from the paternally and maternally inherited alleles, respectively. We have previously shown that another imprinted locus, on human 15q11-q13, contains a large number of tandemly repeated C/D small nucleolar RNA genes (or C/D snoRNAs) only expressed from the paternal allele. Here we show that the region downstream from the GTL2 gene is also characterized by a transcription unit spanning many repeated intron-encoded C/D snoRNA genes, most of them arranged into two tandem arrays of 31 and 9 copies. Intriguingly, these snoRNAs depart from previously reported rRNA or snRNA methylation guides by their tissue-specific expression and by their lack of complementarity to rRNA or snRNA within their sequences. Analysis of the orthologous region in the mouse shows that the previously reported maternally expressed Rian gene, located downstream of Gtl2 on the distal 12 chromosome, encodes at least nine C/D snoRNAs. Through a systematic search in rodents, we could identify other C/D snoRNA genes in this domain. All snoRNAs identified on mouse distal 12 are brain-specific and only expressed from the maternally inherited allele. The human imprinted 14q32 domain therefore shares common genomic features with the imprinted 15q11-q13 loci. This link between tandemly repeated C/D snoRNA genes and genomic imprinting suggests a role for these snoRNAs and/or their host non-coding RNA genes in the evolution and/or mechanism of the epigenetic imprinting process.

Identification of 86 candidates for small non-messenger RNAs from the archaeon Archaeoglobus fulgidus.

In a specialized cDNA library from the archaeon Archaeoglobus fulgidus we have identified a total of 86 different expressed RNA sequences potentially encoding previously uncharacterized small non-messenger RNA (snmRNA) species. Ten of these RNAs resemble eukaryotic small nucleolar RNAs, which guide rRNA 2'-O-methylations (C/D-box type) and pseudouridylations (H/ACA-box type). Thereby, we identified four candidates for H/ACA small RNAs in an archaeal species that are predicted to guide a total of six rRNA pseudouridylations. Furthermore, we have verified the presence of the six predicted pseudouridines experimentally. We demonstrate that 22 snmRNAs are transcribed from a family of short tandem repeats conserved in most archaeal genomes and shown previously to be potentially involved in replicon partitioning. In addition, four snmRNAs derived from the rRNA operon of A. fulgidus were identified and shown to be generated by a splicing/processing pathway of pre-rRNAs. The remaining 50 RNAs could not be assigned to a known class of snmRNAs because of the lack of known structure and/or sequence motifs. Regarding their location on the genome, only nine were located in intergenic regions, whereas 33 were complementary to an ORF, five were overlapping an ORF, and three were derived from the sense orientation within an ORF. Our study further supports the importance of snmRNAs in all three domains of life.

RNomics in Archaea reveals a further link between splicing of archaeal introns and rRNA processing.

The bulge-helix-bulge (BHB) motif recognised by the archaeal splicing endonuclease is also found in the long processing stems of archaeal rRNA precursors in which it is cleaved to generate pre-16S and pre-23S rRNAs. We show that in two species, Archaeoglobus fulgidus and Sulfolobus solfataricus, representatives from the two major archaeal kingdoms Euryarchaeota and Crenarchaeota, respectively, the pre-rRNA spacers cleaved at the BHB motifs surrounding pre-16S and pre-23S rRNAs subsequently become ligated. In addition, we present evidence that this is accompanied by circularization of ribosomal pre-16S and pre-23S rRNAs in both species. These data reveal a further link between intron splicing and pre-rRNA processing in Archaea, which might reflect a common evolutionary origin of the two processes. One spliced RNA species designated 16S-D RNA, resulting from religation at the BHB motif of 16S pre-rRNA, is a highly abundant and stable RNA which folds into a three-stem structure interrupted by two single-stranded regions as assessed by chemical probing. It spans a region of the pre-rRNA 5' external transcribed spacer exhibiting a highly conserved folding pattern in Archaea. Surprisingly, 16S-D RNA contains structural motifs found in archaeal C/D box small RNAs and binds to the L7Ae protein, a core component of archaeal C/D box RNPs. This supports the notion that it might have an important but still unknown role in pre-rRNA biogenesis or might even target RNA molecules other than rRNA.