Next Article in Journal
Genome-Wide Identification of microRNAs in Response to Salt/Alkali Stress in Medicago truncatula through High-Throughput Sequencing
Next Article in Special Issue
Focus on Translation Initiation of the HIV-1 mRNAs
Previous Article in Journal
Roco Proteins and the Parkinson’s Disease-Associated LRRK2
Previous Article in Special Issue
Eukaryotic Initiation Factor 5B (eIF5B) Cooperates with eIF1A and eIF5 to Facilitate uORF2-Mediated Repression of ATF4 Translation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

miRNAmotif—A Tool for the Prediction of Pre-miRNA–Protein Interactions

Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61-704 Poznan, Poland
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2018, 19(12), 4075; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19124075
Submission received: 16 November 2018 / Revised: 12 December 2018 / Accepted: 13 December 2018 / Published: 17 December 2018
(This article belongs to the Special Issue Translational Control)

Abstract

:
MicroRNAs (miRNAs) are short, non-coding post-transcriptional gene regulators. In mammalian cells, mature miRNAs are produced from primary precursors (pri-miRNAs) using canonical protein machinery, which includes Drosha/DGCR8 and Dicer, or the non-canonical mirtron pathway. In plant cells, mature miRNAs are excised from pri-miRNAs by the DICER-LIKE1 (DCL1) protein complex. The involvement of multiple regulatory proteins that bind directly to distinct miRNA precursors in a sequence- or structure-dependent manner adds to the complexity of the miRNA maturation process. Here, we present a web server that enables searches for miRNA precursors that can be recognized by diverse RNA-binding proteins based on known sequence motifs to facilitate the identification of other proteins involved in miRNA biogenesis. The database used by the web server contains known human, murine, and Arabidopsis thaliana pre-miRNAs. The web server can also be used to predict new RNA-binding protein motifs based on a list of user-provided sequences. We show examples of miRNAmotif applications, presenting precursors that contain motifs recognized by Lin28, MCPIP1, and DGCR8 and predicting motifs within pre-miRNA precursors that are recognized by two DEAD-box helicases—DDX1 and DDX17. miRNAmotif is released as an open-source software under the MIT License. The code is available at GitHub (www.github.com/martynaut/mirnamotif). The webserver is freely available at http://mirnamotif.ibch.poznan.pl.

1. Introduction

MicroRNAs (miRNAs) are short, non-coding post-transcriptional gene regulators that control many fundamental cellular processes, such as metabolism, cellular proliferation, apoptosis, immune function, epigenetics, and neurodevelopment [1,2]. The role of microRNA in pathological processes, including human diseases and cancer, is also increasingly acknowledged [3,4,5,6]. Many computational methods are utilized to better understand miRNAs function [7,8,9,10,11,12,13].
Most mammalian miRNAs are produced in the canonical miRNA biogenesis pathway. This pathway includes the excision of mature miRNA from the primary transcript (pri-miRNA) by sequential nuclear (Microprocessor complex-mediated) and cytoplasmic (Dicer-mediated) processing through the intermediate miRNA precursor (pre-miRNA) [14,15,16]. Increasing numbers of other cellular proteins and RNAs (including long non-coding RNAs) are shown to influence miRNA biogenesis and function [17,18,19,20], adding to the complexity of the miRNA pathway. Interestingly, some of the RNA-binding proteins (RBPs) known to bind mRNAs can also interact with miRNA precursors, creating indirect crosstalk between mRNA and miRNA [21].
The characteristic sequence motifs in the terminal loop or stem portion of pre-miRNA can be selectively recognized by RBPs, which can regulate miRNA production and play an important role in RNA metabolism (Figure 1) [22]. Known protein regulators of miRNA biogenesis were recently reviewed by Michlewski et al. [23]. The protein domains that recognize characteristic sequences in the pre-miRNA terminal loop bind mainly to short single-stranded RNA (ssRNA) regions. These proteins perform various functions in miRNA biogenesis, including miRNA precursor relocation, recruitment, or displacement of other proteins with catalytic activity, or the induction of changes in RNA secondary structures [24]. Different proteins may compete to bind to the same terminal loop [25]. Known sequence motifs found in the terminal loop or stem portion of miRNA precursors are shown in Table 1.
Predicted motifs within miRNA precursors may be validated using RNAs containing mutations within potential protein binding sequences [34]; however, the current method of choice is high-throughput crosslinking and immunoprecipitation followed by sequencing (CLIP-seq) to find the exact sequences interacting with candidate proteins [39,40]. For example, DDX17 was demonstrated by CLIP-seq to bind to multiple miRNA precursors [35].
Adding to the complexity of miRNA biogenesis, multiple SNPs and mutations were identified within miRNA precursors influencing their processing and/or function [41,42]. Small sequence changes may influence the binding of proteins due to recognition motif disruption/creation or structure rearrangement [43,44]. Small sequence changes in microRNA precursors may also influence Drosha and Dicer cleavage sites and the homogeneity of produced mature miRNA [45,46].
miRNAmotif is the first application that enables the user (i) to search for specific motifs in the terminal loop, mature miRNA, linking sequence (sequence between mature miRNAs from two arms), or in whole pre-miRNA among sequences in a widely used miRNA database (miRBase) and (ii) to predict enriched motifs within miRNA precursors’ sequences provided by the user. miRNAmotif is a simple and user-friendly application that gives researchers the opportunity to analyze known pre-miRNAs sequences for the presence of specific motifs that can be recognized by protein regulators of miRNA biogenesis.

