asteRIa enables robust interaction modeling between chromatin modifications and epigenetic readers.

Chromatin, the nucleoprotein complex consisting of DNA and histone proteins, plays a crucial role in regulating gene expression by controlling access to DNA. Chromatin modifications are key players in this regulation, as they help to orchestrate DNA transcription, replication, and repair. These modifications recruit epigenetic 'reader' proteins, which mediate downstream events. Most modifications occur in distinctive combinations within a nucleosome, suggesting that epigenetic information can be encoded in combinatorial chromatin modifications. A detailed understanding of how multiple modifications cooperate in recruiting such proteins has, however, remained largely elusive. Here, we integrate nucleosome affinity purification data with high-throughput quantitative proteomics and hierarchical interaction modeling to estimate combinatorial effects of chromatin modifications on protein recruitment. This is facilitated by the computational workflow asteRIa which combines hierarchical interaction modeling, stability-based model selection, and replicate-consistency checks for a stable estimation of Robust Interactions among chromatin modifications. asteRIa identifies several epigenetic reader candidates responding to specific interactions between chromatin modifications. For the polycomb protein CBX8, we independently validate our results using genome-wide ChIP-Seq and bisulphite sequencing datasets. We provide the first quantitative framework for identifying cooperative effects of chromatin modifications on protein binding.

Decoding chromatin states by proteomic profiling of nucleosome readers.

DNA and histone modifications combine into characteristic patterns that demarcate functional regions of the genome1,2. While many 'readers' of individual modifications have been described3-5, how chromatin states comprising composite modification signatures, histone variants and internucleosomal linker DNA are interpreted is a major open question. Here we use a multidimensional proteomics strategy to systematically examine the interaction of around 2,000 nuclear proteins with over 80 modified dinucleosomes representing promoter, enhancer and heterochromatin states. By deconvoluting complex nucleosome-binding profiles into networks of co-regulated proteins and distinct nucleosomal features driving protein recruitment or exclusion, we show comprehensively how chromatin states are decoded by chromatin readers. We find highly distinctive binding responses to different features, many factors that recognize multiple features, and that nucleosomal modifications and linker DNA operate largely independently in regulating protein binding to chromatin. Our online resource, the Modification Atlas of Regulation by Chromatin States (MARCS), provides in-depth analysis tools to engage with our results and advance the discovery of fundamental principles of genome regulation by chromatin states.

Depletion of pyruvate kinase (PK) activity causes glycolytic intermediate imbalances and reveals a PK-TXNIP regulatory axis.

OBJECTIVE: Cancer cells convert more glucose into lactate than healthy cells, what contributes to their growth advantage. Pyruvate kinase (PK) is a key rate limiting enzyme in this process, what makes it a promising potential therapeutic target. However, currently it is still unclear what consequences the inhibition of PK has on cellular processes. Here, we systematically investigate the consequences of PK depletion for gene expression, histone modifications and metabolism. METHODS: Epigenetic, transcriptional and metabolic targets were analysed in different cellular and animal models with stable knockdown or knockout of PK. RESULTS: Depleting PK activity reduces the glycolytic flux and causes accumulation of glucose-6-phosphate (G6P). Such metabolic perturbation results in stimulation of the activity of a heterodimeric pair of transcription factors MondoA and MLX but not in a major reprogramming of the global H3K9ac and H3K4me3 histone modification landscape. The MondoA:MLX heterodimer upregulates expression of thioredoxin-interacting protein (TXNIP) - a tumour suppressor with multifaceted anticancer activity. This effect of TXNIP upregulation extends beyond immortalised cancer cell lines and is applicable to multiple cellular and animal models. CONCLUSIONS: Our work shows that actions of often pro-tumorigenic PK and anti-tumorigenic TXNIP are tightly linked via a glycolytic intermediate. We suggest that PK depletion stimulates the activity of MondoA:MLX transcription factor heterodimers and subsequently, increases cellular TXNIP levels. TXNIP-mediated inhibition of thioredoxin (TXN) can reduce the ability of cells to scavenge reactive oxygen species (ROS) leading to the oxidative damage of cellular structures including DNA. These findings highlight an important regulatory axis affecting tumour suppression mechanisms and provide an attractive opportunity for combination cancer therapies targeting glycolytic activity and ROS-generating pathways.

Nascent RNA antagonizes the interaction of a set of regulatory proteins with chromatin.

A number of regulatory factors are recruited to chromatin by specialized RNAs. Whether RNA has a more general role in regulating the interaction of proteins with chromatin has not been determined. We used proteomics methods to measure the global impact of nascent RNA on chromatin in embryonic stem cells. Surprisingly, we found that nascent RNA primarily antagonized the interaction of chromatin modifiers and transcriptional regulators with chromatin. Transcriptional inhibition and RNA degradation induced recruitment of a set of transcriptional regulators, chromatin modifiers, nucleosome remodelers, and regulators of higher-order structure. RNA directly bound to factors, including BAF, NuRD, EHMT1, and INO80 and inhibited their interaction with nucleosomes. The transcriptional elongation factor P-TEFb directly bound pre-mRNA, and its recruitment to chromatin upon Pol II inhibition was regulated by the 7SK ribonucleoprotein complex. We postulate that by antagonizing the interaction of regulatory proteins with chromatin, nascent RNA links transcriptional output with chromatin composition.

METTL6 is a tRNA m3C methyltransferase that regulates pluripotency and tumor cell growth.

