Integrated genomic sequencing in myeloid blast crisis chronic myeloid leukemia (MBC
Scientific Reports volume 12, Article number: 12816 (2022) Cite this article
1627 Accesses
1 Citations
1 Altmetric
Metrics details
An Author Correction to this article was published on 13 January 2023
This article has been updated
Chronic myeloid leukemia (CML) is a model of leukemogenesis in which the exact molecular mechanisms underlying blast crisis still remained unexplored. The current study identified multiple common and rare important findings in myeloid blast crisis CML (MBC-CML) using integrated genomic sequencing, covering all classes of genes implicated in the leukemogenesis model. Integrated genomic sequencing via Whole Exome Sequencing (WES), Chromosome-seq and RNA-sequencing were conducted on the peripheral blood samples of three CML patients in the myeloid blast crisis. An in-house filtering pipeline was applied to assess important variants in cancer-related genes. Standard variant interpretation guidelines were used for the interpretation of potentially important findings (PIFs) and potentially actionable findings (PAFs). Single nucleotide variation (SNV) and small InDel analysis by WES detected sixteen PIFs affecting all five known classes of leukemogenic genes in myeloid malignancies including signaling pathway components (ABL1, PIK3CB, PTPN11), transcription factors (GATA2, PHF6, IKZF1, WT1), epigenetic regulators (ASXL1), tumor suppressor and DNA repair genes (BRCA2, ATM, CHEK2) and components of spliceosome (PRPF8). These variants affect genes involved in leukemia stem cell proliferation, self-renewal, and differentiation. Both patients No.1 and No.2 had actionable known missense variants on ABL1 (p.Y272H, p.F359V) and frameshift variants on ASXL1 (p.A627Gfs*8, p.G646Wfs*12). The GATA2-L359S in patient No.1, PTPN11-G503V and IKZF1-R208Q variants in the patient No.3 were also PAFs. RNA-sequencing was used to confirm all of the identified variants. In the patient No. 3, chromosome sequencing revealed multiple pathogenic deletions in the short and long arms of chromosome 7, affecting at least three critical leukemogenic genes (IKZF1, EZH2, and CUX1). The large deletion discovered on the short arm of chromosome 17 in patient No. 2 resulted in the deletion of TP53 gene as well. Integrated genomic sequencing combined with RNA-sequencing can successfully discover and confirm a wide range of variants, from SNVs to CNVs. This strategy may be an effective method for identifying actionable findings and understanding the pathophysiological mechanisms underlying MBC-CML, as well as providing further insights into the genetic basis of MBC-CML and its management in the future.
Chronic Myeloid Leukemia (CML) is a unique model for the evolution of cancer and is classified as a triphasic myeloproliferative neoplasm based on clinical and pathological characteristics. Premalignant leukemia stem cells (LSCs) are, in fact, generated in the bone marrow nich by unknown mutagenesis processes1,2. In this step, the disease may be undetectable for a decade or more. Further oncogenic processes cause LSCs to develop into leukemia progenitor cells (LPCs)3. Leukemogenic alterations mainly affect five classes of regulatory proteins: signaling pathway components, transcription factors (TFs), epigenetic regulators (ERs), tumor suppressor genes (TSGs), and components of the spliceosome4.
Without therapeutic interventions or in the case of drug resistance, CML will progress to accelerated phase (AP) and then into an acute leukemia phase that is called blastic phase or blast crisis5. Even in the advent of Tyrosin Kinase Inhibitor therapy, blastic phase survival is less than 12 months. In general, blast crisis in chronic myeloid leukemia is still a fatal disease6,7.
Over the past decade using Massively Parallel Sequencing (MPS) or Next Generation Sequencing (NGS) the multi-dimensional genomic view with higher resolution of molecular profiling has become possible8. It was shown that the combination of exome, genome and RNA Sequencing data9, called integrative sequencing10,11,12 paves the way to better understand the mechanisms and pathogenesis of the cancer and can potentially improve the clinical management and finding new therapeutic targets specially in the advanced phase patients13.
In the present study, we conducted integrated genomic sequencing through an in-house filtering algorithm to uncover the leukemogenic important/actionable findings and investigate the myeloid leukemogenesis model in three Iranian myeloid blast crisis CML patients.
Three patients (from 3 unrelated families) with the clinical diagnosis of blastic phase CML were recruited from the Shariati Hospital in Tehran, Iran. Inclusion criteria were based on clinical presentations and hematological findings in the peripheral blood and bone marrow, according to the WHO 2016-guideline14.
The study was approved by the Ethical Committee of the Iran University of Medical Sciences (Code: IR.IUMS.REC 1395.95-04-87-30235). All subjects and/or their legal guardian(s) gave their informed consent to participate. All methods, including obtaining informed consent, were carried out in compliance with the above-mentioned ethical standards and guidelines.
Whole Exome Sequencing (WES) along with Chromosome-seq and RNA Sequencing were performed on extracted DNA and RNA from peripheral blood samples in the blast crisis phase.
Peripheral blood samples of patients were collected in EDTA containing tubes and used to extract genomic DNA by Blood SV-mini kit (GeneAll Biotechnology Co., LTD, South Korea) according to the manufacturer instruction. The extracted genomic DNA samples were subjected to WES using Agilent SureSelect V6-post Capture/Target Enrichment Kit. For the third patient TWIST Human Core Exome Kit (Twist Bioscience, USA) was used according to the manufacturer’s instruction. The enriched libraries were sequenced on Illumina platforms (Illumina Inc., CA, USA). Table S1 shows a summary report of data analysis metrics from the patients’ WES.
Paired-end sequence reads were mapped to the refseq (NCBI) human reference genome (GRCh38 assembly) by BWA-MEM. Picard software was used to mark duplicate reads. Genome Analysis Tool Kit (GATK) software Mutect2 package was used to call single nucleotide variants (SNVs) and short insertions or deletions (Indels). The intronic variants were removed, however the 100 bp flanking each exon was included in VCF files using interval/bed file.
Several tools were used to annotate and interpret variants including CGI, Varsome, InterVar, CancerVar, and Franklin. The annotated VCF files were filtered to find Tier1 or Tier2 variants according to AMP classification and pathogenic (P), likely pathogenic (LP), or Leaning pathogenic variants based on ACMG classification.
All Tier 1 or Tier 2 variants in AMP classification and the Tier 3 variants that were pathogenic (P), likely pathogenic (LP), or Leaning pathogenic VUS in ACMG classification were selected as potentially important findings (PIFs). PIFs known to be therapeutic, prognostic or diagnostic biomarkers were annotated as potentially actionable findings (PAFs).
RNA from whole blood in PaxGene tube was extracted and processed by the TruSeq Stranded Total RNA with Ribo-Zero Gold Kit (Illumina Inc., CA, USA) using 1 µg of DNAse treated RNA. Then paired-end sequencing was performed on Illumina platforms.
CLC genomics workbench 20 was used for RNA sequencing and fusion gene analyses. RNA sequencing results were used for determining the outlier expression analysis of genes harboring important/actionable variants, analyzing the fusion transcripts as well as confirmation of variants found in WES.
Vista™ Chromosome Sequencing-100 K an NGS based Whole Genome Sequencing was performed in BGI Clinical Laboratories to evaluate copy number variations (≥ 100 Kb).
The identified SNVs and small Indels in WES were visually validated in the mapping read-tracks after importing WES BAM files in the CLC genomics workbench 20. Besides, the validations of identified WES variants were performed using RNASeq BAM files by confirming both variant and gene expression.
The detected CNVs in chromosome-seq data were validated by CNV analysis from WES data using CLC genomics workbench 20.
The study was approved by the Ethical Committee of the Iran University of Medical Sciences (Code: IR.IUMS.REC 1395.95-04-87-30235). Informed consent was obtained from all subjects and/or their legal guardian(s) for participation.
The patient (I), was a 66-year-old female presented by fever, hepatosplenomegaly, leukocytosis, anemia, thrombocytopenia. The blast level in the peripheral blood was 25%. Bone marrow aspiration and immunophenotyping confirmed the myeloid blastic crisis. The patient had a history of thyroid cancer that had been treated many years prior to the development of CML. She was taken to the hospital at the age of 63 with a primary diagnosis of chronic phase CML in terms of leukocytosis with myeloid precursors in the peripheral blood and a positive Philadelphia chromosome in karyotype.
The patient (II), a 55-year old female presented by leukocytosis, anemia, thrombocytopenia. The diagnosis of chronic phase CML was established 15 years before blastic phase at the age of 40 years. There was 36% blast in the peripheral blood. Bone marrow aspiration and immunophenotyping confirmed the myeloid blastic crisis. At blast crisis, the peripheral blood karyotype revealed a clone that was positive for both the philadelphia chromosome and isochromosome 17q, as well as a clone with paracentric inversion in 15q.
Patient (III) was a 45-year old male presented by leukocytosis, anemia and thrombocytopenia. The diagnosis of chronic phase CML was established two years before blastic phase at the age of 43 years. There was 40% blast in the peripheral blood. Bone marrow aspiration and immunophenotyping confirmed the myeloid blastic crisis.
BCR-ABL1 fusion transcript (b13-a2) was also found in the RNA-seq fusion gene study of all patients during the blast crisis phase.
Table 1 summarizes the patients' detailed clinical and paraclinical findings at chronic phase CML diagnosis. The clinical and paraclinical findings at myeloid blast crisis are presented in the Table 2.
WES analysis identified sixteen PIFs (Figs. 1, 2 and 3) in which ABL1-Y272H, ASXL1-A627Gfs*8, GATA2-L359S, ABL1-F359V, ASXL1-G646Wfs*12, PTPN11-G503V and IKZF1-R208Q variants were classified as PAFs. The ABL1-I418T, PIK3CB-R231H, CHEK2:c.-4C > T, BRCA2-P3292L, ABL1-E459G, PRPF8-G1796R, PHF6-R274Q, WT1-R435X and ATM-C117Y variants were potentially important cancer variants which were not actionable (Table 3).
The paired-end mapping read-track displaying reads and detected variants in the first patient. The potentially important findings from WES data analysis is presented in the CLC genomics workbench software. (A) hg38 sequence track with WES mapped reads track; (B) hg38 sequence track with RNA-Seq mapped reads track; (C) Gene expression track; (D) gene, mRNA and CDS tracks; (E) the variant genomic coordinate on hg38 and the variant name. P1: patient No.1.
The paired-end mapping read-track displaying reads and detected variants in the second patient. The potentially important findings from WES data analysis is presented in the CLC genomics workbench software. (A) hg38 sequence track with WES mapped reads track; (B) hg38 sequence track with RNA-Seq mapped reads track; (C) Gene expression track; (D) gene, mRNA and CDS tracks; (E) the variant genomic coordinate on hg38 and the variant name. P2: patient No.2.
The paired-end mapping read-track displaying reads and detected variants in the third patient. The potentially important findings from WES data analysis is presented in the CLC genomics workbench software. (A) hg38 sequence track with WES mapped reads track; (B) hg38 sequence track with RNA-Seq mapped reads track; (C) Gene expression track; (D) gene, mRNA and CDS tracks; (E) the variant genomic coordinate on hg38 and the variant name. P3: patient No.3.
The patient was a BCR-ABL1 positive CML case in the myeloid blast crisis phase. She had abnormal karyotype with Philadelphia chromosome. In chromosome-seq there were no specific copy number alterations.
In WES she had six PIFs which contributed in all classes of leukemogenic genes (Fig. 1) except splicing components (class V). The p.Y272H variant on ABL1 (class I), p.A627Gfs*8 on ASXL1 (class III) and p.L359S on GATA2 (class II) were PAFs.
The patient was a BCR-ABL1 positive CML case in the myeloid blast crisis phase. She had abnormal complex karyotype with Philadelphia chromosome and isochromosome 17q in one clone and paracentric inversion in 15q in the another stemline (46,XX, t(9;22)(q34;q11.2),i(17)(q10)[17] /46,sl,inv(15)(q21,q25)[33]). In chromosome-seq (Fig. S1) there were deletions in 15q and 17p, and duplication in 17q (46,XX,del(15q26.3); del(17p11.2p13.3); dup(17p11.2q25.3)). The TP53 gene is a class II leukemogenic gene (transcription factor) located in 17p13.1 which is deleted in this patient.
In WES she had six PIFs that contribute in all classes of leukemogenic genes (Fig. 2). Considering the TP53 involvement by 17p deletion, this patient has seven PIFs. The p.F359V variant in ABL1 (class I), p.G646Wfs*12 in ASXL1 (class III) and TP53 (class IV) deletion were PAFs.
The patient was a BCR-ABL1 positive myeloid blast crisis CML. In chromosome-seq (Fig. S1) he had several deletions in 7p and 7q (46,XY,del(7p14.3p21.1); del(7p11.2p14.3); del(7q11.23q21.11)). The IKZF1 (7p12.2), EZH2 (7q35-q36) and CUX1 (7q22) are important leukemogenic genes which are located in the deletion regions. EZH2 gene is an epigenetic regulator gene (H3K27 methyltransferase; class III) and CUX1 is a transcription factor (class II) regulating TP53 and ATM.
He also had four PIFs in WES (Fig. 3). In terms of the pathogenic deletions in the short and long arms of chromosome 7 affecting EZH2 and CUX1, this patient has six PIFs which involve four classes of leukemogenic genes. PTPN11:p.G507V (class I) and IKZF1:p.R65Q variants (class II) were PAFs. The WT1:p.R435X (class II) and ATM:p.C117Y (class IV) variants were not PAFs.
IKZF1 variant (p.R65Q) in combination with del(7p11.2p14.3) may affect IKZF1 gene product as an important transcription factor in myeloid malignancies.
The ability to analyze and model the leukemogenesis process may aid in the effective management and treatment of hematologic malignancies. According to the proposed “Slot machine” model of leukemogenesis by Murati et al.4, progression from chronic phase into secondary acute myeloid leukemia (blast crisis) can be modeled by the combination of at least four steps that one classes of leukemogenic genes is affected in each step.
Using integrated genomic sequencing to assess the genes affected by SNVs, Indels and CNVs, in the present study all MBC-CML patients had at least five leukemogenic genes involved. The detected SNV and Indel variants in WES represent all classes of leukemogenic genes (Fig. 4) including signaling pathway gene (ABL1-Y272H, ABL1-I418T, ABL1-F359V, ABL1-E459G, PIK3CB-R231H, PTPN11-G503V), transcription factors (GATA2-L359S, PHF6-R274Q, IKZF1-R208Q, WT1-R435X), epigenetic regulators (ASXL1-A627Gfs*8, ASXL1-G646Wfs*12), tumor suppressor and DNA repair genes (BRCA2-P3292L, ATM-C117Y, CHEK2:c.-4C > T) and components of the spliceosome (PRPF8-G1796R).
The detected SNV and Indel variants in WES related to leukemogenic gene classes. The classes of leukemogenic genes includ signaling pathway genes (class I), transcription factors (class II), epigenetic regulators (class III), tumor suppressor and DNA repair genes (class IV) and components of the spliceosome (class V).
The detected variants in different signaling pathway related genes including ABL1, PIK3CB and PTPN11 in MBC-CML patients may affect ErbB, PI3K-Akt, Ras-MAPK cascade and JAK-STAT signaling pathways15,16,17. ABL1 was involved in the first and second patients. Four important missense variants were detected in ABL1 (p.Y272H, p.I418T, p.F359V, p.E459G). All variants were located in the tyrosine kinase domain of ABL1 and were previously reported in CML (COSM12576, COSM12605, COSM1460549)18. Mutations in the ABL1 kinase domain, in particular, have been identified in CML patients in the blastic phase19. Occurrence of 2 or more mutations in BCR–ABL1 fusion gene following TKI therapy has been reported in CML progression20.
PIK3CB-R231H is an important missense variant in the Ras Binding domain (RB domain). This deleterious variant in PIK3CB gene has not been reported in the hematological malignancies including CML (https://www.mycancergenome.org/content/gene/pik3cb/). This gene plays a critical role in solid tumor development16,21.
The third patient’s PTPN11-G507V variant is located in the Phosphatase (PH) domain. Although PTPN11 mutations have been documented in blastic phase CML, the p.G507V alteration is a novel variant in CML and has only been identified in other myeloid malignancies such as JMML, AML and MDS. (COSM14271).
There were several important variants in different transcription factor genes including GATA2, PHF6, WT1 and IKZF1 in MBC-CML patients. In addition, TP53, IKZF1 and CUX1 genes were affected by chromosomal deletions/CNVs. GATA2-L359S variant in the first patient is located on the zinc finger GATA domain. It has been indicated that GATA2 germline mutations strongly predispose patients to leukemia22. GATA2 mutations were described in MBC-CML; however p.L359S variant has not been previously reported in blastic CML. p.L359S variant was suggested as a predisposing germline likely pathogenic variant in MDS/AML23.
PHF6-R274Q variant in patient No.2 is located in the PHD-like zinc binding domain. Mutations in PHF6 also have been reported in the blastic phase CML while p.R274Q variant is a novel mutation in CML which only reported in Early T-cell precursor (ETP) acute lymphoblastic leukemia (COSM306061).
Biallelic involvement of IKZF1 was observed in the third patient. The detected missense variant was located in the C2H2 zinc finger domain of the IKZF1 (p.R65Q). Because of the 7p deletion, this variant showed a high variant allele fraction (VAF = 100%) in WES analysis. IKZF1 mutations had been found in both chronic and blastic phase CML24.The patient with IKZF1 mutation also had a mutation in PTPN11. It has been shown that in AML patients, mutations in IKZF1 and PTPN11 are associated with aggressive clinical course and primary resistant to chemotherapy25. CUX1 gene involvement was detected in the third patient in the cytogenetic studies in terms of 7q deletion. CUX1 is a transcription factor regulating TP53 and ATM26.
17p deletion in TP53 gene was established by cytogenetic analysis in the second patient. TP53 is a transcription factor located in 17p13.1. TP53 regulate cell cycle, DNA repair and apoptosis. It was observed that 45% of blast phase MPN (MPN-BP) patients have a TP53-related defect such as TP53 gene mutations or haploinsufficiency27. Involvement of TP53 has been described in both chronic and blastic phases of CML28,29.
Additional chromosomal abnormalities (ACAs) have been discovered in the second and the third patients. Patients with high-risk ACA, such as i(17q) and -7/7q-, are known to have a poor response to TKIs and a greater risk of disease progression. -7/7q- are less common in BC than i(17q), although they have a greater unfavorable influence on prognosis30. The patient number 3 with multiple deletions in the chromosome 7 has shorter time to MBC.
WT1-R435X is another significant nonsense variant discovered in the third patient. WT1 is a tumor suppressor gene which has known function as a transcriptional regulator31. This mutation causes loss of function in WT1 by the activation of NMD (Nonsense mediated decay) mechanism. Mutations in WT1 were reported in the blastic phase CML32, Wilms tumor and desmoplastic small round cell tumor (COSM21401). Pathogenic variants in this gene are correlated with poor prognosis in the patients with myelodysplastic syndromes33. The p.R435X germline mutation in WT1 has been classified as pathogenic which is associated to Wilms tumor in the clinvar (RCV000003671).
Two frameshift single nucleotide deletions were identified in ASXL1 (A627Gfs*8, G646Wfs*12). ASXL1 mutations are more frequent in secondary AML and may contribute to the development of CML to AML. Moreover, ASXL1 mutation might be associated with an aggressive phenotype in myeloid malignancies4.
It has been shown that mutations in ABL1 kinase domain, ASXL1, IKZF1, TP53 are common in CML patients who develop BC12,30,34. ASXL1 mutation is the most frequent mutation in poor outcome patients. In addition, individuals with an ASXL1 mutation had a longer median time to develop BC12. In the current study, the first and the second patients with ASXL1 mutation had longer time to MBC than the third patient.
EZH2 gene deletion was established by chromosome-seq in the third patient. EZH2 gene is an epigenetic regulator gene (H3K27 methyltransferase). Loss-of-function mutations in EZH2 have been reported in primary myelofibrosis (PMF), MDS and MDS/MPN overlap syndromes27.
Several important variants have been documented in the tumor suppressor genes specially the genes involved in the DNA damage-recognition and DNA repair pathways including CHEK2, ATM and BRCA2. The CHEK2 c.-4C > T variant is an important Tier3 variant. CHEK2 has a role in DNA damage-recognition and repair, cell cycle and P53 signaling pathway. The c.-4C > T variant is located in the kozak sequence which may affect translation35. This variant has been reported in clinvar as VUS for hereditary cancer-predisposing syndrome (RCV000131287).
ATM-C117Y variant is a Tier2 missense variant. ATM is involved in homologous recombination, cell cycle and P53 signaling pathway36. ATM has a key role in the DNA damage-response pathway and loss of function mutations in ATM may promote acceleration of the blast crisis in CML37. The ATM-CHEK2-p53 axis has been shown to play an important role against cancer initiation by inducing apoptosis, cell cycle arrest or senescence in the cancer cells27.
BRCA2-P3292L variant is a missense Tier2 variant. BRCA2 is a tumor suppressor gene which contributes to homologous recombination38. It has been indicated that 20% of the patients with Therapy-Related Myeloid Neoplasms (t-MN) among breast cancer survivors have a germline mutation in the BRCA1/2, CHEK2, TP53, or PALB2 genes, which are important in DNA repair pathways27.
PRPF8 was the only RNA splicing factor gene with a significant variant. It is a core component of U2-type and U12-type spliceosome complexes39. PRPF8-G1796R variant in the second patient is a missense Tier3 variant. This variant has been classified as driver mutation and has not been previously reported.
The self-renewal, differentiation and proliferation are three main processes which are impaired in malignancies. In the current study, the integrated genomic sequencing detected important variants in twelve different cancer genes. Mutations in signaling genes may induce proliferation of malignant cells while mutations in the epigenetic regulators (ERs) or tumor suppressor genes (TSG) may promote self-renewal in the progenitor cells and transcription factors mutations may affect differentiation. Mutations in the spliceosome component may also involve differentiation and self-renewal of progenitor cells40. Therefore, the alterations in ABL1, PTPN11 and PIK3CB signaling molecules may potentiate the proliferation process in MBC-CML. Variants in the CHEK2, BRCA2 and ATM tumor suppressor genes, as well as the ASXL1 and EZH2 Epigenetic Regulators, may support the self-renewal of the proliferating progenitors. Variants in the GATA2, PHF6, WT1, IKZF1 and CUX1 transcription factors may operate as differentiation blockers. Finally, the variant in the PRPF8 as spliceosome components may act as differentiation blocker or may cause self-renewal of the proliferating progenitor.
Finally, a recent integrated model of blast crisis41 showed that one of the most important mechanism associated to the BC progression is the epigenetic reprogramming related to the abnormal Polycomb Repressive Complex (PRC) function. The components of the PRC complex as an epigenetic complex, repress the expression of certain genes including leukemogenic HOXA cluster. On the other hand, specific gene products interact with PRC components to add/remove repressive marks or recruitment of PRC components to distinct loci. Based on recent studies on BC genome and transcriptome, it has been shown that the components of the PRC complex have mutated or altered in expression respectively. These changes overally affect the gene expression signature of BC cells. ASXL1 gene products interact with PRC2 components including EZH2. IKZF1 also recruits PRC2 to the target gene loci. According to this model, in the first and second patients the frameshift mutations in the ASXL1 may have a critical role in BC progression. In the third patient the large deletions in the chromosome 7 which affect IKZF1 (7p12.2) and EZH2 (7q35-q36) may accelerate BC.
A suggested model of leukemogenesis for the hierarchy of events is presented in the Fig. 5. In this model the initial step is the presence of mildly-deleterious variants (VUS-P) mainly in the tumor suppressor genes (TSG)/DNA damage repair (DDR) genes, transcription factors (TF) and/or components of splisosome. These variants may increase the individual’s susceptibility and predisposition to malignancies. The BCR-ABL1 fusion gene and mutations in the epigenetic regulator (ER)/transcription factor (TF) genes, which mainly influence the function of PRC, are the key genetic discoveries in the second and third steps. The final step in the accelerated and blast crisis phase may be started by the mutations in signaling molecules and/or additional chromosomal abnormalities (ACAs). ABL1 mutations specifically in the kinase domain lead to TKI resistance and delay in TKI switching may affect prognosis and survival by blast crisis. ACAs which involve leukemogenic genes (such as ER, TF or TSG) also may exacerbate the clinical prognosis and enter the disease to the blastic phase.
A suggested model of leukemogenesis. In this model the initial step is the presence of mildly-deleterious variants which may affect the tumor suppressor genes (TSG), transcription factors (TF) and/or components of splisosome and increase susceptibility to cancer. In the second and third steps, the major genetic findings are BCR-ABL1 fusion gene and mutation in the epigenetic regulator (ER)/transcription factor (TF) genes. The final step in the accelerated and blast crisis phase may be started by the mutations in the signaling molecules and/or additional chromosomal abnormalities (ACAs).
Although mutations in ASXL1 or IKZF1 genes may play a key role in BC, they are more often found at the time of diagnosis in the chronic phase, in individuals with poor prognosis. Therefore, mutations in oncogenes such as ABL1 and PTPN11 may act as an accelerator.
In this study, samples related to the different phases of disease progression were not available. Future studies using genomic sequencing in different serial time points may allow detecting the sequence of occurrences related to the disease progression. Furthermore, recent advances in single cell sequencing may help to more accurately clarify the disease mechanisms.
The integrated genomic sequencing in this study effectively identified large spectrum of important and actionable variants from SNVs to CNVs. Mildly-deleterious passenger variants mainly in DNA damage repair (DDR) genes may increase the individual’s predisposition to cancer. BCR-ABL1 fusion, mutations in signaling and epigenetic regulator genes along with ACAs may act as cancer drivers. Epigenetic reprogramming caused by ASXL1 or IKZF1 mutations seems to be one of the most critical mechanisms linked with BC development. ABL1 domain mutations, mutations in other signaling genes and high risk ACAs may promote CML to BC.
Using this strategy we have identified several actionable findings which may improve targeted therapies. These findings may help better understanding the underlying pathophysiological mechanisms and may provide further insights into MBC-CML genetic basis in future. Further studies will shed light on the clinical usefulness of the integrated genomic analysis.
The datasets generated and/or analyzed during the current study are available in the ENA https://www.ebi.ac.uk/ena/browser/view/PRJEB51084.
A Correction to this paper has been published: https://doi.org/10.1038/s41598-023-27872-1
Accelerated phase
Depth of coverage
American college of medical genetics and genomics
American society of clinical oncology
Association for molecular pathology
Benign
Bone marrow aspiration
Chronic myeloid leukemia
Cancer genome interpreter
College of American pathologists
Copy number variation
Diagnosis
Driver
Epigenetic regulators
Genome analysis tool kit
Known driver
Leukemia progenitor cell
Leukemia stem cell
Likely benign
Likely pathogenic
Myelodysplastic syndrome
Myeloid blast crisis
Massively parallel sequencing
Myeloproliferative neoplasms
Next generation sequencing
Not protein affecting
Passenger
Pathogenic
Peripheral blood
Potentially important finding
Potentially actionable finding
Reverse transcriptase-polymerase chain reaction
Single nucleotide variant
Transcription factor
Tumor suppressor gene
Variant of uncertain significance
Variant of uncertain significance leaning benign
Variant of uncertain significance leaning pathogenic
Variant allele fraction
Whole exome sequencing
Houshmand, M. et al. Bone marrow microenvironment: The guardian of leukemia stem cells. World J. Stem Cells 11(8), 476–490 (2019).
Article Google Scholar
Lebecque, B., Bourgne, C., Vidal, V. & Berger, M. G. DNA methylation and intra-clonal heterogeneity: The chronic myeloid leukemia model. Cancers (Basel) 13(14), 3587 (2021).
Article CAS Google Scholar
Valent, P. Imatinib-resistant chronic myeloid leukemia (CML): Current concepts on pathogenesis and new emerging pharmacologic approaches. Biologics 1(4), 433–448 (2007).
CAS Google Scholar
Murati, A. et al. Myeloid malignancies: Mutations, models and management. BMC Cancer 12, 304 (2012).
Article CAS Google Scholar
Flis, S. & Chojnacki, T. Chronic myelogenous leukemia, a still unsolved problem: Pitfalls and new therapeutic possibilities. Drug Des. Dev. Ther. 13, 825–843 (2019).
Article CAS Google Scholar
Apperley, J. F. Chronic myeloid leukaemia. Lancet 385(9976), 1447–1459 (2015).
Article Google Scholar
Bonifacio, M., Stagno, F., Scaffidi, L., Krampera, M. & Di Raimondo, F. Management of chronic myeloid leukemia in advanced phase. Front. Oncol. 9, 1132 (2019).
Article Google Scholar
Avila, M. & Meric-Bernstam, F. Next-generation sequencing for the general cancer patient. Clin. Adv. Hematol. Oncol. 17(8), 447–454 (2019).
Google Scholar
Chudasama, P. et al. Integrative genomic and transcriptomic analysis of leiomyosarcoma. Nat. Commun. 9(1), 144 (2018).
Article ADS Google Scholar
Robinson, D. R. et al. Integrative clinical genomics of metastatic cancer. Nature 548(7667), 297–303 (2017).
Article ADS CAS Google Scholar
Kumar-Sinha, C. & Chinnaiyan, A. M. Precision oncology in the age of integrative genomics. Nat. Biotechnol. 36(1), 46–60 (2018).
Article CAS Google Scholar
Branford, S. et al. Integrative genomic analysis reveals cancer-associated mutations at diagnosis of CML in patients with high-risk disease. Blood 132(9), 948–961 (2018).
Article CAS Google Scholar
Beaubier, N. et al. Integrated genomic profiling expands clinical options for patients with cancer. Nat. Biotechnol. 37(11), 1351–1360 (2019).
Article CAS Google Scholar
Arber, D. A. et al. The 2016 revision to the World Health Organization classification of myeloid neoplasms and acute leukemia. Blood 127(20), 2391–2405 (2016).
Article CAS Google Scholar
Bollu, L. R., Mazumdar, A., Savage, M. I. & Brown, P. H. Molecular pathways: Targeting protein tyrosine phosphatases in cancer. Clin. Cancer Res. 23(9), 2136–2142 (2017).
Article CAS Google Scholar
Yang, J. et al. Targeting PI3K in cancer: mechanisms and advances in clinical trials. Mol. Cancer 18(1), 26 (2019).
Article Google Scholar
Hantschel, O. Structure, regulation, signaling, and targeting of abl kinases in cancer. Genes Cancer 3(5–6), 436–446 (2012).
Article Google Scholar
Jones, D. et al. Uncommon BCR-ABL kinase domain mutations in kinase inhibitor-resistant chronic myelogenous leukemia and Ph+ acute lymphoblastic leukemia show high rates of regression, suggesting weak selective effects. Blood 115(26), 5428–5429 (2010).
Article CAS Google Scholar
Soverini, S. et al. Next-generation sequencing for BCR-ABL1 kinase domain mutation testing in patients with chronic myeloid leukemia: A position paper. J. Hematol. Oncol. 12(1), 131 (2019).
Article Google Scholar
Hochhaus, A. et al. Chronic myeloid leukaemia: ESMO clinical practice guidelines for diagnosis, treatment and follow-up. Ann. Oncol. 28(suppl_4), iv41–iv51 (2017).
Article CAS Google Scholar
Whale, A. D., Colman, L., Lensun, L., Rogers, H. L. & Shuttleworth, S. J. Functional characterization of a novel somatic oncogenic mutation of PIK3CB. Sign. Transduct. Target. Ther. 2, 17063 (2017).
Article Google Scholar
Donadieu, J. et al. Natural history of GATA2 deficiency in a survey of 79 French and Belgian patients. Haematologica 103(8), 1278–1287 (2018).
Article CAS Google Scholar
Gao, J. et al. Heritable GATA2 mutations associated with familial AML-MDS: A case report and review of literature. J. Hematol. Oncol. 7, 36 (2014).
Article Google Scholar
Adnan-Awad, S., Kankainen, M. & Mustjoki, S. Mutational landscape of chronic myeloid leukemia: More than a single oncogene leukemia. Leuk. Lymphoma 62(9), 2064–2078 (2021).
Article CAS Google Scholar
Zhang, X. et al. The specific distribution pattern of IKZF1 mutation in acute myeloid leukemia. J. Hematol. Oncol. 13(1), 140 (2020).
Article Google Scholar
Vadnais, C. et al. CUX1 transcription factor is required for optimal ATM/ATR-mediated responses to DNA damage. Nucleic Acids Res. 40(10), 4483–4495 (2012).
Article CAS Google Scholar
Higgins, A. & Shah, M. V. Genetic and genomic landscape of secondary and therapy-related acute myeloid leukemia. Genes (Basel) 11(7), 749 (2020).
Article CAS Google Scholar
Seifert, H. et al. The prognostic impact of 17p (p53) deletion in 2272 adults with acute myeloid leukemia. Leukemia 23(4), 656–663 (2009).
Article CAS Google Scholar
Nakai, H. & Misawa, S. Chromosome 17 abnormalities and inactivation of the p53 gene in chronic myeloid leukemia and their prognostic significance. Leuk. Lymphoma 19(3–4), 213–221 (1995).
Article CAS Google Scholar
Hehlmann, R. Chronic myeloid leukemia in 2020. Hemasphere 4(5), e468 (2020).
Article Google Scholar
Yang, L., Han, Y., Suarez Saiz, F. & Minden, M. D. A tumor suppressor and oncogene: The WT1 story. Leukemia 21(5), 868–876 (2007).
Article CAS Google Scholar
Perrotti, D., Jamieson, C., Goldman, J. & Skorski, T. Chronic myeloid leukemia: Mechanisms of blastic transformation. J. Clin. Invest. 120(7), 2254–2264 (2010).
Article CAS Google Scholar
Hosoya, N., Miyagawa, K., Mitani, K., Yazaki, Y. & Hirai, H. Mutation analysis of the WT1 gene in myelodysplastic syndromes. Jpn. J. Cancer Res. 89(8), 821–824 (1998).
Article CAS Google Scholar
Grossmann, V. et al. A deep-sequencing study of chronic myeloid leukemia patients in blast crisis (BC-CML) detects mutations in 76.9% of cases. Leukemia 25(3), 557–60 (2011).
Article CAS Google Scholar
Wolf, A. et al. Single base-pair substitutions at the translation initiation sites of human genes as a cause of inherited disease. Hum. Mutat. 32(10), 1137–1143 (2011).
Article CAS Google Scholar
Phan, L. M. & Rezaeian, A. H. ATM: main features, signaling pathways, and its diverse roles in DNA damage response, tumor suppression, and cancer development. Genes (Basel) 12(6), 845 (2021).
Article CAS Google Scholar
Takagi, M. et al. ATM-dependent DNA damage-response pathway as a determinant in chronic myelogenous leukemia. DNA Repair (Amst.) 12(7), 500–507 (2013).
Article CAS Google Scholar
Prakash, R., Zhang, Y., Feng, W. & Jasin, M. Homologous recombination and human health: The roles of BRCA1, BRCA2, and associated proteins. Cold Spring Harb. Perspect. Biol. 7(4), a016600 (2015).
Article Google Scholar
Fica, S. M., Oubridge, C., Wilkinson, M. E., Newman, A. J. & Nagai, K. A human postcatalytic spliceosome structure reveals essential roles of metazoan factors for exon ligation. Science 363(6428), 710–714 (2019).
Article ADS CAS Google Scholar
Aponte, P. M. & Caicedo, A. Stemness in cancer: Stem cells, cancer stem cells, and their microenvironment. Stem Cells Int. 2017, 5619472 (2017).
Article Google Scholar
Ko, T. K. et al. An integrative model of pathway convergence in genetically heterogeneous blast crisis chronic myeloid leukemia. Blood 135(26), 2337–2353 (2020).
Article Google Scholar
Download references
The authors declare no conflict of interest. Our sincere appreciation goes out to Professor Ghavamzadeh, Dr. Shahrbanoo Keyhanian, Dr. Abolghasem Kollaee, Dr. Zarangis Khanpour, Dr. Mohammad Reza Mahdavi and Ms. Hashemabadi who kindly made it possible for us to access the medical data of patients at CML-CP.
This study is part of a PhD thesis funded by the Iran university of medical sciences (Number: 95-04-87-30235). This study was partly supported by the Cellular and Molecular Research Center, Iran university of medical sciences (Number: 97-01-117-33428).
These authors contributed equally: Golnaz Ensieh Kazemi-Sefat and Mohammad Keramatipour.
Department of Molecular Medicine, Faculty of Advanced Technologies in Medicine, Iran University of Medical Sciences, Tehran, Iran
Golnaz Ensieh Kazemi-Sefat
Cellular and Molecular Research Center, Iran University of Medical Sciences, Tehran, Iran
Golnaz Ensieh Kazemi-Sefat & Kazem Mousavizadeh
Department of Medical Genetics, School of Medicine, Tehran University of Medical Sciences, Tehran, Iran
Mohammad Keramatipour
Hematology, Oncology and Stem Cell Transplantation Research Center, Shariati Hospital, Tehran University of Medical Sciences, Tehran, Iran
Mohammad Vaezi, Shahrbano Rostami, Marjan Yaghmaie & Bahram Chahardouli
Oncopathology Research Center, School of Medicine, Iran University of Medical Sciences, Tehran, Iran
Seyed Mohsen Razavi
Laboratory of Complex Biological Systems and Bioinformatics (CBB), Department of Bioinformatics, Institute of Biochemistry and Biophysics (IBB), University of Tehran, Tehran, Iran
Kaveh Kavousi
Department of Medical Genetics, Faculty of Medicine, Mashhad University of Medical Sciences, Mashhad, Iran
Amin Talebi
Department of Medical Genetics, School of Medicine, Iran University of Medical Sciences, Tehran, Iran
Saeed Talebi
Department of Pharmacology, School of Medicine, Iran University of Medical Sciences, Tehran, Iran
Kazem Mousavizadeh
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
G.E.K.S. made contribution to the conception and data analysis, drafted the work and wrote the manuscript. M.K made contribution to the conception and data analysis. K.M. have made contribution to the conception, supervised the project and encouraged the investigation. S.T. suggested the main conceptual ideas and proofed outline and critically reviewed and revised the article for important intellectual content. M.V and S.M.R had role in the clinical management of the patient and contributed to the clinical analysis. K.K contributed in the NGS data analysis. A.T provided critical feedback and contributed in the NGS data analysis. Sh.R, M.Y and B.Ch contributed in performing hematological, cytogenetic and immunophenotyping studies and related data analysis. All authors read and approved the final manuscript. Golnaz Ensieh Kazemi Sefat (G.E.K.S.) and Mohammad Keramatipour (M.K.) contributed equally as the first authors. Saeed Talebi (S.T) and Kazem Mousavizadeh (K.M) contributed equally as corresponding authors. All authors reviewed the manuscript.
Correspondence to Saeed Talebi or Kazem Mousavizadeh.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original online version of this Article was revised: The original version of this Article contained errors in the spelling of the authors Shahrbano Rostami and Marjan Yaghmaie which were incorrectly given as Shahrbanoo Rostami and Marjan Yaghmaei respectively.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Kazemi-Sefat, G.E., Keramatipour, M., Vaezi, M. et al. Integrated genomic sequencing in myeloid blast crisis chronic myeloid leukemia (MBC-CML), identified potentially important findings in the context of leukemogenesis model. Sci Rep 12, 12816 (2022). https://doi.org/10.1038/s41598-022-17232-w
Download citation
Received: 14 March 2022
Accepted: 21 July 2022
Published: 27 July 2022
DOI: https://doi.org/10.1038/s41598-022-17232-w
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.