Total kinetic analysis reveals how combinatorial methylation patterns are established on lysines 27 and 36 of histone H3

Crosstalk within and beyond the Polycomb repressive system

The development of multicellular organisms depends on spatiotemporally controlled differentiation of numerous cell types and their maintenance. To generate such diversity based on the invariant genetic information stored in DNA, epigenetic mechanisms, which are heritable changes in gene function that do not involve alterations to the underlying DNA sequence, are required to establish and maintain unique gene expression programs. Polycomb repressive complexes represent a paradigm of epigenetic regulation of developmentally regulated genes, and the roles of these complexes as well as the epigenetic marks they deposit, namely H3K27me3 and H2AK119ub, have been extensively studied. However, an emerging theme from recent studies is that not only the autonomous functions of the Polycomb repressive system, but also crosstalks of Polycomb with other epigenetic modifications, are important for gene regulation. In this review, we summarize how these crosstalk mechanisms have improved our understanding of Polycomb biology and how such knowledge could help with the design of cancer treatments that target the dysregulated epigenome.

Polycomb repressive complex 2 and its core component EZH2: potential targeted therapeutic strategies for head and neck squamous cell carcinoma

The polycomb group (PcG) comprises a set of proteins that exert epigenetic regulatory effects and play crucial roles in diverse biological processes, ranging from pluripotency and development to carcinogenesis. Among these proteins, enhancer of zeste homolog 2 (EZH2) stands out as a catalytic component of polycomb repressive complex 2 (PRC2), which plays a role in regulating the expression of homologous (Hox) genes and initial stages of x chromosome inactivation. In numerous human cancers, including head and neck squamous cell carcinoma (HNSCC), EZH2 is frequently overexpressed or activated and has been identified as a negative prognostic factor. Notably, EZH2 emerges as a significant gene involved in regulating the STAT3/HOTAIR axis, influencing HNSCC proliferation, differentiation, and promoting metastasis by modulating related oncogenes in oral cancer. Currently, various small molecule compounds have been developed as inhibitors specifically targeting EZH2 and have gained approval for treating refractory tumors. In this review, we delve into the epigenetic regulation mediated by EZH2/PRC2 in HNSCC, with a specific focus on exploring the potential roles and mechanisms of EZH2, its crucial contribution to targeted drug therapy, and its association with cancer markers and epithelial-mesenchymal transition. Furthermore, we aim to unravel its potential as a therapeutic strategy for oral squamous cell carcinoma.

αKG-mediated carnitine synthesis promotes homologous recombination via histone acetylation

Homologous recombination (HR) deficiency enhances sensitivity to DNA damaging agents commonly used to treat cancer. In HR-proficient cancers, metabolic mechanisms driving response or resistance to DNA damaging agents remain unclear. Here we identified that depletion of alpha-ketoglutarate (αKG) sensitizes HR-proficient cells to DNA damaging agents by metabolic regulation of histone acetylation. αKG is required for the activity of αKG-dependent dioxygenases (αKGDDs), and prior work has shown that changes in αKGDD affect demethylases. Using a targeted CRISPR knockout library consisting of 64 αKGDDs, we discovered that Trimethyllysine Hydroxylase Epsilon (TMLHE), the first and rate-limiting enzyme in de novo carnitine synthesis, is necessary for proliferation of HR-proficient cells in the presence of DNA damaging agents. Unexpectedly, αKG-mediated TMLHE-dependent carnitine synthesis was required for histone acetylation, while histone methylation was affected but dispensable. The increase in histone acetylation via αKG-dependent carnitine synthesis promoted HR-mediated DNA repair through site- and substrate-specific histone acetylation. These data demonstrate for the first time that HR-proficiency is mediated through αKG directly influencing histone acetylation via carnitine synthesis and provide a metabolic avenue to induce HR-deficiency and sensitivity to DNA damaging agents.

Epigenetic Causes of Overgrowth Syndromes

Human overgrowth disorders are characterized by excessive prenatal and/or postnatal growth of various tissues. These disorders often present with tall stature, macrocephaly, and/or abdominal organomegaly and are sometimes associated with additional phenotypic abnormalities such as intellectual disability and increased cancer risk. As the genetic etiology of these disorders have been elucidated, a surprising pattern has emerged. Multiple monogenic overgrowth syndromes result from variants in epigenetic regulators: variants in histone methyltransferases NSD1 and EZH2 cause Sotos syndrome and Weaver syndrome, respectively, variants in DNA methyltransferase DNMT3A cause Tatton-Brown-Rahman syndrome, and variants in chromatin remodeler CHD8 cause an autism spectrum disorder with overgrowth. In addition, very recently, a variant in histone reader protein SPIN4 was identified in a new X-linked overgrowth disorder. In this review, we discuss the genetics of these overgrowth disorders and explore possible common underlying mechanisms by which epigenetic pathways regulate human body size.

Discovery of a highly potent, selective, orally bioavailable inhibitor of KAT6A/B histone acetyltransferases with efficacy against KAT6A-high ER+ breast cancer

KAT6A, and its paralog KAT6B, are histone lysine acetyltransferases (HAT) that acetylate histone H3K23 and exert an oncogenic role in several tumor types including breast cancer where KAT6A is frequently amplified/overexpressed. However, pharmacologic targeting of KAT6A to achieve therapeutic benefit has been a challenge. Here we describe identification of a highly potent, selective, and orally bioavailable KAT6A/KAT6B inhibitor CTx-648 (PF-9363), derived from a benzisoxazole series, which demonstrates anti-tumor activity in correlation with H3K23Ac inhibition in KAT6A over-expressing breast cancer. Transcriptional and epigenetic profiling studies show reduced RNA Pol II binding and downregulation of genes involved in estrogen signaling, cell cycle, Myc and stem cell pathways associated with CTx-648 anti-tumor activity in ER-positive (ER+) breast cancer. CTx-648 treatment leads to potent tumor growth inhibition in ER+ breast cancer in vivo models, including models refractory to endocrine therapy, highlighting the potential for targeting KAT6A in ER+ breast cancer.

Iron drives anabolic metabolism through active histone demethylation and mTORC1

All eukaryotic cells require a minimal iron threshold to sustain anabolic metabolism. However, the mechanisms by which cells sense iron to regulate anabolic processes are unclear. Here we report a previously undescribed eukaryotic pathway for iron sensing in which molecular iron is required to sustain active histone demethylation and maintain the expression of critical components of the pro-anabolic mTORC1 pathway. Specifically, we identify the iron-binding histone-demethylase KDM3B as an intrinsic iron sensor that regulates mTORC1 activity by demethylating H3K9me2 at enhancers of a high-affinity leucine transporter, LAT3, and RPTOR. By directly suppressing leucine availability and RAPTOR levels, iron deficiency supersedes other nutrient inputs into mTORC1. This process occurs in vivo and is not an indirect effect by canonical iron-utilizing pathways. Because ancestral eukaryotes share homologues of KDMs and mTORC1 core components, this pathway probably pre-dated the emergence of the other kingdom-specific nutrient sensors for mTORC1.

Likelihood-ratio test statistic for the finite-sample case in nonlinear ordinary differential equation models

