A Multiplexed System for Quantitative Comparisons of Chromatin Landscapes

New frameworks for hematopoiesis derived from single-cell genomics

ABSTRACT

Recent advancements in single-cell genomics have enriched our understanding of hematopoiesis, providing intricate details about hematopoietic stem cell biology, differentiation, and lineage commitment. Technological advancements have highlighted extensive heterogeneity of cell populations and continuity of differentiation routes. Nevertheless, intermediate "attractor" states signify structure in stem and progenitor populations that link state transition dynamics to fate potential. We discuss how innovative model systems quantify lineage bias and how stress accelerates differentiation, thereby reducing fate plasticity compared with native hematopoiesis. We conclude by offering our perspective on the current model of hematopoiesis and discuss how a more precise understanding can translate to strategies that extend healthy hematopoiesis and prevent disease.

Integrated epigenomic exposure signature discovery

Aim: The epigenome influences gene regulation and phenotypes in response to exposures. Epigenome assessment can determine exposure history aiding in diagnosis.Materials & methods: Here we developed and implemented a machine learning algorithm, the exposure signature discovery algorithm (ESDA), to identify the most important features present in multiple epigenomic and transcriptomic datasets to produce an integrated exposure signature (ES).Results: Signatures were developed for seven exposures including Staphylococcus aureus, human immunodeficiency virus, SARS-CoV-2, influenza A (H3N2) virus and Bacillus anthracis vaccinations. ESs differed in the assays and features selected and predictive value.Conclusion: Integrated ESs can potentially be utilized for diagnosis or forensic attribution. The ESDA identifies the most distinguishing features enabling diagnostic panel development for future precision health deployment.

An autoimmune transcriptional circuit drives FOXP3+ regulatory T cell dysfunction

Autoimmune diseases, among the most common disorders of young adults, are mediated by genetic and environmental factors. Although CD4+FOXP3+ regulatory T cells (Tregs) play a central role in preventing autoimmunity, the molecular mechanism underlying their dysfunction is unknown. Here, we performed comprehensive transcriptomic and epigenomic profiling of Tregs in the autoimmune disease multiple sclerosis (MS) to identify critical transcriptional programs regulating human autoimmunity. We found that up-regulation of a primate-specific short isoform of PR domain zinc finger protein 1 (PRDM1-S) induces expression of serum and glucocorticoid-regulated kinase 1 (SGK1) independent from the evolutionarily conserved long PRDM1, which led to destabilization of forkhead box P3 (FOXP3) and Treg dysfunction. This aberrant PRDM1-S/SGK1 axis is shared among other autoimmune diseases. Furthermore, the chromatin landscape profiling in Tregs from individuals with MS revealed enriched activating protein-1 (AP-1)/interferon regulatory factor (IRF) transcription factor binding as candidate upstream regulators of PRDM1-S expression and Treg dysfunction. Our study uncovers a mechanistic model where the evolutionary emergence of PRDM1-S and epigenetic priming of AP-1/IRF may be key drivers of dysfunctional Tregs in autoimmune diseases.

Senescence of human pancreatic beta cells enhances functional maturation through chromatin reorganization and promotes interferon responsiveness

Senescent cells can influence the function of tissues in which they reside, and their propensity for disease. A portion of adult human pancreatic beta cells express the senescence marker p16, yet it is unclear whether they are in a senescent state, and how this affects insulin secretion. We analyzed single-cell transcriptome datasets of adult human beta cells, and found that p16-positive cells express senescence gene signatures, as well as elevated levels of beta-cell maturation genes, consistent with enhanced functionality. Senescent human beta-like cells in culture undergo chromatin reorganization that leads to activation of enhancers regulating functional maturation genes and acquisition of glucose-stimulated insulin secretion capacity. Strikingly, Interferon-stimulated genes are elevated in senescent human beta cells, but genes encoding senescence-associated secretory phenotype (SASP) cytokines are not. Senescent beta cells in culture and in human tissue show elevated levels of cytoplasmic DNA, contributing to their increased interferon responsiveness. Human beta-cell senescence thus involves chromatin-driven upregulation of a functional-maturation program, and increased responsiveness of interferon-stimulated genes, changes that could increase both insulin secretion and immune reactivity.

Loss of CREBBP and KMT2D cooperate to accelerate lymphomagenesis and shape the lymphoma immune microenvironment

Despite regulating overlapping gene enhancers and pathways, CREBBP and KMT2D mutations recurrently co-occur in germinal center (GC) B cell-derived lymphomas, suggesting potential oncogenic cooperation. Herein, we report that combined haploinsufficiency of Crebbp and Kmt2d induces a more severe mouse lymphoma phenotype (vs either allele alone) and unexpectedly confers an immune evasive microenvironment manifesting as CD8+ T-cell exhaustion and reduced infiltration. This is linked to profound repression of immune synapse genes that mediate crosstalk with T-cells, resulting in aberrant GC B cell fate decisions. From the epigenetic perspective, we observe interaction and mutually dependent binding and function of CREBBP and KMT2D on chromatin. Their combined deficiency preferentially impairs activation of immune synapse-responsive super-enhancers, pointing to a particular dependency for both co-activators at these specialized regulatory elements. Together, our data provide an example where chromatin modifier mutations cooperatively shape and induce an immune-evasive microenvironment to facilitate lymphomagenesis.

Disease-associated astrocyte epigenetic memory promotes CNS pathology

Disease-associated astrocyte subsets contribute to the pathology of neurologic diseases, including multiple sclerosis and experimental autoimmune encephalomyelitis1-8 (EAE), an experimental model for multiple sclerosis. However, little is known about the stability of these astrocyte subsets and their ability to integrate past stimulation events. Here we report the identification of an epigenetically controlled memory astrocyte subset that exhibits exacerbated pro-inflammatory responses upon rechallenge. Specifically, using a combination of single-cell RNA sequencing, assay for transposase-accessible chromatin with sequencing, chromatin immunoprecipitation with sequencing, focused interrogation of cells by nucleic acid detection and sequencing, and cell-specific in vivo CRISPR-Cas9-based genetic perturbation studies we established that astrocyte memory is controlled by the metabolic enzyme ATP-citrate lyase (ACLY), which produces acetyl coenzyme A (acetyl-CoA) that is used by histone acetyltransferase p300 to control chromatin accessibility. The number of ACLY+p300+ memory astrocytes is increased in acute and chronic EAE models, and their genetic inactivation ameliorated EAE. We also detected the pro-inflammatory memory phenotype in human astrocytes in vitro; single-cell RNA sequencing and immunohistochemistry studies detected increased numbers of ACLY+p300+ astrocytes in chronic multiple sclerosis lesions. In summary, these studies define an epigenetically controlled memory astrocyte subset that promotes CNS pathology in EAE and, potentially, multiple sclerosis. These findings may guide novel therapeutic approaches for multiple sclerosis and other neurologic diseases.

Epigenetic modulation of myeloid cell functions in HIV and SARS-CoV-2 infection

Myeloid cells play a vital role in innate immune responses as they recognize and phagocytose pathogens like viruses, present antigens, produce cytokines, recruit other immune cells to combat infections, and contribute to the attenuation of immune responses to restore homeostasis. Signal integration by pathogen recognition receptors enables myeloid cells to adapt their functions by a network of transcription factors and chromatin remodelers. This review provides a brief overview of the subtypes of myeloid cells and the main epigenetic regulation mechanisms. Special focus is placed on the epigenomic alterations in viral nucleic acids of HIV and SARS-CoV-2 along with the epigenetic changes in the host's myeloid cell compartment. These changes are important as they lead to immune suppression and promote the progression of the disease. Finally, we highlight some promising examples of 'epidrugs' that modulate the epigenome of immune cells and could be used as therapeutics for viral infections.

