eISSN: 1644-4124
ISSN: 1426-3912
Central European Journal of Immunology
Current issue Archive Manuscripts accepted About the journal Special Issues Editorial board Abstracting and indexing Subscription Contact Instructions for authors Publication charge Ethical standards and procedures
Editorial System
Submit your Manuscript
SCImago Journal & Country Rank
4/2021
vol. 46
 
Share:
Share:
Case report

New genetic variants of TET2 and ASXL1 identified by next generation sequencing and pyrosequencing in a patient with MDS-RS-MLD and secondary acute myeloid leukemia

Monika Małgorzata Adamska
1
,
Ewelina Kowal-Wiśniewska
1, 2
,
Katarzyna Kiwerska
2
,
Adam Ustaszewski
2
,
Joanna Czerwińska-Rybak
1
,
Zuzanna Kanduła
1
,
Marzena Wojtaszewska
1
,
Marta Barańska
1
,
Łukasz Pruchniewski
1
,
Krzysztof Lewandowski
1
,
Małgorzata Jarmuż-Szymczak
1, 2
,
Lidia Gil
1

  1. Department of Hematology and Bone Marrow Transplantation, Poznan University of Medical Sciences, Poland
  2. Institute of Human Genetics, Polish Academy of Sciences, Poznan, Poland
(Cent Eur J Immunol 2021; 46 (4): 524-530)
Online publish date: 2021/12/12
Article file
- New genetic variants.pdf  [0.49 MB]
Get citation
 
PlumX metrics:
 

Introduction

Myelodysplastic syndromes (MDS) represent a heterogeneous group of clonal hematological malignancies characterized by the presence of cytopenias resulting from ineffective hematopoiesis, multilineage dysplasia and high risk of transformation to secondary acute myeloid leukemia (secAML) [1, 2]. Whole genome sequencing identified more than 45 genes involved in MDS molecular pathophysiology, which encode proteins involved in pre-mRNA splicing (SRSF2, SF3B1, U2AF1, ZRSR2), DNA methylation (TET2, DNMT3A, IDH1/2), histone modification (ASXL1, EZH2), signal transduction and the TP53 pathway. Among them also cohesion components (SMC1A, RAD21, STAG1, STAG2, SMC3) and transcription factors (RUNX1, ETV6, GATA2) were found [1, 3]. Various somatic mutations were associated with overall prognosis and specific MDS morphological features [2]. Understanding molecular mechanisms of genetic disease evolution may allow for targeted MDS treatment. Measurable residual disease (MRD) and donor chimerism (DC) monitoring remain essential in the management strategy of MDS patients after allogenic hematopoietic cell transplantation (alloHCT). Next generation sequencing (NGS) application enables detection of sequence variants and their further monitoring. Therefore in this study we intended to evaluate the diagnostic utility of targeted NGS, Sanger sequencing (SSeq) and pyrosequencing in monitoring mutational dynamics with samples collected sequentially at different stages of MDS/secAML.

Material and methods

DNA was isolated from bone marrow samples (BM), peripheral blood samples (PB) and saliva samples (SAL) by a standard phenol : chloroform procedure. Samples collected in 2009-2018 from a patient with MDS/secAML were analyzed in order to examine the mutational profile during the course of the disease (Fig. 1). DNA samples isolated from PB of six healthy controls were used to establish the limit of allele frequency (AF) detection by pyrosequencing. Molecular analysis was performed in accordance with the Declaration of Helsinki and Ethics Committee approval (563/14) with written informed consent obtained from the patient. Both MDS and secAML diagnoses were established according to the 2016 World Health Organization criteria [2]. Treatment response and cytogenetic and genetic risk were evaluated according to the 2017 European LeukemiaNet (ELN) recommendations [4].

Fig. 1

A) Data concerning percentage content of bone marrow blasts, ring sideroblasts, peripheral blood chimerism, implemented treatment, performed molecular analysis and additional information during the course of the patient’s disease. B) Allele frequency results for selected DNA sequence variants (obtained by NGS and pyrosequencing) in BM, PB and SAL as well as information about AF for healthy control and AF cut-off value for selected DNA sequence variants. C) Electrophero- grams presenting results of pyrosequencing (c.1945G>T ASXL1; c.1874G>T SF3B1; c.4044+2dupT TET2; c.4076G>T TET2) performed in: I – May 2009, II – January 2017, III – May 2018

