Perinuclear Anchoring of H3K9-Methylated Chromatin Stabilizes Induced Cell Fate in C. elegans Embryos.

Interphase chromatin is organized in distinct nuclear sub-compartments, reflecting its degree of compaction and transcriptional status. In Caenorhabditis elegans embryos, H3K9 methylation is necessary to silence and to anchor repeat-rich heterochromatin at the nuclear periphery. In a screen for perinuclear anchors of heterochromatin, we identified a previously uncharacterized C. elegans chromodomain protein, CEC-4. CEC-4 binds preferentially mono-, di-, or tri-methylated H3K9 and localizes at the nuclear envelope independently of H3K9 methylation and nuclear lamin. CEC-4 is necessary for endogenous heterochromatin anchoring, but not for transcriptional repression, in contrast to other known H3K9 methyl-binders in worms, which mediate gene repression but not perinuclear anchoring. When we ectopically induce a muscle differentiation program in embryos, cec-4 mutants fail to commit fully to muscle cell fate. This suggests that perinuclear sequestration of chromatin during development helps restrict cell differentiation programs by stabilizing commitment to a specific cell fate. PAPERCLIP.

Step-wise methylation of histone H3K9 positions heterochromatin at the nuclear periphery.

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.

Dynamics of miRNA accumulation during C. elegans larval development.

Temporally and spatially controlled accumulation underlies the functions of microRNAs (miRNAs) in various developmental processes. In Caenorhabditis elegans, this is exemplified by the temporal patterning miRNAs lin-4 and let-7, but for most miRNAs, developmental expression patterns remain poorly resolved. Indeed, experimentally observed long half-lives may constrain possible dynamics. Here, we profile miRNA expression throughout C. elegans postembryonic development at high temporal resolution, which identifies dynamically expressed miRNAs. We use mathematical models to explore the underlying mechanisms. For let-7, we can explain, and experimentally confirm, a striking stepwise accumulation pattern through a combination of rhythmic transcription and stage-specific regulation of precursor processing by the RNA-binding protein LIN-28. By contrast, the dynamics of several other miRNAs cannot be explained by regulation of production rates alone. Specifically, we show that a combination of oscillatory transcription and rhythmic decay drive rhythmic accumulation of miR-235, orthologous to miR-92 in other animals. We demonstrate that decay of miR-235 and additional miRNAs depends on EBAX-1, previously implicated in target-directed miRNA degradation (TDMD). Taken together, our results provide insight into dynamic miRNA decay and establish a resource to studying both the developmental functions of, and the regulatory mechanisms acting on, miRNAs.

Deep quantification of substrate turnover defines protease subsite cooperativity.

Substrate specificity determines protease functions in physiology and in clinical and biotechnological applications, yet quantitative cleavage information is often unavailable, biased, or limited to a small number of events. Here, we develop qPISA (quantitative Protease specificity Inference from Substrate Analysis) to study Dipeptidyl Peptidase Four (DPP4), a key regulator of blood glucose levels. We use mass spectrometry to quantify >40,000 peptides from a complex, commercially available peptide mixture. By analyzing changes in substrate levels quantitatively instead of focusing on qualitative product identification through a binary classifier, we can reveal cooperative interactions within DPP4's active pocket and derive a sequence motif that predicts activity quantitatively. qPISA distinguishes DPP4 from the related C. elegans DPF-3 (a DPP8/9-orthologue), and we relate the differences to the structural features of the two enzymes. We demonstrate that qPISA can direct protein engineering efforts like the stabilization of GLP-1, a key DPP4 substrate used in the treatment of diabetes and obesity. Thus, qPISA offers a versatile approach for profiling protease and especially exopeptidase specificity, facilitating insight into enzyme mechanisms and biotechnological and clinical applications.

Active chromatin marks drive spatial sequestration of heterochromatin in C. elegans nuclei.

