Solary, Éric

PURPOSE Chronic myelomonocytic leukemia (CMML) is a rare myeloid neoplasm characterized by clinical heterogeneity and is associated with poor outcomes. To date, limited molecular information has been incorporated into disease classification and risk stratification. We aimed to integrate genomic features into the clinical decision-making process for CMML. PATIENTS AND METHODS We analyzed a retrospective cohort of 3013 patients with CMML (training set) and a prospective population of 516 patients (validation set). Using an innovative framework for multimodal data analysis, we developed molecular-based disease taxonomy and prognostication. RESULTS Unsupervised clustering identified nine entities with distinct genomic features and outcomes ( P < .001), including splicing machinery, transcription factors, signal transduction and tyrosine kinase pathways aberrations, and high-risk molecular signatures. Notably, 15% of patients showed molecular/clinical overlap with other myeloid neoplasms. We integrated molecular and clinical information to build the international CMML Prognostic Scoring System (iCPSS), incorporating mutations in nine genes together with hematologic parameters and cytogenetic abnormalities. The iCPSS identified five groups with distinct probability of overall and leukemia-free survival in both training and validation cohorts ( P < .001), outperforming existing prognostic models. Importantly, 55% of patients were reassigned to higher or lower risk groups by the iCPSS. Decision analysis demonstrated that iCPSS could refine the optimal timing of allogeneic transplantation at the individual level; compared with conventional prognostic tools, iCPSS-based decision modeling changed transplantation strategy in 31% of cases, resulting in a significant gain-in-life expectancy for eligible patient population ( P < .001). A federated learning platform was implemented to enable continuous, privacy-preserving model update across multiple centers. CONCLUSION Molecular information improves CMML classification and prognostication, supports more effective clinical decision making, and potentially refines the design of clinical trials.

