Mol. Cells

Determinants of Functional MicroRNA Targeting

Hyeonseo Hwang, Hee Ryung Chang

Additional article information

Abstract

MicroRNAs (miRNAs) play cardinal roles in regulating biological pathways and processes, resulting in significant physiological effects. To understand the complex regulatory network of miRNAs, previous studies have utilized massive-scale datasets of miRNA targeting and attempted to computationally predict the functional targets of miRNAs. Many miRNA target prediction tools have been developed and are widely used by scientists from various fields of biology and medicine. Most of these tools consider seed pairing between miRNAs and their mRNA targets and additionally consider other determinants to improve prediction accuracy. However, these tools exhibit limited prediction accuracy and high false positive rates. The utilization of additional determinants, such as RNA modifications and RNA-binding protein binding sites, may further improve miRNA target prediction. In this review, we discuss the determinants of functional miRNA targeting that are currently used in miRNA target prediction and the potentially predictive but unappreciated determinants that may improve prediction accuracy.

Keywords: bioinformatics, microRNA, microRNA target prediction, microRNA targeting, microRNA targeting determinants

CURRENT STATE OF FUNCTIONAL miRNA TARGET PREDICTION

MicroRNAs (miRNAs) regulate a broad range of biological processes and physiological pathways, including tumor suppression and progression (Gregory et al., 2008; Jang et al., 2020; O'Donnell et al., 2005; Peng and Croce, 2016), immune cell development and function (Fontana et al., 2007; Han et al., 2020b; Mehta and Baltimore, 2016; Muljo et al., 2005; O'Connell et al., 2010; Rodriguez et al., 2007; Thai et al., 2007), cardiovascular diseases (Care et al., 2007; Harris et al., 2008), neural development and function (Brennan and Henshall, 2020; Ghosh et al., 2014; Giraldez et al., 2005; Krichevsky et al., 2006; Tan et al., 2013), early embryonic development (Bernstein et al., 2003; Choi et al., 2007; Wienholds et al., 2005), and cytoskeletal dynamics (Fededa et al., 2016; Wu et al., 2014; Xin et al., 2009). Since miRNAs modulate gene expression by binding to their targets, it is crucial to systematically identify and evaluate all the functional targets of miRNAs to better understand complex miRNA regulatory networks (Bracken et al., 2016; Plaisier et al., 2012; Pu et al., 2019).

The generation of a comprehensive list of miRNA targets is challenging because miRNA targets are abundant and widespread across the transcriptome (Friedman et al., 2009; Lim et al., 2005; Selbach et al., 2008). The classical strategy for studying miRNA targeting involves the identification of potential targets using experimental methods, including microarray (Lim et al., 2005), RIP (RNA immunoprecipitation) (Keene et al., 2006), cross-linking immunoprecipitation (CLIP) (Chi et al., 2009), miRNA pull-down (Hassan et al., 2013; Orom and Lund, 2010), and luciferase reporter assays (Ghanbarian et al., 2022; Thomson et al., 2011; Tuschl et al., 1999). The generation of massive amounts of data through these experiments is often limited by a prohibitively large amount of labor, time, and cost. Alternatively, functional miRNA targeting can be identified using computational tools that are much faster and more accessible than conventional experimental methods.

Recently, many computational tools for miRNA targeting prediction have been developed (Nachtigall and Bovolenta, 2022). However, the predictive power of these tools requires further improvement. For instance, even with the state-of-the-art prediction tool (Agarwal et al., 2015; McGeary et al., 2019), the predicted repressive effects show a weak correlation (r2 ≤ 0.2) with the experimental data (Agarwal et al., 2015; McGeary et al., 2019). Here, we review the determinants of functional miRNA targeting that have not been used in computational miRNA target prediction tools. The incorporation of previously unappreciated determinants may improve the performance of miRNA target prediction tools.

PRINCIPLES OF miRNA TARGETING

Mature miRNAs form a protein-RNA complex, termed the RNA-induced silencing complex (RISC), with the Argonaute (AGO) protein (Gregory et al., 2005; Hammond et al., 2000). The current RISC model suggests that four nucleotides (nucleotides 2-5 of guide RNA; g2-5) at the 5′ end of the miRNA are exposed before binding to the target (Schirle et al., 2014). The pairing between this region of a miRNA and its target mRNA induces a conformational change in AGO and exposes additional nucleotides (g6-8 and 13-16) that allow stable pairing between RISC and its target. The sequence of 6 or 7 nucleotides (g2-7 or g2-8) at the 5′ end of the miRNAs, called the seed region, primarily determines the target specificity of the miRNAs (Bartel, 2009). The nucleotides exposed at the 3′ ends of miRNAs may pair with the target and supplement seed pairing (Grimson et al., 2007).

The primary determinant of miRNA targeting efficacy is pairing in the seed region. Mismatches in the seed region reduce the potency of the target site. The miRNA target sites are classified into 'canonical' sites that contain a perfect hexamer match with the seed (g2-7), and the ‘noncanonical’ sites that match imperfectly with the seed (Bartel, 2009) (Fig. 1). The canonical miRNA target sites are further divided into four types according to two criteria: a t8 match with the g8 nucleotide and the presence of adenine at position t1. The 8-mer sites with both a t8 match and a t1A match are the most effective, followed by the 7-mer m8, 7-mer A1, and 6-mer sites (Baek et al., 2008; Bartel, 2009; Lewis et al., 2005). Noncanonical sites generally confer weaker repressive effects than canonical sites, and a subset of the noncanonical sites are only effective in favorable contexts, hence the term 'context-dependent noncanonical site types (CDNSTs)’ (Kim et al., 2016). Currently, offset 7-mer (g3-8 match), offset 6-mer (g3-7 match), 6-mer A1 (g2-6 match with t1A), and four other CDNSTs have been identified as functional noncanonical sites through large-scale bioinformatics analysis (Kim et al., 2016). Many widely used miRNA target prediction tools, such as TargetScan (Agarwal et al., 2015), do not fully account for noncanonical sites.

