Email updates

Keep up to date with the latest news and content from Silence and BioMed Central.

This article is part of the series Non-coding RNA.

Open Access Research

MicroRNAs from the same precursor have different targeting properties

Antonio Marco*, Jamie I MacPherson, Matthew Ronshaugen and Sam Griffiths-Jones*

Author Affiliations

Faculty of Life Sciences, Michael Smith Building, Oxford Road, University of Manchester, Manchester M13 9PT, UK

For all author emails, please log on.

Silence 2012, 3:8  doi:10.1186/1758-907X-3-8


The electronic version of this article is the complete one and can be found online at: http://www.silencejournal.com/content/3/1/8


Received:30 April 2012
Accepted:5 September 2012
Published:27 September 2012

© 2012 Marco et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The processing of a microRNA results in an intermediate duplex of two potential mature products that derive from the two arms (5′ and 3′) of the precursor hairpin. It is often suggested that one of the sequences is degraded and the other is incorporated into the RNA-induced silencing complex. However, both precursor arms may give rise to functional levels of mature microRNA and the dominant product may change from species to species, from tissue to tissue, or between developmental stages. Therefore, both arms of the precursor have the potential to produce functional mature microRNAs.

Results

We have investigated the relationship between predicted mRNA targets of mature sequences derived from the 5′ and 3′ arms of the same pre-microRNAs. Using six state-of-the-art target prediction algorithms, we find that 5′/3′ microRNA pairs target different sites in 3′ untranslated regions of mRNAs. We also find that these pairs do not generally target overlapping sets of genes, or functionally related genes.

Conclusions

We show that alternative mature products produced from the same precursor microRNAs have different targeting properties and therefore different biological functions. These data strongly suggest that developmental or evolutionary changes in arm choice will have significant functional consequences.

Keywords:
Arm switching; Gene regulation; miRNA; Target prediction

Background

MicroRNAs are crucial regulators of gene expression whose biogenesis is tightly controlled by multiple enzymes [1,2]. Primary microRNA transcripts are single-stranded RNA molecules that fold into hairpins, and are cleaved by two RNases producing an approximately 22-nucleotide RNA duplex [1]. In a process called arm-sorting or strand-sorting, one of the sequences of the duplex (derived from one of the arms of the precursor hairpin) associates with the RNA-induced silencing complex (RISC), which will mediate mRNA translational repression or transcript degradation (reviewed in [1]). The other arm (the star sequence or microRNA*) is generally considered a byproduct and is typically degraded [3]. However, it has become clear that both arms of the hairpin may produce functional mature products in many cases [4,5]. Indeed, microRNA* sequences are often highly expressed, evolutionarily conserved, and associated with RISC proteins [5]. For example, the microRNA* product of the Drosophila bantam gene is present at approximately 10-fold greater levels than any other microRNA product across a range of cell types and developmental stages. MicroRNA* sequences have been shown to be loaded into the RISC complex and to repress target genes both in vivo and in vitro[5-7].

Recent studies have shown that precursor microRNAs can change the arm from which the dominant functional mature microRNA is produced. This process, called arm-switching, occurs both in different tissues and developmental stages [4,8,9] and during evolution [10-13]. During microRNA biogenesis, both arms are produced at equal amounts in a given cell and, later on, one of the arms is usually degraded. For that reason, one may expect that microRNAs from the same precursor have similar targeting properties. Alternatively, since opposite arms of the hairpin have different sequences, it is expected that they target different sites. If these sites are in different transcripts, changes in arm usage would have the potential to alter microRNA function. So far, the only studied case is the mir-100/10 family, for which we have shown that opposite arms of precursor microRNAs do not significantly share target genes [13]. The functional consequences of changes in arm usage have not been extensively studied.

MicroRNA target recognition is mediated by complementary base-pairing between the microRNA and the 3′ untranslated regions (UTR) of targeted transcripts [14]. The number of experimentally validated microRNA/target pairs remains limited. However, computational prediction of microRNA targets has been widely used, although these approaches produce high rates of false positives [15]. In spite of this limitation, computational prediction of targets permits the study of general binding properties of a given microRNA. A widely accepted view of microRNA target preferences relies on nucleotides 2 to 7 of a microRNA, the so-called seed sequence, which recognizes binding sites often by perfect complementarity to the targeted transcripts (reviewed in [14]). However, distinct modes of target recognition have been described and they form the basis of distinct prediction algorithms. Since different prediction strategies are based on different assumptions and may give quite different results, it is often useful to apply a variety of algorithms to study the targeting properties of microRNAs.

