Detection of 5'- and 3'-UTR-derived small RNAs and cis-encoded antisense RNAs in Escherichia coli

Evidence is accumulating that small, noncoding RNAs are important regulatory molecules. Computational and experimental searches have led to the identification of ~60 small RNA genes in Escherichia coli. However, most of these studies focused on the intergenic regions and assumed that small RNAs were >50 nt. Thus, the previous screens missed small RNAs encoded on the antisense strand of protein-coding genes and small RNAs of <50 nt. To identify additional small RNAs, we carried out a cloning-based screen focused on RNAs of 30–65 nt. In this screen, we identified RNA species corresponding to fragments of rRNAs, tRNAs and known small RNAs. Several of the small RNAs also corresponded to 5′- and 3′-untranslated regions (UTRs) and internal fragments of mRNAs. Four of the 3′-UTR-derived RNAs were highly abundant and two showed expression patterns that differed from the corresponding mRNAs, suggesting independent functions for the 3′-UTR-derived small RNAs. We also detected three previously unidentified RNAs encoded in intergenic regions and RNAs from the long direct repeat and hok/sok elements. In addition, we identified a few small RNAs that are expressed opposite protein-coding genes and could base pair with 5′ or 3′ ends of the mRNAs with perfect complementarity.     

Small RNAs derived from snoRNAs

Small nucleolar RNAs (snoRNAs) guide RNA modification and are localized in nucleoli and Cajal bodies in eukaryotic cells. Components of the RNA silencing pathway associate with these structures, and two recent reports have revealed that a human and a protozoan snoRNA can be processed into miRNA-like RNAs. Here we show that small RNAs with evolutionary conservation of size and position are derived from the vast majority of snoRNA loci in animals (human, mouse, chicken, fruit fly), Arabidopsis, and fission yeast. In animals, sno-derived RNAs (sdRNAs) from H/ACA snoRNAs are predominantly 20–24 nucleotides (nt) in length and originate from the 3′ end. Those derived from C/D snoRNAs show a bimodal size distribution at ∼17–19 nt and >27 nt and predominantly originate from the 5′ end. SdRNAs are associated with AGO7 in Arabidopsis and Ago1 in fission yeast with characteristic 5′ nucleotide biases and show altered expression patterns in fly loquacious and Dicer-2 and mouse Dicer1 and Dgcr8 mutants. These findings indicate that there is interplay between the RNA silencing and snoRNA-mediated RNA processing systems, and that sdRNAs comprise a novel and ancient class of small RNAs in eukaryotes.     

Molecular characterization of geminivirus-derived small RNAs in different plant species

DNA geminiviruses are thought to be targets of RNA silencing. Here, we characterize small interfering (si) RNAs—the hallmarks of silencing—associated with Cabbage leaf curl begomovirus in Arabidopsis and African cassava mosaic begomovirus in Nicotiana benthamiana and cassava. We detected 21, 22 and 24 nt siRNAs of both polarities, derived from both the coding and the intergenic regions of these geminiviruses. Genetic evidence showed that all the 24 nt and a substantial fraction of the 22 nt viral siRNAs are generated by the dicer-like proteins DCL3 and DCL2, respectively. The viral siRNAs were 5′ end phosphorylated, as shown by phosphatase treatments, and methylated at the 3′-nucleotide, as shown by HEN1 miRNA methylase-dependent resistance to β-elimination. Similar modifications were found in all types of endogenous and transgene-derived siRNAs tested, but not in a major fraction of siRNAs from a cytoplasmic RNA tobamovirus. We conclude that several distinct silencing pathways are involved in DNA virus-plant interactions.     

Small RNAs with 5′-Polyphosphate Termini Associate with a Piwi-Related Protein and Regulate Gene Expression in the Single-Celled Eukaryote Entamoeba histolytica

