Validation of ITD mutations in FLT3 as a therapeutic target in human acute myeloid leukaemia

Effective targeted cancer therapeutic development depends upon distinguishing disease-associated 'driver' mutations, which have causative roles in malignancy pathogenesis, from 'passenger' mutations, which are dispensable for cancer initiation and maintenance. Translational studies of clinically active targeted therapeutics can definitively discriminate driver from passenger lesions and provide valuable insights into human cancer biology. Activating internal tandem duplication (ITD) mutations in FLT3 (FLT3-ITD) are detected in approximately 20% of acute myeloid leukaemia (AML) patients and are associated with a poor prognosis. Abundant scientific and clinical evidence, including the lack of convincing clinical activity of early FLT3 inhibitors, suggests that FLT3-ITD probably represents a passenger lesion. Here we report point mutations at three residues within the kinase domain of FLT3-ITD that confer substantial in vitro resistance to AC220 (quizartinib), an active investigational inhibitor of FLT3, KIT, PDGFRA, PDGFRB and RET; evolution of AC220-resistant substitutions at two of these amino acid positions was observed in eight of eight FLT3-ITD-positive AML patients with acquired resistance to AC220. Our findings demonstrate that FLT3-ITD can represent a driver lesion and valid therapeutic target in human AML. AC220-resistant FLT3 kinase domain mutants represent high-value targets for future FLT3 inhibitor development efforts.     

The clonal and mutational evolution spectrum of primary triple-negative breast cancers

Primary triple-negative breast cancers (TNBCs), a tumour type defined by lack of oestrogen receptor, progesterone receptor and ERBB2 gene amplification, represent approximately 16% of all breast cancers. Here we show in 104 TNBC cases that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some having only a handful of coding somatic aberrations in a few pathways, whereas others contain hundreds of coding somatic mutations. High-throughput RNA sequencing (RNA-seq) revealed that only approximately 36% of mutations are expressed. Using deep re-sequencing measurements of allelic abundance for 2,414 somatic mutations, we determine for the first time-to our knowledge-in an epithelial tumour subtype, the relative abundance of clonal frequencies among cases representative of the population. We show that TNBCs vary widely in their clonal frequencies at the time of diagnosis, with the basal subtype of TNBC showing more variation than non-basal TNBC. Although p53 (also known as TP53), PIK3CA and PTEN somatic mutations seem to be clonally dominant compared to other genes, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal, cell shape and motility proteins occurred at lower clonal frequencies, suggesting that they occurred later during tumour progression. Taken together, our results show that understanding the biology and therapeutic responses of patients with TNBC will require the determination of individual tumour clonal genotypes.     

Clonal expansion of p53 mutant cells is associated with brain tumour progression.

Tumour progression is a fundamental feature of the biology of cancer. Cancers do not arise de novo in their final form, but begin as small, indolent growths, which gradually acquire characteristics associated with malignancy. In the brain, for example, low-grade tumours (astrocytomas) evolve into faster growing, more dysplastic and invasive high-grade tumours (glioblastomas). To define the genetic events underlying brain tumour progression, we analysed the p53 gene in ten primary brain tumour pairs. Seven pairs consisted of tumours that were high grade both at presentation and recurrence (group A) and three pairs consisted of low-grade tumours that had progressed to higher grade tumours (group B). In group A pairs, four of the recurrent tumours contained a p53 gene mutation; in three of them, the same mutation was found in the primary tumour. In group B pairs, progression to high grade was associated with a p53 gene mutation. A subpopulation of cells were present in the low-grade tumours that contained the same p53 gene mutation predominant in the cells of the recurrent tumours that had progressed to glioblastoma. Thus, the histological progression of brain tumours was associated with a clonal expansion of cells that had previously acquired a mutation in the p53 gene, endowing them with a selective growth advantage. These experimental observations strongly support Nowell's clonal evolution model of tumour progression.     

Distant metastasis occurs late during the genetic evolution of pancreatic cancer

