Evolution, biogenesis, expression, and target predictions of a substantially expanded set of Drosophila microRNAs

Characterization of small RNAs in Aplysia reveals a role for miR-124 in constraining synaptic plasticity through CREB

Memory storage and memory-related synaptic plasticity rely on precise spatiotemporal regulation of gene expression. To explore the role of small regulatory RNAs in learning-related synaptic plasticity, we carried out massive parallel sequencing to profile the small RNAs of Aplysia californica. We identified 170 distinct miRNAs, 13 of which were novel and specific to Aplysia. Nine miRNAs were brain enriched, and several of these were rapidly downregulated by transient exposure to serotonin, a modulatory neurotransmitter released during learning. Further characterization of the brain-enriched miRNAs revealed that miR-124, the most abundant and well-conserved brain-specific miRNA, was exclusively present presynaptically in a sensory-motor synapse where it constrains serotonin-induced synaptic facilitation through regulation of the transcriptional factor CREB. We therefore present direct evidence that a modulatory neurotransmitter important for learning can regulate the levels of small RNAs and present a role for miR-124 in long-term plasticity of synapses in the mature nervous system.

Processing of Drosophila endo-siRNAs depends on a specific Loquacious isoform

Drosophila melanogaster expresses three classes of small RNAs, which are classified according to their mechanisms of biogenesis. MicroRNAs are approximately 22-23 nucleotides (nt), ubiquitously expressed small RNAs that are sequentially processed from hairpin-like precursors by Drosha/Pasha and Dcr-1/Loquacious complexes. MicroRNAs usually associate with AGO1 and regulate the expression of protein-coding genes. Piwi-interacting RNAs (piRNAs) of approximately 24-28 nt associate with Piwi-family proteins and can arise from single-stranded precursors. piRNAs function in transposon silencing and are mainly restricted to gonadal tissues. Endo-siRNAs are found in both germline and somatic tissues. These approximately 21-nt RNAs are produced by a distinct Dicer, Dcr-2, and do not depend on Drosha/Pasha complexes. They predominantly bind to AGO2 and target both mobile elements and protein-coding genes. Surprisingly, a subset of endo-siRNAs strongly depend for their production on the dsRNA-binding protein Loquacious (Loqs), thought generally to be a partner for Dcr-1 and a cofactor for miRNA biogenesis. Endo-siRNA production depends on a specific Loqs isoform, Loqs-PD, which is distinct from the one, Loqs-PB, required for the production of microRNAs. Paralleling their roles in the biogenesis of distinct small RNA classes, Loqs-PD and Loqs-PB bind to different Dicer proteins, with Dcr-1/Loqs-PB complexes and Dcr-2/Loqs-PD complexes driving microRNA and endo-siRNA biogenesis, respectively.

Repertoire and evolution of miRNA genes in four divergent nematode species

miRNAs are approximately 22-nt RNA molecules that play important roles in post-transcriptional regulation. We have performed small RNA sequencing in the nematodes Caenorhabditis elegans, C. briggsae, C. remanei, and Pristionchus pacificus, which have diverged up to 400 million years ago, to establish the repertoire and evolutionary dynamics of miRNAs in these species. In addition to previously known miRNA genes from C. elegans and C. briggsae we demonstrate expression of many of their homologs in C. remanei and P. pacificus, and identified in total more than 100 novel expressed miRNA genes, the majority of which belong to P. pacificus. Interestingly, more than half of all identified miRNA genes are conserved at the seed level in all four nematode species, whereas only a few miRNAs appear to be species specific. In our compendium of miRNAs we observed evidence for known mechanisms of miRNA evolution including antisense transcription and arm switching, as well as miRNA family expansion through gene duplication. In addition, we identified a novel mode of miRNA evolution, termed "hairpin shifting," in which an alternative hairpin is formed with up- or downstream sequences, leading to shifting of the hairpin and creation of novel miRNA* species. Finally, we identified 21U-RNAs in all four nematodes, including P. pacificus, where the upstream 21U-RNA motif is more diverged. The identification and systematic analysis of small RNA repertoire in four nematode species described here provides a valuable resource for understanding the evolutionary dynamics of miRNA-mediated gene regulation.

Genomic analysis suggests that mRNA destabilization by the microprocessor is specialized for the auto-regulation of Dgcr8

BACKGROUND

The Microprocessor, containing the RNA binding protein Dgcr8 and RNase III enzyme Drosha, is responsible for processing primary microRNAs to precursor microRNAs. The Microprocessor regulates its own levels by cleaving hairpins in the 5'UTR and coding region of the Dgcr8 mRNA, thereby destabilizing the mature transcript.

METHODOLOGY/PRINCIPAL FINDINGS

To determine whether the Microprocessor has a broader role in directly regulating other coding mRNA levels, we integrated results from expression profiling and ultra high-throughput deep sequencing of small RNAs. Expression analysis of mRNAs in wild-type, Dgcr8 knockout, and Dicer knockout mouse embryonic stem (ES) cells uncovered mRNAs that were specifically upregulated in the Dgcr8 null background. A number of these transcripts had evolutionarily conserved predicted hairpin targets for the Microprocessor. However, analysis of deep sequencing data of 18 to 200nt small RNAs in mouse ES, HeLa, and HepG2 indicates that exonic sequence reads that map in a pattern consistent with Microprocessor activity are unique to Dgcr8.

CONCLUSION/SIGNIFICANCE

We conclude that the Microprocessor's role in directly destabilizing coding mRNAs is likely specifically targeted to Dgcr8 itself, suggesting a specialized cellular mechanism for gene auto-regulation.

Sequence features associated with microRNA strand selection in humans and flies

Background

During microRNA (miRNA) maturation in humans and flies, Drosha and Dicer cut the precursor transcript, thereby producing a short RNA duplex. One strand of this duplex becomes a functional component of the RNA-Induced Silencing Complex (RISC), while the other is eliminated. While thermodynamic asymmetry of the duplex ends appears to play a decisive role in the strand selection process, the details of the selection mechanism are not yet understood.

Results

Here, we assess miRNA strand selection bias in humans and fruit flies by analyzing the sequence composition and relative expression levels of the two strands of the precursor duplex in these species. We find that the sequence elements associated with preferential miRNA strand selection and/or rejection differ between the two species. Further, we identify another feature that distinguishes human and fly miRNA processing machinery: the relative accuracy of the Drosha and Dicer enzymes.

Conclusion

Our result provides clues to the mechanistic aspects of miRNA strand selection in humans and other mammals. Further, it indicates that human and fly miRNA processing pathways are more distinct than currently recognized. Finally, the observed strand selection determinants are instrumental in the rational design of efficient miRNA-based expression regulators.