The execution of developmental programs of gene expression requires an accurate partitioning of the genome into subnuclear compartments, with active euchromatin enriched centrally and silent heterochromatin at the nuclear periphery1. The existence of degenerative diseases linked to lamin A mutations suggests that perinuclear binding of chromatin contributes to cell-type integrity2,3. The methylation of lysine 9 of histone H3 (H3K9me) characterizes heterochromatin and mediates both transcriptional repression and chromatin anchoring at the inner nuclear membrane4. In Caenorhabditis elegans embryos, chromodomain protein CEC-4 bound to the inner nuclear membrane tethers heterochromatin through H3K9me3,5, whereas in differentiated tissues, a second heterochromatin-sequestering pathway is induced. Here we use an RNA interference screen in the cec-4 background and identify MRG-1 as a broadly expressed factor that is necessary for this second chromatin anchor in intestinal cells. However, MRG-1 is exclusively bound to euchromatin, suggesting that it acts indirectly. Heterochromatin detachment in double mrg-1; cec-4 mutants is rescued by depleting the histone acetyltransferase CBP-1/p300 or the transcription factor ATF-8, a member of the bZIP family (which is known to recruit CBP/p300). Overexpression of CBP-1 in cec-4 mutants is sufficient to delocalize heterochromatin in an ATF-8-dependent manner. CBP-1 and H3K27ac levels increase in heterochromatin upon mrg-1 knockdown, coincident with delocalization. This suggests that the spatial organization of chromatin in C. elegans is regulated both by the direct perinuclear attachment of silent chromatin, and by an active retention of CBP-1/p300 in euchromatin. The two pathways contribute differentially in embryos and larval tissues, with CBP-1 sequestration by MRG-1 having a major role in differentiated cells.

Genome-wide Single-Molecule Footprinting Reveals High RNA Polymerase II Turnover at Paused Promoters.

Transcription initiation entails chromatin opening followed by pre-initiation complex formation and RNA polymerase II recruitment. Subsequent polymerase elongation requires additional signals, resulting in increased residence time downstream of the start site, a phenomenon referred to as pausing. Here, we harnessed single-molecule footprinting to quantify distinct steps of initiation in vivo throughout the Drosophila genome. This identifies the impact of promoter structure on initiation dynamics in relation to nucleosomal occupancy. Additionally, perturbation of transcriptional initiation reveals an unexpectedly high turnover of polymerases at paused promoters-an observation confirmed at the level of nascent RNAs. These observations argue that absence of elongation is largely caused by premature termination rather than by stable polymerase stalling. In support of this non-processive model, we observe that induction of the paused heat shock promoter depends on continuous initiation. Our study provides a framework to quantify protein binding at single-molecule resolution and refines concepts of transcriptional pausing.

LIN41 Post-transcriptionally Silences mRNAs by Two Distinct and Position-Dependent Mechanisms.

The RNA-binding protein (RBP) LIN41, also known as LIN-41 or TRIM71, is a key regulator of animal development, but its physiological targets and molecular mechanism of action are largely elusive. Here we find that this RBP has two distinct mRNA-silencing activities. Using genome-wide ribosome profiling, RNA immunoprecipitation, and in vitro-binding experiments, we identify four mRNAs, each encoding a transcription factor or cofactor, as direct physiological targets of C. elegans LIN41. LIN41 silences three of these targets through their 3' UTRs, but it achieves isoform-specific silencing of one target, lin-29A, through its unique 5' UTR. Whereas the 3' UTR targets mab-10, mab-3, and dmd-3 undergo transcript degradation, lin-29A experiences translational repression. Through binding site transplantation experiments, we demonstrate that it is the location of the LIN41-binding site that specifies the silencing mechanism. Such position-dependent dual activity may, when studied more systematically, emerge as a feature shared by other RBPs.

Extensive oscillatory gene expression during C. elegans larval development.

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.

