Large-Scale Screens of miRNA-mRNA Interactions Unveiled That the 3′UTR of a Gene Is Targeted by Multiple miRNAs

Animal microRNA (miRNA) target prediction is still a challenge, although many prediction programs have been exploited. MiRNAs exert their function through partially binding the messenger RNAs (mRNAs; likely at 3′ untranslated regions [3′UTRs]), which makes it possible to detect the miRNA-mRNA interactions in vitro by co-transfection of miRNA and a luciferase reporter gene containing the target mRNA fragment into mammalian cells under a dual-luciferase assay system. Here, we constructed a human miRNA expression library and used a dual-luciferase assay system to perform large-scale screens of interactions between miRNAs and the 3′UTRs of seven genes, which included more than 3,000 interactions with triplicate experiments for each interaction. The screening results showed that the 3′UTR of one gene can be targeted by multiple miRNAs. Among the prediction algorithms, a Bayesian phylogenetic miRNA target identification algorithm and a support vector machine (SVM) presented a relatively better performance (27% for EIMMo and 24.7% for miRDB) against the average precision (17.3%) of the nine prediction programs used here. Additionally, we noticed that a relatively high conservation level was shown at the miRNA 3′ end targeted regions, as well as the 5′ end (seed region) binding sites.     

A non-canonical plant microRNA target site

Plant microRNAs (miRNAs) typically form near-perfect duplexes with their targets and mediate mRNA cleavage. Here, we describe an unconventional miRNA target of miR398 in Arabidopsis, an mRNA encoding the blue copper-binding protein (BCBP). BCBP mRNA carries an miR398 complementary site in its 5′-untranslated region (UTR) with a bulge of six nucleotides opposite to the 5′ region of the miRNA. Despite the disruption of a target site region thought to be especially critical for function, BCBP mRNAs are cleaved by ARGONAUTE1 between nucleotides 10th and 11th, opposite to the miRNA, like conventional plant target sites. Levels of BCBP mRNAs are inversely correlated to levels of miR398 in mutants lacking the miRNA, or transgenic plants overexpressing it. Introducing two mutations that disrupt the miRNA complementarity around the cleavage site renders the target cleavage-resistant. The BCBP site functions outside of the context of the BCBP mRNA and does not depend on 5′-UTR location. Reducing the bulge does not interfere with miR398-mediated regulation and completely removing it increases the efficiency of the slicing. Analysis of degradome data and target predictions revealed that the miR398-BCBP interaction seems to be rather unique. Nevertheless, our results imply that functional target sites with non-perfect pairings in the 5′ region of an ancient conserved miRNA exist in plants.     

Uncoupling of lin-14 mRNA and protein repression by nutrient deprivation in Caenorhabditis elegans

In animals, microRNAs (miRNAs), typically, pair to sites of partial complementarity in the 3′-untranslated regions (3′UTRs) of target genes. Regulation by miRNAs often results in down-regulation of target mRNA and protein expression by mechanisms that are yet to be fully elucidated. Additionally, changes in environmental conditions have been shown to influence miRNA function in some cell culture systems. Here, we report the effect of nutrient deprivation on regulation of an endogenous miRNA target in developing worms. In Caenorhabditis elegans, the lin-4 miRNA recognizes multiple sites in the lin-14 3′UTR and directs mRNA degradation and translational repression, but it is unclear how these processes are coupled. In this study, we demonstrate that nutrient deprivation results in loss of lin-14 mRNA, but not protein, repression. In worms removed from feeding conditions, lin-14 mRNA reaccumulates despite the continued expression of lin-4 miRNA. The relative increase in lin-14 mRNA levels during nutrient deprivation is less pronounced in genetic mutants lacking lin-4 miRNA or the lin-14 3′UTR target sites. In conclusion, regulation of lin-14 at the mRNA and protein levels can be uncoupled by changes in culture conditions, indicating that miRNA function can be modulated by environment in multicellular organisms. The awareness that endogenous miRNA pathways can be sensitive to environment is an important consideration for elucidating the mechanism used by miRNAs to regulate target mRNA and protein expression.     

The microRNAs in an Ancient Protist Repress the Variant-Specific Surface Protein Expression by Targeting the Entire Coding Sequence