2. Results

To show the applicability of the miRNAmotif tool for specific biological problems, we analyzed publicly available data to identify pre-miRNAs that may be regulated by Lin28 (Section 2.1.1), MCPIP1, and DGCR8 (Section 2.1.2), and to find the motifs that are potentially recognized by the DDX1 and DDX17 proteins (Section 2.2.1).

2.1. Searches for Known Motifs within Pre-miRNA Sequences

2.1.1. Human miRNAs with the Lin28 Motif

We used miRNAmotif to identify pre-miRNAs that contain the GGAG motif within their loop or linking sequence recognized by Lin28 in human miRNAs.
First, we searched for the GGAG motif in the 5′→3′ direction only in the terminal loop. We found 70 pre-miRNAs potentially bound by the Lin28 protein (Supplementary Table S1). Next, we searched for the GGAG motif in the linking sequence in the same direction. We found 155 pre-miRNAs that have the GGAG motif within the sequence between mature miRNAs (Supplementary Table S2). As Lin28 is known to regulate the let-7 family of miRNAs and miRNA-9 [47,48], we investigated whether the miRNAmotif algorithm found the motif in the corresponding precursors. In the first search (loop only), we found hsa-let-7a-2 and hsa-mir-9-1. In the linking sequence search, we showed that in addition to these two miRNAs, the GGAG motif was found in other miRNAs from the let-7 family: hsa-let-7a-1, hsa-let-7c, hsa-let-7d, hsa-let-7e, hsa-let-7f-1, hsa-let-7f-2, hsa-let-7g, and hsa-let-7i.
Next, we analyzed the GGAG motif occurrence in whole pre-miRNA sequences of human pre-miRNAs. As expected, we obtained a long list (793) of pre-miRNAs containing the simple four-letter motif in the 5’→3’ orientation (results not shown). Pre-miRNAs found in the loop-based and the linking-sequence-based search were included in the long list obtained in the whole precursor-based search. However, as Lin28 binds to single-strand regions, most of these pre-miRNAs (with the GGAG motif within a stem) are probably false positive.
Lin28 proteins were previously analyzed in the context of pre-miRNA interactions in multiple model organisms, as these highly conserved proteins are involved in the regulation of pluripotency in stem cells. It was shown in Xenopus embryos that Lin28 can also regulate other miRNAs in addition to let-7 family, such as mir-17∼92 and mir-106∼363 cluster miRNAs [49]. One of the regulated miRNAs was xla-mir-363. Analogously, we found hsa-mir-363 to be potentially interacting with Lin28 using miRNAmotif. Another pre-miRNA precursor found in Xenopus research and miRNAmotif results is mir-200b, also expressed in miRNA cluster. These results show that other detected pre-miRNAs, besides the known let-7 family and miRNA-9, may interact with Lin28.
Lin28 proteins are found mainly in testis and placenta tissues [50]. Using miRmine, an miRNA tissue-specific expression database [51], we searched for miRNAs that had expression in testis to further limit the list of found pre-miRNAs that can interact with Lin28. As expected, all let-7 family miRNAs had detectable expression in testis according to miRmine. We found that 40 (Table 2) out of 151 pre-miRNAs having the GGAG motif within a linking sequence had expression in testis. Four of the analyzed pre-miRNAs were not included in the miRmine database. In placenta, we found 36 pre-miRNAs showing detectable expression and the set of two groups (testis and placenta) contained 43 unique pre-miRNAs.

2.1.2. Human miRNAs Containing Highly Similar Motifs

SNPs and mutations are widely known to influence miRNAs processing and function. To show that miRNAmotif can be a useful tool for the precise analysis of motifs differing in one nucleotide, we used motifs recognized by MCPIP1 and DGCR8 proteins (UGC and UGU, respectively) as examples.
First, we searched for both UGC and UGU motifs together with an overlap option in the terminal loop only. We found 244 pre-miRNAs with the UGC motif, which may be recognized by the MCPIP1 protein (Supplementary Table S3a), and 315 pre-miRNAs with the UGU motif, which may be recognized by DGCR8 (Supplementary Table S3b). We also found 47 pre-miRNAs containing both UGC and UGU motifs in their terminal loops (Supplementary Table S3c). Then, using the module to search a single motif, we found 606 pre-miRNAs containing the UGY motif in their terminal loops. This number was equal to that of the first search and, more importantly, miRNAs found by searching for the UGY motif were exactly the same as those found previously using two motifs. This result shows that miRNAmotif is a useful tool for the precise analysis of short motifs with a single-nucleotide difference as well.

2.2. Predicting Protein-Binding Motifs within Pre-miRNA Sequences

2.2.1. Common Motifs in Pre-miRNAs Interacting with DEAD-Box Helicases

