Reproducible features of small RNAs in C. elegans reveal NU RNAs and provide insights into 22G RNAs and 26G RNAs

Small RNAs regulate gene expression and most genes in the worm Caenorhabditis elegans are subject to their regulation. Here, we analyze small RNA data sets and use reproducible features of RNAs present in multiple data sets to discover a new class of small RNAs and to reveal insights into two known classes of small RNAs—22G RNAs and 26G RNAs. We found that reproducibly detected 22-nt RNAs, although are predominantly RNAs with a G at the 5′ end, also include RNAs with A, C, or U at the 5′ end. These RNAs are synthesized downstream from characteristic sequence motifs on mRNA and have U-tailed derivatives. Analysis of 26G RNAs revealed that they are processed from a blunt end of double-stranded RNAs and that production of one 26G RNA generates a hotspot immediately downstream for production of another. To our surprise, analysis of RNAs shorter than 18 nt revealed a new class of RNAs, which we call NU RNAs (pronounced “new RNAs”) because they have a NU bias at the 5′ end, where N is any nucleotide. NU RNAs are antisense to genes and originate downstream from U bases on mRNA. Although many genes have complementary NU RNAs, their genome-wide distribution is distinct from that of previously known classes of small RNAs. Our results suggest that current approaches underestimate reproducibly detected RNAs that are shorter than 18 nt, and theoretical considerations suggest that such shorter RNAs could be used for sequence-specific gene regulation in organisms like C. elegans that have small genomes.     

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.     

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.     

Identification and characterization of endogenous small interfering RNAs from rice

RNA silencing-mediated small interfering RNAs (siRNAs) and microRNAs (miRNAs) have diverse natural roles, ranging from regulation of gene expression and heterochromatin formation to genome defense against transposons and viruses. Unlike miRNAs, endogenous siRNAs are generally not conserved between species; consequently, their identification requires experimental approaches. Thus far, endogenous siRNAs have not been reported from rice, which is a model species for monocotyledonous plants. We identified a large set of putative endogenous siRNAs from root, shoot and inflorescence small RNA cDNA libraries of rice. Most of these siRNAs are from intergenic regions, although a substantial proportion (22%) originates from the introns and exons of protein-coding genes. Northern and RT–PCR analysis revealed that the expression of some of the siRNAs is tissue specific or developmental stage specific. A total of 25 transposons and 21 protein-coding genes were predicted to be cis-targets of some of the siRNAs. Based on sequence homology, we also predicted 111 putative trans-targets for 44 of the siRNAs. Interestingly, ~46% of the predicted trans-targets are transposable elements, which suggests that endogenous siRNAs may play an important role in the suppression of transposon proliferation. Using RNA ligase-mediated-5′ rapid amplification of cDNA end assays, we validated three of the predicted targets and provided evidence for both cis- and trans-silencing of target genes by siRNAs-guided mRNA cleavage.     

Genome-Wide Double-Stranded RNA Sequencing Reveals the Functional Significance of Base-Paired RNAs in Arabidopsis

The functional structure of all biologically active molecules is dependent on intra- and inter-molecular interactions. This is especially evident for RNA molecules whose functionality, maturation, and regulation require formation of correct secondary structure through encoded base-pairing interactions. Unfortunately, intra- and inter-molecular base-pairing information is lacking for most RNAs. Here, we marry classical nuclease-based structure mapping techniques with high-throughput sequencing technology to interrogate all base-paired RNA in Arabidopsis thaliana and identify ∼200 new small (sm)RNA–producing substrates of RNA–DEPENDENT RNA POLYMERASE6. Our comprehensive analysis of paired RNAs reveals conserved functionality within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs, as well as a novel population of functional RNAs, many of which are the precursors of smRNAs. Finally, we identify intra-molecular base-pairing interactions to produce a genome-wide collection of RNA secondary structure models. Although our methodology reveals the pairing status of RNA molecules in the absence of cellular proteins, previous studies have demonstrated that structural information obtained for RNAs in solution accurately reflects their structure in ribonucleoprotein complexes. Furthermore, our identification of RNA–DEPENDENT RNA POLYMERASE6 substrates and conserved functional RNA domains within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs using this approach strongly suggests that RNA molecules are correctly folded into their secondary structure in solution. Overall, our findings highlight the importance of base-paired RNAs in eukaryotes and present an approach that should be widely applicable for the analysis of this key structural feature of RNA.     

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.     

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.     

Intracellular mRNA cleavage by 3' tRNase under the direction of 2'-O-methyl RNA heptamers

