PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of mmrLink to Publisher's site
 
Mol Med Rep. 2017 March; 15(3): 1071–1078.
Published online 2017 January 13. doi:  10.3892/mmr.2017.6117
PMCID: PMC5367367

IsomiR expression patterns in canonical and Dicer-independent microRNAs

Abstract

Multiple microRNA (miRNA) variants, known as isomiRs, are extensively distributed in miRNA loci and predominantly derive from the alternative cleavage of Drosha/Dicer and 3′addition events. The present study aimed to investigate the expression patterns of multiple isomiRs in typical miRNA and Dicer-independent miRNA loci by conducting evolutionary and expression analysis using public datasets. Although different miRNA maturation processes exist, multiple isomiRs can be detected by similar expression distributions. However, isomiR expression in Dicer-independent miRNA loci tends to be at a moderate level, particularly for random distribution in the ends that are split by Dicer in the typical miRNA loci. Compared with the mature miRNA locus (dominant miRNA locus), the non-dominant miRNA locus indicates an expression distribution similar to that of the Dicer-independent miRNA locus. These results increase the understanding of multiple isomiRs in the progression of diseases.

Keywords: microRNA, isomiR, expression, evolution

Introduction

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have been extensively studied as crucial negative regulatory molecules. Canonical miRNA is first processed from a primary miRNA (pri-miRNA) transcript via the Drosha-dependent microprocessor complex (1), and the process generates a precursor miRNA (pre-miRNA) with a stem-loop hairpin structure. Then, it is recognized and cut by Dicer into miRNA-miRNA* duplexes (2) and subsequently pre-miRNA is transported to the cytoplasm by exportin-5 (3). The Dicer protein is widely distributed in plants, metazoans and fungi (4,5), and it is also involved in the processing of other small RNA species, including small nuclear RNAs (6). Recurrent somatic mutations in Drosha can induce changes in miRNA expression (7), and somatic mutations in Drosha and Dicer1 can impair miRNA biogenesis (8). Simultaneously, alternative pathways of miRNA biogenesis have been reported (9), and some miRNAs are generated by a novel processing pathway independent of Dicer (1013). For example, mirtrons can be processed by splicing from precursor hairpins (14,15), and their precursor sequences are shorter than canonical pri-miRNAs because these comprise miRNA-miRNA* duplexes alone. Pre-miRNA hairpins of miRNAs from mirtron-like sources are generated via splicing of short introns and exosome-mediated trimming (16,17).

Typical miRNA biogenesis indicates that mature miRNAs can be incorporated into the RNA-induced silencing complex that then binds to the 3′-untranslated region of the target mRNA to degrade mRNAs or to repress translation (18), while another strand, termed miRNA star (miRNA*), is degraded to an inactive strand. However, accumulating evidence suggests that degraded miRNA* may also serve an important role in gene regulation at the post-transcriptional levels, as well as in the mature miRNA sequence (1923), and miRNA-miRNA* duplexes are also termed miRNA (miR)-#-5p-miR-#-3p duplexes. In the miRNA locus, a series of miRNA variants, called isomiRs, have been widely detected based on high-throughput sequencing datasets (2430). These multiple isomiRs are predominantly derived from the alternative and imprecise cleavage of Drosha and Dicer during pri-miRNA/pre-miRNA processing, and 3′ addition events in miRNA maturation processes (27). In the specific miRNA locus, isomiR profiles are always stable across different samples and different species (26,30), indicating a relatively stable miRNA maturation process.

Drosha and Dicer have important roles in the miRNA maturation process, and their imprecise and alternative cleavage processes largely contribute to the generation of multiple isomiRs. It is known that the isomiR profiles are always stable because of Drosha and Dicer, even across different animal species, however fewer studies have focused on isomiR expression profiles between typical miRNA loci and Dicer-independent miRNA loci. Generally, several dominant isomiRs (typically one to three) can be identified in the canonical miRNA locus, and others always possess lower expression rates (26,29,31). Therefore, according to characteristics of cleavage and expression patterns of multiple isomiRs, the current study attempted to explore the potential associations between isomiR profiles in different miRNA loci with diversity in the maturation processes, and simultaneously discuss evolutionary patterns of Drosha/Dicer. According to miRNA expression profiles, abundantly expressed miR-451 has been identified as an miRNA that is Dicer-independent, and based on previous studies focusing on isomiR expression (3234), relevant expression and evolutionary analysis were performed using public datasets. The present study may provide additional information regarding miRNA/isomiR biogenesis.

