Pancreatic cancer is a highly lethal malignancy with few effective therapies. We performed exome sequencing and copy number analysis to define genomic aberrations in a prospectively accrued clinical cohort (n = 142) of early (stage I and II) sporadic pancreatic ductal adenocarcinoma. Detailed analysis of 99 informative tumours identified substantial heterogeneity with 2,016 non-silent mutations and 1,628 copy-number variations. We define 16 significantly mutated genes, reaffirming known mutations (KRAS, TP53, CDKN2A, SMAD4, MLL3, TGFBR2, ARID1A and SF3B1), and uncover novel mutated genes including additional genes involved in chromatin modification (EPC1 and ARID2), DNA damage repair (ATM) and other mechanisms (ZIM2, MAP2K4, NALCN, SLC16A4 and MAGEA6). Integrative analysis with in vitro functional data and animal models provided supportive evidence for potential roles for these genetic aberrations in carcinogenesis. Pathway-based analysis of recurrently mutated genes recapitulated clustering in core signalling pathways in pancreatic ductal adenocarcinoma, and identified new mutated genes in each pathway. We also identified frequent and diverse somatic aberrations in genes described traditionally as embryonic regulators of axon guidance, particularly SLIT/ROBO signalling, which was also evident in murine Sleeping Beauty transposon-mediated somatic mutagenesis models of pancreatic cancer, providing further supportive evidence for the potential involvement of axon guidance genes in pancreatic carcinogenesis.

Summary In acute myeloid leukemia (AML), the cell of origin, nature and biological consequences of initiating lesions and order of subsequent mutations remain poorly understood, as AML is typically diagnosed without observation of a pre-leukemic phase. Here, highly purified hematopoietic stem cells (HSC), progenitor and mature cell fractions from the blood of AML patients were found to contain recurrent DNMT3a mutations (DNMT3amut) at high allele frequency, but without coincident NPM1 mutations (NPM1c) present in AML blasts. DNMT3amut-bearing HSC exhibited multilineage repopulation advantage over non-mutated HSC in xenografts, establishing their identity as pre-leukemic-HSC (preL-HSC). preL-HSC were found in remission samples indicating that they survive chemotherapy. Thus DNMT3amut arises early in AML evolution, likely in HSC, leading to a clonally expanded pool of preL-HSC from which AML evolves. Our findings provide a paradigm for the detection and treatment of pre-leukemic clones before the acquisition of additional genetic lesions engenders greater therapeutic resistance.

Intratumoral heterogeneity arises through the evolution of genetically diverse subclones during tumor progression. However, it remains unknown whether cells within single genetic clones are functionally equivalent. By combining DNA copy number alteration (CNA) profiling, sequencing, and lentiviral lineage tracking, we followed the repopulation dynamics of 150 single lentivirus-marked lineages from 10 human colorectal cancers through serial xenograft passages in mice. CNA and mutational analysis distinguished individual clones and showed that clones remained stable upon serial transplantation. Despite this stability, the proliferation, persistence, and chemotherapy tolerance of lentivirally marked lineages were variable within each clone. Chemotherapy promoted the dominance of previously minor or dormant lineages. Thus, apart from genetic diversity, tumor cells display inherent functional variability in tumor propagation potential, which contributes to both cancer growth and therapy tolerance.

Pancreas cancer (PC), a highly aggressive tumour type with uniformly poor prognosis, is an exemplar of the classical view of cancer development based on stepwise progression1. The current progression model, based on analyses of precursor lesions termed pancreatic intraepithelial neoplasm (PanINs) lesions, makes two predictions: 1) PC develops through a particular sequence of genetic alterations2–5 (KRAS > CDKN2A > TP53/SMAD4); and 2) the evolutionary trajectory of PC progression is gradual because each alteration is acquired independently. One shortcoming of this nearly two decade old contention is that clonally expanded precursor lesions have been identified that do not always belong to the tumour lineage2,5–9, indicating that the evolutionary trajectory of the tumour lineage and precursor lesions can be divergent. This prevailing view of tumourigenesis has contributed to the clinical notion that PC evolves slowly and presents at a late stage10. However, the propensity for this disease to rapidly metastasize and the inability to improve patient outcomes despite efforts aimed at early detection11, argue that PC progression is anything but gradual. By tracking DNA copy number changes and their associated rearrangements from tumour-enriched genomes using novel informatics tools, we found that PC tumourigenesis neither is gradual nor follows the accepted mutation order. Two-thirds of tumours harbour complex rearrangement patterns associated with mitotic errors, consistent with punctuated equilibrium as the principal evolutionary trajectory12. In a subset of cases, the consequence of such errors was the simultaneous, rather than sequential, knockout of canonical preneoplastic genetic drivers that likely set-off invasive cancer growth. These findings challenge the current model of PC tumourigenesis and provide novel insights into the mutational processes giving rise to these aggressive tumours.

Herein we provide a detailed molecular analysis of the spatial heterogeneity of clinically localized, multifocal prostate cancer to delineate new oncogenes or tumor suppressors. We initially determined the copy number aberration (CNA) profiles of 74 patients with index tumors of Gleason score 7. Of these, 5 patients were subjected to whole-genome sequencing using DNA quantities achievable in diagnostic biopsies, with detailed spatial sampling of 23 distinct tumor regions to assess intraprostatic heterogeneity in focal genomics. Multifocal tumors are highly heterogeneous for single-nucleotide variants (SNVs), CNAs and genomic rearrangements. We identified and validated a new recurrent amplification of MYCL, which is associated with TP53 deletion and unique profiles of DNA damage and transcriptional dysregulation. Moreover, we demonstrate divergent tumor evolution in multifocal cancer and, in some cases, tumors of independent clonal origin. These data represent the first systematic relation of intraprostatic genomic heterogeneity to predicted clinical outcome and inform the development of novel biomarkers that reflect individual prognosis.