Prediction of Conserved Precursors of miRNAs and Their Mature Forms by Integrating Position-Specific Structural Features

MicroRNA (miRNA) precursor hairpins have a unique secondary structure, nucleotide length, and nucleotide content that are in most cases evolutionarily conserved. The aim of this study was to utilize position-specific features of miRNA hairpins to improve their identification. To this end, we defined the evolutionary and structurally conserved features in each position of miRNA hairpins with heuristically derived values, which were successfully integrated using a probabilistic framework. Our method, miRRim2, can not only accurately detect miRNA hairpins, but infer the location of a mature miRNA sequence. To evaluate the accuracy of miRRim2, we designed a cross validation test in which the whole human genome was used for evaluation. miRRim2 could more accurately detect miRNA hairpins than the other computational predictions that had been performed on the human genome, and detect the position of the 5'-end of mature miRNAs with sensitivity and positive predictive value (PPV) above 0.4. To further evaluate miRRim2 on independent data, we applied it to the Ciona intestinalis genome. Our method detected 47 known miRNA hairpins among top 115 candidates, and pinpointed the 5'-end of mature miRNAs with sensitivity and PPV about 0.4. When our results were compared with deep-sequencing reads of small RNA libraries from Ciona intestinalis cells, we found several candidates in which the predicted mature miRNAs were in good accordance with deep-sequencing results.     

Transposable element-associated microRNA hairpins produce 21-nt sRNAs integrated into typical microRNA pathways in rice

microRNAs (miRNAs) are a class of small RNAs (sRNAs) of ～21 nucleotides (nt) in length processed from foldback hairpins by DICER-LIKE1 (DCL1) or DCL4. They regulate the expression of target mRNAs by base pairing through RNA-induced silencing complex (RISC). In the RISC, ARGONAUTE1 (AGO1) is the key protein that cleaves miRNA targets at position ten of a miRNA:target duplex. The authenticity of many annotated rice miRNA hairpins is under debate because of their homology to repeat sequences. Some of them, like miR1884b, have been removed from the current release of miRBase based on incomplete information. In this study, we investigated the association of transposable element (TE)-derived miRNAs with typical miRNA pathways (DCL1/4- and AGO1-dependent) using publicly available deep sequencing datasets. Seven miRNA hairpins with 13 unique sRNAs were specifically enriched in AGO1 immunoprecipitation samples and relatively reduced in DCL1/4 knockdown genotypes. Interestingly, these species are ～21-nt long, instead of 24-nt as annotated in miRBase and the literature. Their expression profiles meet current criteria for functional annotation of miRNAs. In addition, diagnostic cleavage tags were found in degradome datasets for predicted target mRNAs. Most of these miRNA hairpins share significant homology with miniature inverted-repeat transposable elements, one type of abundant DNA transposons in rice. Finally, the root-specific production of a 24-nt miRNA-like sRNA was confirmed by RNA blot for a novel EST that maps to the 3′-UTR of a candidate pseudogene showing extensive sequence homology to miR1884b hairpin. Our data are consistent with the hypothesis that TEs can serve as a driving force for the evolution of some MIRNAs, where co-opting of DICER-LIKE1/4 processing and integration into AGO1 could exapt transcribed TE-associated hairpins into typical miRNA pathways.     

MiRduplexSVM: A High-Performing MiRNA-Duplex Prediction and Evaluation Methodology

We address the problem of predicting the position of a miRNA duplex on a microRNA hairpin via the development and application of a novel SVM-based methodology. Our method combines a unique problem representation and an unbiased optimization protocol to learn from mirBase19.0 an accurate predictive model, termed MiRduplexSVM. This is the first model that provides precise information about all four ends of the miRNA duplex. We show that (a) our method outperforms four state-of-the-art tools, namely MaturePred, MiRPara, MatureBayes, MiRdup as well as a Simple Geometric Locator when applied on the same training datasets employed for each tool and evaluated on a common blind test set. (b) In all comparisons, MiRduplexSVM shows superior performance, achieving up to a 60% increase in prediction accuracy for mammalian hairpins and can generalize very well on plant hairpins, without any special optimization. (c) The tool has a number of important applications such as the ability to accurately predict the miRNA or the miRNA*, given the opposite strand of a duplex. Its performance on this task is superior to the 2nts overhang rule commonly used in computational studies and similar to that of a comparative genomic approach, without the need for prior knowledge or the complexity of performing multiple alignments. Finally, it is able to evaluate novel, potential miRNAs found either computationally or experimentally. In relation with recent confidence evaluation methods used in miRBase, MiRduplexSVM was successful in identifying high confidence potential miRNAs.     

