The spatiotemporal control of gene expression is crucial for the successful completion of animal development. The evolutionary constraints on development are particularly strong for the mid-embryonic stage when body segments are specified, as evidenced by a high degree of morphological and protein-coding gene conservation during this period--a phenomenon known as the developmental hourglass. The discovery of microRNA-mediated gene control revealed an entirely new layer of complexity of the molecular networks that orchestrate development. However, the constraints on microRNA developmental expression and evolution, and the implications for animal evolution are less well understood. To systematically explore the conservation of microRNAs during development, we carried out a genome-wide comparative study of microRNA expression levels throughout the ontogenesis of two divergent fruit flies, Drosophila melanogaster and D. virilis. We show that orthologous microRNAs display highly similar temporal profiles regardless of their mutation rates, suggesting that the timely expression of microRNA genes can be more constrained than their sequence. Furthermore, transitions between key developmental events in the different species are accompanied by conserved shifts in microRNA expression profiles, with the mid-embryonic period between gastrulation and segmentation characterized by the highest similarity of microRNA expression. The conservation of microRNA expression therefore displays an hourglass pattern similar to that observed for protein-coding genes.

Background

Piwi-interacting RNAs (piRNAs) are a class of short (~26-31-nucleotide) non-protein-coding RNAs expressed in the metazoan germline. The piRNA pathway in arthropods is best understood in the ovary of Drosophila melanogaster, where it acts to silence active transposable elements (TEs). Maternal loading of piRNAs in oocytes is further required for the inheritance of piRNA-mediated transposon defence. However, our understanding of the diversity, evolution and function of the piRNA complement beyond drosophilids is limited. The red flour beetle, Tribolium castaneum, is an emerging model organism separated from Drosophila by ~ 350 million years of evolution that displays a number of features ancestral to arthropods, including short germ embryogenesis. Here, we characterize the maternally deposited and zygotically expressed small RNA and mRNA complements throughout T. castaneum embryogenesis.

Results

We find that beetle oocytes and embryos of all stages are abundant in heterogeneous ~ 28-nucleotide RNAs. These small RNAs originate from discrete genomic loci enriched in TE sequences and display the molecular signatures of transposon-derived piRNAs. In addition to the maternally loaded primary piRNAs, Tribolium embryos produce secondary piRNAs by the cleavage of zygotically activated TE transcripts via the ping-pong mechanism. The two Tribolium piRNA pathway effector proteins, Tc-Piwi/Aub and Tc-Ago3, are also expressed throughout the soma of early embryos.

Conclusions

Our results show that the piRNA pathway in Tribolium is not restricted to the germline, but also operates in the embryo and may act to antagonize zygotically activated transposons. Taken together, these data highlight a functional divergence of the piRNA pathway between insects.

MicroRNAs are well-established players in the development of multicellular animals. Most of our understanding of microRNA function in arthropod development comes from studies in Drosophila. Despite their advantages as model systems, the long germband embryogenesis of fruit flies is an evolutionary derived state restricted to several holometabolous insect lineages. MicroRNA evolution and expression across development in animals exhibiting the ancestral and more widespread short germband mode of embryogenesis has not been characterized. We sequenced small RNA libraries of oocytes and successive intervals covering the embryonic development of the short germband model organism, Tribolium castaneum. We analyzed the evolution and temporal expression of the microRNA complement and sequenced libraries of total RNA to investigate the relationships with microRNA target expression. We show microRNA maternal loading and sequence-specific 3' end nontemplate oligoadenylation of maternally deposited microRNAs that is conserved between Tribolium and Drosophila. We further uncover large clusters encoding multiple paralogs from several Tribolium-specific microRNA families expressed during a narrow interval of time immediately after the activation of zygotic transcription. These novel microRNAs, together with several early expressed conserved microRNAs, target a significant number of maternally deposited transcripts. Comparison with Drosophila shows that microRNA-mediated maternal transcript targeting is a conserved process in insects, but the number and sequences of microRNAs involved have diverged. The expression of fast-evolving and species-specific microRNAs in the early blastoderm of T. castaneum is consistent with previous findings in Drosophila and shows that the unique permissiveness for microRNA innovation at this stage is a conserved phenomenon.

