Identification of novel Drosophila melanogaster microRNAs

Epigenetic Molecular Mechanisms in Insects

Insects are the largest animal group on Earth both in biomass and diversity. Their outstanding success has inspired genetics and developmental research, allowing the discovery of dynamic process explaining extreme phenotypic plasticity and canalization. Epigenetic molecular mechanisms (EMMs) are vital for several housekeeping functions in multicellular organisms, regulating developmental, ontogenetic trajectories and environmental adaptations. In Insecta, EMMs are involved in the development of extreme phenotypic divergences such as polyphenisms and eusocial castes. Here, we review the history of this research field and how the main EMMs found in insects help to understand their biological processes and diversity. EMMs in insects confer them rapid response capacity allowing insect either to change with plastic divergence or to keep constant when facing different stressors or stimuli. EMMs function both at intra as well as transgenerational scales, playing important roles in insect ecology and evolution. We discuss on how EMMs pervasive influences in Insecta require not only the control of gene expression but also the dynamic interplay of EMMs with further regulatory levels, including genetic, physiological, behavioral, and environmental among others, as was earlier proposed by the Probabilistic Epigenesis model and Developmental System Theory.

Identification and characterization of novel conserved RNA structures in Drosophila

BACKGROUND

Comparative genomics approaches have facilitated the discovery of many novel non-coding and structured RNAs (ncRNAs). The increasing availability of related genomes now makes it possible to systematically search for compensatory base changes - and thus for conserved secondary structures - even in genomic regions that are poorly alignable in the primary sequence. The wealth of available transcriptome data can add valuable insight into expression and possible function for new ncRNA candidates. Earlier work identifying ncRNAs in Drosophila melanogaster made use of sequence-based alignments and employed a sliding window approach, inevitably biasing identification toward RNAs encoded in the more conserved parts of the genome.

RESULTS

To search for conserved RNA structures (CRSs) that may not be highly conserved in sequence and to assess the expression of CRSs, we conducted a genome-wide structural alignment screen of 27 insect genomes including D. melanogaster and integrated this with an extensive set of tiling array data. The structural alignment screen revealed ∼30,000 novel candidate CRSs at an estimated false discovery rate of less than 10%. With more than one quarter of all individual CRS motifs showing sequence identities below 60%, the predicted CRSs largely complement the findings of sliding window approaches applied previously. While a sixth of the CRSs were ubiquitously expressed, we found that most were expressed in specific developmental stages or cell lines. Notably, most statistically significant enrichment of CRSs were observed in pupae, mainly in exons of untranslated regions, promotors, enhancers, and long ncRNAs. Interestingly, cell lines were found to express a different set of CRSs than were found in vivo. Only a small fraction of intergenic CRSs were co-expressed with the adjacent protein coding genes, which suggests that most intergenic CRSs are independent genetic units.

CONCLUSIONS

This study provides a more comprehensive view of the ncRNA transcriptome in fly as well as evidence for differential expression of CRSs during development and in cell lines.

Identification and characterization of microRNAs in the white-backed planthopper, Sogatella furcifera

MicroRNAs (miRNAs) are a novel class of small, non-coding endogenous RNAs that play critical regulatory roles in many metabolic activities in eukaryotes. Reports of the identification of miRNAs in Sogatella furcifera (white-backed planthopper), the insect that acts as the only confirmed vector of the southern rice black-streaked dwarf virus (SRBSDV), are limited. In this study, a total of 382 miRNAs were identified in S. furcifera, including 106 conserved and 276 novel miRNAs, using high-throughput sequencing based on two small RNA libraries from viruliferous and non-viruliferous S. furcifera, and these miRNAs belonged to 52 conserved miRNA families and 58 S. furcifera-specific families, respectively. Comparison with miRNAs from 26 insect species and five other species in miRBase showed that more than half of the conserved miRNA families are highly conserved in Hexapoda, while other miRNAs are only conserved in non-dipterans. Furthermore, 4 117 target genes predicted for the 382 identified miRNAs could be categorized into 45 functional groups annotated by Gene Ontology. Compared with non-viruliferous cells, eight up-regulated miRNAs and four down-regulated miRNAs were identified in cells inoculated with SRBSDV, among which miR-14 and miR-n98a may be involved in the immune response to SRBSDV infection. Analyses of the identified miRNAs will provide insights into the roles of these miRNAs in the regulation and expression of genes involved in the metabolism, development and viral infection of S. furcifera.

Identification and developmental profiling of conserved and novel microRNAs in Manduca sexta