/f/fulltexts/CEJOI/45715/CEJI-46-45715-g001_min.jpg

Targeted NGS was performed in a BM sample taken in January 2015 at the stage of secAML relapse after the first alloHCT. 67 exons (ex.) (91 amplicons) of nine genes – SF3B1 (ex. 12-16), U2AF1 (ex. 2,6), SRSF2 (ex. 1), ZRSR2 (whole coding region), TET2 (ex. 3-11), ASXL1 (ex. 12), DNMT3A (ex. 2-23), TP53 (ex. 2-11) and RUNX1 (ex. 3-8) – were paired-end sequenced using the targeted NGS (MiSeq system (Illumina)). The sequencing coverage accounted for 31-9095 with mean 250 read length. Bioinformatic analyses comprised mapping reads using the bwa-mem algorithm (BWA, 0.7.10) [5] to the reference genome UCSC Genome Browser (GRCh37/hg19). The detection of variants was performed using GATK (3.5) software [6] with their further annotation completed by snpEff (4.2) [7] based on the dbSNP142 database. For further validation we selected only nonsynonymous substitutions and indels, novel variants or those already annotated with mean AF < 1%. Next, detected variants were confirmed with SSeq and pyrosequencing. Primer sequences, amplicon lengths and localization and PCR conditions are listed in Table 1. Pyrosequencing was carried out according to the manufacturer’s protocol at PyroMark Q48 Autoprep (Qiagen). SSeq was performed according to the protocol described by Kiwerska et al. [8]. The obtained results were analyzed in relation to the patient’s clinical data.

Table 1

Primer sequences, genomic localization of the amplified sequences and reaction conditions

GeneSequence alterationPrimerPrimer sequence 5’-3’Annealing
temperature
Product sizeLocalization (GRCh37/hg19)Sequencing method
TET2c.4044+2dupTFbiotin-CCACTCTTATGGCACCAACATAT55°C87bpchr4:106,182,949-106,183,035Pyrosequencing
RCAATTGCTGCCAATGATTATTTA
STGCTGCCAATGATTATTTAAACSequencing primer
FGGGATTCAAAATGTAAGGG60°C323bpchr4:106,182,817-106,183,139Sanger Sequencing
RTGCAGTGGTTTCAACAATTAAG
c.4076G>TFACTTTCGCATTCACACACACTTT61°C161bpchr4:106,190,734-106,190,894Pyrosequencing
Rbiotin-CCATTCTGCATGTTGTGCAAGTC
STGAACACAGAGCACCASequencing primer
FTCATTCCATTTTGTTTCTGGA55°C400bpchr4:106,190,629-106,191,028Sanger Sequencing
RGCCATGTGGAACTGTGAGTC
ASXL1c.1945G>TFGAGGTCACCACTGCCATAGAGAG55°C135bpchr20:31,022,398-31,022,532Pyrosequencing
Rbiotin-CACAGGCCTCACCACCAT
SGGGGGGGGTGGCCCGSequencing primer
FACCACTGCCATAGAGAGGCG67°C419bpchr20:31,022,404-31,022,822Sanger Sequencing
RGGTTTGGGAGGACAGTAGGG
SF3B1c.1874G>TFbiotin-TTAAGAAGGGCAATAAAGAAGGAA55°C151bpchr2:198,267,412-198,267,562Pyrosequencing
RTTTACATTTTAGGCTGCTGGTCT
SATAGATAACATGGATGAGTASequencing primer
FTTGATTATGGAAAGAAATGGTTG60°C435bpchr2:198,267,189-198,267,623Sanger Sequencing
RAGCCCAAAGGTTTGAGTCC
RUNX1c.509-2A>CFbiotin-TGAAGACAGTGATGGTCAGAGTGA61°C76bpchr21:36,231,843-36,231,918Pyrosequencing
RCCACCAACCTCATTCTGTTTTGTT
SCATTCTGTTTTGTTCTCTATCGTSequencing primer
FGGTAACTTGTGCTGAAGGGC65°C312bpchr21:36,231,673-36,231,984Sanger Sequencing
RGGTTGAACCCAAGGAATCTG