Small interfering RNAs regulate gene expression in diverse biological processes, including heterochromatin formation and DNA elimination, developmental regulation, and cell differentiation. In the single-celled eukaryote Entamoeba histolytica, we have identified a population of small RNAs of 27 nt size that (i) have 5′-polyphosphate termini, (ii) map antisense to genes, and (iii) associate with an E. histolytica Piwi-related protein. Whole genome microarray expression analysis revealed that essentially all genes to which antisense small RNAs map were not expressed under trophozoite conditions, the parasite stage from which the small RNAs were cloned. However, a number of these genes were expressed in other E. histolytica strains with an inverse correlation between small RNA and gene expression level, suggesting that these small RNAs mediate silencing of the cognate gene. Overall, our results demonstrate that E. histolytica has an abundant 27 nt small RNA population, with features similar to secondary siRNAs from C. elegans, and which appear to regulate gene expression. These data indicate that a silencing pathway mediated by 5′-polyphosphate siRNAs extends to single-celled eukaryotic organisms.     

The Putative Nucleic Acid Helicase Sen1p Is Required for Formation and Stability of Termini and for Maximal Rates of Synthesis and Levels of Accumulation of Small Nucleolar RNAs in Saccharomyces cerevisiae

Sen1p from Saccharomyces cerevisiae is a nucleic acid helicase related to DEAD box RNA helicases and type I DNA helicases. The temperature-sensitive sen1-1 mutation located in the helicase motif alters the accumulation of pre-tRNAs, pre-rRNAs, and some small nuclear RNAs. In this report, we show that cells carrying sen1-1 exhibit altered accumulation of several small nucleolar RNAs (snoRNAs) immediately upon temperature shift. Using Northern blotting, RNase H cleavage, primer extension, and base compositional analysis, we detected three forms of the snoRNA snR13 in wild-type cells: an abundant TMG-capped 124-nucleotide (nt) mature form (snR13F) and two less abundant RNAs, including a heterogeneous population of ~1,400-nt 3′-extended forms (snR13R) and a 108-nt 5′-truncated form (snR13T) that is missing 16 nt at the 5′ end. A subpopulation of snR13R contains the same 5′ truncation. Newly synthesized snR13R RNA accumulates with time at the expense of snR13F following temperature shift of sen1-1 cells, suggesting a possible precursor-product relationship. snR13R and snR13T both increase in abundance at the restrictive temperature, indicating that Sen1p stabilizes the 5′ end and promotes maturation of the 3′ end. snR13F contains canonical C and D boxes common to many snoRNAs. The 5′ end of snR13T and the 3′ end of snR13F reside within C₂U₄ sequences that immediately flank the C and D boxes. A mutation in the 5′ C₂U₄ repeat causes underaccumulation of snR13F, whereas mutations in the 3′ C₂U₄ repeat cause the accumulation of two novel RNAs that migrate in the 500-nt range. At the restrictive temperature, double mutants carrying sen1-1 and mutations in the 3′ C₂U₄ repeat show reduced accumulation of the novel RNAs and increased accumulation of snR13R RNA, indicating that Sen1p and the 3′ C₂U₄ sequence act in a common pathway to facilitate 3′ end formation. Based on these findings, we propose that Sen1p and the C₂U₄ repeats that flank the C and D boxes promote maturation of the 3′ terminus and stability of the 5′ terminus and are required for maximal rates of synthesis and levels of accumulation of mature snR13F.     

Minimal-length short hairpin RNAs: The relationship of structure and RNAi activity

Small hairpin RNAs (shRNAs) are widely used in RNAi studies and typically consist of a stem of 19–29 base pairs (bp), a loop of at least 4 nucleotides (nt), and a dinucleotide overhang at the 3′ end. Compared with shRNAs with 21–29 bp stems, we have found that shRNAs with 19-bp or shorter stems (sshRNAs) possess some unique structure–activity features that depend on whether the antisense strand is positioned 5′ or 3′ to the loop (L- or R-type sshRNAs, respectively). L sshRNAs can have IC₅₀s in the very low picomolar range, and sshRNAs with nominal loop sizes of 1 or 4 nt were at least as active as those with longer loops. L sshRNAs remained highly potent even when the 3′ end of the antisense strand was directly linked with the 5′ end of the sense strand. In this case, the sense strand can be shorter than the antisense strand, and the loop can be formed entirely by the 3′ end of the antisense strand. Monomer sshRNAs are not processed by recombinant Dicers in vitro. Although they can form dimers that are sometimes Dicer substrates, their RNAi activity is not dependent on the formation of such structures. Our findings have implications for the mechanism of action of sshRNAs, and the ability to design highly potent shRNAs with minimal length is encouraging for the prospects of the therapeutic use of direct-delivered shRNAs.     

