Characterization of 46 patient-specific BCR-ABL1 fusions and detection of SNPs upstream and downstream the breakpoints in chronic myeloid leukemia using next generation sequencing
© Linhartova et al.; licensee BioMed Central. 2015
Received: 22 October 2014
Accepted: 8 April 2015
Published: 18 April 2015
In chronic myeloid leukemia, the identification of individual BCR-ABL1 fusions is required for the development of personalized medicine approach for minimal residual disease monitoring at the DNA level. Next generation sequencing (NGS) of amplicons larger than 1000 bp simplified and accelerated a process of characterization of patient-specific BCR-ABL1 genomic fusions. NGS of large regions upstream and downstream the individual breakpoints in BCR and ABL1 genes, respectively, also provided information about the sequence variants such are single nucleotide polymorphisms.
BCR-ABL1 transcript level monitoring is a crucial part of therapy response evaluation in patients with chronic myeloid leukemia (CML). Molecular response (MR) in CML is expressed in terms of log reduction in BCR-ABL1 transcript levels from a standardized baseline on an international scale . An increasing number of chronic phase CML patients treated with tyrosine kinase inhibitors (TKIs) achieve sustained deep MR, which has led to the initialization of TKI discontinuation studies [2,3]. Although existing results showing that permanent TKI discontinuation is feasible in a proportion of patients, more than half experience molecular relapse. More sensitive MRD detection might allow to more accurately stratifying patients according to the depth of response, hence to the greater or smaller likelihood of relapse after discontinuation. It was recently shown that MRD could be detected at DNA level despite patients were in stable MR with undetectable BCR-ABL1 transcripts [4,5]. It is supposed that quantification of BCR-ABL1 gene would be useful for MRD monitoring in patients with therapy cessation.
At the genomic level the BCR-ABL1 fusion is unique to each CML patient, since breakpoints are scattered within intron 13 (718 bp) or 14 (2128 bp) of BCR and within intron 1 of ABL1 (140 Kbp). The exact genomic sequence of the BCR-ABL1 fusion must therefore be characterized to design a personalized real-time polymerase chain reaction (PCR) assay. The sensitivity of personalized real-time PCR assay quantifying clone-specific fusions was reported to be down to less than 10−5 (depending on the amount of DNA entering the reaction) when nested PCR was used [6,7]. As amplification of large regions (˃1 Kbp) is needed, the characterization of BCR-ABL1 genomic fusions using conventional sequencing has been a challenging task [8-10]. With the aim to characterize patient-specific BCR-ABL1 fusions, we set up a strategy based on the generation of large amplicons by multiplex Long-Range PCR (mLR-PCR) and on subsequent ultra-wide sequencing using NGS technology. NGS provided sequences of large regions upstream and downstream the BCR-ABL1 fusions, allowing us to assess the frequency of annotated single nucleotide polymorphisms (SNPs), which could be associated with the fusion gene formation.
Next generation sequencing of large amplicons carrying BCR-ABL1 fusions
The patient cohort consisted of 24 males and 22 females with CML. Two BCR-ABL1 positive cell lines K562 and JURL-MK were included as controls of BCR-ABL1 fusion identification. Twelve healthy individuals were used as controls for SNP detection in the BCR region from exon 13 to exon 15.
Patient-specific BCR-ABL1 fusion characterization
SNPs detection upstream and downstream from the BCR-ABL1 fusions
Downstream the individual breakpoints in ABL1, 56 annotated SNPs were detected in 15/46 patients (median 3 SNPs/patient, range 1–11). Only 10 of them recurred in two patients, whereas the remaining 46 were detected only once. This heterogeneity in SNP detection was caused by the extreme length of intron 1 (Figure 1, Figure 2A), where the breakpoints were spread, so that the individual ABL1 regions sequenced could not overlap.
NG_009244.1:g.114249A > G|rs527236141|ss1227131697;NG_009244.1:g.114278G > T|rs527236142|ss1227132188;NG_009244.1:g.115017 T > A|rs527236143|ss1227132282.
The 300 bp region downstream of primers hybridizing to BCR exon 13 was sequenced in all cases (Figure 2B). Four SNPs were identified within this region altogether in 9% of patients with e13a2 mRNA fusion and in 12% of patients with e14a2 mRNA. Other 16 SNPs were detected almost exclusively in patients with e14a2 mRNA fusion, which is in accordance with the fact, that the BCR region sequenced in these patients was larger than in the patients with e13a2 mRNA. The minor allele frequency among the BCR-ABL1 alleles analyzed corresponds to that found in the 1000 GenomePhase Population and in the cohort, although small, of healthy donors (Figure 2B).
Although increasing number of studies has been published in last few years [6-16], analysis of BCR-ABL1 at the DNA level has been limited and difficult so far. Our mLR-PCR-NGS approach streamlined the laboratory workflow for patient-specific BCR-ABL1 fusion characterization allowing personalizing real-time PCR assays for patients with MRD. The large collection of 46 BCR-ABL1 fusion sequences was obtained in the short time and is at disposal in the NCBI Nucleotide database. Moreover, the NGS enabled to detect SNPs upstream and downstream the patient-specific breakpoints of the BCR-ABL1 gene. NGS of this region in large series of patients has the potential to provide insights on the association between particular SNPs or haplotypes and CML. Although we describe the method on BCR-ABL1 in patients with CML, the strategy can be used in the same manner as a base for DNA-based MRD monitoring in BCR-ABL1 positive acute lymphoblastic leukemia (ALL), where significant differences between MRD levels measured by BCR-ABL1 expression and Ig/TCR DNA-based approach were described . Moreover identical approach can be used for breakpoints localization and sequence features calling in other large fusion genes associated with various cancers, e.g. TEL/ABL1 (atypical CML or ALL) , TEL/AML1 (acute lymphoid leukemia) , TMPRSS2/ERG (prostate cancer)  or EML4/ALK (non-small cell lung cancer) .
This work has been supported by the Internal Grant Agency of Ministry of Health of the Czech Republic grant number NT11555, the conceptual development of research organization (00023736) from Ministry of Health of the Czech Republic and ERDF OPPK CZ.2.16/3.1.00/28007. The authors thank also the Czech Leukemia Study Group for Life - CELL for support. LH was supported by GAUK 554214.
- Cross NC, White HE, Müller MC, Saglio G, Hochhaus A. Standardized definitions of molecular response in chronic myeloid leukemia. Leukemia. 2012;26:2172–5.View ArticlePubMedGoogle Scholar
- Falchi L, Kantarjian HM, Wang X, Verma D, Quintás-Cardama A, O'Brien S, et al. Significance of deeper molecular responses in patients with chronic myeloid leukemia in early chronic phase treated with tyrosine kinase inhibitors. Am J Hematol. 2013;88:1024–9.View ArticlePubMedGoogle Scholar
- Hehlmann R, Müller MC, Lauseker M, Hanfstein B, Fabarius A, Schreiber A, et al. Deep molecular response is reached by the majority of patients treated with imatinib, predicts survival, and is achieved more quickly by optimized high-dose imatinib: results from the randomized CML-study IV. J Clin Oncol. 2014;32:415–23.View ArticlePubMedGoogle Scholar
- Ross DM, Branford S, Seymour JF, Schwarer AP, Arthur C, Yeung DT, et al. Safety and efficacy of imatinib cessation for CML patients with stable undetectable minimal residual disease: results from the TWISTER study. Blood. 2013;122:515–22.View ArticlePubMedGoogle Scholar
- Rousselot P, Charbonnier A, Cony-Makhoul P, Agape P, Nicolini FE, Varet B, et al. Loss of major molecular response as a trigger for restarting tyrosine kinase inhibitor therapy in patients with chronic-phase chronic myelogenous leukemia who have stopped imatinib after durable undetectable disease. J Clin Oncol. 2014;32:424–30.View ArticlePubMedGoogle Scholar
- Bartley PA, Ross DM, Latham S, Martin-Harris MH, Budgen B, Wilczek V, et al. Sensitive detection and quantification of minimal residual disease in chronic myeloid leukaemia using nested quantitative PCR for BCR-ABL DNA. Int J Lab Hematol. 2010;32:e222–228.View ArticlePubMedGoogle Scholar
- Morley AA, Latham S, Brisco MJ, Sykes PJ, Sutton R, Hughes E, et al. Sensitive and specific measurement of minimal residual disease in acute lymphoblastic leukemia. J Mol Diagn. 2009;11:201–10.View ArticlePubMed CentralPubMedGoogle Scholar
- Sobrinho-Simões M, Wilczek V, Score J, Cross NC, Apperley JF, Melo JV. In search of the original leukemic clone in chronic myeloid leukemia patients in complete molecular remission after stem cell transplantation or imatinib. Blood. 2010;116:1329–35.View ArticlePubMedGoogle Scholar
- Waller CF, Dennebaum G, Feldmann C, Lange W. Long-template DNA polymerase chain reaction for the detection of the bcr/abl translocation in patients with chronic myelogenous leukemia. Clin Cancer Res. 1999;5:4146–51.PubMedGoogle Scholar
- Krumbholz M, Karl M, Tauer JT, Thiede C, Rascher W, Suttorp M, et al. Genomic BCR-ABL1 breakpoints in pediatric chronic myeloid leukemia. Genes Chromosomes Cancer. 2012;51:1045–53.View ArticlePubMedGoogle Scholar
- Ross DM, Branford S, Seymour JF, Schwarer AP, Arthur C, Bartley PA, et al. Patients with chronic myeloid leukemia who maintain a complete molecular response after stopping imatinib treatment have evidence of persistent leukemia by DNA PCR. Leukemia. 2010;24:1719–24.View ArticlePubMedGoogle Scholar
- Score J, Calasanz MJ, Ottman O, Pane F, Yeh RF, Sobrinho-Simões MA, et al. Analysis of genomic breakpoints in p190 and p210 BCR-ABL indicate distinct mechanisms of formation. Leukemia. 2010;24:1742–50.View ArticlePubMedGoogle Scholar
- Burmeister T, Gröger D, Kühn A, Hoelzer D, Thiel E, Reinhardt R. Fine structure of translocation breakpoints within the major breakpoint region in BCR-ABL1-positive leukemias. DNA Repair (Amst). 2011;10:1131–7.View ArticleGoogle Scholar
- Mattarucchi E, Spinelli O, Rambaldi A, Pasquali F, Lo Curto F, Campiotti L, et al. Molecular monitoring of residual disease in chronic myeloid leukemia by genomic DNA compared with conventional mRNA analysis. J Mol Diagn. 2009;11:482–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Karl M, Krumbholz M, Tauer JT, Jacobs U, Metzler M, Suttorp M. Identification of the genomic BCR-ABL1 fusion sequence from blood specimen stored on filter paper. Leuk Res. 2013;37:117–9.View ArticlePubMedGoogle Scholar
- Shibata Y, Malhotra A, Dutta A. Detection of DNA fusion junctions for BCR-ABL translocations by Anchored ChromPET. Genome Med. 2010; doi:10.1186/gm191.Google Scholar
- Zaliova M, Fronkova E, Krejcikova K, Muzikova K, Mejstrikova E, Stary J, et al. Quantification of fusion transcript reveals a subgroup with distinct biological properties and predicts relapse in BCR/ABL-positive ALL: implications for residual disease monitoring. Leukemia. 2009;23:944–51.View ArticlePubMedGoogle Scholar
- Zuna J, Zaliova M, Muzikova K, Meyer C, Lizcova L, Zemanova Z, et al. Acute leukemias with ETV6/ABL1 (TEL/ABL) fusion: poor prognosis and prenatal origin. Genes Chromosomes Cancer. 2010;49:873–84.View ArticlePubMedGoogle Scholar
- Zuna J, Hrusak O, Kalinova M, Muzikova K, Stary J, Trka J. TEL/AML1 positivity in childhood ALL: average or better prognosis? Czech Paediatric Haematology Working Group. Leukemia. 1999;13:22–4.View ArticlePubMedGoogle Scholar
- Wang J, Cai Y, Ren C, Ittmann M. Expression of variant TMPRSS2/ERG fusion messenger RNAs is associated with aggressive prostate cancer. Cancer Res. 2006;66:8347–51.View ArticlePubMedGoogle Scholar
- Choi YL, Takeuchi K, Soda M, Inamura K, Togashi Y, Hatano S, et al. Identification of novel isoforms of the EML4-ALK transforming gene in non-small cell lung cancer. Cancer Res. 2008;68:4971–6.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.