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.     

Regulation of small RNA stability: methylation and beyond

As central components of RNA silencing, small RNAs play diverse and important roles in many biological processes in eukaryotes. Aberrant reduction or elevation in the levels of small RNAs is associated with many developmental and physiological defects. The in vivo levels of small RNAs are precisely regulated through modulating the rates of their biogenesis and turnover. 2′-O-methylation on the 3′ terminal ribose is a major mechanism that increases the stability of small RNAs. The small RNA methyltransferase HUA ENHANCER1 (HEN1) and its homologs methylate microRNAs and small interfering RNAs (siRNAs) in plants, Piwi-interacting RNAs (piRNAs) in animals, and siRNAs in Drosophila. 3′ nucleotide addition, especially uridylation, and 3′-5′ exonucleolytic degradation are major mechanisms that turnover small RNAs. Other mechanisms impacting small RNA stability include complementary RNAs, cis-elements in small RNA sequences and RNA-binding proteins. Investigations are ongoing to further understand how small RNA stability impacts their accumulation in vivo in order to improve the utilization of RNA silencing in biotechnology and 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.     

Dumbbell-PCR: a method to quantify specific small RNA variants with a single nucleotide resolution at terminal sequences

Recent advances in next-generation sequencing technologies have revealed that cellular functional RNAs are not always expressed as single entities with fixed terminal sequences but as multiple isoforms bearing complex heterogeneity in both length and terminal sequences, such as isomiRs, the isoforms of microRNAs. Unraveling the biogenesis and biological significance of heterogenetic RNA expression requires distinctive analysis of each RNA variant. Here, we report the development of dumbbell PCR (Db-PCR), an efficient and convenient method to distinctively quantify a specific individual small RNA variant. In Db-PCR, 5′- and 3′-stem–loop adapters are specifically hybridized and ligated to the 5′- and 3′-ends of target RNAs, respectively, by T4 RNA ligase 2 (Rnl2). The resultant ligation products with ‘dumbbell-like’ structures are subsequently quantified by TaqMan RT-PCR. We confirmed that high specificity of Rnl2 ligation and TaqMan RT-PCR toward target RNAs assured both 5′- and 3′-terminal sequences of target RNAs with single nucleotide resolution so that Db-PCR specifically detected target RNAs but not their corresponding terminal variants. Db-PCR had broad applicability for the quantification of various small RNAs in different cell types, and the results were consistent with those from other quantification method. Therefore, Db-PCR provides a much-needed simple method for analyzing RNA terminal heterogeneity.     

Distinct substrate specificities of Arabidopsis DCL3 and DCL4

In Arabidopsis thaliana, Dicer-like 3 (DCL3) and Dicer-like 4 (DCL4) cleave long, perfect double-stranded RNAs (dsRNAs) into 24 and 21 nucleotides (nt) small interfering RNAs, respectively, which in turn function in RNA-directed DNA methylation and RNA interference, respectively. To reveal how DCL3 and DCL4 individually recognize long perfect dsRNAs as substrates, we biochemically characterized DCL3 and DCL4 and compared their enzymatic properties. DCL3 preferentially cleaves short dsRNAs with 5′ phosphorylated adenosine or uridine and a 1 nt 3′ overhang, whereas DCL4 cleaves long dsRNAs with blunt ends or with a 1 or 2 nt 3′ overhang with similar efficiency. DCL3 produces 24 nt RNA duplexes with 2 nt 3′ overhangs by the 5′ counting rule. Inorganic phosphate, NaCl and KCl enhance DCL3 activity but inhibit DCL4 activity. These results indicate that plants use DCLs with distinct catalytic profiles to ensure each dsRNA substrate generates only a specific length of siRNAs that trigger a unique siRNA-mediated response.     

The 5′ Terminal Trailer Region of Vesicular Stomatitis Virus Contains a Position-Dependent cis-Acting Signal for Assembly of RNA into Infectious Particles