MicroRNAs (miRNAs) are a group of small RNAs involved in translation inhibition or mRNA degradation. Due to its large size, Manduca sexta has long been used as a model to study insect physiology and biochemistry. While transcriptome studies have greatly enriched our knowledge on M. sexta structural genes, little is known about posttranscriptional regulation by miRNAs in this lepidopteran species. We constructed four small RNA libraries from embryos, 4th instar feeding larvae, pupae, and adults, obtained 21 million reads of 18-31 nucleotides by Illumina sequencing, and found 163 conserved and 13 novel miRNAs. By searching the M. sexta genome assembly, we identified precursors of 82 conserved miRNAs, 76 of which had mapped reads in one or more of these libraries. After normalization, we compared numbers of miRNA and miRNA-star reads in these libraries and observed abundance changes during development. Interestingly, mse-miR-281-star, mse-miR-31-star, mse-miR-965-star, mse-miR-9a-star, mse-miR-9b-star, mse-miR-2a-star, mse-miR-92b-star and mse-miR-279c-star are either more abundant or maintained at similar levels compared to respective mature miRNA strand. Expression profiling of the first set of miRNAs provided insights to their possible involvement in developmental regulation. This study will aid in the annotation of miRNA genes in the genome.

MicroRNAs in Drosophila development

Micro-ribonucleic acids (miRNAs) are small (21-24 nucleotide), endogenously expressed, noncoding RNAs that have emerged as important posttranscriptional regulators of gene expression. MiRNAs have been identified and cloned from diverse eukaryotic organisms where they have been shown to control important physiological and developmental processes such as apoptosis, cell division, and differentiation. A high level of conservation of some miRNAs across phyla further emphasizes their importance as posttranscriptional regulators. Research in a variety of model systems has been instrumental in dissecting the biological functions of miRNAs. In this chapter, we discuss the current literature on the role of miRNAs as developmental regulators in Drosophila.

Computational and experimental identification of mirtrons in Drosophila melanogaster and Caenorhabditis elegans

Mirtrons are intronic hairpin substrates of the dicing machinery that generate functional microRNAs. In this study, we describe experimental assays that defined the essential requirements for entry of introns into the mirtron pathway. These data informed a bioinformatic screen that effectively identified functional mirtrons from the Drosophila melanogaster transcriptome. These included 17 known and six confident novel mirtrons among the top 51 candidates, and additional candidates had limited read evidence in available small RNA data. Our computational model also proved effective on Caenorhabditis elegans, for which the identification of 14 cloned mirtrons among the top 22 candidates more than tripled the number of validated mirtrons in this species. A few low-scoring introns generated mirtron-like read patterns from atypical RNA structures, but their paucity suggests that relatively few such loci were not captured by our model. Unexpectedly, we uncovered examples of clustered mirtrons in both fly and worm genomes, including a <8-kb region in C. elegans harboring eight distinct mirtrons. Altogether, we demonstrate that discovery of functional mirtrons, unlike canonical miRNAs, is amenable to computational methods independent of evolutionary constraint.

Deep annotation of Drosophila melanogaster microRNAs yields insights into their processing, modification, and emergence

Since the initial annotation of miRNAs from cloned short RNAs by the Ambros, Tuschl, and Bartel groups in 2001, more than a hundred studies have sought to identify additional miRNAs in various species. We report here a meta-analysis of short RNA data from Drosophila melanogaster, aggregating published libraries with 76 data sets that we generated for the modENCODE project. In total, we began with more than 1 billion raw reads from 187 libraries comprising diverse developmental stages, specific tissue- and cell-types, mutant conditions, and/or Argonaute immunoprecipitations. We elucidated several features of known miRNA loci, including multiple phased byproducts of cropping and dicing, abundant alternative 5' termini of certain miRNAs, frequent 3' untemplated additions, and potential editing events. We also identified 49 novel genomic locations of miRNA production, and 61 additional candidate loci with limited evidence for miRNA biogenesis. Although these loci broaden the Drosophila miRNA catalog, this work supports the notion that a restricted set of cellular transcripts is competent to be specifically processed by the Drosha/Dicer-1 pathway. Unexpectedly, we detected miRNA production from coding and untranslated regions of mRNAs and found the phenomenon of miRNA production from the antisense strand of known loci to be common. Altogether, this study lays a comprehensive foundation for the study of miRNA diversity and evolution in a complex animal model.

MicroRNA biogenesis via splicing and exosome-mediated trimming in Drosophila

microRNAs (miRNAs) are approximately 22 nucleotide regulatory RNAs derived from hairpins generated either by Drosha cleavage (canonical substrates) or by splicing and debranching of short introns (mirtrons). The 5' end of the highly conserved Drosophila mirtron-like locus mir-1017 is coincident with the splice donor, but a substantial "tail" separates its hairpin from the 3'splice acceptor. Genetic and biochemical studies define a biogenesis pathway involving splicing, lariat debranching, and RNA exosome-mediated "trimming," followed by conventional dicing and loading into AGO1 to yield a miRNA that can repress seed-matched targets. Analysis of cloned small RNAs yielded six additional candidate 3' tailed mirtrons in D. melanogaster. Altogether, these data reveal an unexpected role for the exosome in the biogenesis of miRNAs from hybrid mirtron substrates.

Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development

BACKGROUND

In eukaryotes, microRNAs (miRNAs) have emerged as critical regulators of gene expression. The Silkworm (Bombyx mori L.) is one of the most suitable lepidopteran insects for studying the molecular aspects of metamorphosis because of its large size, availability of mutants and genome sequence. Besides, this insect also has been amply studied from a physiological and biochemical perspective. Deep sequencing of small RNAs isolated from different stages of silkworm is a powerful tool not only for measuring the changes in miRNA profile but also for discovering novel miRNAs.

RESULTS

We generated small RNA libraries from feeding larvae, spinning larvae, pupae and adults of B. mori and obtained approximately 2.5 million reads of 18-30 nt. Sequence analysis identified 14 novel and 101 conserved miRNAs. Most novel miRNAs are preferentially expressed in pupae, whereas more than 95% of the conserved miRNAs are dynamically regulated during different developmental stages. Remarkably, the miRNA-star (miR*) of four miRNAs are expressed at much higher levels than their corresponding miRNAs, and their expression profiles are distinct from their corresponding miRNA profiles during different developmental stages. Additionally, we detected two antisense miRNA loci (miR-263-S and miR-263-AS; miR-306-S and miR-306-AS) that are expressed in sense and antisense directions. Interestingly, miR-263 and miR-306 are preferentially and abundantly expressed in pupae and adults, respectively.

CONCLUSIONS

We identified 101 homologs of conserved miRNAs, 14 species-specific and two antisense miRNAs in the silkworm. Our results provided deeper insights into changes in conserved and novel miRNA and miRNA* accumulation during development.

Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development

BACKGROUND

In eukaryotes, microRNAs (miRNAs) have emerged as critical regulators of gene expression. The Silkworm (Bombyx mori L.) is one of the most suitable lepidopteran insects for studying the molecular aspects of metamorphosis because of its large size, availability of mutants and genome sequence. Besides, this insect also has been amply studied from a physiological and biochemical perspective. Deep sequencing of small RNAs isolated from different stages of silkworm is a powerful tool not only for measuring the changes in miRNA profile but also for discovering novel miRNAs.

RESULTS

We generated small RNA libraries from feeding larvae, spinning larvae, pupae and adults of B. mori and obtained approximately 2.5 million reads of 18-30 nt. Sequence analysis identified 14 novel and 101 conserved miRNAs. Most novel miRNAs are preferentially expressed in pupae, whereas more than 95% of the conserved miRNAs are dynamically regulated during different developmental stages. Remarkably, the miRNA-star (miR*) of four miRNAs are expressed at much higher levels than their corresponding miRNAs, and their expression profiles are distinct from their corresponding miRNA profiles during different developmental stages. Additionally, we detected two antisense miRNA loci (miR-263-S and miR-263-AS; miR-306-S and miR-306-AS) that are expressed in sense and antisense directions. Interestingly, miR-263 and miR-306 are preferentially and abundantly expressed in pupae and adults, respectively.

CONCLUSIONS

We identified 101 homologs of conserved miRNAs, 14 species-specific and two antisense miRNAs in the silkworm. Our results provided deeper insights into changes in conserved and novel miRNA and miRNA* accumulation during development.

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.

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.

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.

MiR-433 and miR-127 arise from independent overlapping primary transcripts encoded by the miR-433-127 locus