Case report

The patient, a 54-year-old man, was diagnosed with MDS with ring sideroblasts and multilineage dysplasia (MDS-RS-MLD) in April 2008 (3% BM blasts, 74% BM ring sideroblasts, normal karyotype and intermediate disease risk according to the Revised International Prognostic Scoring System (IPSS-R)). In November 2010 due to transformation to secAML the patient underwent conventional chemotherapy with daunorubicin and cytarabine (without hypomethylating agents) followed by alloHCT (May 2011). Second alloHCT (preceded by treatment with chemotherapy with daunorubicin and cytarabine) was performed in February 2015 and third alloHCT (preceded by hypomethylating agents) was performed in March 2018 (Fig. 1). Targeted NGS was performed in a BM sample taken at the moment of the resistance to the implemented therapy for the secAML relapse (January 2015; with PB DC = 42% and BM blasts = 6%) and revealed the presence of the following DNA sequence variants: c.4044+2dupT (splicing variant) and c.4076G>T (p.Arg1359Leu) in TET2, c.1945G>T (p.Gly649*) in ASXL1 and c.1874G>T (p.Arg625Leu) in SF3B1, with AF of: 25.3%, 26.6%, 27.9%, 25.3%, respectively. Presence of all detected variants was confirmed using SSeq and pyrosequencing in both BM and SAL. For each DNA sequence variant the sensitivity of the pyrosequencing method (AF percentage above which the tested sample is considered as positive) was established. To do so, twice the standard deviation of AF for healthy controls + the highest value of AF observed in controls was counted. The cut-off values for c.4044+2dupT TET2, c.4076G>T TET2, c.1945G>T ASXL1 and c.1874G>T SF3B1 DNA variants were determined as: 11.5%, 2.4%, 2%, 4.4%, respectively. AF values obtained by pyrosequencing corresponded to those obtained by NGS. Moreover, for analyzed DNA sequence variants, AF levels obtained by pyrosequencing in BM and SAL were very similar (Fig. 1).

In retrospective analysis of BM (May 2009; November 2010) with pyrosequencing and SSeq, c.4044+2dupT and c.4076G>T TET2, c.1945G>T ASXL1 and c.1874G>T SF3B1 DNA sequence variants were detected. In May 2009 at the MDS stage of disease, AF was 49.2%, 34.5%, 44.1%, 44.7%, respectively (Fig. 1). All detected variants were reported in hematological malignancies according to the Catalogue of Somatic Mutations in Cancer (COSMIC) database [9]; however, except for c.1874G>T SF3B1, the remaining three variants had not been identified in MDS patients so far. In November 2010 transformation to secAML occurred; the patient underwent chemotherapy and was qualified for alloHCT. Using pyrosequencing, AF of above-mentioned variants was determined as 47%, 35.5%, 33.6% and 39.6%, respectively (Fig. 1). AlloHCT from an unrelated donor was performed in May 2011 in complete remission (CR) according ELN 2017 treatment response criteria. Since January 2014 loss of PB DC (70%, 56%, 41%) has been observed and the patient was considered for second alloHCT, but relapse of the disease was diagnosed with 51% of BM blasts. The patient received chemotherapy, but he did not achieve CR according to the ELN 2017 treatment response criteria (no CR). Second alloHCT from the same donor was performed in February 2015 with 8.4% BM blasts and PB DC 27%. Two years later (January 2017) pancytopenia, 13.5% of BM blasts and PB DC of 91% were reported. In BM and SAL obtained from the patient in January 2017 none of the previously detected DNA sequence variants were found with SSeq due to the high DC. On the other hand, all analyzed variants were confirmed in BM with pyrosequencing. In SAL DNA sequence variants were undetectable, probably due to concurrent high DC at disease relapse. Based on leukemic evolution theory and the possibility of additional mutation acquisition we decided to check whether any additional mutation occurred in a BM sample after second relapse of the disease. Using the pyrosequencing method, the presence of RUNX1 (c.509-2A>C), SRSF2 (c.284C>A), TET2 (c.4638G>C) and DNMT3A (c.1014+1G>T; c.2390A>G) DNA sequence variants was demonstrated. Pyrosequencing revealed the presence of the RUNX1 (c.509-2A>C) DNA sequence variant with AF 4.87%. Importantly, the RUNX1 variant was not observed at earlier disease stages. Salvage chemotherapy was implemented, resulting in CR (according ELN 2017 treatment response criteria) and the patient was qualified for maintenance treatment with azacitidine. Subsequent relapse with 16% of BM blasts occurred; thus the patient received chemotherapy, resulting in CR (according ELN 2017 treatment response criteria). The patient was qualified for allogenic haploidentical hematopoietic cell transplantation (allo-haploHCT) from his daughter (March 2018). Since May 2018 acute graft-versus-host disease with pancytopenia has been observed. BM biopsy revealed low cellularity and DC of 100%. In BM and PB obtained from the patient in May 2018 none of the previously detected variants were found with SSeq and pyrosequencing. Four months after allo-haploHCT the patient died from severe graft-versus-host disease with secondary cytopenia, refractory to the treatment.