Apprehending multicellularity: regulatory networks, genomics, and evolution

The genomic revolution has provided the first glimpses of the architecture of regulatory networks. Combined with evolutionary information, the "network view" of life processes leads to remarkable insights into how biological systems have been shaped by various forces. This understanding is critical because biological systems, including regulatory networks, are not products of engineering but of historical contingencies. In this light, we attempt a synthetic overview of the natural history of regulatory networks operating in the development and differentiation of multicellular organisms. We first introduce regulatory networks and their organizational principles as can be deduced using ideas from the graph theory. We then discuss findings from comparative genomics to illustrate the effects of lineage-specific expansions, gene-loss, and nonprotein-coding DNA on the architecture of networks. We consider the interaction between expansions of transcription factors, and cis regulatory and more general chromatin state stabilizing elements in the emergence of morphological complexity. Finally, we consider a case study of the Notch subnetwork, which is present throughout Metazoa, to examine how such a regulatory system has been pieced together in evolution from new innovations and pre-existing components that were originally functionally distinct.

A discovery of novel microRNAs in the silkworm (Bombyx mori) genome

MicroRNAs (miRNAs) are pivotal regulators involved in various physiological and pathological processes via their post-transcriptional regulation of gene expressions. We sequenced 14 libraries of small RNAs constructed from samples spanning the life cycle of silkworms, and discovered 50 novel miRNAs previously not known in animals and verified 43 of them using stem-loop RT-PCR. Our genome-wide analyses of 27 species-specific miRNAs suggest they arise from transposable elements, protein-coding genes duplication/transposition and random foldback sequences; which is consistent with the idea that novel animal miRNAs may evolve from incomplete self-complementary transcripts and become fixed in the process of co-adaptation with their targets. Computational prediction suggests that the silkworm-specific miRNAs may have a preference of regulating genes that are related to life-cycle-associated traits, and these genes can serve as potential targets for subsequent studies of the modulating networks in the development of Bombyx mori.

Immunopurification of Ago1 miRNPs selects for a distinct class of microRNA targets

microRNAs comprise a few percent of animal genes and have been recognized as important regulators of a diverse range of biological processes. Understanding the biological functions of miRNAs requires effective means to identify their targets. Combined efforts from computational prediction, miRNA over-expression or depletion, and biochemical purification have identified thousands of potential miRNA-target pairs in cells and organisms. Complementarity to the miRNA seed sequence appears to be a common principle in target recognition. Other features, including miRNA-target duplex stability, binding site accessibility, and local UTR structure might affect target recognition. Yet computational approaches using such contextual features have yielded largely nonoverlapping results and experimental assessment of their impact has been limited. Here, we compare two large sets of miRNA targets: targets identified using an improved Ago1 immunopurification method and targets identified among transcripts up-regulated after Ago1 depletion. We found surprisingly limited overlap between these sets. The two sets showed enrichment for target sites with different molecular, structural and functional properties. Intriguingly, we found a strong correlation between UTR length and other contextual features that distinguish the two groups. This finding was extended to all predicted microRNA targets. Distinct repression mechanisms could have evolved to regulate targets with different contextual features. This study reveals a complex relationship among different features in miRNA-target recognition and poses a new challenge for computational prediction.

Structural determinants of miRNAs for RISC loading and slicer-independent unwinding

MicroRNAs (miRNAs) regulate expression of their target mRNAs through the RNA-induced silencing complex (RISC), which contains an Argonaute (Ago) family protein as a core component. In Drosophila melanogaster, miRNAs are generally sorted into Ago1-containing RISC (Ago1-RISC). We established a native gel system that can biochemically dissect the Ago1-RISC assembly pathway. We found that miRNA-miRNA* duplexes are loaded into Ago1 as double-stranded RNAs in an ATP-dependent fashion. In contrast, unexpectedly, unwinding of miRNA-miRNA* duplexes is a passive process that does not require ATP or slicer activity of Ago1. Central mismatches direct miRNA-miRNA* duplexes into pre-Ago1-RISC, whereas mismatches in the seed or guide strand positions 12-15 promote conversion of pre-Ago1-RISC into mature Ago1-RISC. Our findings show that unwinding of miRNAs is a precise mirror-image process of target recognition, and both processes reflect the unique geometry of RNAs in Ago proteins.

Unlocking the secrets of the genome

Despite the successes of genomics, little is known about how genetic information produces complex organisms. A look at the crucial functional elements of fly and worm genomes could change that.

The National Human Genome Research Institute's modENCODE project (the model organism ENCyclopedia Of DNA Elements) was set up in 2007 with the goal of identifying all the sequence-based functional elements in the genomes of two important experimental organisms, Caenorhabditis elegans and Drosophila melanogaster. Armed with modENCODE data, geneticists will be able to undertake the comprehensive molecular studies of regulatory networks that hold the key to how complex multicellular organisms arise from the list of instructions coded in the genome. In this issue, modENCODE team members outline their plan of campaign. Data from the project are to be made available on http://www.modencode.org and elsewhere as the work progresses.

Evolutionary modeling and prediction of non-coding RNAs in Drosophila

We performed benchmarks of phylogenetic grammar-based ncRNA gene prediction, experimenting with eight different models of structural evolution and two different programs for genome alignment. We evaluated our models using alignments of twelve Drosophila genomes. We find that ncRNA prediction performance can vary greatly between different gene predictors and subfamilies of ncRNA gene. Our estimates for false positive rates are based on simulations which preserve local islands of conservation; using these simulations, we predict a higher rate of false positives than previous computational ncRNA screens have reported. Using one of the tested prediction grammars, we provide an updated set of ncRNA predictions for D. melanogaster and compare them to previously-published predictions and experimental data. Many of our predictions show correlations with protein-coding genes. We found significant depletion of intergenic predictions near the 3' end of coding regions and furthermore depletion of predictions in the first intron of protein-coding genes. Some of our predictions are colocated with larger putative unannotated genes: for example, 17 of our predictions showing homology to the RFAM family snoR28 appear in a tandem array on the X chromosome; the 4.5 Kbp spanned by the predicted tandem array is contained within a FlyBase-annotated cDNA.

Repertoire of bovine miRNA and miRNA-like small regulatory RNAs expressed upon viral infection

