Methylation of cytosines is an essential epigenetic modification in mammalian genomes, yet the rules that govern methylation patterns remain largely elusive. To gain insights into this process, we generated base-pair-resolution mouse methylomes in stem cells and neuronal progenitors. Advanced quantitative analysis identified low-methylated regions (LMRs) with an average methylation of 30%. These represent CpG-poor distal regulatory regions as evidenced by location, DNase I hypersensitivity, presence of enhancer chromatin marks and enhancer activity in reporter assays. LMRs are occupied by DNA-binding factors and their binding is necessary and sufficient to create LMRs. A comparison of neuronal and stem-cell methylomes confirms this dependency, as cell-type-specific LMRs are occupied by cell-type-specific transcription factors. This study provides methylome references for the mouse and shows that DNA-binding factors locally influence DNA methylation, enabling the identification of active regulatory regions.

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.

The factors that sequester transcriptionally repressed heterochromatin at the nuclear periphery are currently unknown. In a genome-wide RNAi screen, we found that depletion of S-adenosylmethionine (SAM) synthetase reduces histone methylation globally and causes derepression and release of heterochromatin from the nuclear periphery in Caenorhabditis elegans embryos. Analysis of histone methyltransferases (HMTs) showed that elimination of two HMTs, MET-2 and SET-25, mimics the loss of SAM synthetase, abrogating the perinuclear attachment of heterochromatic transgenes and of native chromosomal arms rich in histone H3 lysine 9 methylation. The two HMTs target H3K9 in a consecutive fashion: MET-2, a SETDB1 homolog, mediates mono- and dimethylation, and SET-25, a previously uncharacterized HMT, deposits H3K9me3. SET-25 colocalizes with its own product in perinuclear foci, in a manner dependent on H3K9me3, but not on its catalytic domain. This colocalization suggests an autonomous, self-reinforcing mechanism for the establishment and propagation of repeat-rich heterochromatin.

Loss of microRNA (miRNA) pathway components negatively affects differentiation of embryonic stem (ES) cells, but the underlying molecular mechanisms remain poorly defined. Here we characterize changes in mouse ES cells lacking Dicer (Dicer1). Transcriptome analysis of Dicer-/- cells indicates that the ES-specific miR-290 cluster has an important regulatory function in undifferentiated ES cells. Consistently, many of the defects in Dicer-deficient cells can be reversed by transfection with miR-290 family miRNAs. We demonstrate that Oct4 (also known as Pou5f1) silencing in differentiating Dicer-/- ES cells is accompanied by accumulation of repressive histone marks but not by DNA methylation, which prevents the stable repression of Oct4. The methylation defect correlates with downregulation of de novo DNA methyltransferases (Dnmts). The downregulation is mediated by Rbl2 and possibly other transcriptional repressors, potential direct targets of miR-290 cluster miRNAs. The defective DNA methylation can be rescued by ectopic expression of de novo Dnmts or by transfection of the miR-290 cluster miRNAs, indicating that de novo DNA methylation in ES cells is controlled by miRNAs.

