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

MicroRNA-Target Identification: A Combinatorial In Silico Approach

Contemporary computational target prediction tools with their distinctive properties and stringency have been playing a vital role in pursuing putative targets for a solitary miRNA or a subcategory of miRNAs. These tools utilize a defined set of probabilistic algorithms, machine learning techniques, and information of experimentally validated miRNA targets to provide the best selection. However, there are numerous false-positive predictions, and a method to choose an archetypal approach and put the data provided by the prediction tools into context is still lacking. Moreover, sensitivity, specificity, and overall efficiency of a single tool have not yet been achieved. Therefore, a systematic combination of selective online tools combining elementary and advanced factors of miRNA target identification might reinforce the current target prediction regime. The focus of this study was to build a comprehensive workflow by combining six available online tools to facilitate the current understanding of miRNA-target prediction.

A combinatorial in silico approach for microRNA-target identification: Order out of chaos

Contemporary computational microRNA(miRNA)-target prediction tools have been playing a vital role in pursuing putative targets for a solitary miRNA or a group of miRNAs. These tools utilise a set of probabilistic algorithms, machine learning techniques and analyse experimentally validated miRNA targets to identify the potential miRNA-target pairs. Unfortunately, current tools generate a huge number of false-positive predictions. A standardized approach with a single tool or a combination of tools is still lacking. Moreover, sensitivity, specificity and overall efficiency of any single tool are yet to be satisfactory. Therefore, a systematic combination of selective online tools combining the factors regarding miRNA-target identification would be valuable as an miRNA-target prediction scheme. The focus of this study was to develop a theoretical framework by combining six available online tools to facilitate the current understanding of miRNA-target identification.

Inconsistencies and Limitations of Current MicroRNA Target Identification Methods

MicroRNAs and their Argonaute protein partners constitute the RISC complex, which can repress specific target mRNAs. The identification of microRNA targets is of central importance, and various experimental and computational methods have been developed over the last 15 years. Most experimental methods are based on the assumption that mRNAs which interact physically with the RISC complex constitute regulatory targets and, similarly, some computational methods only aim at predicting physical interactors for RISC. Besides specific limitations, which we discuss for each method, the mere concept of assuming a functional role for every detected molecular event is likely to identify many deceptive interactions (i.e., interactions that really exist at the molecular scale, but without controlling any biological function at the macroscopic scale).In order to select biologically important interactions, some computational tools interrogate the phylogenetic conservation of microRNA/mRNA interactions, thus theoretically selecting only biologically relevant targets. Yet even comparative genomics can yield false positives.Conceptual and technical limitations for all these techniques tend to be overlooked by the scientific community. This review sums them up, emphasizing on the implications of these issues on our understanding of microRNA biology.

Genetic variant rs3750625 in the 3'UTR of ADRA2A affects stress-dependent acute pain severity after trauma and alters a microRNA-34a regulatory site

α2A adrenergic receptor (α2A-AR) activation has been shown in animal models to play an important role in regulating the balance of acute pain inhibition vs facilitation after both physical and psychological stress. To our knowledge, the influence of genetic variants in the gene encoding α2A-AR, ADRA2A, on acute pain outcomes in humans experiencing traumatic stress has not been assessed. In this study, we tested whether a genetic variant in the 3'UTR of ADRA2A, rs3750625, is associated with acute musculoskeletal pain (MSP) severity following motor vehicle collision (MVC, n = 948) and sexual assault (n = 84), and whether this influence was affected by stress severity. We evaluated rs3750625 because it is located in the seed binding region of miR-34a, a microRNA (miRNA) known to regulate pain and stress responses. In both cohorts, the minor allele at rs3750625 was associated with increased musculoskeletal pain in distressed individuals (stress*rs3750625 P = 0.043 for MVC cohort and P = 0.007 for sexual assault cohort). We further found that (1) miR-34a binds the 3'UTR of ADRA2A, (2) the amount of repression is greater when the minor (risk) allele is present, (3) miR-34a in the IMR-32 adrenergic neuroblastoma cell line affects ADRA2A expression, (4) miR-34a and ADRA2A are expressed in tissues known to play a role in pain and stress, (5) following forced swim stress exposure, rat peripheral nerve tissue expression changes are consistent with miR-34a regulation of ADRA2A. Together, these results suggest that ADRA2A rs3750625 contributes to poststress musculoskeletal pain severity by modulating miR-34a regulation.

Oxidative stress, microRNAs and cytosolic calcium homeostasis

