SUMOylation of Bonus, the Drosophila homolog of Transcription Intermediary Factor 1, safeguards germline identity by recruiting repressive chromatin complexes to silence tissue-specific 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 Bonus's repressive activity. SUMOylation influences Bonus's 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.

Recognition of double-stranded DNA by the Rhodobacter sphaeroides Argonaute protein

In contrast to eukaryotic Argonaute proteins that act on RNA targets, prokaryotic Argonautes (pAgos) can target DNA, using either small RNA or small DNA guides for its recognition. Since pAgos can recognize only a single strand of DNA and lack a helicase activity, it remains unknown how double-stranded DNA can be bound both in vitro and in vivo. Here, using in vitro reconstitution and footprinting assays we analyze formation of specific complexes with target DNA by a catalytically inactive pAgo, RsAgo from Rhodobacter sphaeroides programmed with small guide RNAs. We showed that RsAgo can recognize a specific site in double-stranded DNA after stepwise reconstitution of the complex from individual oligonucleotides or after prior melting of the DNA target. When bound, RsAgo stabilizes an open DNA bubble corresponding to the length of the guide molecule and protects the target DNA from nuclease cleavage. The results suggest that RsAgo and, possibly, other RNA-guided pAgos cannot directly attack double-stranded DNA and likely require DNA opening by other cellular processes for their action.

Stable Polycomb-dependent transgenerational inheritance of chromatin states in Drosophila

Transgenerational Epigenetic Inheritance (TEI) studies the transmission of alternative functional states through multiple generations in the presence of the same genomic DNA sequence. Very little is known on the principles and the molecular mechanisms governing this type of inheritance. Here, by transiently enhancing 3D chromatin interactions, we established stable and isogenic Drosophila epilines that carry alternative epialleles, defined by differential levels of the Polycomb-dependent H3K27me3 mark. Once established, epialleles can be dominantly transmitted to naïve flies and induce paramutation. Importantly, epilines can be reset to a naïve state by disrupting chromatin interactions. Finally, we show that environmental changes can modulate the expressivity of the epialleles and we extend our paradigm to naturally occurring phenotypes. Our work sheds light on how nuclear organization and Polycomb group proteins contribute to epigenetically inheritable phenotypic variability.

The histone chaperone CAF-1 safeguards somatic cell identity

Cellular differentiation involves profound remodeling of chromatic landscapes, yet the mechanisms by which somatic cell identity is subsequently maintained remain incompletely understood. To further elucidate regulatory pathways that safeguard the somatic state, we performed two comprehensive RNAi screens targeting chromatin factors during transcription factor-mediated reprogramming of mouse fibroblasts to induced pluripotent stem cells (iPSCs). Remarkably, subunits of the chromatin assembly factor-1 (CAF-1) complex emerged as the most prominent hits from both screens, followed by modulators of lysine sumoylation and heterochromatin maintenance. Optimal modulation of both CAF-1 and transcription factor levels increased reprogramming efficiency by several orders of magnitude and facilitated iPSC formation in as little as 4 days. Mechanistically, CAF-1 suppression led to a more accessible chromatin structure at enhancer elements early during reprogramming. These changes were accompanied by a decrease in somatic heterochromatin domains, increased binding of Sox2 to pluripotency-specific targets and activation of associated genes. Notably, suppression of CAF-1 also enhanced the direct conversion of B cells into macrophages and fibroblasts into neurons. Together, our findings reveal the histone chaperone CAF-1 as a novel regulator of somatic cell identity during transcription factor-induced cell fate transitions and provide a potential strategy to modulate cellular plasticity in a regenerative setting.

Preparation of small RNA libraries for high-throughput sequencing

This protocol details the process of small RNA cloning for sequencing on the Illumina/Solexa sequencing platform, but it can be easily modified for use on other next-generation platforms (e.g., SOLiD, 454). This procedure is designed to clone canonical small RNA molecules with 5'-monophosphate and 3'-hydroxyl termini. Modifications, such as the presence of a 2'-O-methyl group, can reduce efficiency, although not sufficiently to negate the utility of the approach. Other termini modifications, such as a 5' triphosphate or a 3' phosphate, can be altered by enzymatic treatment before cloning.

A mammalian microRNA expression atlas based on small RNA library sequencing

MicroRNAs (miRNAs) are small noncoding regulatory RNAs that reduce stability and/or translation of fully or partially sequence-complementary target mRNAs. In order to identify miRNAs and to assess their expression patterns, we sequenced over 250 small RNA libraries from 26 different organ systems and cell types of human and rodents that were enriched in neuronal as well as normal and malignant hematopoietic cells and tissues. We present expression profiles derived from clone count data and provide computational tools for their analysis. Unexpectedly, a relatively small set of miRNAs, many of which are ubiquitously expressed, account for most of the differences in miRNA profiles between cell lineages and tissues. This broad survey also provides detailed and accurate information about mature sequences, precursors, genome locations, maturation processes, inferred transcriptional units, and conservation patterns. We also propose a subclassification scheme for miRNAs for assisting future experimental and computational functional analyses.

A novel class of small RNAs bind to MILI protein in mouse testes