SL1 trans Splicing and 3′-End Formation in a Novel Class of Caenorhabditis elegans Operon

Many Caenorhabditis elegans genes exist in operons in which polycistronic precursors are processed by cleavage at the 3′ ends of upstream genes and trans splicing 100 to 400 nucleotides away, at the 5′ ends of downstream genes, to generate monocistronic messages. Of the two spliced leaders, SL1 is trans spliced to the 5′ ends of upstream genes, whereas SL2 is reserved for downstream genes in operons. However, there are isolated examples of what appears to be a different sort of operon, in which trans splicing is exclusively to SL1 and there is no intercistronic region; the polyadenylation signal is only a few base pairs upstream of the trans-splice site. We have analyzed the processing of an operon of this type by inserting the central part of mes-6/cks-1 into an SL2-type operon. In this novel context, cks-1 is trans spliced only to SL1, and mes-6 3′-end formation occurs normally, demonstrating that this unique mode of processing is indeed intrinsic to this kind of operon, which we herein designate “SL1-type.” An exceptionally long polypyrimidine tract found in the 3′ untranslated regions of the three known SL1-type operons is shown to be required for the accumulation of both upstream and downstream mRNAs. Mutations of the trans-splice and poly(A) signals indicate that the two processes are independent and in competition, presumably due to their close proximity, raising the possibility that production of upstream and downstream mRNAs is mutually exclusive.     

Mining small RNA sequencing data: a new approach to identify small nucleolar RNAs in Arabidopsis

Small nucleolar RNAs (snoRNAs) are noncoding RNAs that direct 2′-O-methylation or pseudouridylation on ribosomal RNAs or spliceosomal small nuclear RNAs. These modifications are needed to modulate the activity of ribosomes and spliceosomes. A comprehensive repertoire of snoRNAs is needed to expand the knowledge of these modifications. The sequences corresponding to snoRNAs in 18–26-nt small RNA sequencing data have been rarely explored and remain as a hidden treasure for snoRNA annotation. Here, we showed the enrichment of small RNAs at Arabidopsis snoRNA termini and developed a computational approach to identify snoRNAs on the basis of this characteristic. The approach successfully uncovered the full-length sequences of 144 known Arabidopsis snoRNA genes, including some snoRNAs with improved 5′- or 3′-end annotation. In addition, we identified 27 and 17 candidates for novel box C/D and box H/ACA snoRNAs, respectively. Northern blot analysis and sequencing data from parallel analysis of RNA ends confirmed the expression and the termini of the newly predicted snoRNAs. Our study especially expanded on the current knowledge of box H/ACA snoRNAs and snoRNA species targeting snRNAs. In this study, we demonstrated that the use of small RNA sequencing data can increase the complexity and the accuracy of snoRNA annotation.     

Characterisation of the U83 and U84 small nucleolar RNAs: two novel 2'-O-ribose methylation guide RNAs that lack complementarities to ribosomal RNAs

In eukaryotic cells, the site-specific 2′-O-ribose methy-lation of ribosomal RNAs (rRNAs) and the U6 spliceosomal small nuclear RNA (snRNA) is directed by small nucleolar RNAs (snoRNAs). The C and D box-containing 2′-O-methylation guide snoRNAs select the correct substrate nucleotide through formation of a long 10–21 bp interaction with the target rRNA and U6 snRNA sequences. Here, we report on the characterisation of two novel mammalian C/D box snoRNAs, called U83 and U84, that contain all the elements that are essential for accumulation and function of 2′-O-methylation guide snoRNAs. However, in contrast to all of the known 2′-O-methylation guide RNAs, the human, mouse and pig U83 and U84 snoRNAs feature no antisense elements complementary to rRNA or U6 snRNA sequences. The human U83 and U84 snoRNAs are not associated with maturing nucleolar pre-ribosomal particles, suggesting that they do not function in rRNA biogenesis. Since artificial substrate RNAs complementary to the evolutionarily conserved putative substrate recognition motifs of the U83 and U84 snoRNAs were correctly 2′-O-methy-lated in the nucleolus of mouse cells, we suggest that the new snoRNAs act as 2′-O-methylation guides for cellular RNAs other then rRNAs and the U6 snRNA.     

