Mutational synergy during leukemia induction remodels chromatin accessibility, histone modifications and three-dimensional DNA topology to alter gene expression

Cross-Species Comparative DNA Methylation Reveals Novel Insights into Complex Trait Genetics among Cattle, Sheep, and Goats

The cross-species characterization of evolutionary changes in the functional genome can facilitate the translation of genetic findings across species and the interpretation of the evolutionary basis underlying complex phenotypes. Yet, this has not been fully explored between cattle, sheep, goats, and other mammals. Here, we systematically characterized the evolutionary dynamics of DNA methylation and gene expression in 3 somatic tissues (i.e. brain, liver, and skeletal muscle) and sperm across 7 mammalian species, including 3 ruminant livestock species (cattle, sheep, and goats), humans, pigs, mice, and dogs, by generating and integrating 160 DNA methylation and transcriptomic data sets. We demonstrate dynamic changes of DNA hypomethylated regions and hypermethylated regions in tissue-type manner across cattle, sheep, and goats. Specifically, based on the phylo-epigenetic model of DNA methylome, we identified a total of 25,074 hypomethylated region extension events specific to cattle, which participated in rewiring tissue-specific regulatory network. Furthermore, by integrating genome-wide association studies of 50 cattle traits, we provided novel insights into the genetic and evolutionary basis of complex phenotypes in cattle. Overall, our study provides a valuable resource for exploring the evolutionary dynamics of the functional genome and highlights the importance of cross-species characterization of multiomics data sets for the evolutionary interpretation of complex phenotypes in cattle livestock.

A new genomic framework to categorize pediatric acute myeloid leukemia

Recent studies on pediatric acute myeloid leukemia (pAML) have revealed pediatric-specific driver alterations, many of which are underrepresented in the current classification schemas. To comprehensively define the genomic landscape of pAML, we systematically categorized 887 pAML into 23 mutually distinct molecular categories, including new major entities such as UBTF or BCL11B, covering 91.4% of the cohort. These molecular categories were associated with unique expression profiles and mutational patterns. For instance, molecular categories characterized by specific HOXA or HOXB expression signatures showed distinct mutation patterns of RAS pathway genes, FLT3 or WT1, suggesting shared biological mechanisms. We show that molecular categories were strongly associated with clinical outcomes using two independent cohorts, leading to the establishment of a new prognostic framework for pAML based on these updated molecular categories and minimal residual disease. Together, this comprehensive diagnostic and prognostic framework forms the basis for future classification of pAML and treatment strategies.

Haploinsufficient Transcription Factors in Myeloid Neoplasms

Many transcription factors (TFs) function as tumor suppressor genes with heterozygous phenotypes, yet haploinsufficiency generally has an underappreciated role in neoplasia. This is no less true in myeloid cells, which are normally regulated by a delicately balanced and interconnected transcriptional network. Detailed understanding of TF dose in this circuitry sheds light on the leukemic transcriptome. In this review, we discuss the emerging features of haploinsufficient transcription factors (HITFs). We posit that: (a) monoallelic and biallelic losses can have distinct cellular outcomes; (b) the activity of a TF exists in a greater range than the traditional Mendelian genetic doses; and (c) how a TF is deleted or mutated impacts the cellular phenotype. The net effect of a HITF is a myeloid differentiation block and increased intercellular heterogeneity in the course of myeloid neoplasia.

In vivo screening characterizes chromatin factor functions during normal and malignant hematopoiesis

Cellular differentiation requires extensive alterations in chromatin structure and function, which is elicited by the coordinated action of chromatin and transcription factors. By contrast with transcription factors, the roles of chromatin factors in differentiation have not been systematically characterized. Here, we combine bulk ex vivo and single-cell in vivo CRISPR screens to characterize the role of chromatin factor families in hematopoiesis. We uncover marked lineage specificities for 142 chromatin factors, revealing functional diversity among related chromatin factors (i.e. barrier-to-autointegration factor subcomplexes) as well as shared roles for unrelated repressive complexes that restrain excessive myeloid differentiation. Using epigenetic profiling, we identify functional interactions between lineage-determining transcription factors and several chromatin factors that explain their lineage dependencies. Studying chromatin factor functions in leukemia, we show that leukemia cells engage homeostatic chromatin factor functions to block differentiation, generating specific chromatin factor-transcription factor interactions that might be therapeutically targeted. Together, our work elucidates the lineage-determining properties of chromatin factors across normal and malignant hematopoiesis.

