Blood Journal
Leading the way in experimental and clinical research in hematology

Unraveling the complexity of tyrosine kinase inhibitor–resistant populations by ultra-deep sequencing of the BCR-ABL kinase domain

  1. Simona Soverini1,
  2. Caterina De Benedittis1,
  3. K. Machova Polakova2,
  4. Adela Brouckova2,
  5. David Horner3,
  6. Michele Iacono4,
  7. Fausto Castagnetti1,
  8. Gabriele Gugliotta1,
  9. Francesca Palandri1,
  10. Cristina Papayannidis1,
  11. Ilaria Iacobucci1,
  12. Claudia Venturi1,
  13. Maria Teresa Bochicchio1,
  14. Hana Klamova2,
  15. Federica Cattina5,
  16. Domenico Russo5,
  17. Paola Bresciani6,
  18. Gianni Binotto7,
  19. Barbara Giannini8,
  20. Alexander Kohlmann9,
  21. Torsten Haferlach9,
  22. Andreas Roller9,
  23. Gianantonio Rosti1,
  24. Michele Cavo1,
  25. Michele Baccarani1, and
  26. Giovanni Martinelli1
  1. 1Department of Experimental, Diagnostic and Specialty Medicine, Institute of Hematology “L. e A. Seràgnoli”, University of Bologna, Bologna, Italy;
  2. 2Institute of Hematology and Blood Transfusion, Prague, Czech Republic;
  3. 3Department of Biosciences, University of Milan, Milan, Italy;
  4. 4Roche Applied Science, Monza, Italy;
  5. 5Unit of Blood Diseases and Stem Cell Transplantation University of Brescia, Brescia, Italy;
  6. 6Department of Medical and Surgical Sciences, Azienda Ospedaliero-Universitaria Policlinico, University of Modena e Reggio Emilia, Modena, Italy;
  7. 7Department of Clinical and Experimental Medicine, Hematology Section, Padua University School of Medicine, Padua, Italy;
  8. 8Laboratorio Unico Area Vasta Romagna, Pievesestina di Cesena, Italy; and
  9. 9Munich Leukemia Laboratory, Munich, Germany

Key Points

  • UDS demonstrated that BCR-ABL KD mutations detectable with conventional methods may just be the tip of the iceberg.

  • The information provided by conventional Sanger sequencing may not always be sufficient to predict responsiveness to a given TKI.

Abstract

In chronic myeloid leukemia and Philadelphia chromosome–positive acute lymphoblastic leukemia, tyrosine kinase inhibitor (TKI) therapy may select for drug-resistant BCR-ABL mutants. We used an ultra-deep sequencing (UDS) approach to resolve qualitatively and quantitatively the complexity of mutated populations surviving TKIs and to investigate their clonal structure and evolution over time in relation to therapeutic intervention. To this purpose, we performed a longitudinal analysis of 106 samples from 33 patients who had received sequential treatment with multiple TKIs and had experienced sequential relapses accompanied by selection of 1 or more TKI-resistant mutations. We found that conventional Sanger sequencing had misclassified or underestimated BCR-ABL mutation status in 55% of the samples, where mutations with 1% to 15% abundance were detected. A complex clonal texture was uncovered by clonal analysis of samples harboring multiple mutations and up to 13 different mutated populations were identified. The landscape of these mutated populations was found to be highly dynamic. The high degree of complexity uncovered by UDS indicates that conventional Sanger sequencing might be an inadequate tool to assess BCR-ABL kinase domain mutation status, which currently represents an important component of the therapeutic decision algorithms. Further evaluation of the clinical usefulness of UDS-based approaches is warranted.

Introduction