Metastasis, the dissemination and growth of neoplastic cells in an organ distinct from that in which they originated, is the most common cause of death in cancer patients. This is particularly true for pancreatic cancers, where most patients are diagnosed with metastatic disease and few show a sustained response to chemotherapy or radiation therapy. Whether the dismal prognosis of patients with pancreatic cancer compared to patients with other types of cancer is a result of late diagnosis or early dissemination of disease to distant organs is not known. Here we rely on data generated by sequencing the genomes of seven pancreatic cancer metastases to evaluate the clonal relationships among primary and metastatic cancers. We find that clonal populations that give rise to distant metastases are represented within the primary carcinoma, but these clones are genetically evolved from the original parental, non-metastatic clone. Thus, genetic heterogeneity of metastases reflects that within the primary carcinoma. A quantitative analysis of the timing of the genetic evolution of pancreatic cancer was performed, indicating at least a decade between the occurrence of the initiating mutation and the birth of the parental, non-metastatic founder cell. At least five more years are required for the acquisition of metastatic ability and patients die an average of two years thereafter. These data provide novel insights into the genetic features underlying pancreatic cancer progression and define a broad time window of opportunity for early detection to prevent deaths from metastatic disease.     

The patterns and dynamics of genomic instability in metastatic pancreatic cancer

Pancreatic cancer is an aggressive malignancy with a five-year mortality of 97-98%, usually due to widespread metastatic disease. Previous studies indicate that this disease has a complex genomic landscape, with frequent copy number changes and point mutations, but genomic rearrangements have not been characterized in detail. Despite the clinical importance of metastasis, there remain fundamental questions about the clonal structures of metastatic tumours, including phylogenetic relationships among metastases, the scale of ongoing parallel evolution in metastatic and primary sites, and how the tumour disseminates. Here we harness advances in DNA sequencing to annotate genomic rearrangements in 13 patients with pancreatic cancer and explore clonal relationships among metastases. We find that pancreatic cancer acquires rearrangements indicative of telomere dysfunction and abnormal cell-cycle control, namely dysregulated G1-to-S-phase transition with intact G2-M checkpoint. These initiate amplification of cancer genes and occur predominantly in early cancer development rather than the later stages of the disease. Genomic instability frequently persists after cancer dissemination, resulting in ongoing, parallel and even convergent evolution among different metastases. We find evidence that there is genetic heterogeneity among metastasis-initiating cells, that seeding metastasis may require driver mutations beyond those required for primary tumours, and that phylogenetic trees across metastases show organ-specific branches. These data attest to the richness of genetic variation in cancer, brought about by the tandem forces of genomic instability and evolutionary selection.     

Tumour evolution inferred by single-cell sequencing

Genomic analysis provides insights into the role of copy number variation in disease, but most methods are not designed to resolve mixed populations of cells. In tumours, where genetic heterogeneity is common, very important information may be lost that would be useful for reconstructing evolutionary history. Here we show that with flow-sorted nuclei, whole genome amplification and next generation sequencing we can accurately quantify genomic copy number within an individual nucleus. We apply single-nucleus sequencing to investigate tumour population structure and evolution in two human breast cancer cases. Analysis of 100 single cells from a polygenomic tumour revealed three distinct clonal subpopulations that probably represent sequential clonal expansions. Additional analysis of 100 single cells from a monogenomic primary tumour and its liver metastasis indicated that a single clonal expansion formed the primary tumour and seeded the metastasis. In both primary tumours, we also identified an unexpectedly abundant subpopulation of genetically diverse 'pseudodiploid' cells that do not travel to the metastatic site. In contrast to gradual models of tumour progression, our data indicate that tumours grow by punctuated clonal expansions with few persistent intermediates.     

DNA sequencing of a cytogenetically normal acute myeloid leukaemia genome

Acute myeloid leukaemia is a highly malignant haematopoietic tumour that affects about 13,000 adults in the United States each year. The treatment of this disease has changed little in the past two decades, because most of the genetic events that initiate the disease remain undiscovered. Whole-genome sequencing is now possible at a reasonable cost and timeframe to use this approach for the unbiased discovery of tumour-specific somatic mutations that alter the protein-coding genes. Here we present the results obtained from sequencing a typical acute myeloid leukaemia genome, and its matched normal counterpart obtained from the same patient's skin. We discovered ten genes with acquired mutations; two were previously described mutations that are thought to contribute to tumour progression, and eight were new mutations present in virtually all tumour cells at presentation and relapse, the function of which is not yet known. Our study establishes whole-genome sequencing as an unbiased method for discovering cancer-initiating mutations in previously unidentified genes that may respond to targeted therapies.     

Genetic variegation of clonal architecture and propagating cells in leukaemia