Discussion

Genetic findings mentioned above indicated that four sequence variants identified with NGS were detected in patient samples since the MDS phase of disease (May 2009) and were persistently reported after transformation to secAML (November 2010, January 2015, January 2017) with the additional RUNX1 c.509-2A>C DNA sequence variant reported since January 2017. Importantly, c.4044+2dupT, c.4076G>T TET2, and c.1945G>T ASXL1 DNA sequence variants were detected for the first time in MDS. Moreover, the presence of c.4044+2dupT and c.4076G>T TET2, c.1945G>T ASXL1 and c.1874G>T SF3B1 DNA sequence variants was confirmed with pyrosequencing and SSeq in both BM and SAL (January 2015). It should be noted that AF detected during pyrosequencing of BM and SAL was very similar (Fig. 1), which highlights the diagnostic utility of SAL testing. However, it is worth noting that the analysis of the mutation presence in SAL is limited by the low amount of leukemic cells as well as presence of high DC in recipients of alloHCT. With clinical suspicion of disease relapse, BM evaluation remains necessary, but our results suggest the clinical usefulness of noninvasive testing with SAL or PB in monitoring the course of disease. Nevertheless, this should be tested on a larger group of patients first.

Comparing pyrosequencing AF results obtained in May 2009, November 2010 and January 2015 and taking into account DC loss (2015 – 42%), we observed a slight AF decrease in tested samples. At the moment of secAML relapse (January 2017) with DC = 91% using the same method, four variants previously identified with targeted NGS, as well as the RUNX1 variant, were detected only in BM but not SAL. Progressive disease in both MDS and secAML is characterized by the emergence of subclones influenced by selection pressure after implemented chemotherapy [10]. In May 2018, after allo-haploHCT (CR, DC = 100%) pyrosequencing results were negative for both BM and PB. Thus, pyrosequencing results corresponded to DC results and additionally introduced information concerning presence of DNA sequence variants in patients after alloHCT. With this approach we are able to confirm the presence of key sequence variants previously detected by NGS. Moreover, the detected DNA sequence variants in selected genes and their prognostic value could influence the therapeutic strategy.

Myelodysplastic syndromes and secAML are very heterogenous diseases and the presence of multiple genetic founding clones translates to different sensitivity to the applied therapy. Moreover, the increasing number of early driver somatic mutations steadily reduces leukemia-free survival [11]. Sequence variants of TET2 and ASXL1 genes are associated with clonal hematopoiesis of indeterminate potential (CHIP). However, CHIP is associated with AF < 10% and presence of single mutations, whereas MDS and AML are associated with higher AF and ≥ 2 mutations [12]. The prognostic significance of TET2 DNA sequence variants in MDS remains unclear [13]. According to Smith et al., TET2 sequence variants in MDS patients do not affect the survival and risk of AML transformation [14]. Furthermore, low expression of the TET2 gene is associated with poor prognosis in MDS patients. Moreover, metaanalysis also revealed that the presence of TET2 DNA sequence variants does not have a significant impact on overall survival. MDS patients with TET2 mutations have increased response rates to hypomethylating agents compared to patients without those mutations. However, treatment with hypomethylating agents does not prolong the overall survival (OS) in MDS patients with TET2 mutations compared to those without. Thus, TET2 status may be used as an additional factor for predicting patients’ response to hypomethylating agents [15]. The SF3B1 variant reported here was found to correlate with ring sideroblasts MDS subtype 2, and this was observed in our case too. Talati et al. associated this variant with favorable prognosis in MDS [16].