Mammalian tRNA 3′ processing endoribonuclease (3′-tRNase) can cleave any RNA at any site under the direction of small guide RNA (sgRNA) in vitro. sgRNAs can be as short as heptamers, which are much smaller than small interfering RNAs of ~21 nt. Together with such flexibility in substrate recognition, the ubiquity and the constitutive expression of 3′-tRNase have suggested that this enzyme can be utilized for specific cleavage of cellular RNAs by introducing appropriate sgRNAs into living cells. Here we demonstrated that the expression of chloramphenicol acetyltransferase can be downregulated by an appropriate sgRNA which is introduced into Madin–Darby canine kidney epithelial cells as an expression plasmid or a synthetic 2′-O-methyl RNA. We also showed that 2′-O-methyl RNA heptamers can attack luciferase mRNAs with a high specificity and induce 3′-tRNase-mediated knock-down of the mRNAs in 293 cells. Furthermore, the MTT cell viability assay suggested that an RNA heptamer can downregulate the endogenous Bcl-2 mRNA in Sarcoma 180 cells. This novel sgRNA/3′-tRNase strategy for destroying specific cellular RNAs may be utilized for therapeutic applications.     

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.     

Identification and characterization of extensive intra-molecular associations between 3'-UTRs and their ORFs

During eukaryotic translation, mRNAs may form intra-molecular interactions between distant domains. The 5′-cap and the polyA tail were shown to interact through their associated proteins, and this can induce physical compaction of the mRNA in vitro. However, the stability of this intra-molecular association in translating mRNAs and whether additional contacts exist in vivo are largely unknown. To explore this, we applied a novel approach in which several endogenous polysomal mRNAs from Saccharomyces cerevisiae were cleaved near their stop codon and the resulting 3′-UTR fragments were tested either for co-sedimentation or co-immunoprecipitation (co-IP) with their ORFs. In all cases a significant fraction of the 3′-UTR fragments sedimented similarly to their ORF-containing fragments, yet the extent of co-sedimentation differed between mRNAs. Similar observations were obtained by a co-IP assay. Interestingly, various treatments that are expected to interfere with the cap to polyA interactions had no effect on the co-sedimentation pattern. Moreover, the 3′-UTR appeared to co-sediment with different regions from within the ORF. Taken together, these results indicate extensive physical associations between 3′-UTRs and their ORFs that vary between genes. This implies that polyribosomal mRNAs are in a compact configuration in vivo.     

Evolutionary Modeling and Prediction of Non-Coding RNAs in Drosophila

We performed benchmarks of phylogenetic grammar-based ncRNA gene prediction, experimenting with eight different models of structural evolution and two different programs for genome alignment. We evaluated our models using alignments of twelve Drosophila genomes. We find that ncRNA prediction performance can vary greatly between different gene predictors and subfamilies of ncRNA gene. Our estimates for false positive rates are based on simulations which preserve local islands of conservation; using these simulations, we predict a higher rate of false positives than previous computational ncRNA screens have reported. Using one of the tested prediction grammars, we provide an updated set of ncRNA predictions for D. melanogaster and compare them to previously-published predictions and experimental data. Many of our predictions show correlations with protein-coding genes. We found significant depletion of intergenic predictions near the 3′ end of coding regions and furthermore depletion of predictions in the first intron of protein-coding genes. Some of our predictions are colocated with larger putative unannotated genes: for example, 17 of our predictions showing homology to the RFAM family snoR28 appear in a tandem array on the X chromosome; the 4.5 Kbp spanned by the predicted tandem array is contained within a FlyBase-annotated cDNA.     

Modified ribosome profiling reveals high abundance of ribosome protected mRNA fragments derived from 3′ untranslated regions

Ribosome profiling identifies ribosome positions on translated mRNAs. A prominent feature of published datasets is the near complete absence of ribosomes in 3′ untranslated regions (3′UTR) although substantial ribosome density can be observed on non-coding RNAs. Here we perform ribosome profiling in cultured Drosophila and human cells and show that different features of translation are revealed depending on the nuclease and the digestion conditions used. Most importantly, we observe high abundance of ribosome protected fragments in 3′UTRs of thousands of genes without manipulation of translation termination. Affinity purification of ribosomes indicates that the 3′UTR reads originate from ribosome protected fragments. Association of ribosomes with the 3′UTR may be due to ribosome migration through the stop codon or 3′UTR mRNA binding to ribosomes on the coding sequence. This association depends primarily on the relative length of the 3′UTR and may be related to translational regulation or ribosome recycling, for which the efficiency is known to inversely correlate with 3′UTR length. Together our results indicate that ribosome profiling is highly dependent on digestion conditions and that ribosomes commonly associate with the 3′UTR, which may have a role in translational regulation.