Inference of miRNA targets using evolutionary conservation and pathway analysis

Background

MicroRNAs have emerged as important regulatory genes in a variety of cellular processes and, in recent years, hundreds of such genes have been discovered in animals. In contrast, functional annotations are available only for a very small fraction of these miRNAs, and even in these cases only partially.

Results

We developed a general Bayesian method for the inference of miRNA target sites, in which, for each miRNA, we explicitly model the evolution of orthologous target sites in a set of related species. Using this method we predict target sites for all known miRNAs in flies, worms, fish, and mammals. By comparing our predictions in fly with a reference set of experimentally tested miRNA-mRNA interactions we show that our general method performs at least as well as the most accurate methods available to date, including ones specifically tailored for target prediction in fly. An important novel feature of our model is that it explicitly infers the phylogenetic distribution of functional target sites, independently for each miRNA. This allows us to infer species-specific and clade-specific miRNA targeting. We also show that, in long human 3' UTRs, miRNA target sites occur preferentially near the start and near the end of the 3' UTR. To characterize miRNA function beyond the predicted lists of targets we further present a method to infer significant associations between the sets of targets predicted for individual miRNAs and specific biochemical pathways, in particular those of the KEGG pathway database. We show that this approach retrieves several known functional miRNA-mRNA associations, and predicts novel functions for known miRNAs in cell growth and in development.

Conclusion

We have presented a Bayesian target prediction algorithm without any tunable parameters, that can be applied to sequences from any clade of species. The algorithm automatically infers the phylogenetic distribution of functional sites for each miRNA, and assigns a posterior probability to each putative target site. The results presented here indicate that our general method achieves very good performance in predicting miRNA target sites, providing at the same time insights into the evolution of target sites for individual miRNAs. Moreover, by combining our predictions with pathway analysis, we propose functions of specific miRNAs in nervous system development, inter-cellular communication and cell growth. The complete target site predictions as well as the miRNA/pathway associations are accessible on the ElMMo web server.     

Novel and mechanical stress-responsive MicroRNAs in Populus trichocarpa that are absent from Arabidopsis

MicroRNAs (miRNAs) are small, noncoding RNAs that can play crucial regulatory roles in eukaryotes by targeting mRNAs for silencing. To test whether miRNAs play roles in the regulation of wood development in tree species, we isolated small RNAs from the developing xylem of Populus trichocarpa stems and cloned 22 miRNAs. They are the founding members of 21 miRNA gene families for 48 miRNA sequences, represented by 98 loci in the Populus genome. A majority of these miRNAs were predicted to target developmental- and stress/defense-related genes and possible functions associated with the biosynthesis of cell wall metabolites. Of the 21 P. trichocarpa miRNA families, 11 have sequence conservation in Arabidopsis thaliana but exhibited species-specific developmental expression patterns, suggesting that even conserved miRNAs may have different regulatory roles in different species. Most unexpectedly, the remaining 10 miRNAs, for which 17 predicted targets were experimentally validated in vivo, are absent from the Arabidopsis genome, suggesting possible roles in tree-specific processes. In fact, the expression of a majority of the cloned miRNAs was upregulated or downregulated in woody stems in a manner consistent with tree-specific corrective growth against tension and compression stresses, two constant mechanical loads in trees. Our results show that plant miRNAs can be induced by mechanical stress and may function in one of the most critical defense systems for structural and mechanical fitness.     

Conservation and evolution of miRNA regulatory programs in plant development

Over the past two years, microarray technologies, large-scale small RNA and whole genome sequencing projects, and data mining have provided a wealth of information about the spectrum of miRNAs and miRNA targets present in different plant species and the alga Chlamydomonas. Such studies have shown that a number of key miRNA regulatory modules for plant development are conserved throughout the plant kingdom, suggesting that these programs were crucial to the colonization of land. New genetic and biochemical studies of miRNA pathways in Arabidopsis, the spatiotemporal expression patterns of several conserved miRNAs and their targets, and the characterization of mutations in Arabidopsis and maize have begun to reveal the functions of these ancient miRNA-regulated developmental programs. In addition to these conserved miRNAs, there are many clade and species-specific miRNAs, which have evolved more recently and whose functions are currently unknown.     

myMIR: a genome-wide microRNA targets identification and annotation tool

