MicroRNAs (miRNAs) are endogenous approximately 22 nt RNAs that can play important regulatory roles in animals and plants by targeting mRNAs for cleavage or translational repression. Although they escaped notice until relatively recently, miRNAs comprise one of the more abundant classes of gene regulatory molecules in multicellular organisms and likely influence the output of many protein-coding genes.

In Caenorhabditis elegans, lin-4 and let-7 encode 22- and 21-nucleotide (nt) RNAs, respectively, which function as key regulators of developmental timing. Because the appearance of these short RNAs is regulated during development, they are also referred to as small temporal RNAs (stRNAs). We show that many 21- and 22-nt expressed RNAs, termed microRNAs, exist in invertebrates and vertebrates and that some of these novel RNAs, similar to let-7 stRNA, are highly conserved. This suggests that sequence-specific, posttranscriptional regulatory mechanisms mediated by small RNAs are more general than previously appreciated.

BACKGROUND

The recent discoveries of microRNA (miRNA) genes and characterization of the first few target genes regulated by miRNAs in Caenorhabditis elegans and Drosophila melanogaster have set the stage for elucidation of a novel network of regulatory control. We present a computational method for whole-genome prediction of miRNA target genes. The method is validated using known examples. For each miRNA, target genes are selected on the basis of three properties: sequence complementarity using a position-weighted local alignment algorithm, free energies of RNA-RNA duplexes, and conservation of target sites in related genomes. Application to the D. melanogaster, Drosophila pseudoobscura and Anopheles gambiae genomes identifies several hundred target genes potentially regulated by one or more known miRNAs.

RESULTS

These potential targets are rich in genes that are expressed at specific developmental stages and that are involved in cell fate specification, morphogenesis and the coordination of developmental processes, as well as genes that are active in the mature nervous system. High-ranking target genes are enriched in transcription factors two-fold and include genes already known to be under translational regulation. Our results reaffirm the thesis that miRNAs have an important role in establishing the complex spatial and temporal patterns of gene activity necessary for the orderly progression of development and suggest additional roles in the function of the mature organism. In addition the results point the way to directed experiments to determine miRNA functions.

CONCLUSIONS

The emerging combinatorics of miRNA target sites in the 3' untranslated regions of messenger RNAs are reminiscent of transcriptional regulation in promoter regions of DNA, with both one-to-many and many-to-one relationships between regulator and target. Typically, more than one miRNA regulates one message, indicative of cooperative translational control. Conversely, one miRNA may have several target genes, reflecting target multiplicity. As a guide to focused experiments, we provide detailed online information about likely target genes and binding sites in their untranslated regions, organized by miRNA or by gene and ranked by likelihood of match. The target prediction algorithm is freely available and can be applied to whole genome sequences using identified miRNA sequences.

A key step in RNA interference (RNAi) is assembly of the RISC, the protein-siRNA complex that mediates target RNA cleavage. Here, we show that the two strands of an siRNA duplex are not equally eligible for assembly into RISC. Rather, both the absolute and relative stabilities of the base pairs at the 5' ends of the two siRNA strands determine the degree to which each strand participates in the RNAi pathway. siRNA duplexes can be functionally asymmetric, with only one of the two strands able to trigger RNAi. Asymmetry is the hallmark of a related class of small, single-stranded, noncoding RNAs, microRNAs (miRNAs). We suggest that single-stranded miRNAs are initially generated as siRNA-like duplexes whose structures predestine one strand to enter the RISC and the other strand to be destroyed. Thus, the common step of RISC assembly is an unexpected source of asymmetry for both siRNA function and miRNA biogenesis.

MicroRNAs are small noncoding RNAs that serve as posttranscriptional regulators of gene expression in higher eukaryotes. Their widespread and important role in animals is highlighted by recent estimates that 20%-30% of all genes are microRNA targets. Here, we report that a large set of genes involved in basic cellular processes avoid microRNA regulation due to short 3'UTRs that are specifically depleted of microRNA binding sites. For individual microRNAs, we find that coexpressed genes avoid microRNA sites, whereas target genes and microRNAs are preferentially expressed in neighboring tissues. This mutually exclusive expression argues that microRNAs confer accuracy to developmental gene-expression programs, thus ensuring tissue identity and supporting cell-lineage decisions.

