The transcriptional repressor B lymphocyte-induced maturation protein-1 (BLIMP1) regulates gene expression and cell fate. The DNA motif bound by BLIMP1 in vitro overlaps with that of interferon regulatory factors (IRFs), which respond to inflammatory/immune signals. At such sites, BLIMP1 and IRFs can antagonistically regulate promoter activity. In vitro motif selection predicts that only a subset of BLIMP1 or IRF sites is subject to antagonistic regulation, but the extent to which antagonism occurs is unknown, since an unbiased assessment of BLIMP1 occupancy in vivo is lacking. To address this, we identified an extended set of promoters occupied by BLIMP1. Motif discovery and enrichment analysis demonstrate that multiple motif variants are required to capture BLIMP1 binding specificity. These are differentially associated with CpG content, leading to the observation that BLIMP1 DNA-binding is methylation sensitive. In occupied promoters, only a subset of BLIMP1 motifs overlap with IRF motifs. Conversely, a distinct subset of IRF motifs is not enriched amongst occupied promoters. Genes linked to occupied promoters containing overlapping BLIMP1/IRF motifs (e.g. AIM2, SP110, BTN3A3) are shown to constitute a dynamic target set which is preferentially activated by BLIMP1 knock-down. These data confirm and extend the competitive model of BLIMP1 and IRF interaction.

Mammalian development is regulated by the interplay of tissue-specific and ubiquitously expressed transcription factors, such as Sp1. Sp1 knockout mice die in utero with multiple phenotypic aberrations, but the underlying molecular mechanism of this differentiation failure has been elusive. Here, we have used conditional knockout mice as well as the differentiation of mouse ES cells as a model with which to address this issue. To this end, we examined differentiation potential, global gene expression patterns and Sp1 target regions in Sp1 wild-type and Sp1-deficient cells representing different stages of hematopoiesis. Sp1(-/-) cells progress through most embryonic stages of blood cell development but cannot complete terminal differentiation. This failure to fully differentiate is not seen when Sp1 is knocked out at later developmental stages. For most Sp1 target and non-target genes, gene expression is unaffected by Sp1 inactivation. However, Cdx genes and multiple Hox genes are stage-specific targets of Sp1 and are downregulated at an early stage. As a consequence, expression of genes involved in hematopoietic specification is progressively deregulated. Our work demonstrates that the early absence of active Sp1 sets a cascade in motion that culminates in a failure of terminal hematopoietic differentiation and emphasizes the role of ubiquitously expressed transcription factors for tissue-specific gene regulation. In addition, our global side-by-side analysis of the response of the transcriptional network to perturbation sheds a new light on the regulatory hierarchy of hematopoietic specification.

Immunological memory is a defining feature of vertebrate physiology, allowing rapid responses to repeat infections. However, the molecular mechanisms required for its establishment and maintenance remain poorly understood. Here, we demonstrated that the first steps in the acquisition of T-cell memory occurred during the initial activation phase of naïve T cells by an antigenic stimulus. This event initiated extensive chromatin remodeling that reprogrammed immune response genes toward a stably maintained primed state, prior to terminal differentiation. Activation induced the transcription factors NFAT and AP-1 which created thousands of new DNase I-hypersensitive sites (DHSs), enabling ETS-1 and RUNX1 recruitment to previously inaccessible sites. Significantly, these DHSs remained stable long after activation ceased, were preserved following replication, and were maintained in memory-phenotype cells. We show that primed DHSs maintain regions of active chromatin in the vicinity of inducible genes and enhancers that regulate immune responses. We suggest that this priming mechanism may contribute to immunological memory in T cells by facilitating the induction of nearby inducible regulatory elements in previously activated T cells.

Acute myeloid leukemia (AML) is characterized by recurrent mutations that affect the epigenetic regulatory machinery and signaling molecules, leading to a block in hematopoietic differentiation. Constitutive signaling from mutated growth factor receptors is a major driver of leukemic growth, but how aberrant signaling affects the epigenome in AML is less understood. Furthermore, AML cells undergo extensive clonal evolution, and the mutations in signaling genes are often secondary events. To elucidate how chronic growth factor signaling alters the transcriptional network in AML, we performed a system-wide multi-omics study of primary cells from patients suffering from AML with internal tandem duplications in the FLT3 transmembrane domain (FLT3-ITD). This strategy revealed cooperation between the MAP kinase (MAPK) inducible transcription factor AP-1 and RUNX1 as a major driver of a common, FLT3-ITD-specific gene expression and chromatin signature, demonstrating a major impact of MAPK signaling pathways in shaping the epigenome of FLT3-ITD AML.