Likelihood ratios are frequently utilized as basis for statistical tests, for model selection criteria and for assessing parameter and prediction uncertainties, e.g. using the profile likelihood. However, translating these likelihood ratios into p-values or confidence intervals requires the exact form of the test statistic's distribution. The lack of knowledge about this distribution for nonlinear ordinary differential equation (ODE) models requires an approximation which assumes the so-called asymptotic setting, i.e. a sufficiently large amount of data. Since the amount of data from quantitative molecular biology is typically limited in applications, this finite-sample case regularly occurs for mechanistic models of dynamical systems, e.g. biochemical reaction networks or infectious disease models. Thus, it is unclear whether the standard approach of using statistical thresholds derived for the asymptotic large-sample setting in realistic applications results in valid conclusions. In this study, empirical likelihood ratios for parameters from 19 published nonlinear ODE benchmark models are investigated using a resampling approach for the original data designs. Their distributions are compared to the asymptotic approximation and statistical thresholds are checked for conservativeness. It turns out, that corrections of the likelihood ratios in such finite-sample applications are required in order to avoid anti-conservative results.

Realistic simulation of time-course measurements in systems biology

In systems biology, the analysis of complex nonlinear systems faces many methodological challenges. For the evaluation and comparison of the performances of novel and competing computational methods, one major bottleneck is the availability of realistic test problems. We present an approach for performing realistic simulation studies for analyses of time course data as they are typically measured in systems biology. Since the design of experiments in practice depends on the process of interest, our approach considers the size and the dynamics of the mathematical model which is intended to be used for the simulation study. To this end, we used 19 published systems biology models with experimental data and evaluated the relationship between model features (e.g., the size and the dynamics) and features of the measurements such as the number and type of observed quantities, the number and the selection of measurement times, and the magnitude of measurement errors. Based on these typical relationships, our novel approach enables suggestions of realistic simulation study designs in the systems biology context and the realistic generation of simulated data for any dynamic model. The approach is demonstrated on three models in detail and its performance is validated on nine models by comparing ODE integration, parameter optimization, and parameter identifiability. The presented approach enables more realistic and less biased benchmark studies and thereby constitutes an important tool for the development of novel methods for dynamic modeling.

Distinct roles for canonical and variant histone H3 lysine-36 in Polycomb silencing

Polycomb complexes regulate cell type-specific gene expression programs through heritable silencing of target genes. Trimethylation of histone H3 lysine 27 (H3K27me3) is essential for this process. Perturbation of H3K36 is thought to interfere with H3K27me3. We show that mutants of Drosophila replication-dependent (H3.2K36R) or replication-independent (H3.3K36R) histone H3 genes generally maintain Polycomb silencing and reach later stages of development. In contrast, combined (H3.3K36RH3.2K36R) mutants display widespread Hox gene misexpression and fail to develop past the first larval stage. Chromatin profiling revealed that the H3.2K36R mutation disrupts H3K27me3 levels broadly throughout silenced domains, whereas these regions are mostly unaffected in H3.3K36R animals. Analysis of H3.3 distributions showed that this histone is enriched at presumptive Polycomb response elements located outside of silenced domains but relatively depleted from those inside. We conclude that H3.2 and H3.3 K36 residues collaborate to repress Hox genes using different mechanisms.

Antecedent chromatin organization determines cGAS recruitment to ruptured micronuclei

Micronuclei (MN) are cytosolic bodies that sequester acentric fragments or mis-segregated chromosomes from the primary nucleus. Spontaneous rupture of the MN envelope allows recognition by the viral receptor cyclic GMP-AMP synthase (cGAS), initiating interferon signaling downstream of DNA damage. Here, we demonstrate that MN rupture is permissive but not sufficient for cGAS localization. Chromatin characteristics such as histone 3, lysine 79 dimethylation (H3K79me2) are present in the nucleus before DNA damage, retained in ruptured MN, and regulate cGAS recruitment. cGAS is further responsive to dynamic intra-MN processes occurring prior to rupture, including transcription. MN chromatin tethering via the nucleosome acidic patch is necessary for cGAS-dependent interferon signaling. Our data suggest that both damage-antecedent nuclear chromatin status and MN-contained chromatin organizational changes dictate cGAS recruitment and the magnitude of the cGAS-driven interferon cascade. Our work defines MN as integrative signaling hubs for the cellular response to genotoxic stress.

Polycomb Repressive Complex 2 in Oncology

Dynamic regulation of the chromatin state by Polycomb Repressive Complex 2 (PRC2) provides an important mean for epigenetic gene control that can profoundly influence normal development and cell lineage specification. PRC2 and PRC2-induced methylation of histone H3 lysine 27 (H3K27) are critically involved in a wide range of DNA-templated processes, which at least include transcriptional repression and gene imprinting, organization of three-dimensional chromatin structure, DNA replication and DNA damage response and repair. PRC2-based genome regulation often goes wrong in diseases, notably cancer. This chapter discusses about different modes-of-action through which PRC2 and EZH2, a catalytic subunit of PRC2, mediate (epi)genomic and transcriptomic regulation. We will also discuss about how alteration or mutation of the PRC2 core or axillary component promotes oncogenesis, how post-translational modification regulates functionality of EZH2 and PRC2, and how PRC2 and other epigenetic pathways crosstalk. Lastly, we will briefly touch on advances in targeting EZH2 and PRC2 dependence as cancer therapeutics.

Efficient computation of adjoint sensitivities at steady-state in ODE models of biochemical reaction networks

Dynamical models in the form of systems of ordinary differential equations have become a standard tool in systems biology. Many parameters of such models are usually unknown and have to be inferred from experimental data. Gradient-based optimization has proven to be effective for parameter estimation. However, computing gradients becomes increasingly costly for larger models, which are required for capturing the complex interactions of multiple biochemical pathways. Adjoint sensitivity analysis has been pivotal for working with such large models, but methods tailored for steady-state data are currently not available. We propose a new adjoint method for computing gradients, which is applicable if the experimental data include steady-state measurements. The method is based on a reformulation of the backward integration problem to a system of linear algebraic equations. The evaluation of the proposed method using real-world problems shows a speedup of total simulation time by a factor of up to 4.4. Our results demonstrate that the proposed approach can achieve a substantial improvement in computation time, in particular for large-scale models, where computational efficiency is critical.

Understanding the Roles of the NSD Protein Methyltransferases in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is the sixth most prevalent non-skin cancer in the world. While immunotherapy has revolutionized the standard of care treatment in patients with recurrent/metastatic HNSCC, more than 70% of patients do not respond to this treatment, making the identification of novel therapeutic targets urgent. Recently, research endeavors have focused on how epigenetic modifications may affect tumor initiation and progression of HNSCC. The nuclear receptor binding SET domain (NSD) family of protein methyltransferases NSD1-NSD3 is of particular interest for HNSCC, with NSD1 and NSD3 being amongst the most commonly mutated or amplified genes respectively in HNSCC. Preclinical studies have identified both oncogenic and tumor-suppressing properties across NSD1, NSD2, and NSD3 within the context of HNSCC. The purpose of this review is to provide a better understanding of the contribution of the NSD family of protein methyltransferases to the pathogenesis of HNSCC, underscoring their promise as novel therapeutic targets in this devastating disease.