Figure F1
Structures of guide RNA-loaded human Argonaute protein 2 (hAGO2) before (top left; PDB ID: 4W5N) and after (top right; PDB ID: 4W5O) binding to the target RNA are shown (Schirle ...

DETERMINANTS OF miRNA TARGETING PREVIOUSLY UTILIZED FOR TARGET PREDICTION

Software tools for miRNA targeting prediction

Each miRNA target prediction tool utilizes a different set of determinants (Table 1). Earlier tools, including PITA (Kertesz et al., 2007), MicroTar (Thadani and Tammi, 2006), RNAhybrid (Kruger and Rehmsmeier, 2006), and PicTar (Krek et al., 2005) rely only on a few determinants, whereas more recent tools, such as TargetScan (Agarwal et al., 2015; McGeary et al., 2019), MIRZA-G (Gumienny and Zavolan, 2015), DIANA-microT-CDS (Reczko et al., 2012), and miRanda-mirSVR (Betel et al., 2010) utilize a larger number of determinants. Most of these tools use seed pairing, thermodynamic stability, evolutionary conservation, and the structural accessibility of target sites. Of the currently available miRNA targeting prediction tools, TargetScan is considered to be the most accurate tool that utilizes the largest number of determinants in the model. Here, we summarize the determinants of miRNA targeting used in these tools (Fig. 2).

Figure F2
Determinants of miRNA targeting that have already been utilized for target prediction include 3′UTR length, minimum distance to the 3′UTR end, structural accessibility, local AU content, evolutionary conservation, ORF target ...
Table 1

Local AU content

Local AU content around the target site is one of the most prominent determinants of miRNA targeting. Effective and conserved target sites are associated with higher local AU content but not global AU content (Grimson et al., 2007; Nielsen et al., 2007). Although the mechanism by which local AU content affects miRNA targeting remains unclear, many target prediction tools, including TargetScan7, miRanda-mirSVR, and DIANA-microT-CDS, utilize this feature. Two hypotheses that are not mutually exclusive have been proposed to explain the correlation between stronger repression and higher local AU content. The first hypothesis is that RISC interacts with AREs (AU-rich elements), and the repressive efficacy of the nearby target is enhanced by this interaction (Jing et al., 2005; Nielsen et al., 2007; Vasudevan and Steitz, 2007). An alternative hypothesis is that high AU content inhibits the formation of secondary structures and increases the accessibility of the target, leading to stronger repression (Grimson et al., 2007; Nielsen et al., 2007).

Structural accessibility

The structural accessibility of the target site, mainly determined by the secondary or tertiary structure of the 3′UTR, positively correlates with targeting efficacy (Kertesz et al., 2007). Higher site accessibility increases the chance of RISC binding to its target, leading to stronger repression. The structural accessibility of the target site can be measured as ΔG or ΔΔG, which is the difference in ΔG between the folded and open 3′UTRs. A previous study demonstrated that ΔΔG had a stronger correlation with the repressive effect than ΔG of the folded 3′UTR (Kertesz et al., 2007). Another study argued that structural accessibility, measured by ΔG, lost its correlation with repression when local AU was controlled (Grimson et al., 2007). A more recent study showed that when structural accessibility was measured in log unpaired probability using RNAplFold (Bernhart et al., 2006), the correlation was present even when other confounding factors, including local AU, were controlled (Agarwal et al., 2015).

Pairing stability and target abundance

The thermodynamic pairing stability between a miRNA and its target is a determinant widely used for miRNA targeting prediction. There are two approaches to calculating pairing stability: (1) using the whole miRNA sequence and (2) using only the seed region of the miRNA, which is termed seed pairing stability. Many miRNA target prediction tools, including DIANA-microT-CDS, MIRZA-G, PITA, PicTar, RNAhybrid, and MicroTar, rely heavily on the predicted thermodynamic pairing stability between the entire mature miRNA and target. However, TargetScan7 utilizes seed pairing stability as a feature (Agarwal et al., 2015). Seed pairing stability is negatively correlated with the abundance of the target site in the 3′UTR of mRNAs (Garcia et al., 2011). Targets with a higher abundance are associated with weaker repressive effects, which may be due to competition between miRNA targets (Arvey et al., 2010). Despite the correlation between target abundance and seed pairing stability, both determinants affect miRNA targeting efficacy independently and globally (Garcia et al., 2011).

Binding affinity

Recently, RNA Bind-n-Seq (RBNS), a technique that measures the affinity between a protein and random RNA sequences by pull-down of RNA-binding protein (RBP) with bound RNA followed by high-throughput sequencing (Lambert et al., 2014), has been employed to measure the binding affinity between miRNA-loaded AGO and 12-nt RNA sequences (McGeary et al., 2019). The AGO-RBNS data have been used to develop a biochemical model of miRNA targeting and are included in TargetScan8 (McGeary et al., 2019).

3′ pairing

When RISC opens up after the initial pairing between the miRNA and the target RNA, two regions of the miRNA open up: nucleotides 2-8 and 13-16 (Schirle et al., 2014). While miRNA targeting is primarily determined by the former, the latter region also affects the repressive effect of miRNAs (Friedman et al., 2009; Grimson et al., 2007). The pairing of the 3′ regions of the miRNA, including nucleotides 13-16, may supplement the canonical seed match or compensate for the noncanonical seed match (Bartel, 2009). About 5% of the conserved target sites have supplementary 3′ pairings, and less than 2% have compensatory 3′ pairings (Friedman et al., 2009). A recent study (McGeary et al., 2022) reported that each miRNA has different preferences for an optimal 3′ pairing, while some miRNAs, such as let-7, can have two 3′ binding modes with different offsets, and nucleotides outside of positions 13-16 could have a significant impact on the 3′ binding affinities.

ORF sites

Although the majority of functional miRNA target sites are located in the 3′UTRs, targets in the ORFs may have weak repressive effects. ORF sites can strengthen the repressive effect of 3′UTR sites in a synergistic or additive manner (Fang and Rajewsky, 2011). The results of DIANA-microT-CDS suggest that the incorporation of ORF site as a determinant may aid miRNA targeting prediction (Reczko et al., 2012).

Distance to 3′UTR ends

Target sites closer to either end of the 3′UTR are associated with higher repressive efficacy (Gaidatzis et al., 2007; Grimson et al., 2007; Majoros and Ohler, 2007). There are two hypotheses to explain the correlation between repression and the minimum distance to the 3′UTR. The first hypothesis is that the target sites near either end of the 3′UTR could be closer to the translation or RNA-processing complexes (Gaidatzis et al., 2007; Grimson et al., 2007). At the 3′ end, mRNA looping may bring the site closer to the protein complex (Grimson et al., 2007). Another possibility is that the sites near the ends have higher structural accessibility because ribosomes and poly(A)-binding proteins inhibit secondary structure formation (Grimson et al., 2007; Majoros and Ohler, 2007). This hypothesis is congruent with a recent study that showed that RBPs that bind near the target site increase targeting efficacy by opening the RNA secondary structure (Kim et al., 2021). While sites closer to the ends of 3′UTRs confer higher targeting efficacies, studies have also reported that conserved miRNA target sites are relatively scarce in the 15- to 20-nt region after the stop codon, likely due to interference by the translation machinery (Gaidatzis et al., 2007; Grimson et al., 2007; Majoros and Ohler, 2007).

3′UTR length

The role of 3′UTR length as a determinant of functional miRNA targeting has been suggested by previous studies (Hausser et al., 2009; Hong et al., 2009; Sandberg et al., 2008; Stark et al., 2005), but the validity and direction of the correlation between 3′UTR length and targeting efficacy remain open to debate (Wen et al., 2011). Nevertheless, regression analyses from recent studies have shown that the 3′UTR length has a significant negative correlation with targeting efficacy (Agarwal et al., 2015; Betel et al., 2010). Three hypotheses have been proposed to explain this correlation. First, a previous study showed a correlation between shorter 3′UTRs and higher structural accessibility through an association analysis, which could improve miRNA targeting (Hong et al., 2009). Second, shorter 3′UTRs have higher target site densities (Hong et al., 2009), which may increase the chance of cooperative action among adjacent target sites (Briskin et al., 2020; Grimson et al., 2007; Saetrom et al., 2007). Third, the higher number of endogenous miRNA target sites in longer 3′UTRs can induce stronger local competition between target sites, leading to derepression and decreased proficiency of target sites for exogenous miRNA (Kim et al., 2014).

Evolutionary conservation

Functionally important regulatory elements are likely to be conserved during evolution, and the same principle applies to miRNA target sites. Because miRNA target sites are frequently under negative selection pressure, more than 60% of human protein-coding genes contain conserved target sites (Friedman et al., 2009). Notably, the degree of target site conservation is positively correlated with a repressive effect (Friedman et al., 2009; Nielsen et al., 2007). Most widely used miRNA prediction tools, including TargetScan, miRanda-mirSVR, MIRZA-G, PicTar, and DIANA-microT-CDS, employ evolutionary conservation as a determinant to predict miRNA targeting.

Affected isoform ratio (AIR)

Fifty-four percent of human genes have multiple polyadenylation (poly(A)) sites, and 51% of human poly(A) sites have heterogeneous cleavage sites (Tian et al., 2005), resulting in transcript isoforms with varying 3′UTR lengths. Such isoforms can confound the 3′UTR-related determinants of miRNA targeting, specifically the distance to the 3′UTR ends and the 3′UTR length (Agarwal et al., 2015). The repressive effect of the site can be reduced if the target site is present in only some of the isoforms, as shown by a study that showed that the isoforms containing the target sites are more strongly repressed by miRNAs than those without the target sites (Legendre et al., 2006). Moreover, the difference in miRNA targeting between cell types has been attributed to the 3′UTR isoform composition (Nam et al., 2014). To account for the 3′UTR isoforms, poly(A)-position profiling by sequencing (3P-seq) (Jan et al., 2011) can be used. TargetScan incorporated 3P-seq by calculating the AIR, which is defined as a fraction of the isoforms that contain a specific target site (Agarwal et al., 2015; Nam et al., 2014).

DETERMINANTS OF miRNA TARGETING THAT HAVE NOT BEEN EMPLOYED IN TARGET PREDICTION

RNA modifications (RMs) in primary miRNAs that affect targeting

RNAs can be chemically modified in more than 100 different ways, and these RMs play crucial roles in the regulation of both coding and noncoding RNAs (Roundtree et al., 2017). It is widely accepted that RMs, such as N6-methyladenosine (m6A), of miRNA precursors are essential for proper processing and maturation (Alarcon et al., 2015b). Recent reports have suggested that various RMs in miRNAs alter their targeting properties and may have significant physiological effects.

One of the first RMs reported to affect miRNA targeting is adenosine-to-inosine (A-to-I) editing. A-to-I editing is induced by adenosine deaminase acting on the RNA (ADAR) family, both in the nucleus and cytoplasm (Bass and Weintraub, 1988; Wagner et al., 1989). Primary miRNAs (pri-miRNAs) can undergo A-to-I editing via ADARs (Blow et al., 2006; Luciano et al., 2004). While A-to-I editing of pri-miRNAs can interfere with miRNA processing (Yang et al., 2006), some edited nucleotides are retained in mature miRNAs and can also affect miRNA targeting (Kawahara et al., 2007). The effect of A-to-I editing on miRNA targeting was first shown for pri-miR-376a, in which single A-to-I editing in the middle of the seed drastically redirected its target specificity (Kawahara et al., 2007). Only 2 of the 78 predicted targets were retained after A-to-I editing of miR-376a. One of the redirected targets of the edited miR-376a is PRPS1, which encodes an enzyme in the uric acid synthesis pathway. Increased levels of uric acid in the brain cortex of ADAR2 knockout mice demonstrated the biological significance of A-to-I editing of miR-376a. In addition, A-to-I editing of the seed sequence of pri-miR-589-3p by ADAR2 has been reported to redirect its targeting, resulting in the inhibition of glioblastoma progression, in part due to the suppression of ADAM12 expression (Cesarini et al., 2018). The redirection of miRNA targeting by A-to-I editing has been reported for at least 14 different miRNAs (Nishikura, 2016). Although A-to-I editing is considered similar to A-to-G substitution because inosine preferentially pairs with cytidine, one study showed that the effect of A-to-I editing of the miRNA seed region on miRNA targeting is different from the effect of A-to-G substitutions (Kume et al., 2014).

RMs present in mature miRNAs

In addition to A-to-I editing, three other classes of RMs present in mature RNAs have been found to affect miRNA targeting. m6A methylation is one of the most prevalent types of RM (Dominissini et al., 2012). m6A modifications are made by writer proteins, including the methyltransferase-like (METTL) family of proteins (Liu et al., 2014), and recognized by reader proteins, including the YTH domain-containing family of proteins (Liao et al., 2018). Although the role of m6A in miRNA biogenesis is well known (Alarcon et al., 2015a; 2015b), its role in miRNA targeting has not been systematically evaluated. A recent study found that the m6A modification of mature miR-200c-3p diminishes its repressive function. Structural prediction of RISC showed that a single m6A modification of the seed of let-7a-5p is sufficient to globally modify the 3-D structure of RISC (Konno et al., 2019).

When guanine is oxidized to 8-oxoguanine (o8G), it can pair with either adenine or cytosine (Michaels et al., 1992). Recent studies have shown that o8G modification of miRNA seeds can redirect its targets, which may result in significant physiological outcomes. For miR-1, when activated adrenergic receptors induce o8G modification in the seed, its targeting is globally redirected (Seok et al., 2020). While miR-1 is known to cause atrophy (Li et al., 2010), the 7o8G miR-1 induces cardiac hypertrophy. It was also reported that the o8G modification of miR-184 redirects the miRNA to Bcl-xL and Bcl-w, key anti-apoptotic genes that are not targeted by unoxidized miR-184 (Wang et al., 2015). Oxidized miR-184 promotes apoptosis and myocardial infarction, highlighting the clinical significance of o8G present in mature miRNAs.

Another abundant type of RM is 5-methylcytosine (m5C), installed by enzymes in the NOL1/NOP2/SUN domain (NSUN) family and DNA methyltransferase 2 (DNMT2) (Bohnsack et al., 2019). A recent study showed that the m5C modification of mature miR-181a-5p abolished its tumor suppressor function through the derepression of Bcl-2-like protein 11 (BIM), a key protein in the initiation of apoptosis (Cheray et al., 2020). In glioblastoma patients, m5C of miR-181a-5p was associated with a worse prognosis, emphasizing the potential physiological importance of m5C modification of mature miRNAs. The study also found that the DNMT3A/AGO4 complex was responsible for m5C modification of miRNAs.

RMs in target mRNAs

RMs at miRNA target sites and adjacent regions affect miRNA targeting. The binding of the m6A reader protein IGF2BP1 near the miRNA target site in the SRF transcript hinders the binding of miR-2 and miR-125 in an m6A-dependent manner (Muller et al., 2019). In contrast, the binding of another m6A reader protein, IGF2BP2, to the target site of miR-133a increases the repressive efficacy of the miRNA by physically interacting with AGO (Qian et al., 2021). m5C is also suggested to play a role in miRNA targeting based on the enrichment of its putative sites in miRNA target sites (Squires et al., 2012). A-to-I editing in the target sites of miR-30b-3p and miR-573 in the 3′UTR of the ARHGAP26 transcript blocks miRNAs from pairing to the target (Wang et al., 2013).

RNA-binding proteins

It is widely accepted that miRNA target sites with higher structural accessibility elicit stronger miRNA-mediated repression (Agarwal et al., 2015). RBPs can open the secondary structures of the 3′UTRs and hence make the target site more accessible for RISC. This phenomenon has been observed for PCBP2 (Lin et al., 2016) and Pumilio (Kedde et al., 2010), inducing a local structural change in the target 3′UTR and thus enhancing miRNA targeting. This hypothesis is supported by a recent study that showed that the number of RBP-binding sites and the proximity of RBP-binding sites to miRNA target sites are positively correlated with increased repression (Kim et al., 2021). Therefore, the locations of the RBP-binding sites in the 3′UTR, which can be accurately determined by enhanced CLIP (eCLIP) (Van Nostrand et al., 2016), could improve miRNA targeting prediction.

Individual RBPs can affect miRNA targeting via mechanisms other than increasing structural accessibility. PTB, an RBP that binds to the polypyrimidine tract, can compete with RISC to bind to miRNA target sites and inhibit miRNA targeting (Xue et al., 2013). HuR, an RBP known to regulate mRNA stability, competes with miR-125b when its binding site is adjacent to its miRNA target site (Ahuja et al., 2016). A recent study showed that UPF1, an essential RBP for NMD (nonsense-mediated decay), directly interacts with Ago2 and enhances target repression through CCR4-NOT-mediated deadenylation, which is another mechanism of miRNA targeting regulation by RBPs (Park et al., 2019). Therefore, utilizing the precise binding site location for each RBP species would allow for accurate miRNA targeting prediction.

Target RNA-directed miRNA degradation

In general, stronger binding between a miRNA and its target leads to better repression. However, when a perfect or near-perfect match occurs in both the seed region and the 3′ regions of the miRNA, the miRNA may be subject to degradation, a phenomenon known as target RNA-directed miRNA degradation (TDMD) (Ameres et al., 2010). TDMD was initially studied using artificial exogenous RNAs (Ameres et al., 2010) and viral RNAs (Cazalla et al., 2010). Some viral RNA transcripts, including those of HVS (herpesvirus saimiri)(Cazalla et al., 2010), MCMV (murine cytomegalovirus) (Libri et al., 2012), and HCMV (human cytomegalovirus) (Lee et al., 2013), utilize TDMD as an evolutionary strategy for immune evasion. However, the physiological role of TDMD in animal cells has been debated (Fuchs Wightman et al., 2018).

Recently, the physiological role of TDMD induced by endogenous miRNA targets has been suggested in three important RNA examples: Cyrano, Nrep, and Serpine1. Cyrano is a long non-coding RNA (lncRNA) broadly conserved in vertebrates and plays a crucial role in early embryonic development (Ulitsky et al., 2011). Cyrano pairs perfectly with miR-7 both at the seed and 3′ regions, inducing effective degradation of miR-7 (Kleaveland et al., 2018). Degradation of miR-7 by Cyrano leads to the derepression and accumulation of the circular RNA (circRNA) Cdr1as, which is associated with neuronal development, neurodegenerative diseases, and cancer metastasis (Memczak et al., 2013). Further studies found that TDMD by Cyrano is mediated by ZSWIM8 polyubiquitin ligase, which induces degradation of AGO by polyubiquitination and exposes the miRNA to destruction (Han et al., 2020a; Shi et al., 2020). Serpine1 plays an important role in the cell cycle re-entry of quiescent cells (Iyer et al., 1999). The serpine1 transcript induces TDMD of miR-30b/c during cell cycle re-entry through miRNA tailing, and the degradation of miR-30b/c is crucial for accelerating the G1/S transition during cell cycle re-entry (Ghini et al., 2018). The Nrep transcript (known as libra in zebrafish) is expressed in vertebrate brains and induces TDMD of miR-29b by 3′ trimming but not tailing. Scrambling the Nrep transcript to inhibit TDMD results in an increase in miR-29b expression in the cerebellum, which in turn leads to behavioral deficits in mice (Bitetti et al., 2018). These examples demonstrate the significant physiological effects of TDMD on endogenous transcripts in mammals.

Competitive endogenous RNAs

It is widely recognized that the abundance of miRNA target sites is negatively correlated with targeting efficacy (Arvey et al., 2010; Garcia et al., 2011). This correlation can be explained by the competitive endogenous RNA (ceRNA) hypothesis that transcriptome-wide competition occurs among target sites, including those in the mRNAs, circRNAs, pseudogenes, and lncRNAs (Salmena et al., 2011; Thomson and Dinger, 2016). Experimental evidence indicates that some transcripts, called miRNA sponges, can effectively sequester miRNAs. For instance, the PTEN pseudogene (PTENP1) (Poliseno et al., 2010) and Braf pseudogene (BRAFP1) (Karreth et al., 2015) act as ceRNAs and increase the expression of PTEN and BRAF through derepression. lncRNA Linc-md1 acts as a ceRNA and regulates muscle differentiation (Cesana et al., 2011). circRNA Cdr1as contains more than 70 miR-7 target sites and can effectively downregulate the repressive efficacy of miR-7 (Hansen et al., 2013; Memczak et al., 2013). However, some studies have suggested that excessively high levels of ceRNA expression are required for their function (Bosson et al., 2014; Denzler et al., 2014). Despite this, ceRNAs can be effective when the concentration of miRNA is sufficiently low by TDMD (Jens and Rajewsky, 2015).

DEEP LEARNING MAY IMPROVE miRNA TARGET PREDICTION

To effectively incorporate previously unappreciated determinants into computational miRNA target prediction tools, high-dimensional input data should be utilized (Fig. 4). For instance, RBP-binding site information generated by eCLIP-seq (Van Nostrand et al., 2016) is represented in a multidimensional score matrix. To account for the RMs in the miRNA and 3′UTR, a large-scale dataset such as m6A-seq (Dominissini et al., 2012) should be utilized. An alternative approach that maps RMs directly from Oxford Nanopore sequencing has been suggested and can be used for miRNA target prediction (Leger et al., 2021). Secondary structures of the 3′UTRs can also be incorporated into the model using selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE)-seq (Lucks et al., 2011) and dimethyl sulfate (DMS)-seq (Rouskin et al., 2014). By combining all these datasets, a high-dimensional input dataset can be prepared for functional miRNA target prediction use.