The first Bcr-Abl tyrosine kinase inhibitor (TKI), imatinib, was introduced in the treatment protocols of chronic myeloid leukemias (CMLs) and Philadelphia chromosome–positive (Ph+) acute lymphoblastic leukemias (ALLs) more than a decade ago.1 Soon after, however, it was observed that BCR-ABL kinase domain (KD) mutated forms with reduced or no sensitivity to imatinib could be selected.2 Second-(dasatinib/nilotinib/bosutinib) and third-generation (ponatinib) TKIs with much fewer insensitive mutations are now already approved or pending approval (Table 13-10). Sequential switch from a TKI to another may rescue response, although further gain of mutations by the same (compound mutations) or different (polyclonal mutations) Ph+ clones is possible.11

View this table:
Table 1

Summary of the BCR-ABL KD amino acid substitutions identified in clinical samples from patients reported to be resistant to the currently approved TKIs

Capillary Sanger sequencing (SS) is the most widely adopted method to assess BCR-ABL KD mutation status3 despite its multiple technical limitations: it cannot robustly identify mutated populations <10% to 15%, it provides only rough estimates of mutated clone abundance, and it cannot discriminate between polyclonal and compound mutations, unless it is preceded by a cumbersome step of cloning. Before the advent of next-generation sequencing (NGS) technologies, however, no method was available that could improve upon these limitations while allowing for scanning of any sequence variant at any position within the KD. NGS carries out thousands to millions of picoliter-scale sequencing reactions simultaneously yielding thousands to millions of sequence reads, each one corresponding to a single, clonally amplified, DNA molecule.12 The use of this approach to sequence a nucleotide position multiple times, thus achieving high sensitivity, is defined as ultra-deep sequencing (UDS). Among NGS technologies commercially available, the Roche 454 Life Sciences is particularly suitable for projects of targeted amplicon resequencing, the average read length being 400 to 600 bp.13 Its application for the UDS of the BCR-ABL KD would allow for: (1) full characterization of the spectrum of minor (<10%-15%) mutated variants; (2) the ability to follow the dynamics of resistant mutations over time; and (3) reconstruction of the clonal architecture of mutated populations in the case of multiple mutations occurring within the same amplicon.

We thus decided to take advantage of a UDS-based approach in order to resolve qualitatively and quantitatively the complexity of mutated Ph+ populations surviving TKIs and to investigate their clonal structure and evolution over time in relation to treatment.

Materials and methods

Patients

We retrospectively selected 33 CML or Ph+ ALL patients who had received sequential treatment with multiple TKIs (2-4 TKIs among imatinib, dasatinib, nilotinib, ponatinib) and had experienced sequential relapses accompanied by selection of 1 or more TKI-resistant mutations. Their main characteristics are presented in Table 2. Written informed consent had been obtained, in accordance with the Declaration of Helsinki. Up to 10 samples were analyzed for each patient, for a total of 106 samples.

View this table:
Table 2

Features of the patients included in the present study

We also randomly selected and analyzed 15 CML patients who had achieved an optimal response to imatinib (according to the European Leukemia Net [ELN] definitions)14 for comparison.

This study was approved by the review boards of the S. Orsola-Malpighi Hospital (Bologna) and of the Institute of Hematology and Blood Transfusion (Prague); it was conducted in accordance with the Declaration of Helsinki.

SS of the BCR-ABL KD

SS of the BCR-ABL KD was performed on an ABI PRISM 3730 (Applied Biosystems) as previously described.15

UDS of the BCR-ABL KD