Disease-associated astrocyte epigenetic memory promotes CNS pathology

Astrocytes play important roles in the central nervous system (CNS) physiology and pathology. Indeed, astrocyte subsets defined by specific transcriptional activation states contribute to the pathology of neurologic diseases, including multiple sclerosis (MS) and its pre-clinical model experimental autoimmune encephalomyelitis (EAE) 1-8 . However, little is known about the stability of these disease-associated astrocyte subsets, their regulation, and whether they integrate past stimulation events to respond to subsequent challenges. Here, we describe the identification of an epigenetically controlled memory astrocyte subset which exhibits exacerbated pro-inflammatory responses upon re-challenge. Specifically, using a combination of single-cell RNA sequencing (scRNA-seq), assay for transposase-accessible chromatin with sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing (ChIP-seq), focused interrogation of cells by nucleic acid detection and sequencing (FIND-seq), and cell-specific in vivo CRISPR/Cas9-based genetic perturbation studies we established that astrocyte memory is controlled by the metabolic enzyme ATP citrate lyase (ACLY), which produces acetyl coenzyme A (acetyl-CoA) used by the histone acetyltransferase p300 to control chromatin accessibility. ACLY + p300 + memory astrocytes are increased in acute and chronic EAE models; the genetic targeting of ACLY + p300 + astrocytes using CRISPR/Cas9 ameliorated EAE. We also detected responses consistent with a pro-inflammatory memory phenotype in human astrocytes in vitro ; scRNA-seq and immunohistochemistry studies detected increased ACLY + p300 + astrocytes in chronic MS lesions. In summary, these studies define an epigenetically controlled memory astrocyte subset that promotes CNS pathology in EAE and, potentially, MS. These findings may guide novel therapeutic approaches for MS and other neurologic diseases.

Omics Analysis and Quality Control Pipelines in a High-Performance Computing Environment

Data quality is often an overlooked feature in the analysis of omics data. This is particularly relevant in studies of chemical and pathogen exposures that can modify an individual's epigenome and transcriptome with persistence over time. Portable, quality control (QC) pipelines for multiple different omics datasets are therefore needed. To meet these goals, portable quality assurance (QA) metrics, metric acceptability criterion, and pipelines to compute these metrics were developed and consolidated into one framework for 12 different omics assays. Performance of these QA metrics and pipelines were evaluated on human data generated by the Defense Advanced Research Projects Agency (DARPA) Epigenetic CHaracterization and Observation (ECHO) program. Twelve analytical pipelines were developed leveraging standard tools when possible. These QC pipelines were containerized using Singularity to ensure portability and scalability. Datasets for these 12 omics assays were analyzed and results were summarized. The quality thresholds and metrics used were described. We found that these pipelines enabled early identification of lower quality datasets, datasets with insufficient reads for additional sequencing, and experimental protocols needing refinements. These omics data analysis and QC pipelines are available as open-source resources as reported and discussed in this article for the omics and life sciences communities.

Dihydrofolate reductase activity controls neurogenic transitions in the developing neocortex

One-carbon/folate (1C) metabolism supplies methyl groups required for DNA and histone methylation, and is involved in the maintenance of self-renewal in stem cells. Dihydrofolate reductase (DHFR), a key enzyme in 1C metabolism, is highly expressed in human and mouse neural progenitors at the early stages of neocortical development. Here, we have investigated the role of DHFR in the developing neocortex and report that reducing its activity in human neural organoids and mouse embryonic neocortex accelerates indirect neurogenesis, thereby affecting neuronal composition of the neocortex. Furthermore, we show that decreasing DHFR activity in neural progenitors leads to a reduction in one-carbon/folate metabolites and correlates with modifications of H3K4me3 levels. Our findings reveal an unanticipated role for DHFR in controlling specific steps of neocortex development and indicate that variations in 1C metabolic cues impact cell fate transitions.

The pioneer factor SOX9 competes for epigenetic factors to switch stem cell fates

During development, progenitors simultaneously activate one lineage while silencing another, a feature highly regulated in adult stem cells but derailed in cancers. Equipped to bind cognate motifs in closed chromatin, pioneer factors operate at these crossroads, but how they perform fate switching remains elusive. Here we tackle this question with SOX9, a master regulator that diverts embryonic epidermal stem cells (EpdSCs) into becoming hair follicle stem cells. By engineering mice to re-activate SOX9 in adult EpdSCs, we trigger fate switching. Combining epigenetic, proteomic and functional analyses, we interrogate the ensuing chromatin and transcriptional dynamics, slowed temporally by the mature EpdSC niche microenvironment. We show that as SOX9 binds and opens key hair follicle enhancers de novo in EpdSCs, it simultaneously recruits co-factors away from epidermal enhancers, which are silenced. Unhinged from its normal regulation, sustained SOX9 subsequently activates oncogenic transcriptional regulators that chart the path to cancers typified by constitutive SOX9 expression.

The ENCODE Uniform Analysis Pipelines