Small RNAs ranging in size between 20 and 30 nucleotides are involved in different types of regulation of gene expression including mRNA degradation, translational repression, and chromatin modification. Here we describe the small RNA profile of Drosophila melanogaster as a function of development. We have cloned and sequenced over 4000 small RNAs, 560 of which have the characteristics of RNase III cleavage products. A nonredundant set of 62 miRNAs was identified. We also isolated 178 repeat-associated small interfering RNAs (rasiRNAs), which are cognate to transposable elements, satellite and microsatellite DNA, and Suppressor of Stellate repeats, suggesting that small RNAs participate in defining chromatin structure. rasiRNAs are most abundant in testes and early embryos, where regulation of transposon activity is critical and dramatic changes in heterochromatin structure occur.

Sequencing of multiple related species followed by comparative genomics analysis constitutes a powerful approach for the systematic understanding of any genome. Here, we use the genomes of 12 Drosophila species for the de novo discovery of functional elements in the fly. Each type of functional element shows characteristic patterns of change, or 'evolutionary signatures', dictated by its precise selective constraints. Such signatures enable recognition of new protein-coding genes and exons, spurious and incorrect gene annotations, and numerous unusual gene structures, including abundant stop-codon readthrough. Similarly, we predict non-protein-coding RNA genes and structures, and new microRNA (miRNA) genes. We provide evidence of miRNA processing and functionality from both hairpin arms and both DNA strands. We identify several classes of pre- and post-transcriptional regulatory motifs, and predict individual motif instances with high confidence. We also study how discovery power scales with the divergence and number of species compared, and we provide general guidelines for comparative studies.

MicroRNA (miRNA) genes give rise to small regulatory RNAs in a wide variety of organisms. We used computational methods to predict miRNAs conserved among Drosophila species and large-scale sequencing of small RNAs from Drosophila melanogaster to experimentally confirm and complement these predictions. In addition to validating 20 of our top 45 predictions for novel miRNA loci, the large-scale sequencing identified many miRNAs that had not been predicted. In total, 59 novel genes were identified, increasing our tally of confirmed fly miRNAs to 148. The large-scale sequencing also refined the identities of previously known miRNAs and provided insights into their biogenesis and expression. Many miRNAs were expressed in particular developmental contexts, with a large cohort of miRNAs expressed primarily in imaginal discs. Conserved miRNAs typically were expressed more broadly and robustly than were nonconserved miRNAs, and those conserved miRNAs with more restricted expression tended to have fewer predicted targets than those expressed more broadly. Predicted targets for the expanded set of microRNAs substantially increased and revised the miRNA-target relationships that appear conserved among the fly species. Insights were also provided into miRNA gene evolution, including evidence for emergent regulatory function deriving from the opposite arm of the miRNA hairpin, exemplified by mir-10, and even the opposite strand of the DNA, exemplified by mir-iab-4.

With their power to shape animal morphology, few genes have captured the imagination of biologists as the evolutionarily conserved members of the Hox clusters have done. Recent research has provided new insight into how Hox proteins cause morphological diversity at the organismal and evolutionary levels. Furthermore, an expanding collection of sequences that are directly regulated by Hox proteins provides information on the specificity of target-gene activation, which might allow the successful prediction of novel Hox-response genes. Finally, the recent discovery of microRNA genes within the Hox gene clusters indicates yet another level of control by Hox genes in development and evolution.

BACKGROUND

MicroRNAs (miRNAs) are a large family of 21-22 nucleotide non-coding RNAs with presumed post-transcriptional regulatory activity. Most miRNAs were identified by direct cloning of small RNAs, an approach that favors detection of abundant miRNAs. Three observations suggested that miRNA genes might be identified using a computational approach. First, miRNAs generally derive from precursor transcripts of 70-100 nucleotides with extended stem-loop structure. Second, miRNAs are usually highly conserved between the genomes of related species. Third, miRNAs display a characteristic pattern of evolutionary divergence.

RESULTS