Phylogenetic shadowing and computational identification of human microRNA genes

We sequenced 122 miRNAs in 10 primate species to reveal conservation characteristics of miRNA genes. Strong conservation is observed in stems of miRNA hairpins and increased variation in loop sequences. Interestingly, a striking drop in conservation was found for sequences immediately flanking the miRNA hairpins. This characteristic profile was employed to predict novel miRNAs using cross-species comparisons. Nine hundred and seventy-six candidate miRNAs were identified by scanning whole-genome human/mouse and human/rat alignments. Most of the novel candidates are conserved also in other vertebrates (dog, cow, chicken, opossum, zebrafish). Northern blot analysis confirmed the expression of mature miRNAs for 16 out of 69 representative candidates. Additional support for the expression of 179 novel candidates can be found in public databases, their presence in gene clusters, and literature that appeared after these predictions were made. Taken together, these results suggest the presence of significantly higher numbers of miRNAs in the human genome than previously estimated.     

MicroRNA Prediction Using a Fixed-Order Markov Model Based on the Secondary Structure Pattern

Predicting miRNAs is an arduous task, due to the diversity of the precursors and complexity of enzyme processes. Although several prediction approaches have reached impressive performances, few of them could achieve a full-function recognition of mature miRNA directly from the candidate hairpins across species. Therefore, researchers continue to seek a more powerful model close to biological recognition to miRNA structure. In this report, we describe a novel miRNA prediction algorithm, known as FOMmiR, using a fixed-order Markov model based on the secondary structural pattern. For a training dataset containing 809 human pre-miRNAs and 6441 human pseudo-miRNA hairpins, the model’s parameters were defined and evaluated. The results showed that FOMmiR reached 91% accuracy on the human dataset through 5-fold cross-validation. Moreover, for the independent test datasets, the FOMmiR presented an outstanding prediction in human and other species including vertebrates, Drosophila, worms and viruses, even plants, in contrast to the well-known algorithms and models. Especially, the FOMmiR was not only able to distinguish the miRNA precursors from the hairpins, but also locate the position and strand of the mature miRNA. Therefore, this study provides a new generation of miRNA prediction algorithm, which successfully realizes a full-function recognition of the mature miRNAs directly from the hairpin sequences. And it presents a new understanding of the biological recognition based on the strongest signal’s location detected by FOMmiR, which might be closely associated with the enzyme cleavage mechanism during the miRNA maturation.     

Integrated sequence-structure motifs suffice to identify microRNA precursors

Background

Upwards of 1200 miRNA loci have hitherto been annotated in the human genome. The specific features defining a miRNA precursor and deciding its recognition and subsequent processing are not yet exhaustively described and miRNA loci can thus not be computationally identified with sufficient confidence.

Results

We rendered pre-miRNA and non-pre-miRNA hairpins as strings of integrated sequence-structure information, and used the software Teiresias to identify sequence-structure motifs (ss-motifs) of variable length in these data sets. Using only ss-motifs as features in a Support Vector Machine (SVM) algorithm for pre-miRNA identification achieved 99.2% specificity and 97.6% sensitivity on a human test data set, which is comparable to previously published algorithms employing combinations of sequence-structure and additional features. Further analysis of the ss-motif information contents revealed strongly significant deviations from those of the respective training sets, revealing important potential clues as to how the sequence and structural information of RNA hairpins are utilized by the miRNA processing apparatus.

Conclusion

Integrated sequence-structure motifs of variable length apparently capture nearly all information required to distinguish miRNA precursors from other stem-loop structures.     

Integration of expressed sequence tag data flanking predicted RNA secondary structures facilitates novel non-coding RNA discovery

Many computational methods have been used to predict novel non-coding RNAs (ncRNAs), but none, to our knowledge, have explicitly investigated the impact of integrating existing cDNA-based Expressed Sequence Tag (EST) data that flank structural RNA predictions. To determine whether flanking EST data can assist in microRNA (miRNA) prediction, we identified genomic sites encoding putative miRNAs by combining functional RNA predictions with flanking ESTs data in a model consistent with miRNAs undergoing cleavage during maturation. In both human and mouse genomes, we observed that the inclusion of flanking ESTs adjacent to and not overlapping predicted miRNAs significantly improved the performance of various methods of miRNA prediction, including direct high-throughput sequencing of small RNA libraries. We analyzed the expression of hundreds of miRNAs predicted to be expressed during myogenic differentiation using a customized microarray and identified several known and predicted myogenic miRNA hairpins. Our results indicate that integrating ESTs flanking structural RNA predictions improves the quality of cleaved miRNA predictions and suggest that this strategy can be used to predict other non-coding RNAs undergoing cleavage during maturation.     