Reactive oxygen species increase cytosolic [Ca(2+)], (Cai), and also modulate the expression of some microRNAs (miRNAs), however the link among oxidative stress, miRNAs and Cai is poorly characterized. In this review we have focused on three groups of miRNAs: (a) miRNAs that are modulated both by ROS and Cai: miR-181a and miR-205; (b) miRNAs that are modulated by ROS and have an effect on Cai: miR-1, miR-21, miR-24, miR-25, miR-185 and miR-214; (c) miRNAs that modulate both ROS and Cai: miR-133; miR-145, miR-495, and we have analyzed their effects on cell signaling and cell function. Finally, in the last section we have examined the role of these miRNAs in the skin, under conditions associated with enhanced oxidative stress, i.e. skin aging, the response to ultraviolet light and two important skin diseases, psoriasis and atopic dermatitis. It is apparent that although some experimental evidence is already available on (a) the role of Cai in miRNAs expression and (b) on the ability of some miRNAs to modulate Cai-dependent intracellular signaling, these research lines are still largely unexplored and represent important areas of future studies.

Inferred miRNA activity identifies miRNA-mediated regulatory networks underlying multiple cancers

MOTIVATION

MicroRNAs (miRNAs) play a key role in regulating tumor progression and metastasis. Identifying key miRNAs, defined by their functional activities, can provide a deeper understanding of biology of miRNAs in cancer. However, miRNA expression level cannot accurately reflect miRNA activity.

RESULTS

We developed a computational approach, ActMiR, for identifying active miRNAs and miRNA-mediated regulatory mechanisms. Applying ActMiR to four cancer datasets in The Cancer Genome Atlas (TCGA), we showed that (i) miRNA activity was tumor subtype specific; (ii) genes correlated with inferred miRNA activities were more likely to enrich for miRNA binding motifs; (iii) expression levels of these genes and inferred miRNA activities were more likely to be negatively correlated. For the four cancer types in TCGA we identified 77-229 key miRNAs for each cancer subtype and annotated their biological functions. The miRNA-target pairs, predicted by our ActMiR algorithm but not by correlation of miRNA expression levels, were experimentally validated. The functional activities of key miRNAs were further demonstrated to be associated with clinical outcomes for other cancer types using independent datasets. For ER(-)/HER2(-) breast cancers, we identified activities of key miRNAs let-7d and miR-18a as potential prognostic markers and validated them in two independent ER(-)/HER2(-) breast cancer datasets. Our work provides a novel scheme to facilitate our understanding of miRNA. In summary, inferred activity of key miRNA provided a functional link to its mediated regulatory network, and can be used to robustly predict patient's survival.

AVAILABILITY AND IMPLEMENTATION

the software is freely available at http://research.mssm.edu/integrative-network-biology/Software.html.

CONTACT

jun.zhu@mssm.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

miRNA-target chimeras reveal miRNA 3'-end pairing as a major determinant of Argonaute target specificity

microRNAs (miRNAs) act as sequence-specific guides for Argonaute (AGO) proteins, which mediate posttranscriptional silencing of target messenger RNAs. Despite their importance in many biological processes, rules governing AGO–miRNA targeting are only partially understood. Here we report a modified AGO HITS-CLIP strategy termed CLEAR (covalent ligation of endogenous Argonaute-bound RNAs)-CLIP, which enriches miRNAs ligated to their endogenous mRNA targets. CLEAR-CLIP mapped ∼130,000 endogenous miRNA–target interactions in mouse brain and ∼40,000 in human hepatoma cells. Motif and structural analysis define expanded pairing rules for over 200 mammalian miRNAs. Most interactions combine seed-based pairing with distinct, miRNA-specific patterns of auxiliary pairing. At some regulatory sites, this specificity confers distinct silencing functions to miRNA family members with shared seed sequences but divergent 3′-ends. This work provides a means for explicit biochemical identification of miRNA sites in vivo, leading to the discovery that miRNA 3′-end pairing is a general determinant of AGO binding specificity.

microRNAs (miRNAs) act as sequence-specific guides for Argonaute (AGO) proteins. By using a modified AGO HITS-CLIP strategy that enriches miRNAs ligated to their endogenous mRNA targets, here the authors show that miRNA 3' end pairing is a general determinant of AGO binding specificity.

Cupid: simultaneous reconstruction of microRNA-target and ceRNA networks

We introduce a method for simultaneous prediction of microRNA–target interactions and their mediated competitive endogenous RNA (ceRNA) interactions. Using high-throughput validation assays in breast cancer cell lines, we show that our integrative approach significantly improves on microRNA–target prediction accuracy as assessed by both mRNA and protein level measurements. Our biochemical assays support nearly 500 microRNA–target interactions with evidence for regulation in breast cancer tumors. Moreover, these assays constitute the most extensive validation platform for computationally inferred networks of microRNA–target interactions in breast cancer tumors, providing a useful benchmark to ascertain future improvements.