Multi-omic profiling of histone variant H3.3 lysine 27 methylation reveals a distinct role from canonical H3 in stem cell differentiation

Histone variants, such as histone H3.3, replace canonical histones within the nucleosome to alter chromatin accessibility and gene expression. Although the biological roles of selected histone post-translational modifications (PTMs) have been extensively characterized, the potential differences in the function of a given PTM on different histone variants is almost always elusive. By applying proteomics and genomics techniques, we investigate the role of lysine 27 tri-methylation specifically on the histone variant H3.3 (H3.3K27me3) in the context of mouse embryonic stem cell pluripotency and differentiation as a model system for development. We demonstrate that while the steady state overall levels of methylation on both H3K27 and H3.3K27 decrease during differentiation, methylation dynamics studies indicate that methylation on H3.3K27 is maintained more than on H3K27. Using a custom-made antibody, we identify a unique enrichment of H3.3K27me3 at lineage-specific genes, such as olfactory receptor genes, and at binding motifs for the transcription factors FOXJ2/3. REST, a predicted FOXJ2/3 target that acts as a transcriptional repressor of terminal neuronal genes, was identified with H3.3K27me3 at its promoter region. H3.3K27A mutant cells confirmed an upregulation of FOXJ2/3 targets upon the loss of methylation at H3.3K27. Thus, while canonical H3K27me3 has been characterized to regulate the expression of transcription factors that play a general role in differentiation, our work suggests H3.3K27me3 is essential for regulating distinct terminal differentiation genes. This work highlights the importance of understanding the effects of PTMs not only on canonical histones but also on specific histone variants, as they may exhibit distinct roles.

A synthetic lethality screen reveals ING5 as a genetic dependency of catalytically dead Set1A/COMPASS in mouse embryonic stem cells

Embryonic stem cells (ESCs) are defined by their ability to self-renew and the potential to differentiate into all tissues of the developing organism. We previously demonstrated that deleting the catalytic SET domain of the Set1A/complex of proteins associated with SET1 histone methyltransferase (Set1A/COMPASS) in mouse ESCs does not impair their viability or ability to self-renew; however, it leads to defects in differentiation. The precise mechanisms by which Set1A executes these functions remain to be elucidated. In this study, we demonstrate that mice lacking the SET domain of Set1A are embryonic lethal at a stage that is unique from null alleles. To gain insight into Set1A function in regulating pluripotency, we conducted a CRISPR/Cas9-mediated dropout screen and identified the MOZ/MORF (monocytic leukaemia zinc finger protein/monocytic leukaemia zinc finger protein-related factor) and HBO1 (HAT bound to ORC1) acetyltransferase complex member ING5 as a synthetic perturbation to Set1A. The loss of Ing5 in Set1AΔSET mouse ESCs decreases the fitness of these cells, and the simultaneous loss of ING5 and in Set1AΔSET leads to up-regulation of differentiation-associated genes. Taken together, our results point toward Set1A/COMPASS and ING5 as potential coregulators of the self-renewal and differentiation status of ESCs.

Precision Methylome and In Vivo Methylation Kinetics Characterization of Klebsiella pneumoniae

Klebsiella pneumoniae (K. pneumoniae) is an important pathogen that can cause severe hospital- and community-acquired infections. To systematically investigate its methylation features, we determined the whole-genome sequences of 14 K. pneumoniae strains covering varying serotypes, multilocus sequence types, clonal groups, viscosity/virulence, and drug resistance. Their methylomes were further characterized using Pacific Biosciences single-molecule real-time and bisulfite technologies. We identified 15 methylation motifs [13 N6-methyladenine (6mA) and two 5-methylcytosine (5mC) motifs], among which eight were novel. Their corresponding DNA methyltransferases were also validated. Additionally, we analyzed the genomic distribution of GATC and CCWGG methylation motifs shared by all strains, and identified differential distribution patterns of some hemi-/un-methylated GATC motifs, which tend to be located within intergenic regions (IGRs). Specifically, we characterized the in vivo methylation kinetics at single-base resolution on a genome-wide scale by simulating the dynamic processes of replication-mediated passive demethylation and MTase-catalyzed re-methylation. The slow methylation of the GATC motifs in the replication origin (oriC) regions and IGRs implicates the epigenetic regulation of replication initiation and transcription. Our findings illustrate the first comprehensive dynamic methylome map of K. pneumoniae at single-base resolution, and provide a useful reference to better understand epigenetic regulation in this and other bacterial species.

Activation of MAT2A-ACSL3 pathway protects cells from ferroptosis in gastric cancer

BACKGROUND

Ferroptosis, a unique form of nonapoptotic-regulated cell death caused by overwhelming lipid peroxidation, represents an emerging tumor suppression mechanism. Growing evidence has demonstrated that cell metabolism plays an important role in the regulation of ferroptosis. Specifically, the association between methionine metabolism and ferroptosis remains undefined.

METHODS

We performed in vitro and in vivo experiments to evaluate the influence of methionine metabolism on ferroptosis sensitivity. Pharmacological and genetic blockade of the methionine cycle was utilized and relevant molecular analyses were performed.

RESULTS

We identified MAT2A as a driver of ferroptosis resistance. Mechanistically, MAT2A mediates the production of S-adenosylmethionine (SAM), which upregulates ACSL3 by increasing the trimethylation of lysine-4 on histone H3 (H3K4me3) at the promoter area, resulting in ferroptosis resistance.

CONCLUSIONS

Collectively, these results established a link between methionine cycle activity and ferroptosis vulnerability in gastric cancer.

Interplay between chromatin marks in development and disease

DNA methylation (DNAme) and histone post-translational modifications (PTMs) have important roles in transcriptional regulation. Although many reports have characterized the functions of such chromatin marks in isolation, recent genome-wide studies reveal surprisingly complex interactions between them. Here, we focus on the interplay between DNAme and methylation of specific lysine residues on the histone H3 tail. We describe the impact of genetic perturbation of the relevant methyltransferases in the mouse on the landscape of chromatin marks as well as the transcriptome. In addition, we discuss the specific neurodevelopmental growth syndromes and cancers resulting from pathogenic mutations in the human orthologues of these genes. Integrating these observations underscores the fundamental importance of crosstalk between DNA and histone H3 methylation in development and disease.

High fidelity epigenetic inheritance: Information theoretic model predicts threshold filling of histone modifications post replication

During cell devision, maintaining the epigenetic information encoded in histone modification patterns is crucial for survival and identity of cells. The faithful inheritance of the histone marks from the parental to the daughter strands is a puzzle, given that each strand gets only half of the parental nucleosomes. Mapping DNA replication and reconstruction of modifications to equivalent problems in communication of information, we ask how well enzymes can recover the parental modifications, if they were ideal computing machines. Studying a parameter regime where realistic enzymes can function, our analysis predicts that enzymes may implement a critical threshold filling algorithm which fills unmodified regions of length at most k. This algorithm, motivated from communication theory, is derived from the maximum à posteriori probability (MAP) decoding which identifies the most probable modification sequence based on available observations. Simulations using our method produce modification patterns similar to what has been observed in recent experiments. We also show that our results can be naturally extended to explain inheritance of spatially distinct antagonistic modifications.