The Molecular Context of Oxidant Stress Response in Cancer Establishes ALDH1A1 as a Critical Target: What This Means for Acute Myeloid Leukemia

The protein family of aldehyde dehydrogenases (ALDH) encompasses nineteen members. The ALDH1 subfamily consists of enzymes with similar activity, having the capacity to neutralize lipid peroxidation products and to generate retinoic acid; however, only ALDH1A1 emerges as a significant risk factor in acute myeloid leukemia. Not only is the gene ALDH1A1 on average significantly overexpressed in the poor prognosis group at the RNA level, but its protein product, ALDH1A1 protects acute myeloid leukemia cells from lipid peroxidation byproducts. This capacity to protect cells can be ascribed to the stability of the enzyme under conditions of oxidant stress. The capacity to protect cells is evident both in vitro, as well as in mouse xenografts of those cells, shielding cells effectively from a number of potent antineoplastic agents. However, the role of ALDH1A1 in acute myeloid leukemia has been unclear in the past due to evidence that normal cells often have higher aldehyde dehydrogenase activity than leukemic cells. This being true, ALDH1A1 RNA expression is significantly associated with poor prognosis. It is hence imperative that ALDH1A1 is methodically targeted, particularly for the acute myeloid leukemia patients of the poor prognosis risk group that overexpress ALDH1A1 RNA.

Loop stacking organizes genome folding from TADs to chromosomes

Although population-level analyses revealed significant roles for CTCF and cohesin in mammalian genome organization, their contributions at the single-cell level remain incompletely understood. Here, we used a super-resolution microscopy approach to measure the effects of removal of CTCF or cohesin in mouse embryonic stem cells. Single-chromosome traces revealed cohesin-dependent loops, frequently stacked at their loop anchors forming multi-way contacts (hubs), bridging across TAD boundaries. Despite these bridging interactions, chromatin in intervening TADs was not intermixed, remaining separated in distinct loops around the hub. At the multi-TAD scale, steric effects from loop stacking insulated local chromatin from ultra-long range (>4 Mb) contacts. Upon cohesin removal, the chromosomes were more disordered and increased cell-cell variability in gene expression. Our data revise the TAD-centric understanding of CTCF and cohesin and provide a multi-scale, structural picture of how they organize the genome on the single-cell level through distinct contributions to loop stacking.

FLT3ITD drives context-specific changes in cell identity and variable interferon dependence during AML initiation

Acute myeloid leukemia (AML) initiation requires multiple rate-limiting mutations to cooperatively reprogram progenitor cell identity. For example, FLT3 internal tandem duplication (FLT3ITD) mutations cooperate with a variety of different initiating mutations to reprogram myeloid progenitor fate. These initiating mutations often skew toward either pediatric or adult AML patient populations, though FLT3ITD itself occurs at similar frequencies in both age groups. This raises the question of whether FLT3ITD might induce distinct transcriptional programs and unmask distinct therapeutic vulnerabilities when paired with pediatric, as opposed to adult AML-initiating mutations. To explore this possibility, we compared AML evolution in mice that carried Flt3ITD/NUP98-HOXD13 (NHD13) or Flt3ITD/Runx1DEL mutation pairs, which are respectively most common in pediatric and adult AML. Single-cell analyses and epigenome profiling revealed distinct interactions between Flt3ITD and its cooperating mutations. Whereas Flt3ITD and Flt3ITD/Runx1DEL caused aberrant expansion of myeloid progenitors, Flt3ITD/NHD13 drove the emergence of a pre-AML population that did not resemble normal hematopoietic progenitors. Differences between Flt3ITD/Runx1DEL and Flt3ITD/NHD13 cooperative target gene expression extended to fully transformed AML as well. Flt3ITD/NHD13 cooperative target genes were enriched in human NUP98-translocated AML. Flt3ITD/NHD13 selectively hijacked type I interferon signaling to drive expansion of the pre-AML population. Blocking interferon signaling delayed AML initiation and extended survival. Thus, common AML driver mutations, such as FLT3ITD, can coopt different mechanisms of transformation in different genetic contexts. Furthermore, pediatric-biased NUP98 fusions convey actionable interferon dependence.

