Skip Navigation

Human Molecular Genetics 2005 14(Review Issue 2):R183-R190; doi:10.1093/hmg/ddi271
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Berezikov, E.
Right arrow Articles by Plasterk, R. H.A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Berezikov, E.
Right arrow Articles by Plasterk, R. H.A.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2005. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oxfordjournals.org

Camels and zebrafish, viruses and cancer: a microRNA update

Eugene Berezikov* and Ronald H.A. Plasterk

Hubrecht Laboratory, Uppsalalaan 8, 3584 CT Utrecht, The Netherlands

* To whom correspondence should be addressed. Tel: +31 302121828; Fax: +31 302516464; Email: berezikov{at}niob.knaw.nl

Received June 9, 2005; Accepted July 14, 2005


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 miRNA MACHINERY
 IDENTIFICATION OF miRNA GENES
 miRNA FUNCTIONS
 miRNAS AND CANCER
 VIRUS CONNECTIONS
 CONCLUSIONS
 REFERENCES
 
MicroRNAs (miRNAs) form an extensive class of RNA molecules that regulate gene expression at post-transcriptional level. In recent years, much progress has been made in dissection of biogenesis and functions of miRNAs. There are at least several hundred miRNA genes in the human genome, and the emerging evidence suggests that miRNAs are broadly implicated in gene regulation. Here, we review some recent advances, and particularly we discuss how comparative genomics helps to identify novel miRNA genes, how studies in zebrafish reveal roles of miRNAs in morphogenesis, how changes in miRNA expression patterns are connected with cancer and how host–virus coevolution exploits miRNA regulatory pathways.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 miRNA MACHINERY
 IDENTIFICATION OF miRNA GENES
 miRNA FUNCTIONS
 miRNAS AND CANCER
 VIRUS CONNECTIONS
 CONCLUSIONS
 REFERENCES
 
Non-coding RNAs (ncRNAs) are RNA molecules that do not encode proteins but have some other functions. The examples of well- and long-known ncRNAs include ribosomal, spliceosomal and transfer RNAs that perform essential housekeeping functions in a cell. Yet, recognition that other ncRNAs, such as introns and non-protein-coding transcripts, may also have crucial cellular functions, has emerged only during the last decade (1Go–3Go). Although only ~1.5% of the human genome encodes proteins, a substantially larger fraction of the genome is transcribed. The extent of genome transcription is best demonstrated in the recent work of Gingeras and colleagues, who used high-density tiling arrays to interrogate transcription of about one-third of the human genome at the 5 bp resolution level, and found that >15% of the examined regions were transcribed, which is an order of magnitude greater than expected from annotated exons and gene predictions (4Go). Interestingly, ~30% of the transcripts originate from unannotated genomic regions and thus are called transcripts of unknown function (TUFs). Two-thirds of these TUFs have a coding capacity of less than 100 amino acids and can very well be ncRNAs.

Information on the diversity of known ncRNAs is available from several databases (5Go–7Go). For example, RNAdb provides a listing of over 800 unique mammalian ncRNAs, including microRNAs (miRNAs), small nucleolar RNAs (snRNAs) and ncRNAs of unknown function but known to be developmentally regulated, disease-associated, imprinted, expressed pseudogenes or antisense transcripts (7Go). Among the variety of ncRNA classes, miRNAs have attracted most attention in recent years due to their abundance, unexpectedly broad involvement in gene regulation and intersection with RNA interference (RNAi) pathway. In this article, we review some recent advances in the miRNA field: from the many angles from which to view the field we have certainly chosen that of our own laboratory and we apologize for omissions or emphases. For information on snRNAs and other ncRNAs, we refer readers to other reviews (8Go–11Go).


    miRNA MACHINERY
 TOP
 ABSTRACT
 INTRODUCTION
 miRNA MACHINERY
 IDENTIFICATION OF miRNA GENES
 miRNA FUNCTIONS
 miRNAS AND CANCER
 VIRUS CONNECTIONS
 CONCLUSIONS
 REFERENCES
 
miRNAs form an extensive class of ncRNAs that regulate the expression of genes at post-transcriptional level. An unusual mechanism of gene regulation, in which a small 22 nt RNA molecule (lin-4) forms a duplex with the 3'UTR of a target mRNA (lin-14) and blocks its translation, was first discovered in Caenorhabditis elegans more than a decade ago (12Go,13Go). As lin-4/lin-14 sequences were conserved only in nematodes, the generality of this regulatory mechanism was initially not recognized, until 7 years later another C. elegans RNA, let-7, was shown to regulate lin-41 gene in the same manner (14Go,15Go). In contrast to lin-4/lin-14 pair, the sequence and expression pattern of let-7 RNA, as well as lin-41 target sites, appeared to be conserved in a wide variety of animals (16Go,17Go), suggesting an ancient origin of this type of gene regulation and the existence of a larger class of regulatory RNAs of this kind. Indeed, efforts of different laboratories for cloning small RNAs soon resulted in identification of more than a hundred additional genes that, like lin-4 and let-7, expressed ~22 nt RNAs potentially processed from stem–loop precursors (18Go–20Go). The newly discovered RNAs were collectively named miRNAs.

Understanding of miRNA biogenesis was greatly facilitated by advances in the field of RNAi, because it appeared that miRNA and siRNA pathways partially overlap (21Go–26Go). miRNA genes are transcribed by RNA polymerase II as large primary miRNA (pri-miRNA) transcripts that have CAP structures and poly(A) tails (27Go,28Go). Some miRNAs are arranged in clusters and transcribed as polycistrons (29Go), whereas ~40% of miRNAs reside in introns of protein-coding genes (30Go) and are presumably co-transcribed with host genes. Pri-miRNA transcripts form characteristic fold-back structures that are recognized and processed into ~70 nt imperfect stem–loop miRNA precursors (pre-miRNA) by a so-called microprocessor complex (31Go,32Go). This complex includes RNAse III enzyme Drosha (33Go) and its cofactor, the DiGregory syndrome critical region gene 8 (DGCR8) protein, also known as Pasha in D. melanogaser and in C. elegans (31Go,32Go,34Go,35Go). It is believed that Drosha needs Pasha/DGCR8 cofactor for recognition of pri-miRNA tertiary structure, which seems to be a primary determinant for substrate specificity (33Go,36Go,37Go).