A dynamic and combinatorial histone code drives malaria parasite asexual and sexual development

Histone posttranslational modifications (PTMs) frequently co-occur on the same chromatin domains or even in the same molecule. It is now established that these “histone codes” are the result of cross talk between enzymes that catalyze multiple PTMs with univocal readout as compared with these PTMs in isolation. Here, we performed a comprehensive identification and quantification of histone codes of the malaria parasite, Plasmodium falciparum. We used advanced quantitative middle-down proteomics to identify combinations of PTMs in both the proliferative, asexual stages and transmissible, sexual gametocyte stages of P. falciparum. We provide an updated, high-resolution compendium of 77 PTMs on H3 and H3.3, of which 34 are newly identified in P. falciparum. Coexisting PTMs with unique stage distinctions were identified, indicating that many of these combinatorial PTMs are associated with specific stages of the parasite life cycle. We focused on the code H3R17me2K18acK23ac for its unique presence in mature gametocytes; chromatin proteomics identified a gametocyte-specific SAGA-like effector complex including the transcription factor AP2-G2, which we tied to this specific histone code, as involved in regulating gene expression in mature gametocytes. Ultimately, this study unveils previously undiscovered histone PTMs and their functional relationship with coexisting partners. These results highlight that investigating chromatin regulation in the parasite using single histone PTM assays might overlook higher-order gene regulation for distinct proliferation and differentiation processes.

•

First middle-down chromatin proteomics compendium of the malaria parasite, Plasmodium falciparum.

•

Novel histone PTMs (including arginine methylation) in both asexual parasites and transmissible gametocytes.

•

Histone PTM cross talk is dynamic life cycle stage stratified.

•

Gametocytes rely on histone PTM connectivity to allow onward transmission.

•

AP2-G2 is an important effector of H3K18acK23ac in mature gametocytes.

Mini-batch optimization enables training of ODE models on large-scale datasets

Quantitative dynamic models are widely used to study cellular signal processing. A critical step in modelling is the estimation of unknown model parameters from experimental data. As model sizes and datasets are steadily growing, established parameter optimization approaches for mechanistic models become computationally extremely challenging. Mini-batch optimization methods, as employed in deep learning, have better scaling properties. In this work, we adapt, apply, and benchmark mini-batch optimization for ordinary differential equation (ODE) models, thereby establishing a direct link between dynamic modelling and machine learning. On our main application example, a large-scale model of cancer signaling, we benchmark mini-batch optimization against established methods, achieving better optimization results and reducing computation by more than an order of magnitude. We expect that our work will serve as a first step towards mini-batch optimization tailored to ODE models and enable modelling of even larger and more complex systems than what is currently possible.

PRC2 Inhibitors Overcome Glucocorticoid Resistance Driven by NSD2 Mutation in Pediatric Acute Lymphoblastic Leukemia

Mutations in epigenetic regulators are common in relapsed pediatric acute lymphoblastic leukemia (ALL). Here, we uncovered the mechanism underlying the relapse of ALL driven by an activating mutation of the NSD2 histone methyltransferase (p.E1099K). Using high-throughput drug screening, we found that NSD2-mutant cells were specifically resistant to glucocorticoids. Correction of this mutation restored glucocorticoid sensitivity. The transcriptional response to glucocorticoids was blocked in NSD2-mutant cells due to depressed glucocorticoid receptor (GR) levels and the failure of glucocorticoids to autoactivate GR expression. Although H3K27me3 was globally decreased by NSD2 p.E1099K, H3K27me3 accumulated at the NR3C1 (GR) promoter. Pretreatment of NSD2 p.E1099K cell lines and patient-derived xenograft samples with PRC2 inhibitors reversed glucocorticoid resistance in vitro and in vivo. PRC2 inhibitors restored NR3C1 autoactivation by glucocorticoids, increasing GR levels and allowing GR binding and activation of proapoptotic genes. These findings suggest a new therapeutic approach to relapsed ALL associated with NSD2 mutation. SIGNIFICANCE: NSD2 histone methyltransferase mutations observed in relapsed pediatric ALL drove glucocorticoid resistance by repression of the GR and abrogation of GR gene autoactivation due to accumulation of K3K27me3 at its promoter. Pretreatment with PRC2 inhibitors reversed resistance, suggesting a new therapeutic approach to these patients with ALL.This article is highlighted in the In This Issue feature, p. 1.

Host Plant Adaptation Drives Changes in Diaphorina citri Proteome Regulation, Proteoform Expression, and Transmission of 'Candidatus Liberibacter asiaticus', the Citrus Greening Pathogen

The Asian citrus psyllid (Diaphorina citri) is a pest of citrus and the primary insect vector of the bacterial pathogen, 'Candidatus Liberibacter asiaticus' (CLas), which is associated with citrus greening disease. The citrus relative Murraya paniculata (orange jasmine) is a host plant of D. citri but is more resistant to CLas compared with all tested Citrus genotypes. The effect of host switching of D. citri between Citrus medica (citron) and M. paniculata plants on the acquisition and transmission of CLas was investigated. The psyllid CLas titer and the proportion of CLas-infected psyllids decreased in the generations after transfer from CLas-infected citron to healthy M. paniculata plants. Furthermore, after several generations of feeding on M. paniculata, pathogen acquisition (20 to 40% reduction) and transmission rates (15 to 20% reduction) in psyllids transferred to CLas-infected citron were reduced compared with psyllids continually maintained on infected citron. Top-down (difference gel electrophoresis) and bottom-up (shotgun MS/MS) proteomics methods were used to identify changes in D. citri protein expression resulting from host plant switching between Citrus macrophylla and M. paniculata. Changes in expression of insect metabolism, immunity, and cytoskeleton proteins were associated with host plant switching. Both transient and sustained feeding on M. paniculata induced distinct patterns of protein expression in D. citri compared with psyllids reared on C. macrophylla. The results point to complex interactions that affect vector competence and may lead to strategies to control the spread of citrus greening disease.

H3K36 trimethylation-mediated biological functions in cancer

Histone modification is an important form of epigenetic regulation. Thereinto, histone methylation is a critical determination of chromatin states, participating in multiple cellular processes. As a conserved histone methylation mark, histone 3 lysine 36 trimethylation (H3K36me3) can mediate multiple transcriptional-related events, such as the regulation of transcriptional activity, transcription elongation, pre-mRNA alternative splicing, and RNA m⁶A methylation. Additionally, H3K36me3 also contributes to DNA damage repair. Given the crucial function of H3K36me3 in genome regulation, the roles of H3K36me3 and its sole methyltransferase SETD2 in pathogenesis, especially malignancies, have been emphasized in many studies, and it is conceivable that disruption of histone methylation regulatory network composed of "writer", "eraser", "reader", and the mutation of H3K36me3 codes have the capacity of powerfully modulating cancer initiation and development. Here we review H3K36me3-mediated biological processes and summarize the latest findings regarding its role in cancers. We highlight the significance of epigenetic combination therapies in cancers.

Natural isotope correction improves analysis of protein modification dynamics