Mutant NPM1 Hijacks Transcriptional Hubs to Maintain Pathogenic Gene Programs in Acute Myeloid Leukemia

Nucleophosmin (NPM1) is a ubiquitously expressed nucleolar protein with a wide range of biological functions. In 30% of acute myeloid leukemia (AML), the terminal exon of NPM1 is often found mutated, resulting in the addition of a nuclear export signal and a shift of the protein to the cytoplasm (NPM1c). AMLs carrying this mutation have aberrant expression of the HOXA/B genes, whose overexpression leads to leukemogenic transformation. Here, for the first time, we comprehensively prove that NPM1c binds to a subset of active gene promoters in NPM1c AMLs, including well-known leukemia-driving genes-HOXA/B cluster genes and MEIS1. NPM1c sustains the active transcription of key target genes by orchestrating a transcription hub and maintains the active chromatin landscape by inhibiting the activity of histone deacetylases. Together, these findings reveal the neomorphic function of NPM1c as a transcriptional amplifier for leukemic gene expression and open up new paradigms for therapeutic intervention.

SIGNIFICANCE

NPM1 mutation is the most common mutation in AML, yet the mechanism of how the mutant protein results in AML remains unclear. Here, for the first time, we prove mutant NPM1 directly binds to active chromatin regions and hijacks the transcription of AML-driving genes. See related article by Uckelmann et al., p. 746. This article is highlighted in the In This Issue feature, p. 517.

Multiomics data integration to reveal chromatin remodeling and reorganization induced by gene mutational synergy

Recurrent gene mutations often cooperate in a predefined stepwise and synergistic manner to alter global transcription, through directly or indirectly remodeling epigenetic landscape on linear and three-dimensional (3D) scales. Here, we present a multiomics data integration approach to investigate the impact of gene mutational synergy on transcription, chromatin states, and 3D chromatin organization in a murine leukemia model. This protocol provides an executable framework to study epigenetic remodeling induced by cooperating gene mutations and to identify the critical regulatory network involved. For complete details on the use and execution of this protocol, please refer to Yun et al. (2021).

Nuclear factor Nrf2 promotes glycosidase OGG1 expression by activating the AKT pathway to enhance leukemia cell resistance to cytarabine

Chemotherapy resistance is the dominant challenge in the treatment of acute myeloid leukemia (AML). Nuclear factor E2-related factor 2 (Nrf2) exerts a vital function in drug resistance of many tumors. Nevertheless, the potential molecular mechanism of Nrf2 regulating the base excision repair pathway that mediates AML chemotherapy resistance remains unclear. Here, in clinical samples, we found that the high expression of Nrf2 and base excision repair pathway gene encoding 8-hydroxyguanine DNA glycosidase (OGG1) was associated with AML disease progression. In vitro, Nrf2 and OGG1 were highly expressed in drug-resistant leukemia cells. Upregulation of Nrf2 in leukemia cells by lentivirus transfection could decrease the sensitivity of leukemia cells to cytarabine, whereas downregulation of Nrf2 in drug-resistant cells could enhance leukemia cell chemosensitivity. Meanwhile, we found that Nrf2 could positively regulate OGG1 expression in leukemia cells. Our chromatin immunoprecipitation assay revealed that Nrf2 could bind to the promoter of OGG1. Furthermore, the use of OGG1 inhibitor TH5487 could partially reverse the inhibitory effect of upregulated Nrf2 on leukemia cell apoptosis. In vivo, downregulation of Nrf2 could increase the sensitivity of leukemia cell to cytarabine and decrease OGG1 expression. Mechanistically, Nrf2-OGG1 axis-mediated AML resistance might be achieved by activating the AKT signaling pathway to regulate downstream apoptotic proteins. Thus, this study reveals a novel mechanism of Nrf2-promoting drug resistance in leukemia, which may provide a potential therapeutic target for the treatment of drug-resistant/refractory leukemia.

A novel network pharmacology approach for leukaemia differentiation therapy using Mogrify®

