- Case report
- Open Access
Genomic trajectory in leukemogenesis of myeloproliferative neoplasms: a case report
BMC Medical Genomics volume 14, Article number: 137 (2021)
We report a patient with Essential Thrombocythemia (ET), subsequently diagnosed with concurrent myeloid and lymphoid leukemia. Generally, the molecular mechanisms underlying leukemic transformation of Philadelphia-negative myeloproliferative neoplasms (Ph-MPN) are poorly understood. Risk of transformation to acute myelogenous leukemia (AML) is low; transformation to both AML and acute lymphoblastic leukemia (ALL) is extremely low. Genetic defects, including allele burden, order of mutation acquisition, clonal heterogeneity and epigenetic mechanisms are important contributors to disease acceleration.
A 78-year-old Caucasian female originally treated for stable ET, underwent disease acceleration and transition to myeloid sarcoma and B-cell ALL. Genomic reconstruction based on targeted sequencing revealed the presence of a large del(5q) in all three malignancies and somatic driver mutations: TET2, TP53, SF3B1, and ASXL1 at high allele frequency. We propose that the combination of genetic and molecular abnormalities led to hematopoietic stem cell (HSC) injury and disease progression through sub-clone branching. We hypothesize that ancestral reconstruction of genomic data is a useful tool to uncover subclonal events leading to transformation.
The use of ancestral reconstruction of genomic data sheds light on the unique clinical scenario described in this case report. By determining the mutational profile of tumors at several timepoints and deducing the most parsimonious relationship between them, we propose a reconstruction of their origin. We propose that blast progression originated from subclonal events with malignant potential, which coexisted with but did not originate from JAK2 p.V617F-positive ET. We conclude that the application of genomic reconstruction enhances our understanding of leukemogenesis by identifying the timing of molecular events, potentially leading to better chemotherapy choices as well as the development of new targeted therapies.
The molecular mechanisms underlying leukemic transformation of Philadelphia-negative myeloproliferative neoplasms (Ph-MPN) are poorly understood. Risk of transformation to acute myelogenous leukemia (AML) is approximately 10–20% for patients with primary myelofibrosis (PMF) and considerably lower for those with Polycythemia Vera (PV) and Essential Thrombocythemia (ET), 5–10% and 2–5% respectively [1, 2]. Transformation of MPN to acute lymphoblastic leukemia (ALL) is rare, with only seventeen cases reported in the literature . Retrospective studies suggest that genetic defects, including allele burden, order of mutation acquisition, clonal heterogeneity and epigenetic mechanisms play an important role in the observed conversion rate, as well as the varied clinical-pathologic entities that evolve [4, 5].
In this report, we describe a series of clonal events in an elderly patient originally diagnosed with JAK2 V617F-positive ET who presented in accelerated phase, with subsequent progression to concurrent myeloid sarcoma (MS) and B-cell ALL. We hypothesize that reconstruction of genomic mutations to build an ancestral tree may be a useful tool to characterize leukemogenesis in patients with high-risk disease.
A 78-year-old Caucasian female was diagnosed with JAK2 V617F-positive ET in 2010 based on laboratory data and molecular profiling of peripheral blood. She was initiated on hydroxyurea (HU) and low-dose aspirin. The patient remained on anticoagulation with coumadin due to a prior history of atrial fibrillation. She demonstrated stable and adequate platelet response to 10 mg/kg of HU.
In 2016, despite compliance with therapy, she presented with platelet count of 1150 k/µl, hemoglobin of 8.2 g/dL, and WBC of 6.0 k/µl with absolute neutrophil count of 2820/µl without detectable peripheral blasts. She underwent a bone marrow biopsy, which revealed 80–90% cellularity with increased myeloblasts (12%), megakaryocyte hyperplasia, dysmegakaryopoiesis and grade 2/4 reticulin fibrosis (Fig. 1a). Cytogenetics revealed del(5)(q22-q33),del(17)(p11.2)/46,XX, with wild-type CALR and MPL. Fluorescent in situ hybridization (FISH) reconfirmed 5q deletion. Next generation sequencing (NGS) revealed the following mutations: JAK2 p.V617F, SF3B1 p.R625C and TET2 p.C1271fs (Table 1). Due to patient preference for an oral agent, lenalidomide was initiated at 10 mg orally daily for 21 days in a 28-day cycle. This intervention resulted in partial remission based on International Working Group-Myeloproliferative Neoplasms Research and Treatment Response Criteria, with clinical improvement of anemia and thrombocytosis (hemoglobin of 12.1 g/dL and platelet count of 344 K/µL, respectively). After seven cycles of lenalidomide, a repeat bone marrow biopsy revealed persistent blasts at 8%, persistent dysmegakaryopoiesis, and stable JAK2 p.V617F mutation. Despite cytogenetic remission, (46 XX in 20 metaphases analyzed), her MDS FISH panel revealed a new 20q12 (PTPRT/MYBL2) loss.
In January 2018, she presented with shortness of breath, abdominal pain, nausea and vomiting. Computed tomography (CT) scan of chest, abdomen and pelvis revealed a large 7.9 cm × 4.9 cm ileocecal mass inducing partial colonic obstruction and moderate left pleural effusion. Pathology of the excised cecal mass showed a cluster of medium-sized MPO-positive myeloblasts with eccentric nuclei and scant cytoplasm, positive for CD33 and CD117 based on flow cytometry and consistent with MS (Fig. 1b). Macroscopic examination of her pleural fluid revealed medium-sized cells with large nuclei and scant basophilic cytoplasm positive for CD19, TdT and CD20 consistent with B-cell ALL (Fig. 1c). NGS of both the colonic mass and pleural fluid revealed SF3B1, TP53 exon 6 and exon 3 and TET2 (Table 1). JAK2-V617F and MPL were negative in both the colonic MS and B-cell ALL (Fig. 2).
Amplicon-based NGS of a panel of 48 genes commonly mutated in hematologic malignancies was performed on the Ion Torrent PGM platform from unfractionated bone marrow DNA.
An inferred phylogenetic tree was constructed from NGS and karyotype data from bone marrow samples obtained during the evolution from ET to MDS, as well as from colon and pleural effusion samples obtained during blast transformation. The corresponding time-scale data points were aligned to reflect the hypothetical evolution of hematopoietic stem cells (HSC) most parsimonious with the sequencing and FISH findings reflected in Table 1. Note that between any two splits of the tree, the order of mutations and chromosomal deletions cannot be determined, and mutations are depicted sequentially for visual convenience (Fig. 2).
Discussion and Conclusions
In this single case report, genomic reconstruction importantly revealed the presence of a large del(5q) in all three malignancies: ET, MS and B-cell ALL, suggesting a common cell of origin. In chronic phase MPN, this finding should inform physicians about the imminent likelihood for disease progression. Episomal reprogramming has identified del(5q) as an early cytogenetic lesion with the capacity to perturb genome stability and differentiation . Larger del(5q) size has been correlated with higher mutation frequency , which in this case included the somatic driver mutations: TET2, TP53, SF3B1, and ASXL1 at high allele frequency (Table 1). We propose that this combination led to HSC injury and disease progression through sub-clone branching (Fig. 2). The shared presence of JAK2-V617F and TET2 mutations in ET suggests that “megakaryocytic branching” originated directly from the HSC as opposed to a more lineage-restricted progenitor . The JAK2-V617F mutation confers a weak proliferative advantage to HSC and its absence in blast phases suggests MS and B-cell ALL did not emerge from this subclone.
Before blast conversion, the patient received lenalidomide, an immune modulator that yields cytogenetic remission by inhibiting growth of del(5q) progenitors without affecting other cells . Clinically, lenalidomide improves survival and reduces transfusion requirements in patients with del(5q) MDS [4, 8]. Some reports suggest an association between lenalidomide therapy and transformation to more aggressive phenotypes, including a review documenting transformation of MPN to ALL . Emerging data suggest that progression to AML in patients treated with lenalidomide is associated with karyotype complexity and clonal selection rather than a drug-mediated transformation [4,9]. Accordingly, in our case, bone marrow biopsies after lenalidomide therapy reveal absence of del(5q) clones, likely as a result of known suppression from drug. They also revealed a new 20q deletion (Table 1), which may portend the malignant potential of uninhibited clones, as evidenced later by its high expression in B-cell ALL.
Based on the ancestral reconstruction of genomic data modeled in Fig. 2, we hypothesize a mechanism for disease acceleration, whereby subclonal events with potential for blast conversion coexisted with but did not originate from JAK2 p.V617F-positive ET. The presence of del(5q) and molecular abnormalities (TP53, KMT2A) in both MS and B-cell ALL (Fig. 2) suggest that “disease progression” originated from diversification of a pluripotent HSC capable of both myeloid and lymphoid differentiation [10, 11], thereby leading to two distinct leukemia initiating cells (LIC): one containing 17p deletion and the other 20q deletion. Interestingly, del(5q) was suppressed in bone marrow but not from extramedullary sites where the blast phase manifested, suggesting variable sensitivity of different clones to lenalidomide. The concurrent presence of TP53 mutation in extramedullary sites  is known to confer a negative impact on survival and drug response to patients with del(5q) MDS treated with lenalidomide [7, 13]. Recent studies suggest that patients with high-risk MDS, characterized by unfavorable-risk cytogenetic abnormalities and/or TP53 mutations, exhibit favorable clinical responses with robust mutation clearance when treated with hypomethylating agents (HMA) . Unfortunately, HMA do not provide durable responses. A combination of azacytidine plus anti-CD47 monoclonal antibody is currently being investigated on TP53-AML with preliminary results showing an objective response in 71% of subjects and 48% complete remission .
We are limited in our ability to confirm the proposed order of pathogenic mutations. We lack banked bone marrow cells at all time points of disease evolution to demonstrate the proposed patterns of clonal progression. As this is a single case report, computer simulation of similar cases of accelerated phase ET would reaffirm our proposed model. However, collecting replicate cases is difficult given the rarity of the events described. Despite the limitations in our analysis, the ancestral tree in this case report highlights how the relative accessibility of NGS continues to improve our understanding of leukemogenesis, specifically, the predictive significance of large del(5q). It also has the capacity to inform therapeutic choices. Notwithstanding the presence of del (5q), recent studies support that HMA is a superior choice to lenalidomide under the clinical scenario described here. HMA combination regimens  currently being investigated may provide durable responses to patients with TP53 mutations.
In conclusion, the advances in NGS technology have made it possible to generate a deep snapshot of the genetic composition of rapidly evolving tumor cells. The development and application of computational algorithms to harness NGS data for ancestral reconstruction could have a dramatic impact in how we practice medicine. Large scale simulation of similar cases to that described in this report would enhance our ability to predict disease progression. It would inform how a specific collection of molecular mutations would favor one therapeutic approach versus another. We propose exploring the use of computational algorithms to resolve patterns of clonal progression to enhance our understanding of leukemogenesis and the contribution of molecular targets, thereby leading to the development of more specific, tumor-tailored therapeutic interventions.
Availability of data and materials
The datasets generated and analyzed during the current study are available in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) repository, accession number PRJNA727835 (http://www.ncbi.nlm.nih.gov/bioproject/727835). Mutational analyses by targeted sequencing panels are provided in the data table.
Acute lymphoblastic leukemia
Acute myelogenous leukemia
Fluorescent in situ hybridization
Hematopoietic stem cell
Leukemia initiating cell
Next generation sequencing
Philadelphia-negative myeloproliferative neoplasm
Rampal R, Ahn J, Abdel-Wahab O, et al. Genomic and functional analysis of leukemic transformation of myeloproliferative neoplasms. Proc Natl Acad Sci U S A. 2014;111:E5401–10.
Iurlo A, Cattaneo D, Gianelli U. Blast transformation in myeloproliferative neoplasms: risk factors, biological findings, and targeted therapeutic options. Int J Mol Sci 2019;20:https://doi.org/10.3390/ijms20081839.
Alhuraiji A, Naqvi K, Huh YO, Ho C, Verstovsek S, Bose P. Acute lymphoblastic leukemia secondary to myeloproliferative neoplasms or after lenalidomide exposure. Clin Case Rep. 2017;6:155–61.
Giagounidis A, Mufti GJ, Fenaux P, Germing U, List A, MacBeth KJ. Lenalidomide as a disease-modifying agent in patients with del(5q) myelodysplastic syndromes: linking mechanism of action to clinical outcomes. Ann Hematol. 2014;93:1–11.
Stengel A, Kern W, Haferlach T, Meggendorfer M, Haferlach C. The 5q deletion size in myeloid malignancies is correlated to additional chromosomal aberrations and to TP53 mutations. Genes Chromosomes Cancer. 2016;55:777–85.
Hsu J, Reilly A, Hayes BJ, et al. Reprogramming identifies functionally distinct stages of clonal evolution in myelodysplastic syndromes. Blood. 2019;134:186–98.
Noetzli LJ, French SL, Machlus KR. New Insights Into the Differentiation of Megakaryocytes From Hematopoietic Progenitors. Arterioscler Thromb Vasc Biol. 2019;39:1288–300.
Lee JH, List A, Sallman DA. Molecular pathogenesis of myelodysplastic syndromes with deletion 5q. Eur J Haematol. 2019;102:203–9.
Pellagatti A, Jadersten M, Forsblom AM, et al. Lenalidomide inhibits the malignant clone and up-regulates the SPARC gene mapping to the commonly deleted region in 5q- syndrome patients. Proc Natl Acad Sci U S A. 2007;104:11406–11.
Xie W, Chen Z, Wang SA, et al. Lymphoblastic leukemia following myelodysplastic syndromes or myelodysplastic/myeloproliferative neoplasms. Leuk Lymphoma. 2019;60:2993–3001.
Janssen JW, Buschle M, Layton M, et al. Clonal analysis of myelodysplastic syndromes: evidence of multipotent stem cell origin. Blood. 1989;73:248–54.
Tang G, DiNardo C, Zhang L, et al. MLL gene amplification in acute myeloid leukemia and myelodysplastic syndromes is associated with characteristic clinicopathological findings and TP53 gene mutation. Hum Pathol. 2015;46:65–73.
Mossner M, Jann JC, Nowak D, et al. Prevalence, clonal dynamics and clinical impact of TP53 mutations in patients with myelodysplastic syndrome with isolated deletion (5q) treated with lenalidomide: results from a prospective multicenter study of the german MDS study group (GMDS). Leukemia. 2016;30:1956–9.
Welch JS, Petti AA, Miller CA, et al. TP53 and Decitabine in Acute Myeloid Leukemia and Myelodysplastic Syndromes. N Engl J Med. 2016;375:2023–36.
Sallman D, Asch A, Kambhampati S, et al. The first-in-class anti-CD47 antibody magrolimab combined with azacitadine is well-tolerated and effective in AML patients: Phase 1b results. In: 62nd American Society of Hematology Meeting 2020;Session 613:Abstract 330 (oral).
The authors thank Cliff Bleiden, MD for preparing photomicrographs representing the histopathology of the patient’s disease.
KYK, MK and YC were supported in part by National Institute of Health (NIH) grants R01HL136333 and R01HL134880.
Ethics approval and consent to participate
Submitted case was approved by the institutional review board [Baylor College of Medicine] approval number H-43327 dated 20th, August 2018. This process is in accordance with the Helsinki declaration.
Consent for publication
Written informed consent was obtained from the patient for publication. A copy of the written consent is available for review.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
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/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Chen, Y., Talukder, R., Merritt, B.Y. et al. Genomic trajectory in leukemogenesis of myeloproliferative neoplasms: a case report. BMC Med Genomics 14, 137 (2021). https://doi.org/10.1186/s12920-021-00986-z
- Essential thrombocythemia
- Clonal evolution
- Myeloproliferative neoplasms
- Case report