RNA was converted to complementary DNA (cDNA) with the Transcriptor High-Fidelity cDNA Synthesis kit (Roche Applied Science). To select for the translocated ABL allele, a first step of amplification was performed by polymerase chain reaction (PCR) with a forward primer either on BCR exon 1a (in case of e1a2 BCR-ABL fusion) or on the border of BCR exons 12-13 (in case of b2a2 or b3a2 BCR-ABL fusions) and a reverse primer on ABL, exon 10. A second amplification step was then performed to generate 4 partly overlapping amplicons covering the KD of ABL (supplemental Figure 1, available on the Blood website), tagged with a 10-base “barcode” sequence (multiplex identifier) for sample pooling. Amplifications were done using the FastStart High-Fidelity PCR System kit (Roche Applied Science). UDS was performed on a Roche GS Junior (454-Life Sciences) according to the manufacturer’s instructions. Primer sequences, PCR, and sequencing protocols are detailed in the supplemental Methods. Amplicon Variant Analyzer (version 2.7; 454-Life Sciences) and Sequence Pilot (version 4.0.1; JSI-Medical Systems) were used to align reads to the reference ABL sequence (GenBank accession no. X16416.1) and to calculate variant frequencies. The presence of all relevant mutations was also manually verified by inspection of individual flowgrams at the corresponding positions, with particular attention to homopolymeric regions. Samples harboring multiple mutations were analyzed to gain further insights into the clonal architecture of mutated populations. To this purpose, 2 further amplicons were designed generating clonal reads for cases harboring a P-loop mutation plus a threonine 315 or phenylalanine 317 mutation, and for cases harboring a threonine 315 or phenylalanine 317 mutation plus a methionine 315, phenylalanine 359 or A-loop mutation: the most frequent mutation combinations (supplemental Figure 1). UDS runs were designed to enable high-sensitivity mutation calling, the target sequence coverage ranging from 3438 to 9976 independent reads for each nucleotide position. The sensitivity and reproducibility of an ultra-deep amplicon sequencing approach based on the 454 technology had been previously explored in the framework of the IRON (Interlaboratory Robustness of Next-generation sequencing)16 consortium. More recently, robustness and clinical utility were addressed in an even greater detail by Grossmann et al.17 Moreover, we further explored and confirmed the sensitivity and reproducibility of our UDS-based BCR-ABL KD mutation screening assay in a series of additional experiments detailed in the supplemental Methods.

Because (1) several studies concordantly demonstrated high reliability and reproducibility of 454 technology for variants ≥1%,18-20 and (2) serial dilution experiments (see supplemental Methods for details) confirmed the ability of our assay to detect mutations as low as 1% (whereas lower levels were not investigated), we decided to reduce the likelihood of false-positive results by filtering-out variants with <1% abundance and to exploit high coverage only for haplotype and clonal evolution analyses.

Phylogenetic analysis

Multiple read alignment and clustering was done with Jalview (version 2.8; www.jalview.org). Phylogenetic reconstruction was performed using the bioneighbor joining method with the Kimura 2-parameter substitution model with 1000 bootstrap replicates as implemented in the program SEAVIEW.21,22 To minimize the impact of sequencing errors, only sequences that represent >0.5% of the reads sequenced for at least 1 time point were considered for phylogenetic analyses.

Results

BCR-ABL KD mutation status may be more complex than SS shows

All the nucleotide substitutions that had been detected by SS (“major” mutations) were also detected by UDS, with fairly good concordance between the percentage of variant reads assessed by UDS and the mutation abundance estimated from the relative peak height in the SS chromatogram (Table 3; response definitions in references 14 and 23). In 58 of 106 (55%) samples, however, UDS revealed that mutations undetectable by SS (from now on referred to as “minor” mutations; abundance between 1% and 15%) were also present (Table 3). The type of minor mutations detected by UDS could frequently be accounted for by TKI exposure history because half (57 of 111, 51%) could be recognized to be poorly sensitive either to the TKI being administered or to the previous TKI received, or both (Table 1 reports the mutations know to be resistant to imatinib, dasatinib, nilotinib). In 4 patients (MBC-11, ALL-23, ALL-32, ALL-33), different Ph+ populations were found to have acquired different nucleotide substitutions leading to the same TKI-resistant amino acid change (eg, g>c and g>t at position 903, both resulting in a Q252H mutation; t>c at 1096 and c>a at 1098, both resulting in an F317L mutation). In the remaining cases, minor mutations were either silent (20 of 111, 18%) or never reported in association with TKI resistance (34 of 111, 28%).

View this table:
Table 3