Acute myeloid leukaemia (AML) is a rapidly fatal blood cancer that is characterised by the accumulation of immature myeloid cells in the blood and bone marrow as a result of blocked differentiation. Methods which identify master transcriptional regulators of AML subtype-specific leukaemia cell states and their combinations could be critical for discovering novel differentiation-inducing therapies. In this proof-of-concept study, we demonstrate a novel utility of the Mogrify® algorithm in identifying combinations of transcription factors (TFs) and drugs, which recapitulate granulocytic differentiation of the NB4 acute promyelocytic leukaemia (APL) cell line, using two different approaches. In the first approach, Connectivity Map (CMAP) analysis of these TFs and their target networks outperformed standard approaches, retrieving ATRA as the top hit. We identify dimaprit and mebendazole as a drug combination which induces myeloid differentiation. In the second approach, we show that genetic manipulation of specific Mogrify®-identified TFs (MYC and IRF1) leads to co-operative induction of APL differentiation, as does pharmacological targeting of these TFs using currently available compounds. We also show that loss of IRF1 blunts ATRA-mediated differentiation, and that MYC represses IRF1 expression through recruitment of PML-RARα, the driver fusion oncoprotein in APL, to the IRF1 promoter. Finally, we demonstrate that these drug combinations can also induce differentiation of primary patient-derived APL cells, and highlight the potential of targeting MYC and IRF1 in high-risk APL. Thus, these results suggest that Mogrify® could be used for drug discovery or repositioning in leukaemia differentiation therapy for other subtypes of leukaemia or cancers.

Dissection of the MEF2D-IRF8 transcriptional circuit dependency in acute myeloid leukemia

Transcriptional dysregulation is a prominent feature in leukemia. Here, we systematically surveyed transcription factor (TF) vulnerabilities in leukemia and uncovered TF clusters that exhibit context-specific vulnerabilities within and between different subtypes of leukemia. Among these TF clusters, we demonstrated that acute myeloid leukemia (AML) with high IRF8 expression was addicted to MEF2D. MEF2D and IRF8 form an autoregulatory loop via direct binding to mutual enhancer elements. One important function of this circuit in AML is to sustain PU.1/MEIS1 co-regulated transcriptional outputs via stabilizing PU.1’s chromatin occupancy. We illustrated that AML could acquire dependency on this circuit through various oncogenic mechanisms that results in the activation of their enhancers. In addition to forming a circuit, MEF2D and IRF8 can also separately regulate gene expression, and dual perturbation of these two TFs leads to a more robust inhibition of AML proliferation. Collectively, our results revealed a TF circuit essential for AML survival.

•

MEF2D is a context-specific vulnerability in IRF8^hi AML

•

MEF2D and IRF8 form a transcriptional circuit via binding to each other’s enhancers

•

MEF2D-IRF8 circuit supports PU.1’s chromatin occupancy and transcriptional output

•

MEF2D and IRF8 can regulate separate gene expression programs alongside the circuit

PU.1-Dependent Enhancer Inhibition Separates Tet2-Deficient Hematopoiesis from Malignant Transformation

Cytosine hypermethylation in and around DNA-binding sites of master transcription factors, including PU.1, occurs in aging hematopoietic stem cells following acquired loss-of-function mutations of DNA methyl-cytosine dioxygenase ten-eleven translocation-2 (TET2), albeit functional relevance has been unclear. We show that Tet2-deficient mouse hematopoietic stem and progenitor cells undergo malignant transformation upon compromised gene regulation through heterozygous deletion of an upstream regulatory region (UREΔ/WT) of the PU.1 gene. Although compatible with multilineage blood formation at young age, Tet2-deficient PU.1 UREΔ/WT mice develop highly penetrant, transplantable acute myeloid leukemia (AML) during aging. Leukemic stem and progenitor cells show hypermethylation at putative PU.1-binding sites, fail to activate myeloid enhancers, and are hallmarked by a signature of genes with impaired expression shared with human AML. Our study demonstrates that Tet2 and PU.1 jointly suppress leukemogenesis and uncovers a methylation-sensitive PU.1-dependent gene network as a unifying molecular vulnerability associated with AML.

SIGNIFICANCE

We identify moderately impaired PU.1 mRNA expression as a biological modality predisposing Tet2-deficient hematopoietic stem and progenitor cells to malignant transformation. Our study furthermore uncovers a methylation-sensitive PU.1 gene network as a common feature of myeloid leukemia potentially allowing for the identification of patients at risk for malignant transformation. See related commentary by Schleicher and Pietras, p. 378. This article is highlighted in the In This Issue feature, p. 369.