Next, we focused on the search of new motifs potentially interacting with DEAD-box helicases. To perform the first analysis, we used the data generated for pre-miRNAs interactions with one of the DEAD-box helicases, DDX1 [34]. A total of 25 pre-miRNA sequences were analyzed (Table 3).
We searched for the motif in the whole pre-miRNA sequence in the 5’→3’ direction. It was shown that DDX proteins (including DDX1) can interact with sequences outside the terminal loop of the miRNA precursors. As a set of negative sequences, we used all human pre-miRNA sequences. Our software detected five significantly enriched motifs—CUAAYACU (p = 1.7 × 10−7), CGYUUU (p = 5.5 × 10−7), CUGKUAA (p = 1.4 × 10−6), CRUCUUAC (p = 1.4 × 10−6), and AAUCRU (p = 1.1 × 10−5) (Supplementary Table S4). Three out of five found motifs harbor the dinucleotide AA detected by Han et al. [34]; however, the miRNAmotif results showed that this short two-nucleotide motif might not be adequate or sufficient for protein binding. The miRNAmotif tool also detected a single significantly enriched motif within the terminal loop in the 5’→3’ direction with all human loops used as negative sequences (UUUAU, p = 9.9 × 10−5).
In the next step, we analyzed which miRNAs outside of the experimental group of miRNA precursors could be recognized by DDX1 based on the most significant predicted motif, CUAAYACU. We found seven human miRNA precursors, of which five were included in the experimental group (Supplementary Table S5).
To validate the miRNAmotif function for longer known motifs, we focused on DDX17. We used the data generated for DDX17 interaction with pre-miRNAs. A total of 153 pre-miRNAs sequences were analyzed [35]. Seven pre-miRNAs from the original list (hsa-miR-3648/3687/6087/7641-1/7641-2/3792/3653) had to be excluded from the analysis, as they are no longer included in the new version of the miRBase database or may relate to multiple records in miRBase.
First, miRNAmotif was used to search for motifs located within the terminal loop. The software did not find any significant motifs enriched in the terminal loops of provided sequences in reference to all human pre-miRNA terminal loops. Similarly, no significantly enriched motifs were found in linking sequences. A second search was performed using whole pre-miRNA sequences. Using these settings, the miRNAmotif algorithm detected two significantly enriched motifs, SSGGG (p = 6.8 × 10−9) and CCSS (p = 4.7 × 10−7). Next, we analyzed which precursors outside of the experimental group of pre-miRNAs contain the most significantly enriched SSGGG motif. We found 26 pre-miRNAs (with three included in the experimental group) with the SSGGG motif in their terminal loops and 410 pre-miRNAs containing the searched motif in the whole precursor sequence. It is worth mentioning that Moy et al. did not find any enriched motif for DDX17 using CLIP-seq data in their study; however, they restricted their analysis to hexamers only [35]. This proves that miRNAmotif may expand previously performed analyses.

3. Discussion

The identification of RNA–protein interactions is a key step towards improving our understanding of the post-transcriptional regulation of gene expression. The number of known motifs recognized by various RNA-binding proteins is increasing due to the rapid development and application of high-throughput sequencing techniques. These motifs are usually considered in the context of mRNA sequences, but can be also ubiquitous in non-coding RNAs, including long non-coding RNAs and miRNA precursors.
Protein binding, specifically to a group of pre-miRNAs, may influence the efficiency of miRNA biogenesis and explain the differences in cellular levels of various miRNAs. Furthermore, specific protein binding can help improve our understanding of tissue- and cell-specific miRNA-mediated mRNA regulation, as was recently highlighted in two protein–miRNA interaction studies [33,52]. With increased knowledge of miRNA and RBP tissue-specific expression levels, in both physiological and pathological states [53], we can better evaluate which interactions may actually happen and in what cell types, as we showed for Lin28 interactions. In addition to cell- and tissue-specificity, in predicting pre-miRNA–protein interactions, the subcellular localization of protein and pre-miRNA should be also taken into account. In particular, non-ubiquitously expressed proteins that are localized in specific cellular compartments where pre-miRNAs are not present would not interact with these RNA molecules even if a recognized motif was found in the miRNA precursor.
We acknowledge that our software has several limitations. In its current form, it does not take into account the subcellular localization of RBPs and tissue-specific expression levels of specific miRNAs. Future improvements of the algorithm may therefore include a tissue-specific version of the tool, as was previously shown in splicing research [54]. Moreover, like all in silico analyses, miRNAmotif results need to be verified experimentally. This applies especially to short (e.g., two-nucleotide) motifs that may be found in the long list of pre-miRNA precursors. The miRNAmotif tool aims to facilitate and ease the search for miRNA precursors with specific sequence motifs for further evaluation.
Taken together, miRNAmotif is a useful and easy-to-employ application that may increase our understanding of miRNA biogenesis by searching for protein-interacting motifs in pre-miRNA sequences or predicting new motifs in a user-provided group of pre-miRNAs. The miRNAmotif software will be updated with the new releases of the miRBase database.

4. Materials and Methods