Furthermore, patients with SF3B1 mutations showed better treatment responses to erythropoiesis-stimulating agents with no differences in response to hypomethylating agents or lenalidomide [17]. It was also demonstrated that the presence of ASXL1 somatic variants translates to adverse outcome as well as disease relapse and shorter survival after alloHCT in both MDS and MDS/AML patients [18]. It is worth noting that the implemented therapy may also affect time to relapse. The use of anti-cancer treatment results in multiple dynamic changes leading to the subclonal architecture of cancer, driving disease progression and treatment resistance, and is called ‘clonal evolution’ [19]. Mutations in the RUNX1 gene are associated with adverse prognosis in MDS [20], so acquisition of the RUNX1 mutation and its co-occurrence with previously detected DNA variants could reflect disease progression and its influence on the response to the therapy. Importantly, based on Ping Wu’s results, co-occurrence of ASXL1 and RUNX1 DNA sequence variants is related to poor response to hypomethylating agents such as azacytidine or decitabine and their influence on shorter survival [21]. Although the rate of relapse or treatment related mortality is high, alloHCT represents the only curative treatment for high-risk MDS, and it remains crucial to carefully monitor the persistence of genetic alterations after alloHCT [22]. Targeted sequencing of 40 genes in 87 MDS patients undergoing alloHCT identified TET2, DNMT3A and TP53 mutations to be associated with reduced overall survival [23]. Both literature data and our results suggest a significant need to screen MDS patients for the presence of DNA sequence variants in genes of proven clinical relevance. This information could help in treatment management and further course of disease monitoring.

Conclusions

The appropriate MRD monitoring in MDS remains a challenge as the quantity of disease clones is lower than in AML [24]. Our results revealed that the variants previously described for AML (c.4044+2dupT TET2, c.1945G>T ASXL1) were present at the MDS stage of disease. It confirms that the molecular landscape of secAML changes with the disease progression (acquiring c.509-2A>C RUNX1 variant). We suggest: 1) NGS to detect the mutational profile of early MDS and 2) tracking changes in the molecular landscape of MDS/secAML during the course of the disease as essential for improving further treatment. Our paper is the first report of c.4044+2dupT, c.4076G>T TET2 and c.1945G>T ASXL1 DNA sequence variants in MDS.

Notes

[1] Conflicts of interest The authors declare no conflict of interest.

References

1 

Kennedy JA, Ebert BL (2017): Clinical implications of genetic mutations in myelodysplastic syndrome. J Clin Oncol 35: 968-974.

2 

Arber DA, Orazi A, Hasserjian R, et al. (2016): The 2016 revision to the World Health Organization classification of myeloid neoplasms and acute leukemia. Blood 127: 2391-2405.

3 

Gill H, Leung AY, Kwong YL (2016): Molecular and cellular mechanisms of myelodysplastic syndrome: implications on targeted therapy. Int J Mol Sci 17: 440.

4 

Döhner H, Estey E, Grimwade D, et al. (2017): Diagnosis and management of AML in adults: 2017 ELN recommendations from an international expert panel. Blood 129: 424-447.

5 

Li H, Durbin R (2010): Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26: 589-595.

6 

Van der Auwera GA, Carneiro MO, Hartl C, et al. (2013): From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Curr Protoc Bioinformatics 43: 11.10.1-11.10.33.

7 

Cingolani P, Platts A, Wang le L, et al. (2012): A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin) 6: 80-92.

