PLCG1 is required for AML1-ETO leukemia stem cell self-renewal.

In an effort to identify novel drugs targeting fusion-oncogene-induced acute myeloid leukemia (AML), we performed high-resolution proteomic analysis. In AML1-ETO (AE)-driven AML, we uncovered a deregulation of phospholipase C (PLC) signaling. We identified PLCgamma 1 (PLCG1) as a specific target of the AE fusion protein that is induced after AE binding to intergenic regulatory DNA elements. Genetic inactivation of PLCG1 in murine and human AML inhibited AML1-ETO dependent self-renewal programs, leukemic proliferation, and leukemia maintenance in vivo. In contrast, PLCG1 was dispensable for normal hematopoietic stem and progenitor cell function. These findings are extended to and confirmed by pharmacologic perturbation of Ca++-signaling in AML1-ETO AML cells, indicating that the PLCG1 pathway poses an important therapeutic target for AML1-ETO+ leukemic stem cells.

Epigenetic regulator genes direct lineage switching in MLL/AF4 leukemia.

The fusion gene MLL/AF4 defines a high-risk subtype of pro-B acute lymphoblastic leukemia. Relapse can be associated with a lineage switch from acute lymphoblastic to acute myeloid leukemia, resulting in poor clinical outcomes caused by resistance to chemotherapies and immunotherapies. In this study, the myeloid relapses shared oncogene fusion breakpoints with their matched lymphoid presentations and originated from various differentiation stages from immature progenitors through to committed B-cell precursors. Lineage switching is linked to substantial changes in chromatin accessibility and rewiring of transcriptional programs, including alternative splicing. These findings indicate that the execution and maintenance of lymphoid lineage differentiation is impaired. The relapsed myeloid phenotype is recurrently associated with the altered expression, splicing, or mutation of chromatin modifiers, including CHD4 coding for the ATPase/helicase of the nucleosome remodelling and deacetylation complex. Perturbation of CHD4 alone or in combination with other mutated epigenetic modifiers induces myeloid gene expression in MLL/AF4+ cell models, indicating that lineage switching in MLL/AF4 leukemia is driven and maintained by disrupted epigenetic regulation.

Identification of gene specific cis-regulatory elements during differentiation of mouse embryonic stem cells: An integrative approach using high-throughput datasets.

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.

Cooperative binding of AP-1 and TEAD4 modulates the balance between vascular smooth muscle and hemogenic cell fate.

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.

A FOXO1-induced oncogenic network defines the AML1-ETO preleukemic program.

Understanding and blocking the self-renewal pathway of preleukemia stem cells could prevent acute myeloid leukemia (AML) relapse. In this study, we show that increased FOXO1 represents a critical mechanism driving aberrant self-renewal in preleukemic cells expressing the t(8;21)-associated oncogene AML1-ETO (AE). Although generally considered as a tumor suppressor, FOXO1 is consistently upregulated in t(8;21) AML. Expression of FOXO1 in human CD34+ cells promotes a preleukemic state with enhanced self-renewal and dysregulated differentiation. The DNA binding domain of FOXO1 is essential for these functions. FOXO1 activates a stem cell molecular signature that is also present in AE preleukemia cells and preserved in t(8;21) patient samples. Genome-wide binding studies show that AE and FOXO1 share the majority of their binding sites, whereby FOXO1 binds to multiple crucial self-renewal genes and is required for their activation. In agreement with this observation, genetic and pharmacological ablation of FOXO1 inhibited the long-term proliferation and clonogenicity of AE cells and t(8;21) AML cell lines. Targeting of FOXO1 therefore provides a potential therapeutic strategy for elimination of stem cells at both preleukemic and leukemic stages.

LMO2 is required for TAL1 DNA binding activity and initiation of definitive haematopoiesis at the haemangioblast stage.

LMO2 is a bridging factor within a DNA binding complex and is required for definitive haematopoiesis to occur. The developmental stage of the block in haematopoietic specification is not known. We show that Lmo2-/- mouse embryonic stem cells differentiated to Flk-1+ haemangioblasts, but less efficiently to haemogenic endothelium, which only produced primitive haematopoietic progenitors. Genome-wide approaches indicated that LMO2 is required at the haemangioblast stage to position the TAL1/LMO2/LDB1 complex to regulatory elements that are important for the establishment of the haematopoietic developmental program. In the absence of LMO2, the target site recognition of TAL1 is impaired. The lack of LMO2 resulted in altered gene expression levels already at the haemangioblast stage, with transcription factor genes accounting for ∼15% of affected genes. Comparison of Lmo2-/- with Tal1-/- Flk-1+ cells further showed that TAL1 was required to initiate or sustain Lmo2 expression.

A crucial role for the ubiquitously expressed transcription factor Sp1 at early stages of hematopoietic specification.

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.

RUNX1 reshapes the epigenetic landscape at the onset of haematopoiesis.