microRNAs (miRNA) have been detected in the deeply branched protist, Giardia lamblia, and shown to repress expression of the family of variant-specific surface proteins (VSPs), only one of which is expressed in Giardia trophozoite at a given time. Three next-generation sequencing libraries of Giardia Argonaute-associated small RNAs were constructed and analyzed. Analysis of the libraries identified a total of 99 new putative miRNAs with a size primarily in the 26 nt range similar to the size previously predicted by the Giardia Dicer crystal structure and identified by our own studies. Bioinformatic analysis identified multiple putative miRNA target sites in the mRNAs of all 73 VSPs. The effect of miRNA target sites within a defined 3′-region were tested on two vsp mRNAs. All the miRNAs showed partial repression of the corresponding vsp expression and were additive when the targeting sites were separately located. But the combined repression still falls short of 100%. Two other relatively short vsp mRNAs with 15 and 11 putative miRNA target sites identified throughout their ORFs were tested with their corresponding miRNAs. The results indicate that; (1) near 100% repression of vsp mRNA expression can be achieved through the combined action of multiple miRNAs on target sites located throughout the ORF; (2) the miRNA machinery could be instrumental in repressing the expression of vsp genes in Giardia; (3) this is the first time that all the miRNA target sites in the entire ORF of a mRNA have been tested and shown to be functional.     

Tissue and Process Specific microRNA–mRNA Co-Expression in Mammalian Development and Malignancy

An association between enrichment and depletion of microRNA (miRNA) binding sites, 3′ UTR length, and mRNA expression has been demonstrated in various developing tissues and tissues from different mature organs; but functional, context-dependent miRNA regulations have yet to be elucidated. Towards that goal, we examined miRNA–mRNA interactions by measuring miRNA and mRNA in the same tissue during development and also in malignant conditions. We identified significant miRNA-mediated biological process categories in developing mouse cerebellum and lung using non-targeted mRNA expression as the negative control. Although miRNAs in general suppress target mRNA messages, many predicted miRNA targets demonstrate a significantly higher level of co-expression than non-target genes in developing cerebellum. This phenomenon is tissue specific since it is not observed in developing lungs. Comparison of mouse cerebellar development and medulloblastoma demonstrates a shared miRNA–mRNA co-expression program for brain-specific neurologic processes such as synaptic transmission and exocytosis, in which miRNA target expression increases with the accumulation of multiple miRNAs in developing cerebellum and decreases with the loss of these miRNAs in brain tumors. These findings demonstrate the context-dependence of miRNA–mRNA co-expression.     

miRTar Hunter: A prediction system for identifying human microRNA target sites

MicroRNAs (miRNAs) are important regulators of gene expression and play crucial roles in many biological processes including apoptosis, differentiation, development, and tumorigenesis. Recent estimates suggest that more than 50% of human protein coding genes may be regulated by miRNAs and that each miRNA may bind to 300–400 target genes. Approximately 1,000 human miRNAs have been identified so far with each having up to hundreds of unique target mRNAs. However, the targets for a majority of these miRNAs have not been identified due to the lack of large-scale experimental detection techniques. Experimental detection of miRNA target sites is a costly and time-consuming process, even though identification of miRNA targets is critical to unraveling their functions in various biological processes. To identify miRNA targets, we developed miRTar Hunter, a novel computational approach for predicting target sites regardless of the presence or absence of a seed match or evolutionary sequence conservation. Our approach is based on a dynamic programming algorithm that incorporates more sequence-specific features and reflects the properties of various types of target sites that determine diverse aspects of complementarities between miRNAs and their targets. We evaluated the performance of our algorithm on 532 known human miRNA:target pairs and 59 experimentally-verified negative miRNA:target pairs, and also compared our method with three popular programs for 481 miRNA:target pairs. miRTar Hunter outperformed three popular existing algorithms in terms of recall and precision, indicating that our unique scheme to quantify the determinants of complementary sites is effective at detecting miRNA targets. miRTar Hunter is now available at http://203.230.194.162/~kbkim.     

MREdictor: a two-step dynamic interaction model that accounts for mRNA accessibility and Pumilio binding accurately predicts microRNA targets

The prediction of pairing between microRNAs (miRNAs) and the miRNA recognition elements (MREs) on mRNAs is expected to be an important tool for understanding gene regulation. Here, we show that mRNAs that contain Pumilio recognition elements (PRE) in the proximity of predicted miRNA-binding sites are more likely to form stable secondary structures within their 3′-UTR, and we demonstrated using a PUM1 and PUM2 double knockdown that Pumilio proteins are general regulators of miRNA accessibility. On the basis of these findings, we developed a computational method for predicting miRNA targets that accounts for the presence of PRE in the proximity of seed-match sequences within poorly accessible structures. Moreover, we implement the miRNA-MRE duplex pairing as a two-step model, which better fits the available structural data. This algorithm, called MREdictor, allows for the identification of miRNA targets in poorly accessible regions and is not restricted to a perfect seed-match; these features are not present in other computational prediction methods.     