After processing by the microprocessor complex, pre-miRNAs are exported by Ran-GTP dependent transporter exportin-5 (38Go–40Go) from the nucleus into cytoplasm, where they are processed by another RNAse III enzyme, Dicer, into imperfect dsRNA duplexes that contain both mature and complementary miRNA strands (41Go–43Go). The thermodynamic energy of the 5' ends of a miRNA–miRNA* duplex is usually different and the strand with lower 5' end stability is preferentially loaded into RNA-induced silencing complex (RISC), thus becoming a mature miRNA, and the complementary miRNA* strand is rapidly degraded (44Go,45Go). RISCs are ribonucleoprotein complexes that exist in different subtypes that probably reflect different functions of RISCs. The core components of RISC are proteins from the extensive Argonaute family (46Go–48Go). For example, the human miRNA-containing RISC contains Argonaute protein eIF2C2 and helicases Gemin3 and Gemin4 (49Go). Additional identified proteins that associate with RISC include nuclease Tudor-SN (50Go) and RNA-binding VIG and fragile-X related protein (51Go,52Go). Despite the extensive characterization of RISC complexes, exact biochemical mechanisms of RISC function remain largely unknown. Good progress in this direction was achieved with recent identification of Argonaute2 as a catalytic center of target mRNA cleavage by siRNA-containing RISC (53Go). As an alternative route to degradation, RISC can impose translational block of a target mRNA, and it is thought that the mode of RISC action is largely determined by the extent of complementarity between a effector RNA (siRNA or miRNA) and a target mRNA, with a perfect identity leading to mRNA cleavage and imperfect matching to translational block (47Go,54Go–56Go). Whereas in plants most miRNAs target mRNA to degradation (57Go–60Go), in animals all known miRNA–mRNA interactions lead to translational block (61Go), and only one case of miRNA-mediated cleavage of target mRNA has been reported so far (62Go). Recently, Lim et al. (63Go) used microarrays to analyze the effects of transfection of human cells with particular miRNAs and found that levels of many direct mRNA targets decreased, thus suggesting that influence of miRNAs on the levels of target mRNAs can be widespread. A more detailed recent reviews on miRNA biogenesis are described earlier (26Go,64Go).


    IDENTIFICATION OF miRNA GENES
 TOP
 ABSTRACT
 INTRODUCTION
 miRNA MACHINERY
 IDENTIFICATION OF miRNA GENES
 miRNA FUNCTIONS
 miRNAS AND CANCER
 VIRUS CONNECTIONS
 CONCLUSIONS
 REFERENCES
 
Information about all known miRNAs is accumulated in the miRNA registry—a database of published miRNA sequences (65Go), the recent release of which (6.0) contains 227 human, 230 mouse and 191 rat miRNA genes. Essentially, there are two approaches for the identification of novel miRNA genes: cloning of size-fractionated (18–25 nt) RNAs and computational prediction based on different structural features of miRNAs. The cloning approach was particularly efficient at the initial stages of miRNA hunting and allowed identification of more than two hundred miRNAs from human, mouse, Drosophila and C. elegans (18Go–20Go,49Go,66Go–73Go). Although powerful, the cloning approach has limitations and depends on the abundance and spatio-temporal expression patterns of miRNAs for their detection. The large number of miRNAs identified by cloning, however, formed a solid data set essential for designing efficient prediction algorithms, and several groups have developed such computational algorithms for identification of miRNAs in different organisms (67Go,74Go–80Go). All these algorithms follow one common theme (81Go) and search for fold-back structures that have characteristics similar to those of known miRNAs and are evolutionarily conserved. On the basis of the correlations between the number of computationally predicted miRNA candidates, the fraction of known miRNAs included in the predictions, frequencies of cloning of known miRNAs and verification rates for novel candidates, Bartel and colleagues (76Go) estimated that the human genome contains at most 255 miRNA genes. Similarly, it was estimated that there are about 120 and 110 miRNAs in C. elegans and in D. melanogaster genomes, respectively (75Go,77Go). An alternative estimate of 140–300 miRNA genes in C. elegans genome is provided by Ruvkun and colleagues (74Go).

We have recently revisited the estimate of a possible number of miRNA genes in the human genome (80Go). Our original intention was to investigate conservation patterns of different miRNA genes and their flanking sequences, and to that end, we sequenced genomic regions encompassing 122 miRNA genes in 10 primate species. This approach, known as phylogenetic shadowing (82Go), provides fine-grain resolution of sequence conservation and allows identification of functional elements that would have been missed in comparisons of evolutionary distant species. The main goal of our analysis was to probe for the presence of conserved elements in sequences immediately flanking miRNA genes. For example, such conserved elements, if found, could play a role in miRNA biogenesis. The obtained data suggested, however, that there are no such conserved elements common to different miRNAs. Instead, we found that immediate flanks of most miRNAs analyzed are not conserved, which lead to the realization that flanks variability may be a good filter for improving existing algorithms of miRNA gene prediction. In its most prominent form, miRNA conservation profiles resembled camel's humps (Fig. 1), and therefore, the typical profile received a code-name ‘camel’. We next developed a computational pipeline for identification of regions in human/rodent whole-genome alignments (WGA) that have ‘camel’-like conservation profiles and are able to form fold-back structures. Besides structural and conservation/variation criteria, our computational approach for identification of miRNA genes included an additional filter—the recently discovered property of miRNAs to have a lower folding free energy than random sequences of the same nucleotide composition (83Go). Analysis of human/mouse and human/rat WGAs identified a total of 978 candidate miRNA regions that satisfied the earlier mentioned criteria and included 158 known miRNAs (>80% of known miRNAs present in initial alignments). Searches in dog, cow, opossum, chicken and zebrafish genomes revealed that 678 of the predicted miRNAs could be also identified in at least one more species besides rodents. In addition, we confirmed a number of novel miRNA candidates by northern blot analysis. The conservative interpretation of our results suggested that there may be between 200 and 300 novel miRNAs in the candidate data set, bringing the total number of miRNA genes in the human genome to around 500. We also speculated that the actual number of miRNAs may be even higher, up to 1000, if highly restricted temporal and spatial expression patterns as well as taxon-specificity of some miRNAs are taken into account. Microarray and cloning experiments are currently underway to interrogate the complete set of predictions.