MicroRNA (miRNA) and other types of small regulatory RNAs play a crucial role in the regulation of gene expression in eukaryotes. Several distinct classes of small regulatory RNAs have been discovered in recent years. To extend the repertoire of small RNAs characterized in mammals and to examine relationship between host miRNA expression and viral infection we used Illumina's ultrahigh throughput sequencing approach. We sequenced three small RNA libraries prepared from cell line derived from the adult bovine kidney under normal conditions and upon infection of the cell line with Bovine herpesvirus 1. We used a bioinformatics approach to distinguish authentic mature miRNA sequences from other classes of small RNAs and short RNA fragments represented in the sequencing data. Using this approach we detected 219 out of 356 known bovine miRNAs and 115 respective miRNA* sequences. In addition we identified five new bovine orthologs of known mammalian miRNAs and discovered 268 new cow miRNAs many of which are not identifiable in other mammalian genomes and thus might be specific to the ruminant lineage. In addition we found seven new bovine mirtron candidates. We also discovered 10 small nucleolar RNA (snoRNA) loci that give rise to small RNA with possible miRNA-like function. Results presented in this study extend our knowledge of the biology and evolution of small regulatory RNAs in mammals and illuminate mechanisms of small RNA biogenesis and function. New miRNA sequences and the original sequencing data have been submitted to miRNA repository (miRBase) and NCBI GEO archive respectively. We envisage that these resources will facilitate functional annotation of the bovine genome and promote further functional and comparative genomics studies of small regulatory RNA in mammals.

Small RNAs derived from snoRNAs

Small nucleolar RNAs (snoRNAs) guide RNA modification and are localized in nucleoli and Cajal bodies in eukaryotic cells. Components of the RNA silencing pathway associate with these structures, and two recent reports have revealed that a human and a protozoan snoRNA can be processed into miRNA-like RNAs. Here we show that small RNAs with evolutionary conservation of size and position are derived from the vast majority of snoRNA loci in animals (human, mouse, chicken, fruit fly), Arabidopsis, and fission yeast. In animals, sno-derived RNAs (sdRNAs) from H/ACA snoRNAs are predominantly 20-24 nucleotides (nt) in length and originate from the 3' end. Those derived from C/D snoRNAs show a bimodal size distribution at approximately 17-19 nt and >27 nt and predominantly originate from the 5' end. SdRNAs are associated with AGO7 in Arabidopsis and Ago1 in fission yeast with characteristic 5' nucleotide biases and show altered expression patterns in fly loquacious and Dicer-2 and mouse Dicer1 and Dgcr8 mutants. These findings indicate that there is interplay between the RNA silencing and snoRNA-mediated RNA processing systems, and that sdRNAs comprise a novel and ancient class of small RNAs in eukaryotes.

Abundant primary piRNAs, endo-siRNAs, and microRNAs in a Drosophila ovary cell line

Piwi proteins, a subclass of Argonaute-family proteins, carry approximately 24-30-nt Piwi-interacting RNAs (piRNAs) that mediate gonadal defense against transposable elements (TEs). We analyzed the Drosophila ovary somatic sheet (OSS) cell line and found that it expresses miRNAs, endogenous small interfering RNAs (endo-siRNAs), and piRNAs in abundance. In contrast to intact gonads, which contain mixtures of germline and somatic cell types that express different Piwi-class proteins, OSS cells are a homogenous somatic cell population that expresses only PIWI and primary piRNAs. Detailed examination of its TE-derived piRNAs and endo-siRNAs revealed aspects of TE defense that do not rely upon ping-pong amplification. In particular, we provide evidence that a subset of piRNA master clusters, including flamenco, are specifically expressed in OSS and ovarian follicle cells. These data indicate that the restriction of certain TEs in somatic gonadal cells is largely mediated by a primary piRNA pathway.

Identification and analysis of miRNAs in human breast cancer and teratoma samples using deep sequencing

Background

MiRNAs play important roles in cellular control and in various disease states such as cancers, where they may serve as markers or possibly even therapeutics. Identifying the whole repertoire of miRNAs and understanding their expression patterns is therefore an important goal.

Methods

Here we describe the analysis of 454 pyrosequencing of small RNA from four different tissues: Breast cancer, normal adjacent breast, and two teratoma cell lines. We developed a pipeline for identifying new miRNAs, emphasizing extracting and retaining as much data as possible from even noisy sequencing data. We investigated differential expression of miRNAs in the breast cancer and normal adjacent breast samples, and systematically examined the mature sequence end variability of miRNA compared to non-miRNA loci.

Results

We identified five novel miRNAs, as well as two putative alternative precursors for known miRNAs. Several miRNAs were differentially expressed between the breast cancer and normal breast samples. The end variability was shown to be significantly different between miRNA and non-miRNA loci.

Conclusion

Pyrosequencing of small RNAs, together with a computational pipeline, can be used to identify miRNAs in tumor and other tissues. Measures of miRNA end variability may in the future be incorporated into the discovery pipeline as a discriminatory feature. Breast cancer samples show a distinct miRNA expression profile compared to normal adjacent breast.

miR-10 in development and cancer

The microRNA (miRNA) miR-10 family has attracted attention because of its conservation and the position of the miR-10 genes within the Hox clusters of developmental regulators. In several species, miR-10 is coexpressed with a set of Hox genes and has been found to regulate the translation of Hox transcripts. In addition, members of the miR-10 family are de-regulated in several cancer forms. Aside from acting in translational repression, miR-10 was recently found to bind a group of transcripts containing a terminal oligo-pyrimidine (TOP) motif and to induce their translation, thereby adding a new function to the miRNA repertoire.

Repertoire of microRNAs in epithelial ovarian cancer as determined by next generation sequencing of small RNA cDNA libraries

Background

MicroRNAs (miRNAs) are small regulatory RNAs that are implicated in cancer pathogenesis and have recently shown promise as blood-based biomarkers for cancer detection. Epithelial ovarian cancer is a deadly disease for which improved outcomes could be achieved by successful early detection and enhanced understanding of molecular pathogenesis that leads to improved therapies. A critical step toward these goals is to establish a comprehensive view of miRNAs expressed in epithelial ovarian cancer tissues as well as in normal ovarian surface epithelial cells.

Methodology

We used massively parallel pyrosequencing (i.e., “454 sequencing”) to discover and characterize novel and known miRNAs expressed in primary cultures of normal human ovarian surface epithelium (HOSE) and in tissue from three of the most common histotypes of ovarian cancer. Deep sequencing of small RNA cDNA libraries derived from normal HOSE and ovarian cancer samples yielded a total of 738,710 high-quality sequence reads, generating comprehensive digital profiles of miRNA expression. Expression profiles for 498 previously annotated miRNAs were delineated and we discovered six novel miRNAs and 39 candidate miRNAs. A set of 124 miRNAs was differentially expressed in normal versus cancer samples and 38 miRNAs were differentially expressed across histologic subtypes of ovarian cancer. Taqman qRT-PCR performed on a subset of miRNAs confirmed results of the sequencing-based study.

Conclusions

This report expands the body of miRNAs known to be expressed in epithelial ovarian cancer and provides a useful resource for future studies of the role of miRNAs in the pathogenesis and early detection of ovarian cancer.