Binding sites of miR-1273 family on the mRNA of target genes

This study examined binding sites of 2,578 miRNAs in the mRNAs of 12,175 human genes using the MirTarget program. It found that the miRNAs of miR-1273 family have between 33 and 1,074 mRNA target genes, with a free hybridization energy of 90% or more of its maximum value. The miR-1273 family consists of miR-1273a, miR-1273c, miR-1273d, miR-1273e, miR-1273f, miR-1273g-3p, miR-1273g-5p, miR-1273h-3p, and miR-1273h-5p. Unique miRNAs (miR-1273e, miR-1273f, and miR-1273g-3p) have more than 400 target genes. We established 99 mRNA nucleotide sequences that contain arranged binding sites for the miR-1273 family. High conservation of each miRNA binding site in the mRNA of the target genes was found. The arranged binding sites of the miR-1273 family are located in the 5′UTR, CDS, or 3′UTR of many mRNAs. Five repeating sites containing some of the miR-1273 family's binding sites were found in the 3′UTR of several target genes. The oligonucleotide sequences of miR-1273 binding sites located in CDSs code for homologous amino acid sequences in the proteins of target genes. The biological role of unique miRNAs was also discussed.

MicroRNA-125 family members exert a similar role in the regulation of murine hematopoiesis

MicroRNAs (miRNAs) are crucial for proper functioning of hematopoietic stem and progenitor cells (HSPCs). Members of the miRNA-125 family (consisting of miR-125a, miR-125b1, and miR-125b2) are known to confer a proliferative advantage on cells upon overexpression, to decrease the rate of apoptosis by targeting proapoptotic genes, and to promote differentiation toward the myeloid lineage in mice. However, many distinct biological effects of the three miR-125 species have been reported as well. In the current study, we set out to assess whether the three miRNA-125s that carry identical seed sequences could be functionally different. Our data show that overexpression of each of the three miR-125 family members preserves HSPCs in a primitive state in vitro, results in a competitive advantage upon serial transplantation, and promotes skewing toward the myeloid lineage. All miR-125 family members decreased the pool of phenotypically defined Lin(-)Sca(+)Kit(+)CD48(-)CD150(+) long-term hematopoietic stem cells, simultaneously increasing the self-renewal activity upon secondary transplantation. The downregulation of miR-125s in hematopoietic stem cells abolishes these effects and impairs long-term contribution to blood cell production. The introduction of a point mutation within the miRNA-125 seed sequence abolishes all abovementioned effects and leads to the restoration of normal hematopoiesis. Our results show that all miR-125 family members are similar in function, they likely operate in a seed-sequence-dependent manner, and they induce a highly comparable hematopoietic phenotype.

Similarity in targets with REST points to neural and glioblastoma related miRNAs

There are groups of genes that need coordinated repression in multiple contexts, for example if they code for proteins that work together in a pathway or in a protein complex. Redundancy of biological regulatory networks implies that such coordinated repression might occur at both the pre- and post-transcriptional level, though not necessarily simultaneously or under the same conditions. Here, we propose that such redundancy in the global regulatory network can be detected by the overlap between the putative targets of a transcriptional repressor, as identified by a ChIP-seq experiment, and predicted targets of a microRNA (miRNA). To test this hypothesis, we used publicly available ChIP-seq data of the neural transcriptional repressor RE1 silencing transcription factor (REST) from 15 different cell samples. We found 20 miRNAs, each of which shares a significant amount of predicted targets with REST. The set of predicted associations between these 20 miRNAs and the overlapping REST targets is enriched in known miRNA targets. Many of the detected miRNAs have functions related to neural identity and glioblastoma, which could be expected from their overlap in targets with REST. We propose that the integration of experimentally determined transcription factor binding sites with miRNA-target predictions provides functional information on miRNAs.

Developing microRNA screening as a functional genomics tool for disease research

Originally discovered as regulators of developmental timing in C. elegans, microRNAs (miRNAs) have emerged as modulators of nearly every cellular process, from normal development to pathogenesis. With the advent of whole genome libraries of miRNA mimics suitable for high throughput screening, it is possible to comprehensively evaluate the function of each member of the miRNAome in cell-based assays. Since the relatively few microRNAs in the genome are thought to directly regulate a large portion of the proteome, miRNAome screening, coupled with the identification of the regulated proteins, might be a powerful new approach to gaining insight into complex biological processes.

