High-efficiency RNA cloning enables accurate quantification of miRNA expression by deep sequencing

Small RNA cloning and sequencing is uniquely positioned as a genome-wide approach to quantify miRNAs with single-nucleotide resolution. However, significant biases introduced by RNA ligation in current protocols lead to inaccurate miRNA quantification by 1000-fold. Here we report an RNA cloning method that achieves over 95% efficiency for both 5′ and 3′ ligations. It achieves accurate quantification of synthetic miRNAs with less than two-fold deviation from the anticipated value and over a dynamic range of four orders of magnitude. Taken together, this high-efficiency RNA cloning method permits accurate genome-wide miRNA profiling from total RNAs.     

Structural bias in T4 RNA ligase-mediated 3'-adapter ligation

T4 RNA ligases are commonly used to attach adapters to RNAs, but large differences in ligation efficiency make detection and quantitation problematic. We developed a ligation selection strategy using random RNAs in combination with high-throughput sequencing to gain insight into the differences in efficiency of ligating pre-adenylated DNA adapters to RNA 3'-ends. After analyzing biases in RNA sequence, secondary structure and RNA-adapter cofold structure, we conclude that T4 RNA ligases do not show significant primary sequence preference in RNA substrates, but are biased against structural features within RNAs and adapters. Specifically, RNAs with less than three unstructured nucleotides at the 3'-end and RNAs that are predicted to cofold with an adapter in unfavorable structures are likely to be poorly ligated. The effect of RNA-adapter cofold structures on ligation is supported by experiments where the ligation efficiency of specific miRNAs was changed by designing adapters to alter cofold structure. In addition, we show that using adapters with randomized regions results in higher ligation efficiency and reduced ligation bias. We propose that using randomized adapters may improve RNA representation in experiments that include a 3'-adapter ligation step.     

Argonaute loading improves the 5' precision of both MicroRNAs and their miRNA* strands in flies

MicroRNAs (miRNAs) are short regulatory RNAs that direct repression of their mRNA targets. The miRNA "seed"-nucleotides 2-7-establishes target specificity by mediating target binding. Accurate processing of the miRNA 5' end is thought to be under strong selective pressure because a shift by just one nucleotide in the 5' end of a miRNA alters its seed sequence, redefining its repertoire of targets (Figure 1). Animal miRNAs are produced by the sequential cleavage of partially double-stranded precursors by the RNase III endonucleases Drosha and Dicer, thereby generating a transitory double-stranded intermediate comprising the miRNA paired to its partially complementary miRNA strand. Here, we report that in flies, the 5' ends of miRNAs and miRNA strands are typically more precisely defined than their 3' ends. Surprisingly, the precision of the 5' ends of both miRNA and miRNA sequences increases after Argonaute2 (Ago2) loading. Our data imply that either many miRNA sequences are under evolutionary pressure to maintain their seed sequences-that is, they have targets-or that secondary constraints, such as the sequence requirements for loading small RNAs into functional Argonaute complexes, narrow the range of miRNA and miRNA 5' ends that accumulate in flies.     

Amplified microRNA detection by templated chemistry

MicroRNAs (miRNAs) are a class of RNAs that play important regulatory roles in the cell. The detection of microRNA has attracted significant interest recently, as abnormal miRNA expression has been linked to cancer and other diseases. Here, we present a straightforward method for isothermal amplified detection of miRNA that involves two separate nucleic acid-templated chemistry steps. The miRNA first templates the cyclization of an oligodeoxynucleotide from a linear precursor containing a 5'-iodide and a 3'-phosphorothioate. The sequence is amplified through rolling circle amplification with 29 DNA polymerase and then detected via a second amplification using fluorogenic templated probes. Tests showed that the cyclization proceeds in ～50% yield over 24 h and is compatible with the conditions required for rolling circle polymerization, unlike enzymatic ligations which required non-compatible buffer conditions. The polymerization yielded 188-fold amplification, and separate experiments showed ～15-fold signal amplification from the templated fluorogenic probes. When all components are combined, results show miRNA detection down to 200 pM in solution, and correlation of the detected signal with the initial concentration of miRNA. The doubly templated double-amplification method demonstrates a new approach to detection of rolling circle products and significant advantages in ease of operation for miRNA detection.     

Identification and characterization of microRNAs from the bovine adipose tissue and mammary gland

