Analysis of 3760 hematologic malignancies reveals rare transcriptomic aberrations of driver genes.

BACKGROUND: Rare oncogenic driver events, particularly affecting the expression or splicing of driver genes, are suspected to substantially contribute to the large heterogeneity of hematologic malignancies. However, their identification remains challenging. METHODS: To address this issue, we generated the largest dataset to date of matched whole genome sequencing and total RNA sequencing of hematologic malignancies from 3760 patients spanning 24 disease entities. Taking advantage of our dataset size, we focused on discovering rare regulatory aberrations. Therefore, we called expression and splicing outliers using an extension of the workflow DROP (Detection of RNA Outliers Pipeline) and AbSplice, a variant effect predictor that identifies genetic variants causing aberrant splicing. We next trained a machine learning model integrating these results to prioritize new candidate disease-specific driver genes. RESULTS: We found a median of seven expression outlier genes, two splicing outlier genes, and two rare splice-affecting variants per sample. Each category showed significant enrichment for already well-characterized driver genes, with odds ratios exceeding three among genes called in more than five samples. On held-out data, our integrative modeling significantly outperformed modeling based solely on genomic data and revealed promising novel candidate driver genes. Remarkably, we found a truncated form of the low density lipoprotein receptor LRP1B transcript to be aberrantly overexpressed in about half of hairy cell leukemia variant (HCL-V) samples and, to a lesser extent, in closely related B-cell neoplasms. This observation, which was confirmed in an independent cohort, suggests LRP1B as a novel marker for a HCL-V subclass and a yet unreported functional role of LRP1B within these rare entities. CONCLUSIONS: Altogether, our census of expression and splicing outliers for 24 hematologic malignancy entities and the companion computational workflow constitute unique resources to deepen our understanding of rare oncogenic events in hematologic cancers.

Comparing malignant monocytosis across the updated WHO and ICC classifications of 2022.

ABSTRACT: The World Health Organization (WHO) classification of hematolymphoid tumors and the International Consensus Classification (ICC) of 2022 introduced major changes to the definition of chronic myelomonocytic leukemia (CMML). To assess its qualitative and quantitative implications for patient care, we started with 3311 established CMML cases (according to WHO 2017 criteria) and included 2130 oligomonocytosis cases fulfilling the new CMML diagnostic criteria. Applying both 2022 classification systems, 356 and 241 of oligomonocytosis cases were newly classified as myelodysplastic (MD)-CMML (WHO and ICC 2022, respectively), most of which were diagnosed as myelodysplastic syndrome (MDS) according to the WHO 2017 classification. Importantly, 1.5 times more oligomonocytosis cases were classified as CMML according to WHO 2022 than based on ICC, because of different diagnostic criteria. Genetic analyses of the newly classified CMML cases showed a distinct mutational profile with strong enrichment of MDS-typical alterations, resulting in a transcriptional subgroup separated from established MD and myeloproliferative CMML. Despite a different cytogenetic, molecular, immunophenotypic, and transcriptional landscape, no differences in overall survival were found between newly classified and established MD-CMML cases. To the best of our knowledge, this study represents the most comprehensive analysis of routine CMML cases to date, both in terms of clinical characterization and transcriptomic analysis, placing newly classified CMML cases on a disease continuum between MDS and previously established CMML.

The clinical and genomic landscape of patients with DDX41 variants identified during diagnostic sequencing.