Cell fate decisions during haematopoiesis are governed by lineage-specific transcription factors, such as RUNX1, SCL/TAL1, FLI1 and C/EBP family members. To gain insight into how these transcription factors regulate the activation of haematopoietic genes during embryonic development, we measured the genome-wide dynamics of transcription factor assembly on their target genes during the RUNX1-dependent transition from haemogenic endothelium (HE) to haematopoietic progenitors. Using a Runx1-/- embryonic stem cell differentiation model expressing an inducible Runx1 gene, we show that in the absence of RUNX1, haematopoietic genes bind SCL/TAL1, FLI1 and C/EBPß and that this early priming is required for correct temporal expression of the myeloid master regulator PU.1 and its downstream targets. After induction, RUNX1 binds to numerous de novo sites, initiating a local increase in histone acetylation and rapid global alterations in the binding patterns of SCL/TAL1 and FLI1. The acquisition of haematopoietic fate controlled by Runx1 therefore does not represent the establishment of a new regulatory layer on top of a pre-existing HE program but instead entails global reorganization of lineage-specific transcription factor assemblies.

Wellington-bootstrap: differential DNase-seq footprinting identifies cell-type determining transcription factors.

BACKGROUND: The analysis of differential gene expression is a fundamental tool to relate gene regulation with specific biological processes. Differential binding of transcription factors (TFs) can drive differential gene expression. While DNase-seq data can provide global snapshots of TF binding, tools for detecting differential binding from pairs of DNase-seq data sets are lacking. RESULTS: In order to link expression changes with changes in TF binding we introduce the concept of differential footprinting alongside a computational tool. We demonstrate that differential footprinting is associated with differential gene expression and can be used to define cell types by their specific TF occupancy patterns. CONCLUSIONS: Our new tool, Wellington-bootstrap, will enable the detection of differential TF binding facilitating the study of gene regulatory systems.

Leukemic stem cells activate lineage inappropriate signalling pathways to promote their growth.

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.

Gene regulatory network analysis predicts cooperating transcription factor regulons required for FLT3-ITD+ AML growth.

AML is a heterogenous disease caused by different mutations. We have previously shown that each mutational sub-type develops its specific gene regulatory network (GRN) with transcription factors interacting with multiple gene modules, many of which are transcription factor genes themselves. Here we hypothesized that highly connected nodes within such networks comprise crucial regulators of AML maintenance. We tested this hypothesis using FLT3-ITD mutated AML as a model and conducted an shRNA drop-out screen informed by this analysis. We show that AML-specific GRNs predict identifying crucial regulatory modules required for AML but not normal cellular 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.

Identification and interrogation of the gene regulatory network of CEBPA-double mutant acute myeloid leukemia.

Acute myeloid leukemia (AML) is a heterogeneous hematological malignancy caused by mutations in genes encoding transcriptional and epigenetic regulators together with signaling genes. It is characterized by a disturbance of differentiation and abnormal proliferation of hematopoietic progenitors. We have previously shown that each AML subtype establishes its own core gene regulatory network (GRN), consisting of transcription factors binding to their target genes and imposing a specific gene expression pattern that is required for AML maintenance. In this study, we integrate gene expression, open chromatin and ChIP data with promoter-capture Hi-C data to define a refined core GRN common to all patients with CEBPA-double mutant (CEBPAN/C) AML. These mutations disrupt the structure of a major regulator of myelopoiesis. We identify the binding sites of mutated C/EBPα proteins in primary cells, we show that C/EBPα, AP-1 factors and RUNX1 colocalize and are required for AML maintenance, and we employ single cell experiments to link important network nodes to the specific differentiation trajectory from leukemic stem to blast cells. Taken together, our study provides an important resource which predicts the specific therapeutic vulnerabilities of this AML subtype in human cells.

Gene regulatory network analysis predicts cooperating transcription factor regulons required for FLT3-ITD+ AML growth.

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.

Transcriptional reprogramming by mutated IRF4 in lymphoma.

Disease-causing mutations in genes encoding transcription factors (TFs) can affect TF interactions with their cognate DNA-binding motifs. Whether and how TF mutations impact upon the binding to TF composite elements (CE) and the interaction with other TFs is unclear. Here, we report a distinct mechanism of TF alteration in human lymphomas with perturbed B cell identity, in particular classic Hodgkin lymphoma. It is caused by a recurrent somatic missense mutation c.295 T > C (p.Cys99Arg; p.C99R) targeting the center of the DNA-binding domain of Interferon Regulatory Factor 4 (IRF4), a key TF in immune cells. IRF4-C99R fundamentally alters IRF4 DNA-binding, with loss-of-binding to canonical IRF motifs and neomorphic gain-of-binding to canonical and non-canonical IRF CEs. IRF4-C99R thoroughly modifies IRF4 function by blocking IRF4-dependent plasma cell induction, and up-regulates disease-specific genes in a non-canonical Activator Protein-1 (AP-1)-IRF-CE (AICE)-dependent manner. Our data explain how a single mutation causes a complex switch of TF specificity and gene regulation and open the perspective to specifically block the neomorphic DNA-binding activities of a mutant TF.