A n A ly s i s RNA-seq experiments generate reads derived not only from mature RNA transcripts but also from pre-mRNA. Here we present a computational approach called exon-intron split analysis (EISA) that measures changes in mature RNA and pre-mRNA reads across different experimental conditions to quantify transcriptional and post-transcriptional regulation of gene expression. We apply EISA to 17 diverse data sets to show that most intronic reads arise from nuclear RNA and changes in intronic read counts accurately predict changes in transcriptional activity. Furthermore, changes in posttranscriptional regulation can be predicted from differences between exonic and intronic changes. EISA reveals both transcriptional and post-transcriptional contributions to expression changes, increasing the amount of information that can be gained from RNA-seq data sets.Cellular RNAs are regulated at multiple stages, including transcription, RNA maturation and degradation. Several analytic methods have been developed to measure these processes on a transcriptomewide scale. For example, global run-on sequencing (GRO-seq) 1 uses incorporation of a nucleotide analog to enrich for nascent RNA. In Nascent-seq 2,3 , newly transcribed RNAs are isolated by purification of their complex with proteins and the DNA template. Cellular fractionation techniques 4 have also been adapted to measure nascent transcripts, which are enriched in the nucleus. mRNA half-lives have been determined, for example, by blockage of transcription followed by transcriptional profiling 5 . RNA sequencing, the most widely used method for transcriptome analysis, has been applied in numerous studies to determine steady-state mRNA levels 6,7 and alternative splicing events 8 and to identify previously unknown transcripts and noncoding RNAs 9-11 . In general, these protocols aim to enrich for mature mRNA by selection of polyadenylated RNA or by depletion of ribosomal RNA.Many computational methods (reviewed in refs. 12,13) have been developed for the analysis of RNA-seq data, to enable spliced alignment 14,15 , transcript assembly 16,17 , transcript quantification 14,18 and differential expression analysis [19][20][21] . Although RNA-seq mostly generates reads that map to exons, it also captures less abundant intronic sequences 6 . However, their interpretation has remained controversial. Some have suggested that they originate from DNA contamination and can thus be used as a quality metric for RNA-seq data 22 (see also RNA-seq guidelines of the Roadmap Epigenomics Consortium, http://www.roadmapepigenomics.org/), whereas others have hypothesized that they stem from unknown exons or intronic enhancers 6,7 . In a study based on exon arrays, probes mapping to introns were used to investigate pre-mRNA dynamics 23 . Three recent studies based on RNA-seq provided evidence that intronic reads might correlate with transcriptional activity. In two of these, the read coverage along introns was related to nascent transcription in combination with co-transcriptional splicing e...

microRNAs (miRNAs) regulate the expression of mRNAs in animals and plants through miRNA-containing ribonucleoprotein particles (RNPs). At the core of these miRNA silencing effector complexes are the Argonaute (AGO) proteins that bind miRNAs and mediate target mRNA recognition. We generated HEK293 cell lines stably expressing epitope-tagged human AGO proteins and other RNA silencing-related proteins and used these cells to purify miRNA-containing RNPs. Mass spectrometric analyses of the proteins associated with different AGO proteins revealed a common set of helicases and mRNA-binding proteins, among them the three trinucleotide repeat containing proteins 6 (TNRC6A,-B,-C). mRNA microarray analyses of these miRNAassociated RNPs revealed that AGO and TNRC6 proteins bind highly similar sets of transcripts enriched in binding sites for highly expressed endogenous miRNAs, indicating that the TNRC6 proteins are a component of the mRNA-targeting miRNA silencing complex. Together with the very similar proteomic composition of each AGO complex, this result suggests substantial functional redundancy within families of human AGO and TNRC6 proteins. Our results further demonstrate that we have developed an effective biochemical approach to identify physiologically relevant human miRNA targets.

Oscillations are a key to achieving dynamic behavior and thus occur in biological systems as diverse as the beating heart, defecating worms, and nascent somites. Here we report pervasive, large-amplitude, and phase-locked oscillations of gene expression in developing C. elegans larvae, caused by periodic transcription. Nearly one fifth of detectably expressed transcripts oscillate with an 8 hr period, and hundreds change >10-fold. Oscillations are important for molting but occur in all phases, implying additional functions. Ribosome profiling reveals that periodic mRNA accumulation causes rhythmic translation, potentially facilitating transient protein accumulation as well as coordinated production of stable, complex structures such as the cuticle. Finally, large-amplitude oscillations in RNA sampled from whole worms indicate robust synchronization of gene expression programs across cells and tissues, suggesting that these oscillations will be a powerful new model to study coordinated gene expression in an animal.

Brain circuits are assembled from a large variety of morphologically and functionally diverse cell types. It is not known how the intermingled cell types of an individual adult brain region differ in their expressed genomes. Here we describe an atlas of cell type transcriptomes in one brain region, the mouse retina. We found that each adult cell type expressed a specific set of genes, including a unique set of transcription factors, forming a 'barcode' for cell identity. Cell type transcriptomes carried enough information to categorize cells into morphological classes and types. Several genes that were specifically expressed in particular retinal circuit elements, such as inhibitory neuron types, are associated with eye diseases. The resource described here allows gene expression to be compared across adult retinal cell types, experimenting with specific transcription factors to differentiate stem or somatic cells to retinal cell types, and predicting cellular targets of newly discovered disease-associated genes.