Single-cell genomic variation induced by mutational processes in cancer.

How cell-to-cell copy number alterations that underpin genomic instability1 in human cancers drive genomic and phenotypic variation, and consequently the evolution of cancer2, remains understudied. Here, by applying scaled single-cell whole-genome sequencing3 to wild-type, TP53-deficient and TP53-deficient;BRCA1-deficient or TP53-deficient;BRCA2-deficient mammary epithelial cells (13,818 genomes), and to primary triple-negative breast cancer (TNBC) and high-grade serous ovarian cancer (HGSC) cells (22,057 genomes), we identify three distinct 'foreground' mutational patterns that are defined by cell-to-cell structural variation. Cell- and clone-specific high-level amplifications, parallel haplotype-specific copy number alterations and copy number segment length variation (serrate structural variations) had measurable phenotypic and evolutionary consequences. In TNBC and HGSC, clone-specific high-level amplifications in known oncogenes were highly prevalent in tumours bearing fold-back inversions, relative to tumours with homologous recombination deficiency, and were associated with increased clone-to-clone phenotypic variation. Parallel haplotype-specific alterations were also commonly observed, leading to phylogenetic evolutionary diversity and clone-specific mono-allelic expression. Serrate variants were increased in tumours with fold-back inversions and were highly correlated with increased genomic diversity of cellular populations. Together, our findings show that cell-to-cell structural variation contributes to the origins of phenotypic and evolutionary diversity in TNBC and HGSC, and provide insight into the genomic and mutational states of individual cancer cells.

Accurate determination of CRISPR-mediated gene fitness in transplantable tumours.

Assessing tumour gene fitness in physiologically-relevant model systems is challenging due to biological features of in vivo tumour regeneration, including extreme variations in single cell lineage progeny. Here we develop a reproducible, quantitative approach to pooled genetic perturbation in patient-derived xenografts (PDXs), by encoding single cell output from transplanted CRISPR-transduced cells in combination with a Bayesian hierarchical model. We apply this to 181 PDX transplants from 21 breast cancer patients. We show that uncertainty in fitness estimates depends critically on the number of transplant cell clones and the variability in clone sizes. We use a pathway-directed allelic series to characterize Notch signaling, and quantify TP53 / MDM2 drug-gene conditional fitness in outlier patients. We show that fitness outlier identification can be mirrored by pharmacological perturbation. Overall, we demonstrate that the gene fitness landscape in breast PDXs is dominated by inter-patient differences.

Results of the phase I CCTG IND.231 trial of CX-5461 in patients with advanced solid tumors enriched for DNA-repair deficiencies.

CX-5461 is a G-quadruplex stabilizer that exhibits synthetic lethality in homologous recombination-deficient models. In this multicentre phase I trial in patients with solid tumors, 40 patients are treated across 10 dose levels (50-650 mg/m2) to determine the recommended phase II dose (primary outcome), and evaluate safety, tolerability, pharmacokinetics (secondary outcomes). Defective homologous recombination is explored as a predictive biomarker of response. CX-5461 is generally well tolerated, with a recommended phase II dose of 475 mg/m2 days 1, 8 and 15 every 4 weeks, and dose limiting phototoxicity. Responses are observed in 14% of patients, primarily in patients with defective homologous recombination. Reversion mutations in PALB2 and BRCA2 are detected on progression following initial response in germline carriers, confirming the underlying synthetic lethal mechanism. In vitro characterization of UV sensitization shows this toxicity is related to the CX-5461 chemotype, independent of G-quadruplex synthetic lethality. These results establish clinical proof-of-concept for this G-quadruplex stabilizer. Clinicaltrials.gov NCT02719977.

Pharmacological systems analysis defines EIF4A3 functions in cell-cycle and RNA stress granule formation.

The RNA helicase EIF4A3 regulates the exon junction complex and nonsense-mediated mRNA decay functions in RNA transcript processing. However, a transcriptome-wide network definition of these functions has been lacking, in part due to the lack of suitable pharmacological inhibitors. Here we employ short-duration graded EIF4A3 inhibition using small molecule allosteric inhibitors to define the transcriptome-wide dependencies of EIF4A3. We thus define conserved cellular functions, such as cell cycle control, that are EIF4A3 dependent. We show that EIF4A3-dependent splicing reactions have a distinct genome-wide pattern of associated RNA-binding protein motifs. We also uncover an unanticipated role of EIF4A3 in the biology of RNA stress granules, which sequester and silence the translation of most mRNAs under stress conditions and are implicated in cell survival and tumour progression. We show that stress granule induction and maintenance is suppressed on the inhibition of EIF4A3, in part through EIF4A3-associated regulation of G3BP1 and TIA1 scaffold protein expression.