The influences of PRG-1 on the expression of small RNAs and mRNAs

Background

In metazoans, Piwi-related Argonaute proteins play important roles in maintaining germline integrity and fertility and have been linked to a class of germline-enriched small RNAs termed piRNAs. Caenorhabditis elegans encodes two Piwi family proteins called PRG-1 and PRG-2, and PRG-1 interacts with the C. elegans piRNAs (21U-RNAs). Previous studies found that mutation of prg-1 causes a marked reduction in the expression of 21U-RNAs, temperature-sensitive defects in fertility and other phenotypic defects.

Results

In this study, we wanted to systematically demonstrate the function of PRG-1 in the regulation of small RNAs and their targets. By analyzing small RNAs and mRNAs with and without a mutation in prg-1 during C. elegans development, we demonstrated that (1) mutation of prg-1 leads to a decrease in the expression of 21U-RNAs, and causes 35 ~ 40% of miRNAs to be down-regulated; (2) in C. elegans, approximately 3% (6% in L4) of protein-coding genes are differentially expressed after mutating prg-1, and 60 ~ 70% of these substantially altered protein-coding genes are up-regulated; (3) the target genes of the down-regulated miRNAs and the candidate target genes of the down-regulated 21U-RNAs are enriched in the up-regulated protein-coding genes; and (4) PRG-1 regulates protein-coding genes by down-regulating small RNAs (miRNAs and 21U-RNAs) that target genes that participate in the development of C. elegans.

Conclusions

In prg-1-mutated C. elegans, the expression of miRNAs and 21U-RNAs was reduced, and the protein-coding targets, which were associated with the development of C. elegans, were up-regulated. This may be the mechanism underlying PRG-1 function.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-321) contains supplementary material, which is available to authorized users.     

Genome-wide analysis of single non-templated nucleotides in plant endogenous siRNAs and miRNAs

Plant small RNAs are subject to various modifications. Previous reports revealed widespread 3′ modifications (truncations and non-templated tailing) of plant miRNAs when the 2′-O-methyltransferase HEN1 is absent. However, non-templated nucleotides in plant heterochromatic siRNAs have not been deeply studied, especially in wild-type plants. We systematically studied non-templated nucleotide patterns in plant small RNAs by analyzing small RNA sequencing libraries from Arabidopsis, tomato, Medicago, rice, maize and Physcomitrella. Elevated rates of non-templated nucleotides were observed at the 3′ ends of both miRNAs and endogenous siRNAs from wild-type specimens of all species. ‘Off-sized’ small RNAs, such as 25 and 23 nt siRNAs arising from loci dominated by 24 nt siRNAs, often had very high rates of 3′-non-templated nucleotides. The same pattern was observed in all species that we studied. Further analysis of 24 nt siRNA clusters in Arabidopsis revealed distinct patterns of 3′-non-templated nucleotides of 23 nt siRNAs arising from heterochromatic siRNA loci. This pattern of non-templated 3′ nucleotides on 23 nt siRNAs is not affected by loss of known small RNA 3′-end modifying enzymes, and may result from modifications added to longer heterochromatic siRNA precursors.     

RNA-Seq of the Nucleolus Reveals Abundant SNORD44-Derived Small RNAs