Comparison between mutations detected by SS and mutations detected by UDS and estimated clonal composition of the samples harboring multiple mutations as assessed by UDS

Samples from 15 CML patients who had achieved stable optimal response to imatinib14 were also analyzed by UDS, for comparison. None were found to harbor point mutations at a cutoff level of ≥1%.

A complex clonal texture is uncovered by clonal analysis of samples harboring multiple mutations by UDS

When multiple mutations fell in KD regions ≤450 bp, it was possible to design “ad hoc” amplicons mapping to those regions (supplemental Figure 1) and take advantage of the clonal nature of the sequence reads generated by UDS to see whether multiple mutations belonged to the same (compound mutations) or to different (polyclonal mutations) BCR-ABL transcripts (hence, defined 1 or multiple Ph+ populations). Supplemental Figures 8 and 9 show 2 representative examples. The great majority of the samples turned out to be a complex mosaic of populations harboring the mutations alone as well as in combination (Table 3). Single mutants were 136 of 274 (49.6%); compound mutants were almost as frequent (Figure 1).

Figure 1

Relative frequency of single as opposed to compound mutants. Compound mutants harboring 2 paired mutations were almost as frequent (38.3%) as single mutants (49.6%) and were by far more frequent than triple and quadruple (they accounted for 105 of the 138 [76%] compound mutants overall identified). Mutated Ph+ populations harboring 3 or 4 mutations in the same BCR-ABL molecule were occasionally detected, but in 1 case only they were found to have achieved clonal dominance over those with 1 or 2 mutations.

The landscape of mutated populations is highly dynamic

Longitudinal quantitative follow-up of mutated populations painted an elaborate picture of how the relative frequency of competing populations can ebb and flow over time and with therapeutic intervention (Table 3). Some representative examples are illustrated in Figure 2. The switch from a TKI to another determined the fall of previously dominant population(s) and the rise of new dominant one(s), not necessarily preexisting at the time of switchover (or, at least, not always detectable at the time of switchover with the level of sensitivity allowed by our experimental approach). A new dominant population could be unrelated to the former, rather arising from an unmutated population (as exemplified in Figure 2A-B), or could result from the acquisition of new mutations by the former, which generated compound mutants (Figure 2B-C). In several cases, the evolution of the pattern of mutated populations suggested that the same mutation could have been acquired in parallel by independent populations (ie, 1 unmutated and 1 already harboring a mutation; Figure 2C-D). Selection/deselection of mutated populations could be strikingly rapid. Some compound mutants (M351T+F317L, Y253H+T315I, Y253H+F317L, F359V+T315I) were observed to have higher selective advantage over the respective single mutants. Other compound mutants (H396R+F317L, H396R+T315I, T315I+F317L, E355G+T315I, G250E+F317L, E255V+Y253H, E255V+T315I) were identified that did not overcome the respective single mutants. The E255K+T315I compound mutant was detected in 6 cases: in 3, it became dominant over the E255K and T315I mutants, whereas in another 3 it did not. The triple and quadruple compound mutants detected fluctuated at low levels and were never able to gain dominance, except in patient LBC-17, suggesting that accumulation of >2 mutations, when tolerated, has almost always limited selective advantage.

Figure 2