Tiny RNAs associated with transcription start sites in animals

It has been reported that relatively short RNAs of heterogeneous sizes are derived from sequences near the promoters of eukaryotic genes. In conjunction with the FANTOM4 project, we have identified tiny RNAs with a modal length of 18 nt that map within -60 to +120 nt of transcription start sites (TSSs) in human, chicken and Drosophila. These transcription initiation RNAs (tiRNAs) are derived from sequences on the same strand as the TSS and are preferentially associated with G+C-rich promoters. The 5' ends of tiRNAs show peak density 10-30 nt downstream of TSSs, indicating that they are processed. tiRNAs are generally, although not exclusively, associated with highly expressed transcripts and sites of RNA polymerase II binding. We suggest that tiRNAs may be a general feature of transcription in metazoa and possibly all eukaryotes.

Identification and characterization of new miRNAs cloned from normal mouse mammary gland

BACKGROUND

MicroRNAs (miRNAs) are small non-coding RNAs that have been found to play important roles in silencing target genes and that are involved in the regulation of various normal cellular processes. Until now their implication in the mammary gland biology was suggested by few studies mainly focusing on pathological situations allowing the characterization of miRNAs as markers of breast cancer tumour classes. If in the normal mammary gland, the expression of known miRNAs has been studied in human and mice but the full repertoire of miRNAs expressed in this tissue is not yet available.

RESULTS

To extend the repertoire of mouse mammary gland expressed miRNAs, we have constructed several libraries of small miRNAs allowing the cloning of 455 sequences. After bioinformatics' analysis, 3 known miRNA (present in miRbase) and 33 new miRNAs were identified. Expression of 24 out of the 33 has been confirmed by RT-PCR. Expression of none of them was found to be mammary specific, despite a tissue-restricted distribution of some of them. No correlation could be established between their expression pattern and evolutionary conservation. Six of them appear to be mouse specific. In several cases, multiple potential precursors of miRNA were present in the genome and we have developed a strategy to determine which of them was able to mature the miRNA.

CONCLUSION

The cloning approach has allowed improving the repertoire of miRNAs in the mammary gland, an evolutionary recent organ. This tissue is a good candidate to find tissue-specific miRNAs and to detect miRNA specific to mammals. We provide evidence for 24 new miRNA. If none of them is mammary gland specific, a few of them are not ubiquitously expressed. For the first time 6 mouse specific miRNA have been identified.

The interplay between transcription factors and microRNAs in genome-scale regulatory networks

Metazoan genomes contain thousands of protein-coding and non-coding RNA genes, most of which are differentially expressed, i.e., at different locations, at different times during development, or in response to environmental signals. Differential gene expression is achieved through complex regulatory networks that are controlled in part by two types of trans-regulators: transcription factors (TFs) and microRNAs (miRNAs). TFs bind to cis-regulatory DNA elements that are often located in or near their target genes, while miRNAs hybridize to cis-regulatory RNA elements mostly located in the 3' untranslated region of their target mRNAs. Here, we describe how these trans-regulators interact with each other in the context of gene regulatory networks to coordinate gene expression at the genome-scale level, and discuss future challenges of integrating these networks with other types of functional networks.

Genome-wide identification of targets of the drosha-pasha/DGCR8 complex

Drosha is a type III RNase, which plays a critical role in miRNA biogenesis. Drosha and its double-stranded RNA-binding partner protein Pasha/DGCR8 likely recognize and cleave miRNA precursor RNAs or pri-miRNA hairpins cotranscriptionally. To identify RNAs processed by Drosha, we used tiling microarrays to examine transcripts after depletion of drosha mRNA with dsRNA in Drosophila Schneider S2 cells. This strategy identified 137 Drosha-regulated RNAs, including 11 putative pri-miRNAs comprising 15 annotated miRNAs. Most of the identified pri-miRNAs seem extremely large, >10 kb as revealed by both the Drosha knock-down strategy and by RNA PolII chromatin IP followed by Drosophila tiling microarrays. Surprisingly, more than a hundred additional RNAs not annotated as miRNAs are under Drosha control and are likely to be direct targets of Drosha action. This is because many of them encode annotated genes, and unlike bona fide pri-miRNAs, they are not affected by depletion of the miRNA processing factor, dicer-1. Moreover, application of the evofold analysis software indicates that at least 25 of the Drosha-regulated RNAs contain evolutionarily conserved hairpins similar to those recognized by the Drosha-Pasha/DGCR8 complex in pri-miRNAs. One of these hairpins is located in the 5' UTR of both pasha and mammalian DGCR8. These observations suggest that a negative feedback loop acting on pasha mRNA may regulate the miRNA-biogenesis pathway: i.e., excess Drosha cleaves pasha/DGCR8 primary transcripts and leads to a reduction in pasha/DGCR8 mRNA levels and Pasha/DGCR8 synthesis.

Dicer-1-dependent Dacapo suppression acts downstream of Insulin receptor in regulating cell division of Drosophila germline stem cells

It is important to understand the regulation of stem cell division because defects in this process can cause altered tissue homeostasis or cancer. The cyclin-dependent kinase inhibitor Dacapo (Dap), a p21/p27 homolog, acts downstream of the microRNA (miRNA) pathway to regulate the cell cycle in Drosophila melanogaster germline stem cells (GSCs). Tissue-extrinsic signals, including insulin, also regulate cell division of GSCs. We report that intrinsic and extrinsic regulators intersect in GSC division control; the Insulin receptor (InR) pathway regulates Dap levels through miRNAs, thereby controlling GSC division. Using GFP-dap 3'UTR sensors in vivo, we show that in GSCs the dap 3'UTR is responsive to Dicer-1, an RNA endonuclease III required for miRNA processing. Furthermore, the dap 3'UTR can be directly targeted by miR-7, miR-278 and miR-309 in luciferase assays. Consistent with this, miR-278 and miR-7 mutant GSCs are partially defective in GSC division and show abnormal cell cycle marker expression, respectively. These data suggest that the GSC cell cycle is regulated via the dap 3'UTR by multiple miRNAs. Furthermore, the GFP-dap 3'UTR sensors respond to InR but not to TGF-beta signaling, suggesting that InR signaling utilizes Dap for GSC cell cycle regulation. We further demonstrate that the miRNA-based Dap regulation may act downstream of InR signaling; Dcr-1 and Dap are required for nutrition-dependent cell cycle regulation in GSCs and reduction of dap partially rescues the cell cycle defect of InR-deficient GSCs. These data suggest that miRNA- and Dap-based cell cycle regulation in GSCs can be controlled by InR signaling.

Current tools for the identification of miRNA genes and their targets