Little is known of the genetic architecture of cancer at the subclonal and single-cell level or in the cells responsible for cancer clone maintenance and propagation. Here we have examined this issue in childhood acute lymphoblastic leukaemia in which the ETV6-RUNX1 gene fusion is an early or initiating genetic lesion followed by a modest number of recurrent or 'driver' copy number alterations. By multiplexing fluorescence in situ hybridization probes for these mutations, up to eight genetic abnormalities can be detected in single cells, a genetic signature of subclones identified and a composite picture of subclonal architecture and putative ancestral trees assembled. Subclones in acute lymphoblastic leukaemia have variegated genetics and complex, nonlinear or branching evolutionary histories. Copy number alterations are independently and reiteratively acquired in subclones of individual patients, and in no preferential order. Clonal architecture is dynamic and is subject to change in the lead-up to a diagnosis and in relapse. Leukaemia propagating cells, assayed by serial transplantation in NOD/SCID IL2Rγ(null) mice, are also genetically variegated, mirroring subclonal patterns, and vary in competitive regenerative capacity in vivo. These data have implications for cancer genomics and for the targeted therapy of cancer.     

Activating mutations in ALK provide a therapeutic target in neuroblastoma

Neuroblastoma, an embryonal tumour of the peripheral sympathetic nervous system, accounts for approximately 15% of all deaths due to childhood cancer. High-risk neuroblastomas are rapidly progressive; even with intensive myeloablative chemotherapy, relapse is common and almost uniformly fatal. Here we report the detection of previously unknown mutations in the ALK gene, which encodes a receptor tyrosine kinase, in 8% of primary neuroblastomas. Five non-synonymous sequence variations were identified in the kinase domain of ALK, of which three were somatic and two were germ line. The most frequent mutation, F1174L, was also identified in three different neuroblastoma cell lines. ALK complementary DNAs encoding the F1174L and R1275Q variants, but not the wild-type ALK cDNA, transformed interleukin-3-dependent murine haematopoietic Ba/F3 cells to cytokine-independent growth. Ba/F3 cells expressing these mutations were sensitive to the small-molecule inhibitor of ALK, TAE684 (ref. 4). Furthermore, two human neuroblastoma cell lines harbouring the F1174L mutation were also sensitive to the inhibitor. Cytotoxicity was associated with increased amounts of apoptosis as measured by TdT-mediated dUTP nick end labelling (TUNEL). Short hairpin RNA (shRNA)-mediated knockdown of ALK expression in neuroblastoma cell lines with the F1174L mutation also resulted in apoptosis and impaired cell proliferation. Thus, activating alleles of the ALK receptor tyrosine kinase are present in primary neuroblastoma tumours and in established neuroblastoma cell lines, and confer sensitivity to ALK inhibition with small molecules, providing a molecular rationale for targeted therapy of this disease.     

Secondary mutations as a mechanism of cisplatin resistance in BRCA2-mutated cancers

Ovarian carcinomas with mutations in the tumour suppressor BRCA2 are particularly sensitive to platinum compounds. However, such carcinomas ultimately develop cisplatin resistance. The mechanism of that resistance is largely unknown. Here we show that acquired resistance to cisplatin can be mediated by secondary intragenic mutations in BRCA2 that restore the wild-type BRCA2 reading frame. First, in a cisplatin-resistant BRCA2-mutated breast-cancer cell line, HCC1428, a secondary genetic change in BRCA2 rescued BRCA2 function. Second, cisplatin selection of a BRCA2-mutated pancreatic cancer cell line, Capan-1 (refs 3, 4), led to five different secondary mutations that restored the wild-type BRCA2 reading frame. All clones with secondary mutations were resistant both to cisplatin and to a poly(ADP-ribose) polymerase (PARP) inhibitor (AG14361). Finally, we evaluated recurrent cancers from patients whose primary BRCA2-mutated ovarian carcinomas were treated with cisplatin. The recurrent tumour that acquired cisplatin resistance had undergone reversion of its BRCA2 mutation. Our results suggest that secondary mutations that restore the wild-type BRCA2 reading frame may be a major clinical mediator of acquired resistance to platinum-based chemotherapy.     

Whole-genome analysis informs breast cancer response to aromatase inhibition