View larger version (21K):
[in this window]
[in a new window]
 
Figure 1. Camel-shaped conservation profile is characteristic for miRNA genes. An average conservation level (middle line) with 95% confidence intervals (upper and lower lines) was calculated for 64 miRNAs. Corresponding regions of the precursor miRNA (top) are shown by arrows. The figure is modified from Berezikov et al. (80Go).

 
A number of candidates from our data set are already proved to be real miRNAs by independent works. Margalit and colleagues (84Go) performed analysis of conservation and clustering patterns of known human miRNAs and used conservation profiles similar to ‘camel’, to predict a number of novel miRNA genes in the vicinity of known miRNAs. They then experimentally confirmed the expression of 14 predicted miRNAs, of which nine were also on our candidates list.

A higher estimate of the total number of miRNA genes in the human genome is also supported by the recent original work of Xie et al. (85Go), who analyzed conservation patterns of promoter and 3' UTR sequences in human, mouse, rat and dog genomes, and found a number of over-represented conserved regulatory motifs. Unexpectedly, motifs found in 3' UTRs showed a length distribution with a prominent peak at an eight base length (72 of 106 conserved motifs), which led authors to speculate that these motifs might be miRNA target sites. In addition, indeed, about half of known miRNAs matched through Watson-Crick pairing to the highly conserved 8mer motifs. The discovered 8mers were next used to search for novel miRNA genes and 242 conserved hairpin sequences were identified that included 113 known miRNAs and 129 novel predictions. Experimental verification confirmed the expression of six out of 12 tested candidates, allowing the authors to conclude that many of their predictions are real miRNAs. It is worth noting that from 129 candidate miRNAs identified by Xie et al., 76 are overlapping with our predictions, further strengthening the notion that there is still a substantial number of miRNAs in mammalian genomes waiting for experimental confirmation.


    miRNA FUNCTIONS
 TOP
 ABSTRACT
 INTRODUCTION
 miRNA MACHINERY
 IDENTIFICATION OF miRNA GENES
 miRNA FUNCTIONS
 miRNAS AND CANCER
 VIRUS CONNECTIONS
 CONCLUSIONS
 REFERENCES
 
The functions of some miRNAs were elucidated by forward genetics screens (61Go). In C. elegans lin-4 and let-7 miRNAs were shown to be involved in developmental timing (12Go,14Go) and lsy-6 and miR-273 in neuronal cell fate (86Go,87Go); bantam and miR-14 in D. melanogaster are responsible for cell death, proliferation and fat storage (88Go,89Go). Cloning of miRNAs from different tissues and/or establishment of miRNA expression patterns may also provide clues for the function of particular miRNAs. For example, miR-181 was cloned from mouse bone marrow and thymus, where it is predominantly expressed and was consequently shown to be involved in B- and T-cell development (90Go). Poy et al. (91Go) found that miR-375, cloned from pancreatic cells, targets myotrophyn gene, involved in insulin secretion.

Functions of miRNAs can also be predicted computationally— the approach actively explored by different groups (58Go,92Go–98Go). However, the initial algorithms for miRNA target predictions were far from perfect, with little or no overlap among top-predicted targets (99Go), which were attributed to the lack of experimental data necessary to formulate clear rules for miRNA target recognition. To fill this gap, several research groups performed systematic analysis of pairing requirements between miRNA and target mRNA and revealed that first eight bases of miRNA starting from its 5' end are most important for miRNA–mRNA interaction (55Go,100Go,101Go). Surprisingly, these eight bases, also known as miRNA seed, can alone, with little or no pairing in the 3' end of the miRNA, be sufficient for a functional miRNA–mRNA duplex. The importance of the miRNA seed in target recognition was also recognized by computational analysis: it appeared that many of the conserved 8mer motifs in 3' UTRs of genes correspond to miRNA seeds (85Go,93Go,102Go). However, perfect matching at the 5' end of a miRNA is not an absolute requirement, because Brennecke et al. (101Go) found that sites with mismatches in the seed can also be functional if there is strong compensatory pairing at the 3' end of a miRNA. Equipped with the better understanding of miRNA target recognition principles, Kerk et al. (98Go) and Lewis et al. (102Go) recently developed improved algorithms for miRNA target prediction and independently came to the similar estimates that in vertebrates miRNAs can target on average 200 transcripts, resulting in regulation of as much as 30% of human genes. Thus, there is a growing evidence of a much broader involvement of miRNAs in gene regulation than previously anticipated.

Information on miRNA expression patterns will be essential for dissection of miRNA-containing gene regulatory networks. We have recently used a combination of microarray and in situ hybridization techniques to establish expression profiles of 115 conserved vertebrate miRNAs in zebrafish embryos (103Go). In situ hybridizations using LNA (locked-nucleic-acid)-modified probes revealed that most (68%) of the investigated miRNAs are expressed only in specific tissues, e.g. muscles, nerves, sensory organs or digestive system (Fig. 2). In addition, the expression of miRNAs was detected starting at segmentation stages of zebrafish embryo development but not at earlier stages. On the basis of these observations, we concluded that miRNAs are not crucial for early patterning but play a role in morphogenesis and maintenance of tissue identity. Similar conclusions came from the work of Giraldez et al. (104Go), who used germline replacement technique to generate Dicer-deficient zebrafish embryos completely void of mature miRNAs. This trick was necessary to circumvent the problem of maternally inherited miRNAs in Dicer mutants (105Go). It appeared that the maternal-zygotic Dicer null mutants display no abnormalities in early development and undergo normal axis formation and regionalization, i.e. form all major subregions and cell types. However, morphogenesis at later stages of development was severely affected. For example, defects in gastrulation, somitogenesis, cardiovascular morphogenesis and brain development were observed. Interestingly, injection of synthetic miRNAs of a single subfamily (miR-430) was able to rescue the brain morphogenesis defects (104Go), demonstrating a fundamental role of miRNAs in organ formation.



View larger version (58K):
[in this window]
[in a new window]
 