Mutated populations rise and fall in dominance over time in relation to therapeutic intervention. Graphical illustration of the kinetics of mutated population abundances in 4 representative cases. Arrows indicate the time points at which UDS was performed. Patient IDs are as in Table 3. (A) At the time of first relapse, 7 distinct imatinib-resistant mutated populations were detected. Dasatinib treatment cleared these mutants as quickly as in 2 months, but just as quickly a pan-resistant T315I mutant was found to have emerged. The patient achieved a transient hematologic response after 1 month but lost it shortly after detection of T315I. The question mark indicates that no T315I had been detected by UDS at the time of switchover (at a coverage of 4527 reads, suggesting that either it was present in <1/5000 transcripts, or that it was acquired some time later). (B) At the time of first relapse, a single imatinib-resistant Y253H mutant that accounted for almost 90% of BCR-ABL–positive cells was detected by SS and UDS. After 3 months of dasatinib therapy, the patient had achieved a complete cytogenetic response (no Ph+ metaphases detectable in the bone marrow by standard chromosome banding analysis) although residual disease remained detectable at the molecular level (as assessed by real-time quantitative [RT-Q]-PCR for BCR-ABL transcript). Neither the Y253H (known to be substantially sensitive to dasatinib) nor other mutants were detectable any longer by UDS. After 9 months on dasatinib, the patient was found to have lost the cytogenetic response. UDS showed the coexistence of 3 distinct compound mutants where Y253H was coupled with 3 well-known dasatinib-resistant mutations (a T315I and an F317L resulting from 2 different nucleotide substitutions). It might be hypothesized that the original Y253H-positive cells were never completely eliminated by dasatinib and persisted at very low levels (undetectable by UDS) until they happened to gain a selective advantage again, although de novo acquisition of mutations by previously unmutated cells cannot be ruled out. Interestingly, Y253H and T315I were already detectable by UDS 1 month before. (C) At the time of second relapse, after 6 months of dasatinib therapy, a dasatinib-resistant F317L mutation was detected. UDS portrayed a complex scenario with 3 distinct populations where the former imatinib-resistant Y253H and the newly acquired dasatinib-resistant F317L were present alone and in combination, although the Y253H+F317L compound mutant quantitatively dominated over the F317L- and Y253H-positive ones. It might be hypothesized that the same mutation was acquired in parallel by independent populations (ie, one unmutated and one already positive for the Y253H). (D) After 6 months of second-line dasatinib treatment, during which the patient achieved a complete cytogenetic response but no molecular response, the former imatinib-resistant H396R mutant plus 2 additional populations, 1 harboring a dasatinib-resistant F317L and 1 harboring H396R+F317L, were detected. Three months later, the F317L had become the dominant one, while H396R+F317L and H396R had declined. During dasatinib therapy, a pan-resistant T315I was also acquired because at the time of switchover to nilotinib, a T315I together with H396R+T315I and F317L+T315I compound mutants were already detectable by UDS (4.44%, 0.76%, and 0.22%, respectively). As quickly as in 2 months, nilotinib treatment selected the T315I mutant that expanded, achieving almost full dominance. IM, imatinib; DAS, dasatinib; NIL, nilotinib.

Further insights into BCR-ABL KD sequence evolution: a phylogenetic approach

Figure 3 shows the neighbor joining bootstrap consensus tree for sequences recovered at 5 time points for patient ALL-30. The tree can be considered as consisting of 3 clusters. Reassuringly, within each cluster, both the timing of first observation of each sequence and the types of amino acid substitutions present are fairly consistent, suggesting that, despite the low bootstrap support for key branches, relevant hypotheses can be generated from the phylogenetic analysis. The tree is consistent with E255V occurring once, along the branch dividing cluster A from the rest of the tree. This substitution, which affords resistance to imatinib and nilotinib, is observed in the majority of the sequences in clusters B and C, and was already frequent in ALL-30-01 after 6 weeks of nilotinib treatment (seq128). Indeed, in ALL-30-02 (after 6 weeks on imatinib), seq128 is the dominant haplotype, although other amino acid substitutions are observed in cluster C combined with E255V. The tree and the timing of emergence of these haplotypes are both consistent with the hypothesis that they evolved from clones represented by seq128. In ALL-30-03, 2 other imatinib- and nilotinib-resistant mutations (E255K and Y253H) are observed. However, the tree indicates that these derive from an unmutated haplotype rather than on the predominant E255V background. The final clinically relevant substitution, the pan-resistant T315I, is observed in cluster B, often in conjunction with E255V (and Q252E whose possible clinical relevance is unknown). T315I thus appears to have arisen at the base of this cluster on an E255V background, rather than from the seq128 haplotype observed in cluster C. This observation suggests that, as in the cases of E255K and Y253H, new TKI-resistant mutations can arise from low frequency, possibly slow-proliferating clones that do not represent the majority of the circulating haplotypes at earlier time points. Interestingly, sequences 31, 38, and 104 do not exhibit E255V, suggesting that after withdrawal of imatinib/nilotinib treatment, reversion of this mutation can occur. Indeed, the relatively high frequency of seq104 might be consistent with the hypothesis that such a reversion is positively selected in the context of T315I under dasatinib treatment. An alternative interpretation of the distribution of T315I would be that it arose on an unmutated background, and that E255V subsequently arose independently in some lines within cluster B (as well as on the branch linking cluster C with clusters A and B).

