HBO1 is required for the maintenance of leukaemia stem cells.

Acute myeloid leukaemia (AML) is a heterogeneous disease characterized by transcriptional dysregulation that results in a block in differentiation and increased malignant self-renewal. Various epigenetic therapies aimed at reversing these hallmarks of AML have progressed into clinical trials, but most show only modest efficacy owing to an inability to effectively eradicate leukaemia stem cells (LSCs)1. Here, to specifically identify novel dependencies in LSCs, we screened a bespoke library of small hairpin RNAs that target chromatin regulators in a unique ex vivo mouse model of LSCs. We identify the MYST acetyltransferase HBO1 (also known as KAT7 or MYST2) and several known members of the HBO1 protein complex as critical regulators of LSC maintenance. Using CRISPR domain screening and quantitative mass spectrometry, we identified the histone acetyltransferase domain of HBO1 as being essential in the acetylation of histone H3 at K14. H3 acetylated at K14 (H3K14ac) facilitates the processivity of RNA polymerase II to maintain the high expression of key genes (including Hoxa9 and Hoxa10) that help to sustain the functional properties of LSCs. To leverage this dependency therapeutically, we developed a highly potent small-molecule inhibitor of HBO1 and demonstrate its mode of activity as a competitive analogue of acetyl-CoA. Inhibition of HBO1 phenocopied our genetic data and showed efficacy in a broad range of human cell lines and primary AML cells from patients. These biological, structural and chemical insights into a therapeutic target in AML will enable the clinical translation of these findings.

Evolutionary routes and KRAS dosage define pancreatic cancer phenotypes.

The poor correlation of mutational landscapes with phenotypes limits our understanding of the pathogenesis and metastasis of pancreatic ductal adenocarcinoma (PDAC). Here we show that oncogenic dosage-variation has a critical role in PDAC biology and phenotypic diversification. We find an increase in gene dosage of mutant KRAS in human PDAC precursors, which drives both early tumorigenesis and metastasis and thus rationalizes early PDAC dissemination. To overcome the limitations posed to gene dosage studies by the stromal richness of PDAC, we have developed large cell culture resources of metastatic mouse PDAC. Integration of cell culture genomes, transcriptomes and tumour phenotypes with functional studies and human data reveals additional widespread effects of oncogenic dosage variation on cell morphology and plasticity, histopathology and clinical outcome, with the highest KrasMUT levels underlying aggressive undifferentiated phenotypes. We also identify alternative oncogenic gains (Myc, Yap1 or Nfkb2), which collaborate with heterozygous KrasMUT in driving tumorigenesis, but have lower metastatic potential. Mechanistically, different oncogenic gains and dosages evolve along distinct evolutionary routes, licensed by defined allelic states and/or combinations of hallmark tumour suppressor alterations (Cdkn2a, Trp53, Tgfß-pathway). Thus, evolutionary constraints and contingencies direct oncogenic dosage gain and variation along defined routes to drive the early progression of PDAC and shape its downstream biology. Our study uncovers universal principles of Ras-driven oncogenesis that have potential relevance beyond pancreatic cancer.

Mutant calreticulin knockin mice develop thrombocytosis and myelofibrosis without a stem cell self-renewal advantage.

Somatic mutations in the endoplasmic reticulum chaperone calreticulin (CALR) are detected in approximately 40% of patients with essential thrombocythemia (ET) and primary myelofibrosis (PMF). Multiple different mutations have been reported, but all result in a +1-bp frameshift and generate a novel protein C terminus. In this study, we generated a conditional mouse knockin model of the most common CALR mutation, a 52-bp deletion. The mutant novel human C-terminal sequence is integrated into the otherwise intact mouse CALR gene and results in mutant CALR expression under the control of the endogenous mouse locus. CALRdel/+ mice develop a transplantable ET-like disease with marked thrombocytosis, which is associated with increased and morphologically abnormal megakaryocytes and increased numbers of phenotypically defined hematopoietic stem cells (HSCs). Homozygous CALRdel/del mice developed extreme thrombocytosis accompanied by features of MF, including leukocytosis, reduced hematocrit, splenomegaly, and increased bone marrow reticulin. CALRdel/+ HSCs were more proliferative in vitro, but neither CALRdel/+ nor CALRdel/del displayed a competitive transplantation advantage in primary or secondary recipient mice. These results demonstrate the consequences of heterozygous and homozygous CALR mutations and provide a powerful model for dissecting the pathogenesis of CALR-mutant ET and PMF.