Deleterious germ line variants in DDX41 are a common cause of genetic predisposition to hematologic malignancies, particularly myelodysplastic neoplasms (MDS) and acute myeloid leukemia (AML). Targeted next-generation sequencing was performed in a large cohort of sequentially recruited patients with myeloid malignancy, covering DDX41 as well as 30 other genes frequently mutated in myeloid malignancy. Whole genome transcriptome sequencing data was analyzed on a separate cohort of patients with a range of hematologic malignancies to investigate the spectrum of cancer predisposition. Altogether, 5737 patients with myeloid malignancies were studied, with 152 different DDX41 variants detected. Multiple novel variants were detected, including synonymous variants affecting splicing as demonstrated by RNA-sequencing. The presence of a somatic DDX41 variant was highly associated with DDX41 germ line variants in patients with MDS and AML, and we developed a statistical approach to incorporate the co-occurrence of a somatic DDX41 variant into germ line variant classification at a very strong level (as per the American College of Medical Genetics and Genomics/Association for Molecular Pathology guidelines). Using this approach, the MDS cohort contained 108 of 2865 (3.8%) patients with germ line likely pathogenic/pathogenic (LP/P) variants, and the AML cohort 106 of 2157 (4.9%). DDX41 LP/P variants were markedly enriched in patients with AML and MDS compared with those in patients with myeloproliferative neoplasms, B-cell neoplasm, and T- or B-cell acute lymphoblastic leukemia. In summary, we have developed a framework to enhance DDX41 variant curation as well as highlighted the importance of assessment of all types of genomic variants (including synonymous and multiexon deletions) to fully detect the landscape of possible clinically relevant DDX41 variants.

Germ line variant GFI1-36N affects DNA repair and sensitizes AML cells to DNA damage and repair therapy.

ABSTRACT: Growth factor independence 1 (GFI1) is a DNA-binding transcription factor and a key regulator of hematopoiesis. GFI1-36N is a germ line variant, causing a change of serine (S) to asparagine (N) at position 36. We previously reported that the GFI1-36N allele has a prevalence of 10% to 15% among patients with acute myeloid leukemia (AML) and 5% to 7% among healthy Caucasians and promotes the development of this disease. Using a multiomics approach, we show here that GFI1-36N expression is associated with increased frequencies of chromosomal aberrations, mutational burden, and mutational signatures in both murine and human AML and impedes homologous recombination (HR)-directed DNA repair in leukemic cells. GFI1-36N exhibits impaired binding to N-Myc downstream-regulated gene 1 (Ndrg1) regulatory elements, causing decreased NDRG1 levels, which leads to a reduction of O6-methylguanine-DNA-methyltransferase (MGMT) expression levels, as illustrated by both transcriptome and proteome analyses. Targeting MGMT via temozolomide, a DNA alkylating drug, and HR via olaparib, a poly-ADP ribose polymerase 1 inhibitor, caused synthetic lethality in human and murine AML samples expressing GFI1-36N, whereas the effects were insignificant in nonmalignant GFI1-36S or GFI1-36N cells. In addition, mice that received transplantation with GFI1-36N leukemic cells treated with a combination of temozolomide and olaparib had significantly longer AML-free survival than mice that received transplantation with GFI1-36S leukemic cells. This suggests that reduced MGMT expression leaves GFI1-36N leukemic cells particularly vulnerable to DNA damage initiating chemotherapeutics. Our data provide critical insights into novel options to treat patients with AML carrying the GFI1-36N variant.

Interface-guided phenotyping of coding variants in the transcription factor RUNX1 with SEUSS.

Understanding the consequences of single amino acid substitutions in cancer driver genes remains an unmet need. Perturb-seq provides a tool to investigate the effects of individual mutations on cellular programs. Here we deploy SEUSS, a Perturb-seq like approach, to generate and assay mutations at physical interfaces of the RUNX1 Runt domain. We measured the impact of 115 mutations on RNA profiles in single myelogenous leukemia cells and used the profiles to categorize mutations into three functionally distinct groups: wild-type (WT)-like, loss-of-function (LOF)-like and hypomorphic. Notably, the largest concentration of functional mutations (non-WT-like) clustered at the DNA binding site and contained many of the more frequently observed mutations in human cancers. Hypomorphic variants shared characteristics with loss of function variants but had gene expression profiles indicative of response to neural growth factor and cytokine recruitment of neutrophils. Additionally, DNA accessibility changes upon perturbations were enriched for RUNX1 binding motifs, particularly near differentially expressed genes. Overall, our work demonstrates the potential of targeting protein interaction interfaces to better define the landscape of prospective phenotypes reachable by amino acid substitutions.

Interplay of TP53 allelic state, blast count, and complex karyotype on survival of patients with AML and MDS.