Abstract Introduction Blast transformation (BT) in CMML occurs in 15-30% of patients (pts) with a median overall survival (OS) &lt;6 months (mo). While clonal evolution can be observed, biological mechanisms of BT remain unknown for many pts. Using paired samples from chronic phase (CP) and at BT, this study defines the cytogenetic and molecular landscape of BT-CMML. Methods BT-CMML pts from Mayo Clinic (n=153, 31%), MD Anderson (n=195, 40%), and 7 centers of the Groupe Francophone des Myélodysplasies (n=145, 29%) were included after IRB approval. ICC 2022 criteria defined CMML and AML. Pt data were collected at initial diagnosis or referral at CP and again at BT. Clonal dominance was inferred from a natural log-fold difference of mutation variant allele fractions (VAF). OS was censored at allogeneic stem cell transplant (alloSCT). Single-cell multiomic analyses used the Tapestri platform (Mission Bio). Statistics used R. Results Among 1,883 cases of CMML with a median follow up of 96 mo, 493 (26%) underwent BT. The median age at CP and BT was 69 and 71 years, respectively. BT-CMML pts were predominantly male (69%) and had proliferative CMML (53%) or CMML-1 (69%) in CP; 39% and 33% were considered high-risk by CPSS-Mol and BLAST-Mol criteria, respectively. The median time to BT was 16 mo (interquartile range 7-35 mo) and was not different in males vs females (p=0.60). Compared to their CP parameters (p&lt;0.001 each), pts at BT had lower median hemoglobin (10.8 vs 9.2 g/dL) and platelets (100 vs 44 x109/L) but higher leukocytes (13.9 vs 16.3 x109/L), peripheral blasts (0% vs 14%), and bone marrow blasts (6% vs 32%). Blasts exhibited mostly monocytic (57%) or immature myeloid (40%) differentiation. Abnormal karyotypes (AK) were more prevalent at BT (25% vs 62%, p&lt;0.001). Among 318 pts with paired karyotype data, 221 (69%) had a normal karyotype (NK) in CP, of which 123 (56%) remained normal at BT and 98 (44%) acquired new abnormalities. Conversely, 97 (31%) pts had an AK in CP, of which 67 (69%) acquired abnormalities or increased clonal metaphases (6%), 23 (24%) remained stable, and 7 (7%) reverted to a NK at BT. Overall, 46% acquired and 8% lost a clonal cytogenetic abnormality, while 46% remained unchanged. Numeric (37%) and structural (30%) changes were most common. Complex (12% vs 3%) and monosomal (5% vs 2%) karyotypes were more frequent at BT than in CP. Only 6 (2%) pts acquired an AML-defining rearrangement at BT. Mutational data were available for 248 paired cases. The number of somatic mutations was higher at BT (median 3 vs 4, p&lt;0.001). Mutations in RUNX1 (27% vs 39%), NRAS (23% vs 33%), KRAS (11% vs 17%), FLT3 (6% vs 15%), NPM1 (0% vs 11%), IDH2 (5% vs 11%), and CEBPA (3% vs 10%) were more frequent at BT (all p&lt;0.05). At BT, the average VAF of the dominant clone increased in 32%, decreased in 35%, and remained unchanged (delta ≤5%) in 33%. The most frequent dominant clones both in CP and at BT had mutations in TET2 (8% vs 4%) or co-dominant mutations involving SRSF2-TET2 (7% vs 3%). Mutational clonal evolution was observed in 72%, with mutations gained in 45% and lost from the dominant clone in 19%. The most gained mutations at BT were in NRAS (15%) and RUNX1 (13%). Single cell proteogenomic analyses in 4 selected CP and BT pts shed light on the impact of NRAS and RUNX1 mutations on clonal hierarchies and differentiation. Among 196 fully evaluable cases, a new cytogenetic and/or molecular aberration was uniquely gained in 106 (54%) cases and lost in 9 (5%), with 40 (20%) both gaining and losing a genetic alteration. Forty-one (21%) patients had no identifiable genetic change between CP and BT, suggesting driver alterations outside of that detected by routine analyses. Median OS was 31 mo from CMML diagnosis and 5.6 mo from BT. Compared to supportive care, treatment in CP did not prolong median time to BT (14 vs 17 mo, p=0.91) or OS (35 vs 29 mo, p=0.05), regardless of intensity. After BT, median OS was 1.2 mo with supportive care (n=83), 8.7 mo with a hypomethylating agent (HMA, n=39), 10 mo with HMA plus venetoclax (n=37), 7.8 mo with chemotherapy (n=91), and 31 mo in pts who underwent alloSCT (n=57). Conclusions In this large BT-CMML cohort, cytogenetic (46% gained, 8% lost) or mutational (45% gained, 19% lost) driver events were identified in 79% of pts, with 21% having no identifiable changes by conventional testing. Whole genome/transcriptome analyses are underway in this subgroup to identify novel drivers of BT.

