Functional siRNAs and miRNAs exhibit strand bias

Both microRNAs (miRNA) and small interfering RNAs (siRNA) share a common set of cellular proteins (Dicer and the RNA-induced silencing complex [RISC]) to elicit RNA interference. In the following work, a statistical analysis of the internal stability of published miRNA sequences in the context of miRNA precursor hairpins revealed enhanced flexibility of miRNA precursors, especially at the 5'-anti-sense (AS) terminal base pair. The same trend was observed in siRNA, with functional duplexes displaying a lower internal stability (Delta0.5 kcal/mol) at the 5'-AS end than nonfunctional duplexes. Average internal stability of siRNA molecules retrieved from plant cells after introduction of long RNA sequences also shows this characteristic thermodynamic signature. Together, these results suggest that the thermodynamic properties of siRNA play a critical role in determining the molecule's function and longevity, possibly biasing the steps involved in duplex unwinding and strand retention by RISC.     

Global profiling of miRNAs and the hairpin precursors: insights into miRNA processing and novel miRNA discovery

MicroRNAs (miRNAs) constitute an important class of small regulatory RNAs that are derived from distinct hairpin precursors (pre-miRNAs). In contrast to mature miRNAs, which have been characterized in numerous genome-wide studies of different organisms, research on global profiling of pre-miRNAs is limited. Here, using massive parallel sequencing, we have performed global characterization of both mouse mature and precursor miRNAs. In total, 87 369 704 and 252 003 sequencing reads derived from 887 mature and 281 precursor miRNAs were obtained, respectively. Our analysis revealed new aspects of miRNA/pre-miRNA processing and modification, including eight Ago2-cleaved pre-miRNAs, eight new instances of miRNA editing and exclusively 5′ tailed mirtrons. Furthermore, based on the sequences of both mature and precursor miRNAs, we developed a miRNA discovery pipeline, miRGrep, which does not rely on the availability of genome reference sequences. In addition to 239 known mouse pre-miRNAs, miRGrep predicted 41 novel ones with high confidence. Similar as known ones, the mature miRNAs derived from most of these novel loci showed both reduced abundance following Dicer knockdown and the binding with Argonaute2. Evaluation on data sets obtained from Caenorhabditis elegans and Caenorhabditis sp.11 demonstrated that miRGrep could be widely used for miRNA discovery in metazoans, especially in those without genome reference sequences.     

Instability of miRNA and cDNAs derivatives in RNA preparations

Micro RNAs (miRNAs) are small RNA molecules, which function as important regulators of gene expression. We found that RNA preparation methods commonly utilized for miRNA expression studies yield highly unstable miRNAs. We studied the stability of four miRNAs belonging to different miRNAs families. A significant degradation of these molecules may be observed already three days after RNA isolation. Moreover, the respective cDNAs are highly unstable as well. Our findings indicate that instability of miRNAs and their cDNAs should be considered when designing miRNA expression studies.     

Study of E. coli Hfq’s RNA annealing acceleration and duplex destabilization activities using substrates with different GC-contents

Folding of RNA molecules into their functional three-dimensional structures is often supported by RNA chaperones, some of which can catalyse the two elementary reactions helix disruption and helix formation. Hfq is one such RNA chaperone, but its strand displacement activity is controversial. Whereas some groups found Hfq to destabilize secondary structures, others did not observe such an activity with their RNA substrates. We studied Hfq’s activities using a set of short RNAs of different thermodynamic stabilities (GC-contents from 4.8% to 61.9%), but constant length. We show that Hfq’s strand displacement as well as its annealing activity are strongly dependent on the substrate’s GC-content. However, this is due to Hfq’s preferred binding of AU-rich sequences and not to the substrate’s thermodynamic stability. Importantly, Hfq catalyses both annealing and strand displacement with comparable rates for different substrates, hinting at RNA strand diffusion and annealing nucleation being rate-limiting for both reactions. Hfq’s strand displacement activity is a result of the thermodynamic destabilization of the RNA through preferred single-strand binding whereas annealing acceleration is independent from Hfq’s thermodynamic influence. Therefore, the two apparently disparate activities annealing acceleration and duplex destabilization are not in energetic conflict with each other.     

Characteristic comparison between two types of miRNA precursors in metazoan species

MicroRNAs (miRNAs) are a class of small non-coding RNAs discovered in recent years, which are found to play important regulatory roles in various organisms. As the number of experimentally validated miRNAs is rapidly increasing, systematic analysis on the characteristics of these known miRNAs is necessary and indispensable, especially for computational prediction of new miRNAs. We extensively analyzed precursor sequences for all experimentally validated mature miRNAs in metazoan species, focusing on the characteristics at the level of primary sequences and secondary structures. An observation over the secondary structures of 2729 miRNA precursors (pre-miRNAs) reveals that these hairpin structures can be approximately classified into two types: one with a hairpin loop, and the other with multiple loops. Interestingly, the two types of pre-miRNAs show significant differences in both sequence and structure characteristics, and our study indicates that separate consideration on each type of pre-miRNAs is more reasonable, especially in computational prediction. Besides, we develop a new criterion called mAMFE which shows robust discriminative power in distinguishing pre-miRNAs against other RNAs, thus can potentially serve as a discriminative feature in prediction of new pre-miRNAs.     

An integrated analysis of miRNA and mRNA expressions in non-small cell lung cancers