miRNAmotif is released as an open-source software under the MIT License. The code is available at GitHub (www.github.com/martynaut/mirnamotif). The web server is freely available at http://mirnamotif.ibch.poznan.pl. miRNAmotif was written in Python programming language using the Pycherry web framework. The software can be used through a user-friendly web interface (Figure 2) or using a command line interface (CLI).
miRNAmotif used sequences of pre-miRNAs extracted from miRBase database version 22 [55]. miRNAs found in humans (hsa), mice (mmu), and Arabidopsis thaliana (ath) were included. The mfold algorithm for structure prediction was used to predict secondary structures of miRNA precursors and define the terminal loop portion with a 2-nt overhang [56]; the same algorithm was used to predict the structures presented in miRBase. Linking sequences (sequences between mature miRNAs from two arms) were extracted based on the genomic coordinates of mature miRNAs and their precursors from miRbase. The mfold algorithm was used to predict secondary structures of miRNA precursors and define the passenger strand for pre-miRNAs with defined mature miRNA from a single arm only.
For motif definition, we used IUPAC (International Union of Pure and Applied Chemistry) nucleic acid nomenclature, where R = G or A, Y = U or C, M = A or C, K = G or U, S = G or C, W = A or U, H = A or C or U, B = G or U or C, V = G or C or A, D = G or U or A, and N = G or U or A or C.
Three independent modules are available to the user. The first module is used to search for sequences containing one user-provided motif. The user may choose the organism database of pre-miRNAs, whether the motif should be located within loop, linking sequence, mature miRNA, or anywhere in the pre-miRNA (default: loop), and whether the motif should be searched in one direction (forward, 5′→3′) or in both the 5′→3′ and reverse directions (default: 5′→3′). The .txt and .rft files with sequences that contain the requested motif (marked in red in the .rtf file) are then sent to the user at the indicated e-mail address and printed out on the webpage.
The second module works analogously to the first, but the user can search for two motifs simultaneously. In addition to the options indicated earlier, the user may exclude miRNAs in which two motifs overlap (default: enable overlapping). As a result, the .txt and .rtf files with pre-miRNA, mature miRNA, and linking or loop sequences that contain a single or both motifs (marked in red in the .rtf file) are sent to the user and printed out on the webpage.
The third module enables the identification of new motifs common to a group of miRNAs that are shown to interact with the chosen protein or that exhibit similar processing. The user needs to paste sequences or the names of miRNAs (in miRBase nomenclature, e.g., hsa-mir-21) into the window or upload a text file with the sequences or names of miRNAs. An exemplary input file is provided for download on the webpage. The motif finder has a similar option as the two above modules, in which the user may find motifs in the 5’→3’ direction or in both the 5’→3’ direction and the reverse direction (default: 5’→3’) and can choose whether the motif needs to be found in the terminal loop, linking sequence, mature miRNA, or in a whole pre-miRNA sequence (default: loop). As a result, the user receives a .txt file with motifs and .eps files with motif representations in Weblogo [57]. For motif prediction, we used the DREME algorithm from the MEME suite [58].

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/19/12/4075/s1.

Author Contributions

Conceptualization, M.O.U.-T. and E.J.; Software, M.O.U.-T.; Supervision, W.J.K.; Validation, M.O.U.-T. and E.J.; Writing—original draft, M.O.U.-T. and E.J.; and Writing—review and editing, M.O.U.-T., E.J., and W.J.K.

Funding

This research was funded by the National Science Centre (2015/17/N/NZ3/03629 to M.O.U.-T., 2012/06/A/NZ1/00094 to W.J.K., and 2015/19/B/NZ2/02453 to W.J.K.); the L’Oreal Poland Women in Science Fellowship (to M.O.U.-T.); the Polish Ministry of Science and Higher Education (statutory funds for young investigators to E.J.); and the Polish Ministry of Science and Higher Education under the KNOW program.

Acknowledgments

We dedicate this work to the recently deceased Wlodzimierz J. Krzyzosiak.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

miRNAmicroRNA
DDXDEAD box proteins
RSArginine/serine-rich domain
RRMRNA recognition motif
RBPRNA-binding protein
IUPACInternational Union of Pure and Applied Chemistry
CLIPCrosslinking immunoprecipitation
CLICommand line interface