Stable isotope labelling in combination with high-resolution mass spectrometry approaches are increasingly used to analyze both metabolite and protein modification dynamics. To enable correct estimation of the resulting dynamics, it is critical to correct the measured values for naturally occurring stable isotopes, a process commonly called isotopologue correction or deconvolution. While the importance of isotopologue correction is well recognized in metabolomics, it has received far less attention in proteomics approaches. Although several tools exist that enable isotopologue correction of mass spectrometry data, the majority is tailored for the analysis of low molecular weight metabolites. We here present PICor which has been developed for isotopologue correction of complex isotope labelling experiments in proteomics or metabolomics and demonstrate the importance of appropriate correction for accurate determination of protein modifications dynamics, using histone acetylation as an example.

Combined Metabolic and Chemical (CoMetChem) Labeling Using Stable Isotopes-a Strategy to Reveal Site-Specific Histone Acetylation and Deacetylation Rates by LC-MS

Histone acetylation is an important, reversible post-translational protein modification and a hallmark of epigenetic regulation. However, little is known about the dynamics of this process, due to the lack of analytical methods that can capture site-specific acetylation and deacetylation reactions. We present a new approach that combines metabolic and chemical labeling (CoMetChem) using uniformly 13C-labeled glucose and stable isotope-labeled acetic anhydride. Thereby, chemically equivalent, fully acetylated histone species are generated, enabling accurate relative quantification of site-specific lysine acetylation dynamics in tryptic peptides using high-resolution mass spectrometry. We show that CoMetChem enables site-specific quantification of the incorporation or loss of lysine acetylation over time, allowing the determination of reaction rates for acetylation and deacetylation. Thus, the CoMetChem methodology provides a comprehensive description of site-specific acetylation dynamics.

DOT1L modulates the senescence-associated secretory phenotype through epigenetic regulation of IL1A

Oncogene-induced senescence (OIS) is a stable cell cycle arrest that occurs in normal cells upon oncogene activation. Cells undergoing OIS express a wide variety of secreted factors that affect the senescent microenvironment termed the senescence-associated secretory phenotype (SASP), which is beneficial or detrimental in a context-dependent manner. OIS cells are also characterized by marked epigenetic changes. We globally assessed histone modifications of OIS cells and discovered an increase in the active histone marks H3K79me2/3. The H3K79 methyltransferase disruptor of telomeric silencing 1-like (DOT1L) was necessary and sufficient for increased H3K79me2/3 occupancy at the IL1A gene locus, but not other SASP genes, and was downstream of STING. Modulating DOT1L expression did not affect the cell cycle arrest. Together, our studies establish DOT1L as an epigenetic regulator of the SASP, whose expression is uncoupled from the senescence-associated cell cycle arrest, providing a potential strategy to inhibit the negative side effects of senescence while maintaining the beneficial inhibition of proliferation.

The language of chromatin modification in human cancers

The genetic information of human cells is stored in the context of chromatin, which is subjected to DNA methylation and various histone modifications. Such a 'language' of chromatin modification constitutes a fundamental means of gene and (epi)genome regulation, underlying a myriad of cellular and developmental processes. In recent years, mounting evidence has demonstrated that miswriting, misreading or mis-erasing of the modification language embedded in chromatin represents a common, sometimes early and pivotal, event across a wide range of human cancers, contributing to oncogenesis through the induction of epigenetic, transcriptomic and phenotypic alterations. It is increasingly clear that cancer-related metabolic perturbations and oncohistone mutations also directly impact chromatin modification, thereby promoting cancerous transformation. Phase separation-based deregulation of chromatin modulators and chromatin structure is also emerging to be an important underpinning of tumorigenesis. Understanding the various molecular pathways that underscore a misregulated chromatin language in cancer, together with discovery and development of more effective drugs to target these chromatin-related vulnerabilities, will enhance treatment of human malignancies.

From the (Epi)Genome to Metabolism and Vice Versa; Examples from Hematologic Malignancy

Hematologic malignancies comprise a heterogeneous group of neoplasms arising from hematopoietic cells or their precursors and most commonly presenting as leukemias, lymphomas, and myelomas. Genetic analyses have uncovered recurrent mutations which initiate or accumulate in the course of malignant transformation, as they provide selective growth advantage to the cell. These include mutations in genes encoding transcription factors and epigenetic regulators of metabolic genes, as well as genes encoding key metabolic enzymes. The resulting alterations contribute to the extensive metabolic reprogramming characterizing the transformed cell, supporting its increased biosynthetic needs and allowing it to withstand the metabolic stress that arises as a consequence of increased metabolic rates and changes in its microenvironment. Interestingly, this cross-talk is bidirectional, as metabolites also signal back to the nucleus and, via their widespread effects on modulating epigenetic modifications, shape the chromatin landscape and the transcriptional programs of the cell. In this article, we provide an overview of the main metabolic changes and relevant genetic alterations that characterize malignant hematopoiesis and discuss how, in turn, metabolites regulate epigenetic events during this process. The aim is to illustrate the intricate interrelationship between the genome (and epigenome) and metabolism and its relevance to hematologic malignancy.

Polycomb Gene Silencing Mechanisms: PRC2 Chromatin Targeting, H3K27me3 'Readout', and Phase Separation-Based Compaction

Modulation of chromatin structure and/or modification by Polycomb repressive complexes (PRCs) provides an important means to partition the genome into functionally distinct subdomains and to regulate the activity of the underlying genes. Both the enzymatic activity of PRC2 and its chromatin recruitment, spreading, and eviction are exquisitely regulated via interactions with cofactors and DNA elements (such as unmethylated CpG islands), histones, RNA (nascent mRNA and long noncoding RNA), and R-loops. PRC2-catalyzed histone H3 lysine 27 trimethylation (H3K27me3) is recognized by distinct classes of effectors such as canonical PRC1 and BAH module-containing proteins (notably BAHCC1 in human). These effectors mediate gene silencing by different mechanisms including phase separation-related chromatin compaction and histone deacetylation. We discuss recent advances in understanding the structural architecture of PRC2, the regulation of its activity and chromatin recruitment, and the molecular mechanisms underlying Polycomb-mediated gene silencing. Because PRC deregulation is intimately associated with the development of diseases, a better appreciation of Polycomb-based (epi)genomic regulation will have far-reaching implications in biology and medicine.

Separation and Characterization of Endogenous Nucleosomes by Native Capillary Zone Electrophoresis-Top-Down Mass Spectrometry

We report a novel platform [native capillary zone electrophoresis-top-down mass spectrometry (nCZE-TDMS)] for the separation and characterization of whole nucleosomes, their histone subunits, and post-translational modifications (PTMs). As the repeating unit of chromatin, mononucleosomes (Nucs) are an ∼200 kDa complex of DNA and histone proteins involved in the regulation of key cellular processes central to human health and disease. Unraveling the covalent modification landscape of histones and their defined stoichiometries within Nucs helps to explain epigenetic regulatory mechanisms. In nCZE-TDMS, online Nuc separation is followed by a three-tier tandem MS approach that measures the intact mass of Nucs, ejects and detects the constituent histones, and fragments to sequence the histone. The new platform was optimized with synthetic Nucs to significantly reduce both sample requirements and cost compared to direct infusion. Limits of detection were in the low-attomole range, with linearity of over ∼3 orders of magnitude. The nCZE-TDMS platform was applied to endogenous Nucs from two cell lines distinguished by overexpression or knockout of histone methyltransferase NSD2/MMSET, where analysis of constituent histones revealed changes in histone abundances over the course of the CZE separation. We are confident the nCZE-TDMS platform will help advance nucleosome-level research in the fields of chromatin and epigenetics.