Recently, covalent modifications of RNA, such as methylation, have emerged as key regulators of all aspects of RNA biology and have been implicated in numerous diseases, for instance, cancer. Here, we undertook a combination of in vitro and in vivo screens to test 78 potential methyltransferases for their roles in hepatocellular carcinoma (HCC) cell proliferation. We identified methyltransferase-like protein 6 (METTL6) as a crucial regulator of tumor cell growth. We show that METTL6 is a bona fide transfer RNA (tRNA) methyltransferase, catalyzing the formation of 3-methylcytidine at C32 of specific serine tRNA isoacceptors. Deletion of Mettl6 in mouse stem cells results in changes in ribosome occupancy and RNA levels, as well as impaired pluripotency. In mice, Mettl6 knockout results in reduced energy expenditure. We reveal a previously unknown pathway in the maintenance of translation efficiency with a role in maintaining stem cell self-renewal, as well as impacting tumor cell growth profoundly.

An Integrated Chemical Proteomics Approach for Quantitative Profiling of Intracellular ADP-Ribosylation.

ADP-ribosylation is integral to a diverse range of cellular processes such as DNA repair, chromatin regulation and RNA processing. However, proteome-wide investigation of its cellular functions has been limited due to numerous technical challenges including the complexity of the poly(ADP-ribose) (PAR) chains, low abundance of the modification and lack of sensitive enrichment methods. We herein show that an adenosine analogue with a terminal alkyne functionality at position 2 of the adenine (2-alkyne adenosine or 2YnAd) is suitable for selective enrichment, fluorescence detection and mass spectrometry proteomics analysis of the candidate ADP-ribosylome in mammalian cells. Although similar labelling profiles were observed via fluorescence imaging for 2YnAd and 6YnAd, a previously reported clickable NAD+ precursor, quantitative mass spectrometry analysis of the two probes in MDA-MB-231 breast cancer cells revealed a significant increase in protein coverage of the 2YnAd probe. To facilitate global enrichment of ADP-ribosylated proteins, we developed a dual metabolic labelling approach that involves simultaneous treatment of live cells with both 2YnAd and 6YnAd. By combining this dual metabolic labelling strategy with highly sensitive tandem mass tag (TMT) isobaric mass spectrometry and hierarchical Bayesian analysis, we have quantified the responses of thousands of endogenous proteins to clinical PARP inhibitors Olaparib and Rucaparib.

H4K20me0 recognition by BRCA1-BARD1 directs homologous recombination to sister chromatids.

Genotoxic DNA double-strand breaks (DSBs) can be repaired by error-free homologous recombination (HR) or mutagenic non-homologous end-joining1. HR supresses tumorigenesis1, but is restricted to the S and G2 phases of the cell cycle when a sister chromatid is present2. Breast cancer type 1 susceptibility protein (BRCA1) promotes HR by antagonizing the anti-resection factor TP53-binding protein 1(53BP1) (refs. 2-5), but it remains unknown how BRCA1 function is limited to the S and G2 phases. We show that BRCA1 recruitment requires recognition of histone H4 unmethylated at lysine 20 (H4K20me0), linking DSB repair pathway choice directly to sister chromatid availability. We identify the ankyrin repeat domain of BRCA1-associated RING domain protein 1 (BARD1)-the obligate BRCA1 binding partner3-as a reader of H4K20me0 present on new histones in post-replicative chromatin6. BARD1 ankyrin repeat domain mutations disabling H4K20me0 recognition abrogate accumulation of BRCA1 at DSBs, causing aberrant build-up of 53BP1, and allowing anti-resection activity to prevail in S and G2. Consequently, BARD1 recognition of H4K20me0 is required for HR and resistance to poly (ADP-ribose) polymerase inhibitors. Collectively, this reveals that BRCA1-BARD1 monitors the replicative state of the genome to oppose 53BP1 function, routing only DSBs within sister chromatids to HR.

MEANS: python package for Moment Expansion Approximation, iNference and Simulation.

MOTIVATION: Many biochemical systems require stochastic descriptions. Unfortunately these can only be solved for the simplest cases and their direct simulation can become prohibitively expensive, precluding thorough analysis. As an alternative, moment closure approximation methods generate equations for the time-evolution of the system's moments and apply a closure ansatz to obtain a closed set of differential equations; that can become the basis for the deterministic analysis of the moments of the outputs of stochastic systems. RESULTS: We present a free, user-friendly tool implementing an efficient moment expansion approximation with parametric closures that integrates well with the IPython interactive environment. Our package enables the analysis of complex stochastic systems without any constraints on the number of species and moments studied and the type of rate laws in the system. In addition to the approximation method our package provides numerous tools to help non-expert users in stochastic analysis. AVAILABILITY AND IMPLEMENTATION: https://github.com/theosysbio/means CONTACTS: m.stumpf@imperial.ac.uk or e.lakatos13@imperial.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DGW: an exploratory data analysis tool for clustering and visualisation of epigenomic marks.

BACKGROUND: Functional genomic and epigenomic research relies fundamentally on sequencing based methods like ChIP-seq for the detection of DNA-protein interactions. These techniques return large, high dimensional data sets with visually complex structures, such as multi-modal peaks extended over large genomic regions. Current tools for visualisation and data exploration represent and leverage these complex features only to a limited extent. RESULTS: We present DGW, an open source software package for simultaneous alignment and clustering of multiple epigenomic marks. DGW uses Dynamic Time Warping to adaptively rescale and align genomic distances which allows to group regions of interest with similar shapes, thereby capturing the structure of epigenomic marks. We demonstrate the effectiveness of the approach in a simulation study and on a real epigenomic data set from the ENCODE project. CONCLUSIONS: Our results show that DGW automatically recognises and aligns important genomic features such as transcription start sites and splicing sites from histone marks. DGW is available as an open source Python package.