References

  1. Ameres, S.L.; Zamore, P.D. Diversifying microRNA sequence and function. Nat. Rev. Mol. Cell Biol. 2013, 14, 475–488. [Google Scholar] [CrossRef] [PubMed]
  2. Bartel, D.P. Metazoan MicroRNAs. Cell 2018, 173, 20–51. [Google Scholar] [CrossRef] [PubMed]
  3. Svoronos, A.A.; Engelman, D.M.; Slack, F.J. OncomiR or Tumor Suppressor? The Duplicity of MicroRNAs in Cancer. Cancer Res. 2016, 76, 3666–3670. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Garo, L.P.; Murugaiyan, G. Contribution of MicroRNAs to autoimmune diseases. Cell. Mol. Life Sci. 2016, 73, 2041–2051. [Google Scholar] [CrossRef] [PubMed]
  5. Zeng, X.; Liu, L.; Lü, L.; Zou, Q. Prediction of potential disease-associated microRNAs using structural perturbation method. Bioinforma. Oxf. Engl. 2018, 34, 2425–2432. [Google Scholar] [CrossRef]
  6. Tang, W.; Wan, S.; Yang, Z.; Teschendorff, A.E.; Zou, Q. Tumor origin detection with tissue-specific miRNA and DNA methylation markers. Bioinforma. Oxf. Engl. 2018, 34, 398–406. [Google Scholar] [CrossRef] [PubMed]
  7. Zou, Q.; Mao, Y.; Hu, L.; Wu, Y.; Ji, Z. miRClassify: An advanced web server for miRNA family classification and annotation. Comput. Biol. Med. 2014, 45, 157–160. [Google Scholar] [CrossRef] [PubMed]
  8. Li, Y.; Zhang, Z. Computational Biology in microRNA. Wiley Interdiscip. Rev. RNA 2015, 6, 435–452. [Google Scholar] [CrossRef] [PubMed]
  9. Jiang, H.; Wang, J.; Li, M.; Lan, W.; Wu, F.; Pan, Y. miRTRS: A Recommendation Algorithm for Predicting miRNA Targets. IEEE/ACM Trans. Comput. Biol. Bioinform. 2018. [Google Scholar] [CrossRef] [PubMed]
  10. Andrés-León, E.; Rojas, A.M. miARma-Seq, a comprehensive pipeline for the simultaneous study and integration of miRNA and mRNA expression data. Methods San Diego Calif. 2018. [Google Scholar] [CrossRef] [PubMed]
  11. Huang, Z.; Shi, J.; Gao, Y.; Cui, C.; Zhang, S.; Li, J.; Zhou, Y.; Cui, Q. HMDD v3.0: A database for experimentally supported human microRNA-disease associations. Nucleic Acids Res. 2018. [Google Scholar] [CrossRef] [PubMed]
  12. Chen, X.; Wang, C.-C.; Yin, J.; You, Z.-H. Novel Human miRNA-Disease Association Inference Based on Random Forest. Mol. Ther. Nucleic Acids 2018, 13, 568–579. [Google Scholar] [CrossRef] [PubMed]
  13. Liu, W.; Cui, Z.; Zan, X. Identifying cancer-related microRNAs based on subpathways. IET Syst. Biol. 2018, 12, 273–278. [Google Scholar] [CrossRef] [PubMed]
  14. Stroynowska-Czerwinska, A.; Fiszer, A.; Krzyzosiak, W.J. The panorama of miRNA-mediated mechanisms in mammalian cells. Cell. Mol. Life Sci. CMLS 2014, 71, 2253–2270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Ha, M.; Kim, V.N. Regulation of microRNA biogenesis. Nat. Rev. Mol. Cell Biol. 2014, 15, 509–524. [Google Scholar] [CrossRef]
  16. Creugny, A.; Fender, A.; Pfeffer, S. Regulation of primary microRNA processing. FEBS Lett. 2018, 592, 1980–1996. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Bayoumi, A.S.; Sayed, A.; Broskova, Z.; Teoh, J.-P.; Wilson, J.; Su, H.; Tang, Y.-L.; Kim, I. Crosstalk between Long Noncoding RNAs and MicroRNAs in Health and Disease. Int. J. Mol. Sci. 2016, 17, 356. [Google Scholar] [CrossRef] [PubMed]
  18. Witkos, T.M.; Krzyzosiak, W.J.; Fiszer, A.; Koscianska, E. A potential role of extended simple sequence repeats in competing endogenous RNA crosstalk. RNA Biol. 2018, 1–11. [Google Scholar] [CrossRef]
  19. Ratnadiwakara, M.; Mohenska, M.; Änkö, M.-L. Splicing factors as regulators of miRNA biogenesis - links to human disease. Semin. Cell Dev. Biol. 2018, 79, 113–122. [Google Scholar] [CrossRef]
  20. Gebert, L.F.R.; MacRae, I.J. Regulation of microRNA function in animals. Nat. Rev. Mol. Cell Biol. 2018. [Google Scholar] [CrossRef]
  21. Treiber, T.; Treiber, N.; Meister, G. Regulation of microRNA biogenesis and its crosstalk with other cellular pathways. Nat. Rev. Mol. Cell Biol. 2018. [Google Scholar] [CrossRef] [PubMed]
  22. Krol, J.; Loedige, I.; Filipowicz, W. The widespread regulation of microRNA biogenesis, function and decay. Nat. Rev. Genet. 2010, 11, 597–610. [Google Scholar] [CrossRef] [PubMed]
  23. Michlewski, G.; Caceres, J.F. Post-transcriptional control of miRNA biogenesis. RNA 2018, rna.068692.118. [Google Scholar] [CrossRef] [PubMed]
  24. Connerty, P.; Ahadi, A.; Hutvagner, G. RNA Binding Proteins in the miRNA Pathway. Int. J. Mol. Sci. 2016, 17, 31. [Google Scholar] [CrossRef] [PubMed]
  25. Castilla-Llorente, V.; Nicastro, G.; Ramos, A. Terminal loop-mediated regulation of miRNA biogenesis: Selectivity and mechanisms. Biochem. Soc. Trans. 2013, 41, 861–865. [Google Scholar] [CrossRef] [PubMed]
  26. Burd, C.G.; Dreyfuss, G. RNA binding specificity of hnRNP A1: Significance of hnRNP A1 high-affinity binding sites in pre-mRNA splicing. EMBO J. 1994, 13, 1197–1204. [Google Scholar] [CrossRef]
  27. Myer, V.E.; Fan, X.C.; Steitz, J.A. Identification of HuR as a protein implicated in AUUUA-mediated mRNA decay. EMBO J. 1997, 16, 2130–2139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Nicastro, G.; García-Mayoral, M.F.; Hollingworth, D.; Kelly, G.; Martin, S.R.; Briata, P.; Gherzi, R.; Ramos, A. Noncanonical G recognition mediates KSRP regulation of let-7 biogenesis. Nat. Struct. Mol. Biol. 2012, 19, 1282–1286. [Google Scholar] [CrossRef] [Green Version]
  29. Nam, Y.; Chen, C.; Gregory, R.I.; Chou, J.J.; Sliz, P. Molecular basis for interaction of let-7 microRNAs with Lin28. Cell 2011, 147, 1080–1091. [Google Scholar] [CrossRef]
  30. Goers, E.S.; Purcell, J.; Voelker, R.B.; Gates, D.P.; Berglund, J.A. MBNL1 binds GC motifs embedded in pyrimidines to regulate alternative splicing. Nucleic Acids Res. 2010, 38, 2467–2484. [Google Scholar] [CrossRef] [Green Version]
  31. Suzuki, H.I.; Arase, M.; Matsuyama, H.; Choi, Y.L.; Ueno, T.; Mano, H.; Sugimoto, K.; Miyazono, K. MCPIP1 ribonuclease antagonizes dicer and terminates microRNA biogenesis through precursor microRNA degradation. Mol. Cell 2011, 44, 424–436. [Google Scholar] [CrossRef] [PubMed]
  32. Nguyen, T.A.; Jo, M.H.; Choi, Y.-G.; Park, J.; Kwon, S.C.; Hohng, S.; Kim, V.N.; Woo, J.-S. Functional Anatomy of the Human Microprocessor. Cell 2015, 161, 1374–1387. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Treiber, T.; Treiber, N.; Plessmann, U.; Harlander, S.; Daiß, J.-L.; Eichner, N.; Lehmann, G.; Schall, K.; Urlaub, H.; Meister, G. A Compendium of RNA-Binding Proteins that Regulate MicroRNA Biogenesis. Mol. Cell 2017, 66, 270–284.e13. [Google Scholar] [CrossRef] [PubMed]
  34. Han, C.; Liu, Y.; Wan, G.; Choi, H.J.; Zhao, L.; Ivan, C.; He, X.; Sood, A.K.; Zhang, X.; Lu, X. The RNA-binding protein DDX1 promotes primary microRNA maturation and inhibits ovarian tumor progression. Cell Rep. 2014, 8, 1447–1460. [Google Scholar] [CrossRef] [PubMed]
  35. Moy, R.H.; Cole, B.S.; Yasunaga, A.; Gold, B.; Shankarling, G.; Varble, A.; Molleston, J.M.; TenOever, B.R.; Lynch, K.W.; Cherry, S. Stem-loop recognition by DDX17 facilitates miRNA processing and antiviral defense. Cell 2014, 158, 764–777. [Google Scholar] [CrossRef] [PubMed]
  36. Chen, Y.; Zubovic, L.; Yang, F.; Godin, K.; Pavelitz, T.; Castellanos, J.; Macchi, P.; Varani, G. Rbfox proteins regulate microRNA biogenesis by sequence-specific binding to their precursors and target downstream Dicer. Nucleic Acids Res. 2016, 44, 4381–4395. [Google Scholar] [CrossRef] [PubMed]
  37. Davis, B.N.; Hilyard, A.C.; Nguyen, P.H.; Lagna, G.; Hata, A. Smad proteins bind a conserved RNA sequence to promote microRNA maturation by Drosha. Mol. Cell 2010, 39, 373–384. [Google Scholar] [CrossRef]
  38. Kim, K.; Nguyen, T.D.; Li, S.; Nguyen, T.A. SRSF3 recruits DROSHA to the basal junction of primary microRNAs. RNA N. Y. 2018, 24, 892–898. [Google Scholar] [CrossRef]
  39. Stork, C.; Zheng, S. Genome-Wide Profiling of RNA-Protein Interactions Using CLIP-Seq. Methods Mol. Biol. 2016, 1421, 137–151. [Google Scholar]
  40. Yang, Y.-C.T.; Di, C.; Hu, B.; Zhou, M.; Liu, Y.; Song, N.; Li, Y.; Umetsu, J.; Lu, Z.J. CLIPdb: A CLIP-seq database for protein-RNA interactions. BMC Genom. 2015, 16, 51. [Google Scholar] [CrossRef]
  41. Króliczewski, J.; Sobolewska, A.; Lejnowski, D.; Collawn, J.F.; Bartoszewski, R. microRNA single polynucleotide polymorphism influences on microRNA biogenesis and mRNA target specificity. Gene 2018, 640, 66–72. [Google Scholar] [CrossRef] [PubMed]
  42. Moszyńska, A.; Gebert, M.; Collawn, J.F.; Bartoszewski, R. SNPs in microRNA target sites and their potential role in human disease. Open Biol. 2017, 7, 170019. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Roden, C.; Gaillard, J.; Kanoria, S.; Rennie, W.; Barish, S.; Cheng, J.; Pan, W.; Liu, J.; Cotsapas, C.; Ding, Y.; et al. Novel determinants of mammalian primary microRNA processing revealed by systematic evaluation of hairpin-containing transcripts and human genetic variation. Genome Res. 2017, 27, 374–384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Fernandez, N.; Cordiner, R.A.; Young, R.S.; Hug, N.; Macias, S.; Cáceres, J.F. Genetic variation and RNA structure regulate microRNA biogenesis. Nat. Commun. 2017, 8, 15114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Starega-Roslan, J.; Galka-Marciniak, P.; Krzyzosiak, W.J. Nucleotide sequence of miRNA precursor contributes to cleavage site selection by Dicer. Nucleic Acids Res. 2015, 43, 10939–10951. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Galka-Marciniak, P.; Olejniczak, M.; Starega-Roslan, J.; Szczesniak, M.W.; Makalowska, I.; Krzyzosiak, W.J. siRNA release from pri-miRNA scaffolds is controlled by the sequence and structure of RNA. Biochim. Biophys. Acta 2016, 1859, 639–649. [Google Scholar] [CrossRef] [PubMed]
  47. Nowak, J.S.; Hobor, F.; Downie Ruiz Velasco, A.; Choudhury, N.R.; Heikel, G.; Kerr, A.; Ramos, A.; Michlewski, G. Lin28a uses distinct mechanisms of binding to RNA and affects miRNA levels positively and negatively. RNA 2017, 23, 317–332. [Google Scholar] [CrossRef] [PubMed]
  48. Ustianenko, D.; Chiu, H.-S.; Treiber, T.; Weyn-Vanhentenryck, S.M.; Treiber, N.; Meister, G.; Sumazin, P.; Zhang, C. LIN28 Selectively Modulates a Subclass of Let-7 MicroRNAs. Mol. Cell 2018, 71, 271–283.e5. [Google Scholar] [CrossRef] [PubMed]
  49. Warrander, F.; Faas, L.; Kovalevskiy, O.; Peters, D.; Coles, M.; Antson, A.A.; Genever, P.; Isaacs, H.V. lin28 proteins promote expression of 17∼92 family miRNAs during amphibian development. Dev. Dyn. 2016, 245, 34–46. [Google Scholar] [CrossRef]
  50. Berglund, L.; Björling, E.; Oksvold, P.; Fagerberg, L.; Asplund, A.; Szigyarto, C.A.-K.; Persson, A.; Ottosson, J.; Wernérus, H.; Nilsson, P.; et al. A genecentric Human Protein Atlas for expression profiles based on antibodies. Mol. Cell. Proteom. MCP 2008, 7, 2019–2027. [Google Scholar] [CrossRef]
  51. Panwar, B.; Omenn, G.S.; Guan, Y. miRmine: A database of human miRNA expression profiles. Bioinformatics 2017, 33, 1554–1560. [Google Scholar] [CrossRef] [PubMed]
  52. Nussbacher, J.K.; Yeo, G.W. Systematic Discovery of RNA Binding Proteins that Regulate MicroRNA Levels. Mol. Cell 2018, 69, 1005–1016.e7. [Google Scholar] [CrossRef] [PubMed]
  53. Sakaue, S.; Hirata, J.; Maeda, Y.; Kawakami, E.; Nii, T.; Kishikawa, T.; Ishigaki, K.; Terao, C.; Suzuki, K.; Akiyama, M.; et al. Integration of genetics and miRNA-target gene network identified disease biology implicated in tissue specificity. Nucleic Acids Res. 2018. [Google Scholar] [CrossRef] [PubMed]
  54. Piva, F.; Giulietti, M.; Burini, A.B.; Principato, G. SpliceAid 2: A database of human splicing factors expression data and RNA target motifs. Hum. Mutat. 2012, 33, 81–85. [Google Scholar] [CrossRef] [PubMed]
  55. Kozomara, A.; Griffiths-Jones, S. miRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014, 42, D68–D73. [Google Scholar] [CrossRef]
  56. Zuker, M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31, 3406–3415. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Crooks, G.E.; Hon, G.; Chandonia, J.-M.; Brenner, S.E. WebLogo: A sequence logo generator. Genome Res. 2004, 14, 1188–1190. [Google Scholar] [CrossRef]
  58. Bailey, T.L. DREME: Motif discovery in transcription factor ChIP-seq data. Bioinforma. Oxf. Engl. 2011, 27, 1653–1659. [Google Scholar] [CrossRef]