Abstract Background The expansion of hematopoietic stem and progenitor cells (HSPCs) in the bone marrow is universally associated with adverse outcomes and disease progression in all chronic myeloid neoplasms, to include chronic myelomonocytic leukemia (CMML). However, HSPCs are a heterogeneous cellular compartment and the impact of its composition on clinical outcomes remains poorly defined. Single cell RNA-sequencing (scRNA-seq) represents an optimal strategy to annotate cellular composition but no large clinically annotated scRNA-seq datasets are currently available. To address this, we have embarked to generate a clinically annotated 500 patient CMML scRNA-seq atlas from CD34+ enriched BMNCs to establish among the largest clinically and genetically annotated single cell references for any chronic myeloid malignancy. Here we report the interim analysis of 113 clinically annotated cases that demonstrates massive reshaping of HSPC composition that is independently a predictor of overall survival. Methods Single-cell RNA expression profiling was performed on 113 CD34-enriched CMML BMNC samples using 10x Genomics GEM-X Universal 3' Gene Expression v4 (4-plex) chemistry. Libraries were sequenced on NovaSeq S4 flow cells in two runs, and data were processed using Cell Ranger to generate gene expression matrices. A total of 502,659 single cells from 113 samples were processed using Seurat. The upper threshold for gene count was defined as the mean number of detected genes plus two standard deviations data normalization was performed using the log-normalization method. Batch correction and integration were conducted with Harmony to mitigate technical variation across pools. Descriptive statistics were used to tabulate clinical and genetic parameters. KM curves, log-rank test, and cox regression were used to calculate survival, test for statistical differences, and perform multivariate analysis, respectively. Results To derive high resolution cell annotations, we employed the SingleR algorithm using a large hematopoietic reference of over 300,000 single cells (Zhang, X. Nature Immunology, 2024). Using this approach a diverse set of cell types were identified that included HSCs, 23 myeloid progenitors, 18 erythroid progenitors, and 18 lymphoid progenitors. This diverse cellular repertoire displayed massive reshaping relative to the normal CD34 enriched reference including a &gt;53x expansion of multilineage GMPs and a &gt;100x depletion of B-cell progenitors consistent with the known myeloid skewing seen in this disease. To understand whether expansion of specific progenitors was associated with outcomes we iteratively evaluated key progenitor populations and observed that expansions of HSC (HR 0.468 p&lt;0.005) and MEP (HR 0.485, p&lt;0.005) were associated with improved survival while GMP (HR 1.729, p=0.029) expansions were associated with inferior survival when using a cut point derived with a maximally selected ranks statistics approach. To achieve a more comprehensive metric of HSPC composition we explored GMP ratios to include GMP/HSC ratios and GMP/MEP ratios and observed that both were associated with inferior survival in univariate and multivariate when considering clinical and genetic features such as the CPSS-molecular prognostic score and FAB classification. When considering both GMP/HSC and GMP/MEP ratios in a multivariate model only GMP/MEP ratios and CPSS molecular remained independently prognostic. As an orthogonal approach, pseudotime projections were carried out using a reference map from Setty et al., 2019 to calculate differentiation potential and branch probabilities. This approach enabled both cell identity assignment and placement of cells along the hematopoietic differentiation trajectory. Samples were classified as monocyte-biased, MEP-biased, or normal-like, as previously described (Ferrall-Fairbanks et al. 2022) and demonstrated that increased GMP ratios are associated with a monocytic biased trajectory and that this trajectory is also associated with overall survival. Conclusion This study represents an interim analysis of among the largest scRNA-seq cohorts generated in chronic myeloid malignancies. This enabled us, for the first time, to demonstrate that HSPC composition is independently prognostic even when considering modern CMML prognostic models. As this cohort continues to increase, validation of these findings and examining the impact of therapy on HPSC composition in serial samples will be prioritized.