Molecular synergy underlies the co-occurrence patterns and phenotype of NPM1-mutant acute myeloid leukemia.

NPM1 mutations define the commonest subgroup of acute myeloid leukemia (AML) and frequently co-occur with FLT3 internal tandem duplications (ITD) or, less commonly, NRAS or KRAS mutations. Co-occurrence of mutant NPM1 with FLT3-ITD carries a significantly worse prognosis than NPM1-RAS combinations. To understand the molecular basis of these observations, we compare the effects of the 2 combinations on hematopoiesis and leukemogenesis in knock-in mice. Early effects of these mutations on hematopoiesis show that compound Npm1cA/+;NrasG12D/+ or Npm1cA;Flt3ITD share a number of features: Hox gene overexpression, enhanced self-renewal, expansion of hematopoietic progenitors, and myeloid differentiation bias. However, Npm1cA;Flt3ITD mutants displayed significantly higher peripheral leukocyte counts, early depletion of common lymphoid progenitors, and a monocytic bias in comparison with the granulocytic bias in Npm1cA/+;NrasG12D/+ mutants. Underlying this was a striking molecular synergy manifested as a dramatically altered gene expression profile in Npm1cA;Flt3ITD , but not Npm1cA/+;NrasG12D/+ , progenitors compared with wild-type. Both double-mutant models developed high-penetrance AML, although latency was significantly longer with Npm1cA/+;NrasG12D/+ During AML evolution, both models acquired additional copies of the mutant Flt3 or Nras alleles, but only Npm1cA/+;NrasG12D/+ mice showed acquisition of other human AML mutations, including IDH1 R132Q. We also find, using primary Cas9-expressing AMLs, that Hoxa genes and selected interactors or downstream targets are required for survival of both types of double-mutant AML. Our results show that molecular complementarity underlies the higher frequency and significantly worse prognosis associated with NPM1c/FLT3-ITD vs NPM1/NRAS-G12D-mutant AML and functionally confirm the role of HOXA genes in NPM1c-driven AML.

Histone deacetylase (HDAC) 1 and 2 are essential for accurate cell division and the pluripotency of embryonic stem cells.

Histone deacetylases 1 and 2 (HDAC1/2) form the core catalytic components of corepressor complexes that modulate gene expression. In most cell types, deletion of both Hdac1 and Hdac2 is required to generate a discernible phenotype, suggesting their activity is largely redundant. We have therefore generated an ES cell line in which Hdac1 and Hdac2 can be inactivated simultaneously. Loss of HDAC1/2 resulted in a 60% reduction in total HDAC activity and a loss of cell viability. Cell death is dependent upon cell cycle progression, because differentiated, nonproliferating cells retain their viability. Furthermore, we observe increased mitotic defects, chromatin bridges, and micronuclei, suggesting HDAC1/2 are necessary for accurate chromosome segregation. Consistent with a critical role in the regulation of gene expression, microarray analysis of Hdac1/2-deleted cells reveals 1,708 differentially expressed genes. Significantly for the maintenance of stem cell self-renewal, we detected a reduction in the expression of the pluripotent transcription factors, Oct4, Nanog, Esrrb, and Rex1. HDAC1/2 activity is regulated through binding of an inositol tetraphosphate molecule (IP4) sandwiched between the HDAC and its cognate corepressor. This raises the important question of whether IP4 regulates the activity of the complex in cells. By rescuing the viability of double-knockout cells, we demonstrate for the first time (to our knowledge) that mutations that abolish IP4 binding reduce the activity of HDAC1/2 in vivo. Our data indicate that HDAC1/2 have essential and pleiotropic roles in cellular proliferation and regulate stem cell self-renewal by maintaining expression of key pluripotent transcription factors.