Small RNAs bound to Argonaute proteins recognize partially or fully complementary nucleic acid targets in diverse gene-silencing processes. A subgroup of the Argonaute proteins--known as the 'Piwi family'--is required for germ- and stem-cell development in invertebrates, and two Piwi members--MILI and MIWI--are essential for spermatogenesis in mouse. Here we describe a new class of small RNAs that bind to MILI in mouse male germ cells, where they accumulate at the onset of meiosis. The sequences of the over 1,000 identified unique molecules share a strong preference for a 5' uridine, but otherwise cannot be readily classified into sequence families. Genomic mapping of these small RNAs reveals a limited number of clusters, suggesting that these RNAs are processed from long primary transcripts. The small RNAs are 26-31 nucleotides (nt) in length--clearly distinct from the 21-23 nt of microRNAs (miRNAs) or short interfering RNAs (siRNAs)--and we refer to them as 'Piwi-interacting RNAs' or piRNAs. Orthologous human chromosomal regions also give rise to small RNAs with the characteristics of piRNAs, but the cloned sequences are distinct. The identification of this new class of small RNAs provides an important starting point to determine the molecular function of Piwi proteins in mammalian spermatogenesis.

Identification of clustered microRNAs using an ab initio prediction method

Background

MicroRNAs (miRNAs) are endogenous 21 to 23-nucleotide RNA molecules that regulate protein-coding gene expression in plants and animals via the RNA interference pathway. Hundreds of them have been identified in the last five years and very recent works indicate that their total number is still larger. Therefore miRNAs gene discovery remains an important aspect of understanding this new and still widely unknown regulation mechanism. Bioinformatics approaches have proved to be very useful toward this goal by guiding the experimental investigations.

Results

In this work we describe our computational method for miRNA prediction and the results of its application to the discovery of novel mammalian miRNAs. We focus on genomic regions around already known miRNAs, in order to exploit the property that miRNAs are occasionally found in clusters. Starting with the known human, mouse and rat miRNAs we analyze 20 kb of flanking genomic regions for the presence of putative precursor miRNAs (pre-miRNAs). Each genome is analyzed separately, allowing us to study the species-specific identity and genome organization of miRNA loci. We only use cross-species comparisons to make conservative estimates of the number of novel miRNAs. Our ab initio method predicts between fifty and hundred novel pre-miRNAs for each of the considered species. Around 30% of these already have experimental support in a large set of cloned mammalian small RNAs. The validation rate among predicted cases that are conserved in at least one other species is higher, about 60%, and many of them have not been detected by prediction methods that used cross-species comparisons. A large fraction of the experimentally confirmed predictions correspond to an imprinted locus residing on chromosome 14 in human, 12 in mouse and 6 in rat. Our computational tool can be accessed on the world-wide-web.

Conclusion

Our results show that the assumption that many miRNAs occur in clusters is fruitful for the discovery of novel miRNAs. Additionally we show that although the overall miRNA content in the observed clusters is very similar across the three considered species, the internal organization of the clusters changes in evolution.

Identification and characterization of small RNAs involved in RNA silencing

Double-stranded RNA (dsRNA) is a potent trigger of sequence-specific gene silencing mechanisms known as RNA silencing or RNA interference. The recognition of the target sequences is mediated by ribonucleoprotein complexes that contain 21- to 28-nucleotide (nt) guide RNAs derived from processing of the trigger dsRNA. Here, we review the experimental and bioinformatic approaches that were used to identify and characterize these small RNAs isolated from cells and tissues. The identification and characterization of small RNAs and their expression patterns is important for elucidating gene regulatory networks.

Clustering and conservation patterns of human microRNAs

MicroRNAs (miRNAs) are ∼22 nt-long non-coding RNA molecules, believed to play important roles in gene regulation. We present a comprehensive analysis of the conservation and clustering patterns of known miRNAs in human. We show that human miRNA gene clustering is significantly higher than expected at random. A total of 37% of the known human miRNA genes analyzed in this study appear in clusters of two or more with pairwise chromosomal distances of at most 3000 nt. Comparison of the miRNA sequences with their homologs in four other organisms reveals a typical conservation pattern, persistent throughout the clusters. Furthermore, we show enrichment in the typical conservation patterns and other miRNA-like properties in the vicinity of known miRNA genes, compared with random genomic regions. This may imply that additional, yet unknown, miRNAs reside in these regions, consistent with the current recognition that there are overlooked miRNAs. Indeed, by comparing our predictions with cloning results and with identified miRNA genes in other mammals, we corroborate the predictions of 18 additional human miRNA genes in the vicinity of the previously known ones. Our study raises the proportion of clustered human miRNAs that are <3000 nt apart to 42%. This suggests that the clustering of miRNA genes is higher than currently acknowledged, alluding to its evolutionary and functional implications.

Human MicroRNA targets

MicroRNAs (miRNAs) interact with target mRNAs at specific sites to induce cleavage of the message or inhibit translation. The specific function of most mammalian miRNAs is unknown. We have predicted target sites on the 3' untranslated regions of human gene transcripts for all currently known 218 mammalian miRNAs to facilitate focused experiments. We report about 2,000 human genes with miRNA target sites conserved in mammals and about 250 human genes conserved as targets between mammals and fish. The prediction algorithm optimizes sequence complementarity using position-specific rules and relies on strict requirements of interspecies conservation. Experimental support for the validity of the method comes from known targets and from strong enrichment of predicted targets in mRNAs associated with the fragile X mental retardation protein in mammals. This is consistent with the hypothesis that miRNAs act as sequence-specific adaptors in the interaction of ribonuclear particles with translationally regulated messages. Overrepresented groups of targets include mRNAs coding for transcription factors, components of the miRNA machinery, and other proteins involved in translational regulation, as well as components of the ubiquitin machinery, representing novel feedback loops in gene regulation. Detailed information about target genes, target processes, and open-source software for target prediction (miRanda) is available at http://www.microrna.org. Our analysis suggests that miRNA genes, which are about 1% of all human genes, regulate protein production for 10% or more of all human genes.