The overall goal of this method is to detect and precisely map microRNA (miRNA) targets in high-throughput. MiRNAs are endogenous non-coding RNAs ~22 nucleotides in length. Following transcription and processing, mature miRNAs are incorporated in a protein complex called the RNA induced silencing complex (RISC). Each miRNA guides the RISC to target elements located primarily in the 3’untranslated regions (3’UTRs) of messenger RNAs (mRNAs), resulting in either translation repression or mRNA cleavage 1. MiRNA recognize target sites based on standard Watson-Crick and G:U wobble base pairing, and are degenerate in nature, containing multiple mismatched base pairs and bulged regions. Many miRNAs are broadly conserved from plants to humans 2,3, where they play a diverse range of biological roles. In metazoans miRNAs can influence multiple biological processes including cell fate decisions 4, developmental timing 5, and frequently exhibit tissue specific expression patterns 6,7. MiRNA misexpression can also result in aberrant gene regulation, which can have substantial influence on cell behaviour based solely on the function of target genes. As such, miRNAs are linked to a wide range of diseases, including neurodegeneration 8,9, diabetes 10 and cancer 11. Bioinformatic and wet-bench approaches suggest that each miRNA may be capable of targeting hundreds to thousands of distinct mRNAs 12-14, indicating that high-throughput or genome wide approaches are required to probe this large pool of potential interactions.
Identifying target genes is a critical component of mechanistically defining miRNA function, and to do so researchers must be able to reveal targets on a large scale. Several approaches have been developed to identify miRNA targets, including bioinformatic prediction algorithms, high-throughput sequencing of targeted mRNAs, and reporter based assays. Each of these approaches has inherent strengths and weaknesses. Given that miRNA targeting is guided by sequence specificity, most notably of nucleotides 2-6 of the miRNA (termed the seed region), several algorithms have been developed to predict miRNA targets throughout the genome of many organisms. These algorithms are trained using the observed base-pairing motifs of validated miRNA targets, and frequently utilize parameters such as stringent seed pairing, site conservation, and/or thermodynamic stability 15. While these filters refine the large number of putative targets with sufficient complementarity to only high confidence targets, they may exclude species specific and non-canonical miRNA target sites, which recent evidence suggests are widespread 16-24. Furthermore, these predictions do not take into account mechanisms of mRNA processing that exclude miRNA target sites, such as alternative polyadenylation 25, RNA editing 26, RNA methylation 27, and cooperative binding. As such, high false positive and false negative rates have been reported for many algorithms 22,24,28. While these algorithms are useful to identify candidate miRNA targets for subsequent experimental validation, these high error rates limit the efficacy of bioinformatic approaches for systematic miRNA target detection.
To systematically probe for interactions between a given miRNA and potentially targeted 3’UTRs we have developed a high-throughput assay called Luminescent Identification of Functional Elements in 3’UTRs (3’LIFE) 24. This assay measures direct interactions and translational repression of the test 3’UTR by a query miRNA using a dual luciferase reporter system. In this system, the 3’UTR of a gene of interest is cloned downstream of the firefly luciferase (fluc) reporter reading frame. The reporter construct is cotransfected with a query miRNA in HEK293T cells. MiRNA targeting is determined by measuring the relative change between the test fluc::3’UTR reporter and a second non-specific Renilla luciferase reporter. Importantly, luciferase assays detect functional miRNA/mRNA interactions that influence the translational output of the reporter. This is a key advantage over traditional methods to detect miRNA regulation, such as RT-qPCR and Western blots, in that this bypasses differences in mRNA degradation and translational repression, as well as changes in protein abundance independent of 3’UTR based regulation.
Luciferase assays are widely utilized to validate direct miRNA targets because of their relative simplicity and sensitivity, yet their use in high-throughput screens is limited by high costs associated with consumable reagents, the lack of 3’UTR libraries from public sources, and the absence of standardized luciferase protocols, leading to difficulties in comparing functional repression across multiple datasets. To facilitate the use of the 3’LIFE assay, we have placed emphasis on simplification of experimental design, utilization of non-commercial transfection 24 and luciferase reagents 29, creating a 3’UTR library which is regularly updated and expanded, and is available through a public plasmid repository 30.
The scalability of the 3’LIFE assay allows screening of a large 3’UTR library for targeting by a given miRNA without biasing the screen towards bioinformatically identified genes. In addition to testing canonical and predicted interactions, this systematic approach allows the identification of novel targets driven via non-canonical and/or species-specific interactions. Importantly, the effect of miRNA targeting on protein production is generally understood to result in modest translational repression 15,31, suggesting that a primary role of miRNA regulation is to fine-tune protein output, protect against aberrant levels of gene expression, and provide robustness to cell specific programs 32,33. The sensitivity of the luciferase assay combined with the inherently large number of negative miRNA/mRNA interactions in the 3’LIFE screen allows the detection of subtle effects of miRNA targeting on a large number of genes, and the identification of multiple components of gene networks that are regulated by a given miRNA 24.
Here we describe the 3'LIFE protocol, and demonstrate it’s feasibility by screening two well characterized miRNAs, miR-10b and let-7c against a panel of 275 human 3'UTRs (Figure 1).