We developed an informatic procedure called 'miRseeker', which analyzed the completed euchromatic sequences of Drosophila melanogaster and D. pseudoobscura for conserved sequences that adopt an extended stem-loop structure and display a pattern of nucleotide divergence characteristic of known miRNAs. The sensitivity of this computational procedure was demonstrated by the presence of 75% (18/24) of previously identified Drosophila miRNAs within the top 124 candidates. In total, we identified 48 novel miRNA candidates that were strongly conserved in more distant insect, nematode, or vertebrate genomes. We verified expression for a total of 24 novel miRNA genes, including 20 of 27 candidates conserved in a third species and 4 of 11 high-scoring, Drosophila-specific candidates. Our analyses lead us to estimate that drosophilid genomes contain around 110 miRNA genes.

CONCLUSIONS

Our computational strategy succeeded in identifying bona fide miRNA genes and suggests that miRNAs constitute nearly 1% of predicted protein-coding genes in Drosophila, a percentage similar to the percentage of miRNAs recently attributed to other metazoan genomes.

During microRNA (miRNA) biogenesis, one strand of a approximately 21-22-nucleotide RNA duplex is preferentially selected for entry into a silencing complex. The other strand, known as the miRNA* species, has typically been assumed to be a carrier strand. Here we show that, although Drosophila melanogaster miRNA* species are less abundant than their partners, they are often present at physiologically relevant levels and can associate with Argonaute proteins. Comparative genomic analyses revealed that >40% of miRNA* sequences resist nucleotide divergence across Drosophilid evolution, and at least half of these well-conserved miRNA* species select for conserved 3' untranslated region seed matches well above background noise. Finally, we validated the inhibitory activity of miRNA* species in both cultured cells and transgenic animals. These data broaden the reach of the miRNA regulatory network and suggest an important mechanism that diversifies miRNA function during evolution.

We present a fluorescence-based, multiplex in situ hybridization method that permits the simultaneous detection of five differently labeled antisense RNA probes and up to seven differ-ent transcripts in a single Drosophila embryo. We also show that it should be possible to increase the number of detected transcripts substantially with nascent transcript multiplex fluorescent in situ hybridization. These multiplex methods fill a current technological gap between high-resolution in situ hybridization with one or two fluorescently labeled probes and low-resolution but genome-wide microarray RNA profiling and should be of great utility in establishing gene networks.

In flies, small silencing RNAs are sorted between Argonaute1 (Ago1), the central protein component of the microRNA (miRNA) pathway, and Argonaute2 (Ago2), which mediates RNA interference. Extensive double-stranded character-as is found in small interfering RNAs (siRNAs)-directs duplexes into Ago2, whereas central mismatches, like those found in miRNA/miRNA* duplexes, direct duplexes into Ago1. Central to this sorting decision is the affinity of the small RNA duplex for the Dcr-2/R2D2 heterodimer, which loads small RNAs into Ago2. Here, we show that while most Drosophila miRNAs are bound to Ago1, miRNA* strands accumulate bound to Ago2. Like siRNA loading, efficient loading of miRNA* strands in Ago2 favors duplexes with a paired central region and requires both Dcr-2 and R2D2. Those miRNA and miRNA* sequences bound to Ago2, like siRNAs diced in vivo from long double-stranded RNA, typically begin with cytidine, whereas Ago1-bound miRNA and miRNA* disproportionately begin with uridine. Consequently, some pre-miRNA generate two or more isoforms from the same side of the stem that differentially partition between Ago1 and Ago2. Our findings provide the first genome-wide test for the idea that Drosophila small RNAs are sorted between Ago1 and Ago2 according to their duplex structure and the identity of their first nucleotide.

microRNAs are short, endogenous non-coding RNAs that act as post-transcriptional modulators of gene expression. Important functions for microRNAs have been found in the regulation of development, cellular proliferation and differentiation, while perturbed miRNA expression patterns have been observed in many human cancers. Here we present a method for specific inhibition of miRNA function through interaction with LNA-modified antisense oligonucleotides and report the specificity of this application. We show that LNA-modified oligonucleotides can inhibit exogenously introduced miRNAs with high specificity using a heterologous reporter assay, and furthermore demonstrate their ability to inhibit an endogenous miRNA in Drosophila melanogaster cells, leading to up-regulation of the cognate target protein. The method shows stoichiometric and reliable inhibition of the targeted miRNA and can thus be applied to studies of miRNA functions and validation of putative target genes.