Capillary Zone Electrophoresis-Tandem Mass Spectrometry As an Alternative to Liquid Chromatography-Tandem Mass Spectrometry for Top-down Proteomics of Histones

Top-down proteomics (TDP) is an ideal approach for deciphering the histone code and it routinely employs reversed-phase liquid chromatography (RPLC)-tandem mass spectrometry (MS/MS). Because of the extreme complexity of histones regarding the number of proteoforms, new analytical tools with high-capacity separation and highly sensitive detection of proteoforms are required for TDP of histones. Here we present capillary zone electrophoresis (CZE)-MS/MS via the electro-kinetically pumped sheath-flow CE-MS interface for large-scale top-down delineation of histone proteoforms. CZE-MS/MS identified a comparable number of proteoforms to RPLC-MS/MS from a calf histone sample with more than 30-fold less sample consumption (75-ng vs. Three μg), indicating its substantially higher sensitivity. We identified about 400 histone proteoforms from the calf histone sample using two-dimensional size-exclusion chromatography (SEC)-CZE-MS/MS with less than 300-ng proteins consumed. We identified histone proteoforms carrying various tentative post-translational modifications (PTMs), for example, acetylation, methylation (mono-, di-, and tri-), phosphorylation, and succinylation. The electrophoretic mobility (μef) of unmodified histone proteoforms can be predicted accurately (R2 = 0.98) with an optimized semiempirical model based on our recent work. The results render CZE-MS/MS as a useful tool for deciphering the histone code in a proteoform-specific manner and on a global scale.

H3K27M in Gliomas Causes a One-Step Decrease in H3K27 Methylation and Reduced Spreading within the Constraints of H3K36 Methylation

The discovery of H3K27M mutations in pediatric gliomas marked a new chapter in cancer epigenomics. Numerous studies have investigated the effect of this mutation on H3K27 trimethylation, but only recently have we started to realize its additional effects on the epigenome. Here, we use isogenic glioma H3K27M^+/− cell lines to investigate H3K27 methylation and its interaction with H3K36 and H3K9 modifications. We describe a “step down” effect of H3K27M on the distribution of H3K27 methylation: me3 is reduced to me2, me2 is reduced to me1, whereas H3K36me2/3 delineates the boundaries for the spread of H3K27me marks. We also observe a replacement of H3K27me2/3 silencing by H3K9me3. Using a computational simulation, we explain our observations by reduced effectiveness of PRC2 and constraints imposed on the deposition of H3K27me by antagonistic H3K36 modifications. Our work further elucidates the effects of H3K27M in gliomas as well as the general principles of deposition in H3K27 methylation.

Harutyunyan et al. use isogenic glioma H3K27M^+/− cell lines to demonstrate the rewiring of the epigenome, specifically H3K27me1/2/3, H3K36me2/3, and H3K9me3. The dynamic deposition of histone marks is simulated by a stochastic model. This work further advances the understanding of the deposition of H3K27 methylation in H3K27M mutant gliomas.

Targeted detection and quantitation of histone modifications from 1,000 cells

Histone post-translational modifications (PTMs) create a powerful regulatory mechanism for maintaining chromosomal integrity in cells. Histone acetylation and methylation, the most widely studied histone PTMs, act in concert with chromatin-associated proteins to control access to genetic information during transcription. Alterations in cellular histone PTMs have been linked to disease states and have crucial biomarker and therapeutic potential. Traditional bottom-up mass spectrometry of histones requires large numbers of cells, typically one million or more. However, for some cell subtype-specific studies, it is difficult or impossible to obtain such large numbers of cells and quantification of rare histone PTMs is often unachievable. An established targeted LC-MS/MS method was used to quantify the abundance of histone PTMs from cell lines and primary human specimens. Sample preparation was modified by omitting nuclear isolation and reducing the rounds of histone derivatization to improve detection of histone peptides down to 1,000 cells. In the current study, we developed and validated a quantitative LC-MS/MS approach tailored for a targeted histone assay of 75 histone peptides with as few as 10,000 cells. Furthermore, we were able to detect and quantify 61 histone peptides from just 1,000 primary human stem cells. Detection of 37 histone peptides was possible from 1,000 acute myeloid leukemia patient cells. We anticipate that this revised method can be used in many applications where achieving large cell numbers is challenging, including rare human cell populations.

Domain Model Explains Propagation Dynamics and Stability of Histone H3K27 and H3K36 Methylation Landscapes

Chromatin states must be maintained during cell proliferation to uphold cellular identity and genome integrity. Inheritance of histone modifications is central in this process. However, the histone modification landscape is challenged by incorporation of new unmodified histones during each cell cycle, and the principles governing heritability remain unclear. We take a quantitative computational modeling approach to describe propagation of histone H3K27 and H3K36 methylation states. We measure combinatorial H3K27 and H3K36 methylation patterns by quantitative mass spectrometry on subsequent generations of histones. Using model comparison, we reject active global demethylation and invoke the existence of domains defined by distinct methylation endpoints. We find that H3K27me3 on pre-existing histones stimulates the rate of de novo H3K27me3 establishment, supporting a read-write mechanism in timely chromatin restoration. Finally, we provide a detailed quantitative picture of the mutual antagonism between H3K27 and H3K36 methylation and propose that it stabilizes epigenetic states across cell division.

NSD1-deposited H3K36me2 directs de novo methylation in the mouse male germline and counteracts Polycomb-associated silencing

De novo DNA methylation (DNAme) in mammalian germ cells is dependent on DNMT3A and DNMT3L. However, oocytes and spermatozoa show distinct patterns of DNAme. In mouse oocytes, de novo DNAme requires the lysine methyltransferase (KMTase) SETD2, which deposits H3K36me3. We show here that SETD2 is dispensable for de novo DNAme in the male germline. Instead, the lysine methyltransferase NSD1, which broadly deposits H3K36me2 in euchromatic regions, plays a critical role in de novo DNAme in prospermatogonia, including at imprinted genes. However, males deficient in germline NSD1 show a more severe defect in spermatogenesis than Dnmt3l^-/- males. Notably, unlike DNMT3L, NSD1 safeguards a subset of genes against H3K27me3-associated transcriptional silencing. In contrast, H3K36me2 in oocytes is predominantly dependent on SETD2 and coincides with H3K36me3. Furthermore, females with NSD1-deficient oocytes are fertile. Thus, the sexually dimorphic pattern of DNAme in mature mouse gametes is orchestrated by distinct profiles of H3K36 methylation.

Proteomics in Non-model Organisms: A New Analytical Frontier