Figure F4
Examples of high-dimensional input data (left), including eCLIP-seq (Van Nostrand et al., 2016), SHAPE-seq (Lucks et al., 2011), Nanopore (Jain et al., 2016), DMS-seq (Rouskin et al., 2014), and PhyloP ...

Deep learning approaches can be used to integrate high-dimensional datasets for accurate miRNA target prediction. Recently, deep learning models, such as CNNs (convolutional neural networks) (Townshend et al., 2021; Zeng et al., 2016), RNNs (recurrent neural networks) (Quang and Xie, 2016; Tasdelen and Sen, 2021), and transformers (Ji et al., 2021; Wang et al., 2021) that utilize nucleic acid sequences and multi-omics data, have generated promising results. Indeed, multiple efforts have been made to employ deep learning for miRNA target prediction (McGeary et al., 2019; Min et al., 2022; Talukder et al., 2022). However, their overall predictive accuracies have been limited, perhaps because these models do not fully integrate the aforementioned determinants of miRNA targeting into the deep learning model. Therefore, a deep learning model that integrates high-dimensional data, accounting for a more comprehensive list of determinants, including those that have not been used previously, may yield significantly improved predictions of functional miRNA targeting.

CONCLUSION

miRNAs are indispensable regulators of gene expression that affect various biological processes, pathways, and diseases by forming complex regulatory networks. Currently, a large number of miRNA target prediction tools have been developed and are widely used. Most of these tools use seed pairing with other well-known determinants of miRNA targeting, including local AU content, structural accessibility, seed-pairing stability, target abundance, and evolutionary conservation. Each of these determinants is significantly correlated with miRNA targeting efficacy and helps to accurately predict functional miRNA targeting. Tools that incorporate a larger number of determinants, such as TargetScan, produce more accurate predictions than those that do not. However, all currently available computational tools exhibit limited prediction accuracy.