The CSR-1 endogenous RNAi pathway ensures accurate transcriptional reprogramming during the oocyte-to-embryo transition in Caenorhabditis elegans.

Endogenous RNAi (endoRNAi) is a conserved mechanism for fine-tuning gene expression. In the nematode Caenorhabditis elegans, several endoRNAi pathways are required for the successful development of reproductive cells. The CSR-1 endoRNAi pathway promotes germ cell development, primarily by facilitating the expression of germline genes. In this study, we report a novel function for the CSR-1 pathway in preventing premature activation of embryonic transcription in the developing oocytes, which is accompanied by a general Pol II activation. This CSR-1 function requires its RNase activity, suggesting that, by controlling the levels of maternal mRNAs, CSR-1-dependent endoRNAi contributes to an orderly reprogramming of transcription during the oocyte-to-embryo transition.

RNAi keeps Atf1-bound stress response genes in check at nuclear pores.

RNAi pathways are prevalent throughout the eukaryotic kingdom and are well known to regulate gene expression on a post-transcriptional level in the cytoplasm. Less is known about possible functions of RNAi in the nucleus. In the fission yeast Schizosaccharomyces pombe, RNAi is crucial to establish and maintain centromeric heterochromatin and functions to repress genome activity by a chromatin silencing mechanism referred to as cotranscriptional gene silencing (CTGS). Mechanistic details and the physiological relevance of CTGS are unknown. Here we show that RNAi components interact with chromatin at nuclear pores to keep stress response genes in check. We demonstrate that RNAi-mediated CTGS represses stress-inducible genes by degrading mRNAs under noninduced conditions. Under chronic heat stress conditions, a Dicer thermoswitch deports Dicer to the cytoplasm, thereby disrupting CTGS and enabling expression of genes implicated in the acquisition of thermotolerance. Taken together, our work highlights a role for nuclear pores and the stress response transcription factor Atf1 in coordinating the interplay between the RNAi machinery and the S. pombe genome and uncovers a novel mode of RNAi regulation in response to an environmental cue.

Potent degradation of neuronal miRNAs induced by highly complementary targets.

MicroRNAs (miRNAs) regulate target mRNAs by silencing them. Reciprocally, however, target mRNAs can also modulate miRNA stability. Here, we uncover a remarkable efficacy of target RNA-directed miRNA degradation (TDMD) in rodent primary neurons. Coincident with degradation, and while still bound to Argonaute, targeted miRNAs are 3' terminally tailed and trimmed. Absolute quantification of both miRNAs and their decay-inducing targets suggests that neuronal TDMD is multiple turnover and does not involve co-degradation of the target but rather competes with miRNA-mediated decay of the target. Moreover, mRNA silencing, but not TDMD, relies on cooperativity among multiple target sites to reach high efficacy. This knowledge can be harnessed for effective depletion of abundant miRNAs. Our findings bring insight into a potent miRNA degradation pathway in primary neurons, whose TDMD activity greatly surpasses that of non-neuronal cells and established cell lines. Thus, TDMD may be particularly relevant for miRNA regulation in the nervous system.

DNA-binding factors shape the mouse methylome at distal regulatory regions.

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.

DNA sequence explains seemingly disordered methylation levels in partially methylated domains of Mammalian genomes.

For the most part metazoan genomes are highly methylated and harbor only small regions with low or absent methylation. In contrast, partially methylated domains (PMDs), recently discovered in a variety of cell lines and tissues, do not fit this paradigm as they show partial methylation for large portions (20%-40%) of the genome. While in PMDs methylation levels are reduced on average, we found that at single CpG resolution, they show extensive variability along the genome outside of CpG islands and DNase I hypersensitive sites (DHS). Methylation levels range from 0% to 100% in a roughly uniform fashion with only little similarity between neighboring CpGs. A comparison of various PMD-containing methylomes showed that these seemingly disordered states of methylation are strongly conserved across cell types for virtually every PMD. Comparative sequence analysis suggests that DNA sequence is a major determinant of these methylation states. This is further substantiated by a purely sequence based model which can predict 31% (R(2)) of the variation in methylation. The model revealed CpG density as the main driving feature promoting methylation, opposite to what has been shown for CpG islands, followed by various dinucleotides immediately flanking the CpG and a minor contribution from sequence preferences reflecting nucleosome positioning. Taken together we provide a reinterpretation for the nucleotide-specific methylation levels observed in PMDs, demonstrate their conservation across tissues and suggest that they are mainly determined by specific DNA sequence features.