Oxidative stress and microRNAs in vascular diseases

Oxidative stress has been demonstrated to play a causal role in different vascular diseases, such as hypertension, diabetic vasculopathy, hypercholesterolemia and atherosclerosis. Indeed, increased reactive oxygen species (ROS) production is known to impair endothelial and vascular smooth muscle cell functions, contributing to the development of cardiovascular diseases. MicroRNAs (miRNAs) are non-coding RNA molecules that modulate the stability and/or the translational efficiency of target messenger RNAs. They have been shown to be modulated in most biological processes, including in cellular responses to redox imbalance. In particular, miR-200 family members play a crucial role in oxidative-stress dependent endothelial dysfunction, as well as in cardiovascular complications of diabetes and obesity. In addition, different miRNAs, such as miR-210, have been demonstrated to play a key role in mitochondrial metabolism, therefore modulating ROS production and sensitivity. In this review, we will discuss miRNAs modulated by ROS or involved in ROS production, and implicated in vascular diseases in which redox imbalance has a pathogenetic role.

Aryl hydrocarbon receptor-mediated induction of the microRNA-132/212 cluster promotes interleukin-17-producing T-helper cell differentiation

Aryl hydrocarbon receptor (AHR) plays critical roles in various autoimmune diseases such as multiple sclerosis by controlling interleukin-17 (IL-17)-producing T-helper (TH17) and regulatory T cells. Although various transcription factors and cytokines have been identified as key participants in TH17 generation, the role of microRNAs in this process is poorly understood. In this study, we found that expression of the microRNA (miR)-132/212 cluster is up-regulated by AHR activation under TH17-inducing, but not regulatory T-inducing conditions. Deficiency of the miR-132/212 cluster prevented the enhancement of TH17 differentiation by AHR activation. We also identified B-cell lymphoma 6, a negative regulator of TH17 differentiation, as a potential target of the miR-212. Finally, we investigated the roles of the miR-132/212 cluster in experimental autoimmune encephalomyelitis, a murine model of multiple sclerosis. Mice deficient in the miR-132/212 cluster exhibited significantly higher resistance to the development of experimental autoimmune encephalomyelitis and lower frequencies of both TH1 and TH17 cells in draining lymph nodes. Our findings reveal a unique mechanism of AHR-dependent TH17 differentiation that depends on the miR-132/212 cluster.

Mapping the human miRNA interactome by CLASH reveals frequent noncanonical binding

MicroRNAs (miRNAs) play key roles in gene regulation, but reliable bioinformatic or experimental identification of their targets remains difficult. To provide an unbiased view of human miRNA targets, we developed a technique for ligation and sequencing of miRNA-target RNA duplexes associated with human AGO1. Here, we report data sets of more than 18,000 high-confidence miRNA-mRNA interactions. The binding of most miRNAs includes the 5' seed region, but around 60% of seed interactions are noncanonical, containing bulged or mismatched nucleotides. Moreover, seed interactions are generally accompanied by specific, nonseed base pairing. 18% of miRNA-mRNA interactions involve the miRNA 3' end, with little evidence for 5' contacts, and some of these were functionally validated. Analyses of miRNA:mRNA base pairing showed that miRNA species systematically differ in their target RNA interactions, and strongly overrepresented motifs were found in the interaction sites of several miRNAs. We speculate that these affect the response of RISC to miRNA-target binding.

New advances of microRNAs in the pathogenesis of rheumatoid arthritis, with a focus on the crosstalk between DNA methylation and the microRNA machinery

Rheumatoid arthritis (RA) is a symmetrical polyarticular disease of unknown aetiology that affects primarily the articular cartilage and bone. Characteristic features of RA pathogenesis are persistent inflammation, synovium hyperplasia and cartilage erosion accompanied by joint swelling and joint destruction. Several lines of evidence have showed a crucial role of activated fibroblast-like synoviocytes (FLS) in the pathogenesis of RA. MicroRNAs (miRNAs) are endogenous, single-stranded, non-coding RNAs with about 21 nucleotides in length and have been detected in a variety of sources, including tissues, serum, and other body fluids, such as saliva. In light of key roles of miRNAs in the regulation of gene expression, miRNAs influence a wide range of physiological and pathological processes. For example, miRNAs are evident in various malignant and nonmalignant diseases, and accumulating evidence also shows that miRNAs have important roles in the pathogenesis of RA. It has been demonstrated that miRNAs can be aberrantly expressed even in the different stages of RA progression, allowing miRNAs to help understand the pathogenesis of the disease, to act as important biomarkers, and to monitor the disease severity and the effects of drug treatment. In addition, miRNAs are emerging as potential targets for new therapeutic strategies of this kind of autoimmune disorders. The ultimate goal is the identification of miRNA targets that could be manipulated through specific therapies, aiming at activation or inhibition of specific miRNAs responsible for the RA development. In this review, the importance of miRNAs in the pathogenesis of RA is discussed systematically, with particular emphasis on the role of the crosstalk between DNA methylation and the microRNA machinery.