Histone deacetylase 1 and 2 are essential for normal T-cell development and genomic stability in mice.

Histone deacetylase 1 and 2 (HDAC1/2) regulate chromatin structure as the catalytic core of the Sin3A, NuRD and CoREST co-repressor complexes. To better understand the key pathways regulated by HDAC1/2 in the adaptive immune system and inform their exploitation as drug targets, we have generated mice with a T-cell specific deletion. Loss of either HDAC1 or HDAC2 alone has little effect, while dual inactivation results in a 5-fold reduction in thymocyte cellularity, accompanied by developmental arrest at the double-negative to double-positive transition. Transcriptome analysis revealed 892 misregulated genes in Hdac1/2 knock-out thymocytes, including down-regulation of LAT, Themis and Itk, key components of the T-cell receptor (TCR) signaling pathway. Down-regulation of these genes suggests a model in which HDAC1/2 deficiency results in defective propagation of TCR signaling, thus blocking development. Furthermore, mice with reduced HDAC1/2 activity (Hdac1 deleted and a single Hdac2 allele) develop a lethal pathology by 3-months of age, caused by neoplastic transformation of immature T cells in the thymus. Tumor cells become aneuploid, express increased levels of c-Myc and show elevated levels of the DNA damage marker, γH2AX. These data demonstrate a crucial role for HDAC1/2 in T-cell development and the maintenance of genomic stability.

Histone deacetylase 1 (HDAC1), but not HDAC2, controls embryonic stem cell differentiation.

Histone deacetylases (HDAC) 1 and 2 are highly similar enzymes that help regulate chromatin structure as the core catalytic components of corepressor complexes. Although tissue-specific deletion of HDAC1 and HDAC2 has demonstrated functional redundancy, germ-line deletion of HDAC1 in the mouse causes early embryonic lethality, whereas HDAC2 does not. To address the unique requirement for HDAC1 in early embryogenesis we have generated conditional knockout embryonic stem (ES) cells in which HDAC1 or HDAC2 genes can be inactivated. Deletion of HDAC1, but not HDAC2, causes a significant reduction in the HDAC activity of Sin3A, NuRD, and CoREST corepressor complexes. This reduced corepressor activity results in a specific 1.6-fold increase in histone H3 K56 acetylation (H3K56Ac), thus providing genetic evidence that H3K56Ac is a substrate of HDAC1. In culture, ES cell proliferation was unaffected by loss of either HDAC1 or HDAC2. Rather, we find that loss of HDAC1 affects ES cell differentiation. ES cells lacking either HDAC1 or HDAC2 were capable of forming embryoid bodies (EBs), which stimulates differentiation into the three primary germ layers. However, HDAC1-deficient EBs were significantly smaller, showed spontaneous rhythmic contraction, and increased expression of both cardiomyocyte and neuronal markers. In summary, our genetic study of HDAC1 and HDAC2 in ES cells, which mimic the embryonic epiblast, has identified a unique requirement for HDAC1 in the optimal activity of HDAC1/2 corepressor complexes and cell fate determination during differentiation.

Optimized whole-genome CRISPR interference screens identify ARID1A-dependent growth regulators in human induced pluripotent stem cells.