miRNA target genes prediction represents a crucial step in miRNAs functional characterization. In this context, the challenging issue remains predictions accuracy and recognition of false positive results. In this article myMIR, a web based system for increasing reliability of miRNAs predicted targets lists, is presented. myMIR implements an integrated pipeline for computing ranked miRNA::target lists and provides annotations for narrowing them down. The system relies on knowledge base data, suitably integrated in order to extend the functional characterization of targeted genes to miRNAs, by highlighting the search on over-represented annotation terms. Validation results show a dramatic reduction in the quantity of predictions and an increase in the sensitivity, when compared to other methods. This improves the predictions accuracy and allows the formulation of novel hypotheses on miRNAs functional involvement.     

A dictionary on microRNAs and their putative target pathways

While in the last decade mRNA expression profiling was among the most popular research areas, over the past years the study of non-coding RNAs, especially microRNAs (miRNAs), has gained increasing interest. For almost 900 known human miRNAs hundreds of pretended targets are known. However, there is only limited knowledge about putative systemic effects of changes in the expression of miRNAs and their regulatory influence. We determined for each known miRNA the biochemical pathways in the KEGG and TRANSPATH database and the Gene Ontology categories that are enriched with respect to its target genes. We refer to these pathways and categories as target pathways of the corresponding miRNA. Investigating target pathways of miRNAs we found a strong relation to disease-related regulatory pathways, including mitogen-activated protein kinase (MAPK) signaling cascade, Transforming growth factor (TGF)-beta signaling pathway or the p53 network. Performing a sophisticated analysis of differentially expressed genes of 13 cancer data sets extracted from gene expression omnibus (GEO) showed that targets of specific miRNAs were significantly deregulated in these sets. The respective miRNA target analysis is also a novel part of our gene set analysis pipeline GeneTrail. Our study represents a comprehensive theoretical analysis of the relationship between miRNAs and their predicted target pathways. Our target pathways analysis provides a ‘miRNA-target pathway’ dictionary, which enables researchers to identify target pathways of differentially regulated miRNAs.     

Isolation and identification of miRNAs in Jatropha curcas

MicroRNAs (miRNAs) are small noncoding RNAs that play crucial regulatory roles by targeting mRNAs for silencing. To identify miRNAs in Jatropha curcas L, a bioenergy crop, cDNA clones from two small RNA libraries of leaves and seeds were sequenced and analyzed using bioinformatic tools. Fifty-two putative miRNAs were found from the two libraries, among them six were identical to known miRNAs and 46 were novel. Differential expression patterns of 15 miRNAs in root, stem, leave, fruit and seed were detected using quantitative real-time PCR. Ten miRNAs were highly expressed in fruit or seed, implying that they may be involved in seed development or fatty acids synthesis in seed. Moreover, 28 targets of the isolated miRNAs were predicted from a jatropha cDNA library database. The miRNA target genes were predicted to encode a broad range of proteins. Sixteen targets had clear BLASTX hits to the Uniprot database and were associated with genes belonging to the three major gene ontology categories of biological process, cellular component, and molecular function. Four targets were identified for JcumiR004. By silencing JcumiR004 primary miRNA, expressions of the four target genes were up-regulated and oil composition were modulated significantly, indicating diverse functions of JcumiR004.     

Interactive Web-based Annotation of Plant MicroRNAs with iwa-miRNA

MicroRNAs (miRNAs) are important regulators of gene expression. The large-scale detection and profiling of miRNAs have been accelerated with the development of high-throughput small RNA sequencing (sRNA-Seq) techniques and bioinformatics tools. However, generating high-quality comprehensive miRNA annotations remains challenging due to the intrinsic complexity of sRNA-Seq data and inherent limitations of existing miRNA prediction tools. Here, we present iwa-miRNA, a Galaxy-based framework that can facilitate miRNA annotation in plant species by combining computational analysis and manual curation. iwa-miRNA is specifically designed to generate a comprehensive list of miRNA candidates, bridging the gap between already annotated miRNAs provided by public miRNA databases and new predictions from sRNA-Seq datasets. It can also assist users in selecting promising miRNA candidates in an interactive mode, contributing to the accessibility and reproducibility of genome-wide miRNA annotation. iwa-miRNA is user-friendly and can be easily deployed as a web application for researchers without programming experience. With flexible, interactive, and easy-to-use features, iwa-miRNA is a valuable tool for the annotation of miRNAs in plant species with reference genomes. We also illustrate the application of iwa-miRNA for miRNA annotation using data from plant species with varying genomic complexity. The source codes and web server of iwa-miRNA are freely accessible at http://iwa-miRNA.omicstudio.cloud/.     

miRDB: a microRNA target prediction and functional annotation database with a wiki interface