Several clinical and genetic factors impact overall survival (OS) in myelodysplastic neoplasms (MDS) and acute myeloid leukemia (AML), including complex karyotype (CK), TP53 allelic state, and blast count. We analyzed the interplay of these factors by performing Cox regression analysis and by determining the frequency of TP53 single-hit (sh) and double-hit (dh) events and OS in MDS (n = 747) with <5% blasts, with ≥5% but <10% blasts, and ≥10% but <20% blasts and AML (n = 772). MDS with <5% blasts showed the best outcome, followed by with ≥5% but <10% blasts, and ≥10% but <20% blasts, and AML (median OS: 75, 54, 27, and 18 months, respectively). The same hierarchy was observed when each subgroup was divided into TP53sh, TP53dh, and without TP53 alterations (alt), revealing a dismal outcome of TP53dh in all subgroups (17, 10, 8, and 1 month[s], respectively). MDS with <5% blasts differed from the other subgroups by showing predominantly TP53sh (76% of TP53alt cases), and by an independent adverse impact of CK on OS (hazard ratio, 5.2; P < .001). The remaining subgroups displayed many similarities, with TP53dh found at high frequencies (67%, 91%, and 71%, respectively) and only TP53alt but not CK independently influencing OS, and TP53dh showing the strongest influence. When the total cohort was split based on TP53 state, only the blast count and not CK had an independent adverse impact on OS in all subgroups. Thus, TP53dh is the strongest prognostic factor, further supporting its integration into risk stratification guidelines and classification as a separate entity. However, the blast count also influences OS independent of TP53 state, whereas CK plays a minor prognostic role.

AML classification in the year 2023: How to avoid a Babylonian confusion of languages.

In parallel to the 5th edition of the World Health Organization Classification of Haematolymphoid Tumours (WHO 2022), an alternative International Consensus Classification (ICC) has been proposed. To evaluate the impact of the new classifications on AML diagnoses and ELN-based risk classification, we analyzed 717 MDS and 734 AML non-therapy-related patients diagnosed according to the revised 4th WHO edition (WHO 2017) by whole genome and transcriptome sequencing. In both new classifications, the purely morphologically defined AML entities decreased from 13% to 5%. Myelodysplasia-related (MR) AML increased from 22% to 28% (WHO 2022) and 26% (ICC). Other genetically-defined AML remained the largest group, and the abandoned AML-RUNX1 was mainly reclassified as AML-MR (WHO 2022: 77%; ICC: 96%). Different inclusion criteria of AML-CEBPA and AML-MR (i.a. exclusion of TP53 mutated cases according to ICC) were associated with differences in overall survival. In conclusion, both classifications focus on more genetics-based definitions with similar basic concepts and a large degree of agreement. The remaining non-comparability (e.g., TP53 mutated AML) needs additional studies to definitely answer open questions on disease categorization in an unbiased way.

SF3B1 mutated MDS: Blast count, genetic co-abnormalities and their impact on classification and prognosis.

Recently, MDS with mutated SF3B1 and blast count <5% was proposed as distinct entity with favorable prognosis by the international working group for the prognosis of MDS (IWG-PM), the 5th edition of the WHO classification and the International Consensus Classification. To further characterize this entity with respect to the genomic landscape, AML transformation rate and clinical outcome, we analyzed 734 MDS patients by whole genome sequencing. SF3B1 mutations were identified in 31% (n = 231), most frequently accompanied by TET2 mutations (29%). 144/231 (62%) SF3B1mut samples fulfilled entity criteria proposed by IWG-PM (SF3B1ent). These cases were associated with longer survival, lower AML transformation rate, normal karyotypes and harbored less accompanying mutations compared to SF3B1mut samples not falling into the proposed SF3B1 entity (SF3B1nent). Of SF3B1mut cases 7% (15/231; SF3B1ent: 3/144 [2%]; SF3B1nent: 12/87 [14%]) progressed to AML compared to 15% SF3B1 wild-type patients (75/503). Of these 15 SF3B1mut cases, 10 (67%) showed RUNX1 mutations at MDS or AML stage. Multivariate analysis revealed that del(5q) and RUNX1 mutations were independent negative prognostic factors for overall survival, while blast count >5% was not. In conclusion, SF3B1mut MDS has a favorable prognosis independent of blast count if karyotype and RUNX1 mutations are considered.