Incorporation of unappreciated miRNA targeting determinants such as RMs in mature miRNAs, RMs in miRNA target mRNAs, target RNA-directed miRNA degradation, ceRNAs, and RBP-binding sites may further improve the accuracy of the computational tools aforementioned. For most of these determinants, their regulatory impacts on miRNA targeting have been validated in various studies, yet some determinants, such as RMs, should be assessed more systematically. To integrate these unappreciated determinants into the prediction model, high-dimensional multiomics datasets should be utilized, and improving computational tools for functional miRNA prediction would help to accurately decipher the complex regulatory network of miRNAs.

Article information

Mol. Cells.Jan 31, 2023; 46(1): 21-32.
Published online 2023-01-4. doi:  10.14348/molcells.2023.2157
School of Biological Sciences, Seoul National University, Seoul 08826, Korea
*Correspondence: baek@snu.ac.kr
Received October 14, 2022; Accepted November 15, 2022.
Articles from Mol. Cells are provided here courtesy of Mol. Cells

References

  • Agarwal, V., Bell, G.W., Nam, J.W., Bartel, D.P. (2015). Predicting effective microRNA target sites in mammalian mRNAs. Elife. 4, e05005.
  • Ahuja, D., Goyal, A., Ray, P.S. (2016). Interplay between RNA-binding protein HuR and microRNA-125b regulates p53 mRNA translation in response to genotoxic stress. RNA Biol.. 13, 1152-1165.
  • Alarcon, C.R., Goodarzi, H., Lee, H., Liu, X., Tavazoie, S., Tavazoie, S.F. (2015a). HNRNPA2B1 is a mediator of m6A-dependent nuclear RNA processing events. Cell. 162, 1299-1308.
  • Alarcon, C.R., Lee, H., Goodarzi, H., Halberg, N., Tavazoie, S.F. (2015b). N6-methyladenosine marks primary microRNAs for processing. Nature. 519, 482-485.
  • Ameres, S.L., Horwich, M.D., Hung, J.H., Xu, J., Ghildiyal, M., Weng, Z., Zamore, P.D. (2010). Target RNA-directed trimming and tailing of small silencing RNAs. Science. 328, 1534-1539.
  • Arvey, A., Larsson, E., Sander, C., Leslie, C.S., Marks, D.S. (2010). Target mRNA abundance dilutes microRNA and siRNA activity. Mol. Syst. Biol.. 6, 363.
  • Baek, D., Villen, J., Shin, C., Camargo, F.D., Gygi, S.P., Bartel, D.P. (2008). The impact of microRNAs on protein output. Nature. 455, 64-71.
  • Bartel, D.P. (2009). MicroRNAs: target recognition and regulatory functions. Cell. 136, 215-233.
  • Bass, B.L., Weintraub, H. (1988). An unwinding activity that covalently modifies its double-stranded-RNA substrate. Cell. 55, 1089-1098.
  • Bernhart, S.H., Hofacker, I.L., Stadler, P.F. (2006). Local RNA base pairing probabilities in large sequences. Bioinformatics. 22, 614-615.
  • Bernstein, E., Kim, S.Y., Carmell, M.A., Murchison, E.P., Alcorn, H., Li, M.Z., Mills, A.A., Elledge, S.J., Anderson, K.V., Hannon, G.J. (2003). Dicer is essential for mouse development. Nat. Genet.. 35, 215-217.
  • Betel, D., Koppal, A., Agius, P., Sander, C., Leslie, C. (2010). Comprehensive modeling of microRNA targets predicts functional non-conserved and non-canonical sites. Genome Biol.. 11, R90.
  • Bitetti, A., Mallory, A.C., Golini, E., Carrieri, C., Carreno Gutierrez, H., Perlas, E., Perez-Rico, Y.A., Tocchini-Valentini, G.P., Enright, A.J., Norton, W.H.J. (2018). MicroRNA degradation by a conserved target RNA regulates animal behavior. Nat. Struct. Mol. Biol.. 25, 244-251.
  • Blow, M.J., Grocock, R.J., van Dongen, S., Enright, A.J., Dicks, E., Futreal, P.A., Wooster, R., Stratton, M.R. (2006). RNA editing of human microRNAs. Genome Biol.. 7, R27.
  • Bohnsack, K.E., Höbartner, C., Bohnsack, M.T. (2019). Eukaryotic 5-methylcytosine (m⁵C) RNA methyltransferases: mechanisms, cellular functions, and links to disease. Genes (Basel). 10, 102.
  • Bosson, A.D., Zamudio, J.R., Sharp, P.A. (2014). Endogenous miRNA and target concentrations determine susceptibility to potential ceRNA competition. Mol. Cell. 56, 347-359.
  • Bracken, C.P., Scott, H.S., Goodall, G.J. (2016). A network-biology perspective of microRNA function and dysfunction in cancer. Nat. Rev. Genet.. 17, 719-732.
  • Brennan, G.P., Henshall, D.C. (2020). MicroRNAs as regulators of brain function and targets for treatment of epilepsy. Nat. Rev. Neurol.. 16, 506-519.
  • Briskin, D., Wang, P.Y., Bartel, D.P. (2020). The biochemical basis for the cooperative action of microRNAs. Proc. Natl. Acad. Sci. U. S. A.. 117, 17764-17774.
  • Care, A., Catalucci, D., Felicetti, F., Bonci, D., Addario, A., Gallo, P., Bang, M.L., Segnalini, P., Gu, Y., Dalton, N.D. (2007). MicroRNA-133 controls cardiac hypertrophy. Nat. Med.. 13, 613-618.
  • Cazalla, D., Yario, T., Steitz, J.A. (2010). Down-regulation of a host microRNA by a Herpesvirus saimiri noncoding RNA. Science. 328, 1563-1566.
  • Cesana, M., Cacchiarelli, D., Legnini, I., Santini, T., Sthandier, O., Chinappi, M., Tramontano, A., Bozzoni, I. (2011). A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell. 147, 358-369.
  • Cesarini, V., Silvestris, D.A., Tassinari, V., Tomaselli, S., Alon, S., Eisenberg, E., Locatelli, F., Gallo, A. (2018). ADAR2/miR-589-3p axis controls glioblastoma cell migration/invasion. Nucleic Acids Res.. 46, 2045-2059.
  • Cheray, M., Etcheverry, A., Jacques, C., Pacaud, R., Bougras-Cartron, G., Aubry, M., Denoual, F., Peterlongo, P., Nadaradjane, A., Briand, J. (2020). Cytosine methylation of mature microRNAs inhibits their functions and is associated with poor prognosis in glioblastoma multiforme. Mol. Cancer. 19, 36.
  • Chi, S.W., Zang, J.B., Mele, A., Darnell, R.B. (2009). Argonaute HITS-CLIP decodes microRNA-mRNA interaction maps. Nature. 460, 479-486.
  • Choi, W.Y., Giraldez, A.J., Schier, A.F. (2007). Target protectors reveal dampening and balancing of Nodal agonist and antagonist by miR-430. Science. 318, 271-274.
  • Denzler, R., Agarwal, V., Stefano, J., Bartel, D.P., Stoffel, M. (2014). Assessing the ceRNA hypothesis with quantitative measurements of miRNA and target abundance. Mol. Cell. 54, 766-776.
  • Dominissini, D., Moshitch-Moshkovitz, S., Schwartz, S., Salmon-Divon, M., Ungar, L., Osenberg, S., Cesarkas, K., Jacob-Hirsch, J., Amariglio, N., Kupiec, M. (2012). Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature. 485, 201-206.
  • Fang, Z., Rajewsky, N. (2011). The impact of miRNA target sites in coding sequences and in 3'UTRs. PLoS One. 6, e18067.
  • Fededa, J.P., Esk, C., Mierzwa, B., Stanyte, R., Yuan, S., Zheng, H., Ebnet, K., Yan, W., Knoblich, J.A., Gerlich, D.W. (2016). MicroRNA-34/449 controls mitotic spindle orientation during mammalian cortex development. EMBO J.. 35, 2386-2398.
  • Fontana, L., Pelosi, E., Greco, P., Racanicchi, S., Testa, U., Liuzzi, F., Croce, C.M., Brunetti, E., Grignani, F., Peschle, C. (2007). MicroRNAs 17-5p-20a-106a control monocytopoiesis through AML1 targeting and M-CSF receptor upregulation. Nat. Cell Biol.. 9, 775-787.
  • Friedman, R.C., Farh, K.K., Burge, C.B., Bartel, D.P. (2009). Most mammalian mRNAs are conserved targets of microRNAs. Genome Res.. 19, 92-105.
  • Fuchs Wightman, F., Giono, L.E., Fededa, J.P., de la Mata, M. (2018). Target RNAs strike back on microRNAs. Front. Genet.. 9, 435.
  • Gaidatzis, D., van Nimwegen, E., Hausser, J., Zavolan, M. (2007). Inference of miRNA targets using evolutionary conservation and pathway analysis. BMC Bioinformatics. 8, 69.
  • Garcia, D.M., Baek, D., Shin, C., Bell, G.W., Grimson, A., Bartel, D.P. (2011). Weak seed-pairing stability and high target-site abundance decrease the proficiency of lsy-6 and other microRNAs. Nat. Struct. Mol. Biol.. 18, 1139-1146.
  • Ghanbarian, H., Yildiz, M.T., Tutar, Y. (2022). MicroRNA targeting. Methods Mol. Biol.. 2257, 105-130.
  • Ghini, F., Rubolino, C., Climent, M., Simeone, I., Marzi, M.J., Nicassio, F. (2018). Endogenous transcripts control miRNA levels and activity in mammalian cells by target-directed miRNA degradation. Nat. Commun.. 9, 3119.
  • Ghosh, T., Aprea, J., Nardelli, J., Engel, H., Selinger, C., Mombereau, C., Lemonnier, T., Moutkine, I., Schwendimann, L., Dori, M. (2014). MicroRNAs establish robustness and adaptability of a critical gene network to regulate progenitor fate decisions during cortical neurogenesis. Cell Rep.. 7, 1779-1788.
  • Giraldez, A.J., Cinalli, R.M., Glasner, M.E., Enright, A.J., Thomson, J.M., Baskerville, S., Hammond, S.M., Bartel, D.P., Schier, A.F. (2005). MicroRNAs regulate brain morphogenesis in zebrafish. Science. 308, 833-838.
  • Gregory, P.A., Bert, A.G., Paterson, E.L., Barry, S.C., Tsykin, A., Farshid, G., Vadas, M.A., Khew-Goodall, Y., Goodall, G.J. (2008). The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nat. Cell Biol.. 10, 593-601.
  • Gregory, R.I., Chendrimada, T.P., Cooch, N., Shiekhattar, R. (2005). Human RISC couples microRNA biogenesis and posttranscriptional gene silencing. Cell. 123, 631-640.
  • Grimson, A., Farh, K.K., Johnston, W.K., Garrett-Engele, P., Lim, L.P., Bartel, D.P. (2007). MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol. Cell. 27, 91-105.
  • Gumienny, R., Zavolan, M. (2015). Accurate transcriptome-wide prediction of microRNA targets and small interfering RNA off-targets with MIRZA-G. Nucleic Acids Res.. 43, 1380-1391.
  • Hammond, S.M., Bernstein, E., Beach, D., Hannon, G.J. (2000). An RNA-directed nuclease mediates post-transcriptional gene silencing in Drosophila cells. Nature. 404, 293-296.
  • Han, J., LaVigne, C.A., Jones, B.T., Zhang, H., Gillett, F., Mendell, J.T. (2020a). A ubiquitin ligase mediates target-directed microRNA decay independently of tailing and trimming. Science. 370, eabc9546.
  • Han, X., Luan, T., Sun, Y., Yan, W., Wang, D., Zeng, X. (2020b). MicroRNA 449c mediates the generation of monocytic myeloid-derived suppressor cells by targeting STAT6. Mol. Cells. 43, 793-803.
  • Hansen, T.B., Jensen, T.I., Clausen, B.H., Bramsen, J.B., Finsen, B., Damgaard, C.K., Kjems, J. (2013). Natural RNA circles function as efficient microRNA sponges. Nature. 495, 384-388.
  • Harris, T.A., Yamakuchi, M., Ferlito, M., Mendell, J.T., Lowenstein, C.J. (2008). MicroRNA-126 regulates endothelial expression of vascular cell adhesion molecule 1. Proc. Natl. Acad. Sci. U. S. A.. 105, 1516-1521.
  • Hassan, T., Smith, S.G., Gaughan, K., Oglesby, I.K., O'Neill, S., McElvaney, N.G., Greene, C.M. (2013). Isolation and identification of cell-specific microRNAs targeting a messenger RNA using a biotinylated anti-sense oligonucleotide capture affinity technique. Nucleic Acids Res.. 41, e71.
  • Hausser, J., Landthaler, M., Jaskiewicz, L., Gaidatzis, D., Zavolan, M. (2009). Relative contribution of sequence and structure features to the mRNA binding of Argonaute/EIF2C-miRNA complexes and the degradation of miRNA targets. Genome Res.. 19, 2009-2020.
  • Hong, X., Hammell, M., Ambros, V., Cohen, S.M. (2009). Immunopurification of Ago1 miRNPs selects for a distinct class of microRNA targets. Proc. Natl. Acad. Sci. U. S. A.. 106, 15085-15090.
  • Iyer, V.R., Eisen, M.B., Ross, D.T., Schuler, G., Moore, T., Lee, J.C., Trent, J.M., Staudt, L.M., Hudson, J., Boguski, M.S. (1999). The transcriptional program in the response of human fibroblasts to serum. Science. 283, 83-87.
  • Jain, M., Olsen, H.E., Paten, B., Akeson, M. (2016). The Oxford Nanopore MinION: delivery of nanopore sequencing to the genomics community. Genome Biol.. 17, 239.
  • Jan, C.H., Friedman, R.C., Ruby, J.G., Bartel, D.P. (2011). Formation, regulation and evolution of Caenorhabditis elegans 3'UTRs. Nature. 469, 97-101.
  • Jang, H., Park, S., Kim, J., Kim, J.H., Kim, S.Y., Cho, S., Park, S.G., Park, B.C., Kim, S., Kim, J.H. (2020). The tumor suppressor, p53, negatively regulates non-canonical NF-κB signaling through miRNA-induced silencing of NF-κB-inducing kinase. Mol. Cells. 43, 23-33.
  • Jens, M., Rajewsky, N. (2015). Competition between target sites of regulators shapes post-transcriptional gene regulation. Nat. Rev. Genet.. 16, 113-126.
  • Ji, Y., Zhou, Z., Liu, H., Davuluri, R.V. (2021). DNABERT: pre-trained Bidirectional Encoder Representations from Transformers model for DNA-language in genome. Bioinformatics. 37, 2112-2120.
  • Jing, Q., Huang, S., Guth, S., Zarubin, T., Motoyama, A., Chen, J., Di Padova, F., Lin, S.C., Gram, H., Han, J. (2005). Involvement of microRNA in AU-rich element-mediated mRNA instability. Cell. 120, 623-634.
  • Karreth, F.A., Reschke, M., Ruocco, A., Ng, C., Chapuy, B., Leopold, V., Sjoberg, M., Keane, T.M., Verma, A., Ala, U. (2015). The BRAF pseudogene functions as a competitive endogenous RNA and induces lymphoma in vivo. Cell. 161, 319-332.
  • Kawahara, Y., Zinshteyn, B., Sethupathy, P., Iizasa, H., Hatzigeorgiou, A.G., Nishikura, K. (2007). Redirection of silencing targets by adenosine-to-inosine editing of miRNAs. Science. 315, 1137-1140.
  • Kedde, M., van Kouwenhove, M., Zwart, W., Oude Vrielink, J.A., Elkon, R., Agami, R. (2010). A Pumilio-induced RNA structure switch in p27-3' UTR controls miR-221 and miR-222 accessibility. Nat. Cell Biol.. 12, 1014-1020.
  • Keene, J.D., Komisarow, J.M., Friedersdorf, M.B. (2006). RIP-Chip: the isolation and identification of mRNAs, microRNAs and protein components of ribonucleoprotein complexes from cell extracts. Nat. Protoc.. 1, 302-307.
  • Kertesz, M., Iovino, N., Unnerstall, U., Gaul, U., Segal, E. (2007). The role of site accessibility in microRNA target recognition. Nat. Genet.. 39, 1278-1284.
  • Kim, D., Kim, J., Baek, D. (2014). Global and local competition between exogenously introduced microRNAs and endogenously expressed microRNAs. Mol. Cells. 37, 412-417.
  • Kim, D., Sung, Y.M., Park, J., Kim, S., Kim, J., Park, J., Ha, H., Bae, J.Y., Kim, S., Baek, D. (2016). General rules for functional microRNA targeting. Nat. Genet.. 48, 1517-1526.
  • Kim, S., Kim, S., Chang, H.R., Kim, D., Park, J., Son, N., Park, J., Yoon, M., Chae, G., Kim, Y.K. (2021). The regulatory impact of RNA-binding proteins on microRNA targeting. Nat. Commun.. 12, 5057.
  • Kleaveland, B., Shi, C.Y., Stefano, J., Bartel, D.P. (2018). A network of noncoding regulatory RNAs acts in the mammalian brain. Cell. 174, 350-362.e17.
  • Konno, M., Koseki, J., Asai, A., Yamagata, A., Shimamura, T., Motooka, D., Okuzaki, D., Kawamoto, K., Mizushima, T., Eguchi, H. (2019). Distinct methylation levels of mature microRNAs in gastrointestinal cancers. Nat. Commun.. 10, 3888.
  • Krek, A., Grun, D., Poy, M.N., Wolf, R., Rosenberg, L., Epstein, E.J., MacMenamin, P., da Piedade, I., Gunsalus, K.C., Stoffel, M. (2005). Combinatorial microRNA target predictions. Nat. Genet.. 37, 495-500.
  • Krichevsky, A.M., Sonntag, K.C., Isacson, O., Kosik, K.S. (2006). Specific microRNAs modulate embryonic stem cell-derived neurogenesis. Stem Cells. 24, 857-864.
  • Kruger, J., Rehmsmeier, M. (2006). RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res.. 34, .
  • Kume, H., Hino, K., Galipon, J., Ui-Tei, K. (2014). A-to-I editing in the miRNA seed region regulates target mRNA selection and silencing efficiency. Nucleic Acids Res.. 42, 10050-10060.
  • Lambert, N., Robertson, A., Jangi, M., McGeary, S., Sharp, P.A., Burge, C.B. (2014). RNA Bind-n-Seq: quantitative assessment of the sequence and structural binding specificity of RNA binding proteins. Mol. Cell. 54, 887-900.
  • Lee, S., Song, J., Kim, S., Kim, J., Hong, Y., Kim, Y., Kim, D., Baek, D., Ahn, K. (2013). Selective degradation of host MicroRNAs by an intergenic HCMV noncoding RNA accelerates virus production. Cell Host Microbe. 13, 678-690.
  • Legendre, M., Ritchie, W., Lopez, F., Gautheret, D. (2006). Differential repression of alternative transcripts: a screen for miRNA targets. PLoS Comput. Biol.. 2, e43.
  • Leger, A., Amaral, P.P., Pandolfini, L., Capitanchik, C., Capraro, F., Miano, V., Migliori, V., Toolan-Kerr, P., Sideri, T., Enright, A.J. (2021). RNA modifications detection by comparative Nanopore direct RNA sequencing. Nat. Commun.. 12, 7198.
  • Lewis, B.P., Burge, C.B., Bartel, D.P. (2005). Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 120, 15-20.
  • Li, Q., Song, X.W., Zou, J., Wang, G.K., Kremneva, E., Li, X.Q., Zhu, N., Lappalainen, P., Yuan, W.J., Qin, Y.W. (2010). Attenuation of microRNA-1 derepresses the cytoskeleton regulatory protein twinfilin-1 to provoke cardiac hypertrophy (vol123, pg2444, 2010). J. Cell Sci.. 123, 2680.
  • Liao, S., Sun, H., Xu, C. (2018). YTH domain: a family of N6-methyladenosine (m6A) readers. Genomics Proteomics Bioinformatics. 16, 99-107.
  • Libri, V., Helwak, A., Miesen, P., Santhakumar, D., Borger, J.G., Kudla, G., Grey, F., Tollervey, D., Buck, A.H. (2012). Murine cytomegalovirus encodes a miR-27 inhibitor disguised as a target. Proc. Natl. Acad. Sci. U. S. A.. 109, 279-284.
  • Lim, L.P., Lau, N.C., Garrett-Engele, P., Grimson, A., Schelter, J.M., Castle, J., Bartel, D.P., Linsley, P.S., Johnson, J.M. (2005). Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 433, 769-773.
  • Lin, X., Yang, B., Liu, W., Tan, X., Wu, F., Hu, P., Jiang, T., Bao, Z., Yuan, J., Qiang, B. (2016). Interplay between PCBP2 and miRNA modulates ARHGDIA expression and function in glioma migration and invasion. Oncotarget. 7, 19483-19498.
  • Liu, J., Yue, Y., Han, D., Wang, X., Fu, Y., Zhang, L., Jia, G., Yu, M., Lu, Z., Deng, X. (2014). A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation. Nat. Chem. Biol.. 10, 93-95.
  • Loughrey, D., Watters, K.E., Settle, A.H., Lucks, J.B. (2014). SHAPE-Seq 2.0: systematic optimization and extension of high-throughput chemical probing of RNA secondary structure with next generation sequencing. Nucleic Acids Res.. 42, e165.
  • Luciano, D.J., Mirsky, H., Vendetti, N.J., Maas, S. (2004). RNA editing of a miRNA precursor. RNA. 10, 1174-1177.
  • Lucks, J.B., Mortimer, S.A., Trapnell, C., Luo, S., Aviran, S., Schroth, G.P., Pachter, L., Doudna, J.A., Arkin, A.P. (2011). Multiplexed RNA structure characterization with selective 2'-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq). Proc. Natl. Acad. Sci. U. S. A.. 108, 11063-11068.
  • Majoros, W.H., Ohler, U. (2007). Spatial preferences of microRNA targets in 3' untranslated regions. BMC Genomics. 8, 152.
  • McGeary, S.E., Bisaria, N., Pham, T.M., Wang, P.Y., Bartel, D.P. (2022). MicroRNA 3'-compensatory pairing occurs through two binding modes, with affinity shaped by nucleotide identity and position. Elife. 11, e69803.
  • McGeary, S.E., Lin, K.S., Shi, C.Y., Pham, T.M., Bisaria, N., Kelley, G.M., Bartel, D.P. (2019). The biochemical basis of microRNA targeting efficacy. Science. 366, eaav1741.
  • Mehta, A., Baltimore, D. (2016). MicroRNAs as regulatory elements in immune system logic. Nat. Rev. Immunol.. 16, 279-294.
  • Memczak, S., Jens, M., Elefsinioti, A., Torti, F., Krueger, J., Rybak, A., Maier, L., Mackowiak, S.D., Gregersen, L.H., Munschauer, M. (2013). Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 495, 333-338.
  • Michaels, M.L., Cruz, C., Grollman, A.P., Miller, J.H. (1992). Evidence that MutY and MutM combine to prevent mutations by an oxidatively damaged form of guanine in DNA. Proc. Natl. Acad. Sci. U. S. A.. 89, 7022-7025.
  • Min, S., Lee, B., Yoon, S. (2022). TargetNet: functional microRNA target prediction with deep neural networks. Bioinformatics. 38, 671-677.
  • Muljo, S.A., Ansel, K.M., Kanellopoulou, C., Livingston, D.M., Rao, A., Rajewsky, K. (2005). Aberrant T cell differentiation in the absence of Dicer. J. Exp. Med.. 202, 261-269.
  • Muller, S., Glass, M., Singh, A.K., Haase, J., Bley, N., Fuchs, T., Lederer, M., Dahl, A., Huang, H., Chen, J. (2019). IGF2BP1 promotes SRF-dependent transcription in cancer in a m6A- and miRNA-dependent manner. Nucleic Acids Res.. 47, 375-390.
  • Nachtigall, P.G., Bovolenta, L.A. (2022). Computational detection of microRNA targets. Methods Mol. Biol.. 2257, 187-209.
  • Nam, J.W., Rissland, O.S., Koppstein, D., Abreu-Goodger, C., Jan, C.H., Agarwal, V., Yildirim, M.A., Rodriguez, A., Bartel, D.P. (2014). Global analyses of the effect of different cellular contexts on microRNA targeting. Mol. Cell. 53, 1031-1043.
  • Nielsen, C.B., Shomron, N., Sandberg, R., Hornstein, E., Kitzman, J., Burge, C.B. (2007). Determinants of targeting by endogenous and exogenous microRNAs and siRNAs. RNA. 13, 1894-1910.
  • Nishikura, K. (2016). A-to-I editing of coding and non-coding RNAs by ADARs. Nat. Rev. Mol. Cell Biol.. 17, 83-96.
  • O'Connell, R.M., Rao, D.S., Chaudhuri, A.A., Baltimore, D. (2010). Physiological and pathological roles for microRNAs in the immune system. Nat. Rev. Immunol.. 10, 111-122.
  • O'Donnell, K.A., Wentzel, E.A., Zeller, K.I., Dang, C.V., Mendell, J.T. (2005). c-Myc-regulated microRNAs modulate E2F1 expression. Nature. 435, 839-843.
  • Orom, U.A., Lund, A.H. (2010). Experimental identification of microRNA targets. Gene. 451, 1-5.
  • Park, J., Seo, J.W., Ahn, N., Park, S., Hwang, J., Nam, J.W. (2019). UPF1/SMG7-dependent microRNA-mediated gene regulation. Nat. Commun.. 10, 4181.
  • Peng, Y., Croce, C.M. (2016). The role of MicroRNAs in human cancer. Signal Transduct. Target. Ther.. 1, 15004.
  • Pettersen, E.F., Goddard, T.D., Huang, C.C., Meng, E.C., Couch, G.S., Croll, T.I., Morris, J.H., Ferrin, T.E. (2021). UCSF ChimeraX: structure visualization for researchers, educators, and developers. Protein Sci.. 30, 70-82.
  • Plaisier, C.L., Pan, M., Baliga, N.S. (2012). A miRNA-regulatory network explains how dysregulated miRNAs perturb oncogenic processes across diverse cancers. Genome Res.. 22, 2302-2314.
  • Poliseno, L., Salmena, L., Zhang, J., Carver, B., Haveman, W.J., Pandolfi, P.P. (2010). A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 465, 1033-1038.
  • Pollard, K.S., Hubisz, M.J., Rosenbloom, K.R., Siepel, A. (2010). Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res.. 20, 110-121.
  • Pu, M., Chen, J., Tao, Z., Miao, L., Qi, X., Wang, Y., Ren, J. (2019). Regulatory network of miRNA on its target: coordination between transcriptional and post-transcriptional regulation of gene expression. Cell. Mol. Life Sci.. 76, 441-451.
  • Qian, B., Wang, P., Zhang, D., Wu, L. (2021). m6A modification promotes miR-133a repression during cardiac development and hypertrophy via IGF2BP2. Cell Death Discov.. 7, 157.
  • Quang, D., Xie, X. (2016). DanQ: a hybrid convolutional and recurrent deep neural network for quantifying the function of DNA sequences. Nucleic Acids Res.. 44, e107.
  • Reczko, M., Maragkakis, M., Alexiou, P., Grosse, I., Hatzigeorgiou, A.G. (2012). Functional microRNA targets in protein coding sequences. Bioinformatics. 28, 771-776.
  • Rodriguez, A., Vigorito, E., Clare, S., Warren, M.V., Couttet, P., Soond, D.R., van Dongen, S., Grocock, R.J., Das, P.P., Miska, E.A. (2007). Requirement of bic/microRNA-155 for normal immune function. Science. 316, 608-611.
  • Roundtree, I.A., Evans, M.E., Pan, T., He, C. (2017). Dynamic RNA modifications in gene expression regulation. Cell. 169, 1187-1200.
  • Rouskin, S., Zubradt, M., Washietl, S., Kellis, M., Weissman, J.S. (2014). Genome-wide probing of RNA structure reveals active unfolding of mRNA structures in vivo. Nature. 505, 701-705.
  • Saetrom, P., Heale, B.S., Snove, O., Aagaard, L., Alluin, J., Rossi, J.J. (2007). Distance constraints between microRNA target sites dictate efficacy and cooperativity. Nucleic Acids Res.. 35, 2333-2342.
  • Salmena, L., Poliseno, L., Tay, Y., Kats, L., Pandolfi, P.P. (2011). A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language?. Cell. 146, 353-358.
  • Sandberg, R., Neilson, J.R., Sarma, A., Sharp, P.A., Burge, C.B. (2008). Proliferating cells express mRNAs with shortened 3' untranslated regions and fewer microRNA target sites. Science. 320, 1643-1647.
  • Schirle, N.T., Sheu-Gruttadauria, J., MacRae, I.J. (2014). Structural basis for microRNA targeting. Science. 346, 608-613.
  • Selbach, M., Schwanhausser, B., Thierfelder, N., Fang, Z., Khanin, R., Rajewsky, N. (2008). Widespread changes in protein synthesis induced by microRNAs. Nature. 455, 58-63.
  • Seok, H., Lee, H., Lee, S., Ahn, S.H., Lee, H.S., Kim, G.D., Peak, J., Park, J., Cho, Y.K., Jeong, Y. (2020). Position-specific oxidation of miR-1 encodes cardiac hypertrophy. Nature. 584, 279-285.
  • Shi, C.Y., Kingston, E.R., Kleaveland, B., Lin, D.H., Stubna, M.W., Bartel, D.P. (2020). The ZSWIM8 ubiquitin ligase mediates target-directed microRNA degradation. Science. 370, eabc9359.
  • Squires, J.E., Patel, H.R., Nousch, M., Sibbritt, T., Humphreys, D.T., Parker, B.J., Suter, C.M., Preiss, T. (2012). Widespread occurrence of 5-methylcytosine in human coding and non-coding RNA. Nucleic Acids Res.. 40, 5023-5033.
  • Stark, A., Brennecke, J., Bushati, N., Russell, R.B., Cohen, S.M. (2005). Animal MicroRNAs confer robustness to gene expression and have a significant impact on 3'UTR evolution. Cell. 123, 1133-1146.
  • Talukder, A., Zhang, W., Li, X., Hu, H. (2022). A deep learning method for miRNA/isomiR target detection. Sci. Rep.. 12, 10618.
  • Tan, C.L., Plotkin, J.L., Veno, M.T., von Schimmelmann, M., Feinberg, P., Mann, S., Handler, A., Kjems, J., Surmeier, D.J., O'Carroll, D. (2013). MicroRNA-128 governs neuronal excitability and motor behavior in mice. Science. 342, 1254-1258.
  • Tasdelen, A., Sen, B. (2021). A hybrid CNN-LSTM model for pre-miRNA classification. Sci. Rep.. 11, 14125.
  • Thadani, R., Tammi, M.T. (2006). MicroTar: predicting microRNA targets from RNA duplexes. BMC Bioinformatics. 7, S20.
  • Thai, T.H., Calado, D.P., Casola, S., Ansel, K.M., Xiao, C., Xue, Y., Murphy, A., Frendewey, D., Valenzuela, D., Kutok, J.L. (2007). Regulation of the germinal center response by microRNA-155. Science. 316, 604-608.
  • Thomson, D.W., Bracken, C.P., Goodall, G.J. (2011). Experimental strategies for microRNA target identification. Nucleic Acids Res.. 39, 6845-6853.
  • Thomson, D.W., Dinger, M.E. (2016). Endogenous microRNA sponges: evidence and controversy. Nat. Rev. Genet.. 17, 272-283.
  • Tian, B., Hu, J., Zhang, H., Lutz, C.S. (2005). A large-scale analysis of mRNA polyadenylation of human and mouse genes. Nucleic Acids Res.. 33, 201-212.
  • Townshend, R.J.L., Eismann, S., Watkins, A.M., Rangan, R., Karelina, M., Das, R., Dror, R.O. (2021). Geometric deep learning of RNA structure. Science. 373, 1047-1051.
  • Tuschl, T., Zamore, P.D., Lehmann, R., Bartel, D.P., Sharp, P.A. (1999). Targeted mRNA degradation by double-stranded RNA in vitro. Genes Dev.. 13, 3191-3197.
  • Ulitsky, I., Shkumatava, A., Jan, C.H., Sive, H., Bartel, D.P. (2011). Conserved function of lincRNAs in vertebrate embryonic development despite rapid sequence evolution. Cell. 147, 1537-1550.
  • Van Nostrand, E.L., Pratt, G.A., Shishkin, A.A., Gelboin-Burkhart, C., Fang, M.Y., Sundararaman, B., Blue, S.M., Nguyen, T.B., Surka, C., Elkins, K. (2016). Robust transcriptome-wide discovery of RNA-binding protein binding sites with enhanced CLIP (eCLIP). Nat. Methods. 13, 508-514.
  • Van Nostrand, E.L., Pratt, G.A., Yee, B.A., Wheeler, E.C., Blue, S.M., Mueller, J., Park, S.S., Garcia, K.E., Gelboin-Burkhart, C., Nguyen, T.B. (2020). Principles of RNA processing from analysis of enhanced CLIP maps for 150 RNA binding proteins. Genome Biol.. 21, 90.
  • Vasudevan, S., Steitz, J.A. (2007). AU-rich-element-mediated upregulation of translation by FXR1 and Argonaute 2. Cell. 128, 1105-1118.
  • Wagner, R.W., Smith, J.E., Cooperman, B.S., Nishikura, K. (1989). A double-stranded RNA unwinding activity introduces structural alterations by means of adenosine to inosine conversions in mammalian cells and Xenopus eggs. Proc. Natl. Acad. Sci. U. S. A.. 86, 2647-2651.
  • Wang, D., Zhang, Z., Jiang, Y., Mao, Z., Wang, D., Lin, H., Xu, D. (2021). DM3Loc: multi-label mRNA subcellular localization prediction and analysis based on multi-head self-attention mechanism. Nucleic Acids Res.. 49, e46.
  • Wang, J.X., Gao, J., Ding, S.L., Wang, K., Jiao, J.Q., Wang, Y., Sun, T., Zhou, L.Y., Long, B., Zhang, X.J. (2015). Oxidative modification of miR-184 enables it to target Bcl-xL and Bcl-w. Mol. Cell. 59, 50-61.
  • Wang, Q., Hui, H., Guo, Z., Zhang, W., Hu, Y., He, T., Tai, Y., Peng, P., Wang, L. (2013). ADAR1 regulates ARHGAP26 gene expression through RNA editing by disrupting miR-30b-3p and miR-573 binding. RNA. 19, 1525-1536.
  • Wen, J., Parker, B.J., Jacobsen, A., Krogh, A. (2011). MicroRNA transfection and AGO-bound CLIP-seq data sets reveal distinct determinants of miRNA action. RNA. 17, 820-834.
  • Wienholds, E., Kloosterman, W.P., Miska, E., Alvarez-Saavedra, E., Berezikov, E., de Bruijn, E., Horvitz, H.R., Kauppinen, S., Plasterk, R.H. (2005). MicroRNA expression in zebrafish embryonic development. Science. 309, 310-311.
  • Wu, J., Bao, J., Kim, M., Yuan, S., Tang, C., Zheng, H., Mastick, G.S., Xu, C., Yan, W. (2014). Two miRNA clusters, miR-34b/c and miR-449, are essential for normal brain development, motile ciliogenesis, and spermatogenesis. Proc. Natl. Acad. Sci. U. S. A.. 111, E2851-E2857.
  • Xin, M., Small, E.M., Sutherland, L.B., Qi, X., McAnally, J., Plato, C.F., Richardson, J.A., Bassel-Duby, R., Olson, E.N. (2009). MicroRNAs miR-143 and miR-145 modulate cytoskeletal dynamics and responsiveness of smooth muscle cells to injury. Genes Dev.. 23, 2166-2178.
  • Xue, Y., Ouyang, K., Huang, J., Zhou, Y., Ouyang, H., Li, H., Wang, G., Wu, Q., Wei, C., Bi, Y. (2013). Direct conversion of fibroblasts to neurons by reprogramming PTB-regulated microRNA circuits. Cell. 152, 82-96.
  • Yang, W., Chendrimada, T.P., Wang, Q., Higuchi, M., Seeburg, P.H., Shiekhattar, R., Nishikura, K. (2006). Modulation of microRNA processing and expression through RNA editing by ADAR deaminases. Nat. Struct. Mol. Biol.. 13, 13-21.
  • Zeng, H., Edwards, M.D., Liu, G., Gifford, D.K. (2016). Convolutional neural network architectures for predicting DNA-protein binding. Bioinformatics. 32, i121-i127.