Background. Study of myeloid neoplasms (MN) has been rapidly transformed by genome characterization. Gene mutations have been reported to influence disease phenotype and risk of disease progression, and the evaluation of mutational status can provide valuable insights to improve the decision-making process. In this international study, we evaluated the clinical impact of mutational screening in patients with CMML, which is characterized by significant clinical and genomic heterogeneity, including a substantial proportion of patients with poor outcomes and unmet needs. Currently, only limited molecular information has been incorporated into CMML classification and prognostication. Methods. We retrospectively collected 3182 CMML patients. We used MOSAIC, an AI-based framework for multimodal analysis in rare cancers (PMID: 38875514) to develop innovative, molecular-based tools for classification and prognostication. An independent validation of the results on 516 prospectively collected patients was performed. Results. Bayesian Networs and Hierarchical Dirichlet Processes were used to identify genomic associations and define a CMML molecular classification. We identified 11 clusters with distinct clinical and genomic features, including splicing machinery, signal transduction and tyrosine kinase pathways aberrations, and high-risk molecular signatures (i.e., AML-like and TP53 mutations). Up to 15% of patients showed clear genomic overlap with other MN. Importantly, extensive multi-color flow cytometry on T lymphocytes, NK and myeloid cells (BD FACSymphony A5 Cell Analyzer) revealed specific immunologic and inflammation signatures associated with molecular subgroups. We then integrated molecular and clinical information to build an international CMML Prognostic Scoring System (iCPSS). We assessed the performance of different statistical and AI-based methods; L1-penalized Cox Model resulted as the best performing method. Selected features included hematological parameters (WBC, Hb, PLT and marrow blasts), cytogenetic abnormalities, and mutations in 10 genes (ASXL1, DNMT3A, EZH2, NRAS, RUNX1, SETBP1, STAG2, TET2, TP53, U2AF1). By Bayesian thresholds optimization, we identified 5 risk classes (i.e., very-low, low, intermediate, high and very-high) with median overall survival (OS) ranging from 99 to 9 months (P&lt;.001), and median leukemia-free survival (LFS) ranging from 210 to 18 months (P&lt;.001). iCPSS provided better patient discrimination across all clinical endpoints compared to currently available prognostic tools, with a Concordance Index (CI) for OS of 0.75 versus 0.62-0.64, respectively. Notably, compared to available scoring systems, up to 40% of patients were reassigned to higher or lower risk classes by the iCPSS. External validation on an independent prospective cohort confirmed that the performance of iCPSS is superior to existing prognostic models (CI for OS 0.71 vs 0.54-0.61, respectively) To demonstrate the clinical utility of iCPSS, we focused on 753 patients who underwent allogeneic stem cell transplantation (HSCT). In multivariable analysis, iCPSS stratified the probability of OS post-HSCT (P&lt;.001) and identified groups of patients with different probabilities of disease relapse, ranging from 9% to 62% (P&lt;.001). We employed a clinical- and genomic-based decision support system (PMID: 38723212) to determine the optimal timing of HSCT. Preliminary analyses indicated that higher-risk patients according to iCPSS (i.e., intermediate, high, and very high) benefit from an immediate HSCT strategy. Lower-risk patients (i.e., very low and low) had life expectancy maximized with delayed HSCT. Modeling decision analysis using iCPSS compared to conventional scores, resulted in a change in transplantation policy for a significant proportion of patients (up to 25%). Conclusion. Molecular information significantly improves the classification of CMML patients, providing a basis for refining diagnostic boundaries with other MN and for a more rational inclusion of patients in clinical trials. The iCPSS demonstrated superior performance over currently available scoring systems and improved HSCT decision making process.

Chronic myelomonocytic leukemia (CMML) is a severe myeloid malignancy with limited therapeutic options. Single-cell analysis of clonal architecture demonstrates early clonal dominance with few residual WT hematopoietic stem cells. Circulating myeloid cells of the leukemic clone and the cytokines they produce generate a deleterious inflammatory climate. Our hypothesis is that therapeutic control of the inflammatory component in CMML could contribute to stepping down disease progression. The present study explored the contribution of immature granulocytes (iGRANs) to CMML progression. iGRANs were detected and quantified in the peripheral blood of patients by spectral and conventional flow cytometry. Their accumulation was a potent and independent poor prognostic factor. These cells belong to the leukemic clone and behaved as myeloid-derived suppressor cells. Bulk and single-cell RNA-Seq revealed a proinflammatory status of iGRAN that secreted multiple cytokines of which CXCL8 was at the highest level. This cytokine inhibited the proliferation of WT but not CMML hematopoietic stem and progenitor cells (HSPCs) in which CXCL8 receptors were downregulated. CXCL8 receptor inhibitors and CXCL8 blockade restored WT HSPC proliferation, suggesting that relieving CXCL8 selective pressure on WT HSPCs is a potential strategy to slow CMML progression and restore some healthy hematopoiesis.

