Renneville, Aline

Abstract Introduction Blast transformation (BT) in CMML occurs in 15-30% of patients (pts) with a median overall survival (OS) <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<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<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<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<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 >53x expansion of multilineage GMPs and a >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<0.005) and MEP (HR 0.485, p<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.

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.