Small non-coding RNAs represent RNA species that are not translated to proteins, but which have diverse and broad functional activities in physiological and pathophysiological states. The knowledge of these small RNAs is rapidly expanding in part through the use of massive parallel (deep) sequencing efforts. We present here the first deep sequencing of small RNomes in subcellular compartments with particular emphasis on small RNAs (sRNA) associated with the nucleolus. The vast majority of the cellular, cytoplasmic and nuclear sRNAs were identified as miRNAs. In contrast, the nucleolar sRNAs had a unique size distribution consisting of 19–20 and 25 nt RNAs, which were predominantly composed of small snoRNA-derived box C/D RNAs (termed as sdRNA). Sequences from 47 sdRNAs were identified, which mapped to both 5′ and 3′ ends of the snoRNAs, and retained conserved box C or D motifs. SdRNA reads mapping to SNORD44 comprised 74% of all nucleolar sdRNAs, and were confirmed by Northern blotting as comprising both 20 and 25 nt RNAs. A novel 120 nt SNORD44 form was also identified. The expression of the SNORD44 sdRNA and 120 nt form was independent of Dicer/Drosha–mediated processing pathways but was dependent on the box C/D snoRNP proteins/sno-ribonucleoproteins fibrillarin and NOP58. The 120 nt SNORD44-derived RNA bound to fibrillarin suggesting that C/D sno-ribonucleoproteins are involved in regulating the stability or processing of SNORD44. This study reveals sRNA cell-compartment specific expression and the distinctive unique composition of the nucleolar sRNAs.     

Common Sequence at the 5′ Ends of the Segmented RNA Genomes of Influenza A and B Viruses

Guanylyl- and methyltransferases, isolated from purified vaccinia virus, were used to specifically label the 5′ ends of the genome RNAs of influenza A and B viruses. All eight segments were labeled with [α-³²P]guanosine 5′-triphosphate or S-adenosyl[methyl-³H]methionine to form “cap” structures of the type m⁷G(5′)pppN^m-, of which unmethylated (p)ppN- represents the original 5′ end. Further analyses indicated that m⁷G(5′)pppA^m, m⁷G(5′)pppA^mpGp, and m⁷G(5′)pppA^mpGpUp were released from total and individual labeled RNA segments by digestion with nuclease P1, RNase T1, and RNase A, respectively. Consequently, the 5′-terminal sequences of most or all individual genome RNAs of influenza A and B viruses were deduced to be (p)ppApGpUp. The presence of identical sequences at the ends of RNA segments of both types of influenza viruses indicates that they have been specifically conserved during evolution.     

The Nematode Eukaryotic Translation Initiation Factor 4E/G Complex Works with a trans-Spliced Leader Stem-Loop To Enable Efficient Translation of Trimethylguanosine-Capped RNAs