The Encyclopedia of DNA elements (ENCODE) project is a collaborative effort to create a comprehensive catalog of functional elements in the human genome. The current database comprises more than 19000 functional genomics experiments across more than 1000 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of the Homo sapiens and Mus musculus genomes. All experimental data, metadata, and associated computational analyses created by the ENCODE consortium are submitted to the Data Coordination Center (DCC) for validation, tracking, storage, and distribution to community resources and the scientific community. The ENCODE project has engineered and distributed uniform processing pipelines in order to promote data provenance and reproducibility as well as allow interoperability between genomic resources and other consortia. All data files, reference genome versions, software versions, and parameters used by the pipelines are captured and available via the ENCODE Portal. The pipeline code, developed using Docker and Workflow Description Language (WDL; https://openwdl.org/) is publicly available in GitHub, with images available on Dockerhub (https://hub.docker.com), enabling access to a diverse range of biomedical researchers. ENCODE pipelines maintained and used by the DCC can be installed to run on personal computers, local HPC clusters, or in cloud computing environments via Cromwell. Access to the pipelines and data via the cloud allows small labs the ability to use the data or software without access to institutional compute clusters. Standardization of the computational methodologies for analysis and quality control leads to comparable results from different ENCODE collections - a prerequisite for successful integrative analyses.

Unveiling the complexity of transcription factor networks in hematopoietic stem cells: implications for cell therapy and hematological malignancies

The functionality and longevity of hematopoietic tissue is ensured by a tightly controlled balance between self-renewal, quiescence, and differentiation of hematopoietic stem cells (HSCs) into the many different blood lineages. Cell fate determination in HSCs is influenced by signals from extrinsic factors (e.g., cytokines, irradiation, reactive oxygen species, O2 concentration) that are translated and integrated by intrinsic factors such as Transcription Factors (TFs) to establish specific gene regulatory programs. TFs also play a central role in the establishment and/or maintenance of hematological malignancies, highlighting the need to understand their functions in multiple contexts. TFs bind to specific DNA sequences and interact with each other to form transcriptional complexes that directly or indirectly control the expression of multiple genes. Over the past decades, significant research efforts have unraveled molecular programs that control HSC function. This, in turn, led to the identification of more than 50 TF proteins that influence HSC fate. However, much remains to be learned about how these proteins interact to form molecular networks in combination with cofactors (e.g. epigenetics factors) and how they control differentiation, expansion, and maintenance of cellular identity. Understanding these processes is critical for future applications particularly in the field of cell therapy, as this would allow for manipulation of cell fate and induction of expansion, differentiation, or reprogramming of HSCs using specific cocktails of TFs. Here, we review recent findings that have unraveled the complexity of molecular networks controlled by TFs in HSCs and point towards possible applications to obtain functional HSCs ex vivo for therapeutic purposes including hematological malignancies. Furthermore, we discuss the challenges and prospects for the derivation and expansion of functional adult HSCs in the near future.

Tracking Targets of Dynamic Super-Enhancers in Vitro to Better Characterize Osteoclastogenesis and to Evaluate the Effect of Diuron on the Maturation of Human Bone Cells

BACKGROUND

Osteoclasts are major actors in the maintenance of bone homeostasis. The full functional maturation of osteoclasts from monocyte lineage cells is essential for the degradation of old/damaged bone matrix. Diuron is one of the most frequently encountered herbicides, particularly in water sources. However, despite a reported delayed ossification in vivo, its impact on bone cells remains largely unknown.

OBJECTIVES

The objectives of this study were to first better characterize osteoclastogenesis by identifying genes that drive the differentiation of CD14+ monocyte progenitors into osteoclasts and to evaluate the toxicity of diuron on osteoblastic and osteoclastic differentiation in vitro.

METHODS

We performed chromatin immunoprecipitation (ChIP) against H3K27ac followed by ChIP-sequencing (ChIP-Seq) and RNA-sequencing (RNA-Seq) at different stages of differentiation of CD14+ monocytes into active osteoclasts. Differentially activated super-enhancers and their potential target genes were identified. Then to evaluate the toxicity of diuron on osteoblasts and osteoclasts, we performed RNA-Seq and functional tests during in vitro osteoblastic and osteoclastic differentiation by exposing cells to different concentrations of diuron.

RESULTS

The combinatorial study of the epigenetic and transcriptional remodeling taking place during differentiation has revealed a very dynamic epigenetic profile that supports the expression of genes vital for osteoclast differentiation and function. In total, we identified 122 genes induced by dynamic super-enhancers at late days. Our data suggest that high concentration of diuron (50μM) affects viability of mesenchymal stem cells (MSCs) in vitro associated with a decrease of bone mineralization. At a lower concentration (1μM), an inhibitory effect was observed in vitro on the number of osteoclasts derived from CD14+ monocytes without affecting cell viability. Among the diuron-affected genes, our analysis suggests a significant enrichment of genes targeted by pro-differentiation super-enhancers, with an odds ratio of 5.12 (ρ=2.59×10-5).

DISCUSSION

Exposure to high concentrations of diuron decreased the viability of MSCs and could therefore affect osteoblastic differentiation and bone mineralization. This pesticide also disrupted osteoclasts maturation by impairing the expression of cell-identity determining genes. Indeed, at sublethal concentrations, differences in the expression of these key genes were mild during the course of in vitro osteoclast differentiation. Taken together our results suggest that high exposure levels of diuron could have an effect on bone homeostasis. https://doi.org/10.1289/EHP11690.

High-throughput methods for the analysis of transcription factors and chromatin modifications: Low input, single cell and spatial genomic technologies

Genome-wide analysis of transcription factors and epigenomic features is instrumental to shed light on DNA-templated regulatory processes such as transcription, cellular differentiation or to monitor cellular responses to environmental cues. Two decades of technological developments have led to a rich set of approaches progressively pushing the limits of epigenetic profiling towards single cells. More recently, disruptive technologies using innovative biochemistry came into play. Assays such as CUT&RUN, CUT&Tag and variations thereof show considerable potential to survey multiple TFs or histone modifications in parallel from a single experiment and in native conditions. These are in the path to become the dominant assays for genome-wide analysis of TFs and chromatin modifications in bulk, single-cell, and spatial genomic applications. The principles together with pros and cons are discussed.

The ENCODE Uniform Analysis Pipelines

The Encyclopedia of DNA elements (ENCODE) project is a collaborative effort to create a comprehensive catalog of functional elements in the human genome. The current database comprises more than 19000 functional genomics experiments across more than 1000 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of the Homo sapiens and Mus musculus genomes. All experimental data, metadata, and associated computational analyses created by the ENCODE consortium are submitted to the Data Coordination Center (DCC) for validation, tracking, storage, and distribution to community resources and the scientific community. The ENCODE project has engineered and distributed uniform processing pipelines in order to promote data provenance and reproducibility as well as allow interoperability between genomic resources and other consortia. All data files, reference genome versions, software versions, and parameters used by the pipelines are captured and available via the ENCODE Portal. The pipeline code, developed using Docker and Workflow Description Language (WDL; https://openwdl.org/) is publicly available in GitHub, with images available on Dockerhub (https://hub.docker.com), enabling access to a diverse range of biomedical researchers. ENCODE pipelines maintained and used by the DCC can be installed to run on personal computers, local HPC clusters, or in cloud computing environments via Cromwell. Access to the pipelines and data via the cloud allows small labs the ability to use the data or software without access to institutional compute clusters. Standardization of the computational methodologies for analysis and quality control leads to comparable results from different ENCODE collections - a prerequisite for successful integrative analyses.

Regulation and Immunotherapeutic Targeting of the Epigenome in Exhausted CD8 T Cell Responses

Exhaustion is a state of CD8 T cell differentiation that occurs in settings of chronic Ag such as tumors, chronic viral infection, and autoimmunity. Cellular differentiation is driven by a series of environmental signals that promote epigenetic landscapes that set transcriptomes needed for function. For CD8 T cells, the epigenome that underlies exhaustion is distinct from effector and memory cell differentiation, suggesting that signals early on set in motion a process where the epigenome is modified to promote a trajectory toward a dysfunctional state. Although we know many signals that promote exhaustion, putting this in the context of the epigenetic changes that occur during differentiation has been less clear. In this review, we aim to summarize the epigenetic changes associated with exhaustion in the context of signals that promote it, highlighting immunotherapeutic studies that support these observations or areas for future therapeutic opportunities.

Dynamic antagonism between key repressive pathways maintains the placental epigenome

DNA and Histone 3 Lysine 27 methylation typically function as repressive modifications and operate within distinct genomic compartments. In mammals, the majority of the genome is kept in a DNA methylated state, whereas the Polycomb repressive complexes regulate the unmethylated CpG-rich promoters of developmental genes. In contrast to this general framework, the extra-embryonic lineages display non-canonical, globally intermediate DNA methylation levels, including disruption of local Polycomb domains. Here, to better understand this unusual landscape's molecular properties, we genetically and chemically perturbed major epigenetic pathways in mouse trophoblast stem cells. We find that the extra-embryonic epigenome reflects ongoing and dynamic de novo methyltransferase recruitment, which is continuously antagonized by Polycomb to maintain intermediate, locally disordered methylation. Despite its disorganized molecular appearance, our data point to a highly controlled equilibrium between counteracting repressors within extra-embryonic cells, one that can seemingly persist indefinitely without bistable features typically seen for embryonic forms of epigenetic regulation.

PRC2-independent actions of H3.3K27M in embryonic stem cell differentiation

The histone H3 variant, H3.3, is localized at specific regions in the genome, especially promoters and active enhancers, and has been shown to play important roles in development. A lysine to methionine substitution in position 27 (H3.3K27M) is a main cause of Diffuse Intrinsic Pontine Glioma (specifically Diffuse Midline Glioma, K27M-mutant), a lethal type of pediatric cancer. H3.3K27M has a dominant-negative effect by inhibiting the Polycomb Repressor Complex 2 (PRC2) activity. Here, we studied the immediate, genome-wide, consequences of the H3.3K27M mutation independent of PRC2 activity. We developed Doxycycline (Dox)-inducible mouse embryonic stem cells (ESCs) carrying a single extra copy of WT-H3.3, H3.3K27M and H3.3K27L, all fused to HA. We performed RNA-Seq and ChIP-Seq at different times following Dox induction in undifferentiated and differentiated ESCs. We find increased binding of H3.3 around transcription start sites in cells expressing both H3.3K27M and H3.3K27L compared with WT, but not in cells treated with PRC2 inhibitors. Differentiated cells carrying either H3.3K27M or H3.3K27L retain expression of ESC-active genes, in expense of expression of genes related to neuronal differentiation. Taken together, our data suggest that a modifiable H3.3K27 is required for proper histone incorporation and cellular maturation, independent of PRC2 activity.

Cooperative super-enhancer inactivation caused by heterozygous loss of CREBBP and KMT2D skews B cell fate decisions and yields T cell-depleted lymphomas

Mutations affecting enhancer chromatin regulators CREBBP and KMT2D are highly co-occurrent in germinal center (GC)-derived lymphomas and other tumors, even though regulating similar pathways. Herein, we report that combined haploinsufficiency of Crebbp and Kmt2d (C+K) indeed accelerated lymphomagenesis. C+K haploinsufficiency induced GC hyperplasia by altering cell fate decisions, skewing B cells away from memory and plasma cell differentiation. C+K deficiency particularly impaired enhancer activation for immune synapse genes involved in exiting the GC reaction. This effect was especially severe at super-enhancers for immunoregulatory and differentiation genes. Mechanistically, CREBBP and KMT2D formed a complex, were highly co-localized on chromatin, and were required for each-other's stable recruitment to enhancers. Notably, C+K lymphomas in mice and humans manifested significantly reduced CD8 + T-cell abundance. Hence, deficiency of C+K cooperatively induced an immune evasive phenotype due at least in part to failure to activate key immune synapse super-enhancers, associated with altered immune cell fate decisions.

SIGNIFICANCE

Although CREBBP and KMT2D have similar enhancer regulatory functions, they are paradoxically co-mutated in lymphomas. We show that their combined loss causes specific disruption of super-enhancers driving immune synapse genes. Importantly, this leads to reduction of CD8 cells in lymphomas, linking super-enhancer function to immune surveillance, with implications for immunotherapy resistance.

Viewing AML through a New Lens: Technological Advances in the Study of Epigenetic Regulation

Epigenetic modifications, such as histone modifications and DNA methylation, are essential for ensuring the dynamic control of gene regulation in every cell type. These modifications are associated with gene activation or repression, depending on the genomic context and specific type of modification. In both cases, they are deposited and removed by epigenetic modifier proteins. In acute myeloid leukemia (AML), the function of these proteins is perturbed through genetic mutations (i.e., in the DNA methylation machinery) or translocations (i.e., MLL-rearrangements) arising during leukemogenesis. This can lead to an imbalance in the epigenomic landscape, which drives aberrant gene expression patterns. New technological advances, such as CRISPR editing, are now being used to precisely model genetic mutations and chromosomal translocations. In addition, high-precision epigenomic editing using dCas9 or CRISPR base editing are being used to investigate the function of epigenetic mechanisms in gene regulation. To interrogate these mechanisms at higher resolution, advances in single-cell techniques have begun to highlight the heterogeneity of epigenomic landscapes and how these impact on gene expression within different AML populations in individual cells. Combined, these technologies provide a new lens through which to study the role of epigenetic modifications in normal hematopoiesis and how the underlying mechanisms can be hijacked in the context of malignancies such as AML.

Microglia states and nomenclature: A field at its crossroads

Microglial research has advanced considerably in recent decades yet has been constrained by a rolling series of dichotomies such as "resting versus activated" and "M1 versus M2." This dualistic classification of good or bad microglia is inconsistent with the wide repertoire of microglial states and functions in development, plasticity, aging, and diseases that were elucidated in recent years. New designations continuously arising in an attempt to describe the different microglial states, notably defined using transcriptomics and proteomics, may easily lead to a misleading, although unintentional, coupling of categories and functions. To address these issues, we assembled a group of multidisciplinary experts to discuss our current understanding of microglial states as a dynamic concept and the importance of addressing microglial function. Here, we provide a conceptual framework and recommendations on the use of microglial nomenclature for researchers, reviewers, and editors, which will serve as the foundations for a future white paper.

Epigenetic reader SP140 loss of function drives Crohn's disease due to uncontrolled macrophage topoisomerases

How mis-regulated chromatin directly impacts human immune disorders is poorly understood. Speckled Protein 140 (SP140) is an immune-restricted PHD and bromodomain-containing epigenetic "reader," and SP140 loss-of-function mutations associate with Crohn's disease (CD), multiple sclerosis (MS), and chronic lymphocytic leukemia (CLL). However, the relevance of these mutations and mechanisms underlying SP140-driven pathogenicity remains unexplored. Using a global proteomic strategy, we identified SP140 as a repressor of topoisomerases (TOPs) that maintains heterochromatin and macrophage fate. In humans and mice, SP140 loss resulted in unleashed TOP activity, de-repression of developmentally silenced genes, and ultimately defective microbe-inducible macrophage transcriptional programs and bacterial killing that drive intestinal pathology. Pharmacological inhibition of TOP1/2 rescued these defects. Furthermore, exacerbated colitis was restored with TOP1/2 inhibitors in Sp140-/- mice, but not wild-type mice, in vivo. Collectively, we identify SP140 as a TOP repressor and reveal repurposing of TOP inhibition to reverse immune diseases driven by SP140 loss.

Polycomb repressive complex 2 shields naïve human pluripotent cells from trophectoderm differentiation

The first lineage choice in human embryo development separates trophectoderm from the inner cell mass. Naïve human embryonic stem cells are derived from the inner cell mass and offer possibilities to explore how lineage integrity is maintained. Here, we discover that polycomb repressive complex 2 (PRC2) maintains naïve pluripotency and restricts differentiation to trophectoderm and mesoderm lineages. Through quantitative epigenome profiling, we found that a broad gain of histone H3 lysine 27 trimethylation (H3K27me3) is a distinct feature of naïve pluripotency. We define shared and naïve-specific bivalent promoters featuring PRC2-mediated H3K27me3 concomitant with H3K4me3. Naïve bivalency maintains key trophectoderm and mesoderm transcription factors in a transcriptionally poised state. Inhibition of PRC2 forces naïve human embryonic stem cells into an 'activated' state, characterized by co-expression of pluripotency and lineage-specific transcription factors, followed by differentiation into either trophectoderm or mesoderm lineages. In summary, PRC2-mediated repression provides a highly adaptive mechanism to restrict lineage potential during early human development.

One genome, many cell states: epigenetic control of innate immunity

A hallmark of the innate immune system is its ability to rapidly initiate short-lived or sustained transcriptional programs in a cell-specific and pathogen-specific manner that is dependent on dynamic chromatin states. Much of the epigenetic landscape is set during cellular differentiation; however, pathogens and other environmental cues also induce changes in chromatin that can either promote tolerance or 'train' innate immune cells for amplified secondary responses. We review chromatin processes that enable innate immune cell differentiation and functional transcriptional responses in naive or experienced cells, in concert with signal transduction and cellular metabolic shifts. We discuss how immune chromatin mechanisms are maladapted in disease and novel therapeutic approaches for cellular reprogramming.

H3K27me3 conditions chemotolerance in triple-negative breast cancer

The persistence of cancer cells resistant to therapy remains a major clinical challenge. In triple-negative breast cancer, resistance to chemotherapy results in the highest recurrence risk among breast cancer subtypes. The drug-tolerant state seems largely defined by nongenetic features, but the underlying mechanisms are poorly understood. Here, by monitoring epigenomes, transcriptomes and lineages with single-cell resolution, we show that the repressive histone mark H3K27me3 (trimethylation of histone H3 at lysine 27) regulates cell fate at the onset of chemotherapy. We report that a persister expression program is primed with both H3K4me3 (trimethylation of histone H3 at lysine 4) and H3K27me3 in unchallenged cells, with H3K27me3 being the lock to its transcriptional activation. We further demonstrate that depleting H3K27me3 enhances the potential of cancer cells to tolerate chemotherapy. Conversely, preventing H3K27me3 demethylation simultaneously to chemotherapy inhibits the transition to a drug-tolerant state, and delays tumor recurrence in vivo. Our results highlight how chromatin landscapes shape the potential of cancer cells to respond to initial therapy.

Slow nucleosome dynamics set the transcriptional speed limit and induce RNA polymerase II traffic jams and bursts

Nucleosomes are recognized as key regulators of transcription. However, the relationship between slow nucleosome unwrapping dynamics and bulk transcriptional properties has not been thoroughly explored. Here, an agent-based model that we call the dynamic defect Totally Asymmetric Simple Exclusion Process (ddTASEP) was constructed to investigate the effects of nucleosome-induced pausing on transcriptional dynamics. Pausing due to slow nucleosome dynamics induced RNAPII convoy formation, which would cooperatively prevent nucleosome rebinding leading to bursts of transcription. The mean first passage time (MFPT) and the variance of first passage time (VFPT) were analytically expressed in terms of the nucleosome rate constants, allowing for the direct quantification of the effects of nucleosome-induced pausing on pioneering polymerase dynamics. The mean first passage elongation rate γ(h_c, h_o) is inversely proportional to the MFPT and can be considered to be a new axis of the ddTASEP phase diagram, orthogonal to the classical αβ-plane (where α and β are the initiation and termination rates). Subsequently, we showed that, for β = 1, there is a novel jamming transition in the αγ-plane that separates the ddTASEP dynamics into initiation-limited and nucleosome pausing-limited regions. We propose analytical estimates for the RNAPII density ρ, average elongation rate v, and transcription flux J and verified them numerically. We demonstrate that the intra-burst RNAPII waiting times t_in follow the time-headway distribution of a max flux TASEP and that the average inter-burst interval tIBI¯ correlates with the index of dispersion D_e. In the limit γ→0, the average burst size reaches a maximum set by the closing rate h_c. When α≪1, the burst sizes are geometrically distributed, allowing large bursts even while the average burst size NB¯ is small. Last, preliminary results on the relative effects of static and dynamic defects are presented to show that dynamic defects can induce equal or greater pausing than static bottle necks.

To perform specific functions, cells must express specific genes by copying the information in DNA into RNA via transcription. Structural proteins called nucleosomes are spaced every 200 base pairs along the length of a strand of DNA and play a crucial function in the regulation of gene activity by tightly binding DNA strands and condensing them into heterochromatin, preventing transcription by RNA polymerase II (RNAPII). Even on active genes where nucleosomes are loosely attached to DNA strands, the wrapping and unwrapping of nucleosomes pause transcription as RNAPII passes by. Previous mathematical models of transcription have compared this biological process to traffic on a one lane highway without obstructions. In contrast, our proposed model simulates transcription like traffic in a grid system where nucleosomes can be thought of as pedestrians or other vehicles crossing the road at regularly spaced intersections. Just as side street traffic and pedestrian crossings can cause cars to form convoys and cause jams limiting the max speed in an area, nucleosomes can cause RNAPII to form convoys that lead to bursts of mRNA production and limit the average polymerase flux through the gene.

Epigenetic basis of oncogenic-Kras-mediated epithelial-cellular proliferation and plasticity

Oncogenic Kras induces a hyper-proliferative state that permits cells to progress to neoplasms in diverse epithelial tissues. Depending on the cell of origin, this also involves lineage transformation. Although a multitude of downstream factors have been implicated in these processes, the precise chronology of molecular events controlling them remains elusive. Using mouse models, primary human tissues, and cell lines, we show that, in Kras-mutant alveolar type II cells (AEC2), FOSL1-based AP-1 factor guides the mSWI/SNF complex to increase chromatin accessibility at genomic loci controlling the expression of genes necessary for neoplastic transformation. We identified two orthogonal processes in Kras-mutant distal airway club cells. The first promoted their transdifferentiation into an AEC2-like state through NKX2.1, and the second controlled oncogenic transformation through the AP-1 complex. Our results suggest that neoplasms retain an epigenetic memory of their cell of origin through cell-type-specific transcription factors. Our analysis showed that a cross-tissue-conserved AP-1-dependent chromatin remodeling program regulates carcinogenesis.

ETV6-NCOA2 fusion induces T/myeloid mixed-phenotype leukemia through transformation of nonthymic hematopoietic progenitor cells

Mixed-phenotype acute leukemia is a rare subtype of leukemia in which both myeloid and lymphoid markers are co-expressed on the same malignant cells. The pathogenesis is largely unknown, and the treatment is challenging. We previously reported the specific association of the recurrent t(8;12)(q13;p13) chromosomal translocation that creates the ETV6-NCOA2 fusion with T/myeloid leukemias. Here we report that ETV6-NCOA2 initiates T/myeloid leukemia in preclinical models; ectopic expression of ETV6-NCOA2 in mouse bone marrow hematopoietic progenitors induced T/myeloid lymphoma accompanied by spontaneous Notch1-activating mutations. Similarly, cotransduction of human cord blood CD34+ progenitors with ETV6-NCOA2 and a nontransforming NOTCH1 mutant induced T/myeloid leukemia in immunodeficient mice; the immunophenotype and gene expression pattern were similar to those of patient-derived ETV6-NCOA2 leukemias. Mechanistically, we show that ETV6-NCOA2 forms a transcriptional complex with ETV6 and the histone acetyltransferase p300, leading to derepression of ETV6 target genes. The expression of ETV6-NCOA2 in human and mouse nonthymic hematopoietic progenitor cells induces transcriptional dysregulation, which activates a lymphoid program while failing to repress the expression of myeloid genes such as CSF1 and MEF2C. The ETV6-NCOA2 induced arrest at an early immature T-cell developmental stage. The additional acquisition of activating NOTCH1 mutations transforms the early immature ETV6-NCOA2 cells into T/myeloid leukemias. Here, we describe the first preclinical model to depict the initiation of T/myeloid leukemia by a specific somatic genetic aberration.

Profiling Histone Methylation in Low Numbers of Cells

Chromatin immunoprecipitation (ChIP) enables the study of DNA-protein interactions. When coupled with high-throughput sequencing (ChIP-seq), this method allows the generation of genome-wide profiles of the distribution of specific proteins in a given cellular context. Typical ChIP-seq experiments require millions of cells as input material and thus are not ideal to study many in vivo cell populations. Here, we describe an ultra-low-input native ChIP-seq method, ULI-NChIP-seq, to profile histone modification patterns in as low as 150 cells.

High-throughput single-cell epigenomic profiling by targeted insertion of promoters (TIP-seq)

Chromatin profiling in single cells has been extremely challenging and almost exclusively limited to histone proteins. In cases where single-cell methods have shown promise, many require highly specialized equipment or cell type-specific protocols and are relatively low throughput. Here, we combine the advantages of tagmentation, linear amplification, and combinatorial indexing to produce a high-throughput single-cell DNA binding site mapping method that is simple, inexpensive, and capable of multiplexing several independent samples per experiment. Targeted insertion of promoters sequencing (TIP-seq) uses Tn5 fused to proteinA to insert a T7 RNA polymerase promoter adjacent to a chromatin protein of interest. Linear amplification of flanking DNA with T7 polymerase before sequencing library preparation provides ∼10-fold higher unique reads per single cell compared with other methods. We applied TIP-seq to map histone modifications, RNA polymerase II (RNAPII), and transcription factor CTCF binding sites in single human and mouse cells.

FACT-seq: profiling histone modifications in formalin-fixed paraffin-embedded samples with low cell numbers

The majority of biopsies in both basic research and translational cancer studies are preserved in the format of archived formalin-fixed paraffin-embedded (FFPE) samples. Profiling histone modifications in archived FFPE tissues is critically important to understand gene regulation in human disease. The required input for current genome-wide histone modification profiling studies from FFPE samples is either 10-20 tissue sections or whole tissue blocks, which prevents better resolved analyses. But it is desirable to consume a minimal amount of FFPE tissue sections in the analysis as clinical tissues of interest are limited. Here, we present FFPE tissue with antibody-guided chromatin tagmentation with sequencing (FACT-seq), the first highly sensitive method to efficiently profile histone modifications in FFPE tissues by combining a novel fusion protein of hyperactive Tn5 transposase and protein A (T7-pA-Tn5) transposition and T7 in vitro transcription. FACT-seq generates high-quality chromatin profiles from different histone modifications with low number of FFPE nuclei. We proved a very small piece of FFPE tissue section containing ∼4000 nuclei is sufficient to decode H3K27ac modifications with FACT-seq. H3K27ac FACT-seq revealed disease-specific super enhancers in the archived FFPE human colorectal and human glioblastoma cancer tissue. In summary, FACT-seq allows decoding the histone modifications in archival FFPE tissues with high sensitivity and help researchers to better understand epigenetic regulation in cancer and human disease.

Establishment, maintenance, and recall of inflammatory memory

Known for nearly a century but through mechanisms that remain elusive, cells retain a memory of inflammation that equips them to react quickly and broadly to diverse secondary stimuli. Using murine epidermal stem cells as a model, we elucidate how cells establish, maintain, and recall inflammatory memory. Specifically, we landscape and functionally interrogate temporal, dynamic changes to chromatin accessibility, histone modifications, and transcription factor binding that occur during inflammation, post-resolution, and in memory recall following injury. We unearth an essential, unifying role for the general stress-responsive transcription factor FOS, which partners with JUN and cooperates with stimulus-specific STAT3 to establish memory; JUN then remains with other homeostatic factors on memory domains, facilitating rapid FOS re-recruitment and gene re-activation upon diverse secondary challenges. Extending our findings, we offer a comprehensive, potentially universal mechanism behind inflammatory memory and less discriminate recall phenomena with profound implications for tissue fitness in health and disease.

Chromatin Remodeling of Colorectal Cancer Liver Metastasis is Mediated by an HGF-PU.1-DPP4 Axis

Colorectal cancer (CRC) metastasizes mainly to the liver, which accounts for the majority of CRC-related deaths. Here it is shown that metastatic cells undergo specific chromatin remodeling in the liver. Hepatic growth factor (HGF) induces phosphorylation of PU.1, a pioneer factor, which in turn binds and opens chromatin regions of downstream effector genes. PU.1 increases histone acetylation at the DPP4 locus. Precise epigenetic silencing by CRISPR/dCas9KRAB or CRISPR/dCas9HDAC revealed that individual PU.1-remodeled regulatory elements collectively modulate DPP4 expression and liver metastasis growth. Genetic silencing or pharmacological inhibition of each factor along this chromatin remodeling axis strongly suppressed liver metastasis. Therefore, microenvironment-induced epimutation is an important mechanism for metastatic tumor cells to grow in their new niche. This study presents a potential strategy to target chromatin remodeling in metastatic cancer and the promise of repurposing drugs to treat metastasis.

Recent advance in single-cell epigenomics

The analysis of gene expression regulation, or the epigenome analysis, at the single-cell level is at the forefront of genomics research. To elucidate the mechanisms that regulate gene expression, chromatin immunoprecipitation has been conventionally used for determining the binding sites of DNA-binding proteins, such as histones and transcription factors. Now several new approaches have been emerged to reveal epigenome states at the single-cell level. Instead of using immunoprecipitation of fragmented chromatin, in situ reactions using cells or nuclei, combining with transposase tagging and other methods, have enabled single-cell analysis. Furthermore, single-cell multiomics techniques to simultaneously profiling transcriptome and open chromatin or histone modification have been developed. These single-cell analyses have the potential to identify different cell types in a cell population and reveal the dynamic changes of gene regulation, although those technologies have not yet reached a level for general application.

Epigenetic silencing by SETDB1 suppresses tumour intrinsic immunogenicity

Epigenetic dysregulation is a defining feature of tumorigenesis that is implicated in immune escape^1,2. Here, to identify factors that modulate the immune sensitivity of cancer cells, we performed in vivo CRISPR-Cas9 screens targeting 936 chromatin regulators in mouse tumour models treated with immune checkpoint blockade. We identified the H3K9 methyltransferase SETDB1 and other members of the HUSH and KAP1 complexes as mediators of immune escape^3-5. We also found that amplification of SETDB1 (1q21.3) in human tumours is associated with immune exclusion and resistance to immune checkpoint blockade. SETDB1 represses broad domains, primarily within the open genome compartment. These domains are enriched for transposable elements (TEs) and immune clusters associated with segmental duplication events, a central mechanism of genome evolution⁶. SETDB1 loss derepresses latent TE-derived regulatory elements, immunostimulatory genes, and TE-encoded retroviral antigens in these regions, and triggers TE-specific cytotoxic T cell responses in vivo. Our study establishes SETDB1 as an epigenetic checkpoint that suppresses tumour-intrinsic immunogenicity, and thus represents a candidate target for immunotherapy.

Transcription factor chromatin profiling genome-wide using uliCUT&RUN in single cells and individual blastocysts

Determining chromatin-associated protein localization across the genome has provided insight into the functions of DNA-binding proteins and their connections to disease. However, established protocols requiring large quantities of cell or tissue samples currently limit applications for clinical and biomedical research in this field. Furthermore, most technologies have been optimized to assess abundant histone protein localization, prohibiting the investigation of nonhistone protein localization in low cell numbers. We recently described a protocol to profile chromatin-associated protein localization in as low as one cell: ultra-low-input cleavage under targets and release using nuclease (uliCUT&RUN). Optimized from chromatin immunocleavage and CUT&RUN, uliCUT&RUN is a tethered enzyme-based protocol that utilizes a combination of recombinant protein, antibody recognition and stringent purification to selectively target proteins of interest and isolate the associated DNA. Performed in native conditions, uliCUT&RUN profiles protein localization to chromatin with low input and high precision. Compared with other profiling technologies, uliCUT&RUN can determine nonhistone protein chromatin occupancies in low cell numbers, permitting the investigation into the molecular functions of a range of DNA-binding proteins within rare samples. From sample preparation to sequencing library submission, the uliCUT&RUN protocol takes <2 d to perform, with the accompanying data analysis timeline dependent on experience level.

A plug and play microfluidic platform for standardized sensitive low-input chromatin immunoprecipitation

Epigenetic profiling by chromatin immunoprecipitation followed by sequencing (ChIP-seq) has become a powerful tool for genome-wide identification of regulatory elements, for defining transcriptional regulatory networks, and for screening for biomarkers. However, the ChIP-seq protocol for low-input samples is laborious and time-consuming and suffers from experimental variation, resulting in poor reproducibility and low throughput. Although prototypic microfluidic ChIP-seq platforms have been developed, these are poorly transferable as they require sophisticated custom-made equipment and in-depth microfluidic and ChIP expertise, while lacking parallelization. To enable standardized, automated ChIP-seq profiling of low-input samples, we constructed microfluidic PDMS-based plates capable of performing 24 sensitive ChIP reactions within 30 min of hands-on time and 4.5 h of machine-running time. These disposable plates can be conveniently loaded into a widely available controller for pneumatics and thermocycling. In light of the plug and play (PnP) ChIP plates and workflow, we named our procedure PnP-ChIP-seq. We show high-quality ChIP-seq on hundreds to a few thousand of cells for all six post-translational histone modifications that are included in the International Human Epigenome Consortium set of reference epigenomes. PnP-ChIP-seq robustly detects epigenetic differences on promoters and enhancers between naive and more primed mouse embryonic stem cells (mESCs). Furthermore, we used our platform to generate epigenetic profiles of rare subpopulations of mESCs that resemble the two-cell stage of embryonic development. PnP-ChIP-seq allows nonexpert laboratories worldwide to conveniently run robust, standardized ChIP-seq, whereas its high throughput, consistency, and sensitivity pave the way toward large-scale profiling of precious sample types such as rare subpopulations of cells or biopsies.

Smc3 dosage regulates B cell transit through germinal centers and restricts their malignant transformation

During the germinal center (GC) reaction, B cells undergo extensive redistribution of cohesin complex and three-dimensional reorganization of their genomes. Yet, the significance of cohesin and architectural programming in the humoral immune response is unknown. Herein we report that homozygous deletion of Smc3, encoding the cohesin ATPase subunit, abrogated GC formation, while, in marked contrast, Smc3 haploinsufficiency resulted in GC hyperplasia, skewing of GC polarity and impaired plasma cell (PC) differentiation. Genome-wide chromosomal conformation and transcriptional profiling revealed defects in GC B cell terminal differentiation programs controlled by the lymphoma epigenetic tumor suppressors Tet2 and Kmt2d and failure of Smc3-haploinsufficient GC B cells to switch from B cell- to PC-defining transcription factors. Smc3 haploinsufficiency preferentially impaired the connectivity of enhancer elements controlling various lymphoma tumor suppressor genes, and, accordingly, Smc3 haploinsufficiency accelerated lymphomagenesis in mice with constitutive Bcl6 expression. Collectively, our data indicate a dose-dependent function for cohesin in humoral immunity to facilitate the B cell to PC phenotypic switch while restricting malignant transformation.

Evolution of the Epigenetic Landscape in Childhood B Acute Lymphoblastic Leukemia and Its Role in Drug Resistance

Although B-cell acute lymphoblastic leukemia (B-ALL) is the most common malignancy in children and while highly curable, it remains a leading cause of cancer-related mortality. The outgrowth of tumor subclones carrying mutations in genes responsible for resistance to therapy has led to a Darwinian model of clonal selection. Previous work has indicated that alterations in the epigenome might contribute to clonal selection, yet the extent to which the chromatin state is altered under the selective pressures of therapy is unknown. To address this, we performed chromatin immunoprecipitation, gene expression analysis, and enhanced reduced representation bisulfite sequencing on a cohort of paired diagnosis and relapse samples from individual patients who all but one relapsed within 36 months of initial diagnosis. The chromatin state at diagnosis varied widely among patients, while the majority of peaks remained stable between diagnosis and relapse. Yet a significant fraction was either lost or newly gained, with some patients showing few differences and others showing massive changes of the epigenetic state. Evolution of the epigenome was associated with pathways previously linked to therapy resistance as well as novel candidate pathways through alterations in pyrimidine biosynthesis and downregulation of polycomb repressive complex 2 targets. Three novel, relapse-specific superenhancers were shared by a majority of patients including one associated with S100A8, the top upregulated gene seen at relapse in childhood B-ALL. Overall, our results support a role of the epigenome in clonal evolution and uncover new candidate pathways associated with relapse. SIGNIFICANCE: This study suggests a major role for epigenetic mechanisms in driving clonal evolution in B-ALL and identifies novel pathways associated with drug resistance.

DOT1L-mediated murine neuronal differentiation associates with H3K79me2 accumulation and preserves SOX2-enhancer accessibility

During neuronal differentiation, the transcriptional profile and the epigenetic context of neural committed cells is subject to significant rearrangements, but a systematic quantification of global histone modification changes is still missing. Here, we show that H3K79me2 increases and H3K27ac decreases globally during in-vitro neuronal differentiation of murine embryonic stem cells. DOT1L mediates all three degrees of methylation of H3K79 and its enzymatic activity is critical to modulate cellular differentiation and reprogramming. In this context, we find that inhibition of DOT1L in neural progenitor cells biases the transcriptional state towards neuronal differentiation, resulting in transcriptional upregulation of genes marked with H3K27me3 on the promoter region. We further show that DOT1L inhibition affects accessibility of SOX2-bound enhancers and impairs SOX2 binding in neural progenitors. Our work provides evidence that DOT1L activity gates differentiation of progenitors by allowing SOX2-dependent transcription of stemness programs.

Targeting DNA Methylation Depletes Uterine Leiomyoma Stem Cell-enriched Population by Stimulating Their Differentiation

Uterine leiomyoma (LM) is the most common tumor in women and can cause severe morbidity. Leiomyoma growth requires the maintenance and proliferation of a stem cell population. Dysregulated deoxyribonucleic acid (DNA) methylation has been reported in LM, but its role in LM stem cell regulation remains unclear. Here, we fluorescence-activated cell sorting (FACS)-sorted cells from human LM tissues into 3 populations: LM stem cell-like cells (LSC, 5%), LM intermediate cells (LIC, 7%), and differentiated LM cells (LDC, 88%), and we analyzed the transcriptome and epigenetic landscape of LM cells at different differentiation stages. Leiomyoma stem cell-like cells harbored a unique methylome, with 8862 differentially methylated regions compared to LIC and 9444 compared to LDC, most of which were hypermethylated. Consistent with global hypermethylation, transcript levels of TET1 and TET3 methylcytosine dioxygenases were lower in LSC. Integrative analyses revealed an inverse relationship between methylation and gene expression changes during LSC differentiation. In LSC, hypermethylation suppressed the genes important for myometrium- and LM-associated functions, including muscle contraction and hormone action, to maintain stemness. The hypomethylating drug, 5'-Aza, stimulated LSC differentiation, depleting the stem cell population and inhibiting tumor initiation. Our data suggest that DNA methylation maintains the pool of LSC, which is critical for the regeneration of LM tumors.

Chromatin integration labeling for mapping DNA-binding proteins and modifications with low input

Cell identity is determined by the selective activation or silencing of specific genes via transcription factor binding and epigenetic modifications on the genome. Chromatin immunoprecipitation (ChIP) has been the standard technique for mapping the sites of transcription factor binding and histone modification. Recently, alternative methods to ChIP have been developed for addressing the increasing demands for low-input epigenomic profiling. Chromatin integration labeling (ChIL) followed by sequencing (ChIL-seq) has been demonstrated to be particularly useful for epigenomic profiling of low-input samples or even single cells because the technique amplifies the target genomic sequence before cell lysis. After labeling the target protein or modification in situ with an oligonucleotide-conjugated antibody (ChIL probe), the nearby genome sequence is amplified by Tn5 transposase-mediated transposition followed by T7 RNA polymerase-mediated transcription. ChIL-seq enables the detection of the antibody target localization under a fluorescence microscope and at the genomic level. Here we describe the detailed protocol of ChIL-seq with assessment methods for the key steps, including ChIL probe reaction, transposition, in situ transcription and sequencing library preparation. The protocol usually takes 3 d to prepare the sequencing library, including overnight incubations for the ChIL probe reaction and in situ transcription. The ChIL probe can be separately prepared and stored for several months, and its preparation and evaluation protocols are also documented in detail. An optional analysis for multiple targets (multitarget ChIL-seq) is also described. We anticipate that the protocol presented here will make the ChIL technique more widely accessible for analyzing precious samples and facilitate further applications.

A physical basis for quantitative ChIP-sequencing

ChIP followed by next-generation sequencing (ChIP-Seq) is a key technique for mapping the distribution of histone posttranslational modifications (PTMs) and chromatin-associated factors across genomes. There is a perceived challenge to define a quantitative scale for ChIP-Seq data, and as such, several approaches making use of exogenous additives, or "spike-ins," have recently been developed. Herein, we report on the development of a quantitative, physical model defining ChIP-Seq. The quantitative scale on which ChIP-Seq results should be compared emerges from the model. To test the model and demonstrate the quantitative scale, we examine the impacts of an EZH2 inhibitor through the lens of ChIP-Seq. We report a significant increase in immunoprecipitation of presumed off-target histone PTMs after inhibitor treatment, a trend predicted by the model but contrary to spike-in-based indications. Our work also identifies a sensitivity issue in spike-in normalization that has not been considered in the literature, placing limitations on its utility and trustworthiness. We call our new approach the sans-spike-in method for quantitative ChIP-sequencing (siQ-ChIP). A number of changes in community practice of ChIP-Seq, data reporting, and analysis are motivated by this work.

Epigenetic fine-mapping: identification of causal mechanisms for autoimmunity

Genome-wide association studies (GWAS) have identified genetic susceptibility loci for a variety of autoimmune and inflammatory diseases. These studies confirm the fundamental genetic basis of individual autoimmune diseases, and also point to shared etiological mechanisms across the spectrum of autoimmunity. While hundreds of genetic loci have been implicated in autoimmune diseases, the translation of individual susceptibility loci into specific molecular mechanisms for individual diseases remains difficult. This review highlights recent advances in the genetics of autoimmune disease, and the emerging use of epigenetic techniques to identify pathogenic cell types and causal molecular mechanisms of autoimmunity.

Persistence of a regeneration-associated, transitional alveolar epithelial cell state in pulmonary fibrosis

Stem cells undergo dynamic changes in response to injury to regenerate lost cells. However, the identity of transitional states and the mechanisms that drive their trajectories remain understudied. Using lung organoids, multiple in vivo repair models, single-cell transcriptomics and lineage tracing, we find that alveolar type-2 epithelial cells undergoing differentiation into type-1 cells acquire pre-alveolar type-1 transitional cell state (PATS) en route to terminal maturation. Transitional cells undergo extensive stretching during differentiation, making them vulnerable to DNA damage. Cells in the PATS show an enrichment of TP53, TGFβ, DNA-damage-response signalling and cellular senescence. Gain and loss of function as well as genomic binding assays revealed a direct transcriptional control of PATS by TP53 signalling. Notably, accumulation of PATS-like cells in human fibrotic lungs was observed, suggesting persistence of the transitional state in fibrosis. Our study thus implicates a transient state associated with senescence in normal epithelial tissue repair and its abnormal persistence in disease conditions.

AutoRELACS: automated generation and analysis of ultra-parallel ChIP-seq

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a method used to profile protein-DNA interactions genome-wide. Restriction Enzyme-based Labeling of Chromatin in Situ (RELACS) is a recently developed ChIP-seq protocol that deploys a chromatin barcoding strategy to enable standardized and high-throughput generation of ChIP-seq data. The manual implementation of RELACS is constrained by human processivity in both data generation and data analysis. To overcome these limitations, we have developed AutoRELACS, an automated implementation of the RELACS protocol using the liquid handler Biomek i7 workstation. We match the unprecedented processivity in data generation allowed by AutoRELACS with the automated computation pipelines offered by snakePipes. In doing so, we build a continuous workflow that streamlines epigenetic profiling, from sample collection to biological interpretation. Here, we show that AutoRELACS successfully automates chromatin barcode integration, and is able to generate high-quality ChIP-seq data comparable with the standards of the manual protocol, also for limited amounts of biological samples.

A parallelized, automated platform enabling individual or sequential ChIP of histone marks and transcription factors

Despite its popularity, chromatin immunoprecipitation followed by sequencing (ChIP-seq) remains a tedious (>2 d), manually intensive, low-sensitivity and low-throughput approach. Here, we combine principles of microengineering, surface chemistry, and molecular biology to address the major limitations of standard ChIP-seq. The resulting technology, FloChIP, automates and miniaturizes ChIP in a beadless fashion while facilitating the downstream library preparation process through on-chip chromatin tagmentation. FloChIP is fast (<2 h), has a wide dynamic range (from 10⁶ to 500 cells), is scalable and parallelized, and supports antibody- or sample-multiplexed ChIP on both histone marks and transcription factors. In addition, FloChIP's interconnected design allows for straightforward chromatin reimmunoprecipitation, which allows this technology to also act as a microfluidic sequential ChIP-seq system. Finally, we ran FloChIP for the transcription factor MEF2A in 32 distinct human lymphoblastoid cell lines, providing insights into the main factors driving collaborative DNA binding of MEF2A and into its role in B cell-specific gene regulation. Together, our results validate FloChIP as a flexible and reproducible automated solution for individual or sequential ChIP-seq.

Loci specific epigenetic drug sensitivity

Therapeutic targeting of epigenetic modulators offers a novel approach to the treatment of multiple diseases. The cellular consequences of chemical compounds that target epigenetic regulators (epi-drugs) are complex. Epi-drugs affect global cellular phenotypes and cause local changes to gene expression due to alteration of a gene chromatin environment. Despite increasing use in the clinic, the mechanisms responsible for cellular changes are unclear. Specifically, to what degree the effects are a result of cell-wide changes or disease related locus specific effects is unknown. Here we developed a platform to systematically and simultaneously investigate the sensitivity of epi-drugs at hundreds of genomic locations by combining DNA barcoding, unique split-pool encoding, and single cell expression measurements. Internal controls are used to isolate locus specific effects separately from any global consequences these drugs have. Using this platform we discovered wide-spread loci specific sensitivities to epi-drugs for three distinct epi-drugs that target histone deacetylase, DNA methylation and bromodomain proteins. By leveraging ENCODE data on chromatin modification, we identified features of chromatin environments that are most likely to be affected by epi-drugs. The measurements of loci specific epi-drugs sensitivities will pave the way to the development of targeted therapy for personalized medicine.

Canonical Wnt/β-catenin activity and differential epigenetic marks direct sexually dimorphic regulation of Irx3 and Irx5 in developing mouse gonads

Members of the Iroquois B (IrxB) homeodomain cluster genes, specifically Irx3 and Irx5, are crucial for heart, limb and bone development. Recently, we reported their importance for oocyte and follicle survival within the developing ovary. Irx3 and Irx5 expression begins after sex determination in the ovary but remains absent in the fetal testis. Mutually antagonistic molecular signals ensure ovary versus testis differentiation with canonical Wnt/β-catenin signals paramount for promoting the ovary pathway. Notably, few direct downstream targets have been identified. We report that Wnt/β-catenin signaling directly stimulates Irx3 and Irx5 transcription in the developing ovary. Using in silico analysis of ATAC- and ChIP-Seq databases in conjunction with mouse gonad explant transfection assays, we identified TCF/LEF-binding sequences within two distal enhancers of the IrxB locus that promote β-catenin-responsive ovary expression. Meanwhile, Irx3 and Irx5 transcription is suppressed within the developing testis by the presence of H3K27me3 on these same sites. Thus, we resolved sexually dimorphic regulation of Irx3 and Irx5 via epigenetic and β-catenin transcriptional control where their ovarian presence promotes oocyte and follicle survival vital for future ovarian health.