Here, we use multiple target prediction algorithms to predict targets of human and fly microRNAs. We assess whether pairs of mature sequences derived from the 5′ and 3′ arms of the same precursor target identical sites (Figure  1A), different sites in the same gene transcripts (Figure  1B) and different genes in the same functional pathways (Figure  1C).

thumbnailFigure 1. Possible targeting properties of 5′/3′ microRNA pairs. (A) Both 5′ and 3′ products bind to the same target. (B) MicroRNA products bind to different sites in the same transcript. (C) MicroRNA products bind to different transcripts that act in the same functional pathway.

Results

Mature microRNAs from the same precursor have distinct target sites

We tested whether alternative mature microRNA products derived from the 5′ and 3′ arms of the same precursor share predicted target sites (Figure  1A). We predicted all canonical seed targets for all microRNAs in Drosophila melanogaster and human [14] and counted how many target sites have pairs of microRNAs from the same precursor in common. We observed that not a single predicted site was shared between the pairs of mature microRNAs from Drosophila. In humans, only one 5′/3′ microRNA pair, derived from mir-3648, had common targets, sharing 61 predicted sites out of a total of 569 and 455 sites predicted for the 5′ and 3′ microRNAs respectively. This is explained by the fact that both mature sequences are GC rich, and both seed sixmers are identical: GCCGCG. A closer inspection of the patterns of deep sequencing reads mapped to the mir-3648 locus (as shown in miRBase; [16]) suggests that mir-3648 may not be a bona fide microRNA, since it does not show a read pattern compatible with small RNA processing. In general, mature microRNAs from opposite arms have different sequences, therefore their propensity to target different sites is expected.

5′/3′ microRNA pairs target non-overlapping gene lists

UTRs may contain multiple target sites for different microRNAs. Therefore 5′/3′ pairs of microRNAs may target sites in the same transcript (Figure  1B). To test whether 5′/3′ microRNA pairs target common genes, we predicted regulated genes using six different and complementary methods: canonical seeds, miRanda, PITA, Diana-microT, RNAhybrid and TargetScan (with conservation - see Methods). For each 5′/3′ microRNA pair we compared the overlap between the predicted target lists and the expected overlap for random pairs of microRNAs (see Methods).

Canonical seeds, PITA, Diana-microT, RNAhybrid and TargetScan methods consistently showed that the overlap between genes targeted by 5′/3′ microRNA pairs is not statistically different from random expectation (Figure  2A). Only the miRanda algorithm suggests a significant overlap of genes targeted by 5′/3′ microRNA pairs (see below). The number of microRNA pairs with overlapping target gene predictions in the human dataset is about twice that for Drosophila. This is likely due to the fact that human 3′ UTRs are longer than those from Drosophila, and therefore the number of microRNAs predicted to target each transcript is significantly larger. Therefore, we performed a second analysis in the human set using a more stringent set of parameters (see Methods). These strict predictions yielded smaller overlapping values, but the overall findings remain robust to the parameter changes: only the miRanda set showed significant differences between the observed and the expected overlap values (Figure  2A).

thumbnailFigure 2. Distribution of target overlaps between 5′/3′ microRNA pairs. Average overlap of target predictions for multiple algorithms (grey boxes) and expected overlap based on random sampling (white boxes). Error bars depict standard errors of the means. Statistical differences between distributions underlying the plotted data were assessed by one-tailed Kolmogorov-Smirnov tests with Bonferroni correction. Only P-values below 0.05 are shown. (A) Overlap between lists of targeted genes from microRNA pairs from the same precursor. The human datasets include additional values for strict sets of predictions for four of the algorithms. (B) Term overlap between the targets of 5′/3′ pairs of microRNAs. Human term overlap values were calculated, when available, for the strict target prediction sets.