The cis-acting genomic RNA requirements for the assembly of vesicular stomatitis virus (VSV) ribonucleocapsids into infectious particles were investigated. Using a biological assay based on particle infectivity, we demonstrated that subgenomic replicons that contained all four possible combinations of the natural genomic termini, the 3′ leader (Le) and 5′ trailer (Tr) regions, were replication competent; however, a 3′ copyback replicon (3′CB), containing the natural 3′ terminus but having the 5′ Tr replaced by a sequence complementary to the 3′ Le for 46 nucleotides, was unable to assemble infectious particles, despite efficient replication. When a copy of Tr was inserted 51 nucleotides from the 5′ end of 3′CB, infectious particles were produced. However, analysis of the replication products of these particles showed that the 51 nucleotides which corresponded to the Le complement sequences at the 5′ terminus were removed during RNA replication, thus restoring the wild-type 5′ Tr to the exact 5′ terminus. These data showed that a cis-acting signal was necessary for assembly of VSV RNAs into infectious particles and that this signal was supplied by Tr when located at the 5′ end. The regions within Tr required for assembly were analyzed by a series of deletions and exchanges for Le complement sequences, which demonstrated that the 5′ terminal 29 nucleotides of Tr allowed assembly of infectious particles but that the 5′ terminal 22 nucleotides functioned poorly. Deletions in Tr also altered the balance between negative- and positive-strand genomic RNA and affected levels of replication. RNAs that retained fewer than 45 but at least 22 nucleotides of the 5′ terminus could replicate but were impaired in RNA replication, and RNAs that retained only 14 nucleotides of the 5′ terminus were severely reduced in ability to replicate. These data define the VSV Tr as a position-dependent, cis-acting element for the assembly of RNAs into infectious particles, and they delineate RNA sequences that are essential for negative-strand RNA synthesis. These observations are consistent with, and offer an explanation for, the absence of 3′ copyback defective interfering particles in nature.     

Genome 3'-end repair in dengue virus type 2

Genomes of RNA viruses encounter a continual threat from host cellular ribonucleases. Therefore, viruses have evolved mechanisms to protect the integrity of their genomes. To study the mechanism of 3′-end repair in dengue virus-2 in mammalian cells, a series of 3′-end deletions in the genome were evaluated for virus replication by detection of viral antigen NS1 and by sequence analysis. Limited deletions did not cause any delay in the detection of NS1 within 5 d. However, deletions of 7–10 nucleotides caused a delay of 9 d in the detection of NS1. Sequence analysis of RNAs from recovered viruses showed that at early times, virus progenies evolved through RNA molecules of heterogeneous lengths and nucleotide sequences at the 3′ end, suggesting a possible role for terminal nucleotidyl transferase activity of the viral polymerase (NS5). However, this diversity gradually diminished and consensus sequences emerged. Template activities of 3′-end mutants in the synthesis of negative-strand RNA in vitro by purified NS5 correlate well with the abilities of mutant RNAs to repair and produce virus progenies. Using the Mfold program for RNA structure prediction, we show that if the 3′ stem–loop (3′ SL) structure was abrogated by mutations, viruses eventually restored the 3′ SL structure. Taken together, these results favor a two-step repair process: non-template-based nucleotide addition followed by evolutionary selection of 3′-end sequences based on the best-fit RNA structure that can support viral replication.     

Aberration-corrected ultrafine analysis of miRNA reads at single-base resolution: a k-mer lattice approach