Figure 2. Example of a tissue-specific miRNA expression pattern. (A) Ventral and (B) lateral view of miR-137 whole-mount in situ expression in brain structures and spinal cord of zebrafish embryos. The pictures are kindly provided by Erno Wienholds (103Go).

 
Differential gene expression is a principal driving force behind developmental processes such as morphogenesis. Lim et al. (63Go) has recently shown that a single miRNA can influence the entire transcriptome profile of a cell. The authors transfected human cells with miR-124 and miR-1 miRNAs, which are specifically expressed in brain and muscles, respectively, and demonstrated by microarray analysis that in case of miR-124, the expression profile was shifted towards that of brain, whereas miR-1 transfection resulted in the expression profile characteristic for muscles. These experiments provide additional support for the notion that miRNAs help to establish and maintain the identity of different cell types.


    miRNAS AND CANCER
 TOP
 ABSTRACT
 INTRODUCTION
 miRNA MACHINERY
 IDENTIFICATION OF miRNA GENES
 miRNA FUNCTIONS
 miRNAS AND CANCER
 VIRUS CONNECTIONS
 CONCLUSIONS
 REFERENCES
 
A connection between miRNAs and cancer was suggested not long after the discovery of miRNAs. First, Calin et al. (106Go) found that miRNAs mir-15 and mir-16 are located in a locus 13q14, which is frequently deleted in B-cell chronic lymphocytic leukemias. Extension of that study led to the discovery that, in fact, about half of miRNA genes (98 of 186 studied) are situated in cancer-associated genomic regions (107Go) and expression profiling revealed that both up- and down-regulation of many miRNAs occurs in cancer tissues (108Go–110Go). For example, miRNAs from the let-7 family appeared to be poorly expressed in lung cancers (107Go,111Go) but not in other types of cancer. In addition, Johnson et al. (112Go) demonstrated that let-7 family members regulate RAS oncogenes in human. He et al. (113Go) found that the mir-17-19 cluster of miRNAs is overexpressed in certain human lymphomas and demonstrated that its increased expression, mediated by c-MYC oncogene, accelerates tumor development. O'Donnell et al. (114Go) also showed that c-Myc activates expression of mir-17-19 miRNA cluster, whereas two miRNAs from this cluster, miR-17-5p and miR-20a, negatively regulate E2F1 transcription factor (TF) involved in cell cycle progression. Interestingly, E2F1 and c-Myc are also known to activate each other. Thus the role of c-Myc-regulated miRNAs can be in preventing a positive-feedback loop between c-Myc and E2F1, which otherwise could lead to overexpression of both genes and failure in cell cycle regulation (115Go).

Altered expression of particular miRNAs in different types of cancer can possibly be used for cancer diagnosis and prognosis. Brown et al. (116Go) profiled the expression of miRNAs in more than 60 patients with lung, colon, breast, bladder, pancreatic, prostate or thymus cancer and found that different types of tumor have different miRNA expression signatures. For example, in 70% of the patients with lung cancer, six miRNAs (miR-30a, -126, -143, -146, -188, and -331) were expressed at significantly lower levels, whereas three (miR-21, -189 and -200b) were significantly over-expressed. Similarly, colon cancers had a characteristic set of five over-expressed and four under-expressed miRNAs. Notably, there are many known oncogenes among putative targets of differentially expressed miRNAs. Lu et al. (117Go) demonstrated recently that miRNAs can indeed be developed into potent cancer markers. Expression profiling of 217 mammalian miRNA using bead-based flow cytometry revealed that most miRNAs are downregulated in tumor tissues and that different types of tumors can be distinguished by miRNA expression patterns. Moreover, Lu et al. showed that miRNA profiles can be used for highly reliable diagnosis of tumors of histologically uncertain cellular origin.

Taken together, there is growing evidence that miRNAs can act as both oncogenes and tumor suppressors. Further studies are certainly required to reveal the full breadth of miRNA involvement in cancerogenesis.


    VIRUS CONNECTIONS
 TOP
 ABSTRACT
 INTRODUCTION
 miRNA MACHINERY
 IDENTIFICATION OF miRNA GENES
 miRNA FUNCTIONS
 miRNAS AND CANCER
 VIRUS CONNECTIONS
 CONCLUSIONS
 REFERENCES
 
Besides animals and plants, miRNAs were also recently discovered in viruses (118Go–120Go). Pfeffer et al. (118Go) performed a small RNA profile of human cells infected by Epstein-Barr virus, a member of the human herpesvirus family, to study the role that RNA silencing can play during viral infection in animals. Besides host miRNAs they identified five miRNAs originating from the virus. They next developed a computational algorithm to identify miRNAs in virus genomes and surveyed a wide range of viruses (119Go). Experimental verification of predictions in herpes viruses resulted in identification of 28 novel miRNAs from Kaposi sarcoma-associated virus, mouse gamma-herpes virus and human cytomegalovirus. Unexpectedly, virus miRNAs appeared to have no homology between each other or with human miRNAs. Positions of miRNAs within different virus genomes were also not conserved. Lack of sequence and positional conservation suggests that miRNAs are likely to be recent acquisitions in viral genomes evolved to perform very specific functions. The first example of such a specific function of viral miRNA in a polyomavirus was recently published by Sullivan et al. (120Go), who discovered that simian virus 40 (SV40) encodes a miRNA that regulates one of the viral genes called T antigene. The authors found that SV40 miRNA accumulates in late stages of infection and targets a early viral T antigene mRNA for degradation, thus reducing expression of viral T antigens. Importantly, besides involvement in virus replication, T antigenes also trigger the cytotoxic response of T lymphocytes. Therefore, reduced levels of T antigene expression enhance chances of the virus to escape lysis by T-cells. Although the exact roles of many other viral miRNAs remain to be established, the available evidence suggests that DNA viruses evolved to exploit miRNA silencing pathway for their advantage.