Figure 1. RNA-binding proteins (RBPs) may interact with microRNA (miRNA) precursors through sequence motifs found within the terminal loop (like competing LIN28 and MBNL1) or other elements (like SRSF3). RBPs binding to pre-miRNA precursors may positively or negatively regulate miRNA biogenesis resulting in increased or decreased mature miRNA levels, respectively. Alternatively, protein binding may influence the secondary structure of the miRNA precursor, leading to changed specificity of the Drosha or Dicer cleavage.
Figure 1. RNA-binding proteins (RBPs) may interact with microRNA (miRNA) precursors through sequence motifs found within the terminal loop (like competing LIN28 and MBNL1) or other elements (like SRSF3). RBPs binding to pre-miRNA precursors may positively or negatively regulate miRNA biogenesis resulting in increased or decreased mature miRNA levels, respectively. Alternatively, protein binding may influence the secondary structure of the miRNA precursor, leading to changed specificity of the Drosha or Dicer cleavage.
Ijms 19 04075 g001
Figure 2. Web server interface enabling three modules available to the user.
Figure 2. Web server interface enabling three modules available to the user.
Ijms 19 04075 g002
Table 1. Known motifs in miRNA precursors recognized by protein regulators.
Table 1. Known motifs in miRNA precursors recognized by protein regulators.
ProteinMotifProtein Domain TypeReference
Terminal LoophnRNPA1UAGGGAWRRM[26]
HuRAUUUUUAUUUURRM[27]
KSRPGGGUKH[28]
Lin28GGAGCSD&CCHC-ZnF[29]
MBNL1YGCYCCCH-ZnF[30]
MCPIP1UGCCCCH-ZnF[31]
DGCR8UGURhed[32]
MATR3AUCUUnd[33]
ZC3H7SMTANYnd
YBX1CAUCnd
TRIM71UAUAAnd
PTBP1/3UUUUUCCNUCUUUnd
StemDDX1AAnd[34]
DDX17VCAUCHDEAD box[35]
RbfoxGCAUGRRM[36]
SMADCAGACMH1[37]
CELF1/2UGUNNNNNNNUGUnd[33]
ZC3H10GCAGCGCnd
SRSF3CNNCRS/RRM[38]
nd—no data.
Table 2. Forty pre-miRNAs having the GGAG motif within a linking sequence that are expressed in testis.
Table 2. Forty pre-miRNAs having the GGAG motif within a linking sequence that are expressed in testis.
pre-miRNA
hsa-let-7a-1hsa-miR-204
hsa-let-7a-2hsa-miR-30e
hsa-let-7chsa-miR-320a
hsa-let-7dhsa-miR-324
hsa-let-7ehsa-miR-363
hsa-let-7f-1hsa-miR-3655
hsa-let-7f-2hsa-miR-378i
hsa-let-7ghsa-miR-4286
hsa-let-7ihsa-miR-4784
hsa-miR-107hsa-miR-5006
hsa-miR-1236hsa-miR-5579
hsa-miR-132hsa-miR-5690
hsa-miR-1323hsa-miR-629
hsa-miR-139hsa-miR-6775
hsa-miR-142hsa-miR-6777
hsa-miR-143hsa-miR-6834
hsa-miR-149hsa-miR-7154
hsa-miR-152hsa-miR-7161
hsa-miR-200bhsa-miR-7847
hsa-miR-200chsa-miR-7975
Table 3. Twenty-five pre-miRNAs interacting with DDX1 based on a work by Han et al. [34].
Table 3. Twenty-five pre-miRNAs interacting with DDX1 based on a work by Han et al. [34].
pre-miRNA
hsa-miR-82hsa-miR-410
hsa-miR-96hsa-miR-429
hsa-miR-101hsa-miR-449a
hsa-miR-129hsa-miR-487b
hsa-miR-138hsa-miR-490
hsa-miR-141hsa-miR-495
hsa-miR-146bhsa-miR-499
hsa-miR-155hsa-miR-518e
hsa-miR-200ahsa-miR-524
hsa-miR-200bhsa-miR-539
hsa-miR-200chsa-miR-542
hsa-miR-376hsa-miR-590
hsa-miR-376a

Share and Cite

MDPI and ACS Style

Urbanek-Trzeciak, M.O.; Jaworska, E.; Krzyzosiak, W.J. miRNAmotif—A Tool for the Prediction of Pre-miRNA–Protein Interactions. Int. J. Mol. Sci. 2018, 19, 4075. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19124075

AMA Style

Urbanek-Trzeciak MO, Jaworska E, Krzyzosiak WJ. miRNAmotif—A Tool for the Prediction of Pre-miRNA–Protein Interactions. International Journal of Molecular Sciences. 2018; 19(12):4075. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19124075

Chicago/Turabian Style

Urbanek-Trzeciak, Martyna O., Edyta Jaworska, and Wlodzimierz J. Krzyzosiak. 2018. "miRNAmotif—A Tool for the Prediction of Pre-miRNA–Protein Interactions" International Journal of Molecular Sciences 19, no. 12: 4075. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19124075

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop