Clonal dynamics and copy number variants by single‐cell analysis in leukemic evolution of myeloproliferative neoplasms

Transformation from chronic (CP) to blast phase (BP) in myeloproliferative neoplasm (MPN) remains poorly characterized, and no specific mutation pattern has been highlighted. BP-MPN represents an unmet need, due to its refractoriness to treatment and dismal outcome. Taking advantage of the granularity provided by single-cell sequencing (SCS), we analyzed paired samples of CP and BP in 10 patients to map clonal trajectories and interrogate target copy number variants (CNVs). Already at diagnosis, MPN present as oligoclonal diseases with varying ratio of mutated and wild-type cells, including cases where normal hematopoiesis was entirely surmised by mutated clones. BP originated from increasing clonal complexity, either on top or independent of a driver mutation, through acquisition of novel mutations as well as accumulation of clones harboring multiple mutations, that were detected at CP by SCS but were missed by bulk sequencing. There were progressive copy-number imbalances from CP to BP, that configured distinct clonal profiles and identified recurrences in genes including NF1, TET2, and BCOR, suggesting an additional level of complexity and contribution to leukemic transformation. EZH2 emerged as the gene most frequently affected by single nucleotide and CNVs, that might result in EZH2/PRC2-mediated transcriptional deregulation, as supported by combined scATAC-seq and snRNA-seq analysis of the leukemic clone in a representative case. Overall, findings provided insights into the pathogenesis of MPN-BP, identified CNVs as a hitherto poorly characterized mechanism and point to EZH2 dysregulation as target. Serial assessment of clonal dynamics might potentially allow early detection of impending disease transformation, with therapeutic implications.