Materials and methods

Source data

Both Drosha and Dicer nucleotide and protein sequences of eight vertebrates, including Amphibia (Xenopus tropicalis), Aves (Gallus gallus; Taeniopygia guttata), Mammalia (Homo sapiens; Sus scrofa; Mus musculus), Pisces (Danio rerio) and Sauria (Anolis carolinensis) were collected from GenBank (Table I). Simultaneously, available expression data of miRNA/isomiR in breast cancer (BC) samples (n=683) and normal samples (n=87) were acquired from the Cancer Genome Atlas pilot project (http://cancergenome.nih.gov). Then, due to the fact that the selected BC samples were derived from female patients, samples from other human diseases from the Cancer Genome Atlas pilot project (http://cancergenome.nih.gov) were selected: Thyroid carcinoma (tumor, n=507; normal, n=59) and prostate adenocarcinoma (tumor, n=498; normal, n=52). The selection generated more information about isomiR expression with potential gender differences.

Table I.
Drosha and Dicer in eight vertebrates.

Sequence and expression analyses

The amino acid sequences of Drosha and Dicer were aligned using Clustal X 2.0 (http://www.clustal.org/clustal2) (35), and phylogenetic associations were reconstructed using the MEGA 5.10 software (http://www.megasoftware.net) (36) based on the neighbor-joining method. Genetic distance was simultaneously estimated using MEGA. Functional domains of Drosha and Dicer were identified using the Pfam database (http://pfam.xfam.org) (37), and inferred domains were predicted using SMART (http://smart.embl-heidelberg.de) (38).

According to biological characteristics of isomiRs, the relative expression rate of each isomiR was estimated based on all the isomiRs in the miRNA locus, and then the average expression rate was calculated across different individuals. Herein, due to the larger differences in sample sizes between tumor and normal samples, expression patterns were first estimated using different sample sizes, and then were estimated using equal sample sizes by randomly selecting 87 tumor samples in BC. Simultaneously, relevant analysis was further performed using other samples. According to miRNA/isomiR expression profiles and typical isomiR expression patterns, the simultaneously dominantly expressed miR-#-5p and miR-#-3p from hsa-miR-21, −30a and −30e were collected, and another miRNA that was Dicer-independent (miR-451) was selected for further analysis.

Statistical analysis

The expression levels of the isomiRs in multiple individuals, particularly those of dominant isomiRs, were expressed as the mean ± standard deviation. Differences in genetic distance were estimated using the t-test between Drosha and Dicer, and P<0.05 was considered to indicate a statistically significant difference. Simultaneously, the 95% confidence interval was also estimated. The degree of variation of relative expression of isomiRs was presented using a box plot. Relevant statistical description or analysis was conducted using Stata software, version 11 (StataCorp LP, College Station, TX, USA).

Results and discussion

The typical miRNA maturation process depends on Drosha and Dicer, and the miRNA-#-5p:miRNA-#-3p duplex is generated through cleavage of pri-miRNA and pre-miRNA (Fig. 1A). Although single miRNA sequences have been widely studied, the miRNA locus generates a series of isomiRs with diverse sequences and expression levels that are predominantly derived from the imprecise and alternative cleavage of Drosha/Dicer and 3′ addition events (2430). As a specific sequence in multiple isomiRs, the canonical miRNA sequence is not always the most dominantly expressed sequence (such as hsa-miR-30a-5p in the miRNA locus; Fig. 1) (29). Although sample sizes differed between tumor and normal samples, similar expression patterns were obtained using equally sized samples (Fig. 1). The isomiR expression profiles were always stable across different samples, and this feature was also observed in different cells, tissues and animal species, although abnormal isomiR expression patterns were also detected, particularly in certain diseased samples. The current study attempted to compare isomiR expression profiles between typical miRNA loci and miRNA loci with independent of Dicer using public sequencing datasets. Compared with recent studies on the expression and evolutionary patterns of isomiRs (29,3941), the present study examined the expression profile of a specific miRNA gene (hsa-miR-30a), that can generate two kinds of mature miRNAs, and a dominantly expressed miRNA independent of Dicer (hsa-miR-451), which has been updated as miR-451a, and was selected due to the fact that fewer miRNAs independent of Dicer are dominantly expressed.

Figure 1.
An example of multiple isomiRs and their expression patterns. Percentages (%) indicate the percentage of a certain isomiR compared with all detected isomiRs in the miRNA locus. (A) Hsa-mir-30a yields miR-30a-5p and miR-30a-3p with a series of multiple ...

As crucial factors that lead to multiple isomiRs in miRNA maturation, Drosha and Dicer are well-conserved and have similar evolutionary rates across different vertebrates (Fig. 2 and Table II). These similar evolutionary patterns may contribute to recognition and cleavage of Drosha and Dicer in the miRNA maturation process, and further contribute to stable isomiR expression profiles, particularly for the majority of isomiRs that are 3′ isomiRs with the same 5′ ends. In pre-miRNA with a stem-loop structure, dominantly expressed miRNA may be located in the 5p or 3p arm, however all the miRNA loci generate dominant 3′ isomiRs and rare 5′ isomiRs. Notably, more pre-miRNAs have been identified to generate mature and functional miR-#-5p and miR-#-3p. These observations suggest the relative precise cleavage of Drosha and Dicer on the 5′ ends during miRNA maturation process, although they serve important roles in different regions in cells. Divergence of 3′ ends among isomiRs from a given miRNA locus and among miRNAs in different species may be predominantly derived from the imprecise cleavage of Drosha and Dicer on 3′ ends, in addition to modification events following miRNA maturation.

Figure 2.
Phylogenetic trees of Drosha and Dicer. Similar phylogenetic associations are detected based on Drosha and Dicer sequences. hsa, Homo sapiens; ssc, Sus scrofa; mmu, Mus musculus; aca, Anolis carolinensis; gga, Gallus gallus; tgu, Taeniopygia guttata; ...
Table II.
Evolutionary distance between Drosha and Dicer.

However, compared with typical miRNAs that are generated by cleavage of Drosha and Dicer, certain miRNAs are generated by alternative pathways during miRNA biogenesis (9), including a novel processing pathway independent of Dicer (1012). For example, miR-451 is generated via the Dicer-independent pathway, and then contributes to the regulatory network as a functional miRNA (13). Although different maturation processes occur in miR-451 and canonical miRNA loci (including miR-21, 30a and 30e), multiple isomiR sequences are also present in the Dicer-independent miRNA locus (Figs. 35 and Table III). Similar to typical miRNA loci (26,29,42), several dominant isomiRs are dominantly expressed, and the majority of isomiRs are 3′ isomiRs. The highly conserved 5′ ends ensure the consistency of the function between different isomiR sequences, which in turn facilitate the co-regulation of target mRNAs. However, the most dominant isomiR only possesses moderate expression (approximately 45% of the total expression), and other dominant isomiRs have similar expression levels (Fig. 3A). The moderate expression patterns differ from the most typical miRNAs that are prone to have predominant isomiRs with absolute abundant expression (such as isomiRs from the miR-30a-5p locus).

Figure 3.
Multiple isomiRs in miR-451 and variations across different samples. Percentage (%) in the Y-axis indicates the percentage of a certain isomiR compared with all detected isomiRs in the miRNA locus. (A) hsa-miR-451 was identified to undergo Dicer-independent ...
Figure 5.
Variations in isomiRs across different samples using equal sample sizes (n=87). Percentage (%) in the Y-axis indicates the percentage of a specific isomiR compared with all detected isomiRs in the miRNA locus. miR, microRNA; BC, breast cancer.
Table III.
Dominantly expressed isomiR sequences in Figs. 2 and and33.

Another non-dominant miRNA strand such as the miR-30a-3p locus induces moderate isomiR expression profiles, similar to that observed in Dicer-independent miR-451 locus (40). Notably, although several dominant isomiRs are present in the passenger strand, distinct differences were detected between dominant isomiRs and other rare isomiRs, whereas in the miR-451 locus, no evident boundary was detected. The moderate level of expression may have resulted from the random cleavage of the hairpin, whereas isomiRs in most typical miRNAs are strictly controlled. The mir-451 hairpin could not be recognized by Dicer because of the insufficient amount of duplexes formed, and therefore it was directly cleaved by Ago2. The Dicer-independent maturation process may lead to relative random cleavage sites at the 3′ ends of miRNAs, which is different from the non-random 3′ ends of multiple isomiRs derived from typical miRNA locus. Generally, in typical miRNA loci, dominant isomiRs are evident, and there are always one to three types of abundantly expressed isomiRs that can be used to infer dominant cleavage sites in Drosha and Dicer. These observations indicate the cleavage bias of Drosha and Dicer, which further contributes to the distinct expression profiles of dominant isomiRs.

The current study also observed that although isomiR expression profiles are always stable across different samples, including between tumor and normal samples, and variations were also detected between diseased and normal samples (39). In diseased samples, a larger degree of variation was observed compared with that among normal samples, which implicates that the miRNA maturation process may be affected in the abnormal micro-environment (Figs. 3 and and5).5). The dispersed expression patterns of isomiRs may be derived from changes in Drosha and Dicer cleavage, modification or regulation following miRNA maturation in tumor samples. These distinct characteristics may furthermore provide additional information on the dynamic expression profiles of the coding-non-coding RNA regulatory network.

In conclusion, expression analysis of isomiRs demonstrated that Dicer-independent miRNA involves a moderate level of isomiR expression compared with that observed in typical miRNA loci. Notably, a similar expression profile can be detected between non-dominant miRNA and Dicer-independent miRNA loci, thereby suggesting more complex miRNA maturation processes, particularly at the isomiR levels. Dynamic miRNA/isomiR expression profiles further enrich the regulatory network, particularly those involving miRNA-miRNA and miRNA-mRNA interactions. Additional studies on multiple isomiRs, particularly their potential versatile biological roles, may prove beneficial.

Figure 4.
Expression patterns of multiple isomiRs in different miRNA loci and samples. Percentage (%) in the Y-axis indicates the percentage of a specific isomiR compared with all detected isomiRs in the miRNA locus and the X-axis indicates the isomiR types. A ...

Acknowledgements

The current study was supported by the National Natural Science Foundation of China (grant nos. 31271261, 31401009 and 61301251), the Program for the Top Young Talents by the Organization Department of the CPC Central Committee, the National Natural Science Foundation of Jiangsu (grant no. BK20130885), the Key Project of Social Development in Jiangsu Province (grant no. BE2016773), the Research Fund for the Doctoral Program of Higher Education of China (grant no. 20133234120009), Shandong Provincial Key Laboratory of Functional Macromolecular Biophysics, the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), sponsored by NUPTSF (grant no. NY215068).