Engineered in-vitro cell line mixtures and robust evaluation of computational methods for clonal decomposition and longitudinal dynamics in cancer.

Characterization and quantification of tumour clonal populations over time via longitudinal sampling are essential components in understanding and predicting the response to therapeutic interventions. Computational methods for inferring tumour clonal composition from deep-targeted sequencing data are ubiquitous, however due to the lack of a ground truth biological data, evaluating their performance is difficult. In this work, we generate a benchmark data set that simulates tumour longitudinal growth and heterogeneity by in vitro mixing of cancer cell lines with known proportions. We apply four different algorithms to our ground truth data set and assess their performance in inferring clonal composition using different metrics. We also analyse the performance of these algorithms on breast tumour xenograft samples. We conclude that methods that can simultaneously analyse multiple samples while accounting for copy number alterations as a factor in allelic measurements exhibit the most accurate predictions. These results will inform future functional genomics oriented studies of model systems where time series measurements in the context of therapeutic interventions are becoming increasingly common. These studies will need computational models which accurately reflect the multi-factorial nature of allele measurement in cancer including, as we show here, segmental aneuploidies.

The histone methyltransferase EZH2 is a therapeutic target in small cell carcinoma of the ovary, hypercalcaemic type.

Small cell carcinoma of the ovary, hypercalcaemic type (SCCOHT) is a rare but aggressive and untreatable malignancy affecting young women. We and others recently discovered that SMARCA4, a gene encoding the ATPase of the SWI/SNF chromatin-remodelling complex, is the only gene recurrently mutated in the majority of SCCOHT. The low somatic complexity of SCCOHT genomes and the prominent role of the SWI/SNF chromatin-remodelling complex in transcriptional control of genes suggest that SCCOHT cells may rely on epigenetic rewiring for oncogenic transformation. Herein, we report that approximately 80% (19/24) of SCCOHT tumour samples have strong expression of the histone methyltransferase EZH2 by immunohistochemistry, with the rest expressing variable amounts of EZH2. Re-expression of SMARCA4 suppressed the expression of EZH2 in SCCOHT cells. In comparison to other ovarian cell lines, SCCOHT cells displayed hypersensitivity to EZH2 shRNAs and two selective EZH2 inhibitors, GSK126 and EPZ-6438. EZH2 inhibitors induced cell cycle arrest, apoptosis, and cell differentiation in SCCOHT cells, along with the induction of genes involved in cell cycle regulation, apoptosis, and neuron-like differentiation. EZH2 inhibitors suppressed tumour growth and improved the survival of mice bearing SCCOHT xenografts. Therefore, our data suggest that loss of SMARCA4 creates a dependency on the catalytic activity of EZH2 in SCCOHT cells and that pharmacological inhibition of EZH2 is a promising therapeutic strategy for treating this disease. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Atrophin controls developmental signaling pathways via interactions with Trithorax-like.

Mutations in human Atrophin1, a transcriptional corepressor, cause dentatorubral-pallidoluysian atrophy, a neurodegenerative disease. Drosophila Atrophin (Atro) mutants display many phenotypes, including neurodegeneration, segmentation, patterning and planar polarity defects. Despite Atro's critical role in development and disease, relatively little is known about Atro's binding partners and downstream targets. We present the first genomic analysis of Atro using ChIP-seq against endogenous Atro. ChIP-seq identified 1300 potential direct targets of Atro including engrailed, and components of the Dpp and Notch signaling pathways. We show that Atro regulates Dpp and Notch signaling in larval imaginal discs, at least partially via regulation of thickveins and fringe. In addition, bioinformatics analyses, sequential ChIP and coimmunoprecipitation experiments reveal that Atro interacts with the Drosophila GAGA Factor, Trithorax-like (Trl), and they bind to the same loci simultaneously. Phenotypic analyses of Trl and Atro clones suggest that Atro is required to modulate the transcription activation by Trl in larval imaginal discs. Taken together, these data indicate that Atro is a major Trl cofactor that functions to moderate developmental gene transcription.