Figure 3

Neighbor joining bootstrap consensus tree (K2P distances) for nucleotide sequences recovered at 5 time points for patient ALL-30 (samples ALL-30-01 through ALL-30-05 in Table 3). The patient was initially enrolled in a clinical study testing the use of nilotinib and imatinib administered in rotation in 6-week cycles. Sample ALL-30-01 (time point 1) was collected after the first (nilotinib) cycle. Sample ALL-30-02 (time point 2) was collected after the second (imatinib) cycle. Sample ALL-30-03 (time point 3) was collected after the fourth cycle. Sample ALL-30-04 was collected after the patient had discontinued the nilotinib-imatinib rotation schedule for hematologic relapse and had received 2 weeks of dasatinib treatment. Sample ALL-30-05 was collected after 6 weeks of dasatinib treatment and a few days before the exitus. Branches with <25% support are collapsed to polytomies. Amino acid substitutions are mapped onto the tree and shown for each sequence along with the percentage of reads that each sequence represented at each time point. “Germline” indicates the unmutated sequence. The scale bar shows substitutions per site. The colors identify the 3 sequence clusters that can be recognized (A-C; see “Further insights into BCR-ABL KD sequence evolution: a phylogenetic approach” for details).

Discussion

UDS has revolutionized the way we can approach the study of drug-resistant cellular populations. Virology is one of the fields that has benefited most from the possibility of highlighting heretofore-undetectable minor mutated variants and performing haplotype analysis, thus allowing for characterization and monitoring of population diversity in HIV and hepatitis viruses.24 There seem to be striking similarities between viral populations and Ph+ leukemia populations, in that they both tend to accumulate mutations to escape antiviral or TKI therapy, respectively. In both scenarios, testing for resistance-associated mutations is important to guide selection of the most appropriate treatment regimen. We have now found that there may be a high degree of heterogeneity, recalling that of viral quasispecies, in BCR-ABL KD sequences from patients failing multiple sequential TKIs. In 55% of the samples, “major” mutations (detectable by SS) were found to be only “the tip of the iceberg”: UDS revealed that additional “minor” (<10%-15%) mutations might be present, even by using a quite “conservative” lower abundance cutoff of 1%. In 51% of the cases, minor mutations could be recognized as poorly sensitive either to the TKI being administered or to the previous TKI received. They most likely corresponded either to outgrowing mutations anticipating an imminent relapse (as for the pan-resistant T315I in samples CP-01-02, AP-10-01, ALL-26-01, ALL-28-01, ALL-29-04) or, more rarely, to “withdrawing” mutants not (yet) entirely deselected by the change in TKI (as in CP-03-06, CP-04-03, CP-06-02, and CP-06-03). In other cases, they could rather be seen as the result of “secondary route(s)” toward resistance followed by some Ph+ cells as an alternative to the one(s) leading to the dominant population(s) (as in AP-10-02, LBC-14-01, ALL-27-03, ALL-32-04, ALL-33-03). In a not negligible proportion of cases, however, minor mutations were either silent or never reported in association with TKI resistance. One would expect such mutations to be always colocalized (“passengers”) on BCR-ABL molecules already harboring a TKI-resistant mutation (“driver”; as in sample CP-03-05); however, this was not always the case, suggesting that, depending on the specific context, a mutant that per se would not be so markedly insensitive to treatment may somehow, at least temporarily, survive and expand to a certain extent, although it will never be able to achieve dominance.