Perturbing expression is a powerful way to understand the role of individual genes, but can be challenging in important models. CRISPR-Cas screens in human induced pluripotent stem cells (iPSCs) are of limited efficiency due to DNA break-induced stress, while the less stressful silencing with an inactive Cas9 has been considered less effective so far. Here, we developed the dCas9-KRAB-MeCP2 fusion protein for screening in iPSCs from multiple donors. We found silencing in a 200 bp window around the transcription start site in polyclonal pools to be as effective as using wild-type Cas9 for identifying essential genes, but with much reduced cell numbers. Whole-genome screens to identify ARID1A-dependent dosage sensitivity revealed the PSMB2 gene, and enrichment of proteasome genes among the hits. This selective dependency was replicated with a proteasome inhibitor, indicating a targetable drug-gene interaction. Many more plausible targets in challenging cell models can be efficiently identified with our approach.

Transcriptional variability accelerates preleukemia by cell diversification and perturbation of protein synthesis.

Transcriptional variability facilitates stochastic cell diversification and can in turn underpin adaptation to stress or injury. We hypothesize that it may analogously facilitate progression of premalignancy to cancer. To investigate this, we initiated preleukemia in mouse cells with enhanced transcriptional variability due to conditional disruption of the histone lysine acetyltransferase gene Kat2a. By combining single-cell RNA sequencing of preleukemia with functional analysis of transformation, we show that Kat2a loss results in global variegation of cell identity and accumulation of preleukemic cells. Leukemia progression is subsequently facilitated by destabilization of ribosome biogenesis and protein synthesis, which confer a transient transformation advantage. The contribution of transcriptional variability to early cancer evolution reflects a generic role in promoting cell fate transitions, which, in the case of well-adapted malignancies, contrastingly differentiates and depletes cancer stem cells. That is, transcriptional variability confers forward momentum to cell fate systems, with differential multistage impact throughout cancer evolution.

The CADM1 tumor suppressor gene is a major candidate gene in MDS with deletion of the long arm of chromosome 11.

Myelodysplastic syndromes (MDS) represent a heterogeneous group of clonal hematopoietic stem cell disorders characterized by ineffective hematopoiesis leading to peripheral cytopenias and in a substantial proportion of cases to acute myeloid leukemia. The deletion of the long arm of chromosome 11, del(11q), is a rare but recurrent clonal event in MDS. Here, we detail the largest series of 113 cases of MDS and myelodysplastic syndromes/myeloproliferative neoplasms (MDS/MPN) harboring a del(11q) analyzed at clinical, cytological, cytogenetic, and molecular levels. Female predominance, a survival prognosis similar to other MDS, a low monocyte count, and dysmegakaryopoiesis were the specific clinical and cytological features of del(11q) MDS. In most cases, del(11q) was isolated, primary and interstitial encompassing the 11q22-23 region containing ATM, KMT2A, and CBL genes. The common deleted region at 11q23.2 is centered on an intergenic region between CADM1 (also known as Tumor Suppressor in Lung Cancer 1) and NXPE2. CADM1 was expressed in all myeloid cells analyzed in contrast to NXPE2. At the functional level, the deletion of Cadm1 in murine Lineage-Sca1+Kit+ cells modifies the lymphoid-to-myeloid ratio in bone marrow, although not altering their multilineage hematopoietic reconstitution potential after syngenic transplantation. Together with the frequent simultaneous deletions of KMT2A, ATM, and CBL and mutations of ASXL1, SF3B1, and CBL, we show that CADM1 may be important in the physiopathology of the del(11q) MDS, extending its role as tumor-suppressor gene from solid tumors to hematopoietic malignancies.

Minimal genome-wide human CRISPR-Cas9 library.

CRISPR guide RNA libraries have been iteratively improved to provide increasingly efficient reagents, although their large size is a barrier for many applications. We design an optimised minimal genome-wide human CRISPR-Cas9 library (MinLibCas9) by mining existing large-scale gene loss-of-function datasets, resulting in a greater than 42% reduction in size compared to other CRISPR-Cas9 libraries while preserving assay sensitivity and specificity. MinLibCas9 provides backward compatibility with existing datasets, increases the dynamic range of CRISPR-Cas9 screens and extends their application to complex models and assays.

