Myelodysplastic syndromes and related disorders (myelodysplasia) are a heterogeneous group of myeloid neoplasms showing deregulated blood cell production with evidence of myeloid dysplasia and a predisposition to acute myeloid leukaemia, whose pathogenesis is only incompletely understood. Here we report whole-exome sequencing of 29 myelodysplasia specimens, which unexpectedly revealed novel pathway mutations involving multiple components of the RNA splicing machinery, including U2AF35, ZRSR2, SRSF2 and SF3B1. In a large series analysis, these splicing pathway mutations were frequent (∼45 to ∼85%) in, and highly specific to, myeloid neoplasms showing features of myelodysplasia. Conspicuously, most of the mutations, which occurred in a mutually exclusive manner, affected genes involved in the 3'-splice site recognition during pre-mRNA processing, inducing abnormal RNA splicing and compromised haematopoiesis. Our results provide the first evidence indicating that genetic alterations of the major splicing components could be involved in human pathogenesis, also implicating a novel therapeutic possibility for myelodysplasia.

Contemporary treatment of pediatric acute myeloid leukemia (AML) requires the assignment of patients to specific risk groups. To explore whether expression profiling of leukemic blasts could accurately distinguish between the known risk groups of AML, we analyzed 130 pediatric and 20 adult AML diagnostic bone marrow or peripheral blood samples using the Affymetrix U133A microarray. Class discriminating genes were identified for each of the major prognostic subtypes of pediatric AML, including t(15;17) [ PML-RAR␣], t(8;21)[AML1-ETO], inv 16 [CBF␤-MYH11], MLL chimeric fusion genes, and cases classified as FAB-M7. When subsets of these genes were used in supervised learning algorithms, an overall classification accuracy of more than 93% was achieved. Moreover, we were able to use the expression signatures generated from the pediatric samples to accurately classify adult de novo AMLs with the same genetic lesions. The class discriminating genes also provided novel insights into the molecular pathobiology of these leukemias. Finally, using a combined pediatric data set of 130 AMLs and 137 acute lymphoblastic leukemias, we identified an expression signature for cases with MLL chimeric fusion genes irrespective of lineage. Surprisingly, AMLs containing partial tandem duplications of MLL failed to cluster with MLL chimeric fusion gene cases, suggesting a significant difference in their underlying mechanism of transformation. IntroductionAcute myeloid leukemia (AML) is a relatively rare malignancy in the pediatric population, comprising only 15% to 20% of the acute leukemias diagnosed in this age group. 1 Nevertheless, it remains a challenging disease with an inferior treatment outcome compared with pediatric acute lymphoblastic leukemia (ALL). Despite the introduction of new drugs, the aggressive use of allogeneic and autologous bone marrow transplantation, and improvements in supportive care, overall cure rates of AML in most contemporary treatment protocols remain below 60%. [2][3][4][5] Further improvements in cure rates are likely to come from a better understanding of both the molecular abnormalities responsible for the formation and growth of the leukemic cells, and the mechanisms underlying drug resistance.Increasingly, contemporary treatment protocols are incorporating methods for both accurate diagnosis and subsequent risk stratification. To achieve this requires not only distinguishing myeloblasts from lymphoblasts, but also assessing the extent of lineage commitment and differentiation, as well as the presence of specific molecular lesions or chromosomal abnormalities. Efforts over the last several decades have revealed AML to be a heterogeneous disease, with marked differences in cure rates between various genetic subtypes. [6][7][8][9] Acute promyelocytic leukemia was the first clear example of a clinically distinct AML subtype, being characterized by FAB-M3 morphology and expression of the t(15;17)-encoded promyelocytic leukemia-retinoic acid receptor alpha (PML-RAR␣) fusion protein. [10][11][12][13][14] ...

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).

To elucidate differential roles of mutations in myelodysplastic syndromes (MDS), we investigated clonal dynamics using whole-exome and/or targeted sequencing of 699 patients, of whom 122 were analyzed longitudinally. Including the results from previous reports, we assessed a total of 2,250 patients for mutational enrichment patterns. During progression, the number of mutations, their diversity and clone sizes increased, with alterations frequently present in dominant clones with or without their sweeping previous clones. Enriched in secondary acute myeloid leukemia (sAML; in comparison to high-risk MDS), FLT3, PTPN11, WT1, IDH1, NPM1, IDH2 and NRAS mutations (type 1) tended to be newly acquired, and were associated with faster sAML progression and a shorter overall survival time. Significantly enriched in high-risk MDS (in comparison to low-risk MDS), TP53, GATA2, KRAS, RUNX1, STAG2, ASXL1, ZRSR2 and TET2 mutations (type 2) had a weaker impact on sAML progression and overall survival than type-1 mutations. The distinct roles of type-1 and type-2 mutations suggest their potential utility in disease monitoring.