Acute myeloid leukemia (AML) is a highly heterogeneous cancer associated with different patterns of gene expression determined by the nature of their DNA mutations. These mutations mostly act to deregulate gene expression by various mechanisms at the level of the nucleus. By performing genome-wide epigenetic profiling of cis-regulatory elements, we found that AML encompasses different mutation-specific subclasses associated with the rewiring of the gene regulatory networks that drive differentiation into different directions away from normal myeloid development. By integrating epigenetic profiles with gene expression and chromatin conformation data, we defined pathways within gene regulation networks that were differentially rewired within each mutation-specific subclass of AML. This analysis revealed 2 major classes of AML: one class defined by mutations in signaling molecules that activate AP-1 via the mitogen-activated protein (MAP) kinase pathway and a second class defined by mutations within genes encoding transcription factors such as RUNX1/CBFβ and C/EBPα. By identifying specific DNA motifs protected from DNase I digestion at cis-regulatory elements, we were able to infer candidate transcription factors bound to these motifs. These integrated analyses allowed the identification of AML subtype-specific core regulatory networks that are required for AML development and maintenance, which could now be targeted in personalized therapies.

This protocol offers a detailed procedure for the in vitro differentiation of human pluripotent stem cells (hPSCs) to multipotent hematopoietic progenitors that arise from SOX17+ hemogenic endothelium, mimicking intra-embryonic, HOXA-positive, aorta-gonad mesonephros (AGM) hematopoiesis. The generated endothelium displays transcriptional similarities to cells sorted from human 5-week AGM, and CD45+CD34+RUNX1C+ progenitors share an accessible chromatin profile with adult hematopoietic stem cells and multipotent progenitors. Therefore, this protocol is suitable for the mechanistic study of human multipotent progenitor development and for modeling childhood leukemias. For complete details on the use and execution of this protocol, please refer to Nafria et al. (2020).

Acute Myeloid Leukemia (AML) is caused by multiple mutations which dysregulate growth and differentiation of myeloid cells. Cells adopt different gene regulatory networks specific to individual mutations, maintaining a rapidly proliferating blast cell population with fatal consequences for the patient if not treated. The most common treatment option is still chemotherapy which targets such cells. However, patients harbour a population of quiescent leukemic stem cells (LSCs) which can emerge from quiescence to trigger relapse after therapy. The processes that allow such cells to re-grow remain unknown. Here, we examine the well characterised t(8;21) AML sub-type as a model to address this question. Using four primary AML samples and a novel t(8;21) patient-derived xenograft model, we show that t(8;21) LSCs aberrantly activate the VEGF and IL-5 signalling pathways. Both pathways operate within a regulatory circuit consisting of the driver oncoprotein RUNX1::ETO and an AP-1/GATA2 axis allowing LSCs to re-enter the cell cycle while preserving self-renewal capacity.

Intra-tumor heterogeneity caused by clonal evolution is a major problem in cancer treatment. To address this problem, we performed label-free quantitative proteomics on primary acute myeloid leukemia (AML) samples. We identified 50 leukemia-enriched plasma membrane proteins enabling the prospective isolation of genetically distinct subclones from individual AML patients. Subclones differed in their regulatory phenotype, drug sensitivity, growth, and engraftment behavior, as determined by RNA sequencing, DNase I hypersensitive site mapping, transcription factor occupancy analysis, in vitro culture, and xenograft transplantation. Finally, we show that these markers can be used to identify and longitudinally track distinct leukemic clones in patients in routine diagnostics. Our study describes a strategy for a major improvement in stratifying cancer diagnosis and treatment.