The discovery of microRNAs (miRNAs), almost 10 years ago, changed dramatically our perspective on eukaryotic gene expression regulation. However, the broad and important functions of these regulators are only now becoming apparent. The expansion of our catalogue of miRNA genes and the identification of the genes they regulate owe much to the development of sophisticated computational tools that have helped either to focus or interpret experimental assays. In this article, we review the methods for miRNA gene finding and target identification that have been proposed in the last few years. We identify some problems that current approaches have not yet been able to overcome and we offer some perspectives on the next generation of computational methods.

Small silencing RNAs: an expanding universe

Since the discovery in 1993 of the first small silencing RNA, a dizzying number of small RNA classes have been identified, including microRNAs (miRNAs), small interfering RNAs (siRNAs) and Piwi-interacting RNAs (piRNAs). These classes differ in their biogenesis, their modes of target regulation and in the biological pathways they regulate. There is a growing realization that, despite their differences, these distinct small RNA pathways are interconnected, and that small RNA pathways compete and collaborate as they regulate genes and protect the genome from external and internal threats.

MicroRNAs: target recognition and regulatory functions

MicroRNAs (miRNAs) are endogenous approximately 23 nt RNAs that play important gene-regulatory roles in animals and plants by pairing to the mRNAs of protein-coding genes to direct their posttranscriptional repression. This review outlines the current understanding of miRNA target recognition in animals and discusses the widespread impact of miRNAs on both the expression and evolution of protein-coding genes.

MicroRNAs: target recognition and regulatory functions

MicroRNAs (miRNAs) are endogenous approximately 23 nt RNAs that play important gene-regulatory roles in animals and plants by pairing to the mRNAs of protein-coding genes to direct their posttranscriptional repression. This review outlines the current understanding of miRNA target recognition in animals and discusses the widespread impact of miRNAs on both the expression and evolution of protein-coding genes.

MicroRNAs: target recognition and regulatory functions

MicroRNAs (miRNAs) are endogenous approximately 23 nt RNAs that play important gene-regulatory roles in animals and plants by pairing to the mRNAs of protein-coding genes to direct their posttranscriptional repression. This review outlines the current understanding of miRNA target recognition in animals and discusses the widespread impact of miRNAs on both the expression and evolution of protein-coding genes.

A distinct class of small RNAs arises from pre-miRNA-proximal regions in a simple chordate

MicroRNAs (miRNAs) have been implicated in various cellular processes. They are thought to function primarily as inhibitors of gene activity by attenuating translation or promoting mRNA degradation. A typical miRNA gene produces a predominant approximately 21-nucleotide (nt) RNA (the miRNA) along with a less abundant miRNA(*) product. We sought to identify miRNAs from the simple chordate Ciona intestinalis through comprehensive sequencing of small RNA libraries created from different developmental stages. Unexpectedly, half of the identified miRNA loci encode up to four distinct, stable small RNAs. The additional RNAs, miRNA-offset RNAs (moRs), are generated from sequences immediately adjacent to the predicted approximately 60-nt pre-miRNA. moRs seem to be produced by RNAse III-like processing, are approximately 20 nt long and, like miRNAs, are observed at specific developmental stages. We present evidence suggesting that the biogenesis of moRs results from an intrinsic property of the miRNA processing machinery in C. intestinalis.

Characterization and comparative profiling of the small RNA transcriptomes in two phases of locust

High-throughput sequencing of the small RNA transcriptome of locust reveals differences in post-transcriptional regulation between solitary and swarming phases and provides insights into the evolution of insect small RNAs.

Background

All the reports on insect small RNAs come from holometabolous insects whose genome sequence data are available. Therefore, study of hemimetabolous insect small RNAs could provide more insights into evolution and function of small RNAs in insects. The locust is an important, economically harmful hemimetabolous insect. Its phase changes, as a phenotypic plasticity, result from differential gene expression potentially regulated at both the post-transcriptional level, mediated by small RNAs, and the transcriptional level.

Results

Here, using high-throughput sequencing, we characterize the small RNA transcriptome in the locust. We identified 50 conserved microRNA families by similarity searching against miRBase, and a maximum of 185 potential locust-specific microRNA family candidates were identified using our newly developed method independent of locust genome sequence. We also demonstrate conservation of microRNA*, and evolutionary analysis of locust microRNAs indicates that the generation of miRNAs in locusts is concentrated along three phylogenetic tree branches: bilaterians, coelomates, and insects. Our study identified thousands of endogenous small interfering RNAs, some of which were of transposon origin, and also detected many Piwi-interacting RNA-like small RNAs. Comparison of small RNA expression patterns of the two phases showed that longer small RNAs were expressed more abundantly in the solitary phase and that each category of small RNAs exhibited different expression profiles between the two phases.

Conclusions

The abundance of small RNAs in the locust might indicate a long evolutionary history of post-transcriptional gene expression regulation, and differential expression of small RNAs between the two phases might further disclose the molecular mechanism of phase changes.

The evolution of RNAi as a defence against viruses and transposable elements

RNA interference (RNAi) is an important defence against viruses and transposable elements (TEs). RNAi not only protects against viruses by degrading viral RNA, but hosts and viruses can also use RNAi to manipulate each other's gene expression, and hosts can encode microRNAs that target viral sequences. In response, viruses have evolved a myriad of adaptations to suppress and evade RNAi. RNAi can also protect cells against TEs, both by degrading TE transcripts and by preventing TE expression through heterochromatin formation. The aim of our review is to summarize and evaluate the current data on the evolution of these RNAi defence mechanisms. To this end, we also extend a previous analysis of the evolution of genes of the RNAi pathways. Strikingly, we find that antiviral RNAi genes, anti-TE RNAi genes and viral suppressors of RNAi all evolve rapidly, suggestive of an evolutionary arms race between hosts and parasites. Over longer time scales, key RNAi genes are repeatedly duplicated or lost across the metazoan phylogeny, with important implications for RNAi as an immune defence.

The evolution and functional diversification of animal microRNA genes

microRNAs (miRNAs) are an abundant class of approximately 22 nucleotide (nt) regulatory RNAs that are pervasive in higher eukaryotic genomes. In order to fully understand their prominence in genomes, it is necessary to elucidate the molecular mechanisms that can diversify miRNA activities. In this review, we describe some of the many strategies that allow novel miRNA functions to emerge, with particular emphasis on how miRNA genes evolve in animals. These mechanisms include changes in their sequence, processing, or expression pattern; acquisition of miRNA* functionality or antisense processing; and de novo gene birth. The facility and versatility of miRNAs to evolve and change likely underlies how they have become dominant constituents of higher genomes.

A double-edged sword to force posterior dominance of Hox genes