PCRPi: Presaging Critical Residues in Protein interfaces, a new computational tool to chart hot spots in protein interfaces.

Protein-protein interactions (PPIs) are ubiquitous in Biology, and thus offer an enormous potential for the discovery of novel therapeutics. Although protein interfaces are large and lack defining physiochemical traits, is well established that only a small portion of interface residues, the so-called hot spot residues, contribute the most to the binding energy of the protein complex. Moreover, recent successes in development of novel drugs aimed at disrupting PPIs rely on targeting such residues. Experimental methods for describing critical residues are lengthy and costly; therefore, there is a need for computational tools that can complement experimental efforts. Here, we describe a new computational approach to predict hot spot residues in protein interfaces. The method, called Presaging Critical Residues in Protein interfaces (PCRPi), depends on the integration of diverse metrics into a unique probabilistic measure by using Bayesian Networks. We have benchmarked our method using a large set of experimentally verified hot spot residues and on a blind prediction on the protein complex formed by HRAS protein and a single domain antibody. Under both scenarios, PCRPi delivered consistent and accurate predictions. Finally, PCRPi is able to handle cases where some of the input data is either missing or not reliable (e.g. evolutionary information).

Isoform-specific and signaling-dependent propagation of acute myeloid leukemia by Wilms tumor 1.

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.

RUNX1/RUNX1T1 mediates alternative splicing and reorganises the transcriptional landscape in leukemia.

The fusion oncogene RUNX1/RUNX1T1 encodes an aberrant transcription factor, which plays a key role in the initiation and maintenance of acute myeloid leukemia. Here we show that the RUNX1/RUNX1T1 oncogene is a regulator of alternative RNA splicing in leukemic cells. The comprehensive analysis of RUNX1/RUNX1T1-associated splicing events identifies two principal mechanisms that underlie the differential production of RNA isoforms: (i) RUNX1/RUNX1T1-mediated regulation of alternative transcription start site selection, and (ii) direct or indirect control of the expression of genes encoding splicing factors. The first mechanism leads to the expression of RNA isoforms with alternative structure of the 5'-UTR regions. The second mechanism generates alternative transcripts with new junctions between internal cassettes and constitutive exons. We also show that RUNX1/RUNX1T1-mediated differential splicing affects several functional groups of genes and produces proteins with unique conserved domain structures. In summary, this study reveals alternative splicing as an important component of transcriptome re-organization in leukemia by an aberrant transcriptional regulator.

RUNX1/ETO and mutant KIT both contribute to programming the transcriptional and chromatin landscape in t(8;21) acute myeloid leukemia.

Acute myeloid leukemia development occurs in a stepwise fashion whereby an original driver mutation is followed by additional mutations. The first type of mutations tends to be in genes encoding members of the epigenetic/transcription regulatory machinery (i.e., RUNX1, DNMT3A, TET2), while the secondary mutations often involve genes encoding members of signaling pathways that cause uncontrolled growth of such cells such as the growth factor receptors c-KIT of FLT3. Patients usually present with both types of mutations, but it is currently unclear how both mutational events shape the epigenome in developing AML cells. To this end we generated an in vitro model of t(8;21) AML by expressing its driver oncoprotein RUNX1-ETO with or without a mutated (N822K) KIT protein. Expression of N822K-c-KIT strongly increases the self-renewal capacity of RUNX1-ETO-expressing cells. Global analysis of gene expression changes and alterations in the epigenome revealed that N822K-c-KIT expression profoundly influences the open chromatin landscape and transcription factor binding. However, our experiments also revealed that double mutant cells still differ from their patient-derived counterparts, highlighting the importance of studying patient cells to obtain a true picture of how gene regulatory networks have been reprogrammed during tumorigenesis.

Rewiring of the Transcription Factor Network in Acute Myeloid Leukemia.

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.

Global long terminal repeat activation participates in establishing the unique gene expression programme of classical Hodgkin lymphoma.

Long terminal repeat (LTR) elements are wide-spread in the human genome and have the potential to act as promoters and enhancers. Their expression is therefore under tight epigenetic control. We previously reported in classical Hodgkin Lymphoma (cHL) that a member of the THE1B class of LTR elements acted as a promoter for the proto-oncogene and growth factor receptor gene CSF1R and that expression of this gene is required for cHL tumour survival. However, to which extent and how such elements participate in globally shaping the unique cHL gene expression programme is unknown. To address this question we mapped the genome-wide activation of THE1-LTRs in cHL cells using a targeted next generation sequencing approach (RACE-Seq). Integration of these data with global gene expression data from cHL and control B cell lines showed a unique pattern of LTR activation impacting on gene expression, including genes associated with the cHL phenotype. We also show that global LTR activation is induced by strong inflammatory stimuli. Together these results demonstrate that LTR activation provides an additional layer of gene deregulation in classical Hodgkin lymphoma and highlight the potential impact of genome-wide LTR activation in other inflammatory diseases.