We investigated whether the observed overlap for miRanda predictions of gene targets of 5′/3′ microRNAs pairs could be explained by sequence composition biases. In particular, programs that use hybrid stability to detect microRNA targets (such as miRanda) may be biased by variable GC content [17]. We therefore studied the potential effect of composition bias on predicted microRNA targets in humans. We find that the number of predicted gene targets is highly correlated with the GC content of the microRNA (R2 = 0.72, P <0.001). There is also a positive correlation between the microRNA duplex GC content and the overlap between the genes targeted (R2 = 0.58, P <0.001). After removing those microRNAs with high GC content (defined as greater than 67% as in [17]), the overlap between target genes of human 5′/3′ microRNA pairs was still significant (P = 0.003). The overlap between miRanda predictions for 5′/3′ microRNAs pairs is therefore robust to sequence bias.

5′/3′ microRNA pairs do not target genes in the same functional classes

Different genes targeted by different microRNAs may have related functions or be involved in related pathways (Figure  1C). The functional similarity of two genes can be quantified by assessing the similarity of their annotation, for example using Gene Ontology (GO) terms [18,19]. This class of methods is known as semantic similarity measures. Semantic similarity using GO term annotation has been widely applied in genomics to compare functional similarity between pairs of genes (for example, [19,20]). Here we use a measure called average term overlap (TO) to estimate the functional similarity between lists of genes (see Methods). Values for average TO were calculated for the lists of genes targeted by 5′/3′ pairs of microRNAs. We did not observe any significant overlap in the functions of genes targeted by 5′/3′ pairs of microRNAs based on GO annotations with any of the algorithms. A slight bias (although not significant) for 5′/3′ microRNAs to target genes with related functions using miRanda (Figure  2B) is explained by the significant overlap of targeted genes discussed above (Figure  2A). From these analyses, we conclude that alternative microRNAs from the same precursor have significantly different targeting properties.

Cases in which 5′/3′ pairs have similar targets

We have shown that miRanda predictions suggest that some 5′/3′ microRNA pairs tend to target common genes. We explored whether the relative amount of microRNA produced from each arm of the hairpin precursor is associated with the targeting properties for the human dataset. In Figure  3 we plot the average gene overlap for different levels of arm usage bias. Arm usage bias reflects the number of reads from deep sequencing experiments that map to one arm with respect to the other (see Methods), and was calculated only for microRNAs that have reads associated with both arms. The impact of arm usage bias in the targeting properties of human microRNAs is shown in Table  1. Where pairs of alternate microRNAs from the same hairpin are produced at ratios of at least 10:1 (that is, a mature product from one arm dominates), we find that the 5′/3′ pairs of microRNAs do not bind to overlapping lists of genes. MicroRNAs with low or no arm usage bias produce pairs of mature sequences that do bind to overlapping lists of genes (Table  1). By contrast, mature 5′/3′ microRNA pairs that are expressed at similar levels tend to bind more similar lists of genes (Figure  3). We observe a similar pattern in Drosophila: the subset of microRNAs producing mature sequences approximately equally from both arms share more targets than expected by chance (Table  1), although the differences are not statistically significant. The set of human microRNAs that produce similar amounts of mature products from each arm (ratio less than 3:1, and a minimum of 10 reads mapping to either arm) is shown in Table  2. Three out of the 11 pairs have a target overlap above the expected value (>0.071). We therefore show that the significant overlap of predicted gene targets of 5′/3′ microRNA pairs can be attributed to microRNAs that produce approximately equal amounts of mature sequences from both arms.

thumbnailFigure 3. Effect of microRNA arm usage on targeting properties. Sliding-window plot showing (A) mean target list overlap and (B) mean term (function) overlap with respect to the microRNA arm usage bias (see Methods).

Table 1. Effect of arm usage bias on gene overlap of miRanda predictions of 5/3microRNA pairs

Table 2. Human microRNAs with low arm usage bias

Discussion

