Need Help?

An integrative model of pathway convergence in genetically heterogeneous blast crisis chronic myeloid leukemia (CML)

Targeted therapies against the BCR-ABL1 kinase have revolutionized treatment of chronic phase (CP) chronic myeloid leukemia (CML). In contrast, management of blast crisis (BC) CML remains challenging because BC cells acquire complex molecular alterations that confer stemness features to progenitor populations and resistance to BCR-ABL1 tyrosine kinase inhibitors. Comprehensive models of BC transformation have proved elusive because of the rarity and genetic heterogeneity of BC, but are important for developing biomarkers predicting BC progression and effective therapies. To better understand BC, we performed an integrated multiomics analysis of 74 CP and BC samples using whole-genome and exome sequencing, transcriptome and methylome profiling, and chromatin immunoprecipitation followed by high-throughput sequencing. Employing pathway-based analysis, we found the BC genome was significantly enriched for mutations affecting components of the polycomb repressive complex (PRC) pathway. While transcriptomically, BC progenitors were enriched and depleted for PRC1- and PRC2-related gene sets respectively. By integrating our data sets, we determined that BC progenitors undergo PRC-driven epigenetic reprogramming toward a convergent transcriptomic state. Specifically, PRC2 directs BC DNA hypermethylation, which in turn silences key genes involved in myeloid differentiation and tumor suppressor function via so-called epigenetic switching, whereas PRC1 represses an overlapping and distinct set of genes, including novel BC tumor suppressors. On the basis of these observations, we developed an integrated model of BC that facilitated the identification of combinatorial therapies capable of reversing BC reprogramming (decitabine+PRC1 inhibitors), novel PRC-silenced tumor suppressor genes (NR4A2), and gene expression signatures predictive of disease progression and drug resistance in CP.

Click on a Dataset ID in the table below to learn more, and to find out who to contact about access to these data

Dataset ID Description Technology Samples
EGAD00001002053 29
EGAD00001002054 24
EGAD00001002057 Illumina HiSeq 2000 29
EGAD00001002058 HiSeq X Ten Illumina HiSeq 2000 33
EGAD00001002059 Illumina HiSeq 2000 24
EGAD00001002060 Illumina HiSeq 2000 24
EGAD00001002061 Illumina HiSeq 2500 3
EGAD00001002062 Illumina HiSeq 2500 3
EGAD00010002209 Illumina Beadchip HT12v4 34
EGAD00010002210 Illumina Beadchip HT12v4 44
EGAD00010002211 Illumina Beadchip HT12v4 48
EGAD00010002216 Illumina 450K 65