Potent degradation of neuronal miRNAs induced by highly complementary targets

MicroRNAs (miRNAs) regulate target mRNAs by silencing them. Reciprocally, however, target mRNAs can also modulate miRNA stability. Here, we uncover a remarkable efficacy of target RNA-directed miRNA degradation (TDMD) in rodent primary neurons. Coincident with degradation, and while still bound to Argonaute, targeted miRNAs are 3′ terminally tailed and trimmed. Absolute quantification of both miRNAs and their decay-inducing targets suggests that neuronal TDMD is multiple turnover and does not involve co-degradation of the target but rather competes with miRNA-mediated decay of the target. Moreover, mRNA silencing, but not TDMD, relies on cooperativity among multiple target sites to reach high efficacy. This knowledge can be harnessed for effective depletion of abundant miRNAs. Our findings bring insight into a potent miRNA degradation pathway in primary neurons, whose TDMD activity greatly surpasses that of non-neuronal cells and established cell lines. Thus, TDMD may be particularly relevant for miRNA regulation in the nervous system.     

The Impact of 3′UTR Variants on Differential Expression of Candidate Cancer Susceptibility Genes

Variants in regulatory regions are predicted to play an important role in disease susceptibility of common diseases. Polymorphisms mapping to microRNA (miRNA) binding sites have been shown to disrupt the ability of miRNAs to target genes resulting in differential mRNA and protein expression. Skin tumor susceptibility 5 (Skts5) was identified as a locus conferring susceptibility to chemically-induced skin cancer in NIH/Ola by SPRET/Outbred F1 backcrosses. To determine if polymorphisms between the strains which mapped to putative miRNA binding sites in the 3′ untranslated region (3′UTR) of genes at Skts5 influenced expression, we conducted a systematic evaluation of 3′UTRs of candidate genes across this locus. Nine genes had polymorphisms in their 3′UTRs which fit the linkage data and eight of these contained polymorphisms suspected to interfere with or introduce miRNA binding. 3′UTRs of six genes, Bcap29, Dgkb, Hbp1, Pik3cg, Twistnb, and Tspan13 differentially affected luciferase expression, but did not appear to be differentially regulated by the evaluated miRNAs predicted to bind to only one of the two isoforms. 3′UTRs from four additional genes chosen from the locus that fit less stringent criteria were evaluated. Ifrd1 and Etv1 showed differences and contained polymorphisms predicted to disrupt or create miRNA binding sites but showed no difference in regulation by the miRNAs tested. In summary, multiple 3′UTRs with putative functional variants between susceptible and resistant strains of mice influenced differential expression independent of predicted miRNA binding.     

Anti-Argonaute RIP-Chip shows that miRNA transfections alter global patterns of mRNA recruitment to microribonucleoprotein complexes

MicroRNAs (miRNAs) play key roles in gene expression regulation by guiding Argonaute (AGO)-containing microribonucleoprotein (miRNP) effector complexes to target polynucleotides. There are still uncertainties about how miRNAs interact with mRNAs. Here we employed a biochemical approach to isolate AGO-containing miRNPs from human H4 tumor cells by co-immunoprecipitation (co-IP) with a previously described anti-AGO antibody. Co-immunoprecipitated (co-IPed) RNAs were subjected to downstream Affymetrix Human Gene 1.0 ST microarray analysis. During rigorous validation, the “RIP-Chip” assay identified target mRNAs specifically associated with AGO complexes. RIP-Chip was performed after transfecting brain-enriched miRNAs (miR-107, miR-124, miR-128, and miR-320) and nonphysiologic control miRNA to identify miRNA targets. As expected, the miRNA transfections altered the mRNA content of the miRNPs. Specific mRNA species recruited to miRNPs after miRNA transfections were moderately in agreement with computational target predictions. In addition to recruiting mRNA targets into miRNPs, miR-107 and to a lesser extent miR-128, but not miR-124 or miR-320, caused apparent exclusion of some mRNAs that are normally associated with miRNPs. MiR-107 and miR-128 transfections also result in decreased AGO mRNA and protein levels. However, AGO mRNAs were not recruited to miRNPs after either miR-107 or miR-128 transfection, confirming that miRNAs may alter gene expression without stable association between particular mRNAs and miRNPs. In summary, RIP-Chip assays constitute an optimized, validated, direct, and high-throughput biochemical assay that provides data about specific miRNA:mRNA interactions, as well as global patterns of regulation by miRNAs.