In this work, we have shown that, in general, 5′/3′ mature microRNA sequences derived from the same microRNA precursor target non-overlapping lists of genes. The only exceptions derive from predictions made with the miRanda algorithm [21] of targets of mature sequences produced in equal concentrations from both arms of the precursor. miRanda takes into account hybrid stability of the target and the microRNA, as well as strong sequence complementarity in the seed region [21]. We envisage two possible explanations for the different result from miRanda predictions. On the one hand, the relaxation of the requirement for perfect complementarity in the seed region may allow miRanda to detect targets and trends that escape other prediction algorithms (probably at the expense of prediction specificity). Indeed, a small number of cases of 5′/3′ microRNA pairs binding to the same transcript have been described (for example, [22]). On the other hand, miRanda predictions may be susceptible to unknown biases such that the observed pattern is an artifact of the algorithm (although we rule out the effects of GC bias here). Nevertheless, all six different algorithms with two different sets of parameters, covering the spectrum of most existing target prediction algorithms [23], concur that 5′/3′ mature microRNA pairs do not target the same genes or pathways when the precursor produces functional products primarily from one of the arms.

Early experiments suggested that the thermodynamic properties of the microRNA duplex determine the sequence that is incorporated into the RISC, and hence, which arm is functional [5,6,24]. However, we recently proved that identical duplex sequences in Drosophila melanogaster and the beetle Tribolium castaneum can produce functional microRNAs from opposite arms [13]. Moreover, the dominant arm can change within the same species in different developmental stages or tissues [9-11]. This suggests that arm sorting can be determined by signals outside the mature microRNA duplex. Thus, changes in arm usage may occur without changing the nucleotide sequences of mature microRNAs, such that the potential targeting properties of both arms are unchanged (see also [5,6]). We have described five instances of arm switching between Drosophila and Tribolium microRNAs [12]: mir-10, mir-33, mir-275, mir-929 and mir-993. These microRNAs are highly expressed and, in each case, mature sequences are produced in ratios of around 10:1 [25]. In this work, we provide evidence that the targeting properties of 5′/3′ microRNA products are not similar when one mature product dominates. Therefore, arm-switching events in these five microRNAs [12] are predicted to lead to functional changes, as we previously suggested for mir-10 in Drosophila and Tribolium[13].

Conclusions

Alternative mature products from the same precursor microRNA have different targeting properties. Exceptions to this rule are observed for microRNAs from which both arms produce significant amounts of mature products using miRanda gene predictions. We therefore strongly suggest that microRNA arm preferences have important functional consequences. Comparative analysis of regulatory networks accounting for microRNA arm usage will be slightly more complex, yet biologically more meaningful.

Methods

We extracted all fly (D. melanogaster) and human (Homo sapiens) microRNAs from miRBase (version 16; [16]). This version of miRBase does not index 5′ and 3′ mature sequences for all microRNAs. Where a single mature sequence from a microRNA precursor is reported, we selected as the miR* sequence the most abundant read from the appropriate arm from high-throughput sequencing data displayed in miRBase (December 2010; [16]) and discarded sequences with no evidence for a miR* sequence. This resulted in a total of 163 and 426 pre-microRNAs in fly and human respectively. The expression datasets used in this analysis are listed in Additional file 1: Table S1.

Additional file 1. Table S1. Gene expression datasets.

Format: DOC Size: 41KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

We used six different algorithms to detect potential targets of mature microRNA sequences: canonical seeds as described in [14]; miRanda [21], a method based on hybrid energy and stability; PITA [26], which takes into account the site accessibility at 3′ UTRs; Diana-microT [27], a predictor that combines multiple features; RNAhybrid [28], which detects stable RNA-RNA duplexes; and TargetScan [29,30], a canonical seed detection program that also takes into account conservation of microRNAs and target sites. We ran TargetScan to identify target sites conserved in at least two species in the 3′ UTR alignments available from their webpage [30]. We generated target prediction datasets for each algorithm using default parameters. We also generated a second prediction set for human microRNAs (called the strict set) using each algorithm with the following parameter modifications: at least two sites in canonical seed predictions; miRanda targets with a score above 1,000, to reduce the number of targets to a tenth of the original predictions; PITA-predicted targets of a size of 7 to 8, with no mismatches or wobble positions; Diana-microT predictions with an MRE score above 0.6 as suggested by the authors.