Oncogenic transcription factors such as RUNX1/ETO, which is generated by the chromosomal translocation t(8;21), subvert normal blood cell development by impairing differentiation and driving malignant self-renewal. Here, we use digital footprinting and chromatin immunoprecipitation sequencing (ChIP-seq) to identify the core RUNX1/ETO-responsive transcriptional network of t(8;21) cells. We show that the transcriptional program underlying leukemic propagation is regulated by a dynamic equilibrium between RUNX1/ETO and RUNX1 complexes, which bind to identical DNA sites in a mutually exclusive fashion. Perturbation of this equilibrium in t(8;21) cells by RUNX1/ETO depletion leads to a global redistribution of transcription factor complexes within preexisting open chromatin, resulting in the formation of a transcriptional network that drives myeloid differentiation. Our work demonstrates on a genome-wide level that the extent of impaired myeloid differentiation in t(8;21) is controlled by the dynamic balance between RUNX1/ETO and RUNX1 activities through the repression of transcription factors that drive differentiation.

Acute myeloid leukemia (AML) is a heterogeneous disease caused by mutations in transcriptional regulator genes, but how different mutant regulators shape the chromatin landscape is unclear. Here, we compared the transcriptional networks of two types of AML with chromosomal translocations of the RUNX1 locus that fuse the RUNX1 DNA-binding domain to different regulators, the t(8;21) expressing RUNX1-ETO and the t(3;21) expressing RUNX1-EVI1. Despite containing the same DNA-binding domain, the two fusion proteins display distinct binding patterns, show differences in gene expression and chromatin landscape, and are dependent on different transcription factors. RUNX1-EVI1 directs a stem cell-like transcriptional network reliant on GATA2, whereas that of RUNX1-ETO-expressing cells is more mature and depends on RUNX1. However, both types of AML are dependent on the continuous expression of the fusion proteins. Our data provide a molecular explanation for the differences in clinical prognosis for these types of AML.

Haematopoietic cells arise from endothelial cells within the dorsal aorta of the embryo via a process called the endothelial-haematopoietic transition (EHT). This process crucially depends on the transcription factor RUNX1 which rapidly activates the expression of genes essential for haematopoietic development. Using an inducible version of RUNX1 in a mouse embryonic stem cell differentiation model we showed that prior to the EHT, haematopoietic genes are primed by the binding of the transcription factor FLI1. Once expressed, RUNX1 relocates FLI1 towards its binding sites. However, the nature of the transcription factor assemblies recruited by RUNX1 to reshape the chromatin landscape and initiate mRNA synthesis are unclear. Here, we performed genome-wide analyses of RUNX1-dependent binding of factors associated with transcription elongation to address this question. We demonstrate that RUNX1 induction moves FLI1 from distal ETS/GATA sites to RUNX1/ETS sites and recruits the basal transcription factors CDK9, BRD4, the Mediator complex and the looping factor LDB1. Our study explains how the expression of a single transcription factor can drive rapid and replication independent transitions in cellular shape which are widely observed in development and disease.

Acute myeloid leukemia (AML) is a hematopoietic malignancy caused by recurrent mutations in genes encoding transcriptional, chromatin, and/or signaling regulators. The t(8;21) translocation generates the aberrant transcription factor RUNX1-ETO (RUNX1-RUNX1T1), which by itself is insufficient to cause disease. t(8;21) AML patients show extensive chromatin reprogramming and have acquired additional mutations. Therefore, the genomic and developmental effects directly and solely attributable to RUNX1-ETO expression are unclear. To address this, we employ a human embryonic stem cell differentiation system capable of forming definitive myeloid progenitor cells to express RUNX1-ETO in an inducible fashion. Induction of RUNX1-ETO causes extensive chromatin reprogramming by interfering with RUNX1 binding, blocks differentiation, and arrests cellular growth, whereby growth arrest is reversible following RUNX1-ETO removal. Single-cell gene expression analyses show that RUNX1-ETO induction alters the differentiation of early myeloid progenitors, but not of other progenitor types, indicating that oncoprotein-mediated transcriptional reprogramming is highly target cell specific.