Abstract TET2 mutations (mTET2) are common genetic events in myeloid malignancies and clonal hematopoiesis. These mutations arise in the founding clone and are implicated in many clinical sequelae associated with oncogenic feedforward inflammatory circuits. However, the direct downstream effector of mTET2 responsible for the potentiation of these inflammatory circuits is unknown. To address this, we performed scRNA-seq and scATAC-seq in patients with COVID-19 with and without TET2-mutated clonal hematopoiesis reasoning that inflammation from COVID-19 may highlight critical downstream transcriptional targets of mTET2. Using this approach, we identified metastasis-associated lung adenoma transcript 1 (MALAT1), a therapeutically tractable lncRNA, as a central downstream effector of mTET2 that is both necessary and sufficient to induce the oncogenic proinflammatory features of mTET2 in vivo. We also elucidate the mechanism by which mTET2 upregulate MALAT1 and describe an interaction between MALAT1 and p65, which leads to RNA “shielding” from protein phosphatase 2A dephosphorylation, thus preventing resolution of inflammatory signaling. Significance: This work identifies MALAT1 as a requisite downstream effector of oncogenic feedforward inflammatory circuits necessary for the development of TET2-mutated CH and fulminant myeloid malignancy. We elucidate a novel mechanism by which MALAT1 “shields” p65 from dephosphorylation to potentiate this circuit and nominate MALAT1 inhibition as a future therapeutic strategy.