An even higher degree of complexity emerged when we tried to reconstruct the different haplotypes in the samples harboring multiple mutations (Table 3). In the early days of the second-generation TKI era, one study had suggested that compound mutations might be particularly insidious because they may be associated with enhanced oncogenic potential and TKI insensitivity than the separate mutants would exhibit.11 However, the clinical relevance of compound as against polyclonal mutations in patients failing multiple TKIs has long been underestimated. This is mainly because SS, the most widely used method for routine BCR-ABL KD mutation screening,3 precludes the determination of whether multiple mutations are colocated on the same BCR-ABL molecule: only SS of an appropriate number of bacterial colonies with cloned-in BCR-ABL KD would bypass this limitation, but the labor-intensiveness of this approach limits its use to exploratory studies of small patient cohorts.25 Our UDS approach revealed that compound and polyclonal mutations are not 2 mutually exclusive scenarios. This reminds the pattern observed for some patients with myeloproliferative disorders in which JAK2 and TET2 mutations were observed in the same clone and also in different clones.26 Thus, sequential changes in TKI-selective pressure result in heterogeneous mosaics of Ph+ populations harboring different mutations or mutation combinations. Longitudinal observation of the dynamics of these populations in vivo in relation to TKI treatment might suggest that some compound mutants (M351T+F317L, Y253H+T315I, Y253H+F317L) are selectively at an advantage over single mutants, whereas others (H396R+F317L, H396R+T315I, T315I+F317L) are not. On the other hand, the same E255K+T315I (the most frequent compound mutant identified in our samples) was found to achieve dominance in some cases but not in others, suggesting that the patient-specific context, including number and features of coexisting populations, may shape the fitness of a compound mutant and a general rule cannot easily be inferred.

Our results on the complexity of mutated populations and on their dynamics under sequential TKI treatment concur to depict a model in which evolution of BCR-ABL–positive cells is mainly shaped by TKI-selective pressure (whether microenvironment and/or the immune surveillance may also be playing a role is currently unknown) and the fitness of each mutated population is the net result of an “absolute” fitness (the ability to survive treatment depending on the intrinsic sensitivity to the specific TKI administered) and of a “relative” fitness (the ability to survive the competition with all other coexisting populations). For the majority of arising mutants, the fitness will luckily not be high enough to sustain clonal expansion (at least not over a certain threshold), and the population will rather face extinction. Other mutants will expand, although the switch to another TKI may later turn a “fit” into an “unfit” population and lead to an evolutionary “dead-end”. Rare cells within existing mutated populations will gain additional mutations (giving rise to compound mutants), and this will be either an advantage or a disadvantage depending on the type of mutation. Recent whole-genome sequencing studies have challenged the longstanding model of linear evolution of tumors through stepwise accumulation of genetic hits in a founding clone.27 In multiple myeloma, for example, tumor progression has been shown to proceed in a branching rather than in a linear manner, leading to substantial clonal diversity and coexistence of wide genetic heterogeneity.28,29 On a smaller scale, the same seems to apply to drug-resistant BCR-ABL mutants.