We report the identification of bovine miRNAs by cloning small RNAs from adipose tissue and the mammary gland. Fifty-nine distinct miRNAs were identified, five of them were not homologous to known mammalian miRNAs, and many of them had 3' and/or 5' end variants. Ribonuclease protection assays indicated that miR-23a and miR-24, whose genes are closely located on the same chromosome, were co-expressed in different tissues. The assays also suggested a role for several miRNAs in the mammary gland and a role for miR-133, a previously known skeletal and cardiac muscle-specific miRNA, in the rumen, an organ unique to the ruminant.     

Improved annotation of C. elegans microRNAs by deep sequencing reveals structures associated with processing by Drosha and Dicer

MicroRNAs (miRNAs) are small regulatory RNAs that are essential in all studied metazoans. Research has focused on the prediction and identification of novel miRNAs, while little has been done to validate, annotate, and characterize identified miRNAs. Using Illumina sequencing, ～20 million small RNA sequences were obtained from Caenorhabditis elegans. Of the 175 miRNAs listed on the miRBase database, 106 were validated as deriving from a stem-loop precursor with hallmark characteristics of miRNAs. This result suggests that not all sequences identified as miRNAs belong in this category of small RNAs. Our large data set of validated miRNAs facilitated the determination of general sequence and structural characteristics of miRNAs and miRNA precursors. In contrast to previous observations, we did not observe a preference for the 5' nucleotide of the miRNA to be unpaired compared to the 5' nucleotide of the miRNA*, nor a preference for the miRNA to be on either the 5' or 3' arm of the miRNA precursor stem-loop. We observed that steady-state pools of miRNAs have fairly homogeneous termini, especially at their 5' end. Nearly all mature miRNA-miRNA* duplexes had two nucleotide 3' overhangs, and there was a preference for a uracil in the first and ninth position of the mature miRNA. Finally, we observed that specific nucleotides and structural distortions were overrepresented at certain positions adjacent to Drosha and Dicer cleavage sites. Our study offers a comprehensive data set of C. elegans miRNAs and their precursors that significantly decreases the uncertainty associated with the identity of these molecules in existing databases.     

Next-generation sequencing of the Chinese hamster ovary microRNA transcriptome: Identification, annotation and profiling of microRNAs as targets for cellular engineering

Chinese hamster ovary (CHO) cells are the predominant cell factory for the production of recombinant therapeutic proteins. Nevertheless, the lack in publicly available sequence information is severely limiting advances in CHO cell biology, including the exploration of microRNAs (miRNA) as tools for CHO cell characterization and engineering. In an effort to identify and annotate both conserved and novel CHO miRNAs in the absence of a Chinese hamster genome, we deep-sequenced small RNA fractions of 6 biotechnologically relevant cell lines and mapped the resulting reads to an artificial reference sequence consisting of all known miRNA hairpins. Read alignment patterns and read count ratios of 5' and 3' mature miRNAs were obtained and used for an independent classification into miR/miR* and 5p/3p miRNA pairs and discrimination of miRNAs from other non-coding RNAs, resulting in the annotation of 387 mature CHO miRNAs. The quantitative content of next-generation sequencing data was analyzed and confirmed using qPCR, to find that miRNAs are markers of cell status. Finally, cDNA sequencing of 26 validated targets of miR-17-92 suggests conserved functions for miRNAs in CHO cells, which together with the now publicly available sequence information sets the stage for developing novel RNAi tools for CHO cell engineering.     

Methylation protects miRNAs and siRNAs from a 3'-end uridylation activity in Arabidopsis

Small RNAs of 21-25 nucleotides (nt), including small interfering RNAs (siRNAs) and microRNAs (miRNAs), act as guide RNAs to silence target-gene expression in a sequence-specific manner. In addition to a Dicer homolog, DCL1, the biogenesis of miRNAs in Arabidopsis requires another protein, HEN1. miRNAs are reduced in abundance and increased in size in hen1 mutants. We found that HEN1 is a miRNA methyltransferase that adds a methyl group to the 3'-most nucleotide of miRNAs, but the role of miRNA methylation was unknown. Here, we show that siRNAs from sense transgenes, hairpin transgenes, and transposons or repeat sequences, as well as a new class of siRNAs known as trans-acting siRNAs, are also methylated in vivo by HEN1. In addition, we show that the size increase of small RNAs in the hen1-1 mutant is due to the addition of one to five U residues to the 3' ends of the small RNAs. Therefore, a novel uridylation activity targets the 3' ends of unmethylated miRNAs and siRNAs in hen1 mutants. We conclude that 3'-end methylation is a common step in miRNA and siRNA metabolism and likely protects the 3' ends of the small RNAs from the uridylation activity.     

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.     