Abstract Background: COVID-19 causes significant morbidity and mortality, albeit with considerable heterogeneity among affected individuals. It remains unclear which host factors determine disease severity and survival. Given the propensity of clonal hematopoiesis (CH) to promote inflammation in healthy individuals, we investigated its effect on COVID-19 outcomes. Methods: We performed a multi-omics interrogation of the genome, epigenome, transcriptome, and proteome of peripheral blood mononuclear cells from COVID-19 patients (n=227). We obtained clinical data, laboratory studies, and survival outcomes. We determined CH status and TET2-related DNA methylation. We performed single-cell proteogenomics to understand clonal composition in relation to cell phenotype. We interrogated single-cell gene expression in isolation and in conjunction with DNA accessibility. We integrated these multi-omics data to understand the effect of CH on clonal composition, gene expression, methylation of cis-regulatory elements, and lineage commitment in COVID-19 patients. We performed shRNA knockdowns to validate the effect of one candidate transcription factor in myeloid cell lines. Results: The presence of CH was strongly associated with COVID-19 severity and all-cause mortality, independent of age (HR 3.48, 95% CI 1.45-8.36, p=0.005). Differential methylation of promoters and enhancers was prevalent in TET2-mutant, but not DNMT3A-mutant CH. TET2-mutant CH was associated with enhanced classical/intermediate monocytosis and single-cell proteogenomics confirmed an enrichment of TET2 mutations in these cell types. We identified cell-type specific gene expression changes associated with TET2 mutations in 102,072 single cells (n=34). Single-cell RNA-seq confirmed the skewing of hematopoiesis towards classical and intermediate monocytes and demonstrated the downregulation of EGR1 (a transcription factor important for monocyte differentiation) along with up-regulation of the lncRNA MALAT1 in monocytes. Combined scRNA-/scATAC-seq in 43,160 single cells (n=18) confirmed the skewing of hematopoiesis and up-regulation of MALAT1 in monocytes along with decreased accessibility of EGR1 motifs in known cis-regulatory elements. Using myeloid cell lines for functional validation, shRNA knockdowns of EGR1 confirmed the up-regulation of MALAT1 (in comparison to wildtype controls). Conclusions: CH is an independent prognostic factor in COVID-19 and skews hematopoiesis towards monocytosis. TET2-mutant CH is characterized by differential methylation and accessibility of enhancers binding myeloid transcriptions factors including EGR1. The ensuing loss of EGR1 expression in monocytes causes MALAT1 overexpression, a factor known to promote monocyte differentiation and inflammation. These data provide a mechanistic insight to the adverse prognostic impact of CH in COVID-19. Citation Format: Moritz Binder, Terra L. Lasho, Wazim Mohammed Ismail, Nana A. Ben-Crentsil, Jenna A. Fernandez, Minsuk Kim, Susan M. Geyer, Amelia Mazzone, Christy M. Finke, Abhishek A. Mangaonkar, Jeong-Heon Lee, Kwan Hyun Kim, Vernadette A. Simon, Fariborz Rakhshan Rohakthar, Amik Munankarmy, Susan M. Schwager, Jonathan J. Harrington, Melissa R. Snyder, Nathalie M. Droin, Eric Solary, Keith D. Robertson, Eric D. Wieben, Eric Padron, Nicholas Chia, Alexandre Gaspar-Maia, Mrinal M. Patnaik. Enhancer deregulation inTET2-mutant clonal hematopoiesis is associated with increased COVID-19 severity and mortality [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 806.

Abstract To determine the survival benefit of allogeneic hematopoietic cell transplantation (allo-HCT) in chronic myelomonocytic leukemias (CMML), we assembled a retrospective cohort of CMML patients 18-70 years old diagnosed between 2000 and 2014 from an international CMML dataset (n = 730) and the EBMT registry (n = 384). The prognostic impact of allo-HCT was analyzed through univariable and multivariable time-dependent models and with a multistate model, accounting for age, sex, CMML prognostic scoring system (low or intermediate-1 grouped as lower-risk, intermediate-2 or high as higher-risk) at diagnosis, and AML transformation. In univariable analysis, lower-risk CMMLs had a 5-year overall survival (OS) of 20% with allo-HCT vs 42% without allo-HCT (P &lt; .001). In higher-risk patients, 5-year OS was 27% with allo-HCT vs 15% without allo-HCT (P = .13). With multistate models, performing allo-HCT before AML transformation reduced OS in patients with lower-risk CMML, and a survival benefit was predicted for men with higher-risk CMML. In a multivariable analysis of lower-risk patients, performing allo-HCT before transformation to AML significantly increased the risk of death within 2 years of transplantation (hazard ratio [HR], 3.19; P &lt; .001), with no significant change in long-term survival beyond this time point (HR, 0.98; P = .92). In higher-risk patients, allo-HCT significantly increased the risk of death in the first 2 years after transplant (HR 1.46; P = .01) but not beyond (HR, 0.60; P = .09). Performing allo-HCT before AML transformation decreases life expectancy in lower-risk patients but may be considered in higher-risk patients.

Abstract Mutations in chromatin regulator ASXL1 are frequently identified in myeloid malignancies, in particular ∼40% of patients with chronic myelomonocytic leukemia (CMML). ASXL1 mutations are associated with poor prognosis in CMML and significantly co-occur with NRAS mutations. Here, we show that concurrent ASXL1 and NRAS mutations defined a population of CMML patients who had shorter leukemia-free survival than those with ASXL1 mutation only. Corroborating this human data, Asxl1−/− accelerated CMML progression and promoted CMML transformation to acute myeloid leukemia (AML) in NrasG12D/+ mice. NrasG12D/+;Asxl1−/− (NA) leukemia cells displayed hyperactivation of MEK/ERK signaling, increased global levels of H3K27ac, upregulation of Flt3. Moreover, we find that NA-AML cells overexpressed all the major inhibitory immune checkpoint ligands: programmed death-ligand 1 (PD-L1)/PD-L2, CD155, and CD80/CD86. Among them, overexpression of PD-L1 and CD86 correlated with upregulation of AP-1 transcription factors (TFs) in NA-AML cells. An AP-1 inhibitor or short hairpin RNAs against AP-1 TF Jun decreased PD-L1 and CD86 expression in NA-AML cells. Once NA-AML cells were transplanted into syngeneic recipients, NA-derived T cells were not detectable. Host-derived wild-type T cells overexpressed programmed cell death protein 1 (PD-1) and T-cell immunoreceptor with immunoglobulin and ITIM domains (TIGIT) receptors, leading to a predominant exhausted T-cell phenotype. Combined inhibition of MEK and BET resulted in downregulation of Flt3 and AP-1 expression, partial restoration of the immune microenvironment, enhancement of CD8 T-cell cytotoxicity, and prolonged survival in NA-AML mice. Our study suggests that combined targeted therapy and immunotherapy may be beneficial for treating secondary AML with concurrent ASXL1 and NRAS mutations.