To correlate the variable clinical features of oestrogen-receptor-positive breast cancer with somatic alterations, we studied pretreatment tumour biopsies accrued from patients in two studies of neoadjuvant aromatase inhibitor therapy by massively parallel sequencing and analysis. Eighteen significantly mutated genes were identified, including five genes (RUNX1, CBFB, MYH9, MLL3 and SF3B1) previously linked to haematopoietic disorders. Mutant MAP3K1 was associated with luminal A status, low-grade histology and low proliferation rates, whereas mutant TP53 was associated with the opposite pattern. Moreover, mutant GATA3 correlated with suppression of proliferation upon aromatase inhibitor treatment. Pathway analysis demonstrated that mutations in MAP2K4, a MAP3K1 substrate, produced similar perturbations as MAP3K1 loss. Distinct phenotypes in oestrogen-receptor-positive breast cancer are associated with specific patterns of somatic mutations that map into cellular pathways linked to tumour biology, but most recurrent mutations are relatively infrequent. Prospective clinical trials based on these findings will require comprehensive genome sequencing.     

A murine lung cancer co-clinical trial identifies genetic modifiers of therapeutic response

Targeted therapies have demonstrated efficacy against specific subsets of molecularly defined cancers. Although most patients with lung cancer are stratified according to a single oncogenic driver, cancers harbouring identical activating genetic mutations show large variations in their responses to the same targeted therapy. The biology underlying this heterogeneity is not well understood, and the impact of co-existing genetic mutations, especially the loss of tumour suppressors, has not been fully explored. Here we use genetically engineered mouse models to conduct a 'co-clinical' trial that mirrors an ongoing human clinical trial in patients with KRAS-mutant lung cancers. This trial aims to determine if the MEK inhibitor selumetinib (AZD6244) increases the efficacy of docetaxel, a standard of care chemotherapy. Our studies demonstrate that concomitant loss of either p53 (also known as Tp53) or Lkb1 (also known as Stk11), two clinically relevant tumour suppressors, markedly impaired the response of Kras-mutant cancers to docetaxel monotherapy. We observed that the addition of selumetinib provided substantial benefit for mice with lung cancer caused by Kras and Kras and p53 mutations, but mice with Kras and Lkb1 mutations had primary resistance to this combination therapy. Pharmacodynamic studies, including positron-emission tomography (PET) and computed tomography (CT), identified biological markers in mice and patients that provide a rationale for the differential efficacy of these therapies in the different genotypes. These co-clinical results identify predictive genetic biomarkers that should be validated by interrogating samples from patients enrolled on the concurrent clinical trial. These studies also highlight the rationale for synchronous co-clinical trials, not only to anticipate the results of ongoing human clinical trials, but also to generate clinically relevant hypotheses that can inform the analysis and design of human studies.     

Impaired hydroxylation of 5-methylcytosine in myeloid cancers with mutant TET2

TET2 is a close relative of TET1, an enzyme that converts 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) in DNA. The gene encoding TET2 resides at chromosome 4q24, in a region showing recurrent microdeletions and copy-neutral loss of heterozygosity (CN-LOH) in patients with diverse myeloid malignancies. Somatic TET2 mutations are frequently observed in myelodysplastic syndromes (MDS), myeloproliferative neoplasms (MPN), MDS/MPN overlap syndromes including chronic myelomonocytic leukaemia (CMML), acute myeloid leukaemias (AML) and secondary AML (sAML). We show here that TET2 mutations associated with myeloid malignancies compromise catalytic activity. Bone marrow samples from patients with TET2 mutations displayed uniformly low levels of 5hmC in genomic DNA compared to bone marrow samples from healthy controls. Moreover, small hairpin RNA (shRNA)-mediated depletion of Tet2 in mouse haematopoietic precursors skewed their differentiation towards monocyte/macrophage lineages in culture. There was no significant difference in DNA methylation between bone marrow samples from patients with high 5hmC versus healthy controls, but samples from patients with low 5hmC showed hypomethylation relative to controls at the majority of differentially methylated CpG sites. Our results demonstrate that Tet2 is important for normal myelopoiesis, and suggest that disruption of TET2 enzymatic activity favours myeloid tumorigenesis. Measurement of 5hmC levels in myeloid malignancies may prove valuable as a diagnostic and prognostic tool, to tailor therapies and assess responses to anticancer drugs.     

Defining the mode of tumour growth by clonal analysis