Gene expression governs cell fate, and is regulated via a complex interplay of transcription factors and molecules that change chromatin structure. Advances in sequencing-based assays have enabled investigation of these processes genome-wide, leading to large datasets that combine information on the dynamics of gene expression, transcription factor binding and chromatin structure as cells differentiate. While numerous studies focus on the effects of these features on broader gene regulation, less work has been done on the mechanisms of gene-specific transcriptional control. In this study, we have focussed on the latter by integrating gene expression data for the in vitro differentiation of murine ES cells to macrophages and cardiomyocytes, with dynamic data on chromatin structure, epigenetics and transcription factor binding. Combining a novel strategy to identify communities of related control elements with a penalized regression approach, we developed individual models to identify the potential control elements predictive of the expression of each gene. Our models were compared to an existing method and evaluated using the existing literature and new experimental data from embryonic stem cell differentiation reporter assays. Our method is able to identify transcriptional control elements in a gene specific manner that reflect known regulatory relationships and to generate useful hypotheses for further testing.

The transmission of extracellular signals into the nucleus involves inducible transcription factors, but how different signalling pathways act in a cell type-specific fashion is poorly understood. Here, we studied the regulatory role of the AP-1 transcription factor family in blood development using embryonic stem cell differentiation coupled with genome-wide transcription factor binding and gene expression analyses. AP-1 factors respond to MAP kinase signalling and comprise dimers of FOS, ATF and JUN proteins. To examine genes regulated by AP-1 and to examine how it interacts with other inducible transcription factors, we abrogated its global DNA-binding activity using a dominant-negative FOS peptide. We show that FOS and JUN bind to and activate a specific set of vascular genes and that AP-1 inhibition shifts the balance between smooth muscle and hematopoietic differentiation towards blood. Furthermore, AP-1 is required for de novo binding of TEAD4, a transcription factor connected to Hippo signalling. Our bottom-up approach demonstrates that AP-1- and TEAD4-associated cis-regulatory elements form hubs for multiple signalling-responsive transcription factors and define the cistrome that regulates vascular and hematopoietic development by extrinsic signals.

