N6-methyladenosine (m6A) is an abundant nucleotide modification in mRNA that is required for the differentiation of mouse embryonic stem cells. However, it remains unknown whether m6A controls differentiation of normal and/or malignant myeloid hematopoietic cells. Here we show that shRNA-mediated depletion of the m6A-forming enzyme METTL3 in human hematopoietic stem/progenitor cells promotes differentiation coupled with reduced proliferation. Conversely, overexpression of wild-type METTL3, but not the catalytic-dead form of METTL3, inhibits differentiation and increases cell growth. METTL3 mRNA and protein is expressed more abundantly in acute myeloid leukemia (AML) cells compared to healthy hematopoietic stem/progenitor cells and other types of tumors. Furthermore, METTL3 depletion in humanmyeloid leukemia cell lines induces differentiation and apoptosis and delays leukemia in recipient mice in vivo. Single-nucleotide resolution mapping of m6A coupled with ribosome profiling reveals that m6A promotes the translation of c-MYC, BCL2 and PTEN mRNAs in human myeloid leukemia MOLM13 cells. Moreover, loss of METTL3 leads to increased levels of pAKT, which contributes to the differentiation effects of METTL3 depletion. Overall these results provide a rationale for therapeutic targeting of METTL3 in myeloid leukemia.

Acquired mutations are pervasive across normal tissues. However, our understanding of the processes that drive transformation of certain clones to cancer is limited. Here we study this phenomenon in the context of clonal hematopoiesis (CH) and the development of therapy-related myeloid neoplasms (tMN). We find mutations are selected differentially based on exposures. Mutations in ASXL1 are enriched in current or former smokers, whereas cancer therapy with radiation, platinum and topoisomerase II inhibitors preferentially selects for mutations in DNA damage response (DDR) genes ( TP53, PPM1D, CHEK2 ). Sequential sampling provides definitive evidence that DDR clones outcompete other clones when exposed to certain therapies. Among cases where CH was previously detected, the CH mutation was present at tMN diagnosis. We identify the molecular characteristics of CH that increase risk of tMN. The increasing implementation of clinical sequencing at diagnosis provides an opportunity to identify patients at risk of tMN for prevention strategies.

Tumor protein p53 (TP53) is the most frequently mutated gene in cancer 1,2. In patients with myelodysplastic syndromes (MDS), TP53 mutations are associated with high-risk disease 3,4 , rapid transformation to acute myeloid leukemia (AML) 5 , resistance to conventional therapies 6-8 and dismal outcomes 9. Consistent with the tumor-suppressive role of TP53, patients harbor both mono-and biallelic mutations 10. However, the biological and clinical implications of TP53 allelic state have not been fully investigated in MDS or any other cancer type. We analyzed 3,324 patients with MDS for TP53 mutations and allelic imbalances and delineated two subsets of patients with distinct phenotypes and outcomes. One-third of TP53-mutated patients had monoallelic mutations whereas two-thirds had multiple hits (multi-hit) consistent with biallelic targeting. Established associations with complex karyotype, few co-occurring mutations, high-risk presentation and poor outcomes were specific to multi-hit patients only. TP53 multi-hit state predicted risk of death and leukemic transformation independently of the Revised International Prognostic Scoring System (IPSS-R) 11. Surprisingly, monoallelic patients did not differ from TP53 wild-type patients in outcomes and response to therapy. This study shows that consideration of TP53 allelic state is critical for diagnostic and prognostic precision in MDS as well as in future correlative studies of treatment response. In collaboration with the International Working Group for Prognosis in MDS (Supplementary Table 1), we assembled a cohort of 3,324 peridiagnostic and treatment-naive patients with MDS or closely related myeloid neoplasms (Extended Data Fig. 1 and Supplementary Fig. 1). Genetic profiling included conventional G-banding analyses (CBA) and tumor-only, capture-based, next-generation sequencing (NGS) of a panel of genes recurrently mutated in MDS, as well as genome-wide copy number probes. Allele-specific copy number profiles were generated from NGS data using the CNACS algorithm 7 (see Methods and Code availability). An additional 1,120 samples derived from the Japanese MDS consortium (Extended Data Fig. 2) were used as a validation cohort. To study the effect of TP53 allelic state on genome stability, clinical presentation, outcome and response to therapy, we performed a detailed characterization of alterations at the TP53 locus. First, we assessed genome-wide allelic imbalances in the cohort of 3,324 patients, to include arm-level or focal (~3 Mb) ploidy alterations and regions of copy-neutral loss of heterozygosity (cnLOH) (Extended Data Fig. 3, Supplementary Figs. 2-4 and Methods).

Summary Myeloid malignancies, including acute myeloid leukemia (AML), arise from the expansion of hematopoietic stem/progenitor cells which acquire somatic mutations. Bulk molecular profiling suggests step-wise mutation acquisition, where mutant genes with high variant allele frequencies (VAFs) occur early in leukemogenesis and mutations with lower VAFs are thought to be acquired later 1 – 3 . Although bulk sequencing informs leukemia biology and prognostication, it cannot distinguish which mutations occur in the same clone(s), accurately measure clonal complexity, or definitively elucidate mutational order. To delineate the clonal framework of myeloid malignancies, we performed single cell mutational profiling on 146 samples from 123 patients. We found AML is dominated by a small number of clones, which frequently harbor co-occurring mutations in epigenetic regulators. Conversely, mutations in signaling genes often occur more than once in distinct subclones consistent with increasing clonal diversity. We next mapped clonal trajectories for each sample and uncovered mutation combinations that synergized to promote clonal expansion and dominance. Finally, we combined protein expression with mutational analysis to map somatic genotype and clonal architecture with immunophenotype. Our studies of single cell clonal architecture provides novel insights into the pathogenesis of myeloid transformation and how clonal complexity evolves with disease progression.

Somatic mutations in the isocitrate dehydrogenase 2 gene (IDH2) contribute to the pathogenesis of acute myeloid leukaemia (AML) through the production of the oncometabolite 2-hydroxyglutarate (2HG). Enasidenib (AG-221) is an allosteric inhibitor that binds to the IDH2 dimer interface and blocks the production of 2HG by IDH2 mutants. In a phase I/II clinical trial, enasidenib inhibited the production of 2HG and induced clinical responses in relapsed or refractory IDH2-mutant AML. Here we describe two patients with IDH2-mutant AML who had a clinical response to enasidenib followed by clinical resistance, disease progression, and a recurrent increase in circulating levels of 2HG. We show that therapeutic resistance is associated with the emergence of second-site IDH2 mutations in trans, such that the resistance mutations occurred in the IDH2 allele without the neomorphic R140Q mutation. The in trans mutations occurred at glutamine 316 (Q316E) and isoleucine 319 (I319M), which are at the interface where enasidenib binds to the IDH2 dimer. The expression of either of these mutant disease alleles alone did not induce the production of 2HG; however, the expression of the Q316E or I319M mutation together with the R140Q mutation in trans allowed 2HG production that was resistant to inhibition by enasidenib. Biochemical studies predicted that resistance to allosteric IDH inhibitors could also occur via IDH dimer-interface mutations in cis, which was confirmed in a patient with acquired resistance to the IDH1 inhibitor ivosidenib (AG-120). Our observations uncover a mechanism of acquired resistance to a targeted therapy and underscore the importance of 2HG production in the pathogenesis of IDH-mutant malignancies.