Rare germline alterations of myeloperoxidase predispose to myeloid neoplasms.

Myeloperoxidase (MPO) gene alterations with variable clinical penetrance have been found in hereditary MPO deficiency, but their leukemia association in patients and carriers has not been established. Germline MPO alterations were found to be significantly enriched in myeloid neoplasms: 28 pathogenic/likely pathogenic variants were identified in 100 patients. The most common alterations were c.2031-2 A > C, R569W, M519fs* and Y173C accounting for about half of the cases. While functional experiments showed that the marrow stem cell pool of Mpo-/- mice was not increased, using competitive repopulation demonstrated that Mpo-/- grafts gained growth advantage over MPO wild type cells. This finding also correlated with increased clonogenic potential after serial replating in the setting of H2O2-induced oxidative stress. Furthermore, we demonstrated that H2O2-induced DNA damage and activation of error-prone DNA repair may result in secondary genetic damage potentially predisposing to leukemia leukemic evolution. In conclusion, our study for the first time demonstrates that germline MPO variants may constitute risk alleles for MN evolution.

Pathophysiologic and clinical implications of molecular profiles resultant from deletion 5q.

BACKGROUND: Haploinsufficiency (HI) resulting from deletion of the long arm of chromosome 5 [del(5q)] and the accompanied loss of heterozygosity are likely key pathogenic factors in del(5q) myeloid neoplasia (MN) although the consequences of del(5q) have not been yet clarified. METHODS: Here, we explored mutations, gene expression and clinical phenotypes of 388 del(5q) vs. 841 diploid cases with MN [82% myelodysplastic syndromes (MDS)]. FINDINGS: Del(5q) resulted as founder (better prognosis) or secondary hit (preceded by TP53 mutations). Using Bayesian prediction analyses on 57 HI marker genes we established the minimal del(5q) gene signature that distinguishes del(5q) from diploid cases. Clusters of diploid cases mimicking the del(5q) signature support the overall importance of del(5q) genes in the pathogenesis of MDS in general. Sub-clusters within del(5q) patients pointed towards the inherent intrapatient heterogeneity of HI genes. INTERPRETATION: The underlying clonal expansion drive results from a balance between the "HI-driver" genes (e.g., CSNK1A1, CTNNA1, TCERG1) and the proapoptotic "HI-anti-drivers" (e.g., RPS14, PURA, SIL1). The residual essential clonal expansion drive allows for selection of accelerator mutations such as TP53 (denominating poor) and CSNK1A1 mutations (with a better prognosis) which overcome pro-apoptotic genes (e.g., p21, BAD, BAX), resulting in a clonal expansion. In summary, we describe the complete picture of del(5q) MN identifying the crucial genes, gene clusters and clonal hierarchy dictating the clinical course of del(5q) patients. FUNDING: Torsten Haferlach Leukemia Diagnostics Foundation. US National Institute of Health (NIH) grants R35 HL135795, R01HL123904, R01 HL118281, R01 HL128425, R01 HL132071, and a grant from Edward P. Evans Foundation.

Circulating microbial content in myeloid malignancy patients is associated with disease subtypes and patient outcomes.

Although recent work has described the microbiome in solid tumors, microbial content in hematological malignancies is not well-characterized. Here we analyze existing deep DNA sequence data from the blood and bone marrow of 1870 patients with myeloid malignancies, along with healthy controls, for bacterial, fungal, and viral content. After strict quality filtering, we find evidence for dysbiosis in disease cases, and distinct microbial signatures among disease subtypes. We also find that microbial content is associated with host gene mutations and with myeloblast cell percentages. In patients with low-risk myelodysplastic syndrome, we provide evidence that Epstein-Barr virus status refines risk stratification into more precise categories than the current standard. Motivated by these observations, we construct machine-learning classifiers that can discriminate among disease subtypes based solely on bacterial content. Our study highlights the association between the circulating microbiome and patient outcome, and its relationship with disease subtype.