With its significant contribution to cancer mortality globally, advanced colorectal cancer (CRC) requires new treatment strategies. However, despite recent good results for mismatch repair (MMR)-deficient CRC and other malignancies, such as melanoma, the vast majority of MMR-proficient CRCs are resistant to checkpoint inhibitor (CKI) therapy. MMR-proficient CRCs commonly develop from precursor adenomas with enhanced Wnt-signalling due to adenomatous polyposis coli (APC) mutations. In melanomas with enhanced Wnt signalling due to stabilized β-catenin, immune anergy and resistance to CKI therapy has been observed, which is dependent on micro-environmental myelomonocytic (MM) cell depletion in melanoma models. However, MM populations of colorectal adenomas or CRC have not been studied. To characterize resident intestinal MM cell populations during the early stages of tumorigenesis, the present study utilized the ApcMin/+ mouse as a model of MMR-proficient CRC, using enhanced green fluorescent protein (EGFP) expression in the mouse lysozyme (M-lys) lys-EGFP/+ mouse as a pan-myelomonocytic cell marker and a panel of murine macrophage surface markers. Total intestinal lamina propria mononuclear cell (LPMNC) numbers significantly decreased with age (2.32±1.39×107 [n=4] at 33 days of age vs. 1.06±0.24×107 [n=8] at 109 days of age) during intestinal adenoma development in ApcMin/+ mice (P=0.05; unpaired Student's t-test), but not in wild-type littermates (P=0.35). Decreased total LPMNC numbers were associated with atrophy of intestinal lymphoid follicles and the absence of MM/lymphoid cell aggregates in ApcMin/+ mouse intestine, but not spleen, compared with wild-type mice. Furthermore, during the early stage of intestinal adenoma development, there was a two-fold reduction of M-lys expressing cells (P=0.05) and four-fold reduction of ER-HR3 (macrophage sub-set) expressing cells (P=0.05; two tailed Mann-Whitney U test) in mice with reduced total intestinal LPMNCs (n=3). Further studies are necessary to determine the relevance of these findings to immune-surveillance of colorectal adenomas or MMR-proficient CRC CKI therapy resistance.

Acute myeloid leukemia (AML) is caused by recurrent mutations in members of the gene regulatory and signaling machinery that control hematopoietic progenitor cell growth and differentiation. Here, we show that the transcription factor WT1 forms a major node in the rewired mutation-specific gene regulatory networks of multiple AML subtypes. WT1 is frequently either mutated or upregulated in AML, and its expression is predictive for relapse. The WT1 protein exists as multiple isoforms. For two main AML subtypes, we demonstrate that these isoforms exhibit differential patterns of binding and support contrasting biological activities, including enhanced proliferation. We also show that WT1 responds to oncogenic signaling and is part of a signaling-responsive transcription factor hub that controls AML growth. WT1 therefore plays a central and widespread role in AML biology.

Acute myeloid leukemia (AML) is a heterogeneous disease caused by different mutations. Previously, we showed that each mutational subtype develops its specific gene regulatory network (GRN) with transcription factors interacting within multiple gene modules, many of which are transcription factor genes themselves. Here, we hypothesize that highly connected nodes within such networks comprise crucial regulators of AML maintenance. We test this hypothesis using FLT3-ITD-mutated AML as a model and conduct an shRNA drop-out screen informed by this analysis. We show that AML-specific GRNs predict crucial regulatory modules required for AML growth. Furthermore, our work shows that all modules are highly connected and regulate each other. The careful multi-omic analysis of the role of one (RUNX1) module by shRNA and chemical inhibition shows that this transcription factor and its target genes stabilize the GRN of FLT3-ITD+ AML and that its removal leads to GRN collapse and cell death.

Runx proteins are essential for hematopoiesis and play an important role in T cell development by regulating key target genes, such as CD4 and CD8 as well as lymphokine genes, during the specialization of naive CD4 T cells into distinct T helper subsets. In regulatory T (T reg) cells, the signature transcription factor Foxp3 interacts with and modulates the function of several other DNA binding proteins, including Runx family members, at the protein level. We show that Runx proteins also regulate the initiation and the maintenance of Foxp3 gene expression in CD4 T cells. Full-length Runx promoted the de novo expression of Foxp3 during inducible T reg cell differentiation, whereas the isolated dominant-negative Runt DNA binding domain antagonized de novo Foxp3 expression. Foxp3 expression in natural T reg cells remained dependent on Runx proteins and correlated with the binding of Runx/core-binding factor beta to regulatory elements within the Foxp3 locus. Our data show that Runx and Foxp3 are components of a feed-forward loop in which Runx proteins contribute to the expression of Foxp3 and cooperate with Foxp3 proteins to regulate the expression of downstream target genes.

Transcription-factor-induced somatic cell conversions are highly relevant for both basic and clinical research yet their mechanism is not fully understood and it is unclear whether they reflect normal differentiation processes. Here we show that during pre-B-cell-to-macrophage transdifferentiation, C/EBPα binds to two types of myeloid enhancers in B cells: pre-existing enhancers that are bound by PU.1, providing a platform for incoming C/EBPα; and de novo enhancers that are targeted by C/EBPα, acting as a pioneer factor for subsequent binding by PU.1. The order of factor binding dictates the upregulation kinetics of nearby genes. Pre-existing enhancers are broadly active throughout the hematopoietic lineage tree, including B cells. In contrast, de novo enhancers are silent in most cell types except in myeloid cells where they become activated by C/EBP factors. Our data suggest that C/EBPα recapitulates physiological developmental processes by short-circuiting two macrophage enhancer pathways in pre-B cells.

Acute myeloid leukemia (AML) is associated with mutations in transcriptional and epigenetic regulator genes impairing myeloid differentiation. The t(8;21)(q22;q22) translocation generates the RUNX1-ETO fusion protein, which interferes with the hematopoietic master regulator RUNX1. We previously showed that the maintenance of t(8;21) AML is dependent on RUNX1-ETO expression. Its depletion causes extensive changes in transcription factor binding, as well as gene expression, and initiates myeloid differentiation. However, how these processes are connected within a gene regulatory network is unclear. To address this question, we performed Promoter-Capture Hi-C assays, with or without RUNX1-ETO depletion and assigned interacting cis-regulatory elements to their respective genes. To construct a RUNX1-ETO-dependent gene regulatory network maintaining AML, we integrated cis-regulatory element interactions with gene expression and transcription factor binding data. This analysis shows that RUNX1-ETO participates in cis-regulatory element interactions. However, differential interactions following RUNX1-ETO depletion are driven by alterations in the binding of RUNX1-ETO-regulated transcription factors.