Spatially and temporally restricted expression of Hox genes requires multiple mechanisms at both the transcriptional and the post-transcriptional levels. New insight into this precise expression mechanism comes from recent findings of a novel sense-antisense miRNA combination from the bithorax complex of Drosophila melanogaster. These two miRNAs encoded from the same locus target 3' untranslated regions of anterior hox genes, Antp, Ubx and abd-A to establish the dominance of posterior hox gene Abd-B in its expression domain. Such double-edge tools, sense-antisense miRNA combinations, also operate at multiple loci in the genome implicating their wider impact on the post-transcriptional gene regulation in eukaryotes.

Mouse ES cells express endogenous shRNAs, siRNAs, and other Microprocessor-independent, Dicer-dependent small RNAs

Canonical microRNAs (miRNAs) require two processing steps: the first by the Microprocessor, a complex of DGCR8 and Drosha, and the second by a complex of TRBP and Dicer. dgcr8Delta/Delta mouse embryonic stem cells (mESCs) have less severe phenotypes than dicer1Delta/Delta mESCs, suggesting a physiological role for Microprocessor-independent, Dicer-dependent small RNAs. To identify these small RNAs with unusual biogenesis, we performed high-throughput sequencing from wild-type, dgcr8Delta/Delta, and dicer1Delta/Delta mESCs. Several of the resulting DGCR8-independent, Dicer-dependent RNAs were noncanonical miRNAs. These derived from mirtrons and a newly identified subclass of miRNA precursors, which appears to be the endogenous counterpart of shRNAs. Our analyses also revealed endogenous siRNAs resulting from Dicer cleavage of long hairpins, the vast majority of which originated from one genomic locus with tandem, inverted short interspersed nuclear elements (SINEs). Our results extend the known diversity of mammalian small RNA-generating pathways and show that mammalian siRNAs exist in cell types other than oocytes.

Biological principles of microRNA-mediated regulation: shared themes amid diversity

Regulation of gene activity by microRNAs is critical to myriad aspects of eukaryotic development and physiology. Amidst an extensive regulatory web that is predicted to involve thousands of transcripts, emergent themes are now beginning to illustrate how microRNAs have been incorporated into diverse settings. These include potent inhibition of individual key targets, fine-tuning of target activity, the coordinated regulation of target batteries, and the reversibility of some aspects of microRNA-mediated repression. Such themes may reflect some of the inherent advantages of exploiting microRNA control in biological circuits, and provide insight into the consequences of microRNA dysfunction in disease.

Genome-scale spatiotemporal analysis of Caenorhabditis elegans microRNA promoter activity

The Caenorhabditis elegans genome encodes more than 100 microRNAs (miRNAs). Genetic analyses of miRNA deletion mutants have only provided limited insights into miRNA function. To gain insight into the function of miRNAs, it is important to determine their spatiotemporal expression pattern. Here, we use miRNA promoters driving the expression of GFP as a proxy for miRNA expression. We describe a set of 73 transgenic C. elegans strains, each expressing GFP under the control of a miRNA promoter. Together, these promoters control the expression of 89 miRNAs (66% of all predicted miRNAs). We find that miRNA promoters drive GFP expression in a variety of tissues and that, overall, their activity is similar to that of protein-coding gene promoters. However, miRNAs are expressed later in development, which is consistent with functions after initial body-plan specification. We find that miRNA members belonging to families are more likely to be expressed in overlapping tissues than miRNAs that do not belong to the same family, and provide evidence that intronic miRNAs may be controlled by their own, rather than a host gene promoter. Finally, our data suggest that post-transcriptional mechanisms contribute to differential miRNA expression. The data and strains described here will provide a valuable guide and resource for the functional analysis of C. elegans miRNAs.

Genome-wide identification and analysis of small RNAs originated from natural antisense transcripts in Oryza sativa

Natural antisense transcripts (NATs) have been shown to play important roles in post-transcriptional regulation through the RNA interference pathway. We have combined pyrophosphate-based high-throughput sequencing and computational analysis to identify and analyze, in genome scale, cis-NAT and trans-NAT small RNAs that are derived under normal conditions and in response to drought and salt stresses in the staple plant Oryza sativa. Computationally, we identified 344 cis-NATs and 7,142 trans-NATs that are formed by protein-coding genes. From the deep sequencing data, we found 108 cis-NATs and 7,141 trans-NATs that gave rise to small RNAs from their overlapping regions. Consistent with early findings, the majority of these 108 cis-NATs seem to be associated with specific conditions or developmental stages. Our analyses also revealed several interesting results. The overlapping regions of the cis-NATs and trans-NATs appear to be more enriched with small RNA loci than non-overlapping regions. The small RNAs generated from cis-NATs and trans-NATs have a length bias of 21 nt, even though their lengths spread over a large range. Furthermore, >40% of the small RNAs from cis-NATs and trans-NATs carry an A as their 5'-terminal nucleotides. A substantial portion of the transcripts are involved in both cis-NATs and trans-NATs, and many trans-NATs can form many-to-many relationships, indicating that NATs may form complex regulatory networks in O. sativa. This study is the first genome-wide investigation of NAT-derived small RNAs in O. sativa. It reveals the importance of NATs in biogenesis of small RNAs and broadens our understanding of the roles of NAT-derived small RNAs in gene regulation, particularly in response to environmental stimuli.

Most mammalian mRNAs are conserved targets of microRNAs

MicroRNAs (miRNAs) are small endogenous RNAs that pair to sites in mRNAs to direct post-transcriptional repression. Many sites that match the miRNA seed (nucleotides 2-7), particularly those in 3' untranslated regions (3'UTRs), are preferentially conserved. Here, we overhauled our tool for finding preferential conservation of sequence motifs and applied it to the analysis of human 3'UTRs, increasing by nearly threefold the detected number of preferentially conserved miRNA target sites. The new tool more efficiently incorporates new genomes and more completely controls for background conservation by accounting for mutational biases, dinucleotide conservation rates, and the conservation rates of individual UTRs. The improved background model enabled preferential conservation of a new site type, the "offset 6mer," to be detected. In total, >45,000 miRNA target sites within human 3'UTRs are conserved above background levels, and >60% of human protein-coding genes have been under selective pressure to maintain pairing to miRNAs. Mammalian-specific miRNAs have far fewer conserved targets than do the more broadly conserved miRNAs, even when considering only more recently emerged targets. Although pairing to the 3' end of miRNAs can compensate for seed mismatches, this class of sites constitutes less than 2% of all preferentially conserved sites detected. The new tool enables statistically powerful analysis of individual miRNA target sites, with the probability of preferentially conserved targeting (P(CT)) correlating with experimental measurements of repression. Our expanded set of target predictions (including conserved 3'-compensatory sites), are available at the TargetScan website, which displays the P(CT) for each site and each predicted target.

