Deciphering intratumor heterogeneity and temporal acquisition of driver events to refine precision medicine
© Hiley et al.; licensee BioMed Central 2014
Published: 27 August 2014
Skip to main content
© Hiley et al.; licensee BioMed Central 2014
Published: 27 August 2014
The presence of multiple subclones within tumors mandates understanding of longitudinal and spatial subclonal dynamics. Resolving the spatial and temporal heterogeneity of subclones with cancer driver events may offer insight into therapy response, tumor evolutionary histories and clinical trial design.
The identification of somatic events and mutational processes that drive a cancer is increasingly important for precision cancer diagnosis and therapy. To date, sequencing efforts have identified several hundred cancer-driver mutations and genomic aberrations across multiple cancer types [1–4]. Sequencing studies have also shed light on the extent of tumor diversity, not only among tumors from different patients (intertumor heterogeneity) but also within individual tumors (intratumor heterogeneity). Pathologists have long recognized heterogeneity within tumors at the morphological level, and heterogeneity at the genetic level was first shown several decades ago by cytogenetic analyses (as reviewed by Navin and Hicks ), but more recent sequencing studies have provided deeper insights into the full extent of intertumor and intratumor heterogeneity. It is increasingly recognized that tumors consist of multiple genetically distinct subclones that often evolve following a pattern of branched evolution.
There is a need, therefore, not only to determine which driver events occur in a tumor but also to understand their relative timing during tumor evolution. Moreover, our understanding of how changes in the prevalence of different subclones over time impact upon therapeutic response and clinical outcome remains limited. There is an increased realization of the need to understand a tumor’s evolutionary history using both spatial and longitudinal genomic information and to identify driver events and mutational processes that contribute to tumor initiation, maintenance, progression and subclonal diversification. Here, we summarize recent findings on the relevance of subclonal driver events. We also describe how subclonal diversity might contribute to the limitations of targeted therapies and how it can be leveraged to study the evolutionary history of a tumor and to optimize clinical trial design.
The subclonal diversity within a tumor if viewed as a snapshot, rather than longitudinally, provides little information about the future evolutionary routes that subclonal populations might take. Has the dominant subclone within a tumor outcompeted less fit minor subclones or is a new fitter subclonal population emerging (Figure 1)? A greater understanding of the evolutionary timings and ‘life histories’ of tumors might shed light on the most clinically significant subclones and reveal common rules that govern tumor evolution both within and across cancer subtypes.
In breast and pancreatic cancers, the majority of known mutational and copy number driver events are relatively early events [15,16,27,29]. In ccRCC, however, the majority of identified driver mutations were found to be subclonal. In fact, inactivating mutations in the Von Hippel-Lindau tumor suppressor gene and loss of heterozygosity at chromosome 3p were the only somatic events occurring ‘early’ in ten ccRCC tumors analyzed . Conversely, mutations in certain driver genes, including BAP1 (BRCA1 associated protein-1), PTEN (phosphatase and tensin homolog), PIK3CA (phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha), SETD2 (SET domain containing 2) and TP53 (encoding tumor protein p53), were always subclonal, and thus probably involved in ccRCC progression. Selection for mutations in these genes during later stages of ccRCC development is evidenced by the observation that different subclones acquire mutations in the same gene in parallel.
Driver events are not always clonally dominant across all cancer subtypes. For example, loss of the tumor suppressor PTEN has been identified, by sequencing and fluorescence in situ hybridization analyses, as a subclonal event in prostate tumors but often as a clonal event in triple negative breast cancer [14,29,30]. Similarly, TP53 mutations were usually identified as early events in triple negative breast tumors and cutaneous squamous cell carcinomas [28,29] but predominantly as subclonal in CLL and ccRCC [18,22]. Nevertheless, to define a clonal mutation truly as an initiating event, single cell analysis would be required to reconstruct the evolutionary lineage of a sample. For example, the etiological significance in MM of the chromosomal translocation t(11;14) was further evidenced by the fact that this was the sole abnormality in some subclones and that genetic variegation of this mutation in an initial clone had resulted in the tumor heterogeneity .
Deciphering genetic intratumor heterogeneity sheds light not only on the temporal acquisition of somatic driver events, but also on the temporal dynamics of mutational processes. Analysis of breast cancers found that early mutations were dominated by C-to-T transitions, predominantly in a CpG context, probably reflecting spontaneous deamination of methylated cytosines . By contrast, in ccRCC, C-to-T transitions at CpG sites were enriched in late mutations . Later in breast cancer evolution, novel mutational processes, such as mutagenesis by APOBEC (apolipoprotein B mRNA-editing, enzyme-catalytic, polypeptide-like 3G) cytidine deaminases, were found to dominate in certain samples . Intriguingly, in MM, the contribution of APOBEC-edited mutagenesis was found to either increase or decrease over time, depending on the sample , whereas in bladder cancer, the contribution of APOBEC remained relatively stable between pairs of superficial noninvasive and mucosa- or muscle-invasive tumors from two patients .
Taken together, these studies highlight the diversity in tumor evolutionary processes, with variation in both the temporal acquisition of driver mutations and the mutational signatures themselves observed both within and between tumor types. As more samples are analyzed, it will be important to validate which driver genes are always clonal (that is, on the trunk of a tumor’s evolutionary tree) and which may occur later in tumor evolution, driving subclonal expansions. Early drivers may serve as optimal therapeutic targets in future drug-development strategies . Importantly, whereas a clonal driver mutation in a tumor is fixed and will remain clonal unless subject to copy number loss later in tumor evolution, a subclonal driver mutation may become clonal, evade detection or disappear entirely at a later stage of tumor evolution. Thus, a subclonal driver is more dynamic than a clonal driver. Longitudinal studies are needed to shed light on whether there are epistatic relationships between driver events and to explore the possibility of sets of evolutionary rules that determine how tumor cell populations change over time within each cancer type. The existence of such evolutionary rules is supported by evidence of the parallel evolution of distinct subclones within the same tumor, each harboring distinct somatic events that affect the same gene or signal transduction pathway, and by pairwise associations between different driver events [22,31,34]. Conceivably, if rules regarding the temporal acquisition of somatic driver events can be defined - taking into account the tumor microenvironment, the host genome and early somatic events in tumor evolution - they could inform therapeutic strategies. Similarly, a greater understanding of whether mutational processes, such as APOBEC-mediated mutagenesis, generally occur transiently or accumulate gradually over time may shed light on how specific cancers should be monitored and treated .
Heterogeneity between primary and metastatic lesions has profound implications for approaches to genomic research and patient care, as does the heterogeneity of clones within a single sample that changes over time because of cell-intrinsic mechanisms such as genomic instability or selective pressures from tumor-directed therapy.
Spatial heterogeneity in solid tumors can result in significant sampling bias. In high-grade serous ovarian cancer, multi-region sampling of six patients prior to treatment demonstrated the diversity of somatic mutation, copy number and gene expression within each patient . Other than TP53, few driver genes were ubiquitous in the multiple sampled regions from each patient. In ccRCC, multiple biopsies are required to better define the true extent of genomic heterogeneity and clinically relevant mutations. For example, mutations in the mammalian target of rapamycin kinase that confer resistance to everolimus can be found in some but not all tumor regions [22,23].
Sequencing studies have demonstrated heterogeneity of driver events between primary tumor and metastatic sites. Similar heterogeneity has also been demonstrated for clinically relevant biomarkers [12,16,23]. A retrospective review of HER2 (v-erb-b2 avian erythroblastic leukemia viral oncogene homolog 2) status in primary breast cancer and metastatic relapses showed significant discordance . Patients with stage IV breast cancer who received HER2-directed therapy due to having had HER2-positive early breast cancer and who on retrospective analysis were found to have a HER2-negative metastatic relapse had shorter overall survival compared to those who had a true HER2-positive metastatic relapse. Sequencing of recurrent high-grade gliomas, after earlier surgical resection and sequencing of low grade lesions, showed that almost half of the high-grade relapses did not come from the previously resected low-grade glioma but from an ancestral clone that predated the low-grade component . Many potentially therapeutically targetable driver mutations, such as BRAF (B-Raf proto-oncogene, serine/threonine kinase) V600E, that were present in the primary low-grade lesion, were not present in the recurrent high-grade relapse. Multi-region sequencing and expression analysis from glioblastoma multiforme (GBM) patients have also shown that heterogeneity results in the presence of multiple tumor subtypes, as identified by gene expression classifiers, within the same tumor [37,38]. This questions the utility of such classification systems and gene expression signatures to define individual GBM subtypes.
The sequencing of two temporally separated CLL samples from treated and untreated patients showed that many of these tumors underwent clonal evolution and that the presence of a subclonal driver was an independent risk factor for disease progression . Conversely, in myelodysplastic syndrome, driver mutations had a similar prognostic significance whether they were clonal or subclonal; the absolute number of driver mutations rather than clonality had the biggest implications for outcome . The impact of intratumor heterogeneity and the clonality of driver mutations on prognosis and response to precision medicine has not been studied prospectively. A UK-based longitudinal observational study, Tracking Non-small Cell Lung Cancer Evolution Through Therapy (TRACERx), has been launched to assess this ; 842 patients will have whole-exome sequencing (WES) of multiple regions of their resected primary tumor, as well as of cell-free DNA (cfDNA) and circulating tumor cells (CTCs) at multiple time points throughout follow-up. Patients who suffer disease recurrence will be consented for repeat tumor sampling in order to define the evolutionary routes of individual tumors. This national study may provide insight into the role of intratumor heterogeneity and the (sub)clonality of driver events on outcome, reveal the origins of the lethal subclone and begin to define selection pressures initiated by therapy.
A number of studies have demonstrated that treatment can act as a selection pressure in some malignancies, driving clonal evolution and selecting for certain subclones. This emphasizes the importance of longitudinal tumor sampling strategies to depict tumor genomic landscapes. In acute lymphoblastic leukemia, copy number abnormalities (CNA) were strikingly different between samples taken at diagnosis and after relapse following chemotherapy . On retrospective analysis, the cells responsible for relapse were present as a subclonal population at diagnosis; chemotherapy had selected for a population with CNA in genes involved in the regulation of the cell cycle and B-cell development. In CLL, cancers were more likely to have undergone clonal evolution in patients treated with chemotherapy than in untreated patients . In non-small cell lung cancer (NSCLC), patients who relapsed with MET-amplified epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor (TKI)-resistant disease following treatment with an EGFR inhibitor harbored a low-frequency subclone (<1% of cells) with MET amplification prior to treatment, which was selected for during therapy . Notably, a subclonal population with MET amplification at such a low frequency would be difficult to identify in a heterogeneous biopsy sample. Conceivably, these patients may have benefited from combination EGFR TKI and MET inhibition to forestall selection of the drug-resistant subclone.
Similarly, the presence of the EGFR T790M mutation, which is associated with resistance to EGFR TKI therapy, has been demonstrated prior to treatment with EGFR inhibition in patients with NSCLC . In this study, matrix-assisted laser desorption/ionization mass spectrometry and next-generation sequencing (NGS) were used to detect the presence of low-frequency (<5% of cells) T790M mutations in pre-treatment samples that were not detected using standard Sanger sequencing. They found that a greater prevalence of the T790M subclonal population was detected in post-treatment biopsies and, in the context of selection due to EGFR TKI, the T790M mutation acted as a driver of subclonal expansion.
In a study of patients with colorectal cancer, multiple somatic mutations in KRAS (Kirsten rat sarcoma viral oncogene homolog), which are associated with resistance to anti-EGFR antibody therapy, could be detected non-invasively through cfDNA analysis during the acquisition of drug resistance, and were predicted to be present in a subclone prior to treatment . By contrast, a recent sequencing study of five patients with RAF/MEK-inhibitor-resistant BRAF melanoma found no evidence of a pre-existing resistant subclonal population, suggesting that resistant tumor cells had developed de novo on treatment, or that resistant subclones were present in the pre-treatment tumor at frequencies that were below the limits of detection .
Cancer therapies not only can act as the selection pressure to drive tumor evolution along a certain lineage if pre-existing subclones possess genotypes that are associated with a drug-resistant phenotype [18,19] but also can generate new subclonal driver events. For example, temozolamide, the standard first-line therapy for GBM, induces mutations in tumor DNA. Some are deleterious for the cell and result in death, others neutral and act as passenger mutations, but others such as mutations in mismatch repair (MMR) genes are potentially advantageous for tumor cells. Some GBM tumors treated with temozolamide exhibited a mutator phenotype, resulting from mutations in MMR genes, and were found to harbor driver mutations in RB1 (encoding retinoblastoma 1), PIK3CA and PTEN that bore the signature of temozolamide-induced mutagenesis [21,45].
Mechanisms of resistance to common cytotoxic chemotherapies and precision medicines
Decreased mismatch repair proficiency (e.g.↓MLH1 & ↓MSH2)
Increased efficiency of other modes of DNA repair (e.g. nucleotide excision repair – ↑ERCC1, trans lesion synthesis – ↑POLH, homologous recombination – BRCA1/2 restoration)
Drug efflux via increased expression of MDR-1
Changes in microtubule structure (e.g. mutations in β-tubulin and overexpression of βIII-tubulin
EGFR TK domain
Resistance mutation, e.g. T790M
Transformation to small cell lung cancer
EGFR monoclonal antibody
EGFR extracellular domain
Acquired KRAS or NRAS mutation
Activation of PIK3CA/PTEN pathway
Inhibition of cetuximab binding, e.g. EGFR-S492R
Elevated BRAF/CRAF/COT1 expression
Acquired mutation in other elements of the MAPK pathway
Persistent activation of other kinases, e.g. EGFR and PDGFRβ
Secondary EML4-ALK mutations or rearrangements
In some cases, it is becoming clear that resistance to therapy can be mediated by more than one resistant subclone. In colorectal cancer, different mutations in KRAS (exon 12 and 13), associated with resistance to EGFR monoclonal antibody therapy, were found in the same patient . In anaplastic lymphoma kinase (ALK)-positive, crizotinib-treated NSCLC patients, resistance was driven by multiple mechanisms, such as secondary ALK mutations, amplification of KIT or the EML4-ALK fusion gene itself and EGFR pathway activation, with multiple resistance mechanisms sometimes found simultaneously within the same tumor . Paradoxically, precision medicines may have a detrimental effect in the presence of polyclonal disease with subclonal driver events. BRAF inhibition has been shown to have significant anti-tumor efficacy in BRAFV600E mutant cancers. In BRAF wild-type cells, however, activation of extracellular-signal-regulated kinase signaling with the use of BRAF inhibitors can promote tumorigenesis . In the context of a polyclonal tumor where a BRAF-activating mutation is subclonal in nature with the presence of BRAF wild-type subclones, BRAF inhibition might promote the growth of the BRAF wild-type population, particularly if these subclones harbor KRAS or NRAS (neuroblastoma RAS viral (v-ras) oncogene homolog) driver mutations . The potential emergence of polyclonal resistance and the presence of subclonal drivers should be considered when designing clinical trials to forestall resistance to targeted agents.
The extent to which the presence of subclonal drivers and intratumor heterogeneity impacts upon the different responses witnessed with both systemic chemotherapy and precision medicines remains unclear [23,54–57]. Mixed responses to drug treatments are common and there is no consensus in clinical decision-making in this context . Should the treating physician switch systemic therapy, add a second systemic therapy to combat resistance, advocate locoregional approaches with surgery or radiotherapy, or monitor clinically insignificant progression and continue therapy with the caveat that non-responding lesions might act as an evolutionary sink that later contribute to more widespread progression? Further cataloguing of clonal and subclonal drivers and common mechanisms for resistance to treatment, together with adapting clinical trial design to the challenges of tumor evolution may improve patient care in the future.
The identification prior to therapy and subsequent management of low-frequency subclones harboring driver events that influence clinical outcome is clearly a major challenge. Such subclones may be present at a low frequency within one biopsy, could be spatially separated within a primary tumor or might be differently distributed in the primary tumor and its metastases. In recent years, clinical trial designs have moved from stratification that is based on histology to classification by molecular subtype [59,60]. The next generation of clinical trial design has seen the development of basket trials (MATCH, I-SPY, FOCUS4 and MATRIX) that contain multiple molecular subgroups, each treated with a different therapy, that are based on pre-specified somatic aberrations. These trials often have an adaptive design that allows removal of poorly performing arms during the course of the trial . A potential limitation of this approach is that, in the metastatic setting, molecular biomarker analysis for patient stratification is often performed on archival material, which may not reflect the current state of disease. This strategy may readily identify clonally dominant drivers that occur earlier in tumor evolution, but its ability to detect subclonal drivers will depend on the prevalence of that subclone within the tissue taken for analysis.
Without the use of single-cell approaches, addressing whether subclonal populations that have resistance to therapy are present prior to treatment or are selected during therapy is hindered by the limit of detection of low-frequency cancer cell populations. Conventional Sanger sequencing can detect variants down to the level of approximately 10% prevalence. Some studies have used mutation-specific PCR strategies followed by targeted sequencing or mass spectrometry to enrich low-frequency populations within archival material, allowing detection of subclonal events at 0.1% prevalence [62,63]. Massive parallel sequencing offers the potential to identify unknown subclonal populations. Nevertheless, the optimal detection of low-frequency subclones might in some cases require multi-region analysis, or analysis of circulating tumor DNA, to identify spatially separated populations; alternatively, high sequencing depth might be used to identify low frequency but potentially clinically significant events within single samples. In addition, the error rates of NGS, resulting from PCR errors during library preparation, instrument errors during cluster amplification and image analysis, or errors from base-calling algorithms, introduce a significant false-positive rate. Using molecular barcodes to build read families from sequencing data will help to reduce the errors generated during sequencing but not PCR errors introduced during library preparation [64,65]. The challenge of detecting significant subclonal driver mutations within the statistical noise generated by high-coverage NGS will need to be addressed to allow clinical application.
The advent of NGS has permitted greater resolution in determining the extent of heterogeneity within individual tumors. Distinguishing early and late somatic events may allow a better understanding of the genes and mutational processes that are involved in tumor initiation in comparison with those involved in maintenance and metastasis, which in turn might inform new therapeutic avenues. Moreover, the distinction of clonally dominant events from subclonal driver events might allow the acceleration of drug development towards targeted early, truncal drivers of disease. Although targeting clonally dominant drivers may make intuitive sense, some trunk drivers may only be relevant for tumor initiation and targeting these drivers after clonal diversification might not be efficacious.Furthermore, the origins of the ‘lethal’ subclone or subclones might be determined by somatic events that occur later in tumor evolution and not in the dominant clone. Therefore, efforts to limit disease progression might require a greater understanding and optimal targeting of subclonal driver events.
High-throughput functional assessment of validated mutations will be important to assess the significance of these mutations and to avoid wasting resources on the further investigation of sequencing artifact or passenger mutations. The wealth and depth of data may allow us to find unanticipated gene-gene interactions and might reveal new unappreciated cancer drivers that are involved in epistatic interactions . Clinical translation of these findings will be key. Novel clinical trial approaches that consider clonal evolution in the context of cancer diversity may shed light on the efficacy of unexpected combination therapies and provide evidence for adaptive therapy to avoid the selection of drug-resistant subclones [39,67–69].
Anaplastic lymphoma kinase
Acute myeloid leukemia
Apolipoprotein B mRNA-editing, enzyme-catalytic, polypeptide-like 3G
BRCA1 associated protein-1
B-Raf proto-oncogene, serine/threonine kinase
clear cell renal cell carcinoma
Chronic lymphocytic leukaemia
Copy number abnormalities
Circulating tumour cell
Deciphering Anti-tumour response and evolution with INtratumour heterogeneity
Epidermal growth factor receptor
v-erb-b2 avian erythroblastic leukemia viral oncogene homolog 2
Kirsten rat sarcoma viral oncogene homolog
Non-small cell lung cancer
Polymerase chain reaction
Phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha
Phosphatase and tensin homolog
Encoding retinoblastoma 1
SET domain containing 2
Tyrosine kinase inhibitor
Encoding tumor protein p53
Tracking non-small cell lung cancer evolution through therapy
ECdB is a Rosetrees Trust fellow, CH is funded by the National Institute for Health Research Academic Clinical Fellowship and CS is funded by Cancer Research UK, the Rosetrees Trust, EU FP7 (projects PREDICT and RESPONSIFY, ID:259303), the Prostate Cancer Foundation, the European Research Council, the Breast Cancer Research Foundation, and is supported by the National Institute for Health Research University College London Hospital Biomedical Research Centre.
The licensee has exclusive rights to distribute this article, in any medium, for 12 months following its publication. After this time, the article is available under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.