Recent studies using the isolation of a subpopulation of tumour cells followed by their transplantation into immunodeficient mice provide evidence that certain tumours, including squamous skin tumours, contain cells with high clonogenic potential that have been referred to as cancer stem cells (CSCs). Until now, CSC properties have only been investigated by transplantation assays, and their existence in unperturbed tumour growth is unproven. Here we make use of clonal analysis of squamous skin tumours using genetic lineage tracing to unravel the mode of tumour growth in vivo in its native environment. To this end, we used a genetic labelling strategy that allows individual tumour cells to be marked and traced over time at different stages of tumour progression. Surprisingly, we found that the majority of labelled tumour cells in benign papilloma have only limited proliferative potential, whereas a fraction has the capacity to persist long term, giving rise to progeny that occupy a significant part of the tumour. As well as confirming the presence of two distinct proliferative cell compartments within the papilloma, mirroring the composition, hierarchy and fate behaviour of normal tissue, quantitative analysis of clonal fate data indicates that the more persistent population has stem-cell-like characteristics and cycles twice per day, whereas the second represents a slower cycling transient population that gives rise to terminally differentiated tumour cells. Such behaviour is shown to be consistent with double-labelling experiments and detailed clonal fate characteristics. By contrast, measurements of clone size and proliferative potential in invasive squamous cell carcinoma show a different pattern of behaviour, consistent with geometric expansion of a single CSC population with limited potential for terminal differentiation. This study presents the first experimental evidence for the existence of CSCs during unperturbed solid tumour growth.     

Whole-genome sequencing identifies recurrent mutations in chronic lymphocytic leukaemia

Chronic lymphocytic leukaemia (CLL), the most frequent leukaemia in adults in Western countries, is a heterogeneous disease with variable clinical presentation and evolution. Two major molecular subtypes can be distinguished, characterized respectively by a high or low number of somatic hypermutations in the variable region of immunoglobulin genes. The molecular changes leading to the pathogenesis of the disease are still poorly understood. Here we performed whole-genome sequencing of four cases of CLL and identified 46 somatic mutations that potentially affect gene function. Further analysis of these mutations in 363 patients with CLL identified four genes that are recurrently mutated: notch 1 (NOTCH1), exportin 1 (XPO1), myeloid differentiation primary response gene 88 (MYD88) and kelch-like 6 (KLHL6). Mutations in MYD88 and KLHL6 are predominant in cases of CLL with mutated immunoglobulin genes, whereas NOTCH1 and XPO1 mutations are mainly detected in patients with unmutated immunoglobulins. The patterns of somatic mutation, supported by functional and clinical analyses, strongly indicate that the recurrent NOTCH1, MYD88 and XPO1 mutations are oncogenic changes that contribute to the clinical evolution of the disease. To our knowledge, this is the first comprehensive analysis of CLL combining whole-genome sequencing with clinical characteristics and clinical outcomes. It highlights the usefulness of this approach for the identification of clinically relevant mutations in cancer.     

PML targeting eradicates quiescent leukaemia-initiating cells

The existence of a small population of 'cancer-initiating cells' responsible for tumour maintenance has been firmly demonstrated in leukaemia. This concept is currently being tested in solid tumours. Leukaemia-initiating cells, particularly those that are in a quiescent state, are thought to be resistant to chemotherapy and targeted therapies, resulting in disease relapse. Chronic myeloid leukaemia is a paradigmatic haematopoietic stem cell disease in which the leukaemia-initiating-cell pool is not eradicated by current therapy, leading to disease relapse on drug discontinuation. Here we define the critical role of the promyelocytic leukaemia protein (PML) tumour suppressor in haematopoietic stem cell maintenance, and present a new therapeutic approach for targeting quiescent leukaemia-initiating cells and possibly cancer-initiating cells by pharmacological inhibition of PML.     

Evolution of human BCR-ABL1 lymphoblastic leukaemia-initiating cells

Many tumours are composed of genetically diverse cells; however, little is known about how diversity evolves or the impact that diversity has on functional properties. Here, using xenografting and DNA copy number alteration (CNA) profiling of human BCR-ABL1 lymphoblastic leukaemia, we demonstrate that genetic diversity occurs in functionally defined leukaemia-initiating cells and that many diagnostic patient samples contain multiple genetically distinct leukaemia-initiating cell subclones. Reconstructing the subclonal genetic ancestry of several samples by CNA profiling demonstrated a branching multi-clonal evolution model of leukaemogenesis, rather than linear succession. For some patient samples, the predominant diagnostic clone repopulated xenografts, whereas in others it was outcompeted by minor subclones. Reconstitution with the predominant diagnosis clone was associated with more aggressive growth properties in xenografts, deletion of CDKN2A and CDKN2B, and a trend towards poorer patient outcome. Our findings link clonal diversity with leukaemia-initiating-cell function and underscore the importance of developing therapies that eradicate all intratumoral subclones.