Figure 1


Structures of guide RNA-loaded human Argonaute protein 2 (hAGO2) before (top left; PDB ID: 4W5N) and after (top right; PDB ID: 4W5O) binding to the target RNA are shown ( Schirle et al., 2014). The structures were rendered using UCSF ChimeraX ( Pettersen et al., 2021). Each nucleotide in the target RNA is numbered according to the position of the corresponding nucleotide in the miRNA. O denotes a Watson-Crick base pairing, Ø denotes a mismatch, B denotes non-A bases, and N denotes any base.

Figure 2


Determinants of miRNA targeting that have already been utilized for target prediction include 3′UTR length, minimum distance to the 3′UTR end, structural accessibility, local AU content, evolutionary conservation, ORF target site, target site abundance, seed pairing stability, thermodynamic pairing stability, the binding affinity between miRNA-loaded AGO and target, 3′UTR isoforms, and 3′ pairing. RISC, RNA-induced silencing complex; AGO, Argonaute; RBNS, RNA Bind-n-Seq.

Figure 3


Determinants that have not been employed in currently available miRNA target prediction tools include TDMD, RNA modifications (RMs) in miRNAs, RMs in target mRNAs, RBP binding sites, and ceRNAs. RISC, RNA-induced silencing complex.

Figure 4