QuasR: quantification and annotation of short reads in R.

UNLABELLED: QuasR is a package for the integrated analysis of high-throughput sequencing data in R, covering all steps from read preprocessing, alignment and quality control to quantification. QuasR supports different experiment types (including RNA-seq, ChIP-seq and Bis-seq) and analysis variants (e.g. paired-end, stranded, spliced and allele-specific), and is integrated in Bioconductor so that its output can be directly processed for statistical analysis and visualization. AVAILABILITY AND IMPLEMENTATION: QuasR is implemented in R and C/C++. Source code and binaries for major platforms (Linux, OS X and MS Windows) are available from Bioconductor (www.bioconductor.org/packages/release/bioc/html/QuasR.html). The package includes a 'vignette' with step-by-step examples for typical work flows. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Engineering of a conditional allele reveals multiple roles of XRN2 in Caenorhabditis elegans development and substrate specificity in microRNA turnover.

Although XRN2 proteins are highly conserved eukaryotic 5'→3' exonucleases, little is known about their function in animals. Here, we characterize Caenorhabditis elegans XRN2, which we find to be a broadly and constitutively expressed nuclear protein. An xrn-2 null mutation or loss of XRN2 catalytic activity causes a molting defect and early larval arrest. However, by generating a conditionally mutant xrn-2ts strain de novo through an approach that may be also applicable to other genes of interest, we reveal further functions in fertility, during embryogenesis and during additional larval stages. Consistent with the known role of XRN2 in controlling microRNA (miRNA) levels, we can demonstrate that loss of XRN2 activity stabilizes some rapidly decaying miRNAs. Surprisingly, however, other miRNAs continue to decay rapidly in xrn-2ts animals. Thus, XRN2 has unanticipated miRNA specificity in vivo, and its diverse developmental functions may relate to distinct substrates. Finally, our global analysis of miRNA stability during larval stage 1 reveals that miRNA passenger strands (miR*s) are substantially less stable than guide strands (miRs), supporting the notion that the former are mostly byproducts of biogenesis rather than a less abundant functional species.

Functional characterization of C. elegans Y-box-binding proteins reveals tissue-specific functions and a critical role in the formation of polysomes.

The cold shock domain is one of the most highly conserved motifs between bacteria and higher eukaryotes. Y-box-binding proteins represent a subfamily of cold shock domain proteins with pleiotropic functions, ranging from transcription in the nucleus to translation in the cytoplasm. These proteins have been investigated in all major model organisms except Caenorhabditis elegans. In this study, we set out to fill this gap and present a functional characterization of CEYs, the C. elegans Y-box-binding proteins. We find that, similar to other organisms, CEYs are essential for proper gametogenesis. However, we also report a novel function of these proteins in the formation of large polysomes in the soma. In the absence of the somatic CEYs, polysomes are dramatically reduced with a simultaneous increase in monosomes and disomes, which, unexpectedly, has no obvious impact on animal biology. Because transcripts that are enriched in polysomes in wild-type animals tend to be less abundant in the absence of CEYs, our findings suggest that large polysomes might depend on transcript stabilization mediated by CEY proteins.

Transcription factor occupancy can mediate active turnover of DNA methylation at regulatory regions.