SETBP1 overexpression acts in the place of class-defining mutations to drive FLT3-ITD-mutant AML.

Advances in cancer genomics have revealed genomic classes of acute myeloid leukemia (AML) characterized by class-defining mutations, such as chimeric fusion genes or in genes such as NPM1, MLL, and CEBPA. These class-defining mutations frequently synergize with internal tandem duplications in FLT3 (FLT3-ITDs) to drive leukemogenesis. However, ∼20% of FLT3-ITD-positive AMLs bare no class-defining mutations, and mechanisms of leukemic transformation in these cases are unknown. To identify pathways that drive FLT3-ITD mutant AML in the absence of class-defining mutations, we performed an insertional mutagenesis (IM) screening in Flt3-ITD mice, using Sleeping Beauty transposons. All mice developed acute leukemia (predominantly AML) after a median of 73 days. Analysis of transposon insertions in 38 samples from Flt3-ITD/IM leukemic mice identified recurrent integrations at 22 loci, including Setbp1 (20/38), Ets1 (11/38), Ash1l (8/38), Notch1 (8/38), Erg (7/38), and Runx1 (5/38). Insertions at Setbp1 led exclusively to AML and activated a transcriptional program similar, but not identical, to those of NPM1-mutant and MLL-rearranged AMLs. Guide RNA targeting of Setbp1 was highly detrimental to Flt3ITD/+/Setbp1IM+, but not to Flt3ITD/+/Npm1cA/+, AMLs. Also, analysis of RNA-sequencing data from hundreds of human AMLs revealed that SETBP1 expression is significantly higher in FLT3-ITD AMLs lacking class-defining mutations. These findings propose that SETBP1 overexpression collaborates with FLT3-ITD to drive a subtype of human AML. To identify genetic vulnerabilities of these AMLs, we performed genome-wide CRISPR-Cas9 screening in Flt3ITD/+/Setbp1IM+ AMLs and identified potential therapeutic targets, including Kdm1a, Brd3, Ezh2, and Hmgcr. Our study gives new insights into epigenetic pathways that can drive AMLs lacking class-defining mutations and proposes therapeutic approaches against such cases.

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

Altered transcription is a cardinal feature of acute myeloid leukemia (AML); however, exactly how mutations synergize to remodel the epigenetic landscape and rewire three-dimensional DNA topology is unknown. Here, we apply an integrated genomic approach to a murine allelic series that models the two most common mutations in AML: Flt3-ITD and Npm1c. We then deconvolute the contribution of each mutation to alterations of the epigenetic landscape and genome organization, and infer how mutations synergize in the induction of AML. Our studies demonstrate that Flt3-ITD signals to chromatin to alter the epigenetic environment and synergizes with mutations in Npm1c to alter gene expression and drive leukemia induction. These analyses also allow the identification of long-range cis-regulatory circuits, including a previously unknown superenhancer of Hoxa locus, as well as larger and more detailed gene-regulatory networks, driven by transcription factors including PU.1 and IRF8, whose importance we demonstrate through perturbation of network members.

Author Correction: The long non-coding RNA HOXB-AS3 regulates ribosomal RNA transcription in NPM1-mutated acute myeloid leukemia.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Publisher Correction: The long non-coding RNA HOXB-AS3 regulates ribosomal RNA transcription in NPM1-mutated acute myeloid leukemia.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

SOCS2 is part of a highly prognostic 4-gene signature in AML and promotes disease aggressiveness.