For the last century we have relied on model organisms to help understand fundamental biological processes. Now, with advancements in genome sequencing, assembly, and annotation, non-model organisms may be studied with the same advanced bioanalytical toolkit as model organisms. Proteomics is one such technique, which classically relies on predicted protein sequences to catalog and measure complex proteomes across tissues and biofluids. Applying proteomics to non-model organisms can advance and accelerate biomimicry studies, biomedical advancements, veterinary medicine, agricultural research, behavioral ecology, and food safety. In this postmodel organism era, we can study almost any species, meaning that many non-model organisms are, in fact, important emerging model organisms. Herein we specifically focus on eukaryotic organisms and discuss the steps to generate sequence databases, analyze proteomic data with or without a database, and interpret results as well as future research opportunities. Proteomics is more accessible than ever before and will continue to rapidly advance in the coming years, enabling critical research and discoveries in non-model organisms that were hitherto impossible.

Benchmark problems for dynamic modeling of intracellular processes

Motivation

Dynamic models are used in systems biology to study and understand cellular processes like gene regulation or signal transduction. Frequently, ordinary differential equation (ODE) models are used to model the time and dose dependency of the abundances of molecular compounds as well as interactions and translocations. A multitude of computational approaches, e.g. for parameter estimation or uncertainty analysis have been developed within recent years. However, many of these approaches lack proper testing in application settings because a comprehensive set of benchmark problems is yet missing.

Results

We present a collection of 20 benchmark problems in order to evaluate new and existing methodologies, where an ODE model with corresponding experimental data is referred to as problem. In addition to the equations of the dynamical system, the benchmark collection provides observation functions as well as assumptions about measurement noise distributions and parameters. The presented benchmark models comprise problems of different size, complexity and numerical demands. Important characteristics of the models and methodological requirements are summarized, estimated parameters are provided, and some example studies were performed for illustrating the capabilities of the presented benchmark collection.

Availability and implementation

The models are provided in several standardized formats, including an easy-to-use human readable form and machine-readable SBML files. The data is provided as Excel sheets. All files are available at https://github.com/Benchmarking-Initiative/Benchmark-Models, including step-by-step explanations and MATLAB code to process and simulate the models.

Supplementary information

Supplementary data are available at Bioinformatics online.

Mutant EZH2 Induces a Pre-malignant Lymphoma Niche by Reprogramming the Immune Response

Follicular lymphomas (FLs) are slow-growing, indolent tumors containing extensive follicular dendritic cell (FDC) networks and recurrent EZH2 gain-of-function mutations. Paradoxically, FLs originate from highly proliferative germinal center (GC) B cells with proliferation strictly dependent on interactions with T follicular helper cells. Herein, we show that EZH2 mutations initiate FL by attenuating GC B cell requirement for T cell help and driving slow expansion of GC centrocytes that become enmeshed with and dependent on FDCs. By impairing T cell help, mutant EZH2 prevents induction of proliferative MYC programs. Thus, EZH2 mutation fosters malignant transformation by epigenetically reprograming B cells to form an aberrant immunological niche that reflects characteristic features of human FLs, explaining how indolent tumors arise from GC B cells.

Using antagonistic pleiotropy to design a chemotherapy-induced evolutionary trap to target drug resistance in cancer

Local adaptation directs populations towards environment-specific fitness maxima through acquisition of positively selected traits. However, rapid environmental changes can identify hidden fitness trade-offs that turn adaptation into maladaptation, resulting in evolutionary traps. Cancer, a disease that is prone to drug resistance, is in principle susceptible to such traps. We therefore performed pooled CRISPR-Cas9 knockout screens in acute myeloid leukemia (AML) cells treated with various chemotherapies to map the drug-dependent genetic basis of fitness trade-offs, a concept known as antagonistic pleiotropy (AP). We identified a PRC2-NSD2/3-mediated MYC regulatory axis as a drug-induced AP pathway whose ability to confer resistance to bromodomain inhibition and sensitivity to BCL-2 inhibition templates an evolutionary trap. Across diverse AML cell-line and patient-derived xenograft models, we find that acquisition of resistance to bromodomain inhibition through this pathway exposes coincident hypersensitivity to BCL-2 inhibition. Thus, drug-induced AP can be leveraged to design evolutionary traps that selectively target drug resistance in cancer.

Comparative Proteomic Analysis of Histone Modifications upon Acridone Derivative -Induced CCRF-CEM Cells by Data Independent Acquisition

The lead compound acridone derivative 8a showed potent antiproliferative activity by inducing DNA damage through direct stacking with DNA bases and triggering ROS in CCRF-CEM cells. To define the chromatin alterations during DNA damage sensing and repair, a detailed quantitative map of single and coexisting histone post-translational modifications (PTMs) in CCRF-CEM cells affected by 8a was performed by the Data Independent Acquisition (DIA) method on QE-plus. A total of 79 distinct and 164 coexisting histone PTMs were quantified, of which 16 distinct histone PTMs were significantly altered when comparing 8a-treated cells with vehicle control cells. The changes in histone PTMs were confirmed by Western blotting analysis for three H3 and one H4 histone markers. The up-regulated dimethylation on H3K9, H3K36, and H4K20 implied that CCRF-CEM cells might accelerate DNA damage repair to counteract the DNA lesion induced by 8a, which was verified by an increment in the 53BP1 foci localization at the damaged DNA. Most of the significantly altered PTMs were involved in transcriptional regulation, including down-regulated acetylation on H3K18, H3K27, and H3K122, and up-regulated di- and trimethylation on H3K9 and H3K27. This transcription-silencing phenomenon was associated with G2/M cell cycle arrest after 8a treatment by flow cytometry. This study shows that the DIA proteomics strategy provides a sensitive and accurate way to characterize the coexisting histone PTMs changes and their cross-talk in CCRF-CEM cells after 8a treatment. Specifically, histone PTMs rearrange transcription-silencing, and cell cycle arrest DNA damage repair may contribute to the mechanism of epigenetic response affected by 8a.

Coupling Fluorescence-Activated Cell Sorting and Targeted Analysis of Histone Modification Profiles in Primary Human Leukocytes

Histone posttranslational modifications (PTMs) are essential for regulating chromatin and maintaining gene expression throughout cell differentiation. Despite the deep level of understanding of immunophenotypic differentiation pathways in hematopoietic cells, few studies have investigated global levels of histone PTMs required for differentiation and maintenance of these distinct cell types. Here, we describe an approach to couple fluorescence-activated cell sorting (FACS) with targeted mass spectrometry to define global "epi-proteomic" signatures for primary leukocytes. FACS was used to sort closely and distantly related leukocytes from normal human peripheral blood for quantitation of histone PTMs with a multiple reaction monitoring LC-MS/MS method measuring histone PTMs on histones H3 and H4. We validate cell sorting directly into H₂SO₄ for immediate histone extraction to decrease time and number of steps after FACS to analyze histone PTMs. Relative histone PTM levels vary in T cells across healthy donors, and the majority of PTMs remain stable up to 2 days following initial blood draw. Large differences in the levels of histone PTMs are observed across the mature lymphoid and myeloid lineages, as well as between different types within the same lineage, though no differences are observed in closely related T cell subtypes. The results show a streamlined approach for quantifying global changes in histone PTMs in cell types separated by FACS that is poised for clinical deployment.

Mass spectrometric analysis of PTM dynamics using stable isotope labeled metabolic precursors in cell culture

Biological organisms represent highly dynamic systems, which are continually exposed to environmental factors and always strive to restore steady-state homeostasis. Posttranslational modifications are key regulators with which biological systems respond to external stimuli. To understand how homeostasis is restored, it is important to study the kinetics of posttranslational modifications. In this review we discuss proteomic approaches using stable isotope labeled metabolic precursors to study dynamics of posttranslational modifications in cell culture.