Distal regulatory elements, including enhancers, play a critical role in regulating gene activity. Transcription factor binding to these elements correlates with Low Methylated Regions (LMRs) in a process that is poorly understood. Here we ask whether and how actual occupancy of DNA-binding factors is linked to DNA methylation at the level of individual molecules. Using CTCF as an example, we observe that frequency of binding correlates with the likelihood of a demethylated state and sites of low occupancy display heterogeneous DNA methylation within the CTCF motif. In line with a dynamic model of binding and DNA methylation turnover, we find that 5-hydroxymethylcytosine (5hmC), formed as an intermediate state of active demethylation, is enriched at LMRs in stem and somatic cells. Moreover, a significant fraction of changes in 5hmC during differentiation occurs at these regions, suggesting that transcription factor activity could be a key driver for active demethylation. Since deletion of CTCF is lethal for embryonic stem cells, we used genetic deletion of REST as another DNA-binding factor implicated in LMR formation to test this hypothesis. The absence of REST leads to a decrease of hydroxymethylation and a concomitant increase of DNA methylation at its binding sites. These data support a model where DNA-binding factors can mediate turnover of DNA methylation as an integral part of maintenance and reprogramming of regulatory regions.

A quantitative RNA code for mRNA target selection by the germline fate determinant GLD-1.

RNA-binding proteins (RBPs) are critical regulators of gene expression. To understand and predict the outcome of RBP-mediated regulation a comprehensive analysis of their interaction with RNA is necessary. The signal transduction and activation of RNA (STAR) family of RBPs includes developmental regulators and tumour suppressors such as Caenorhabditis elegans GLD-1, which is a key regulator of germ cell development. To obtain a comprehensive picture of GLD-1 interactions with the transcriptome, we identified GLD-1-associated mRNAs by RNA immunoprecipitation followed by microarray detection. Based on the computational analysis of these mRNAs we generated a predictive model, where GLD-1 association with mRNA is determined by the strength and number of 7-mer GLD-1-binding motifs (GBMs) within UTRs. We verified this quantitative model both in vitro, by competition GLD-1/GBM-binding experiments to determine relative affinity, and in vivo, by 'transplantation' experiments, where 'weak' and 'strong' GBMs imposed translational repression of increasing strength on a non-target mRNA. This study demonstrates that transcriptome-wide identification of RBP mRNA targets combined with quantitative computational analysis can generate highly predictive models of post-transcriptional regulatory networks.

The mammalian TRIM-NHL protein TRIM71/LIN-41 is a repressor of mRNA function.

TRIM-NHL proteins are conserved regulators of development and differentiation but their molecular function has remained largely elusive. Here, we report an as yet unrecognized activity for the mammalian TRIM-NHL protein TRIM71 as a repressor of mRNAs. We show that TRIM71 is associated with mRNAs and that it promotes translational repression and mRNA decay. We have identified Rbl1 and Rbl2, two transcription factors whose down-regulation is important for stem cell function, as TRIM71 targets in mouse embryonic stem cells. Furthermore, one of the defining features of TRIM-NHL proteins, the NHL domain, is necessary and sufficient to target TRIM71 to RNA, while the RING domain that confers ubiquitin ligase activity is dispensable for repression. Our results reveal strong similarities between TRIM71 and Drosophila BRAT, the best-studied TRIM-NHL protein and a well-documented translational repressor, suggesting that BRAT and TRIM71 are part of a family of mRNA repressors regulating proliferation and differentiation.

Identification of active regulatory regions from DNA methylation data.

We have recently shown that transcription factor binding leads to defined reduction in DNA methylation, allowing for the identification of active regulatory regions from high-resolution methylomes. Here, we present MethylSeekR, a computational tool to accurately identify such footprints from bisulfite-sequencing data. Applying our method to a large number of published human methylomes, we demonstrate its broad applicability and generalize our previous findings from a neuronal differentiation system to many cell types and tissues. MethylSeekR is available as an R package at www.bioconductor.org.