Raw sequencing reads of miRNAs contain machine-made substitution errors, or even insertions and deletions (indels). Although the error rate can be low at 0.1%, precise rectification of these errors is critically important because isoform variation analysis at single-base resolution such as novel isomiR discovery, editing events understanding, differential expression analysis, or tissue-specific isoform identification is very sensitive to base positions and copy counts of the reads. Existing error correction methods do not work for miRNA sequencing data attributed to miRNAs’ length and per-read-coverage properties distinct from DNA or mRNA sequencing reads. We present a novel lattice structure combining kmers, (k – 1)mers and (k + 1)mers to address this problem. The method is particularly effective for the correction of indel errors. Extensive tests on datasets having known ground truth of errors demonstrate that the method is able to remove almost all of the errors, without introducing any new error, to improve the data quality from every-50-reads containing one error to every-1300-reads containing one error. Studies on experimental miRNA sequencing datasets show that the errors are often rectified at the 5′ ends and the seed regions of the reads, and that there are remarkable changes after the correction in miRNA isoform abundance, volume of singleton reads, overall entropy, isomiR families, tissue-specific miRNAs, and rare-miRNA quantities.     

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.     

Deep Sequencing of the Small RNAs Derived from Two Symptomatic Variants of a Chloroplastic Viroid: Implications for Their Genesis and for Pathogenesis

Northern-blot hybridization and low-scale sequencing have revealed that plants infected by viroids, non-protein-coding RNA replicons, accumulate 21–24 nt viroid-derived small RNAs (vd-sRNAs) similar to the small interfering RNAs, the hallmarks of RNA silencing. These results strongly support that viroids are elicitors and targets of the RNA silencing machinery of their hosts. Low-scale sequencing, however, retrieves partial datasets and may lead to biased interpretations. To overcome this restraint we have examined by deep sequencing (Solexa-Illumina) and computational approaches the vd-sRNAs accumulating in GF-305 peach seedlings infected by two molecular variants of Peach latent mosaic viroid (PLMVd) inciting peach calico (albinism) and peach mosaic. Our results show in both samples multiple PLMVd-sRNAs, with prevalent 21-nt (+) and (−) RNAs presenting a biased distribution of their 5′ nucleotide, and adopting a hotspot profile along the genomic (+) and (−) RNAs. Dicer-like 4 and 2 (DCL4 and DCL2, respectively), which act hierarchically in antiviral defense, likely also mediate the genesis of the 21- and 22-nt PLMVd-sRNAs. More specifically, because PLMVd replicates in plastids wherein RNA silencing has not been reported, DCL4 and DCL2 should dice the PLMVd genomic RNAs during their cytoplasmic movement or the PLMVd-dsRNAs generated by a cytoplasmic RNA-dependent RNA polymerase (RDR), like RDR6, acting in concert with DCL4 processing. Furthermore, given that vd-sRNAs derived from the 12–14-nt insertion containing the pathogenicity determinant of peach calico are underrepresented, it is unlikely that symptoms may result from the accidental targeting of host mRNAs by vd-sRNAs from this determinant guiding the RNA silencing machinery.     

Selenium derivatization and crystallization of DNA and RNA oligonucleotides for X-ray crystallography using multiple anomalous dispersion

We report here the solid phase synthesis of RNA and DNA oligonucleotides containing the 2′-selenium functionality for X-ray crystallography using multiwavelength anomalous dispersion. We have synthesized the novel 2′-methylseleno cytidine phosphoramidite and improved the accessibility of the 2′-methylseleno uridine phosphoramidite for the synthesis of many selenium-derivatized DNAs and RNAs in large scales. The yields of coupling these Se-nucleoside phosphoramidites into DNA or RNA oligonucleotides were over 99% when 5-(benzylmercapto)-1H-tetrazole was used as the coupling reagent. The UV melting study of A-form dsDNAs indicated that the 2′-selenium derivatization had no effect on the stability of the duplexes with the 3′-endo sugar pucker. Thus, the stems of functional RNA molecules with the same 3′-endo sugar pucker appear to be the ideal sites for the selenium derivatization with 2′-Se-C and 2′-Se-U. Crystallization of the selenium-derivatized oligonucleotides is also reported here. The results demonstrate that this 2′-selenium functionality is suitable for RNA and A-form DNA derivatization in X-ray crystallography.