Human neurodegenerative diseases have the temporal hallmark of afflicting the elderly population. Ageing is one of the most prominent factors to influence disease onset and progression, yet little is known about the molecular pathways that connect these processes. To understand this connection it is necessary to identify the pathways that functionally integrate ageing, chronic maintenance of the brain and modulation of neurodegenerative disease. MicroRNAs (miRNA) are emerging as critical factors in gene regulation during development; however, their role in adult-onset, age-associated processes is only beginning to be revealed. Here we report that the conserved miRNA miR-34 regulates age-associated events and long-term brain integrity in Drosophila, providing a molecular link between ageing and neurodegeneration. Fly mir-34 expression exhibits adult-onset, brain-enriched and age-modulated characteristics. Whereas mir-34 loss triggers a gene profile of accelerated brain ageing, late-onset brain degeneration and a catastrophic decline in survival, mir-34 upregulation extends median lifespan and mitigates neurodegeneration induced by human pathogenic polyglutamine disease protein. Some of the age-associated effects of miR-34 require adult-onset translational repression of Eip74EF, an essential ETS domain transcription factor involved in steroid hormone pathways. Our studies indicate that miRNA-dependent pathways may have an impact on adult-onset, age-associated events by silencing developmental genes that later have a deleterious influence on adult life cycle and disease, and highlight fly miR-34 as a key miRNA with a role in this process.

lin-4 and let-7 are founding members of an extensive family of genes that produce small transcripts, termed microRNAs (miRNAs). In Caenorhabditis elegans, lin-4 and let-7 control the timing of postembryonic events by translational repression of target genes, permitting progression from early to late developmental programs. To identify Drosophila melanogaster miRNAs that could play similar roles in the control of developmental timing, we characterized the developmental expression profile of 24 miRNAs in Drosophila, and found 7 miRNAs that are either upregulated or downregulated in conjunction with metamorphosis. The upregulation of three of these miRNAs (mir-100, mir-125, and let-7), and the downregulation of a fourth (mir-34) requires the hormone ecdysone (Ecd) and the activity of the Ecd-inducible gene Broad-Complex. Interestingly, mir-125 is a putative homologue of lin-4. mir-100, -125, and let-7 are clustered within an 800-bp region on chromosome 2L, suggesting that these three miRNAs may be coordinately regulated via common cis-acting elements during metamorphosis. In S2 cells, Ecd and the juvenile hormone analog methoprene exert opposite effects on the expression of these four miRNAs, indicating the participation of both these hormones in the temporal regulation of mir-34, -100, -125, and let-7 expression in vivo.

MicroRNAs (miRNAs) are an extensive class of regulatory RNA whose specific functions in animals are generally unknown. Although computational methods have identified many potential targets of miRNAs, elucidating the spatial expression patterns of miRNAs is necessary to identify the sites of miRNA action. Here, we report the spatial patterns of miRNA transcription during Drosophila embryonic development, as revealed by in situ hybridization to nascent miRNA transcripts. We detect expression of 15 "stand-alone" miRNA loci and 9 intronic miRNA loci, which collectively represent 38 miRNA genes. We observe great variety in the spatial patterns of miRNA transcription, including preblastoderm stripes, in aspects of the central and peripheral nervous systems, and in cellular subsets of the mesoderm and endoderm. We also describe an intronic miRNA (miR-7) whose expression pattern is distinct from that of its host mRNA (bancal), which demonstrates that intronic miRNAs can be subject to independent cis-regulatory control. Intriguingly, the expression patterns of several fly miRNAs are analogous to those of their vertebrate counterparts, suggesting that these miRNAs may have ancient roles in animal patterning.