Acute myeloid leukemia (AML) is a heterogeneous disease with respect to its genetic and molecular basis and to patients´ outcome. Clinical, cytogenetic, and mutational data are used to classify patients into risk groups with different survival, however, within-group heterogeneity is still an issue. Here, we used a robust likelihood-based survival modeling approach and publicly available gene expression data to identify a minimal number of genes whose combined expression values were prognostic of overall survival. The resulting gene expression signature (4-GES) consisted of 4 genes (SOCS2, IL2RA, NPDC1, PHGDH), predicted patient survival as an independent prognostic parameter in several cohorts of AML patients (total, 1272 patients), and further refined prognostication based on the European Leukemia Net classification. An oncogenic role of the top scoring gene in this signature, SOCS2, was investigated using MLL-AF9 and Flt3-ITD/NPM1c driven mouse models of AML. SOCS2 promoted leukemogenesis as well as the abundance, quiescence, and activity of AML stem cells. Overall, the 4-GES represents a highly discriminating prognostic parameter in AML, whose clinical applicability is greatly enhanced by its small number of genes. The newly established role of SOCS2 in leukemia aggressiveness and stemness raises the possibility that the signature might even be exploitable therapeutically.

The long non-coding RNA HOXB-AS3 regulates ribosomal RNA transcription in NPM1-mutated acute myeloid leukemia.

Long non-coding RNAs (lncRNAs) are important regulatory molecules that are implicated in cellular physiology and pathology. In this work, we dissect the functional role of the HOXB-AS3 lncRNA in patients with NPM1-mutated (NPM1mut) acute myeloid leukemia (AML). We show that HOXB-AS3 regulates the proliferative capacity of NPM1mut AML blasts in vitro and in vivo. HOXB-AS3 is shown to interact with the ErbB3-binding protein 1 (EBP1) and guide EBP1 to the ribosomal DNA locus. Via this mechanism, HOXB-AS3 regulates ribosomal RNA transcription and de novo protein synthesis. We propose that in the context of NPM1 mutations, HOXB-AS3 overexpression acts as a compensatory mechanism, which allows adequate protein production in leukemic blasts.

UTX-mediated enhancer and chromatin remodeling suppresses myeloid leukemogenesis through noncatalytic inverse regulation of ETS and GATA programs.

The histone H3 Lys27-specific demethylase UTX (or KDM6A) is targeted by loss-of-function mutations in multiple cancers. Here, we demonstrate that UTX suppresses myeloid leukemogenesis through noncatalytic functions, a property shared with its catalytically inactive Y-chromosome paralog, UTY (or KDM6C). In keeping with this, we demonstrate concomitant loss/mutation of KDM6A (UTX) and UTY in multiple human cancers. Mechanistically, global genomic profiling showed only minor changes in H3K27me3 but significant and bidirectional alterations in H3K27ac and chromatin accessibility; a predominant loss of H3K4me1 modifications; alterations in ETS and GATA-factor binding; and altered gene expression after Utx loss. By integrating proteomic and genomic analyses, we link these changes to UTX regulation of ATP-dependent chromatin remodeling, coordination of the COMPASS complex and enhanced pioneering activity of ETS factors during evolution to AML. Collectively, our findings identify a dual role for UTX in suppressing acute myeloid leukemia via repression of oncogenic ETS and upregulation of tumor-suppressive GATA programs.

Preventing chemotherapy-induced myelosuppression by repurposing the FLT3 inhibitor quizartinib.

We describe an approach to inhibit chemotherapy-induced myelosuppression. We found that short-term exposure of mice to the FLT3 inhibitor quizartinib induced the transient quiescence of multipotent progenitors (MPPs). This property of quizartinib conferred marked protection to MPPs in mice receiving fluorouracil or gemcitabine. The protection resulted in the rapid recovery of bone marrow and blood cellularity, thus preventing otherwise lethal myelosuppression. A treatment strategy involving quizartinib priming that protected wild-type bone marrow progenitors, but not leukemic cells, from fluorouracil provided a more effective treatment than conventional induction therapy in mouse models of acute myeloid leukemia. This strategy has the potential to be extended for use in other cancers where FLT3 inhibition does not adversely affect the effectiveness of chemotherapy. Thus, the addition of quizartinib to cancer treatment regimens could markedly improve cancer patient survival and quality of life.