MicroRNAs (miRNAs) are short noncoding RNAs that are involved in the regulation of thousands of gene targets. Recent studies indicate that miRNAs are likely to be master regulators of many important biological processes. Due to their functional importance, miRNAs are under intense study at present, and many studies have been published in recent years on miRNA functional characterization. The rapid accumulation of miRNA knowledge makes it challenging to properly organize and present miRNA function data. Although several miRNA functional databases have been developed recently, this remains a major bioinformatics challenge to miRNA research community. Here, we describe a new online database system, miRDB, on miRNA target prediction and functional annotation. Flexible web search interface was developed for the retrieval of target prediction results, which were generated with a new bioinformatics algorithm we developed recently. Unlike most other miRNA databases, miRNA functional annotations in miRDB are presented with a primary focus on mature miRNAs, which are the functional carriers of miRNA-mediated gene expression regulation. In addition, a wiki editing interface was established to allow anyone with Internet access to make contributions on miRNA functional annotation. This is a new attempt to develop an interactive community-annotated miRNA functional catalog. All data stored in miRDB are freely accessible at http://mirdb.org.     

miRSystem: An Integrated System for Characterizing Enriched Functions and Pathways of MicroRNA Targets

Background

Many prediction tools for microRNA (miRNA) targets have been developed, but inconsistent predictions were observed across multiple algorithms, which can make further analysis difficult. Moreover, the nomenclature of human miRNAs changes rapidly. To address these issues, we developed a web-based system, miRSystem, for converting queried miRNAs to the latest annotation and predicting the function of miRNA by integrating miRNA target gene prediction and function/pathway analyses.

Results

First, queried miRNA IDs were converted to the latest annotated version to prevent potential conflicts resulting from multiple aliases. Next, by combining seven algorithms and two validated databases, potential gene targets of miRNAs and their functions were predicted based on the consistency across independent algorithms and observed/expected ratios. Lastly, five pathway databases were included to characterize the enriched pathways of target genes through bootstrap approaches. Based on the enriched pathways of target genes, the functions of queried miRNAs could be predicted.

Conclusions

MiRSystem is a user-friendly tool for predicting the target genes and their associated pathways for many miRNAs simultaneously. The web server and the documentation are freely available at http://mirsystem.cgm.ntu.edu.tw/.     

Identification of miRNAs and Their Target Genes in Peach (Prunus persica L.) Using High-Throughput Sequencing and Degradome Analysis

MicroRNAs play critical roles in various biological and metabolic processes. The function of miRNAs has been widely studied in model plants such as Arabidopsis and rice. However, the number of identified miRNAs and related miRNA targets in peach (Prunus persica) is limited. To understand further the relationship between miRNAs and their target genes during tissue development in peach, a small RNA library and three degradome libraries were constructed from three tissues for deep sequencing. We identified 117 conserved miRNAs and 186 novel miRNA candidates in peach by deep sequencing and 19 conserved miRNAs and 13 novel miRNAs were further evaluated for their expression by RT-qPCR. The number of gene targets that were identified for 26 conserved miRNA families and 38 novel miRNA candidates, were 172 and 87, respectively. Some of the identified miRNA targets were abundantly represented as conserved miRNA targets in plant. However, some of them were first identified and showed important roles in peach development. Our study provides information concerning the regulatory network of miRNAs in peach and advances our understanding of miRNA functions during tissue development.     

Sequence variation of MicroRNAs and their binding sites in Arabidopsis

Major differences exist between plants and animals both in the extent of microRNA (miRNA)-based gene regulation and the sequence complementarity requirements for miRNA-messenger RNA pairing. Whether these differences affect how these sites evolve at the molecular level is unknown. To determine the extent of sequence variation at miRNAs and their targets in a plant species, we resequenced 16 miRNA families (66 miRNAs in total) and all 52 of the characterized binding sites for these miRNAs in the plant model Arabidopsis (Arabidopsis thaliana), accounting for around 50% of the known miRNAs and binding sites in this species. As has been shown previously in humans, we find that both miRNAs and their target binding sites have very low nucleotide variation and divergence compared to their flanking sequences in Arabidopsis, indicating strong purifying selection on these sites in this species. Sequence data flanking the mature miRNAs, however, exhibit normal levels of polymorphism for the accessions in this study and, in some cases, nonneutral evolution or subtle effects on predicted pre-miRNA secondary structure, suggesting that there is raw material for the differential function of miRNA alleles. Overall, our results show that despite differences in the architecture of miRNA-based regulation, miRNAs and their targets are similarly constrained in both plants and animals.