Drosophila neuroblasts and ovarian stem cells are well characterized models for stem cell biology. In both cell types, one daughter cell self-renews continuously while the other undergoes a limited number of divisions, stops to proliferate mitotically and differentiates. Whereas neuroblasts segregate the Trim-NHL (tripartite motif and Ncl-1, HT2A and Lin-41 domain)-containing protein Brain tumour (Brat) into one of the two daughter cells, ovarian stem cells are regulated by an extracellular signal from the surrounding stem cell niche. After division, one daughter cell looses niche contact. It undergoes 4 transit-amplifying divisions to form a cyst of 16 interconnected cells that reduce their rate of growth and stop to proliferate mitotically. Here we show that the Trim-NHL protein Mei-P26 (refs 7, 8) restricts growth and proliferation in the ovarian stem cell lineage. Mei-P26 expression is low in stem cells but is strongly induced in 16-cell cysts. In mei-P26 mutants, transit-amplifying cells are larger and proliferate indefinitely leading to the formation of an ovarian tumour. Like brat, mei-P26 regulates nucleolar size and can induce differentiation in Drosophila neuroblasts, suggesting that these genes act through the same pathway. We identify Argonaute-1, a component of the RISC complex, as a common binding partner of Brat and Mei-P26, and show that Mei-P26 acts by inhibiting the microRNA pathway. Mei-P26 and Brat have a similar domain composition that is also found in other tumour suppressors and might be a defining property of a new family of microRNA regulators that act specifically in stem cell lineages.

MicroRNAs (miRNAs) are short regulatory RNAs that inhibit target genes by complementary binding in 3' untranslated regions (3' UTRs). They are one of the most abundant classes of regulators, targeting a large fraction of all genes, making their comprehensive study a requirement for understanding regulation and development. Here we use 12 Drosophila genomes to define structural and evolutionary signatures of miRNA hairpins, which we use for their de novo discovery. We predict >41 novel miRNA genes, which encompass many unique families, and 28 of which are validated experimentally. We also define signals for the precise start position of mature miRNAs, which suggest corrections of previously known miRNAs, often leading to drastic changes in their predicted target spectrum. We show that miRNA discovery power scales with the number and divergence of species compared, suggesting that such approaches can be successful in human as dozens of mammalian genomes become available. Interestingly, for some miRNAs sense and anti-sense hairpins score highly and mature miRNAs from both strands can indeed be found in vivo. Similarly, miRNAs with weak 5' end predictions show increased in vivo processing of multiple alternate 5' ends and have fewer predicted targets. Lastly, we show that several miRNA star sequences score highly and are likely functional. For mir-10 in particular, both arms show abundant processing, and both show highly conserved target sites in Hox genes, suggesting a possible cooperation of the two arms, and their role as a master Hox regulator.

MicroRNAs (miRNAs) modulate transcript stability and translation. Functional mature miRNAs are processed from one or both arms of the hairpin precursor. The miR-100/10 family has undergone three independent evolutionary events that have switched the arm from which the functional miRNA is processed. The dominant miR-10 sequences in the insects Drosophila melanogaster and Tribolium castaneum are processed from opposite arms. However, the duplex produced by Dicer cleavage has an identical sequence in fly and beetle. Expression of the Tribolium miR-10 sequence in Drosophila S2 cells recapitulates the native beetle pattern. Thus, arm usage is encoded in the primary miRNA sequence, but outside the mature miRNA duplex. We show that the predicted messenger RNA targets and inferred function of sequences from opposite arms differ significantly. Arm switching is likely to be general, and provides a fundamental mechanism to evolve the function of a miRNA locus and target gene network.

MicroRNAs (miRNAs) are endogenous noncoding small RNAs with important roles in many biological pathways; their generation and activity are under precise regulation [1-3]. Emerging evidence suggests that miRNA pathways are precisely modulated with controls at the level of transcription [4-8], processing [9-11], and stability [12, 13], with miRNA deregulation linked with diseases [14] and neurodegenerative disorders [15]. In the Drosophila miRNA biogenesis pathway, long primary miRNA transcripts undergo sequential cleavage [16-18] to release the embedded miRNAs. Mature miRNAs are then loaded into Argonaute1 (Ago1) within the RNA-induced silencing complex (RISC) [19, 20]. Intriguingly, we found that Drosophila miR-34 displays multiple isoforms that differ at the 3' end, suggesting a novel biogenesis mechanism involving 3' end processing. To define the cellular factors responsible, we performed an RNA interference (RNAi) screen and identified a putative 3'→5' exoribonuclease CG9247/nibbler essential for the generation of the smaller isoforms of miR-34. Nibbler (Nbr) interacts with Ago1 and processes miR-34 within RISC. Deep sequencing analysis revealed a larger set of multi-isoform miRNAs that are controlled by nibbler. These findings suggest that Nbr-mediated 3' end processing represents a critical step in miRNA maturation that impacts miRNA diversity.

