- Open Access
Translating insights into tumor evolution to clinical practice: promises and challenges
Genome Medicinevolume 11, Article number: 20 (2019)
Accelerating technological advances have allowed the widespread genomic profiling of tumors. As yet, however, the vast catalogues of mutations that have been identified have made only a modest impact on clinical medicine. Massively parallel sequencing has informed our understanding of the genetic evolution and heterogeneity of cancers, allowing us to place these mutational catalogues into a meaningful context. Here, we review the methods used to measure tumor evolution and heterogeneity, and the potential and challenges for translating the insights gained to achieve clinical impact for cancer therapy, monitoring, early detection, risk stratification, and prevention. We discuss how tumor evolution can guide cancer therapy by targeting clonal and subclonal mutations both individually and in combination. Circulating tumor DNA and circulating tumor cells can be leveraged for monitoring the efficacy of therapy and for tracking the emergence of resistant subclones. The evolutionary history of tumors can be deduced for late-stage cancers, either directly by sampling precursor lesions or by leveraging computational approaches to infer the timing of driver events. This approach can identify recurrent early driver mutations that represent promising avenues for future early detection strategies. Emerging evidence suggests that mutational processes and complex clonal dynamics are active even in normal development and aging. This will make discriminating developing malignant neoplasms from normal aging cell lineages a challenge. Furthermore, insight into signatures of mutational processes that are active early in tumor evolution may allow the development of cancer-prevention approaches. Research and clinical studies that incorporate an appreciation of the complex evolutionary patterns in tumors will not only produce more meaningful genomic data, but also better exploit the vulnerabilities of cancer, resulting in improved treatment outcomes.
Over time, the therapeutic approach to cancer is evolving from targeting the clinical phenotype (tumor size, location, stage, histological type, and grade), to targeting a molecular phenotype (such as surface receptor status or the presence of activating or sensitizing mutations) [1, 2]. The clinical phenotype can be targeted spatially with surgery and radiotherapy or systemically using cytotoxic chemotherapies. The molecular phenotype has been targeted by both direct and indirect endocrine manipulation, by an array of small molecule inhibitors, and by monoclonal antibody therapies. Both approaches typically consider the target to be static (to be treated until clinical failure) and homogeneous (one sample represents all tumor cells).
The application of evolutionary concepts to cancer was proposed several decades ago by Peter Nowell . Reliable exploration of the degree of variation within and between cancers has only become possible with the increasing availability of next generation sequencing and associated computational analysis [4,5,6].
All of the cells within a tumor are unique, comprising different somatic variants and epigenetic and transcriptomic states. Even normal cells are likely to accrue approximately three somatic mutations every cell cycle [7, 8]. Most of these changes will have no functional impact and are ‘passengers’ on the cells’ evolutionary journey (Box 1). Somatic mutations (or epigenetic changes) that have an advantageous functional impact are ‘drivers’ and will allow a cell to expand clonally and outcompete its neighbors. When a clonal expansion goes to completion, the entire population will be ‘clonally’ descended from that founder cell, or clone. The last complete clonal expansion will have arisen from the most recent common ancestor (MRCA), defined as the most recent individual cell from which all existing cancer cells in a cancer sample are descendants. If a clonal expansion or sweep is incomplete, the expanded population is subclonal, comprising only a fraction of the tumor cells. Diverging subclones with mutually exclusive mutations can co-exist within a tumor . Intra-tumor heterogeneity, or the presence of subclones possessing private mutations within a tumor, has been observed across many cancer types and seems to be nearly ubiquitous [10, 11].
The dynamics of evolution in cancer are still not fully understood . Traditionally, mutation and selection are thought to be slow iterative processes that occur throughout a cancer’s lifetime, a process of gradual evolution. The patterns of mutations observed in some tumors, however, suggest that mutations can also be acquired in sudden bursts, leading to punctuated evolutionary steps [13,14,15,16,17,18,19].
An emerging wealth of cancer genome sequencing data is informing our understanding of tumor evolution, and will cause a fundamental paradigm shift in our approach to cancer. This will impact all aspects of cancer management, including cancer therapy, monitoring, early detection, and prevention (Table 1).
Measuring intra-tumor heterogeneity and tumor evolution
Implicit in the heterogeneity of tumor cells and essential for evolution is variation in either the genome or the epigenome [20,21,22]. Although epigenetic heterogeneity has been shown to have prognostic utility [23,24,25,26] and is the subject of intense study, genetic heterogeneity is better understood at present, and is the focus of this review.
Intra-tumor heterogeneity and evolution can be inferred from the pattern of mutations that is detected. Clonal mutations, which are common to all cells within a tumor, were present in the tumor cells’ most recent common ancestor, whereas subclonal mutations were acquired later and are therefore only found in a proportion of tumor cells (Box 1). The frequency of a mutation in sequencing data (the variant allele frequency (VAF)) can be used to establish its clonality. VAF is influenced by both the proportion of cells that possess the mutation and the number of both mutated and un-mutated copies of that DNA locus. Mutation frequencies can be estimated by sampling, which has intrinsic spatial, genomic, and statistical limitations (Fig. 1). Intra-tumor heterogeneity has been extensively explored using exome or genome sequencing of multiple regions of resected primary tumors [9, 12, 19, 27,28,29]. Paired primary–metastasis studies and post-mortem studies have allowed detailed insight into the evolution and patterns of spread of metastases [30,31,32,33]. Intra-tumor heterogeneity has been shown to be prognostic across cancer types [10, 34, 35], and is predominantly associated with the degree and heterogeneity of aneuploidy. It has also been shown to impact therapy: potentially targetable driver mutations can be subclonal, suggesting that treatment would only be partially effective .
Describing tumor evolution requires measurement over time. Models of tumorigenesis, such as the ‘Vogelgram’, were created by sampling different stages of cancer progression across a population [37, 38]. The initial Vogelgram in colorectal cancer was established by probing a limited number of putative oncogenes, identified from hereditary cases, across the histologically defined spectrum of disease . Mutations that are found across different stages of disease are assumed to arise early in tumor evolution, whereas those found only in established invasive cancers can be assumed to occur later in tumor evolution. Rarely, in individuals with predisposing risk factors such as those who have colitis-associated colon cancer, the whole spectrum of tumor progression can be observed simultaneously . Despite applying modern genomic techniques, models of progression can remain elusive if the genome is already markedly aberrated in pre-invasive lesions, as in the precursors of lung squamous cell carcinoma . This modeling approach also relies on the assumption that cancers of the same histology have a highly stereotyped genetic progression that is common to different tumors.
Computational approaches have been developed to infer the history of an individual tumor that is already established from its own genome, as recently reviewed [41, 42]. Although these approaches typically allow only partial reconstruction of a tumor’s evolutionary history, from a single biopsy, aggregating results across multiple tumors can be a powerful approach . Taking multiple samples from the same tumor over time or across space can also significantly increase the power of these reconstruction approaches [41, 42]. In metastatic solid organ tumors, repeated sampling over time is challenging, so hematological malignancies have been studied most extensively in this context [43,44,45,46]. Circulating tumor DNA (ctDNA) and cells shed from solid tumors offer the potential to track subclonal mutations, albeit with limited sensitivity and specificity.
Most DNA sequencing has been performed on pooled DNA from multiple cells and, consequently, ambiguity can remain as to whether mutations co-occur in the same cell. Single-cell sequencing can overcome this, albeit at higher cost and at the expense of substantial sequencing artifacts [47,48,49,50]. High-throughput techniques have been developed for analyzing large numbers of single cells, although these methods are most advanced for transcriptome sequencing . Single-cell sequencing of other ‘omic layers is currently relatively costly and available for fewer cells , but exciting high-throughput approaches are now emerging . Techniques to analyze multiple layers simultaneously have also been developed recently [54,55,56], but these are currently costly and lower throughput. These ‘multi-omic’ approaches are likely to significantly improve the interpretation of non-genetic cellular heterogeneity. Such interpretation is also confounded by heterogeneity among non-tumor cells that results from the variety of cell types and states within a tumor [57, 58].
Future approaches for measuring tumor heterogeneity that could be used clinically would need to satisfy the following criteria: (i) sampling should be minimally invasive or performed as part of tumor resection; (ii) sampling of the tumor should be as comprehensive as possible, ideally without any spatial biases; (iii) sample handling and preservation will need to be simple and readily available in the clinic; (iv) simple proxy biomarkers need to be available to assay heterogeneity reliably; and (v) assays need to be rapid and cost-effective.
Recently, a conceptual consideration of how evolution and heterogeneity could be summarized was explored in a consensus statement by Maley et al. . They proposed binary divisions of the degree of heterogeneity (diversity, D) and evolution (rate of change, ∆) that could be combined in a single four-level Evo-Index. As yet, it is not clear how these scores would be generated or whether such a simple binary system is informative.
Can tumor evolution guide cancer therapy?
The rational design of cancer therapies based on genomic data has to date, with a few notable exceptions, been expensive and has delivered limited benefit to patients . Even therapies specifically targeting prevalent tumor mutations, such as the BRAF V600E mutation in melanoma  and a variety of EGFR point mutations in lung cancer , only lead to relatively short-lived tumor responses. Understanding the heterogeneity that exists within tumors and their ability to evolve in response to therapy may allow more optimized treatment strategies (Table 1).
Individual clonal therapies
The simplest conceivable therapeutic approach is to target individual clonal mutations. By targeting mutations that are present in all tumor cells, the entire tumor could in theory be eradicated. Previous targeted therapies have, to some degree, implicitly relied on the presumption that mutations that are highly prevalent in different tumors are probably early events in tumorigenesis and therefore likely to be clonal.
In most cases, single clonal mutations, which are thought to be functionally relevant driver mutations, have been targeted directly. In established cancers, this invariably results in the acquisition of treatment resistance. The simplest examples are the resistance to endocrine therapy in metastatic breast and prostate cancer. The mechanisms of these resistance phenomena are now relatively well understood. Many breast cancers depend on estrogen signaling and are initially sensitive to therapies that reduce the level of circulating estrogen or that target the cellular estrogen receptor, such as aromatase inhibitors or selective estrogen receptor modulators, respectively. Treatment resistance frequently arises when tumor cells develop constitutive activity in the estrogen receptor through mutation of its gene, ESR1 . Likewise, prostate cancers are almost ubiquitously driven by androgen signaling, sensitizing them to chemical or surgical castration. Prostate cancer cells compensate for medically depleted circulating androgen levels through a number of different mechanisms, including amplification of the androgen receptor . Gundem et al.  demonstrated that multiple separate tumor cell populations, across distinct metastatic sites, can develop unique androgen receptor amplifications—a demonstration of convergent evolution. The widespread evolution of resistance suggests that clonal monotherapies are unlikely to achieve permanent tumor control or cure. For those with slow-paced advanced disease, or those who would not tolerate more intensive therapy, individual therapies will continue to play an important role. Most responses to targeted therapies, however, are both incomplete and short-lived and require improvement (Fig. 2a).
Even when a mutation is not treated directly, tumors can develop resistance. Synthetic lethality is a treatment approach that exploits a cellular vulnerability exposed by a clonal driver mutation. BRCA mutations in breast and ovarian cancer, both inherited or acquired, increase genomic instability due to disruption of the repair of double-strand DNA breaks, which not only produces variation during tumorigenesis but also increases the reliance of these tumors on other DNA-repair mechanisms. This is exploited for therapy by inhibiting the single-stranded DNA repair PARP enzymes [65, 66]. PARP inhibition causes the accumulation of lethal DNA damage specifically in tumor cells. BRCA mutations can, however, undergo somatic reversal in multiple tumor subclones, leading to resistance to PARP inhibition [67,68,69].
Resistance to therapy typically results from mutations, which may pre-exist or can appear subsequent to the therapy, or from non-genetic factors. Mutations that exist prior to treatment exposure might be rare, and therefore undetectable by present assays. Once treatment creates selective pressure, resistant cells carrying these mutations will persist and become apparent. It is possible, and perhaps likely in larger tumors, that most resistance mutations exist prior to therapy exposure, even for conventional cytotoxic therapies . Resistance mutations may also occur de novo after treatment exposure, perhaps having been induced by iatrogenic mutagenesis [70, 71]. Unless these mutations are of a distinct type, known to be induced by therapy, it is difficult to exclude the possibility that they did not exist prior to treatment at a very low and undetectable frequency. Resistance may also be non-genetic and either related to cell state plasticity or to a specific molecular resistance pathway, such as aurora kinase activation in anti-EGFR-treated lung cancer [72, 73], but further understanding of these non-heritable resistance mechanisms is needed. There are broad principles of treatment resistance that are common between cancer and infectious diseases : like tumor cell populations, pathogens can be also be genetically heterogeneous , and as observed in the examples of HIV and Mycobacterium tuberculosis, they rarely have prolonged responses to monotherapy.
In principle, individual clonal therapies may still be used curatively if employed very early in tumor evolution, as proposed by Mitchell et al. . Clear cell kidney cancers were modeled to have deleted VHL (on chromosome 3p) several decades prior to a second mutational hit to the remaining VHL allele. Proliferation and tumorigenesis only accelerate after both alleles of VHL are mutated. Therefore, the pool of mutated cells is probably only a few hundred cells for a prolonged period. Depleting this small cell population—even marginally—with a therapy, sensitized by 3p loss, would reduce the probability of a cell with a second hit mutation ever arising. This would have to be achieved decades before these cells become detectable, and therefore would most likely involve the preemptive treatment of healthy individuals. However, considering that most tissues may harbor equivalent cell populations [8, 77, 78], such prophylactic management may not be clinically, economically, or ethically feasible.
Combined clonal therapies
Predictions of a tumor’s evolutionary response to a therapy can allow pre-emptive measures to prevent resistance. For example, ABL1 inhibition in chronic myeloid leukemia (CML), characterized by clonal BCR–ABL1 fusions, has revolutionized therapy for this disease, yet the development of resistance remains a challenge in a proportion of patients. Combining different classes of ABL1 inhibitors with mutually exclusive profiles of resistance mutations can preempt the emergence of resistant subclones (Fig. 2b). Preclinical application of this approach has resulted in durable responses .
Combining different clonal therapies might also reduce the emergence of resistance. Many breast cancers are thought to have cell-cycle dysregulation related to the cyclin-CDK-Rb pathway, in addition to estrogen sensitivity . The addition of CDK4/6 inhibition to aromatase inhibition does indeed prolong the response in patients with metastatic disease. This delays the need for conventional cytotoxic therapy, but at the price of increased toxicity compared to endocrine therapy alone .
The development of effective combination therapies requires a comprehensive understanding of mutation clonality and resistance mechanisms. Metastatic melanomas frequently have activating mutations in the MAPK pathway, and resistance to BRAF inhibitors was thought to result from downstream MEK activation [82, 83]. Trials combining MEK and BRAF inhibition in melanoma have demonstrated modest clinical benefit [84, 85]; however, resistant tumors often have multiple different detectable MAPK mutations, suggesting convergent evolution .
Ideally, larger numbers of drivers could be targeted simultaneously or sequentially, depending on the pace and nature of the evolutionary response of the tumor. Such combination therapies will impact toxicity management, although not always detrimentally. In fact, the addition of MEK inhibition to BRAF inhibitors reduces the cutaneous side-effects that are associated with BRAF inhibitors. The toxicities resulting from combination treatments may, however, require complex pharmacological adjustments that have implications for trial design, drug licensing and healthcare economic assessments.
Many tumors have only few clonal driver mutations and will require alternative strategies [87,88,89]. Effective pharmacological options for targeting driver mutations are also relatively limited. Some driver mutations may be treatable indirectly, either by collateral lethality, whereby susceptibilities created by the loss of genes adjacent to deleted tumor suppressors are harnessed, or by synthetic lethality [90, 91]. Alternatively, immunotherapy exploits the antigenicity of mutations, regardless of their driver status and without relying on the recurrence of mutations in different patients. Where durable clinical responses to immunotherapy have been seen, they are probably brought about by the simultaneous targeting of multiple clonal mutations. Indeed, one of the potential predictive markers of response to immune checkpoint blockade in non-small-cell lung cancer and melanoma is the clonal neoantigen load [92, 93]. If a common mechanism of resistance to an immunotherapy can occur, (epi) genetic variation and selection could drive tumors towards it, even when multi-pronged approaches are used. These mechanisms of immune editing are still a subject of intense study. They include an ability of tumors to reduce their antigen-presenting capability. In melanoma, lung, and ovarian cancer, these changes have been shown to result in part from either somatic (often subclonal) or germline loss of heterozygosity of the HLA locus [94,95,96]. Equivalent loss of expression of class II MHC may also result in treatment failure after allogeneic bone marrow transplant for acute myeloid leukemia . Without a full and diverse HLA repertoire, many neoantigens cannot be successfully presented on the surface of tumor cells and therefore are not recognized by an adaptive immune response.
Targeting subclonal mutations
The detection of subclonal mutations is still an active research topic and therefore potential strategies for their therapeutic use are only conceptual at present. The simplest approach is to target a combination of multiple subclonal mutations, probably coupled with a clonal therapy. In rare circumstances, such as those recently suggested in pediatric brain tumors, subclonal populations might be highly functionally interdependent . In these circumstances, even subclonal population depletion might have a profound effect on the tumor as a whole. Alternatively, if the relative importance and the clinical impact of different subclonal populations can be measured, then those causing the greatest symptomatic burden could be prioritized. Implicit in this more strategic approach is the acceptance that other cell populations that cause lower symptomatic burden will not be eradicated, representing a shift to managing cancer as a chronic disease without the intent to cure .
A combination of conventional cross-sectional imaging with the monitoring of circulating markers could be used to identify spatially or mutationally distinct metastases. If lesions are spatially segregated, they may be amenable to local therapies: surgery, cryotherapy, focused ultrasound, or stereotactic radiotherapy. If they are characterized by treatable mutations, additional systemic therapies could be used. At present, proofs of this concept are yet to emerge.
Finally, the concept of adaptive therapy has also been proposed [100, 101]. Each of the subclones present in a tumor may be either sensitive or insensitive to a potential therapy. They compete for survival within the tumor environment and a mutation that confers resistance to a treatment, possibly through the loss or alteration of an oncogenic driver, might result in a growth disadvantage when that treatment agent is not present. With an adaptive approach, sensitive subclones can be treated to the point at which tumor size is reduced or growth is suppressed to achieve symptomatic benefit. Response may conceivably be monitored with a non-invasive surrogate biomarker, such as serum prostate-specific antigen (PSA) in prostate cancer. Thereafter, treatment can be reduced or withdrawn to allow the competitive suppression of resistant subclones (Fig. 2c). This approach is currently under evaluation in metastatic prostate cancer with the use of individualized PSA thresholds to guide the use of abiraterone, a CYP17A1 inhibitor . To date, only small numbers of patients have been treated, albeit with good clinical outcome and reduced cumulative exposure to medication. It is worth noting that adaptive therapy is not the same as intermittent therapy, in which treatment may also be used discontinuously and with the monitoring of a biomarker, but without any individualization of treatment duration on the basis of response dynamics. For example, intermittent hormonal therapy has been attempted in prostate cancer. Crucially, trials such as TAP22 used fixed PSA thresholds rather than individualized thresholds [103, 104]. This could result in the depletion of treatment-sensitive clones, reducing their ability to suppress their treatment-resistant cousins.
Therapy monitoring: circulating tumor DNA and circulating tumor cells
Liquid biopsies sample more readily available body fluids, mainly blood, for cellular or genomic material that has been shed from the tumor. They are heralded for reducing the invasiveness of clinical assays used for diagnosis [105, 106], prognosis [107, 108], molecular profiling , and response assessment [110,111,112,113,114]. Monitoring the treatment of more advanced disease may be substantially enhanced by monitoring the dynamics of different tumor cell populations.
The therapeutic approaches discussed above, particularly combination subclonal targeting and adaptive therapy, rely on accurate information about the relative importance of different subclonal populations in space and time. Liquid biopsies allow non-invasive assays that can easily be repeated over time. In particular, ctDNA is relatively stable and simple to handle, and its sequence content can be analyzed using a variety of approaches  (Fig. 1).
The detection of early subclinical relapse or minimal residual disease after attempted curative therapy has relied on detecting clonal mutations in circulation. Somatic structural variants are particularly amenable to highly disease-specific PCR-based approaches. Canonical disease-defining genomic rearrangements, such as the BCR–ABL1 fusion in chronic myeloid leukemia, are routinely monitored in hematological malignancies to assess treatment response [116, 117]. Solid organ malignancies have fewer disease-defining rearrangements, but frequently possess unique somatic rearrangements that can be used to define bespoke monitoring panels [118, 119].
Monitoring of subclonal evolution has focused on evaluating somatic point mutations. Murtaza et al.  demonstrated that a dominant subclone, which was responsible for the progression of a chest wall breast cancer metastasis, was detectable by the increasing level of mutations private to that subclone. O’Leary et al. were able to use ctDNA in a small proportion of metastatic breast cancer patients, who were treated with the addition of the CDK4/6 inhibitor palbociclib, to both predict longer progression-free intervals  and detect emerging resistant subclones . Furthermore, Abbosh et al.  showed that ctDNA was detectable 10–346 days (median 70 days) prior to clinical detection of relapsed lung cancer.
There are several challenges to the adoption of this approach. Clearly, bespoke ctDNA monitoring is costly. Abbosh et al.  estimated that even a limited bespoke monitoring panel, based on detected mutations from a single primary tumor region, would cost USD 1750 per patient. In addition, current analyses have only explored minimal numbers of detectable subclones and give an incomplete picture of their number and range. Whether there are substantial biases in the tumor cells that contribute circulating DNA is currently not known. It is likely that highly vascular and necrotic tumors will contribute more to ctDNA than tumors in cryptic sites, such as the central nervous system . The use of other sources of cell-free DNA, such as stool , urine, cerebrospinal fluid, and effusions, may in part compensate for this [124, 125]. There are also likely to be genomic biases because cell-free DNA is predominantly thought to be generated by apoptotic nuclease activity which produces nucleosome-associated DNA fragments [126, 127], resulting in distinct chromatin-associated patterns. These patterns and the degree of apoptosis are likely to vary across tumor cell populations, and result in a bias in circulating tumor DNA.
The detection of subclonal mutations is also limited by the sensitivity of detection assays. Next-generation sequencing approaches that seek to gain an unbiased view of all detectable variants in circulation cannot identify rare subclonal mutations. In the Murtaza et al.  study, even clonal mutations had variant allele fractions of 3.8–34.9%. To compensate for this, most approaches, as exemplified by Abbosh et al. , use a specific amplification method based on fixed expected mutations that are detected in a sequenced primary tumor. This, by definition, means that de novo mutations that arose subsequent to the sampling of the primary tumor will not be detectable in circulation.
Circulating tumor cells (CTCs) can be analyzed using single-cell sequencing approaches. In a study by Carter et al. , the copy number profile of circulating tumor cells at the time of diagnosis of small-cell lung cancer predicted the duration of response to chemotherapy. Cellular approaches are less likely to be confounded by the genomic aberrations that arise in other cells than the index tumor [129, 130]. As a result of the rarity of these tumor cells, they require significant enrichment which is likely to introduce biases, resulting in low sensitivity even for clonal tumor populations [131, 132]. Interestingly, Kwan et al.  demonstrated that after some initial filtration, an RNA expression-based signature can be used to detect breast cancer CTCs, and that the presence of these cells carried prognostic information in the setting of neoadjuvant chemotherapy.
Can insight into tumor evolution improve early diagnosis, risk stratification, and cancer prevention?
In order to improve cancer outcomes, it is essential to alter tumor evolution. This can be achieved throughout the evolutionary timeline by preventing etiological factors, screening cell populations on the path to cancer, or stratifying cancers that will pose the greatest threat.
Cancer screening aims to reduce cancer mortality by increasing detection at a curable stage . This needs to be carefully managed, however, as overtreatment of incidental findings causes unnecessary cost, harm, and anxiety . This problem has beset the introduction of a prostate cancer screening strategy, as many low-grade prostate cancers can be managed with observation alone [136, 137]. Reliable predictive biomarkers of progression in detected lesions could increase the utility of screening programs. To date, risk stratification has relied almost exclusively on histological staging and grading.
Methods are being developed that recapitulate the early evolution of cancers using sequencing information from later-stage cancers alone, as recently reviewed . In general terms, these methods utilize the number of copies of mutations on gained chromosomal segments to infer whether these mutations happened before or after that gain. For example, if a whole chromosome has been duplicated and there are two copies of a mutation found on that chromosome, then it is likely that the mutation occurred first and was duplicated with the chromosomal gain. By analysis of whole-genome sequencing data from primary and metastatic prostate cancers, Wedge et al.  have been able to retrospectively identify chromosomal changes that developed earlier in tumorigenesis. These findings, such as the early gain of chromosome 8q, recapitulated those previously found in prostate intraepithelial neoplasia (PIN), which is thought to be a precursor of prostate adenocarcinoma .
These approaches have also been applied to invasive cancers, which have less well characterized precursor lesions [76, 89]. Recently, the Pan-Cancer Analysis of Whole Genomes (PCAWG) initiative leveraged whole-genome sequencing data to infer evolutionary timelines across cancer types . This work reproduced and refined classic models of mutational progression such as for colorectal cancer, in which APC mutations precede KRAS and TP53 mutations. This information could define mutations that can be used to risk-stratify those pre-malignant or early invasive lesions that require intervention and those that do not. In addition, large datasets and novel computational methods [141, 142] may be able to detect stereotyped evolutionary patterns and trajectories in cancer evolution that may inform early diagnosis or risk-stratification approaches.
Predicting tumor evolution—implications for risk stratification
A deeper and more comprehensive understanding of tumor evolution should allow us to understand how a cancer will behave in the future. This has specific implications for the risk stratification of established cancers. Incidental findings, such as small renal lesions that are often found during investigations for other conditions, are a clinical challenge because definitive resection is morbid but radiological and histological criteria are unreliable for prognostication . In clear cell renal cell carcinoma, Turajlic et al.  have modeled that analyses of two biopsies can allow the quantification of intra-tumor copy number heterogeneity. This can discriminate lesions of higher and lower risk of progression, thereby potentially assisting in the decision-making process for small renal lesions. In a companion study, the same authors also suggested that richer information gleaned from more thorough tumor sampling can identify evolutionary profiles that are more likely to be associated with the development of metastatic disease . In other cancer types, patterns of heterogeneity, such as copy number diversity in lung cancer  and pan-mutational diversity (so-called regional ‘explosions’) in childhood cancers , have also been shown to carry prognostic information. More transformative change to cancer prediction strategies will require the development of more complex computational tools and models [141, 142]. Much as weather forecasting models require vast amounts of measured data from the real world, cancer evolution models will require the sequence-based profiling of the evolution of many more cancers. Ultimately, this will allow these forecasts to guide the optimal management for each patient.
Prevention of key early mutagenic processes
The identification of predisposing factors for cancer, whether heritable, environmental, or infectious, has previously relied on a combination of epidemiological and biological evidence. A deeper understanding of tumor evolution can lead to new insights into the impact of these factors on the genome.
Two clear examples of direct impact on the genome are ultraviolet (UV) radiation exposure for sun-induced cancers, such as cutaneous squamous cell cancers, and exposure to tobacco smoke carcinogens for smoking-related airway cancers. The epidemiological evidence for both has long been established, although its popular acceptance took some time . Mechanisms of mutation as a result of each exposure have been identified: misrepair by transcription-coupled nucleotide excision repair of UV-induced pyrimidine photodimers  and misrepair of guanine damage by the same mechanism , respectively. These specific mutational types can now be detected across the genome as mutational signatures [70, 148], and this allows estimation of the contributions of each mutational signature (and potentially the level of mutagen exposure) in any individual tumor .
The accrual of mutations over time can now be explored retrospectively in a whole-genome-sequenced tumor. Nik-Zainal et al.  used a mutation timing approach to study changes in mutational processes over the life history of breast cancers. By leveraging the power of a large cohort of tumor samples, it becomes possible to identify mutational processes that act early or late in tumor evolution. In lung cancer, the proportion of mutations bearing a smoking signature declines later in tumor evolution, despite ongoing smoke exposure [150, 151]. Conversely, mutagenesis that is related to the activity of the APOBEC family of cytidine deaminases increases later in lung tumor evolution. As expected, inherited defects in DNA repair, such as the deficient mismatch repair seen in Lynch syndrome, can lead to steady and ongoing mutation throughout a tumor’s lifetime .
Many mutational signatures do not have identified etiologies, but direct genomic evidence can provide an objective starting point for both epidemiological and biological study. Identifying causative environmental exposures may suggest preventative measures, akin to smoking cessation and UV protection.
The challenge of somatic variation in normal tissues
The challenge in identifying mutations that are acquired early in tumorigenesis is that many canonical driver mutations, which are thought to be specific and relevant to cancer, may also occur in populations of phenotypically normal cells (Table 1).
Martincorena et al. [77, 129] identified multiple clonal expansions of cells, containing mutations in TP53, NOTCH1, and other known cancer genes, in both sun-exposed normal eyelids and in aging normal esophagus. Interestingly, mutations were much more common in NOTCH1 than in TP53 in normal esophagus, the inverse of the pattern seen in esophageal cancer, suggesting that early NOTCH1 mutations may protect against cancer development. Demeulemeester et al.  analyzed epithelial cells found in bone marrow aspirates of breast cancer patients, identifying cells with copy number aberrations that were completely distinct from the primary breast cancer and therefore from an unknown origin. Gao et al.  also detected similar aberrant cells in tissue adjacent to breast tumors that were once again unrelated to tumor cells. Finally, clonal expansions of hematopoietic cells containing leukemia-associated mutations are reported in the circulation of otherwise healthy adults [154, 155]. These confer an increased risk of the subsequent development of a hematological malignancy, but clearly many do not progress [156, 157].
In order to truly reveal the early evolution of cancer, we will need to understand the frequency of these mutational events in the normal tissues in which cancers arise. Cataloguing mutational events in normal tissues, at rare frequencies, will help to identify the cells of origin of cancer as well as the early mutational steps that occur in these cells .
Conclusions and future perspectives
Intra-tumor heterogeneity and the ability of cancers to evolve continuously has proved a major challenge to the implementation of precision anti-cancer medicine. Molecular therapies, predicted to be effective on the basis of the presence of a sensitizing mutation in a single sample, may be of limited clinical benefit. Driver mutations may be subclonal and resistance mechanisms can evolve rapidly [31, 89]. Deeper understanding of this complexity will allow the development of more robust therapeutic strategies. Without doubt, the complexity of tumor evolution is still far from being fully understood, and on an individual basis, tumors will always make unanticipated moves to evade even our best efforts. The recognition that cancer is an evolving system offers a framework on which to hang our clinical and research observations of cancer behavior and biology. We have discussed the more immediate opportunities for translating knowledge of tumor evolution here, but it seems likely that deeper insight will open additional unforeseen avenues.
Insight into the full spectrum of evolutionary paths that cancers can take may lead to the stratification of subsets of cancers that follow specific evolutionary paths. Potentially, the earliest steps or the rate-limiting steps in tumor evolution could be interrupted, either by the identification of preventable etiological factors or by timely medical interventions. These strategies may lead to a significant reduction in the incidence of some cancers or to a high cure rate in early diagnosed cancers, respectively. In addition, once diagnosed, treatment pathways may be matched according to the anticipated evolutionary path of the cancer, as opposed to classification based on traditional histological tumor subtyping. Patients with indolent tumors may be spared therapy altogether. As future therapies emerge, insight into tumor evolution is likely to inform their further development and maximize their impact. Immune checkpoint blockade is possibly the first class of therapy to emerge in this context, reaping the reward of a better understanding of the spectrum of clinical response [92, 94, 159]. Many cancers will probably need an armory of affordable, effective, and tolerable therapies that can be used safely in combination and sequentially. It is likely that conventional therapies—surgery, radiotherapy, and cytotoxic chemotherapy—will continue to have crucial roles in these treatment paradigms, but with a better understanding of the disease, these conventional therapies could be rationally combined with approaches informed by (epi) genomic insights into tumor evolution to achieve improved outcomes for cancer patients.
Circulating tumor cell
Circulating tumor DNA
Most recent common ancestor
Yeo B, Turner NC, Jones A. An update on the medical management of breast cancer. BMJ. 2014;348:g3608.
Rosell R, Karachaliou N. Optimizing lung cancer treatment approaches. Nat Rev Clin Oncol. 2014;12:75–6.
Nowell PC. The clonal evolution of tumor cell populations. Science. 1976;194:23–8.
Meyerson M, Gabriel S, Getz G. Advances in understanding cancer genomes through second-generation sequencing. Nat Rev Genet. 2010;11:685–96.
Yates LR, Campbell PJ. Evolution of the cancer genome. Nat Rev Genet. 2012;13:795–806.
Baslan T, Hicks J. Unravelling biology and shifting paradigms in cancer with single-cell sequencing. Nat Rev Cancer. 2017;17:557–69.
Lynch M. Rate, molecular spectrum, and consequences of human mutation. Proc Natl Acad Sci U S A. 2010;107:961–8.
Martincorena I, Campbell PJ. Somatic mutation in cancer and normal cells. Science. 2015;349:1483–9.
Gerlinger M, Rowan AJ, Horswell S, Larkin J, Endesfelder D, Gronroos E, et al. Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. New Engl J Med. 2012;366:883–92.
Andor N, Graham TA, Jansen M, Xia LC, Aktipis CA, Petritsch C, et al. Pan-cancer analysis of the extent and consequences of intratumor heterogeneity. Nat Med. 2016;22:105–13.
Dentro SC, Leshchiner I, Haase K, Tarabichi M, Wintersinger J, Deshwar AG, et al. Portraits of genetic intra-tumour heterogeneity and subclonal selection across cancer types. bioRxiv. 2018; https://doi.org/10.1101/312041.
McGranahan N, Swanton C. Clonal heterogeneity and tumor evolution: past, present, and the future. Cell. 2017;168:613–28.
Nik-Zainal S, Van Loo P, Wedge DC, Alexandrov LB, Greenman CD, Lau KW, et al. The life history of 21 breast cancers. Cell. 2012;149:994–1007.
Stephens PJ, Greenman CD, Fu B, Yang F, Bignell GR, Mudie LJ, et al. Massive genomic rearrangement acquired in a single catastrophic event during cancer development. Cell. 2011;144:27–40.
Alexandrov LB, Jones PH, Wedge DC, Sale JE, Campbell PJ, Nik-Zainal S, et al. Clock-like mutational processes in human somatic cells. Nat Genet. 2015;47:1402–7.
Baca SC, Prandi D, Lawrence MS, Mosquera JM, Romanel A, Drier Y, et al. Punctuated evolution of prostate cancer genomes. Cell. 2013;153:666–77.
Markowetz F. A saltationist theory of cancer evolution. Nat Genet. 2016;48:1102–3.
Nik-Zainal S, Alexandrov LB, Wedge DC, Van Loo P, Greenman CD, Raine K, et al. Mutational processes molding the genomes of 21 breast cancers. Cell. 2012;149:979–93.
Sottoriva A, Kang H, Ma Z, Graham TA, Salomon MP, Zhao J, et al. A big bang model of human colorectal tumor growth. Nat Genet. 2015;47:209–16.
Feinberg AP, Ohlsson R, Henikoff S. The epigenetic progenitor origin of human cancer. Nat Rev Genet. 2006;7:21–33.
Prensner JR, Chinnaiyan AM. The emergence of lncRNAs in cancer biology. Cancer Discov. 2011;1:391–407.
Hanash SM, Pitteri SJ, Faca VM. Mining the plasma proteome for cancer biomarkers. Nature. 2008;452:571–9.
Sheffield NC, Pierron G, Klughammer J, Datlinger P, Schönegger A, Schuster M, et al. DNA methylation heterogeneity defines a disease spectrum in Ewing sarcoma. Nat Med. 2017;23:386–95.
Li S, Garrett-Bakelman FE, Chung SS, Sanders MA, Hricik T, Rapaport F, et al. Distinct evolution and dynamics of epigenetic and genetic heterogeneity in acute myeloid leukemia. Nat Med. 2016;22:792–9.
Pan H, Jiang Y, Boi M, Tabbò F, Redmond D, Nie K, et al. Epigenomic evolution in diffuse large B-cell lymphomas. Nat Commun. 2015;6:6921.
Landau DA, Clement K, Ziller MJ, Boyle P, Fan J, Gu H, et al. Locally disordered methylation forms the basis of intratumor methylome variation in chronic lymphocytic leukemia. Cancer Cell. 2014;26:813–25.
Baker A-M, Cross W, Curtius K, Al Bakir I, Choi C-HR, Davis HL, et al. Evolutionary history of human colitis-associated colorectal cancer. Gut. 2018. https://doi.org/10.1136/gutjnl-2018-31619.
Angelova M, Mlecnik B, Vasaturo A, Bindea G, Fredriksen T, Lafontaine L, et al. Evolution of metastases in space and time under immune selection. Cell. 2018;175:751–65.
Cooper CS, Eeles R, Wedge DC, Van Loo P, Gundem G, Alexandrov LB, et al. Analysis of the genetic phylogeny of multifocal prostate cancer identifies multiple independent clonal expansions in neoplastic and morphologically normal prostate tissue. Nat Genet. 2015;47:367–72.
Abbosh C, Birkbak NJ, Wilson GA, Jamal-Hanjani M, Constantin T, Salari R, et al. Phylogenetic ctDNA analysis depicts early-stage lung cancer evolution. Nature. 2017;545:446–51.
Gundem G, Van Loo P, Kremeyer B, Alexandrov LB, Tubio JM, Papaemmanuil E, et al. The evolutionary history of lethal metastatic prostate cancer. Nature. 2015;520:353–7.
Campbell PJ, Yachida S, Mudie LJ, Stephens PJ, Pleasance ED, Stebbings LA, et al. The patterns and dynamics of genomic instability in metastatic pancreatic cancer. Nature. 2010;467:1109–13.
Yates LR, Knappskog S, Wedge D, Farmery JHR, Gonzalez S, Martincorena I, et al. Genomic evolution of breast cancer metastasis and relapse. Cancer Cell. 2017;32:169–84.
Landau DA, Carter SL, Stojanov P, McKenna A, Stevenson K, Lawrence MS, et al. Evolution and impact of subclonal mutations in chronic lymphocytic leukemia. Cell. 2013;152:714–26.
Rocco JW. Mutant allele tumor heterogeneity (MATH) and head and neck squamous cell carcinoma. Head Neck Pathol. 2015;9:1–5.
Yates LR, Gerstung M, Knappskog S, Desmedt C, Gundem G, Van Loo P, et al. Subclonal diversification of primary breast cancer revealed by multiregion sequencing. Nat Med. 2015;21:751–9.
Fearon ER, Vogelstein B. A genetic model for colorectal tumorigenesis. Cell. 1990;61:759–67.
Shain AH, Joseph NM, Yu R, Benhamida J, Liu S, Prow T, et al. Genomic and transcriptomic analysis reveals incremental disruption of key signaling pathways during melanoma evolution. Cancer Cell. 2018;34:45–55.
Vogelstein B, Fearon ER, Hamilton SR, Kern SE, Preisinger AC, Leppert M, et al. Genetic alterations during colorectal-tumor development. New Engl J Med. 1988;319:525–32.
Teixeira VH, Pipinikas CP, Pennycuick A, Lee-Six H, Chandrasekharan D, Beane J, et al. Deciphering the genomic, epigenomic, and transcriptomic landscapes of pre-invasive lung cancer lesions. Nat Med. 2019;25:517–25.
Dentro SC, Wedge DC, Van Loo P. Principles of reconstructing the subclonal architecture of cancers. Cold Spring Harb Perspect Med. 2017;7:a026625.
Jolly C, Van Loo P. Timing somatic events in the evolution of cancer. Genome Biol. 2018;19:95.
Ding L, Ley TJ, Larson DE, Miller CA, Koboldt DC, Welch JS, et al. Clonal evolution in relapsed acute myeloid leukaemia revealed by whole-genome sequencing. Nature. 2012;481:506–10.
Cortes JE, Talpaz M, Giles F, O'Brien S, Rios MB, Shan J, et al. Prognostic significance of cytogenetic clonal evolution in patients with chronic myelogenous leukemia on imatinib mesylate therapy. Blood. 2003;101:3794–800.
Jan M, Snyder TM, Corces-Zimmerman MR, Vyas P, Weissman IL, Quake SR, et al. Clonal evolution of preleukemic hematopoietic stem cells precedes human acute myeloid leukemia. Sci Transl Med. 2012;4:149ra118.
Schuh A, Becq J, Humphray S, Alexa A, Burns A, Clifford R, et al. Monitoring chronic lymphocytic leukemia progression by whole genome sequencing reveals heterogeneous clonal evolution patterns. Blood. 2012;120:4191–6.
Wang Y, Navin NE. Advances and applications of single-cell sequencing technologies. Mol Cell. 2015;58:598–609.
Navin N, Hicks J. Future medical applications of single-cell sequencing in cancer. Genome Med. 2011;3:31.
Van Loo P, Voet T. Single cell analysis of cancer genomes. Curr Opin Genet Dev. 2014;24:82–91.
Navin NE. Cancer genomics: one cell at a time. Genome Biol. 2014;15:452.
Zheng GX, Terry JM, Belgrader P, Ryvkin P, Bent ZW, Wilson R, et al. Massively parallel digital transcriptional profiling of single cells. Nat Commun. 2017;8:14049.
Kim C, Gao R, Sei E, Brandt R, Hartman J, Hatschek T, et al. Chemoresistance evolution in triple-negative breast cancer delineated by single-cell sequencing. Cell. 2018;173:879–93.
Zahn H, Steif A, Laks E, Eirew P, VanInsberghe M, Shah SP, et al. Scalable whole-genome single-cell library preparation without preamplification. Nat Methods. 2017;14:167–73.
Macaulay IC, Haerty W, Kumar P, Li YI, Hu TX, Teng MJ, et al. G&T-seq: parallel sequencing of single-cell genomes and transcriptomes. Nat Methods. 2015;12:519–22.
Macaulay IC, Ponting CP, Voet T. Single-cell multiomics: multiple measurements from single cells. Trends Genet. 2017;33:155–68.
Dey SS, Kester L, Spanjaard B, Bienko M, van Oudenaarden A. Integrated genome and transcriptome sequencing of the same cell. Nat Biotechnol. 2015;33:285–9.
Ryser MD, Yu M, Grady W, Siegmund K, Shibata D. Epigenetic heterogeneity in human colorectal tumors reveals preferential conservation and evidence of immune surveillance. Sci Rep. 2018;8:17292.
Bian S, Hou Y, Zhou X, Li X, Yong J, Wang Y, et al. Single-cell multiomics sequencing and analyses of human colorectal cancer. Science. 2018;362:1060–3.
Maley CC, Aktipis A, Graham TA, Sottoriva A, Boddy AM, Janiszewska M, et al. Classifying the evolutionary and ecological features of neoplasms. Nat Rev Cancer. 2017;17:605–19.
Davis C, Naci H, Gurpinar E, Poplavska E, Pinto A, Aggarwal A. Availability of evidence of benefits on overall survival and quality of life of cancer drugs approved by European medicines agency: retrospective cohort study of drug approvals 2009–13. BMJ. 2017;359:j4530.
McArthur GA, Chapman PB, Robert C, Larkin J, Haanen JB, Dummer R, et al. Safety and efficacy of vemurafenib in BRAF(V600E) and BRAF(V600K) mutation-positive melanoma (BRIM-3): extended follow-up of a phase 3, randomised, open-label study. Lancet Oncol. 2014;15:323–32.
Rosell R, Carcereny E, Gervais R, Vergnenegre A, Massuti B, Felip E, et al. Erlotinib versus standard chemotherapy as first-line treatment for European patients with advanced EGFR mutation-positive non-small-cell lung cancer (EURTAC): a multicentre, open-label, randomised phase 3 trial. Lancet Oncol. 2012;13:239–46.
Jeselsohn R, Buchwalter G, De Angelis C, Brown M, Schiff R. ESR1 mutations—a mechanism for acquired endocrine resistance in breast cancer. Nat Rev Clin Oncol. 2015;12:573–83.
Watson PA, Arora VK, Sawyers CL. Emerging mechanisms of resistance to androgen receptor inhibitors in prostate cancer. Nat Rev Cancer. 2015;15:701–11.
Lord CJ, Ashworth A. BRCAness revisited. Nat Rev Cancer. 2016;16:110–20.
Ashworth A. A synthetic lethal therapeutic approach: poly (ADP) ribose polymerase inhibitors for the treatment of cancers deficient in DNA double-strand break repair. J Clin Oncol. 2008;26:3785–90.
Patch AM, Christie EL, Etemadmoghadam D, Garsed DW, George J, Fereday S, et al. Whole-genome characterization of chemoresistant ovarian cancer. Nature. 2015;521:489–94.
Weigelt B, Comino-Mendez I, de Bruijn I, Tian L, Meisel JL, Garcia-Murillas I, et al. Diverse BRCA1 and BRCA2 reversion mutations in circulating cell-free DNA of therapy-resistant breast or ovarian cancer. Clin Cancer Res. 2017;23:6708–20.
Christie EL, Fereday S, Doig K, Pattnaik S, Dawson SJ, Bowtell DDL. Reversion of BRCA1/2 germline mutations detected in circulating tumor DNA from patients with high-grade serous ovarian cancer. J Clin Oncol. 2017;35:1274–80.
Alexandrov LB, Nik-Zainal S, Wedge DC, Aparicio SA, Behjati S, Biankin AV, et al. Signatures of mutational processes in human cancer. Nature. 2013;500:415–21.
Szikriszt B, Póti Á, Pipek O, Krzystanek M, Kanu N, Molnár J, et al. A comprehensive survey of the mutagenic impact of common cancer cytotoxics. Genome Biol. 2016;17:99.
Sharma A, Cao EY, Kumar V, Zhang X, Leong HS, Wong AML, et al. Longitudinal single-cell RNA sequencing of patient-derived primary cells reveals drug-induced infidelity in stem cell hierarchy. Nat Commun. 2018;9:4931.
Shah KN, Bhatt R, Rotow J, Rohrberg J, Olivas V, Wang VE, et al. Aurora kinase a drives the evolution of resistance to third-generation EGFR inhibitors in lung cancer. Nat Med. 2019;25:111–8.
Glickman MS, Sawyers CL. Converting cancer therapies into cures: lessons from infectious diseases. Cell. 2012;148:1089–98.
Stewart PS, Franklin MJ. Physiological heterogeneity in biofilms. Nat Rev Microbiol. 2008;6:199–210.
Mitchell TJ, Turajlic S, Rowan A, Nicol D, Farmery JHR, O'Brien T, et al. Timing the landmark events in the evolution of clear cell renal cell cancer: TRACERx renal. Cell. 2018;173:611–23.
Martincorena I, Fowler JC, Wabik A, Lawson ARJ, Abascal F, Hall MWJ, et al. Somatic mutant clones colonize the human esophagus with age. Science. 2018;362:911–7.
Yokoyama A, Kakiuchi N, Yoshizato T, Nannya Y, Suzuki H, Takeuchi Y, et al. Age-related remodelling of oesophageal epithelia by mutated cancer drivers. Nature. 2019;565:312–7.
Wylie AA, Schoepfer J, Jahnke W, Cowan-Jacob SW, Loo A, Furet P, et al. The allosteric inhibitor ABL001 enables dual targeting of BCR–ABL1. Nature. 2017;543:733–7.
Dean JL, Thangavel C, McClendon AK, Reed CA, Knudsen ES. Therapeutic CDK4/6 inhibition in breast cancer: key mechanisms of response and failure. Oncogene. 2010;29:4018–32.
Cristofanilli M, Turner NC, Bondarenko I, Ro J, Im SA, Masuda N, et al. Fulvestrant plus palbociclib versus fulvestrant plus placebo for treatment of hormone-receptor-positive, HER2-negative metastatic breast cancer that progressed on previous endocrine therapy (PALOMA-3): final analysis of the multicentre, double-blind, phase 3 randomised controlled trial. Lancet Oncol. 2016;17:425–39.
Nazarian R, Shi H, Wang Q, Kong X, Koya RC, Lee H, et al. Melanomas acquire resistance to B-RAF(V600E) inhibition by RTK or N-RAS upregulation. Nature. 2010;468:973–7.
Johannessen CM, Boehm JS, Kim SY, Thomas SR, Wardwell L, Johnson LA, et al. COT drives resistance to RAF inhibition through MAP kinase pathway reactivation. Nature. 2010;468:968–72.
Long GV, Stroyakovskiy D, Gogas H, Levchenko E, de Braud F, Larkin J, et al. Combined BRAF and MEK inhibition versus BRAF inhibition alone in melanoma. New Engl J Med. 2014;371:1877–88.
Moriceau G, Hugo W, Hong A, Shi H, Kong X, Yu CC, et al. Tunable-combinatorial mechanisms of acquired resistance limit the efficacy of BRAF/MEK cotargeting but result in melanoma drug addiction. Cancer Cell. 2015;27:240–56.
Wagle N, Van Allen EM, Treacy DJ, Frederick DT, Cooper ZA, Taylor-Weiner A, et al. MAP kinase pathway alterations in BRAF-mutant melanoma patients with acquired resistance to combined RAF/MEK inhibition. Cancer Discov. 2014;4:61–8.
Turajlic S, Xu H, Litchfield K, Rowan A, Horswell S, Chambers T, et al. Deterministic evolutionary trajectories influence primary tumor growth: TRACERx renal. Cell. 2018;173:595–610.
Turajlic S, Xu H, Litchfield K, Rowan A, Chambers T, Lopez JI, et al. Tracking cancer evolution reveals constrained routes to metastases: TRACERx renal. Cell. 2018;173:581–94.
Jamal-Hanjani M, Wilson GA, McGranahan N, Birkbak NJ, Watkins TBK, Veeriah S, et al. Tracking the evolution of non–small-cell lung cancer. New Engl J Med. 2017;376:2109–21.
Muller FL, Aquilanti EA, DePinho RA. Collateral lethality: a new therapeutic strategy in oncology. Trends Cancer. 2015;1:161–73.
Muller FL, Colla S, Aquilanti E, Manzo VE, Genovese G, Lee J, et al. Passenger deletions generate therapeutic vulnerabilities in cancer. Nature. 2012;488:337–42.
McGranahan N, Furness AJ, Rosenthal R, Ramskov S, Lyngaa R, Saini SK, et al. Clonal neoantigens elicit T cell immunoreactivity and sensitivity to immune checkpoint blockade. Science. 2016;351:1463–9.
Snyder A, Makarov V, Merghoub T, Yuan J, Zaretsky JM, Desrichard A, et al. Genetic basis for clinical response to CTLA-4 blockade in melanoma. New Engl J Med. 2014;371:2189–99.
McGranahan N, Rosenthal R, Hiley CT, Rowan AJ, Watkins TBK, Wilson GA, et al. Allele-specific HLA loss and immune escape in lung cancer evolution. Cell. 2017;171:1259–71.
Chowell D, Morris LGT, Grigg CM, Weber JK, Samstein RM, Makarov V, et al. Patient HLA class I genotype influences cancer response to checkpoint blockade immunotherapy. Science. 2018;359:582–7.
Zhang AW, McPherson A, Milne K, Kroeger DR, Hamilton PT, Miranda A, et al. Interfaces of malignant and immunologic clonal dynamics in ovarian cancer. Cell. 2018;173:1755–69.
Christopher MJ, Petti AA, Rettig MP, Miller CA, Chendamarai E, Duncavage EJ, et al. Immune escape of relapsed AML cells after allogeneic transplantation. New Engl J Med. 2018;379:2330–41.
Vinci M, Burford A, Molinari V, Kessler K, Popov S, Clarke M, et al. Functional diversity and cooperativity between subclonal populations of pediatric glioblastoma and diffuse intrinsic pontine glioma cells. Nat Med. 2018;24:1204–15.
Beck S, Ng T. C2c: turning cancer into chronic disease. Genome Med. 2014;6:38.
Enriquez-Navas PM, Kam Y, Das T, Hassan S, Silva A, Foroutan P, et al. Exploiting evolutionary principles to prolong tumor control in preclinical models of breast cancer. Sci Transl Med. 2016;8:327ra24.
Gatenby RA, Silva AS, Gillies RJ, Frieden BR. Adaptive therapy. Cancer Res. 2009;69:4894–903.
Zhang J, Cunningham JJ, Brown JS, Gatenby RA. Integrating evolutionary dynamics into treatment of metastatic castrate-resistant prostate cancer. Nat Commun. 2017;8:1816.
Mottet N, Van Damme J, Loulidi S, Russel C, Leitenberger A, Wolff JM, et al. Intermittent hormonal therapy in the treatment of metastatic prostate cancer: a randomized trial. BJU Int. 2012;110:1262–9.
Magnan S, Zarychanski R, Pilote L, Bernier L, Shemilt M, Vigneault E, et al. Intermittent vs continuous androgen deprivation therapy for prostate cancer: a systematic review and meta-analysis. JAMA Oncol. 2015;1:1261–9.
Newman AM, Bratman SV, To J, Wynne JF, Eclov NC, Modlin LA, et al. An ultrasensitive method for quantitating circulating tumor DNA with broad patient coverage. Nat Med. 2014;20:548–54.
Beaver JA, Jelovac D, Balukrishna S, Cochran R, Croessmann S, Zabransky DJ, et al. Detection of cancer DNA in plasma of patients with early-stage breast cancer. Clin Cancer Res. 2014;20:2643–50.
Bettegowda C, Sausen M, Leary RJ, Kinde I, Wang Y, Agrawal N, et al. Detection of circulating tumor DNA in early- and late-stage human malignancies. Sci Transl Med. 2014;6:224ra24.
Khan KH, Cunningham D, Werner B, Vlachogiannis G, Spiteri I, Heide T, et al. Longitudinal liquid biopsy and mathematical modeling of clonal evolution forecast time to treatment failure in the PROSPECT-C phase II colorectal Cancer clinical trial. Cancer Discov. 2018;8:1270–85.
Oxnard GR, Thress KS, Alden RS, Lawrance R, Paweletz CP, Cantarini M, et al. Association between plasma genotyping and outcomes of treatment with osimertinib (AZD9291) in advanced non-small-cell lung cancer. J Clin Oncol. 2016;34:3375–82.
Dawson SJ, Tsui DW, Murtaza M, Biggs H, Rueda OM, Chin SF, et al. Analysis of circulating tumor DNA to monitor metastatic breast cancer. New Engl J Med. 2013;368:1199–209.
O'Leary B, Hrebien S, Morden JP, Beaney M, Fribbens C, Huang X, et al. Early circulating tumor DNA dynamics and clonal selection with palbociclib and fulvestrant for breast cancer. Nat Commun. 2018;9:896.
Spina V, Bruscaggin A, Cuccaro A, Martini M, Di Trani M, Forestieri G, et al. Circulating tumor DNA reveals genetics, clonal evolution, and residual disease in classical Hodgkin lymphoma. Blood. 2018;131:2413–25.
Barault L, Amatu A, Siravegna G, Ponzetti A, Moran S, Cassingena A, et al. Discovery of methylated circulating DNA biomarkers for comprehensive non-invasive monitoring of treatment response in metastatic colorectal cancer. Gut. 2018;67:1995–2005.
Agarwal R, Chan Y-C, Tam CS, Hunter T, Vassiliadis D, Teh CE, et al. Dynamic molecular monitoring reveals that SWI–SNF mutations mediate resistance to ibrutinib plus venetoclax in mantle cell lymphoma. Nat Med. 2019;25:119–29.
Henao Diaz E, Yachnin J, Gronberg H, Lindberg J. The in vitro stability of circulating tumour DNA. PLoS One. 2016;11:e0168153.
Marin D, Kaeda J, Szydlo R, Saunders S, Fleming A, Howard J, et al. Monitoring patients in complete cytogenetic remission after treatment of CML in chronic phase with imatinib: patterns of residual leukaemia and prognostic factors for cytogenetic relapse. Leukemia. 2005;19:507–12.
Branford S, Rudzki Z, Parkinson I, Grigg A, Taylor K, Seymour JF, et al. Real-time quantitative PCR analysis can be used as a primary screen to identify patients with CML treated with imatinib who have BCR-ABL kinase domain mutations. Blood. 2004;104:2926–32.
McBride DJ, Orpana AK, Sotiriou C, Joensuu H, Stephens PJ, Mudie LJ, et al. Use of cancer-specific genomic rearrangements to quantify disease burden in plasma from patients with solid tumors. Genes Chromosomes Cancer. 2010;49:1062–9.
Olsson E, Winter C, George A, Chen Y, Howlin J, Tang MH, et al. Serial monitoring of circulating tumor DNA in patients with primary breast cancer for detection of occult metastatic disease. EMBO Mol Med. 2015;7:1034–47.
Murtaza M, Dawson SJ, Pogrebniak K, Rueda OM, Provenzano E, Grant J, et al. Multifocal clonal evolution characterized using circulating tumour DNA in a case of metastatic breast cancer. Nat Commun. 2015;6:8760.
O'Leary B, Cutts RJ, Liu Y, Hrebien S, Huang X, Fenwick K, et al. The genetic landscape and clonal evolution of breast cancer resistance to palbociclib plus fulvestrant in the PALOMA-3 trial. Cancer Discov. 2018;8:1390–403.
De Mattos-Arruda L, Mayor R, Ng CK, Weigelt B, Martinez-Ricarte F, Torrejon D, et al. Cerebrospinal fluid-derived circulating tumour DNA better represents the genomic alterations of brain tumours than plasma. Nat Commun. 2015;6:8839.
Imperiale TF, Ransohoff DF, Itzkowitz SH. Multitarget stool DNA testing for colorectal-cancer screening. New Engl J Med. 2014;371:187–8.
Peng M, Chen C, Hulbert A, Brock MV, Yu F. Non-blood circulating tumor DNA detection in cancer. Oncotarget. 2017;8:69162–73.
Miller AM, Shah RH, Pentsova EI, Pourmaleki M, Briggs S, Distefano N, et al. Tracking tumour evolution in glioma through liquid biopsies of cerebrospinal fluid. Nature. 2019;565:654–8.
Thierry AR, Mouliere F, Gongora C, Ollier J, Robert B, Ychou M, et al. Origin and quantification of circulating DNA in mice with human colorectal cancer xenografts. Nucleic Acids Res. 2010;38:6159–75.
Giacona MB, Ruben GC, Iczkowski KA, Roos TB, Porter DM, Sorenson GD. Cell-free DNA in human blood plasma: length measurements in patients with pancreatic cancer and healthy controls. Pancreas. 1998;17:89–97.
Carter L, Rothwell DG, Mesquita B, Smowton C, Leong HS, Fernandez-Gutierrez F, et al. Molecular analysis of circulating tumor cells identifies distinct copy-number profiles in patients with chemosensitive and chemorefractory small-cell lung cancer. Nat Med. 2017;23:114–9.
Martincorena I, Roshan A, Gerstung M, Ellis P, Van Loo P, McLaren S, et al. Tumor evolution. High burden and pervasive positive selection of somatic mutations in normal human skin. Science. 2015;348:880–6.
Demeulemeester J, Kumar P, Moller EK, Nord S, Wedge DC, Peterson A, et al. Tracing the origin of disseminated tumor cells in breast cancer using single-cell sequencing. Genome Biol. 2016;17:250.
Alix-Panabieres C, Pantel K. Challenges in circulating tumour cell research. Nat Rev Cancer. 2014;14:623–31.
Krebs MG, Metcalf RL, Carter L, Brady G, Blackhall FH, Dive C. Molecular analysis of circulating tumour cells-biology and biomarkers. Nat Rev Clin Oncol. 2014;11:129–44.
Kwan TT, Bardia A, Spring LM, Giobbie-Hurder A, Kalinich M, Dubash T, et al. A digital RNA signature of circulating tumor cells predicting early therapeutic response in localized and metastatic breast cancer. Cancer Discov. 2018;8:1286–99.
Gray JA. New concepts in screening. Br J Gen Pract. 2004;54:292–8.
Esserman LJ, Thompson IM Jr, Reid B. Overdiagnosis and overtreatment in cancer: an opportunity for improvement. JAMA. 2013;310:797–8.
Eggener SE, Cifu AS, Nabhan C. Prostate cancer screening. JAMA. 2015;314:825–6.
Hugosson J, Carlsson S, Aus G, Bergdahl S, Khatami A, Lodding P, et al. Mortality results from the Goteborg randomised population-based prostate-cancer screening trial. Lancet Oncol. 2010;11:725–32.
Wedge DC, Gundem G, Mitchell T, Woodcock DJ, Martincorena I, Ghori M, et al. Sequencing of prostate cancers identifies new cancer genes, routes of progression and drug targets. Nat Genet. 2018;50:682–92.
Jung SH, Shin S, Kim MS, Baek IP, Lee JY, Lee SH, et al. Genetic progression of high grade prostatic intraepithelial neoplasia to prostate cancer. Eur Urol. 2016;69:823–30.
Gerstung M, Jolly C, Leshchiner I, Dentro SC, Gonzalez S, Mitchell TJ, et al. The evolutionary history of 2,658 cancers. bioRxiv. 2018; doi: https://doi.org/10.1101/161562.
Caravagna G, Giarratano Y, Ramazzotti D, Tomlinson I, Graham TA, Sanguinetti G, et al. Detecting repeated cancer evolution from multi-region tumor sequencing data. Nat Methods. 2018;15:707–14.
Gerhauser C, Favero F, Risch T, Simon R, Feuerbach L, Assenov Y, et al. Molecular evolution of early-onset prostate cancer identifies molecular risk markers and clinical trajectories. Cancer Cell. 2018;34:996–1011.
Ljungberg B, Bensalah K, Canfield S, Dabestani S, Hofmann F, Hora M, et al. EAU guidelines on renal cell carcinoma: 2014 update. Eur Urol. 2015;67:913–24.
Karlsson J, Valind A, Holmquist Mengelbier L, Bredin S, Cornmark L, Jansson C, et al. Four evolutionary trajectories underlie genetic intratumoral variation in childhood cancer. Nat Genet. 2018;50:944–50.
Doll R, Hill AB. Smoking and carcinoma of the lung; preliminary report. BMJ. 1950;2:739–48.
Brash DE, Haseltine WA. UV-induced mutation hotspots occur at DNA damage hotspots. Nature. 1982;298:189–92.
Denissenko MF, Pao A, Tang M, Pfeifer GP. Preferential formation of benzo [a] pyrene adducts at lung cancer mutational hotspots in P53. Science. 1996;274:430–2.
Alexandrov LB, Ju YS, Haase K, Van Loo P, Martincorena I, Nik-Zainal S, et al. Mutational signatures associated with tobacco smoking in human cancer. Science. 2016;354:618–22.
Nik-Zainal S, Davies H, Staaf J, Ramakrishna M, Glodzik D, Zou X, et al. Landscape of somatic mutations in 560 breast cancer whole-genome sequences. Nature. 2016;534:47–54.
Swanton C, McGranahan N, Starrett GJ, Harris RS. APOBEC enzymes: mutagenic fuel for cancer evolution and heterogeneity. Cancer Discov. 2015;5:704–12.
de Bruin EC, McGranahan N, Mitter R, Salm M, Wedge DC, Yates L, et al. Spatial and temporal diversity in genomic instability processes defines lung cancer evolution. Science. 2014;346:251–6.
Campbell BB, Light N, Fabrizio D, Zatzman M, Fuligni F, De Borja R, et al. Comprehensive analysis of hypermutation in human cancer. Cell. 2017;171:1042–56.
Gao R, Davis A, McDonald TO, Sei E, Shi X, Wang Y, et al. Punctuated copy number evolution and clonal stasis in triple-negative breast cancer. Nat Genet. 2016;48:1119–30.
Xie M, Lu C, Wang J, McLellan MD, Johnson KJ, Wendl MC, et al. Age-related mutations associated with clonal hematopoietic expansion and malignancies. Nat Med. 2014;20:1472–8.
Genovese G, Kahler AK, Handsaker RE, Lindberg J, Rose SA, Bakhoum SF, et al. Clonal hematopoiesis and blood-cancer risk inferred from blood DNA sequence. New Engl J Med. 2014;371:2477–87.
Welch JS, Ley TJ, Link DC, Miller CA, Larson DE, Koboldt DC, et al. The origin and evolution of mutations in acute myeloid leukemia. Cell. 2012;150:264–78.
Gibson CJ, Steensma DP. New insights from studies of clonal hematopoiesis. Clin Cancer Res. 2018;24:4633–42.
Young MD, Mitchell TJ, Vieira Braga FA, Tran MGB, Stewart BJ, Ferdinand JR, et al. Single-cell transcriptomes from human kidneys reveal the cellular identity of renal tumors. Science. 2018;361:594–9.
Snyder A, Nathanson T, Funt SA, Ahuja A, Buros Novik J, Hellmann MD, et al. Contribution of systemic and somatic factors to clinical response and resistance to PD-L1 blockade in urothelial cancer: an exploratory multi-omic analysis. PLoS Med. 2017;14:e1002309.
Research in the PVL lab is supported by the Francis Crick Institute, which receives its core funding from Cancer Research UK (FC001202), the UK Medical Research Council (FC001202), and the Wellcome Trust (FC001202).
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.