Massive parallel sequencing, the foundation of next generation sequencing (NGS), allows us to sequence the entire exome (the coding sequences of the genome) of leukemia or lymphoma cells, and can be combined with RNA-Seq to evaluate the transcriptome. Using these techniques, one can search for mutations, indels, gene fusions, copy number alterations, alternative splicing, and gene expression profiles from a blood sample from a person with leukemia. Importantly, these analyses can be performed in just a few days and at modest cost. Given these advances, we might expect data from NGS to be quickly integrated into therapy decisions to effect so-called precision or personalized medicine.21 In fact, cancer treatment was envisioned as one of the most promising applications of PM. In particular, haematological neoplasias, and especially leukemias, were seen as the most direct candidates, given the accessibility of neoplastic cells. Interestingly, and surprisingly, this is not yet widely used, although it continues to generate interest and debates.3 In contrast, our current situation is similar to having many or most pieces of a puzzle, but limited ability to put them in the right place to complete the picture. Several important limitations and challenges associated with using these data to precisely treat persons with leukemia or lymphoma have emerged and are summarized below along with possible solutions. We do not deal here with germline mutations; only with those somatic.
Distinguishing the wheat from the chaff
Many mutations identified by NGS are present before a cell is transformed and are unrelated to leukemia development. These mutations are termed passenger mutations. Passenger mutations are distinguished from driver mutations which cause leukemia transformation. Passenger mutations increase over time and are more frequent in older persons.54 They are also present in persons without leukemia or any blood disorder.6 One way to distinguish passenger from driver mutations is to consider drivers as only those recurrently identified in a substantial proportion of persons with leukemia or clonal hematopoiesis,87 and not age-adjusted normals. However, this approach risks ignoring infrequent mutations (black holes), which may be important in a specific person’s leukemia. In an attempt to identify mutations which may be drivers, we developed OncoScore (revised version submitted for publication in Scientific Reports; available as Bioconductor software at: https://bioconductor.org/packages/release/bioc/html/OncoScore.html and as a web tool at: http://www.galseq.com/oncoscore.html). OncoScore is dynamic software which tracks medical literature in real time, and suggests a numerical score of the probability a mutation is a driver mutation. We tested the possible clinical value of OncoScore in a cohort of 23 persons with Chronic Myeloid Leukemia at diagnosis and found it was better correlated with sustained response to tyrosine kinase-inhibitors than the total number of mutations or Sokal score.9 Further work and analysis will be needed to ascertain the real value of this software.
If the impact of a specific mutation on transformation is unknown, it is also possible to try to predict its impact on the encoded protein function with software such as PolyPhen, DAVID or PROVEAN.1310 These software tools help to indicate whether the observed mutation is likely to cause perturbation in the protein function. Alternatively, animal models could be used to characterize the functional significance of a particular mutation; however, these models require a substantial amount of time, and thus are frequently incompatible with the dynamics and time frames allowed in clinical medicine.
All these tools are imperfect and require improvement and complementation with additional decision making instruments. However, they represent a first step in the direction of differentiating driver and passenger mutations.
Determining mutation hierarchy
Once mutations in a persons’ leukemia cells are identified, we must then add a further dimension: the temporal order in which mutations are acquired. Data from diseases such as MPN-associated myelofibrosis indicated different sequences of mutation acquisition results in different phenotypes despite a similar genotype.14 Reconstruction of the order of acquiring mutations is also important in other settings. Some leukemias first acquire important driver mutations (e.g., BCR/ABL1, NPM/ALK, PML/RARalpha) with substantial tranforming ability such that additional mutations are dispensable. In other diseases, such as myelodysplastic syndrome (MDS), the first mutations are only weakly transforming, and additional genetic alterations are needed for the fully transformed phenotype.
Targeting the earliest driver mutation(s) holds the greatest therapeutic promise when they carry a relevant transforming potential.1615 In this setting (CML, APL, ALK+ lymphomas), PM can change disease prognosis.17 Alternatively, therapies targeting several different driver mutations are a potential therapeutic strategy when the transforming potential of the initial mutation is low. Here, the available evidence for a benefit to patients is more limited, although some promising data are emerging.18 Knowing the mutation hierarchy could also reveal why the same type of mutation, such as ALK containing fusion genes, have different therapeutic implications in different cancers; in lymphoma versus lung cancer, for example, the same drug (crizotinib) obtains quite different therapeutic responses. Hierarchical variant reconstruction is possible but requires sequencing many individual leukemia colonies or sequential studies. In addition, this is presently feasible for myeloid neoplasms, but less so for other cancers. Alternative strategies such as single cell exome/RNA-Seq analysis are being developed to facilitate reconstruction of clonal hierarchy and to eliminate the need to sequence colonies arising from single cells. Statistical methods to infer the order of acquisition of multiple mutations in a cancer were recently suggested by Papaemmanuil et al.19 and Caravagna et al.20 These methods are important as they show that cancer progression follows a defined trajectory, not a random pattern. However they are valid for groups of patients, but cannot assess the order of development of mutations inside a single leukemia.14
All these considerations are valid when the therapeutic strategy is “functional” targeting, e.g., blocking the enzymatic activity of the product of a mutated gene such as a tyrosine kinase. Paradoxically, passenger mutations could represent ideal targets for immune therapy since they are present in all cancer cells, before the driver mutation.21
Signal transduction pathways
An ideal initial driver mutation is one which carries most of the leukemogenic activity and which can be directly targeted. BCR/ABL1 in CML is an example. However, most leukemias and lymphomas are more complex, with a median of >10 mutations/cases and several sub-clones at the time of diagnosis.22 Several software packages, for example DAVID,12 are designed to address this complexity using inputs such as lists of mutated genes or Differentially Expressed Genes (DEG) derived from RNA-Seq analyses. The output can identify the pathway(s) used by the leukemia or lymphoma which could be targeted. Clearly, therapeutic specificity is reduced with this approach as targeting is focused on a pathway used by many normal functions, not solely by the driver mutation. Nonetheless, this approach holds promise when direct targeting is not yet available or feasible.
Which NGS strategy is best?
Whole genome or exome NGS is attractive, but targeted sequencing is gaining favor.19 The strategy is to sequence specific genes (or mutations) identified as recurrent in a specific tumor type and which are actionable, i.e., can be directly targeted with current drugs. Advantages of this approach are: (1) lower cost; (2) higher coverage (mean number of times each nucleotide of the target region is sequenced), which decreases the risk that some targeted loci are insufficiently covered for reliable variant-calling (despite improvements in NGS, it is still common to find parts of targeted genes insufficiently covered); and (3) less complex bioinformatics. Using this approach, one can interrogate <1,000 instead of >13,000 genes. The obvious potential drawback is missing non-conventional mutations (black holes). However, the use of panel-based sequencing is, in some regards, contrary to the original goal of NGS-based PM, namely, characterizing the universe of genetic abnormalities in an individual cancer. Panel-based sequencing is, instead, a standardized approach to PM, which can be unable to reconstruct mutation hierarchy in a person with leukemia or cancer.23 It is difficult to foresee which strategy, comprehensive or panel-based, will be most useful in the future.
In medio stat virtus
Virtue, as is often the case, likely lies between these alternative strategies. Presently, use of pre-defined panels is probably clinically the most feasible one, and will allow physicians to familiarize themselves with the complexities of NGS without being overwhelmed by a mass of omics data. However, in the near future, complete unbiased NGS will, in our opinion, predominate. This evolution may result from improvements in sequencing technology and reduced cost, thus solving, for example, the black holes problem. New, user-friendly bioinformatics tools are likely to be developed, and physicians and researchers will become more familiar with them; new communication skills will also be needed to convey this complex information to patients and their families. Importantly, new clinical trial designs are needed to test the clinical relevance of data from NGS (see below).
A new type of clinical trial is needed
Proving clinical benefit from any therapy intervention requires rigorous methodology, best exemplified by randomized clinical trials. NGS, however, makes it increasingly difficult to identify homogeneous cohorts of persons for study using this trial design. Consequently, new types of clinical trials are needed which preserve the value and rigorous approach of the controlled study, but also take into consideration the NGS-based molecular profile of the tumor.262418 A possible solution could reside in trials in which the strategy of using PM data or not to inform treatment decisions is evaluated in a controlled way, rather than the single therapeutic intervention.
Conclusions
NGS has increased our understanding of leukemia/lymphoma development, and must be translated into better therapy. Consequently, NGS will likely change the way physicians treat lymphomas and leukemias.27 In this scenario, panel-based sequencing will be a bridging-technology to whole-exome or even whole-genome sequencing and RNA-Seq. We are entering an exciting era of PM in general, and leukemia therapy specifically. However, we must also solve the important logistical problems that the use of NGS will inevitably cause. We must find new ways to evaluate therapy strategies, and off-label use of approved drugs must be streamlined and simplified. Finally, we must develop new types of clinical studies and new ways to render the complexity of PM information understandable and useful to patients.
Acknowledgements
This work was supported by Associazione Italiana Ricerca sul Cancro 2013 (IG-14249 to C.G.P.), Associazione Italiana Ricerca sul Cancro 2015 (IG-17727 to R.P.) and by European Union’s Horizon 2020 Marie Skłodowska-Curie Innovative Training Networks (ITN-ETN) under grant agreement No.: 675712CGP, CGP is a member of the European Research Initiative for ALK-Related Malignancies (www.erialcl.net).
The authors regret to be unable to cite many appropriate references, given the limitations allowed by the journal.
References
- Mardis ER, Ding L, Dooling DJ. Recurring mutations found by sequencing an acute myeloid leukemia genome. The New England journal of medicine. 2009; 361(11):1058-1066. PubMedhttps://doi.org/10.1056/NEJMoa0903840Google Scholar
- Collins FS, Varmus H. A new initiative on precision medicine. N Engl J Med. 2015; 372(9):793-795. PubMedhttps://doi.org/10.1056/NEJMp1500523Google Scholar
- Hunter DJ. Uncertainty in the Era of Precision Medicine. N Engl J Med. 2016; 375(8):711-713. PubMedhttps://doi.org/10.1056/NEJMp1608282Google Scholar
- Milholland B, Auton A, Suh Y, Vijg J. Age-related somatic mutations in the cancer genome. Oncotarget. 2015; 6(28):24627-24635. Google Scholar
- Jiang L, Gu ZH, Yan ZX. Exome sequencing identifies somatic mutations of DDX3X in natural killer/T-cell lymphoma. Nat Genet. 2015; 47(9):1061-1066. PubMedhttps://doi.org/10.1038/ng.3358Google Scholar
- Martin GM, Ogburn CE, Colgin LM, Gown AM, Edland SD, Monnat RJ. Somatic mutations are frequent and increase with age in human kidney epithelial cells. Hum Mol Genet. 1996; 5(2):215-221. PubMedhttps://doi.org/10.1093/hmg/5.2.215Google Scholar
- Genovese G, Kahler AK, Handsaker RE. Clonal hematopoiesis and blood-cancer risk inferred from blood DNA sequence. N Engl J Med. 2014; 371(26):2477-2487. PubMedhttps://doi.org/10.1056/NEJMoa1409405Google Scholar
- Jaiswal S, Fontanillas P, Flannick J. Age-related clonal hematopoiesis associated with adverse outcomes. N Engl J Med. 2014; 371(26):2488-2498. PubMedhttps://doi.org/10.1056/NEJMoa1408617Google Scholar
- Piazza R, Mologni L, Ramazzotti D.Paper presented at: ; Google Scholar
- Adzhubei IA, Schmidt S, Peshkin L. A method and server for predicting damaging missense mutations. Nat Methods. 2010; 7(4):248-249. PubMedhttps://doi.org/10.1038/nmeth0410-248Google Scholar
- Choi Y, Sims GE, Murphy S, Miller JR, Chan AP. Predicting the functional effect of amino acid substitutions and indels. PloS One. 2012; 7(10):e46688. PubMedhttps://doi.org/10.1371/journal.pone.0046688Google Scholar
- Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009; 4(1):44-57. PubMedhttps://doi.org/10.1038/nprot.2008.211Google Scholar
- Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009; 37(1):1-13. PubMedhttps://doi.org/10.1093/nar/gkn923Google Scholar
- Ortmann CA, Kent DG, Nangalia J. Effect of mutation order on myeloproliferative neoplasms. N Engl J Med. 2015; 372(7):601-612. PubMedhttps://doi.org/10.1056/NEJMoa1412098Google Scholar
- Kantarjian H, Sawyers C, Hochhaus A. Hematologic and cytogenetic responses to imatinib mesylate in chronic myelogenous leukemia. N Engl J Med. 2002; 346(9):645-652. PubMedhttps://doi.org/10.1056/NEJMoa011573Google Scholar
- Gambacorti-Passerini C, Messa C, Pogliani EM. Crizotinib in anaplastic large-cell lymphoma. N Engl J Med. 2011; 364(8):775-776. PubMedhttps://doi.org/10.1056/NEJMc1013224Google Scholar
- Gambacorti-Passerini C, Piazza R. Imatinib–A New Tyrosine Kinase Inhibitor for First-Line Treatment of Chronic Myeloid Leukemia in 2015. JAMA Oncol. 2015; 1(2):143-144. Google Scholar
- Mody RJ, Wu YM, Lonigro RJ. Integrative Clinical Sequencing in the Management of Refractory or Relapsed Cancer in Youth. JAMA. 2015; 314(9):913-925. PubMedhttps://doi.org/10.1001/jama.2015.10080Google Scholar
- Papaemmanuil E, Gerstung M, Bullinger L. Genomic Classification and Prognosis in Acute Myeloid Leukemia. N Engl J Med. 2016; 374(23):2209-2221. PubMedhttps://doi.org/10.1056/NEJMoa1516192Google Scholar
- Caravagna G, Graudenzi A, Ramazzotti D. Algorithmic methods to infer the evolutionary trajectories in cancer progression. Proc Natl Acad Sci U S A. 2016; 113(28):E4025-4034. PubMedhttps://doi.org/10.1073/pnas.1520213113Google Scholar
- Snyder A, Makarov V, Merghoub T. Genetic basis for clinical response to CTLA-4 blockade in melanoma. N Engl J Med. 2014; 371(23):2189-2199. PubMedhttps://doi.org/10.1056/NEJMoa1406498Google Scholar
- Walter MJ, Shen D, Ding L. Clonal architecture of secondary acute myeloid leukemia. N Engl J Med. 2012; 366(12):1090-1098. PubMedhttps://doi.org/10.1056/NEJMoa1106968Google Scholar
- Xu X, Hou Y, Yin X. Single-cell exome sequencing reveals single-nucleotide mutation characteristics of a kidney tumor. Cell. 2012; 148(5):886-895. PubMedhttps://doi.org/10.1016/j.cell.2012.02.025Google Scholar
- Willyard C. ‘Basket studies’ will hold intricate data for cancer drug approvals. Nat Med. 2013; 19(6):655. PubMedhttps://doi.org/10.1038/nm0613-655Google Scholar
- Patterson SE, Liu R, Statz CM, Durkin D, Lakshminarayana A, Mockus SM. The clinical trial landscape in oncology and connectivity of somatic mutational profiles to targeted therapies. Hum Genomics. 2016; 10:4. Google Scholar
- Schwaederle M, Zhao M, Lee J. Impact of precision medicine in refractory malignancies: A meta-analysis of 13,203 patients in phase I clinical trials. J Clin Oncol. 2016; 34(suppl)Google Scholar
- Jameson JL, Longo DL. Precision medicine–personalized, problematic, and promising. N Engl J Med. 2015; 372(23):2229-2234. PubMedhttps://doi.org/10.1056/NEJMsb1503104Google Scholar