Examples of high-dimensional input data (left), including eCLIP-seq ( Van Nostrand et al., 2016), SHAPE-seq ( Lucks et al., 2011), Nanopore ( Jain et al., 2016), DMS-seq ( Rouskin et al., 2014), and PhyloP ( Pollard et al., 2010), are shown. DMS-seq data was retrieved from GEO GSM1297493 ( Rouskin et al., 2014), eCLIP-seq data from ENCODE ENCFF171MEK ( Van Nostrand et al., 2020), and SHAPE-seq data from RMDB HCIRES_1M7_0001 ( Loughrey et al., 2014). Examples of a deep learning model (middle) and a cartoon representation of a miRNA regulatory network (right) are shown. eCLIP, enhanced cross-linking immunoprecipitation; seq, sequencing; SHAPE, selective 2′-hydroxyl acylation analyzed by primer extension; DMS, dimethyl sulfate.

Table 1

A table of representative computational tools for miRNA target prediction and the determinants they use

Model Seed TPS EC SA Dist. AU Len. 3Sup. TA ORFS
TargetScan7 O SPS O O O O O O O 8m
miRanda-mirSVR O X O O O O O O X X
DIANA-microT-CDS O O O O O O X X X O
MIRZA-G O O O O O X X X X X
PITA Opt. O X O X X X X X X
PicTar O O O X X X X X X X
RNAhybrid Opt. O X X X X X X X X
MicroTar O O X X X X X X X X