PHF19 promotes multiple myeloma tumorigenicity through PRC2 activation and broad H3K27me3 domain formation

Dysregulation of polycomb repressive complex 2 (PRC2) promotes oncogenesis partly through its enzymatic function for inducing trimethylation of histone H3 lysine 27 (H3K27me3). However, it remains to be determined how PRC2 activity is regulated in normal and diseased settings. We here report a PRC2-associated cofactor, PHD finger protein 19 (PHF19; also known as polycomb-like 3), as a crucial mediator of tumorigenicity in multiple myeloma (MM). Overexpression and/or genomic amplification of PHF19 is found associated with malignant progression of MM and plasma cell leukemia, correlating to worse treatment outcomes. Using various MM models, we demonstrated a critical requirement of PHF19 for tumor growth in vitro and in vivo. Mechanistically, PHF19-mediated oncogenic effect relies on its PRC2-interacting and chromatin-binding functions. Chromatin immunoprecipitation followed by sequencing profiling showed a critical role for PHF19 in maintaining the H3K27me3 landscape. PHF19 depletion led to loss of broad H3K27me3 domains, possibly due to impaired H3K27me3 spreading from cytosine guanine dinucleotide islands, which is reminiscent to the reported effect of an "onco"-histone mutation, H3K27 to methionine (H3K27M). RNA-sequencing-based transcriptome profiling in MM lines also demonstrated a requirement of PHF19 for optimal silencing of PRC2 targets, which include cell cycle inhibitors and interferon-JAK-STAT signaling genes critically involved in tumor suppression. Correlation studies using patient sample data sets further support a clinical relevance of the PHF19-regulated pathways. Lastly, we show that MM cells are generally sensitive to PRC2 inhibitors. Collectively, this study demonstrates that PHF19 promotes MM tumorigenesis through enhancing H3K27me3 deposition and PRC2's gene-regulatory functions, lending support for PRC2 blockade as a means for MM therapeutics.

PRC2 functions in development and congenital disorders

Polycomb repressive complex 2 (PRC2) is a conserved chromatin regulator that is responsible for the methylation of histone H3 lysine 27 (H3K27). PRC2 is essential for normal development and its loss of function thus results in a range of developmental phenotypes. Here, we review the latest advances in our understanding of mammalian PRC2 activity and present an updated summary of the phenotypes associated with its loss of function in mice. We then discuss recent studies that have highlighted regulatory interplay between the modifications laid down by PRC2 and other chromatin modifiers, including NSD1 and DNMT3A. Finally, we propose a model in which the dysregulation of these modifications at intergenic regions is a shared molecular feature of genetically distinct but highly phenotypically similar overgrowth syndromes in humans.

PyBioNetFit and the Biological Property Specification Language

In systems biology modeling, important steps include model parameterization, uncertainty quantification, and evaluation of agreement with experimental observations. To help modelers perform these steps, we developed the software PyBioNetFit, which in addition supports checking models against known system properties and solving design problems. PyBioNetFit introduces Biological Property Specification Language (BPSL) for the formal declaration of system properties. BPSL allows qualitative data to be used alone or in combination with quantitative data. PyBioNetFit performs parameterization with parallelized metaheuristic optimization algorithms that work directly with existing model definition standards: BioNetGen Language (BNGL) and Systems Biology Markup Language (SBML). We demonstrate PyBioNetFit's capabilities by solving various example problems, including the challenging problem of parameterizing a 153-parameter model of cell cycle control in yeast based on both quantitative and qualitative data. We demonstrate the model checking and design applications of PyBioNetFit and BPSL by analyzing a model of targeted drug interventions in autophagy signaling.

A clinically relevant murine model unmasks a "two-hit" mechanism for reactivation and dissemination of cytomegalovirus after kidney transplant

Reactivation of latent cytomegalovirus remains an important complication after transplant. Although immunosuppression (IS) has been implicated as a primary cause, we have previously shown that the implantation response of a kidney allograft can lead to early transcriptional activation of latent murine cytomegalovirus (MCMV) genes in an immune-competent host and to MCMV reactivation and dissemination to other organs in a genetically immune-deficient recipient. We now describe a model that allows us to separately analyze the impact of the implantation effect vs that of a clinically relevant IS regimen. Treatment with IS of latently infected mice alone does not induce viral reactivation, but transplant of latently infected allogeneic kidneys combined with IS facilitates MCMV reactivation in the graft and dissemination to other organs. The IS regimen effectively dampens allo-immune inflammatory pathways and depletes recipient anti-MCMV but does not affect ischemia-reperfusion injury pathways. MCMV reactivation similar to that seen in allogeneic transplants combined with also occurs after syngeneic transplants. Thus, our data strongly suggest that while ischemia-reperfusion injury of the implanted graft is sufficient and necessary to initiate transcriptional reactivation of latent MCMV ("first hit"), IS is permissive to the first hit and facilitates dissemination to other organs ("second hit").

PRC2 is high maintenance

As the process that silences gene expression ensues during development, the stage is set for the activity of Polycomb-repressive complex 2 (PRC2) to maintain these repressed gene profiles. PRC2 catalyzes a specific histone posttranslational modification (hPTM) that fosters chromatin compaction. PRC2 also facilitates the inheritance of this hPTM through its self-contained "write and read" activities, key to preserving cellular identity during cell division. As these changes in gene expression occur without changes in DNA sequence and are inherited, the process is epigenetic in scope. Mutants of mammalian PRC2 or of its histone substrate contribute to the cancer process and other diseases, and research into these aberrant pathways is yielding viable candidates for therapeutic targeting. The effectiveness of PRC2 hinges on its being recruited to the proper chromatin sites; however, resolving the determinants to this process in the mammalian case was not straightforward and thus piqued the interest of many in the field. Here, we chronicle the latest advances toward exposing mammalian PRC2 and its high maintenance.

Understanding histone H3 lysine 36 methylation and its deregulation in disease

Methylation of histone H3 lysine 36 (H3K36) plays crucial roles in the partitioning of chromatin to distinctive domains and the regulation of a wide range of biological processes. Trimethylation of H3K36 (H3K36me3) demarcates body regions of the actively transcribed genes, providing signals for modulating transcription fidelity, mRNA splicing and DNA damage repair; and di-methylation of H3K36 (H3K36me2) spreads out within large intragenic regions, regulating distribution of histone H3 lysine 27 trimethylation (H3K27me3) and possibly DNA methylation. These H3K36 methylation-mediated events are biologically crucial and controlled by different classes of proteins responsible for either 'writing', 'reading' or 'erasing' of H3K36 methylation marks. Deregulation of H3K36 methylation and related regulatory factors leads to pathogenesis of disease such as developmental syndrome and cancer. Additionally, recurrent mutations of H3K36 and surrounding histone residues are detected in human tumors, further highlighting the importance of H3K36 in biology and medicine. This review will elaborate on current advances in understanding H3K36 methylation and related molecular players during various chromatin-templated cellular processes, their crosstalks with other chromatin factors, as well as their deregulations in the diseased contexts.