CX-5461 is a DNA G-quadruplex stabilizer with selective lethality in BRCA1/2 deficient tumours.

G-quadruplex DNAs form four-stranded helical structures and are proposed to play key roles in different cellular processes. Targeting G-quadruplex DNAs for cancer treatment is a very promising prospect. Here, we show that CX-5461 is a G-quadruplex stabilizer, with specific toxicity against BRCA deficiencies in cancer cells and polyclonal patient-derived xenograft models, including tumours resistant to PARP inhibition. Exposure to CX-5461, and its related drug CX-3543, blocks replication forks and induces ssDNA gaps or breaks. The BRCA and NHEJ pathways are required for the repair of CX-5461 and CX-3543-induced DNA damage and failure to do so leads to lethality. These data strengthen the concept of G4 targeting as a therapeutic approach, specifically for targeting HR and NHEJ deficient cancers and other tumours deficient for DNA damage repair. CX-5461 is now in advanced phase I clinical trial for patients with BRCA1/2 deficient tumours (Canadian trial, NCT02719977, opened May 2016).

CLK-dependent exon recognition and conjoined gene formation revealed with a novel small molecule inhibitor.

CDC-like kinase phosphorylation of serine/arginine-rich proteins is central to RNA splicing reactions. Yet, the genomic network of CDC-like kinase-dependent RNA processing events remains poorly defined. Here, we explore the connectivity of genomic CDC-like kinase splicing functions by applying graduated, short-exposure, pharmacological CDC-like kinase inhibition using a novel small molecule (T3) with very high potency, selectivity, and cell-based stability. Using RNA-Seq, we define CDC-like kinase-responsive alternative splicing events, the large majority of which monotonically increase or decrease with increasing CDC-like kinase inhibition. We show that distinct RNA-binding motifs are associated with T3 response in skipped exons. Unexpectedly, we observe dose-dependent conjoined gene transcription, which is associated with motif enrichment in the last and second exons of upstream and downstream partners, respectively. siRNA knockdown of CLK2-associated genes significantly increases conjoined gene formation. Collectively, our results reveal an unexpected role for CDC-like kinase in conjoined gene formation, via regulation of 3'-end processing and associated splicing factors.The phosphorylation of serine/arginine-rich proteins by CDC-like kinase is a central regulatory mechanism for RNA splicing reactions. Here, the authors synthesize a novel small molecule CLK inhibitor and map CLK-responsive alternative splicing events and discover an effect on conjoined gene transcription.

Histological Transformation and Progression in Follicular Lymphoma: A Clonal Evolution Study.

BACKGROUND: Follicular lymphoma (FL) is an indolent, yet incurable B cell malignancy. A subset of patients experience an increased mortality rate driven by two distinct clinical end points: histological transformation and early progression after immunochemotherapy. The nature of tumor clonal dynamics leading to these clinical end points is poorly understood, and previously determined genetic alterations do not explain the majority of transformed cases or accurately predict early progressive disease. We contend that detailed knowledge of the expansion patterns of specific cell populations plus their associated mutations would provide insight into therapeutic strategies and disease biology over the time course of FL clinical histories. METHODS AND FINDINGS: Using a combination of whole genome sequencing, targeted deep sequencing, and digital droplet PCR on matched diagnostic and relapse specimens, we deciphered the constituent clonal populations in 15 transformation cases and 6 progression cases, and measured the change in clonal population abundance over time. We observed widely divergent patterns of clonal dynamics in transformed cases relative to progressed cases. Transformation specimens were generally composed of clones that were rare or absent in diagnostic specimens, consistent with dramatic clonal expansions that came to dominate the transformation specimens. This pattern was independent of time to transformation and treatment modality. By contrast, early progression specimens were composed of clones that were already present in the diagnostic specimens and exhibited only moderate clonal dynamics, even in the presence of immunochemotherapy. Analysis of somatic mutations impacting 94 genes was undertaken in an extension cohort consisting of 395 samples from 277 patients in order to decipher disrupted biology in the two clinical end points. We found 12 genes that were more commonly mutated in transformed samples than in the preceding FL tumors, including TP53, B2M, CCND3, GNA13, S1PR2, and P2RY8. Moreover, ten genes were more commonly mutated in diagnostic specimens of patients with early progression, including TP53, BTG1, MKI67, and XBP1. CONCLUSIONS: Our results illuminate contrasting modes of evolution shaping the clinical histories of transformation and progression. They have implications for interpretation of evolutionary dynamics in the context of treatment-induced selective pressures, and indicate that transformation and progression will require different clinical management strategies.

Divergent modes of clonal spread and intraperitoneal mixing in high-grade serous ovarian cancer.

We performed phylogenetic analysis of high-grade serous ovarian cancers (68 samples from seven patients), identifying constituent clones and quantifying their relative abundances at multiple intraperitoneal sites. Through whole-genome and single-nucleus sequencing, we identified evolutionary features including mutation loss, convergence of the structural genome and temporal activation of mutational processes that patterned clonal progression. We then determined the precise clonal mixtures comprising each tumor sample. The majority of sites were clonally pure or composed of clones from a single phylogenetic clade. However, each patient contained at least one site composed of polyphyletic clones. Five patients exhibited monoclonal and unidirectional seeding from the ovary to intraperitoneal sites, and two patients demonstrated polyclonal spread and reseeding. Our findings indicate that at least two distinct modes of intraperitoneal spread operate in clonal dissemination and highlight the distribution of migratory potential over clonal populations comprising high-grade serous ovarian cancers.

Clonal genotype and population structure inference from single-cell tumor sequencing.

Single-cell DNA sequencing has great potential to reveal the clonal genotypes and population structure of human cancers. However, single-cell data suffer from missing values and biased allelic counts as well as false genotype measurements owing to the sequencing of multiple cells. We describe the Single Cell Genotyper (https://bitbucket.org/aroth85/scg), an open-source software based on a statistical model coupled with a mean-field variational inference method, which can be used to address these problems and robustly infer clonal genotypes.

Weaver Syndrome-Associated EZH2 Protein Variants Show Impaired Histone Methyltransferase Function In Vitro.

Weaver syndrome (WS) is a rare congenital disorder characterized by generalized overgrowth, macrocephaly, specific facial features, accelerated bone age, intellectual disability, and susceptibility to cancers. De novo mutations in the enhancer of zeste homolog 2 (EZH2) have been shown to cause WS. EZH2 is a histone methyltransferase that acts as the catalytic agent of the polycomb-repressive complex 2 (PRC2) to maintain gene repression via methylation of lysine 27 on histone H3 (H3K27). Functional studies investigating histone methyltransferase activity of mutant EZH2 from various cancers have been reported, whereas WS-associated mutations remain poorly characterized. To investigate the role of EZH2 in WS, we performed functional studies using artificially assembled PRC2 complexes containing mutagenized human EZH2 that reflected the codon changes predicted from patients with WS. We found that WS-associated amino acid alterations reduce the histone methyltransferase function of EZH2 in this in vitro assay. Our results support the hypothesis that WS is caused by constitutional mutations in EZH2 that alter the histone methyltransferase function of PRC2. However, histone methyltransferase activities of different EZH2 variants do not appear to correlate directly with the phenotypic variability between WS patients and individuals with a common c.553G>C (p.Asp185His) polymorphism in EZH2.

The oncogenic roles of DICER1 RNase IIIb domain mutations in ovarian Sertoli-Leydig cell tumors.

DICER1, an endoribonuclease required for microRNA (miRNA) biogenesis, is essential for embryogenesis and the development of many organs including ovaries. We have recently identified somatic hotspot mutations in RNase IIIb domain of DICER1 in half of ovarian Sertoli-Leydig cell tumors, a rare class of sex-cord stromal cell tumors in young women. These hotspot mutations lost IIIb cleavage activity of DICER1 in vitro and failed to produce 5p-derived miRNAs in mouse Dicer1-null ES cells. However, the oncogenic potential of these hotspot DICER1 mutations has not been studied. Here, we further revealed that the global expression of 5p-derived miRNAs was dramatically reduced in ovarian Sertoli-Leydig cell tumors carrying DICER1 hotspot mutations compared with those without DICER1 hotspot mutation. The miRNA production defect was associated with the deregulation of genes controlling cell proliferation and the cell fate. Using an immortalized human granulosa cell line, SVOG3e, we determined that the D1709N-DICER1 hotspot mutation failed to produce 5p-derived miRNAs, deregulated the expression of several genes that control gonadal differentiation and cell proliferation, and promoted cell growth. Re-expression of let-7 significantly inhibited the growth of D1709N-DICER1 SVOG3e cells, accompanied by the suppression of key regulators of cell cycle control and ovarian gonad differentiation. Taken together, our data revealed that DICER1 hotspot mutations cause systemic loss of 5p-miRNAs that can both drive pseudodifferentiation of testicular elements and cause oncogenic transformation in the ovary.

A co-culture genome-wide RNAi screen with mammary epithelial cells reveals transmembrane signals required for growth and differentiation.

INTRODUCTION: The extracellular signals regulating mammary epithelial cell growth are of relevance to understanding the pathophysiology of mammary epithelia, yet they remain poorly characterized. In this study, we applied an unbiased approach to understanding the functional role of signalling molecules in several models of normal physiological growth and translated these results to the biological understanding of breast cancer subtypes. METHODS: We developed and utilized a cytogenetically normal clonal line of hTERT immortalized human mammary epithelial cells in a fibroblast-enhanced co-culture assay to conduct a genome-wide small interfering RNA (siRNA) screen for evaluation of the functional effect of silencing each gene. Our selected endpoint was inhibition of growth. In rigorous postscreen validation processes, including quantitative RT-PCR, to ensure on-target silencing, deconvolution of pooled siRNAs and independent confirmation of effects with lentiviral short-hairpin RNA constructs, we identified a subset of genes required for mammary epithelial cell growth. Using three-dimensional Matrigel growth and differentiation assays and primary human mammary epithelial cell colony assays, we confirmed that these growth effects were not limited to the 184-hTERT cell line. We utilized the METABRIC dataset of 1,998 breast cancer patients to evaluate both the differential expression of these genes across breast cancer subtypes and their prognostic significance. RESULTS: We identified 47 genes that are critically important for fibroblast-enhanced mammary epithelial cell growth. This group was enriched for several axonal guidance molecules and G protein-coupled receptors, as well as for the endothelin receptor PROCR. The majority of genes (43 of 47) identified in two dimensions were also required for three-dimensional growth, with HSD17B2, SNN and PROCR showing greater than tenfold reductions in acinar formation. Several genes, including PROCR and the neuronal pathfinding molecules EFNA4 and NTN1, were also required for proper differentiation and polarization in three-dimensional cultures. The 47 genes identified showed a significant nonrandom enrichment for differential expression among 10 molecular subtypes of breast cancer sampled from 1,998 patients. CD79A, SERPINH1, KCNJ5 and TMEM14C exhibited breast cancer subtype-independent overall survival differences. CONCLUSION: Diverse transmembrane signals are required for mammary epithelial cell growth in two-dimensional and three-dimensional conditions. Strikingly, we define novel roles for axonal pathfinding receptors and ligands and the endothelin receptor in both growth and differentiation.

Dynamics of genomic clones in breast cancer patient xenografts at single-cell resolution.

Human cancers, including breast cancers, comprise clones differing in mutation content. Clones evolve dynamically in space and time following principles of Darwinian evolution, underpinning important emergent features such as drug resistance and metastasis. Human breast cancer xenoengraftment is used as a means of capturing and studying tumour biology, and breast tumour xenografts are generally assumed to be reasonable models of the originating tumours. However, the consequences and reproducibility of engraftment and propagation on the genomic clonal architecture of tumours have not been systematically examined at single-cell resolution. Here we show, using deep-genome and single-cell sequencing methods, the clonal dynamics of initial engraftment and subsequent serial propagation of primary and metastatic human breast cancers in immunodeficient mice. In all 15 cases examined, clonal selection on engraftment was observed in both primary and metastatic breast tumours, varying in degree from extreme selective engraftment of minor (<5% of starting population) clones to moderate, polyclonal engraftment. Furthermore, ongoing clonal dynamics during serial passaging is a feature of tumours experiencing modest initial selection. Through single-cell sequencing, we show that major mutation clusters estimated from tumour population sequencing relate predictably to the most abundant clonal genotypes, even in clonally complex and rapidly evolving cases. Finally, we show that similar clonal expansion patterns can emerge in independent grafts of the same starting tumour population, indicating that genomic aberrations can be reproducible determinants of evolutionary trajectories. Our results show that measurement of genomically defined clonal population dynamics will be highly informative for functional studies using patient-derived breast cancer xenoengraftment.

TITAN: inference of copy number architectures in clonal cell populations from tumor whole-genome sequence data.

The evolution of cancer genomes within a single tumor creates mixed cell populations with divergent somatic mutational landscapes. Inference of tumor subpopulations has been disproportionately focused on the assessment of somatic point mutations, whereas computational methods targeting evolutionary dynamics of copy number alterations (CNA) and loss of heterozygosity (LOH) in whole-genome sequencing data remain underdeveloped. We present a novel probabilistic model, TITAN, to infer CNA and LOH events while accounting for mixtures of cell populations, thereby estimating the proportion of cells harboring each event. We evaluate TITAN on idealized mixtures, simulating clonal populations from whole-genome sequences taken from genomically heterogeneous ovarian tumor sites collected from the same patient. In addition, we show in 23 whole genomes of breast tumors that the inference of CNA and LOH using TITAN critically informs population structure and the nature of the evolving cancer genome. Finally, we experimentally validated subclonal predictions using fluorescence in situ hybridization (FISH) and single-cell sequencing from an ovarian cancer patient sample, thereby recapitulating the key modeling assumptions of TITAN.

Impact of MLL5 expression on decitabine efficacy and DNA methylation in acute myeloid leukemia.

Hypomethylating agents are widely used in patients with myelodysplastic syndromes and unfit patients with acute myeloid leukemia. However, it is not well understood why only some patients respond to hypomethylating agents. We found previously that the effect of decitabine on hematopoietic stem cell viability differed between Mll5 wild-type and null cells. We, therefore, investigated the role of MLL5 expression levels on outcome of acute myeloid leukemia patients who were treated with decitabine. MLL5 above the median expression level predicted longer overall survival independent of DNMT3A mutation status in bivariate analysis (median overall survival for high vs. low MLL5 expression 292 vs. 167 days; P=0.026). In patients who received three or more courses decitabine, high MLL5 expression and wild-type DNMT3A independently predicted improved overall survival (median overall survival for high vs. low MLL5 expression 468 vs. 243 days; P=0.012). In transformed murine cells, loss of Mll5 was associated with resistance to low-dose decitabine, less global DNA methylation in promoter regions, and reduced DNA demethylation upon decitabine treatment. Together, these data support our clinical observation of improved outcome in decitabine-treated patients who express MLL5 at high levels, and suggest a mechanistic role of MLL5 in the regulation of DNA methylation.

A transgenic mouse model demonstrating the oncogenic role of mutations in the polycomb-group gene EZH2 in lymphomagenesis.

The histone methyltransferase EZH2 is frequently mutated in germinal center-derived diffuse large B-cell lymphoma and follicular lymphoma. To further characterize these EZH2 mutations in lymphomagenesis, we generated a mouse line where EZH2(Y641F) is expressed from a lymphocyte-specific promoter. Spleen cells isolated from the transgenic mice displayed a global increase in trimethylated H3K27, but the mice did not show an increased tendency to develop lymphoma. As EZH2 mutations often coincide with other mutations in lymphoma, we combined the expression of EZH2(Y641F) by crossing these transgenic mice with Eµ-Myc transgenic mice. We observed a dramatic acceleration of lymphoma development in this combination model of Myc and EZH2(Y641F). The lymphomas show histologic features of high-grade disease with a shift toward a more mature B-cell phenotype, increased cycling and gene expression, and epigenetic changes involving important pathways in B-cell regulation and function. Furthermore, they initiate disease in secondary recipients. In summary, EZH2(Y641F) can collaborate with Myc to accelerate lymphomagenesis demonstrating a cooperative role of EZH2 mutations in oncogenesis. This murine lymphoma model provides a new tool to study global changes in the epigenome caused by this frequent mutation and a promising model system for testing novel treatments.

PyClone: statistical inference of clonal population structure in cancer.

We introduce PyClone, a statistical model for inference of clonal population structures in cancers. PyClone is a Bayesian clustering method for grouping sets of deeply sequenced somatic mutations into putative clonal clusters while estimating their cellular prevalences and accounting for allelic imbalances introduced by segmental copy-number changes and normal-cell contamination. Single-cell sequencing validation demonstrates PyClone's accuracy.