Two clinical caveats derive from our findings. We here show that information provided by SS may not always be sufficient to predict responsiveness to a TKI. SS is the currently recommended method for BCR-ABL KD mutation analysis, and its results do influence TKI selection3,30-33: for example, if an imatinib-resistant patient was found to harbor a Y253H by SS, switching to dasatinib rather than to nilotinib should be considered, in light of the fact that this mutant has been shown both in vitro and in vivo to be fully sensitive to dasatinib but poorly sensitive to nilotinib.6,7 If the scenario turned out to be more complex with additional minor mutated populations, clinical efficacy of the selected TKI might be compromised, or only transient.34,35 Such a case should sound as a warning to the physician and a closer monitoring could be beneficial. We acknowledge that in very few instances the information added by UDS would alter TKI selection: in case one of these minor mutated populations were a T315I-positive one, for example, ponatinib might become preferable over dasatinib.4,36 However, if multiple low-level mutations were detected, it would become difficult to predict their role in relapse and progression and to place this information into a therapeutic decision algorithm. We also acknowledge that the clinical significance of low-level mutation detection remains, at present, unclear. Methods with sensitivities ranging between 1% and 10% have been set up and proposed by several authors.37-40 Nevertheless, the optimal lower detection limit allowing for identification of only those mutants that will outgrow and trigger relapse has not been established yet,41,42 although we believe that prospective application of UDS in a large series of patients might ultimately clarify this issue.

Another “warning” from our observations is that sensitivity of a single or a compound mutant to a TKI in vivo might be dictated by more complex factors than the mere in vitro IC50 (half maximal inhibitory concentration) value. If the in vivo sensitivity of any single or compound mutant is influenced by competition with other coexisting mutants, physicians should be warned against using the in vitro IC50 values to estimate the degree of insensitivity of a specific Bcr-Abl mutant protein to a specific TKI.43,44

Many foresee that UDS will soon find routine diagnostic application, especially in the field of hematologic malignancies where the pivotal IRON study has recently provided first evidence of technical feasibility and high concordance of results across multiple laboratories.16 Although it is premature to predict whether UDS will replace SS as the gold standard for BCR-ABL KD mutation analysis, further evaluation of this technology is highly warranted. Currently, the ongoing IRON II study45 is setting-up and testing UDS strategies for the detection of a wide spectrum of mutations, including the BCR-ABL KD screening strategy herein presented.

Authorship

Contribution: S.S. designed the research, performed experiments, analyzed and interpreted results, and wrote the paper; C.D.B. performed experiments and analyzed and interpreted results; K.M.P., A.B., I.I., C.V., M.T.B., and F. Cattina performed experiments; D.H. and M.I. performed phylogenetic analyses; F. Castagnetti, C.P., G.G., F.P., H.K., D.R., P.B., G.B., B.G., and G.R. provided patient samples and clinical data; A.K., T.H., and A.R. provided vital technical and bioinformatics support for the development of the UDS assay; M.C., M.B., and G.M. coordinated the clinical and research team activities and supervised the study; and all authors gave final approval for submission.

Conflict-of-interest disclosure: S.S. was a consultant for and received honoraria from Novartis, Bristol-Myers Squibb, and Ariad. F. Castagnetti, G.G., and G.R. were consultants for and received honoraria from Novartis and Bristol-Myers Squibb. A.K. received honoraria from Roche Diagnostics. T.H. and K.M.P. received research funding from Roche Diagnostics. M.B. and G.M. were consultants for and received honoraria from Novartis, Bristol-Myers Squibb, Ariad, and Pfizer. M.I. is a Roche employee. The remaining authors declare no competing financial interests.

Correspondence: Simona Soverini, Institute of Hematology “L. e A. Seràgnoli”, Via Massarenti 9, 40138 Bologna, Italy; e-mail: simona.soverini{at}unibo.it.

Acknowledgments

The authors thank the Interlaboratory Robustness of Next-generation sequencing (IRON) Phase II Study Group members for helpful discussions on the UDS assay.

This work was supported by European Leukemia Net, Fondazione Cassa di Risparmio in Bologna, Associazione Italiana contro le Leucemie, Programmi di Ricerca di Rilevante Interesse Nazionale 2009 (prot. 2009JSMKY), International Grant Agency Ministry of Health of the Czech Republic NT11555.

Footnotes

  • The online version of this article contains a data supplement.

  • There is an Inside Blood commentary on this article in this issue.

  • The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

  • Submitted March 1, 2013.
  • Accepted June 14, 2013.

References

View Abstract