Interestingly, miRNAs were found in different DNA viruses but not in small-genome retroviruses (119Go). In addition, no viral siRNAs were detected in human cells infected with different RNA viruses, whereas RNAi machinery of the cells remained unaffected, suggesting that RNA viruses do not interfere with RNA silencing pathways in mammalian cells (119Go). This is in contrast to plants and insects, where RNA silencing is used as an immune response to RNA viruses (121Go,122Go). It was discovered recently that mammalian cells actually do use RNA silencing against retroviruses but utilize miRNA pathway instead of RNAi. Lecellier et al. (123Go) reported that in cells expressing a suppressor of RNA silencing, accumulation of a retrovirus called primate foamy virus type 1 (PFV-1) was strongly enhanced, indicating involvement of siRNAs or miRNAs in the control of virus replication. However, similar to the results of other researchers, no viral siRNAs were detected in infected cells. Instead, computational analysis suggested that human miRNA miR-32 can target one of viral ORF sequences, and experimental verification proved that predicted miR-32 target sequence is functional. In addition, depletion of miR-32 from infected cells resulted in an increased accumulation of viral mRNA, proving antiviral properties of miR-32. Leicellier et al. (123Go) also demonstrated that one of PFV-1 proteins, Tas, is a general suppressor of RNA silencing that can offset the miR-32 effect. Thus, a picture of host–virus interactions is emerging in which cellular miRNAs play an important role in limiting virus replication, and viruses counteract this cellular defense by hampering RNA silencing pathways. The authors suggested that miRNAs might be broadly implicated in viral infection of mammalian cells, and virtually every miRNA, disregarding its primary cellular function, can have a fortuitous antiviral potential (123Go). Indeed, as matching of eight bases in the seed of a miRNA to a target sequence is often sufficient for a functional interaction (101Go), host's miRNAs relations to viruses may be seen as a mammalian version of the bacterial restriction enzyme defense system against phages. To test this hypothesis and to establish the extent of host's miRNA involvement in virus regulation, a survey of a broad range of viruses will be required.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 miRNA MACHINERY
 IDENTIFICATION OF miRNA GENES
 miRNA FUNCTIONS
 miRNAS AND CANCER
 VIRUS CONNECTIONS
 CONCLUSIONS
 REFERENCES
 
It is no exaggeration to say that we are evidencing a paradigm shift in our understanding of the complexity in gene regulatory networks, provoked by the discovery of miRNAs. The recent findings that miRNAs provide a new layer in gene regulation, additional to TFs and other regulatory systems, constitute the essence of this shift. The similarity in the logics of action of TFs and miRNAs was discussed by Hobert (124Go), who noticed that different combinations of miRNAs in different cell types, or ‘miRNA codes’, can have roles conceptually similar to ‘TF codes’ in establishing these cell types. Bartel and Chen (125Go) speculated that differential expression of miRNAs in different cell types provides a ‘miRNA milieu’ that influences expression of thousands of mRNAs. Recent experimental data confirm these concepts (63Go,103Go,104Go). Tissue-specific expression of many miRNAs, their necessity for morphogenesis but not for early development, and ability to shift expression profiles of cells to that of specific tissues, suggest that establishment and maintenance of cells in specific differentiated states is the main function of many miRNAs. Thus, by analogy to housekeeping genes that provide basic functions in all the organisms, miRNAs can be called roomkeeping genes, because they are to a large extent responsible for keeping differentiated cells in order, by dumping expression of genes that should not be expressed in particular cell types. Connection between miRNAs and cancer immediately becomes obvious from this analogy: perturbations in miRNA expression patterns may lead to cell dedifferentiation and can be a cause as well as a consequence of tumorogenesis. Apparently, roomkeeping is not the only function of miRNAs. The recently discovered involvement of miRNAs in host–virus coevolution is spectacular: at least some miRNAs evolved to play defensive roles against viruses, while viruses also evolved to use miRNA machinery for their advantage.

Although a steady progress in understanding the miRNA world has been made in recent years, there is still a long way ahead before miRNAs reveal all their secrets. The issue on the number of miRNAs in the human genome is not settled and continuing cloning and other experimental efforts are required to identify all miRNA genes. Knowledge of the complete microRNome is essential for experimental and computational dissection of miRNA-containing gene regulatory networks. Knockout and knockdown experiments are required to reveal specific functions of particular miRNA genes, and technologies for this type of experiments are available (126Go). Although miRNA knockouts/knockdowns will provide loss of function phenotypes, computational approaches will be necessary to assist in identification of mRNA targets behind these phenotypes and further improvements in miRNA target prediction algorithms can be expected. As computational predictions require an experimental confirmation, a high-throughput technology for miRNA target verification, if developed, will greatly benefit the field. Finally, we foresee a substantial increase in volume of research on connection between miRNAs and human disease.