Global analysis of non-coding small RNAs in Arabidopsis in response to jasmonate treatment by deep sequencing technology

In plants, non-coding small RNAs play a vital role in plant development and stress responses. To explore the possible role of non-coding small RNAs in the regulation of the jasmonate (JA) pathway, we compared the non-coding small RNAs between the JA-deficient aos mutant and the JA-treated wild type Arabidopsis via high-throughput sequencing. Thirty new miRNAs and 27 new miRNA candidates were identified through bioinformatics approach. Forty-nine known miRNAs (belonging to 24 families), 15 new miRNAs and new miRNA candidates (belonging to 11 families) and 3 tasiRNA families were induced by JA, whereas 1 new miRNA, 1 tasiRNA family and 22 known miRNAs (belonging to 9 families) were repressed by JA.     

Complexity of murine cardiomyocyte miRNA biogenesis, sequence variant expression and function

microRNAs (miRNAs) are critical to heart development and disease. Emerging research indicates that regulated precursor processing can give rise to an unexpected diversity of miRNA variants. We subjected small RNA from murine HL-1 cardiomyocyte cells to next generation sequencing to investigate the relevance of such diversity to cardiac biology. ～40 million tags were mapped to known miRNA hairpin sequences as deposited in miRBase version 16, calling 403 generic miRNAs as appreciably expressed. Hairpin arm bias broadly agreed with miRBase annotation, although 44 miR* were unexpectedly abundant (>20% of tags); conversely, 33 -5p/-3p annotated hairpins were asymmetrically expressed. Overall, variability was infrequent at the 5' start but common at the 3' end of miRNAs (5.2% and 52.3% of tags, respectively). Nevertheless, 105 miRNAs showed marked 5' isomiR expression (>20% of tags). Among these was miR-133a, a miRNA with important cardiac functions, and we demonstrated differential mRNA targeting by two of its prevalent 5' isomiRs. Analyses of miRNA termini and base-pairing patterns around Drosha and Dicer cleavage regions confirmed the known bias towards uridine at the 5' most position of miRNAs, as well as supporting the thermodynamic asymmetry rule for miRNA strand selection and a role for local structural distortions in fine tuning miRNA processing. We further recorded appreciable expression of 5 novel miR*, 38 extreme variants and 8 antisense miRNAs. Analysis of genome-mapped tags revealed 147 novel candidate miRNAs. In summary, we revealed pronounced sequence diversity among cardiomyocyte miRNAs, knowledge of which will underpin future research into the mechanisms involved in miRNA biogenesis and, importantly, cardiac function, disease and therapy.     

Recapitulation of short RNA-directed translational gene silencing in vitro

microRNAs (miRNAs) are a large class of endogenous short RNAs that repress gene expression. Many miRNAs are conserved throughout evolution, and dysregulation of miRNA pathways has been correlated with an increasing number of human diseases. In animals, miRNAs typically bind to the 3' untranslated region (3'UTR) of target mRNAs with imperfect sequence complementarity and repress translation. Despite their importance in regulating biological processes in numerous organisms, the mechanisms of miRNA function are largely unknown. Here, we report in vitro reactions for miRNA-directed translational gene silencing. These reactions faithfully recapitulate known in vivo hallmarks of mammalian miRNA function, including a requirement for a 5' phosphate and perfect complementarity to the mRNA target in the 5' seed region. Translational gene silencing by miRNAs in vitro requires target mRNAs to possess a 7-methyl G cap and a polyA tail, whereas increasing polyA tail length alone can increase miRNA silencing activity.     

Functional diversity of microRNA variants in plants

MicroRNAs (miRNAs) are small noncoding RNAs that play a crucial role in plant growth, development, and stress responses by regulating target gene expression. With the development of high-throughput sequencing technology, it has been facilitated to identify new miRNAs, as well as diversity and variability of miRNA variants. miRNA variants share the sequences with other closely related miRNAs and contain length and/or sequence variations at the 5’-, 3’-ends, as well as internal positions. They originate from the same miRNA precursor or from the diversity of members in the same miRNA family. It has been reported that tissue- or condition-specific variation in the relative abundance of different miRNA variants could contribute to differential functions of those in development or stress responses. In addition, the diversity of miRNA variants affects stability, loading efficiency onto AGO, and target selection. In this review, the diversification of miRNA sequences and the evidences of the distinct functional role of miRNA variants will be discussed.