We used as potential targets the largest 3′UTR available for each gene in Drosophila in Flybase (genome version BDGP 5.25 [31]) and in human from ENSEMBL (assembly 60 [32]). For each pair of mature products from a precursor microRNA, we identified potential targets with all six methods, and we calculated for each method the overlap between the lists of target sites as the number of commonly targeted sites divided by the total number of sites targeted by both arms (Jaccard similarity; [33]). Similarly, the overlap between lists of target genes was calculated as the number of commonly targeted genes divided by the total number of genes targeted by both arms. The expected distributions of values were calculated by selecting 10,000 random pairs of microRNA arms and calculating the target overlap for each pair.

To assess whether two lists of genes have a similar functional annotation, we cross-compared all gene pairs between the two lists and calculated semantic similarity using the term overlap (TO) measure [19,34] for the ‘biological process’ domain of Gene Ontology [35]. Average TO values for pairs of gene lists are defined as:

T O = i = 1 n j = 1 m T G i , G j n m

where T{Gi,Gj} is the number of common GO terms to which genes Gi and Gj are annotated. The TO analysis in humans was performed only for the strict target prediction sets. Expected average TO values were calculated by generating 1,000 randomized pairs.

Arm usage is defined as the relative production of mature products from one arm with respect to the other arm, and it is calculated as described in [12]. Only microRNAs with reads in both arms were included. An arm usage of 0 means that both arms produce the same amount of product. Each unit above 0 indicates a two-fold increase in the biased production of one of the arms.

Abbreviations

GO: gene ontology; RISC: RNA-induced silencing complex; TO: term overlap; UTR: untranslated region.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AM and SGJ conceived the project. AM and JIM performed the analyses. AM, MR and SGJ interpreted the results and wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We are grateful to Artemis Hatzigeorgiou for constructive discussion, and to AH and Manolis Maragkakis for providing Diana-microT predictions for both arms of each microRNA. We also thank Marc Rehmsmeier for help with the RNAhybrid program. This work was supported by the Biotechnology and Biological Sciences Research Council (BB/G011346/1 and BB/H017801/1), the Wellcome Trust (097820/Z/11/Z), and the University of Manchester (fellowships to SGJ and MR). JIM was funded by the BBSRC and Pfizer Global Research and Development.