Conflict of Interest statement. None declared.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 miRNA MACHINERY
 IDENTIFICATION OF miRNA GENES
 miRNA FUNCTIONS
 miRNAS AND CANCER
 VIRUS CONNECTIONS
 CONCLUSIONS
 REFERENCES
 

  1. Mattick, J.S. (2001) Non-coding RNAs: the architects of eukaryotic complexity. EMBO Rep., 2, 986–991.[CrossRef][ISI][Medline]

  2. Dennis, C. (2002) The brave new world of RNA. Nature, 418, 122–124.[CrossRef][Medline]

  3. Mattick, J.S. (2003) Challenging the dogma: the hidden layer of non-protein-coding RNAs in complex organisms. Bioessays, 25, 930–939.[CrossRef][ISI][Medline]

  4. Cheng, J., Kapranov, P., Drenkow, J., Dike, S., Brubaker, S., Patel, S., Long, J., Stern, D., Tammana, H., Helt, G. et al. (2005) Transcriptional maps of 10 human chromosomes at 5-nucleotide resolution. Science, 308, 1149–1154.[Abstract/Free Full Text]

  5. Griffiths-Jones, S., Moxon, S., Marshall, M., Khanna, A., Eddy, S.R. and Bateman, A. (2005) Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids Res., 33, D121–D124.[Abstract/Free Full Text]

  6. Liu, C., Bai, B., Skogerbo, G., Cai, L., Deng, W., Zhang, Y., Bu, D., Zhao, Y. and Chen, R. (2005) NONCODE: an integrated knowledge database of non-coding RNAs. Nucleic Acids Res., 33, D112–D115.[Abstract/Free Full Text]

  7. Pang, K.C., Stephen, S., Engstrom, P.G., Tajul-Arifin, K., Chen, W., Wahlestedt, C., Lenhard, B., Hayashizaki, Y. and Mattick, J.S. (2005) RNAdb-a comprehensive mammalian noncoding RNA database. Nucleic Acids Res., 33, D125–D130.[Abstract/Free Full Text]

  8. Bachellerie, J.P., Cavaille, J. and Huttenhofer, A. (2002) The expanding snoRNA world. Biochimie, 84, 775–790.[Medline]

  9. Kiss, T. (2002) Small nucleolar RNAs: an abundant group of noncoding RNAs with diverse cellular functions. Cell, 109, 145–148.[CrossRef][ISI][Medline]

  10. Mattick, J.S. and Makunin, I.V. (2005) Small regulatory RNAs in mammals. Hum. Mol. Genet., 14, Spec no 1, R121–R132.[Abstract/Free Full Text]

  11. Huttenhofer, A., Schattner, P. and Polacek, N. (2005) Non-coding RNAs: hope or hype?. Trends Genet., 21, 289–297.[CrossRef][ISI][Medline]

  12. Lee, R.C., Feinbaum, R.L. and Ambros, V. (1993) The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell, 75, 843–854.[CrossRef][ISI][Medline]

  13. Wightman, B., Ha, I. and Ruvkun, G. (1993) Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans. Cell, 75, 855–862.[CrossRef][ISI][Medline]

  14. Reinhart, B.J., Slack, F.J., Basson, M., Pasquinelli, A.E., Bettinger, J.C., Rougvie, A.E., Horvitz, H.R. and Ruvkun, G. (2000) The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature, 403, 901–906.[CrossRef][Medline]

  15. Ruvkun, G., Wightman, B., Ha, I., Arasu, P. et al. (2004) The 20 years it took to recognize the importance of tiny RNAs. Cell, S93–S96.

  16. Pasquinelli, A.E., Reinhart, B.J., Slack, F., Martindale, M.Q., Kuroda, M.I., Maller, B., Hayward, D.C., Ball, E.E., Degnan, B., Muller, P. et al. (2000) Conservation of the sequence and temporal expression of let-7 heterochronic regulatory RNA. Nature, 408, 86–89.[CrossRef][Medline]

  17. Slack, F.J., Basson, M., Liu, Z., Ambros, V., Horvitz, H.R. and Ruvkun, G. (2000) The lin-41 RBCC gene acts in the C. elegans heterochronic pathway between the let-7 regulatory RNA and the LIN-29 transcription factor. Mol. Cell, 5, 659–669.[CrossRef][ISI][Medline]

  18. Lagos-Quintana, M., Rauhut, R., Lendeckel, W. and Tuschl, T. (2001) Identification of novel genes coding for small expressed RNAs. Science, 294, 853–858.[Abstract/Free Full Text]

  19. Lau, N.C., Lim, L.P., Weinstein, E.G. and Bartel, D.P. (2001) An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans. Science, 294, 858–862.[Abstract/Free Full Text]

  20. Lee, R.C. and Ambros, V. (2001) An extensive class of small RNAs in Caenorhabditis elegans. Science, 294, 862–864.[Abstract/Free Full Text]

  21. Bartel, D.P. (2004) MicroRNAs: genomics, biogenesis, mechanism, and function. Cell, 116, 281–297.[CrossRef][ISI][Medline]

  22. He, L. and Hannon, G.J. (2004) MicroRNAs: small RNAs with a big role in gene regulation. Nat. Rev. Genet., 5, 522–531.[CrossRef][Medline]

  23. Mallory, A.C. and Vaucheret, H. (2004) MicroRNAs: something important between the genes. Curr. Opin. Plant Biol., 7, 120–125.[CrossRef][ISI][Medline]

  24. Murchison, E.P. and Hannon, G.J. (2004) miRNAs on the move: miRNA biogenesis and the RNAi machinery. Curr. Opin. Cell Biol., 16, 223–229.[CrossRef][ISI][Medline]

  25. Ruvkun, G., Wightman, B. and Ha, I. (2004) The 20 years it took to recognize the importance of tiny RNAs. Cell, 116, S93–S96.[CrossRef][Medline]

  26. Kim, V.N. (2005) MicroRNA biogenesis: coordinated cropping and dicing. Nat. Rev. Mol. Cell Biol., 6, 376–385.[CrossRef][ISI][Medline]

  27. Cai, X., Hagedorn, C.H. and Cullen, B.R. (2004) Human microRNAs are processed from capped, polyadenylated transcripts that can also function as mRNAs. RNA, 10, 1957–1966.[Abstract/Free Full Text]

  28. Lee, Y., Kim, M., Han, J., Yeom, K.H., Lee, S., Baek, S.H. and Kim, V.N. (2004) MicroRNA genes are transcribed by RNA polymerase II. EMBO J., 23, 4051–4060.[CrossRef][ISI][Medline]

  29. Lee, Y., Jeon, K., Lee, J.T., Kim, S. and Kim, V.N. (2002) MicroRNA maturation: stepwise processing and subcellular localization. EMBO J., 21, 4663–4670.[CrossRef][ISI][Medline]

  30. Rodriguez, A., Griffiths-Jones, S., Ashurst, J.L. and Bradley, A. (2004) Identification of mammalian microRNA host genes and transcription units. Genome Res., 14, 1902–1910.[Abstract/Free Full Text]

  31. Denli, A.M., Tops, B.B., Plasterk, R.H., Ketting, R.F. and Hannon, G.J. (2004) Processing of primary microRNAs by the microprocessor complex. Nature, 432, 231–235.[CrossRef][Medline]

  32. Gregory, R.I., Yan, K.P., Amuthan, G., Chendrimada, T., Doratotaj, B., Cooch, N. and Shiekhattar, R. (2004) The microprocessor complex mediates the genesis of microRNAs. Nature, 432, 235–240.[CrossRef][Medline]

  33. Lee, Y., Ahn, C., Han, J., Choi, H., Kim, J., Yim, J., Lee, J., Provost, P., Radmark, O., Kim, S., et al. (2003) The nuclear RNase III Drosha initiates microRNA processing. Nature, 425, 415–419.[CrossRef][Medline]

  34. Han, J., Lee, Y., Yeom, K.H., Kim, Y.K., Jin, H. and Kim, V.N. (2004) The Drosha-DGCR8 complex in primary microRNA processing. Genes Dev., 18, 3016–3027.[Abstract/Free Full Text]

  35. Landthaler, M., Yalcin, A. and Tuschl, T. (2004) The human DiGeorge syndrome critical region gene 8 and its D. melanogaster homolog are required for miRNA biogenesis. Curr. Biol., 14, 2162–2167.[CrossRef][ISI][Medline]

  36. Zeng, Y. and Cullen, B.R. (2003) Sequence requirements for micro RNA processing and function in human cells. RNA, 9, 112–123.[Abstract/Free Full Text]

  37. Zeng, Y., Yi, R. and Cullen, B.R. (2005) Recognition and cleavage of primary microRNA precursors by the nuclear processing enzyme Drosha. EMBO J., 24, 138–148.[CrossRef][ISI][Medline]

  38. Yi, R., Qin, Y., Macara, I.G. and Cullen, B.R. (2003) Exportin-5 mediates the nuclear export of pre-microRNAs and short hairpin RNAs. Genes Dev., 17, 3011–3016.[Abstract/Free Full Text]

  39. Bohnsack, M.T., Czaplinski, K. and Gorlich, D. (2004) Exportin 5 is a RanGTP-dependent dsRNA-binding protein that mediates nuclear export of pre-miRNAs. RNA, 10, 185–191.[Abstract/Free Full Text]

  40. Lund, E., Guttinger, S., Calado, A., Dahlberg, J.E. and Kutay, U. (2004) Nuclear export of microRNA precursors. Science, 303, 95–98.[Abstract/Free Full Text]

  41. Grishok, A., Pasquinelli, A.E., Conte, D., Li, N., Parrish, S., Ha, I., Baillie, D.L., Fire, A., Ruvkun, G. and Mello, C.C. (2001) Genes and mechanisms related to RNA interference regulate expression of the small temporal RNAs that control C. elegans developmental timing. Cell, 106, 23–34.[CrossRef][ISI][Medline]

  42. Hutvagner, G., McLachlan, J., Pasquinelli, A.E., Balint, E., Tuschl, T. and Zamore, P.D. (2001) A cellular function for the RNA-interference enzyme Dicer in the maturation of the let-7 small temporal RNA. Science, 293, 834–838.[Abstract/Free Full Text]

  43. Ketting, R.F., Fischer, S.E., Bernstein, E., Sijen, T., Hannon, G.J. and Plasterk, R.H. (2001) Dicer functions in RNA interference and in synthesis of small RNA involved in developmental timing in C. elegans. Genes Dev., 15, 2654–2659.[Abstract/Free Full Text]

  44. Khvorova, A., Reynolds, A. and Jayasena, S.D. (2003) Functional siRNAs and miRNAs exhibit strand bias. Cell, 115, 209–216.[CrossRef][ISI][Medline]

  45. Schwarz, D.S., Hutvagner, G., Du, T., Xu, Z., Aronin, N. and Zamore, P.D. (2003) Asymmetry in the assembly of the RNAi enzyme complex. Cell, 115, 199–208.[CrossRef][ISI][Medline]

  46. Hammond, S.M., Boettcher, S., Caudy, A.A., Kobayashi, R. and Hannon, G.J. (2001) Argonaute2, a link between genetic and biochemical analyses of RNAi. Science, 293, 1146–1150.[Abstract/Free Full Text]

  47. Hutvagner, G. and Zamore, P.D. (2002) A microRNA in a multiple-turnover RNAi enzyme complex. Science, 297, 2056–2060.[Abstract/Free Full Text]

  48. Martinez, J., Patkaniowska, A., Urlaub, H., Luhrmann, R. and Tuschl, T. (2002) Single-stranded antisense siRNAs guide target RNA cleavage in RNAi. Cell, 110, 563–574.[CrossRef][ISI][Medline]

  49. Mourelatos, Z., Dostie, J., Paushkin, S., Sharma, A., Charroux, B., Abel, L., Rappsilber, J., Mann, M. and Dreyfuss, G. (2002) miRNPs: a novel class of ribonucleoproteins containing numerous microRNAs. Genes Dev., 16, 720–728.[Abstract/Free Full Text]

  50. Caudy, A.A., Ketting, R.F., Hammond, S.M., Denli, A.M., Bathoorn, A.M., Tops, B.B., Silva, J.M., Myers, M.M., Hannon, G.J. and Plasterk, R.H. (2003) A micrococcal nuclease homologue in RNAi effector complexes. Nature, 425, 411–414.[CrossRef][Medline]

  51. Caudy, A.A., Myers, M., Hannon, G.J. and Hammond, S.M. (2002) Fragile X-related protein and VIG associate with the RNA interference machinery. Genes Dev., 16, 2491–2496.[Abstract/Free Full Text]

  52. Jin, P., Zarnescu, D.C., Ceman, S., Nakamoto, M., Mowrey, J., Jongens, T.A., Nelson, D.L., Moses, K. and Warren, S.T. (2004) Biochemical and genetic interaction between the fragile X mental retardation protein and the microRNA pathway. Nat. Neurosci., 7, 113–117.[CrossRef][ISI][Medline]

  53. Liu, J., Carmell, M.A., Rivas, F.V., Marsden, C.G., Thomson, J.M., Song, J.J., Hammond, S.M., Joshua-Tor, L. and Hannon, G.J. (2004) Argonaute2 is the catalytic engine of mammalian RNAi. Science, 305, 1437–1441.[Abstract/Free Full Text]

  54. Zeng, Y., Wagner, E.J. and Cullen, B.R. (2002) Both natural and designed micro RNAs can inhibit the expression of cognate mRNAs when expressed in human cells. Mol. Cell, 9, 1327–1333.[CrossRef][ISI][Medline]

  55. Doench, J.G., Petersen, C.P. and Sharp, P.A. (2003) siRNAs can function as miRNAs. Genes Dev., 17, 438–442.[Abstract/Free Full Text]

  56. Zeng, Y., Yi, R. and Cullen, B.R. (2003) MicroRNAs and small interfering RNAs can inhibit mRNA expression by similar mechanisms. Proc. Natl Acad. Sci. USA, 100, 9779–9784.[Abstract/Free Full Text]

  57. Llave, C., Xie, Z., Kasschau, K.D. and Carrington, J.C. (2002) Cleavage of scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA. Science, 297, 2053–2056.[Abstract/Free Full Text]

  58. Rhoades, M.W., Reinhart, B.J., Lim, L.P., Burge, C.B., Bartel, B. and Bartel, D.P. (2002) Prediction of plant microRNA targets. Cell, 110, 513–520.[CrossRef][ISI][Medline]

  59. Kasschau, K.D., Xie, Z., Allen, E., Llave, C., Chapman, E.J., Krizan, K.A. and Carrington, J.C. (2003) P1/HC-Pro, a viral suppressor of RNA silencing, interferes with Arabidopsis development and miRNA function. Dev. Cell, 4, 205–217.[CrossRef][ISI][Medline]

  60. Tang, G., Reinhart, B.J., Bartel, D.P. and Zamore, P.D. (2003) A biochemical framework for RNA silencing in plants. Genes Dev., 17, 49–63.[Abstract/Free Full Text]

  61. Ambros, V. (2004) The functions of animal microRNAs. Nature, 431, 350–355.[CrossRef][Medline]

  62. Yekta, S., Shih, I.H. and Bartel, D.P. (2004) MicroRNA-directed cleavage of HOXB8 mRNA. Science, 304, 594–596.[Abstract/Free Full Text]

  63. Lim, L.P., Lau, N.C., Garrett-Engele, P., Grimson, A., Schelter, J.M., Castle, J., Bartel, D.P., Linsley, P.S. and Johnson, J.M. (2005) Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature, 433, 769–773.[CrossRef][Medline]

  64. Filipowicz, W., Jaskiewicz, L., Kolb, F.A. and Pillai, R.S. (2005) Post-transcriptional gene silencing by siRNAs and miRNAs. Curr. Opin. Struct. Biol., 15, 331–341.[CrossRef][ISI][Medline]

  65. Griffiths-Jones, S. (2004) The microRNA Registry. Nucleic Acids Res., 32 Database issue, D109–11.[Abstract/Free Full Text]

  66. Lagos-Quintana, M., Rauhut, R., Yalcin, A., Meyer, J., Lendeckel, W. and Tuschl, T. (2002) Identification of tissue-specific microRNAs from mouse. Curr. Biol., 12, 735–739.[CrossRef][ISI][Medline]

  67. Ambros, V., Lee, R.C., Lavanway, A., Williams, P.T. and Jewell, D. (2003) MicroRNAs and other tiny endogenous RNAs in C. elegans. Curr. Biol., 13, 807–818.[CrossRef][ISI][Medline]

  68. Aravin, A.A., Lagos-Quintana, M., Yalcin, A., Zavolan, M., Marks, D., Snyder, B., Gaasterland, T., Meyer, J. and Tuschl, T. (2003) The small RNA profile during Drosophila melanogaster development. Dev. Cell, 5, 337–350.[CrossRef][ISI][Medline]

  69. Dostie, J., Mourelatos, Z., Yang, M., Sharma, A. and Dreyfuss, G. (2003) Numerous microRNPs in neuronal cells containing novel microRNAs. RNA, 9, 180–186.[Abstract/Free Full Text]

  70. Houbaviy, H.B., Murray, M.F. and Sharp, P.A. (2003) Embryonic stem cell-specific MicroRNAs. Dev. Cell, 5, 351–358.[CrossRef][ISI][Medline]

  71. Lagos-Quintana, M., Rauhut, R., Meyer, J., Borkhardt, A. and Tuschl, T. (2003) New microRNAs from mouse and human. RNA, 9, 175–179.[Abstract/Free Full Text]

  72. Kim, J., Krichevsky, A., Grad, Y., Hayes, G.D., Kosik, K.S., Church, G.M. and Ruvkun, G. (2004) Identification of many microRNAs that copurify with polyribosomes in mammalian neurons. Proc. Natl Acad. Sci. USA, 101, 360–365.[Abstract/Free Full Text]

  73. Suh, M.R., Lee, Y., Kim, J.Y., Kim, S.K., Moon, S.H., Lee, J.Y., Cha, K.Y., Chung, H.M., Yoon, H.S., Moon, S.Y. et al. (2004) Human embryonic stem cells express a unique set of microRNAs. Dev. Biol., 270, 488–498.[CrossRef][ISI][Medline]

  74. Grad, Y., Aach, J., Hayes, G.D., Reinhart, B.J., Church, G.M., Ruvkun, G. and Kim, J. (2003) Computational and experimental identification of C. elegans microRNAs. Mol. Cell, 11, 1253–1263.[CrossRef][ISI][Medline]

  75. Lai, E.C., Tomancak, P., Williams, R.W. and Rubin, G.M. (2003) Computational identification of Drosophila microRNA genes. Genome Biol., 4, R42.[CrossRef][Medline]

  76. Lim, L.P., Glasner, M.E., Yekta, S., Burge, C.B. and Bartel, D.P. (2003) Vertebrate microRNA genes. Science, 299, 1540.[Free Full Text]

  77. Lim, L.P., Lau, N.C., Weinstein, E.G., Abdelhakim, A., Yekta, S., Rhoades, M.W., Burge, C.B. and Bartel, D.P. (2003) The microRNAs of Caenorhabditis elegans. Genes Dev., 17, 991–1008.[Abstract/Free Full Text]

  78. Jones-Rhoades, M.W. and Bartel, D.P. (2004) Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol. Cell, 14, 787–799.[CrossRef][ISI][Medline]

  79. Wang, X.J., Reyes, J.L., Chua, N.H. and Gaasterland, T. (2004) Prediction and identification of Arabidopsis thaliana microRNAs and their mRNA targets. Genome Biol., 5, R65.[CrossRef][Medline]

  80. Berezikov, E., Guryev, V., van de, BeJ., Wienholds, E., Plasterk, R.H. and Cuppen, E. (2005) Phylogenetic shadowing and computational identification of human microRNA genes. Cell, 120, 21–24.[CrossRef][ISI][Medline]

  81. Brennecke, J. and Cohen, S.M. (2003) Towards a complete description of the microRNA complement of animal genomes. Genome Biol., 4, 228.[CrossRef][Medline]