Transposon-based RNAi delivery system for generating knockdown cell lines

RNA interference is rapidly becoming a powerful tool for genetic analyses in mammalian systems. A potential drawback to transient small inhibitory RNA silencing is the short duration of downregulation it confers, usually only 24-72h. Viral-based vector systems for the long-term delivery of RNA hairpins have been developed, yet they require expertise in viral production and transduction. Here we describe a simple plasmid-based system for the generation of long-term gene knockdown utilizing RNA interference combined with the gene delivery capabilities of the mammalian Tc1-like transposon Sleeping Beauty. Designated Maleficent, this system is shown to downregulate exogenous expression of GFP in a constitutively positive cell line. In addition, targeting of the endogenously expressed lamin A gene results in long-term silencing with significant reduction in protein levels (> 95%). Maleficent therefore provides a relatively easy, efficient, and stable means of delivering RNAi hairpins to generate long-term gene-specific knockdown cell lines.     

Comprehensive analyses of microRNA gene evolution in paleopolyploid soybean genome

miRNA genes are thought to undergo quick birth and death processes in genomes and the emergence of a MIRNA‐like hairpin provides the base for functional miRNA gene formation. However, the factors affecting the formation of an active miRNA gene from a MIRNA‐like hairpin within a genome remain unclear. We performed a genome‐wide investigation of MIRNA‐like hairpin accumulation, expression, structural changes and relationships with annotated genomic features in the paleopolyploid soybean genome. Our results showed that adjacent gene and transposable element content, rates of genetic recombination at location of emergence, along with its own gene structure divergence greatly affected miRNA gene evolution. Further investigation suggested that miRNA genes from different duplication sources followed distinct evolutionary trajectories and that the accumulation of MIRNA‐like hairpins might be an important factor in causing long terminal repeat retrotransposons to lose activity during genome evolution.     

An efficient method to enhance gene silencing by using precursor microRNA designed small hairpin RNAs

Gene silencing can be mediated by small interfering RNA (siRNA) and microRNA (miRNA). To investigate the potential application of using a precursor microRNA (pre-miRNA) backbone for gene silencing, we studied the inhibition efficiency of exogenous GFP and endogenous GAPDH by conventional shRNA- and pre-miRNA-designed hairpins, respectively. In this study, the conventional shRNA-, pre-miRNA-30-, and pre-miRNA-155-designed hairpins targeting either GFP or GAPDH were transfected into the HEK293 cells that were mediated by the pSilencer-4.1-neo vector, which carries a modified RNA polymerase II-type CMV promoter. Comparisons with conventional GFP shRNA showed that GFP levels were reduced markedly by pre-miRNA-30- and pre-miRNA-155-designed GFP shRNAs by fluorescence microscopy. The consistent results from semi-quantitative RT-PCR and Western blot analysis revealed that pre-miRNA-30- and pre-miRNA-155-designed GFP shRNAs could suppress GFP expression significantly. As for endogenous GAPDH, the results from semi-quantitative RT-PCR and Western blot analysis showed that pre-miRNA-30- and pre-miRNA-155-designed GAPDH shRNAs could suppress GAPDH expression even more efficiently than conventional GAPDH shRNA. Together, this study confirmed the efficiency of gene silencing mediated by pre-miRNA-30- and pre-miRNA-155-designed shRNAs, demonstrating that pre-miRNA-designed hairpins are a good strategy for gene silencing.     

Structural basis of microRNA length variety

The biogenesis of human microRNAs (miRNAs) includes two RNA cleavage steps in which the activities of the RNases Drosha and Dicer are involved. miRNAs of diverse lengths are generated from different genes, and miRNAs that are heterogeneous in length are produced from a single miRNA gene. We determined the solution structures of many miRNA precursors and analysed the structural basis of miRNA length diversity using a new measure: the weighted average length of diced RNA (WALDI). We found that asymmetrical structural motifs present in precursor hairpins are primarily responsible for the length diversity of miRNAs generated by Dicer. High-resolution northern blots of miRNAs and their precursors revealed that both Dicer and Drosha cleavages of imperfect specificity contributed to the miRNA length heterogeneity. The relevance of these findings to the dynamics of the dicing complex, mRNA regulation by miRNA, RNA interference and miRNA technologies are discussed.