References

  1. Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function.

    Cell 2004, 116:281-297. PubMed Abstract | Publisher Full Text OpenURL

  2. Krol J, Loedige I, Filipowicz W: The widespread regulation of microRNA biogenesis, function and decay.

    Nat Rev Genet 2010, 11:597-610. PubMed Abstract | Publisher Full Text OpenURL

  3. Schwarz DS, Hutvágner G, Du T, Xu Z, Aronin N, Zamore PD: Asymmetry in the assembly of the RNAi enzyme complex.

    Cell 2003, 115:199-208. PubMed Abstract | Publisher Full Text OpenURL

  4. Glazov EA, Cottee PA, Barris WC, Moore RJ, Dalrymple BP, Tizard ML: A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach.

    Genome Res 2008, 18:957-964. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Okamura K, Phillips MD, Tyler DM, Duan H, Chou Y, 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Yang JS, Phillips MD, Betel D, Mu P, Ventura A, Siepel AC, Chen KC, Lai EC: Widespread regulatory activity of vertebrate microRNA* species.

    RNA 2011, 17:312-326. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  9. Chiang HR, Schoenfeld LW, Ruby JG, Auyeung VC, Spies N, Baek D, Johnston WK, Russ C, Luo S, Babiarz JE, Blelloch R, Schroth GP, Nusbaum C, Bartel DP: Mammalian microRNAs: experimental evaluation of novel and previously annotated genes.

    Genes Dev 2010, 24:992-1009. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Ruby JG, Stark A, Johnston WK, Kellis M, Bartel DP, Lai EC: Evolution, biogenesis, expression, and target predictions of a substantially expanded set of Drosophila microRNAs.

    Genome Res 2007, 17:1850-1864. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. de Wit E, Linsen SEV, Cuppen E, Berezikov E: Repertoire and evolution of miRNA genes in four divergent nematode species.

    Genome Res 2009, 19:2064-2074. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Marco A, Hui JHL, Ronshaugen M, Griffiths-Jones S: Functional shifts in insect microRNA evolution.

    Genome Biol Evol 2010, 2:686-696. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Griffiths-Jones S, Hui JHL, Marco A, Ronshaugen M: MicroRNA evolution by arm switching.

    EMBO Rep 2011, 12:172-177. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Bartel DP: MicroRNAs: target recognition and regulatory functions.

    Cell 2009, 136:215-233. PubMed Abstract | Publisher Full Text OpenURL

  15. Alexiou P, Maragkakis M, Papadopoulos GL, Reczko M, Hatzigeorgiou AG: Lost in translation: an assessment and perspective for computational microRNA target identification.

    Bioinformatics 2009, 25:3049-3055. PubMed Abstract | Publisher Full Text OpenURL

  16. Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data.

    Nucleic Acids Res 2011, 39:D152-D157. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Davis N, Biddlecom N, Hecht D, Fogel GB: On the relationship between GC content and the number of predicted microRNA binding sites by MicroInspector.

    Comput Biol Chem 2008, 32:222-226. PubMed Abstract | Publisher Full Text OpenURL

  18. Lord PW, Stevens RD, Brass A, Goble CA: Semantic similarity measures as tools for exploring the gene ontology.

    Pacific Symp Biocomput 2003, 8:601-612. OpenURL

  19. Mistry M, Pavlidis P: Gene Ontology term overlap as a measure of gene functional similarity.

    BMC Bioinformatics 2008, 9:327. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  20. Marco A, Konikoff C, Karr TL, Kumar S: Relationship between gene co-expression and sharing of transcription factor binding sites in Drosophila melanogaster.

    Bioinformatics 2009, 25:2473-2477. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Enright A, John B, Gaul U, Tuschl T, Sander C, Marks DS: MicroRNA targets in Drosophila.

    Genome Biol 2003, 5:R1. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  22. Davis E, Caiment F, Tordoir X, Cavaillé J, Ferguson-Smith A, Cockett N, Georges M, Charlier C: RNAi-mediated allelic trans-interaction at the imprinted Rtl1/Peg11 locus.

    Curr Biol 2005, 15:743-749. PubMed Abstract | Publisher Full Text OpenURL

  23. Thomas M, Lieberman J, Lal A: Desperately seeking microRNA targets.

    Nat Struct Mol Biol 2010, 17:1169-1174. PubMed Abstract | Publisher Full Text OpenURL

  24. Ghildiyal M, Xu J, Seitz H, Weng Z, Zamore PD: Sorting of Drosophila small silencing RNAs partitions microRNA* strands into the RNA interference pathway.

    RNA 2010, 16:43-56. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Chung W-J, Okamura K, Martin R, Lai EC: Endogenous RNA interference provides a somatic defense against Drosophila transposons.

    Curr Biol 2008, 18:795-802. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E: The role of site accessibility in microRNA target recognition.

    Nat Genet 2007, 39:1278-1284. PubMed Abstract | Publisher Full Text OpenURL

  27. Maragkakis M, Alexiou P, Papadopoulos GL, Reczko M, Dalamagas T, Giannopoulos G, Goumas G, Koukis E, Kourtis K, Simossis VA, Sethupathy P, Vergoulis T, Koziris N, Sellis T, Tsanakas P, Hatzigeorgiou AG: Accurate microRNA target prediction correlates with protein repression levels.

    BMC Bioinformatics 2009, 10:295. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  28. Kruger J, Rehmsmeier M: RNAhybrid: microRNA target prediction easy, fast and flexible.

    Nucl Acids Res 2006, 34:W451-W454. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets.

    Cell 2005, 120:15-20. PubMed Abstract | Publisher Full Text OpenURL

  30. TargetScan 6.2.

    http://www.targetscan.org webcite

    OpenURL

  31. Flybase Homepage.

    http://flybase.org/ webcite

    OpenURL

  32. ENSEMBL Homepage.

    http://www.ensembl.org/ webcite

    OpenURL

  33. Gordon AD: Classification. 2nd edition. Chapman and Hall/CRC Press, Boca Ratón; 1999. OpenURL

  34. Lee HK, Hsu AK, Sajdak J, Qin J, Pavlidis P: Coexpression analysis of human genes across many microarray data sets.

    Genome Res 2004, 14:1085-1094. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology.

    Nat Genet 2000, 25:25-29. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL