1One of the pervasive features of cancer is the diversity of mutations found in malignant 2 cells within the same tumor; a phenomenon called clonal diversity or intratumor heterogeneity. 3Clonal diversity allows tumors to adapt to the selective pressure of treatment and likely 4 contributes to the development of treatment resistance and cancer recurrence. Thus, the ability to 5 precisely delineate the clonal substructure of a tumor, including the evolutionary history of its 6 development and the co-occurrence of its mutations, is necessary to understand and overcome 7 treatment resistance. However, DNA sequencing of bulk tumor samples cannot accurately 8 resolve complex clonal architectures. Here, we performed high-throughput single-cell DNA 9 sequencing to quantitatively assess the clonal architecture of acute myeloid leukemia (AML). 10We sequenced a total of 556,951 cells from 77 patients with AML for 19 genes known to be 11 recurrently mutated in AML. The data revealed clonal relationship among AML driver mutations 12 and identified mutations that often co-occurred (e.g., NPM1/FLT3-ITD, DNMT3A/NPM1, 13 SRSF2/IDH2, and WT1/FLT3-ITD) and those that were mutually exclusive (e.g., NRAS/KRAS, 14 FLT3-D835/ITD, and IDH1/IDH2) at single-cell resolution. Reconstruction of the tumor 15 phylogeny uncovered history of tumor development that is characterized by linear and branching 16 clonal evolution patterns with latter involving functional convergence of separately evolved 17 clones. Analysis of longitudinal samples revealed remodeling of clonal architecture in response 18to therapeutic pressure that is driven by clonal selection. Furthermore, in this AML cohort, 19 higher clonal diversity (≥4 subclones) was associated with significantly worse overall survival. 20These data portray clonal relationship, architecture, and evolution of AML driver genes with 21 unprecedented resolution, and illuminate the role of clonal diversity in therapeutic resistance, 22relapse and clinical outcome in AML. 23
Main 1A growing body of evidence supports the role of clonal diversity in therapeutic 2 resistance, recurrence, and poor outcomes in cancer 1 . Clonal diversity also reflects the history of 3 the accumulation of somatic mutations within a tumor. Thus, a precise characterization of clonal 4 diversity reveals not only the extent of a tumor's clonal complexity but also the evolutionary 5 history of the tumor's development. Much of the work characterizing the clonal architecture of 6 tumors has been done by computational inference using variant allele fraction (VAF) data from 7 massively parallel DNA sequencing of bulk tumor samples 2,3 . However, the ability to infer 8 clonal heterogeneity and tumor phylogeny from bulk sequencing data is inherently limited, 9 because bulk sequencing techniques cannot reliably infer mutation co-occurrences and hence 10