MicroRNAs are short non-protein-coding RNAs that regulate gene expression at the post-transcriptional level and are essential for the embryonic development of multicellular animals. Comparative genome-scale analyses have revealed that metazoan evolution is accompanied by the continuous acquisition of novel microRNA genes. This suggests that novel microRNAs may promote innovation and diversity in development. We determined the evolutionary origins of extant Drosophila microRNAs and estimated the sequence divergence between the 130 orthologous microRNAs in Drosophila melanogaster and Drosophila virilis, separated by 63 million years of evolution. We then generated small RNA sequencing data sets covering D. virilis development and explored the relationship between microRNA conservation and expression in a developmental context. We find that late embryonic, larval, and adult stages are dominated by conserved microRNAs. This pattern, however, does not hold for the early embryo, where rapidly evolving microRNAs are uniquely present at high levels in both species. The group of fast-evolving microRNAs that are highly expressed in the early embryo belong to two Drosophilid lineage-specific clusters: mir-310 ∼ 313 and mir-309 ∼ 6. These clusters have particularly complex evolutionary histories of duplication, gain, and loss. Our analyses suggest that the early embryo is a more permissive environment for microRNA changes and innovations. Fast-evolving microRNAs, therefore, have the opportunity to become preferentially integrated in early developmental processes, and may impact the evolution of development. The relationship between microRNA conservation and expression throughout the development of Drosophila differs from that previously observed for protein-coding genes.

The conserved family of Transcription Intermediary Factors (TIF1) proteins consists of key transcriptional regulators that control transcription of target genes by modulating chromatin state. Unlike mammals that have four TIF1 members, Drosophila only encodes one member of the family, Bonus. Bonus has been implicated in embryonic development and organogenesis and shown to regulate several signaling pathways, however, its targets and mechanism of action remained poorly understood. We found that knockdown of Bonus in early oogenesis results in severe defects in ovarian development and in ectopic expression of genes that are normally repressed in the germline, demonstrating its essential function in the ovary. Recruitment of Bonus to chromatin leads to silencing associated with accumulation of the repressive H3K9me3 mark. We show that Bonus associates with the histone methyltransferase SetDB1 and the chromatin remodeler NuRD and depletion of either component releases Bonus-induced repression. We further established that Bonus is SUMOylated at a single site at its N-terminus that is conserved among insects and this modification is indispensable for Bonuss repressive activity. SUMOylation influences Bonuss subnuclear localization, its association with chromatin and interaction with SetDB1. Finally, we showed that Bonus SUMOylation is mediated by the SUMO E3-ligase Su(var)2-10, revealing that although SUMOylation of TIF1 proteins is conserved between insects and mammals, both the mechanism and specific site of modification is different in the two taxa. Together, our work identified Bonus as a regulator of tissue-specific gene expression and revealed the importance of SUMOylation as a regulator of complex formation in the context of transcriptional repression.

MicroRNAs are small RNAs that regulate protein levels. It is commonly assumed that the expression level of a microRNA is directly correlated with its repressive activity - that is, highly expressed microRNAs will repress their target mRNAs more. Here we investigate the quantitative relationship between endogenous microRNA expression and repression for 32 mature microRNAs in Drosophila melanogaster S2 cells. In general, we find that more abundant microRNAs repress their targets to a greater degree. However, the relationship between expression and repression is nonlinear, such that a 10-fold greater microRNA concentration produces only a 10% increase in target repression. The expression/repression relationship is the same for both dominant guide microRNAs and minor mature products (so-called passenger strands/microRNA* sequences). However, we find examples of microRNAs whose cellular concentrations differ by several orders of magnitude, yet induce similar repression of target mRNAs. Likewise, microRNAs with similar expression can have very different repressive abilities. We show that the association of microRNAs with Argonaute proteins does not explain this variation in repression. The observed relationship is consistent with the limiting step in target repression being the association of the microRNA/RISC complex with the target site. These findings argue that modest changes in cellular microRNA concentration will have minor effects on repression of targets.

Members of the conserved Argonaute (Ago) protein family provide defence against invading nucleic acids in eukaryotes in the process of RNA interference. Many prokaryotes also contain Ago proteins that are predicted to be active nucleases; however, their functional activities in host cells remain poorly understood. Here, we characterize the in vitro and in vivo properties of the SeAgo protein from the mesophilic cyanobacterium Synechococcus elongatus. We show that SeAgo is a DNA-guided nuclease preferentially acting on single-stranded DNA targets, with non-specific guide-independent activity observed for double-stranded substrates. The SeAgo gene is steadily expressed in S. elongatus; however, its deletion or overexpression does not change the kinetics of cell growth. When purified from its host cells or from heterologous E. coli, SeAgo is loaded with small guide DNAs whose formation depends on the endonuclease activity of the argonaute protein. SeAgo co-purifies with SSB proteins suggesting that they may also be involved in DNA processing. The SeAgo-associated small DNAs are derived from diverse genomic locations, with certain enrichment for the proposed sites of chromosomal replication initiation and termination, but show no preference for an endogenous plasmid. Therefore, promiscuous genome sampling by SeAgo does not have great effects on cell physiology and plasmid maintenance.

Ribosomal RNAs (rRNAs) are essential components of the ribosome and are among the most abundant macromolecules in the cell. To ensure high rRNA level, eukaryotic genomes contain dozens to hundreds of rDNA genes, however, only a fraction of the rRNA genes seems to be active, while others are transcriptionally silent. We found that individual rDNA genes have high level of cell-to-cell heterogeneity in their expression in Drosophila melanogaster. Insertion of heterologous sequences into rDNA leads to repression associated with reduced expression in individual cells and decreased number of cells expressing rDNA with insertions. We found that SUMO (Small Ubiquitin-like Modifier) and SUMO ligase Ubc9 are required for efficient repression of interrupted rDNA units and variable expression of intact rDNA. Disruption of the SUMO pathway abolishes discrimination of interrupted and intact rDNAs and removes cell-to-cell heterogeneity leading to uniformly high expression of individual rDNA in single cells. Our results suggest that the SUMO pathway is responsible for both repression of interrupted units and control of intact rDNA expression.

The conserved THO/TREX (transcription/export) complex is critical for pre-mRNA processing and mRNA nuclear export. In metazoa, TREX is loaded on nascent RNA transcribed by RNA polymerase II in a splicing-dependent fashion; however, how TREX functions is poorly understood. Here we show that Thoc5 and other TREX components are essential for the biogenesis of piRNA, a distinct class of small noncoding RNAs that control expression of transposable elements (TEs) in the Drosophila germline. Mutations in TREX lead to defects in piRNA biogenesis, resulting in derepression of multiple TE families, gametogenesis defects, and sterility. TREX components are enriched on piRNA precursors transcribed from dual-strand piRNA clusters and colocalize in distinct nuclear foci that overlap with sites of piRNA transcription. The localization of TREX in nuclear foci and its loading on piRNA precursor transcripts depend on Cutoff, a protein associated with chromatin of piRNA clusters. Finally, we show that TREX is required for accumulation of nascent piRNA precursors. Our study reveals a novel splicing-independent mechanism for TREX loading on nascent RNA and its importance in piRNA biogenesis.

MicroRNAs are small (∼22 nt) noncoding RNAs that repress translation and therefore regulate the production of proteins from specific target mRNAs. microRNAs have been found to function in diverse aspects of gene regulation within animal development and many other processes. Among invertebrates, both conserved and novel, lineage specific, microRNAs have been extensively studied predominantly in holometabolous insects such as Drosophila melanogaster However little is known about microRNA repertoires in other arthropod lineages such as the chelicerates. To understand the evolution of microRNAs in this poorly sampled subphylum, we characterized the microRNA repertoire expressed during embryogenesis of the common house spider Parasteatoda tepidariorum We identified a total of 148 microRNAs in P. tepidariorum representing 66 families. Approximately half of these microRNA families are conserved in other metazoans, while the remainder are specific to this spider. Of the 35 conserved microRNAs families 15 had at least two copies in the P. tepidariorum genome. A BLAST-based approach revealed a similar pattern of duplication in other spiders and a scorpion, but not among other chelicerates and arthropods, with the exception of a horseshoe crab. Among the duplicated microRNAs we found examples of lineage-specific tandem duplications, and the duplication of entire microRNA clusters in three spiders, a scorpion, and in a horseshoe crab. Furthermore, we found that paralogs of many P. tepidariorum microRNA families exhibit arm switching, which suggests that duplication was often followed by sub- or neofunctionalization. Our work shows that understanding the evolution of microRNAs in the chelicerates has great potential to provide insights into the process of microRNA duplication and divergence and the evolution of animal development.