Warts is required for PI3K-regulated growth arrest, autophagy, and autophagic cell death in Drosophila

BACKGROUND

Cell growth arrest and autophagy are required for autophagic cell death in Drosophila. Maintenance of growth by expression of either activated Ras, Dp110, or Akt is sufficient to inhibit autophagy and cell death in Drosophila salivary glands, but the mechanism that controls growth arrest is unknown. Although the Warts (Wts) tumor suppressor is a critical regulator of tissue growth in animals, it is not clear how this signaling pathway controls cell growth.

RESULTS

Here, we show that genes in the Wts pathway are required for salivary gland degradation and that wts mutants have defects in cell growth arrest, caspase activity, and autophagy. Expression of Atg1, a regulator of autophagy, in salivary glands is sufficient to rescue wts mutant salivary gland destruction. Surprisingly, expression of Yorkie (Yki) and Scalloped (Sd) in salivary glands fails to phenocopy wts mutants. By contrast, misexpression of the Yki target bantam was able to inhibit salivary gland cell death, even though mutations in bantam fail to suppress the wts mutant salivary gland-persistence phenotype. Significantly, wts mutant salivary glands possess altered phosphoinositide signaling, and decreased function of the class I PI3K-pathway genes chico and TOR suppressed wts defects in cell death.

CONCLUSIONS

Although we have previously shown that salivary gland degradation requires genes in the Wts pathway, this study provides the first evidence that Wts influences autophagy. Our data indicate that the Wts-pathway components Yki, Sd, and bantam fail to function in salivary glands and that Wts regulates salivary gland cell death in a PI3K-dependent manner.

MicroRNAs in the Hox network: an apparent link to posterior prevalence

Homeobox (Hox) transcription factors confer anterior-posterior (AP) axial coordinates to vertebrate embryos. Hox genes are found in clusters that also contain genes for microRNAs (miRNAs). Our analysis of predicted miRNA targets indicates that Hox cluster-embedded miRNAs preferentially target Hox mRNAs. Moreover, the presumed Hox target genes are predominantly situated on the 3' side of each Hox miRNA locus. These results suggest that Hox miRNAs help repress more anterior programmes, thereby reinforcing posterior prevalence, which is the hierarchical dominance of posterior over anterior Hox gene function that is observed in bilaterians. In this way, miRNA-mediated regulation seems to recapitulate interactions at other levels of gene expression, some more ancestral, within a network under stabilizing selection.

Early origins and evolution of microRNAs and Piwi-interacting RNAs in animals

In bilaterian animals, such as humans, flies and worms, hundreds of microRNAs (miRNAs), some conserved throughout bilaterian evolution, collectively regulate a substantial fraction of the transcriptome. In addition to miRNAs, other bilaterian small RNAs, known as Piwi-interacting RNAs (piRNAs), protect the genome from transposons. Here we identify small RNAs from animal phyla that diverged before the emergence of the Bilateria. The cnidarian Nematostella vectensis (starlet sea anemone), a close relative to the Bilateria, possesses an extensive repertoire of miRNA genes, two classes of piRNAs and a complement of proteins specific to small-RNA biology comparable to that of humans. The poriferan Amphimedon queenslandica (sponge), one of the simplest animals and a distant relative of the Bilateria, also possesses miRNAs, both classes of piRNAs and a full complement of the small-RNA machinery. Animal miRNA evolution seems to have been relatively dynamic, with precursor sizes and mature miRNA sequences differing greatly between poriferans, cnidarians and bilaterians. Nonetheless, miRNAs and piRNAs have been available as classes of riboregulators to shape gene expression throughout the evolution and radiation of animal phyla.

Specific alignment of structured RNA: stochastic grammars and sequence annealing

MOTIVATION

Whole-genome screens suggest that eukaryotic genomes are dense with non-coding RNAs (ncRNAs). We introduce a novel approach to RNA multiple alignment which couples a generative probabilistic model of sequence and structure with an efficient sequence annealing approach for exploring the space of multiple alignments. This leads to a new software program, Stemloc-AMA, that is both accurate and specific in the alignment of multiple related RNA sequences.

RESULTS

When tested on the benchmark datasets BRalibase II and BRalibase 2.1, Stemloc-AMA has comparable sensitivity to and better specificity than the best competing methods. We use a large-scale random sequence experiment to show that while most alignment programs maximize sensitivity at the expense of specificity, even to the point of giving complete alignments of non-homologous sequences, Stemloc-AMA aligns only sequences with detectable homology and leaves unrelated sequences largely unaligned. Such accurate and specific alignments are crucial for comparative-genomics analysis, from inferring phylogeny to estimating substitution rates across different lineages.

AVAILABILITY

Stemloc-AMA is available from http://biowiki.org/StemLocAMA as part of the dart software package for sequence analysis.

Endogenous small interfering RNAs in animals

Until recently, only nematodes among animals had a well-defined endogenous small interfering RNA (endo-siRNA) pathway. This has changed dramatically with the recent discovery of diverse intramolecular and intermolecular substrates that generate endo-siRNAs in Drosophila melanogaster and mice. These findings suggest broad and possibly conserved roles for endogenous RNA interference in regulating host-gene expression and transposable element transcripts. They also raise many questions regarding the biogenesis and function of small regulatory RNAs in animals.

The long and short of inverted repeat genes in animals: microRNAs, mirtrons and hairpin RNAs

MicroRNAs (miRNAs) are endogenous transcripts that contain intramolecular double stranded RNA (dsRNA) and are processed by Dicer. Their mature products are approximately 21-24 nucleotides in length, and they collectively regulate a broad network of endogenous transcripts. A subset of animal miRNAs are produced from mirtrons, short hairpin introns whose splicing bypasses the normal nuclear processing of canonical miRNAs. Recent studies revealed novel, extended intramolecular dsRNA produced by defined transcription units in flies and mammals, termed hairpin RNAs (hpRNAs). Detailed biogenesis studies in Drosophila showed that hpRNAs are not merely "long" miRNAs, but are actually processed by a distinct biogenesis pathway that is related to the canonical RNA interference pathway. We compare and contrast the miRNA and hpRNA pathways in this review, and describe some of the key questions that the recognition of this novel pathway raises.

The silkworm (Bombyx mori) microRNAs and their expressions in multiple developmental stages

Background

MicroRNAs (miRNAs) play crucial roles in various physiological processes through post-transcriptional regulation of gene expressions and are involved in development, metabolism, and many other important molecular mechanisms and cellular processes. The Bombyx mori genome sequence provides opportunities for a thorough survey for miRNAs as well as comparative analyses with other sequenced insect species.

Methodology/Principal Findings