References

1. Gregory RI, Yan KP, Amuthan G, Chendrimada T, Doratotaj B, Cooch N, Shiekhattar R. The Microprocessor complex mediates the genesis of microRNAs. Nature. 2004;432:235–240. doi: 10.1038/nature03120. [PubMed] [Cross Ref]
2. He L, Hannon GJ. MicroRNAs: Small RNAs with a big role in gene regulation. Nat Rev Genet. 2004;5:522–531. doi: 10.1038/nrg1415. [PubMed] [Cross Ref]
3. Lund E, Güttinger S, Calado A, Dahlberg JE, Kutay U. Nuclear export of microRNA precursors. Science. 2004;303:95–98. doi: 10.1126/science.1090599. [PubMed] [Cross Ref]
4. Cerutti H, Casas-Mollano JA. On the origin and functions of RNA-mediated silencing: From protists to man. Curr Genet. 2006;50:81–99. doi: 10.1007/s00294-006-0078-x. [PMC free article] [PubMed] [Cross Ref]
5. Shabalina SA, Koonin EV. Origins and evolution of eukaryotic RNA interference. Trends Ecol Evol. 2008;23:578–587. doi: 10.1016/j.tree.2008.06.005. [PMC free article] [PubMed] [Cross Ref]
6. Langenberger D, Bermudez-Santana CI, Stadler PF, Hoffmann S. Identification and classification of small RNAs in transcriptome sequence data. Pac Symp Biocomput. 2010:80–87. [PubMed]
7. Torrezan GT, Ferreira EN, Nakahata AM, Barros BD, Castro MT, Correa BR, Krepischi AC, Olivieri EH, Cunha IW, Tabori U, et al. Recurrent somatic mutation in DROSHA induces microRNA profile changes in Wilms tumour. Nat Commun. 2014;5:4039. doi: 10.1038/ncomms5039. [PMC free article] [PubMed] [Cross Ref]
8. Rakheja D, Chen KS, Liu Y, Shukla AA, Schmid V, Chang TC, Khokhar S, Wickiser JE, Karandikar NJ, Malter JS, et al. Somatic mutations in DROSHA and DICER1 impair microRNA biogenesis through distinct mechanisms in Wilms tumours. Nat Commun. 2014;2:4802. doi: 10.1038/ncomms5802. [PMC free article] [PubMed] [Cross Ref]
9. Yang JS, Lai EC. Alternative miRNA biogenesis pathways and the interpretation of core miRNA pathway mutants. Mol Cell. 2011;43:892–903. doi: 10.1016/j.molcel.2011.07.024. [PMC free article] [PubMed] [Cross Ref]
10. Cheloufi S, Dos Santos CO, Chong MM, Hannon GJ. A dicer-independent miRNA biogenesis pathway that requires Ago catalysis. Nature. 2010;465:584–589. doi: 10.1038/nature09092. [PMC free article] [PubMed] [Cross Ref]
11. Cifuentes D, Xue H, Taylor DW, Patnode H, Mishima Y, Cheloufi S, Ma E, Mane S, Hannon GJ, Lawson ND, et al. A novel miRNA processing pathway independent of Dicer requires Argonaute2 catalytic activity. Science. 2010;328:1694–1698. doi: 10.1126/science.1190809. [PMC free article] [PubMed] [Cross Ref]
12. Miyoshi K, Miyoshi T, Siomi H. Many ways to generate microRNA-like small RNAs: Non-canonical pathways for microRNA production. Mol Genet Genomics. 2010;284:95–103. doi: 10.1007/s00438-010-0556-1. [PubMed] [Cross Ref]
13. Yang JS, Maurin T, Robine N, Rasmussen KD, Jeffrey KL, Chandwani R, Papapetrou EP, Sadelain M, O'Carroll D, Lai EC. Conserved vertebrate mir-451 provides a platform for Dicer-independent, Ago2-mediated microRNA biogenesis. Proc Natl Acad Sci USA. 2010;107:15163–15168. doi: 10.1073/pnas.1006432107. [PubMed] [Cross Ref]
14. Okamura K, Hagen JW, Duan H, Tyler DM, Lai EC. The mirtron pathway generates microRNA-class regulatory RNAs in Drosophila. Cell. 2007;130:89–100. doi: 10.1016/j.cell.2007.06.028. [PMC free article] [PubMed] [Cross Ref]
15. Ruby JG, Jan CH, Bartel DP. Intronic microRNA precursors that bypass Drosha processing. Nature. 2007;448:83–86. doi: 10.1038/nature05983. [PMC free article] [PubMed] [Cross Ref]
16. Flynt AS, Greimann JC, Chung WJ, Lima CD, Lai EC. MicroRNA biogenesis via splicing and exosome-mediated trimming in Drosophila. Mol Cell. 2010;38:900–907. doi: 10.1016/j.molcel.2010.06.014. [PMC free article] [PubMed] [Cross Ref]
17. Chong MM, Zhang G, Cheloufi S, Neubert TA, Hannon GJ, Littman DR. Canonical and alternate functions of the microRNA biogenesis machinery. Genes Dev. 2010;24:1951–1960. doi: 10.1101/gad.1953310. [PubMed] [Cross Ref]
18. Kim VN, Han J, Siomi MC. Biogenesis of small RNAs in animals. Nat Rev Mol Cell Bio. 2009;10:126–139. doi: 10.1038/nrm2632. [PubMed] [Cross Ref]
19. Okamura K, Phillips MD, Tyler DM, Duan H, Chou YT, Lai EC. The regulatory activity of microRNA* species has substantial influence on microRNA and 3′ UTR evolution. Nat Struct Mol Biol. 2008;15:354–363. doi: 10.1038/nsmb.1409. [PMC free article] [PubMed] [Cross Ref]
20. Czech B, Zhou R, Erlich Y, Brennecke J, Binari R, Villalta C, Gordon A, Perrimon N, Hannon GJ. Hierarchical rules for Argonaute loading in Drosophila. Mol Cell. 2009;36:445–456. doi: 10.1016/j.molcel.2009.09.028. [PMC free article] [PubMed] [Cross Ref]
21. Okamura K, Liu N, Lai EC. Distinct mechanisms for microRNA strand selection by Drosophila Argonautes. Mol Cell. 2009;36:431–444. doi: 10.1016/j.molcel.2009.09.027. [PMC free article] [PubMed] [Cross Ref]
22. Guo L, Lu Z. The fate of miRNA* strand through evolutionary analysis: Implication for degradation as merely carrier strand or potential regulatory molecule? PLoS One. 2010;5:e11387. doi: 10.1371/journal.pone.0011387. [PMC free article] [PubMed] [Cross Ref]
23. Jagadeeswaran G, Zheng Y, Sumathipala N, Jiang H, Arrese EL, Soulages JL, Zhang W, Sunkar R. Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development. BMC Genomics. 2010;11:52. doi: 10.1186/1471-2164-11-52. [PMC free article] [PubMed] [Cross Ref]
24. Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, Pfeffer S, Rice A, Kamphorst AO, Landthaler M, et al. A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007;129:1401–1414. doi: 10.1016/j.cell.2007.04.040. [PMC free article] [PubMed] [Cross Ref]
25. Morin RD, Aksay G, Dolgosheina E, Ebhardt HA, Magrini V, Mardis ER, Sahinalp SC, Unrau PJ. Comparative analysis of the small RNA transcriptomes of Pinus contorta and Oryza sativa. Genome Res. 2008;18:571–584. doi: 10.1101/gr.6897308. [PubMed] [Cross Ref]
26. Guo L, Yang Q, Lu J, Li H, Ge Q, Gu W, Bai Y, Lu Z. A comprehensive survey of mirna repertoire and 3′ addition events in the placentas of patients with pre-eclampsia from high-throughput sequencing. PLoS One. 2011;6:e21072. doi: 10.1371/journal.pone.0021072. [PMC free article] [PubMed] [Cross Ref]
27. Neilsen CT, Goodall GJ, Bracken CP. IsomiRs-the overlooked repertoire in the dynamic microRNAome. Trends Genet. 2012;28:544–549. doi: 10.1016/j.tig.2012.07.005. [PubMed] [Cross Ref]
28. Lee LW, Zhang S, Etheridge A, Ma L, Martin D, Galas D, Wang K. Complexity of the microRNA repertoire revealed by next generation sequencing. RNA. 2010;16:2170–2180. doi: 10.1261/rna.2225110. [PubMed] [Cross Ref]
29. Guo L, Chen F. A Challenge for miRNA: Multiple isomiRs in miRNAomics. Gene. 2014;544:1–7. doi: 10.1016/j.gene.2014.04.039. [PubMed] [Cross Ref]
30. Burroughs AM, Ando Y, de Hoon MJ, Tomaru Y, Nishibu T, Ukekawa R, Funakoshi T, Kurokawa T, Suzuki H, Hayashizaki Y, Daub CO. A comprehensive survey of 3′ animal miRNA modification events and a possible role for 3′ adenylation in modulating miRNA targeting effectiveness. Genome Res. 2010;20:1398–1410. doi: 10.1101/gr.106054.110. [PubMed] [Cross Ref]
31. Guo L, Li H, Lu J, Yang Q, Ge Q, Gu W, Bai Y, Lu Z. Tracking miRNA precursor metabolic products and processing sites through completely analyzing high-throughput sequencing data. Mol Biol Rep. 2012;39:2031–2038. doi: 10.1007/s11033-011-0950-8. [PubMed] [Cross Ref]
32. Loher P, Londin ER, Rigoutsos I. Isomir expression profiles in human lymphoblastoid cell lines exhibit population and gender dependencies. Oncotarget. 2014;5:8790–8802. doi: 10.18632/oncotarget.2405. [PMC free article] [PubMed] [Cross Ref]
33. Llorens F, Hummel M, Pantano L, Pastor X, Vivancos A, Castillo E, Mattlin H, Ferrer A, Ingham M, Noguera M, et al. Microarray and deep sequencing cross-platform analysis of the mirRNome and isomiR variation in response to epidermal growth factor. BMC Genomics. 2013;14:371. doi: 10.1186/1471-2164-14-371. [PMC free article] [PubMed] [Cross Ref]
34. Guo L, Yu J, Liang T, Zou Q. miR-isomiRExp: A web-server for the analysis of expression of miRNA at the miRNA/isomiR levels. Sci Rep. 2016;6:23700. doi: 10.1038/srep23700. [PMC free article] [PubMed] [Cross Ref]
35. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al. Clustal W and clustal X version 2.0. Bioinformatics. 2007;23:2947–2948. doi: 10.1093/bioinformatics/btm404. [PubMed] [Cross Ref]
36. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–2729. doi: 10.1093/molbev/mst197. [PMC free article] [PubMed] [Cross Ref]
37. Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, et al. Pfam: The protein families database. Nucleic Acids Res. 2014;42:D222–D230. doi: 10.1093/nar/gkt1223. [PMC free article] [PubMed] [Cross Ref]
38. Letunic I, Doerks T, Bork P. SMART: Recent updates, new developments and status in 2015. Nucleic Acids Res. 2014;43:257–260. doi: 10.1093/nar/gku949. [PMC free article] [PubMed] [Cross Ref]
39. Guo L, Yu J, Yu H, Zhao Y, Chen S, Xu C, Chen F. Evolutionary and expression analysis of miR-#-5p and miR-#-3p at the miRNAs/isomiRs levels. Biomed Res Int. 2015;2015:168358. doi: 10.1155/2015/168358. [PMC free article] [PubMed] [Cross Ref]
40. Guo L, Zhang H, Zhao Y, Yang S, Chen F. Selected isomiR expression profiles via arm switching? Gene. 2014;533:149–155. doi: 10.1016/j.gene.2013.09.102. [PubMed] [Cross Ref]
41. Guo L, Zhao Y, Zhang H, Yang S, Chen F. Close association between paralogous multiple isomiRs and paralogous/orthologues miRNA sequences implicates dominant sequence selection across various animal species. Gene. 2013;527:624–629. doi: 10.1016/j.gene.2013.06.083. [PubMed] [Cross Ref]
42. Guo L, Chen F, Lu Z. Multiple IsomiRs and Diversity of miRNA Sequences Unveil Evolutionary Roles and Functional Relationships Across Animals. MicroRNA and Non-Coding RNA: Technology, Developments and Applications. 2013:127–144.

Articles from Molecular Medicine Reports are provided here courtesy of Spandidos Publications