Dissecting the Genetic and Non-Genetic Heterogeneity of Acute Myeloid Leukemia Using Next-Generation Sequencing and In Vivo Models

Simple Summary

Acute myeloid leukemia (AML) is an extremely aggressive form of blood cancer with high rates of treatment failure. AML arises from the stepwise acquisition of genetic aberrations and is a highly heterogeneous disorder. Recent research has shown that individual AML samples often contain several clones that are defined by a distinct combination of genetic lesions, epigenetic patterns and cell surface marker expression profiles. A better understanding of the clonal dynamics of AML is required to develop novel treatment strategies against this disease. In this review, we discuss the recent developments that have further deepened our understanding of clonal evolution and heterogeneity in AML.

Abstract

Acute myeloid leukemia (AML) is an extremely aggressive and heterogeneous disorder that results from the transformation of hematopoietic stem cells. Although our understanding of the molecular pathology of AML has greatly improved in the last few decades, the overall and relapse free survival rates among AML patients remain quite poor. This is largely due to evolution of the disease and selection of the fittest, treatment-resistant leukemic clones. There is increasing evidence that most AMLs possess a highly complex clonal architecture and individual leukemias are comprised of genetically, phenotypically and epigenetically distinct clones, which are continually evolving. Advances in sequencing technologies as well as studies using murine AML models have provided further insights into the heterogeneity of leukemias. We will review recent advances in the field of genetic and non-genetic heterogeneity in AML.

Cell Origin-Dependent Cooperativity of Mutant Dnmt3a and Npm1 in Clonal Hematopoiesis and Myeloid Malignancy

In adult acute myeloid leukemia (AML), acquisition of driver somatic mutations may be preceded by a benign state termed clonal hematopoiesis (CH). To develop therapeutic strategies to prevent leukemia development from CH, it is important to understand the mechanisms by which CH-driving and AML-driving mutations cooperate. Here, we use mice with inducible mutant alleles common in human CH (DNMT3AR882; mouse Dnmt3aR878H) and AML (NPM1c; mouse Npm1cA). We find that Dnmt3aR878H/+ hematopoietic stem cells (HSCs), but not multipotent progenitor cell (MPP) subsets, have reduced expression of cytokine and pro-inflammatory transcriptional signatures and a functional competitive advantage over their wild-type counterparts. Dnmt3aR878H/+ HSCs are the most potent cell type transformed by Npm1cA, generating myeloid malignancies in which few additional cooperating somatic mutation events were detected. At a molecular level, Npm1cA in cooperation with Dnmt3aR878H acutely increased accessibility of a distinct set of promoters in HSCs compared to MPP cells. These promoters were enriched for cell cycling, PI3K/AKT/mTOR signaling, stem cell signatures, and targets of transcription factors including NFAT and the chromatin binding factor HMGB1, which have been implicated in human AML. These results demonstrate cooperativity between pre-existing Dnmt3aR878H and Npm1cA at the chromatin level, where specific loci altered in accessibility by Npm1cA are dependent on cell context as well as Dnmt3a mutation status. These findings have implications for biological understanding and therapeutic intervention into transformation from CH to AML.

Murine Models of Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is a phenotypically and genetically heterogeneous hematologic malignancy. Extensive sequencing efforts have mapped the genomic landscape of adult and pediatric AML revealing a number of biologically and prognostically relevant driver lesions. Beyond identifying recurrent genetic aberrations, it is of critical importance to fully delineate the complex mechanisms by which they contribute to the initiation and evolution of disease to ultimately facilitate the development of targeted therapies. Towards these aims, murine models of AML are indispensable research tools. The rapid evolution of genetic engineering techniques over the past 20 years has greatly advanced the use of murine models to mirror specific genetic subtypes of human AML, define cell-intrinsic and extrinsic disease mechanisms, study the interaction between co-occurring genetic lesions, and test novel therapeutic approaches. This review summarizes the mouse model systems that have been developed to recapitulate the most common genomic subtypes of AML. We will discuss the strengths and weaknesses of varying modeling strategies, highlight major discoveries emanating from these model systems, and outline future opportunities to leverage emerging technologies for mechanistic and preclinical investigations.