Eukaryotic mRNA translation begins with recruitment of the 40S ribosome complex to the mRNA 5′ end through the eIF4F initiation complex binding to the 5′ m⁷G-mRNA cap. Spliced leader (SL) RNA trans splicing adds a trimethylguanosine (TMG) cap and a sequence, the SL, to the 5′ end of mRNAs. Efficient translation of TMG-capped mRNAs in nematodes requires the SL sequence. Here we define a core set of nucleotides and a stem-loop within the 22-nucleotide nematode SL that stimulate translation of mRNAs with a TMG cap. The structure and core nucleotides are conserved in other nematode SLs and correspond to regions of SL1 required for early Caenorhabditis elegans development. These SL elements do not facilitate translation of m⁷G-capped RNAs in nematodes or TMG-capped mRNAs in mammalian or plant translation systems. Similar stem-loop structures in phylogenetically diverse SLs are predicted. We show that the nematode eukaryotic translation initiation factor 4E/G (eIF4E/G) complex enables efficient translation of the TMG-SL RNAs in diverse in vitro translation systems. TMG-capped mRNA translation is determined by eIF4E/G interaction with the cap and the SL RNA, although the SL does not increase the affinity of eIF4E/G for capped RNA. These results suggest that the mRNA 5′ untranslated region (UTR) can play a positive and novel role in translation initiation through interaction with the eIF4E/G complex in nematodes and raise the issue of whether eIF4E/G-RNA interactions play a role in the translation of other eukaryotic mRNAs.Cap-dependent translation initiation in eukaryotes is a complex process involving many factors and serves as the primary mechanism for eukaryotic translation (37, 44). The first step in the initiation process, recruitment of the m⁷G (7-methylguanosine)-capped mRNA to the ribosome, is widely considered the rate-limiting step. It begins with recognition of and binding to the m⁷G cap at the 5′ end of the mRNA by the eukaryotic translation initiation factor 4F (eIF4F) complex, which contains three proteins: eIF4E (a cap-binding protein), eIF4G (a scaffold protein with RNA binding sites), and eIF4A (an RNA helicase). eIF4G''s interaction with eIF3, itself a multisubunit complex that interacts with the 40S ribosome, facilitates the actual recruitment of capped RNA to the ribosome. With the help of several other initiation factors, the small ribosomal subunit scans the mRNA from 5′ to 3′ until a translation initiation codon (AUG) in appropriate context is identified and an 80S ribosomal complex is formed, after which the first peptide bond is formed, thus ending the initiation process (37, 44). The AUG context can play an important role in the efficiency of translation initiation (23, 44). The length, structure, and presence of AUGs or open reading frames in the mRNA 5′ untranslated region (UTR) can negatively affect cap-dependent translation and ribosomal scanning. In general, long and highly structured 5′ UTRs, as well as upstream AUGs leading to short open reading frames, can impede ribosome scanning and lead to reduced translation (23, 44). In addition, 5′ UTRs less than 10 nucleotides (nt) in length are thought to be too short to enable preinitiation complex assembly and scanning (24). Thus, several attributes of the mRNA 5′ UTR are known to negatively affect translation initiation, whereas only the AUG context and the absence of negative elements are known to have a positive effect on translation initiation (44).Two of the important mRNA features associated with cap-dependent translation, the cap and the 5′ UTR, are significantly altered by an RNA processing event known as spliced leader (SL) trans splicing (3, 8, 17, 26, 36, 47). This takes place in members of a diverse group of eukaryotic organisms, including some protozoa, sponges, cnidarians, chaetognaths, flatworms, nematodes, rotifers, crustaceans, and tunicates (17, 28, 39, 55, 56). In SL trans splicing, a separately transcribed small exon (16 to 51 nucleotides [nt]) with its own cap gets added to the 5′ end of pre-mRNAs. This produces mature mRNAs with a unique cap and a conserved sequence in the 5′ UTR. In metazoa, the m⁷G cap is replaced with a trimethylguanosine (TMG) cap (m^2,2,7GpppN) (27, 30, 46, 49). In nematodes, ∼70% of all mRNAs are trans spliced and therefore have a TMG cap and an SL (2). In general, eukaryotic eIF4E proteins do not effectively recognize the TMG cap (35). This raises the issues of how the translation machinery in trans-splicing metazoa effectively recognizes TMG-capped trans-spliced mRNAs, what role the SL sequence plays in translation initiation, and how the conserved translation initiation machinery has adapted to effectively translate trans-spliced mRNAs.Previous work has shown that efficient translation of TMG-capped messages in nematodes requires the SL sequence (22 nt) immediately downstream of the cap (5, 25, 29). In the current studies, we sought to understand the manner in which the SL enhanced the translation of TMG-capped mRNAs. Using a cell-free nematode in vitro translation system, we carried out mutational analyses that define the specific sequences in the SL that are required and sufficient for efficient translation of TMG-capped mRNAs. These analyses led to the discovery of a small, discrete stem-loop immediately adjacent to the TMG cap in trans-spliced messages required for efficient translation. Notably, the sequences involved in the base pairing of the stem are highly conserved in alternative SL sequences found in nematodes. We further show that the nematode eIF4E/G complex plays a major role in facilitating the SL enhancement of TMG-capped mRNA that likely occurs after the initial cap-binding step. The results demonstrate the importance of specific enhancing elements in the 5′ UTR and adaptation in the eIF4F complex necessary for optimal cap-dependent translation.     

Dicer independent small RNAs associate with telomeric heterochromatin

Small RNAs play important roles in the establishment and maintenance of heterochromatin structures. We show the presence of telomere specific small RNAs (tel-sRNAs) in mouse embryonic stem cells that are ∼24 nucleotides in length, Dicer-independent, and 2′-O-methylated at the 3′ terminus. The tel-sRNAs are asymmetric with specificity toward telomere G-rich strand, and evolutionarily conserved from protozoan to mammalian cells. Furthermore, tel-sRNAs are up-regulated in cells that carry null mutation of H3K4 methyltransferase MLL (Mll^(−/−)) and down-regulated in cells that carry null mutations of histone H3K9 methyltransferase SUV39H (Suv39h1/h2^(−/−)), suggesting that they are subject to epigenetic regulation. These results support that tel-sRNAs are heterochromatin associated pi-like small RNAs.