Using DNA microarrays, we generated both mRNA and miRNA expression data from 6 non-small cell lung cancer (NSCLC) tissues and their matching normal control from adjacent tissues to identify potential miRNA markers for diagnostics. We demonstrated that hsa-miR-96 is significantly and consistently up-regulated in all 6 NSCLCs. We validated this result in an independent set of 35 paired tumors and their adjacent normal tissues, as well as their sera that are collected before surgical resection or chemotherapy, and the results suggested that hsa-miR-96 may play an important role in NSCLC development and has great potential to be used as a noninvasive marker for diagnosing NSCLC. We predicted potential miRNA target mRNAs based on different methods (TargetScan and miRanda). Further classification of miRNA regulated genes based on their relationship with miRNAs revealed that hsa-miR-96 and certain other miRNAs tend to down-regulate their target mRNAs in NSCLC development, which have expression levels permissive to direct interaction between miRNAs and their target mRNAs. In addition, we identified a significant correlation of miRNA regulation with genes coincide with high density of CpG islands, which suggests that miRNA may represent a primary regulatory mechanism governing basic cellular functions and cell differentiations, and such mechanism may be complementary to DNA methylation in repressing or activating gene expression.     

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.     

Focus on RNA isolation: obtaining RNA for microRNA (miRNA) expression profiling analyses of neural tissue

MicroRNAs (miRNAs) are present in all known plant and animal tissues and appear to be somewhat concentrated in the mammalian nervous system. Many different miRNA expression profiling platforms have been described. However, relatively little research has been published to establish the importance of 'upstream' variables in RNA isolation for neural miRNA expression profiling. We tested whether apparent changes in miRNA expression profiles may be associated with tissue processing, RNA isolation techniques, or different cell types in the sample. RNA isolation was performed on a single brain sample using eight different RNA isolation methods, and results were correlated using a conventional miRNA microarray and then cross-referenced to Northern blots. Differing results were seen between samples obtained using different RNA isolation techniques and between microarray and Northern blot results. Another complication of miRNA microarrays is tissue-level heterogeneity of cellular composition. To investigate this phenomenon, miRNA expression profiles were determined and compared between highly-purified primary cerebral cortical cell preparations of rat primary E15-E18 neurons versus rat primary E15-E18 astrocytes. Finally, to assess the importance of dissecting human brain gray matter from subjacent white matter in cerebral cortical studies, miRNA expression profiles were compared between gray matter and immediately contiguous white matter. The results suggest that for microarray studies, cellular composition is important, and dissecting white matter from gray matter improves the specificity of the results. Based on these data, recommendations for miRNA expression profiling in neural tissues, and considerations worthy of further study, are discussed.     

Tissue-dependent paired expression of miRNAs

It is believed that depending on the thermodynamic stability of the 5′-strand and the 3′-strand in the stem-loop structure of a precursor microRNA (pre-miRNA), cells preferentially select the less stable one (called the miRNA or guide strand) and destroy the other one (called the miRNA* or passenger strand). However, our expression profiling analyses revealed that both strands could be co-accumulated as miRNA pairs in some tissues while being subjected to strand selection in other tissues. Our target prediction and validation assays demonstrated that both strands of a miRNA pair could target equal numbers of genes and that both were able to suppress the expression of their target genes. Our finding not only suggests that the numbers of miRNAs and their targets are much greater than what we previously thought, but also implies that novel mechanisms are involved in the tissue-dependent miRNA biogenesis and target selection process.     

miRTour: Plant miRNA and target prediction tool

MicroRNAs (miRNAs) are important negative regulators of gene expression in plant and animals, which are endogenously produced from their own genes. Computational comparative approach based on evolutionary conservation of mature miRNAs has revealed a number of orthologs of known miRNAs in different plant species. The homology-based plant miRNA discovery, followed by target prediction, comprises several steps, which have been done so far manually. Here, we present the bioinformatics pipeline miRTour which automates all the steps of miRNA similarity search, miRNA precursor selection, target prediction and annotation, each of them performed with the same set of input sequences. AVAILABILITY: The database is available for free at http://bio2server.bioinfo.uni-plovdiv.bg/miRTour/     

A wide repertoire of miRNA binding sites: prediction and functional implications

MOTIVATION: Over the past decade, deciphering the roles of microRNAs (miRNAs) has relied heavily upon the identification of their targets. Most of the targets that were computationally and experimentally characterized were evolutionarily conserved 'seed' targets, containing a perfect 6-8 nt match between the miRNA 5(')-region and the messenger RNA (mRNA). Gradually, it has become evident that other types of miRNA binding can confer target regulation, but their characterization has been lagging behind. RESULTS: Here, we complement the putative evolutionarily-conserved seed-containing targets by a wide repertoire of putative targets exhibiting a variety of miRNA binding patterns, predicted by our algorithm RepTar. These include non-conserved sites, 'seed' binding sites with G:U-wobbles within the seed, '3(') compensatory' sites and 'centered' sites. Apart from the centered sites, we demonstrate the functionality of these sites and characterize the target profile of a miRNA by the types of binding sites predicted in its target 3(') UTRs. We find that different miRNAs have individual target profiles, with some more inclined to seed binding and others more inclined to binding through 3(') compensatory sites. This diversity in targeting patterns is also evident within several miRNA families (defined by common seed sequences), leading to divergence in the target sets of members of the same family. The prediction of non-conventional miRNA targets is also beneficial in the search for targets of the non-conserved viral miRNAs. Analyzing the cellular targets of viral miRNAs, we show that viral miRNAs use various binding patterns to exploit cellular miRNA binding sites and suggest roles for these targets in virus-host interactions.