MicroRNAs play significant roles in development, metabolism and carcinogenesis, however, limited information is available about their primary transcripts and the transcriptional regulation of the microRNA genes. We report here the cloning of two primary miRNAs (pri-miR-433 and pri-miR-127) encoded by the miR-433-127 locus. Using both database mining and experimental methods, we isolated the full-length primary transcripts of the mouse miR-433 and miR-127 and demonstrated that they overlapped in a 5'-3' unidirectional way. These two miRNA genes are expressed in a compact space by using overlapping coding regions. This is the first report to identify an economical structure for miR-433 and miR-127 genes, which may be a novel way of miRNA gene to maximize the genetic information in order to fit the complex physiological function of mammalian organism.

Transcriptional mechanism for the paired miR-433 and miR-127 genes by nuclear receptors SHP and ERRgamma

MicroRNAs (miRNAs, miRs) are genomically encoded small ∼22 nt RNA molecules that have been shown to mediate translational repression of target mRNAs involved in cellular proliferation, differentiation and death. Despite intensive studies on their physiological and pathological functions, the molecular mechanism of how miRNA gene transcription is regulated remains largely unknown. Microarray profiling revealed 21 miRNAs clustered on chromosome 12, including miR-433 and miR-127, that were co-upregulated in small heterodimer partner (SHP, NR0B2) SHP knockouts (SHP^–/–) liver. Gene cloning revealed that the 3′-coding region of pri-miR-433 served as the promoter region of pri-miR-127. Estrogen related receptor (ERRγ, NR3B3) robustly activated miR-433 and miR-127 promoter reporters through ERRE, which was transrepressed by SHP. The strong elevation of miR-433 and miR-127 in Hepa-1 cells correlated with the down-regulation of SHP and up-regulation of ERRγ. Ectopic expression of ERRγ induced miR-433 and miR-127 expression, which was repressed by SHP coexpression. In contrast, knockdown ERRγ decreased miR-433 and miR-127 expression. In addition, the ERRγ agonist GSK4716 induced miR-433 and miR-127 expression both in vitro and in vivo, respectively. In summary, the coupled miR-433 and miR-127 genes were transcribed from independent promoters regulated by nuclear receptors ERRγ/SHP in a compact space by using overlapping genomic regions.

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.

Dinucleotide controlled null models for comparative RNA gene prediction

Background

Comparative prediction of RNA structures can be used to identify functional noncoding RNAs in genomic screens. It was shown recently by Babak et al. [BMC Bioinformatics. 8:33] that RNA gene prediction programs can be biased by the genomic dinucleotide content, in particular those programs using a thermodynamic folding model including stacking energies. As a consequence, there is need for dinucleotide-preserving control strategies to assess the significance of such predictions. While there have been randomization algorithms for single sequences for many years, the problem has remained challenging for multiple alignments and there is currently no algorithm available.

Results

We present a program called SISSIz that simulates multiple alignments of a given average dinucleotide content. Meeting additional requirements of an accurate null model, the randomized alignments are on average of the same sequence diversity and preserve local conservation and gap patterns. We make use of a phylogenetic substitution model that includes overlapping dependencies and site-specific rates. Using fast heuristics and a distance based approach, a tree is estimated under this model which is used to guide the simulations. The new algorithm is tested on vertebrate genomic alignments and the effect on RNA structure predictions is studied. In addition, we directly combined the new null model with the RNAalifold consensus folding algorithm giving a new variant of a thermodynamic structure based RNA gene finding program that is not biased by the dinucleotide content.

Conclusion

SISSIz implements an efficient algorithm to randomize multiple alignments preserving dinucleotide content. It can be used to get more accurate estimates of false positive rates of existing programs, to produce negative controls for the training of machine learning based programs, or as standalone RNA gene finding program. Other applications in comparative genomics that require randomization of multiple alignments can be considered.

Availability

SISSIz is available as open source C code that can be compiled for every major platform and downloaded here: .

The effect of central loops in miRNA:MRE duplexes on the efficiency of miRNA-mediated gene regulation

MicroRNAs (miRNAs) guide posttranscriptional repression of mRNAs. Hundreds of miRNAs have been identified but the target identification of mammalian mRNAs is still a difficult task due to a poor understanding of the interaction between miRNAs and the miRNA recognizing element (MRE). In recent research, the importance of the 5' end of the miRNA:MRE duplex has been emphasized and the effect of the tail region addressed, but the role of the central loop has largely remained unexplored. Here we examined the effect of the loop region in miRNA:MRE duplexes and found that the location of the central loop is one of the important factors affecting the efficiency of gene regulation mediated by miRNAs. It was further determined that the addition of a loop score combining both location and size as a new criterion for predicting MREs and their cognate miRNAs significantly decreased the false positive rates and increased the specificity of MRE prediction.