Context-specific microRNA function in developmental complexity

Since their discovery, microRNAs (miRNA) have been implicated in a vast array of biological processes in animals, from fundamental developmental functions including cellular proliferation and differentiation, to more complex and specialized roles such as long-term potentiation and synapse-specific modifications in neurons. This review recounts the history behind this paradigm shift, which has seen small non-coding RNA molecules coming to the forefront of molecular biology, and introduces their role in establishing developmental complexity in animals. The fundamental mechanisms of miRNA biogenesis and function are then considered, leading into a discussion of recent discoveries transforming our understanding of how these molecules regulate gene network behaviour throughout developmental and pathophysiological processes. The emerging complexity of this mechanism is also examined with respect to the influence of cellular context on miRNA function. This discussion highlights the absolute imperative for experimental designs to appreciate the significance of context-specific factors when determining what genes are regulated by a particular miRNA. Moreover, by establishing the timing, location, and mechanism of these regulatory events, we may ultimately understand the true biological function of a specific miRNA in a given cellular environment.

ROD1 is a seedless target gene of hypoxia-induced miR-210

Most metazoan microRNA (miRNA) target sites have perfect pairing to the “seed” sequence, a highly conserved region centering on miRNA nucleotides 2–7. Thus, complementarity to this region is a necessary requirement for target prediction algorithms. However, also non-canonical miRNA binding can confer target regulation. Here, we identified a seedless target of miR-210, a master miRNA of the hypoxic response. We analyzed 20 genes that were inversely correlated to miR-210 expression and did not display any complementarity with miR-210 seed sequence. We validated ROD1 (Regulator of Differentiation 1, also named PTBP3, Polypyrimidine Tract Binding protein 3) as a miR-210 seedless transcript enriched in miR-210-containing RNA-induced silencing complexes. ROD1 was not indirectly targeted by a miR-210-induced miRNA. Conversely, we identified a “centered” miR-210 binding site in ROD1 involving 10 consecutive bases in the central portion of miR-210. Reporter assays showed that miR-210 inhibited ROD1 by the direct binding to this sequence, demonstrating that ROD1 is a bona fide seedless target of miR-210. As expected, both ROD1 mRNA and protein were down-modulated upon hypoxia in a miR-210 dependent manner. ROD1 targeting by miR-210 was biologically significant: the rescue of ROD1 inhibition significantly increased hypoxia-induced cell death. These data highlight the importance of ROD1 regulation by miR-210 for cell homeostasis.

A least angle regression model for the prediction of canonical and non-canonical miRNA-mRNA interactions

microRNAs (miRNAs) are short non-coding RNAs with regulatory functions in various biological processes including cell differentiation, development and oncogenic transformation. They can bind to mRNA transcripts of protein-coding genes and repress their translation or lead to mRNA degradation. Conversely, the transcription of miRNAs is regulated by proteins including transcription factors, co-factors, and messenger molecules in signaling pathways, yielding a bidirectional regulatory network of gene and miRNA expression. We describe here a least angle regression approach for uncovering the functional interplay of gene and miRNA regulation based on paired gene and miRNA expression profiles. First, we show that gene expression profiles can indeed be reconstructed from the expression profiles of miRNAs predicted to be regulating the specific gene. Second, we propose a two-step model where in the first step, sequence information is used to constrain the possible set of regulating miRNAs and in the second step, this constraint is relaxed to find regulating miRNAs that do not rely on perfect seed binding. Finally, a bidirectional network comprised of miRNAs regulating genes and genes regulating miRNAs is built from our previous regulatory predictions. After applying the method to a human cancer cell line data set, an analysis of the underlying network reveals miRNAs known to be associated with cancer when dysregulated are predictors of genes with functions in apoptosis. Among the predicted and newly identified targets that lack a classical miRNA seed binding site of a specific oncomir, miR-19b-1, we found an over-representation of genes with functions in apoptosis, which is in accordance with the previous finding that this miRNA is the key oncogenic factor in the mir-17-92 cluster. In addition, we found genes involved in DNA recombination and repair that underline its importance in maintaining the integrity of the cell.