microRNAs act in a prevalent and conserved post-transcriptional gene regulatory mechanism that impacts development, homeostasis and disease, yet biological functions for the vast majority of miRNAs remain unknown. Given the power of invertebrate genetics to promote rapid evaluation of miRNA function, recently expanded miRNA identifications (miRBase 10.1), and the importance of assessing potential functional redundancies within and between species, we evaluated miRNA sequence relationships by 5' end match and overall homology criteria to compile a snapshot overview of miRNA families within the C. elegans and D. melanogaster genomes that includes their identified human counterparts. This compilation expands literature documentation of both the number of families and the number of family members, within and between nematode and fly models, and highlights sequences conserved between species pairs or among nematodes, flies and humans. Themes that emerge include the substantial potential for functional redundancy of miRNA sequences within species (84/139 C. elegans miRNAs and 70/152 D. melanogaster miRNAs share significant homology with other miRNAs encoded by their respective genomes), and the striking extent to which miRNAs are conserved across species--over half (73/139) C. elegans miRNAs share sequence homology with miRNAs encoded also in both fly and human genomes. This summary analysis of mature miRNA sequence relationships provides a quickly accessible resource that should facilitate functional and evolutionary analyses of miRNAs and miRNA families.

microRNAs (miRNAs) are endogenous short RNAs that mediate vast networks of post-transcriptional gene regulation. Although computational searches and experimental profiling provide evidence for hundreds of functional targets for individual miRNAs, such data rarely provide clear insight into the phenotypic consequences of manipulating miRNAs in vivo. We describe a genome-wide collection of 165 Drosophila miRNA transgenes and find that a majority induced specific developmental defects, including phenocopies of mutants in myriad cell-signaling and patterning genes. Such connections allowed us to validate several likely targets for miRNA-induced phenotypes. Importantly, few of these phenotypes could be predicted from computationally predicted target lists, thus highlighting the value of whole-animal readouts of miRNA activities. Finally, we provide an example of the relevance of these data to miRNA loss-of-function conditions. Whereas misexpression of several K box miRNAs inhibited Notch pathway activity, reciprocal genetic interaction tests with miRNA sponges demonstrated endogenous roles of the K box miRNA family in restricting Notch signaling. In summary, we provide extensive evidence that misexpression of individual miRNAs often induces specific mutant phenotypes that can guide their functional study. By extension, these data suggest that the deregulation of individual miRNAs in other animals may frequently yield relatively specific phenotypes during disease conditions.

Genetic linkage may result in the expression of multiple products from a polycistronic transcript, under the control of a single promoter. In animals, protein-coding polycistronic transcripts are rare. However, microRNAs are frequently clustered in the genomes of animals, and these clusters are often transcribed as a single unit. The evolution of microRNA clusters has been the subject of much speculation, and a selective advantage of clusters of functionally related microRNAs is often proposed. However, the origin of microRNA clusters has not been so far explored. Here, we study the evolution of microRNA clusters in Drosophila melanogaster. We observed that the majority of microRNA clusters arose by the de novo formation of new microRNA-like hairpins in existing microRNA transcripts. Some clusters also emerged by tandem duplication of a single microRNA. Comparative genomics show that these clusters are unlikely to split or undergo rearrangements. We did not find any instances of clusters appearing by rearrangement of pre-existing microRNA genes. We propose a model for microRNA cluster evolution in which selection over one of the microRNAs in the cluster interferes with the evolution of the other linked microRNAs. Our analysis suggests that the study of microRNAs and small RNAs must consider linkage associations.