We identified 114 non-redundant conserved miRNAs and 148 novel putative miRNAs from the B. mori genome with an elaborate computational protocol. We also sequenced 6,720 clones from 14 developmental stage-specific small RNA libraries in which we identified 35 unique miRNAs containing 21 conserved miRNAs (including 17 predicted miRNAs) and 14 novel miRNAs (including 11 predicted novel miRNAs). Among the 114 conserved miRNAs, we found six pairs of clusters evolutionarily conserved cross insect lineages. Our observations on length heterogeneity at 5′ and/or 3′ ends of nine miRNAs between cloned and predicted sequences, and three mature forms deriving from the same arm of putative pre-miRNAs suggest a mechanism by which miRNAs gain new functions. Analyzing development-related miRNAs expression at 14 developmental stages based on clone-sampling and stem-loop RT PCR, we discovered an unusual abundance of 33 sequences representing 12 different miRNAs and sharply fluctuated expression of miRNAs at larva-molting stage. The potential functions of several stage-biased miRNAs were also analyzed in combination with predicted target genes and silkworm's phenotypic traits; our results indicated that miRNAs may play key regulatory roles in specific developmental stages in the silkworm, such as ecdysis.

Conclusions/Significance

Taking a combined approach, we identified 118 conserved miRNAs and 151 novel miRNA candidates from the B. mori genome sequence. Our expression analyses by sampling miRNAs and real-time PCR over multiple developmental stages allowed us to pinpoint molting stages as hotspots of miRNA expression both in sorts and quantities. Based on the analysis of target genes, we hypothesized that miRNAs regulate development through a particular emphasis on complex stages rather than general regulatory mechanisms.

MicroRNAs in development and disease

Since the discovery of microRNAs (miRNAs) in Caenorhabditis elegans, mounting evidence illustrates the important regulatory roles for miRNAs in various developmental, differentiation, cell proliferation, and apoptosis pathways of diverse organisms. We are just beginning to elucidate novel aspects of RNA mediated gene regulation and to understand how heavily various molecular pathways rely on miRNAs for their normal function. miRNAs are small non-protein-coding transcripts that regulate gene expression post-transcriptionally by targeting messenger RNAs (mRNAs). While individual miRNAs have been specifically linked to critical developmental pathways, the deregulated expression of many miRNAs also has been shown to have functional significance for multiple human diseases, such as cancer. We continue to discover novel functional roles for miRNAs at a rapid pace. Here, we summarize some of the key recent findings on miRNAs, their mode of action, and their roles in both normal development and in human pathology. A better understanding of how miRNAs operate during the normal life of a cell as well as in the pathogenesis of disease when deregulated will provide new avenues for diagnostic, prognostic, and therapeutic applications.

Lessons from microRNA mutants in worms, flies and mice

It is apparent that microRNAs (miRNAs) are important components in the regulation of genetic networks in many biological contexts. Based on computational analysis, typical miRNAs are inferred to have tens to hundreds of conserved targets. Many miRNA-target interactions have been validated by various means, including heterologous tests in cultured cells and gain-of-function approaches that can yield striking phenotypes in whole animals. However, these strategies do not report on the endogenous importance of such miRNA activities. Likewise, studies of miRNA pathway mutants can suggest an endogenous role for miRNAs in a given setting, but do not identify roles for specific miRNAs. Therefore, these approaches must be complemented with the analysis of miRNA mutant alleles. In this review, we describe some of the lessons learned from studying miRNA gene deletions in worms, flies and mice, and discuss their implications for the control of endogenous regulatory networks.

Genomic analysis of the relationship between gene expression variation and DNA polymorphism in Drosophila simulans

BACKGROUND

Understanding how DNA sequence polymorphism relates to variation in gene expression is essential to connecting genotypic differences with phenotypic differences among individuals. Addressing this question requires linking population genomic data with gene expression variation.

RESULTS

Using whole genome expression data and recent light shotgun genome sequencing of six Drosophila simulans genotypes, we assessed the relationship between expression variation in males and females and nucleotide polymorphism across thousands of loci. By examining sequence polymorphism in gene features, such as untranslated regions and introns, we find that genes showing greater variation in gene expression between genotypes also have higher levels of sequence polymorphism in many gene features. Accordingly, X-linked genes, which have lower sequence polymorphism levels than autosomal genes, also show less expression variation than autosomal genes. We also find that sex-specifically expressed genes show higher local levels of polymorphism and divergence than both sex-biased and unbiased genes, and that they appear to have simpler regulatory regions.

CONCLUSION

The gene-feature-based analyses and the X-to-autosome comparisons suggest that sequence polymorphism in cis-acting elements is an important determinant of expression variation. However, this relationship varies among the different categories of sex-biased expression, and trans factors might contribute more to male-specific gene expression than cis effects. Our analysis of sex-specific gene expression also shows that female-specific genes have been overlooked in analyses that only point to male-biased genes as having unusual patterns of evolution and that studies of sexually dimorphic traits need to recognize that the relationship between genetic and expression variation at these traits is different from the genome as a whole.

Sequence relationships among C. elegans, D. melanogaster and human microRNAs highlight the extensive conservation of microRNAs in biology

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.

MicroRNA discovery and profiling in human embryonic stem cells by deep sequencing of small RNA libraries

We used massively parallel pyrosequencing to discover and characterize microRNAs (miRNAs) expressed in human embryonic stem cells (hESC). Sequencing of small RNA cDNA libraries derived from undifferentiated hESC and from isogenic differentiating cultures yielded a total of 425,505 high-quality sequence reads. A custom data analysis pipeline delineated expression profiles for 191 previously annotated miRNAs, 13 novel miRNAs, and 56 candidate miRNAs. Further characterization of a subset of the novel miRNAs in Dicer-knockdown hESC demonstrated Dicer-dependent expression, providing additional validation of our results. A set of 14 miRNAs (9 known and 5 novel) was noted to be expressed in undifferentiated hESC and then strongly downregulated with differentiation. Functional annotation analysis of predicted targets of these miRNAs and comparison with a null model using non-hESC-expressed miRNAs identified statistically enriched functional categories, including chromatin remodeling and lineage-specific differentiation annotations. Finally, integration of our data with genome-wide chromatin immunoprecipitation data on OCT4, SOX2, and NANOG binding sites implicates these transcription factors in the regulation of nine of the novel/candidate miRNAs identified here. Comparison of our results with those of recent deep sequencing studies in mouse and human ESC shows that most of the novel/candidate miRNAs found here were not identified in the other studies. The data indicate that hESC express a larger complement of miRNAs than previously appreciated, and they provide a resource for additional studies of miRNA regulation of hESC physiology. Disclosure of potential conflicts of interest is found at the end of this article.