8 

Kiwerska K, Szaumkessel M, Paczkowska J, et al. (2017): Combined deletion and DNA methylation result in silencing of FAM107A gene in laryngeal tumors. Sci Rep 7: 5386.

9 

Tate JG, Bamford S, Jubb HC, et al. (2019): COSMIC: the Catalogue Of Somatic Mutations In Cancer. Nucleic Acids Res 47(D1): D941-D947.

10 

Jacoby MA, Duncavage EJ, Chang GS, et al. (2018): Subclones dominate at MDS progression following allogeneic hematopoietic cell transplant. JCI Insight 3: e98962.

11 

Papaemmanuil E, Gerstung M, Malcovati L, et al. (2013): Clinical and biological implications of driver mutations in myelodysplastic syndromes. Blood 122: 3616-3627; quiz 3699.

12 

Jaiswal S, Fontanillas P, Flannick J, et al. (2014): Age-related clonal hematopoiesis associated with adverse outcomes. N Engl J Med 371: 2488-2498.

13 

Visconte V, Tiu RV, Rogers HJ (2017): Pathogenesis of myelodysplastic syndromes: an overview of molecular and non-molecular aspects of the disease. Blood Res 49: 216-227.

14 

Smith AE, Mohamedali AM, Kulasekararaj A, et al. (2010): Next-generation sequencing of the TET2 gene in 355 MDS and CMML patients reveals low-abundance mutant clones with early origins, but indicates no definite prognostic value. Blood 116: 3923-3932.

15 

Lin Y, Lin Z, Cheng K, et al. (2017): Prognostic role of TET2 deficiency in myelodysplastic syndromes: A meta-analysis. Oncotarget 8: 43295-43305.

16 

Talati C, Sallman D, List A (2017): Lenalidomide: myelodysplastic syndromes with del(5q) and beyond. Semin Hematol 54: 159-166.

17 

Migdady Y, Barnard J, Al Ali N, et al. (2018): Clinical outcomes with ring sideroblasts and SF3B1 mutations in myelodysplastic syndromes: MDS clinical research consortium analysis. Clin Lymphoma Myeloma Leuk 18: 528-532.

18 

Della Porta MG, Galli A, Bacigalupo A, et al. (2016): Clinical effects of driver somatic mutations on the outcomes of patients with myelodysplastic syndromes treated with allogeneic hematopoietic stem-cell transplantation. J Clin Oncol 34: 3627-3637.

19 

Greaves M, Maley CC (2012): Clonal evolution in cancer. Nature 481: 306-313.

20 

Haferlach T, Nagata Y, Grossmann V, et al. (2014): Landscape of genetic lesions in 944 patients with myelodysplastic syndromes. Leukemia 28: 241-247.

21 

Wu P, Weng J, Li M, et al. (2019): Co-occurrence of RUNX1 and ASXL1 mutations underlie poor response and outcome for MDS patients treated with HMAs. Am J Transl Res 11: 3651-3658.

22 

Anthias C, Dignan FL, Morilla R, et al. (2014): Pre-transplant MRD predicts outcome following reduced-intensity and myeloablative allogeneic hemopoietic SCT in AML. Bone Marrow Transplant 49: 679-683.

23 

Bejar R, Stevenson KE, Caughey B, et al. (2014): Somatic mutations predict poor outcome in patients with myelodysplastic syndrome after hematopoietic stem-cell transplantation. J Clin Oncol 32: 2691-2698.

24 

Nassereddine S, Nishihori T, Padron E, et al. (2017): Integrating genomics in myelodysplastic syndrome to predict outcomes after allogeneic hematopoietic cell transplantation. Clin Lymphoma Myeloma Leuk 17: 7-13.

Copyright: © 2021 Polish Society of Experimental and Clinical Immunology This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International (CC BY-NC-SA 4.0) License (http://creativecommons.org/licenses/by-nc-sa/4.0/), allowing third parties to copy and redistribute the material in any medium or format and to remix, transform, and build upon the material, provided the original work is properly cited and states its license.
 
Quick links
© 2024 Termedia Sp. z o.o.
Developed by Bentus.