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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Somatic mutations at EZH2 Y641 act dominantly through a mechanism of selectively altered PRC2 catalytic activity, to increase H3K27 trimethylation.

Next-generation sequencing of follicular lymphoma and diffuse-large B-cell lymphoma has revealed frequent somatic, heterozygous Y641 mutations in the histone methyltransferase EZH2. Heterozygosity and the presence of equal quantities of both mutant and wild-type mRNA and expressed protein suggest a dominant mode of action. Surprisingly, B-cell lymphoma cell lines and lymphoma samples harboring heterozygous EZH2(Y641) mutations have increased levels of histone H3 Lys-27-specific trimethylation (H3K27me3). Expression of EZH2(Y641F/N) mutants in cells with EZH2(WT) resulted in an increase of H3K27me3 levels in vivo. Structural modeling of EZH2(Y641) mutants suggests a "Tyr/Phe switch" model whereby structurally neutral, nontyrosine residues at position 641 would decrease affinity for unmethylated and monomethylated H3K27 substrates and potentially favor trimethylation. We demonstrate, using in vitro enzyme assays of reconstituted PRC2 complexes, that Y641 mutations result in a decrease in monomethylation and an increase in trimethylation activity of the enzyme relative to the wild-type enzyme. This represents the first example of a disease-associated gain-of-function mutation in a histone methyltransferase, whereby somatic EZH2 Y641 mutations in lymphoma act dominantly to increase, rather than decrease, histone methylation. The dominant mode of action suggests that allele-specific EZH2 inhibitors should be a future therapeutic strategy for this disease.

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.

The testosterone-dependent and independent transcriptional networks in the hypothalamus of Gpr54 and Kiss1 knockout male mice are not fully equivalent.

BACKGROUND: Humans and mice with loss of function mutations in GPR54 (KISS1R) or kisspeptin do not progress through puberty, caused by a failure to release GnRH. The transcriptional networks regulated by these proteins in the hypothalamus have yet to be explored by genome-wide methods. RESULTS: We show here, using 1 million exon mouse arrays (Exon 1.0 Affymetrix) and quantitative polymerase chain reaction (QPCR) validation to analyse microdissected hypothalamic tissue from Gpr54 and Kiss1 knockout mice, the extent of transcriptional regulation in the hypothalamus. The sensitivity to detect important transcript differences in microdissected RNA was confirmed by the observation of counter-regulation of Kiss1 expression in Gpr54 knockouts and confirmed by immunohistochemistry (IHC). Since Gpr54 and Kiss1 knockout animals are effectively pre-pubertal with low testosterone (T) levels, we also determined which of the validated transcripts were T-responsive and which varied according to genotype alone. We observed four types of transcriptional regulation (i) genotype only dependent regulation, (ii) T only dependent regulation, (iii) genotype and T-dependent regulation with interaction between these variables, (iv) genotype and T-dependent regulation with no interaction between these variables. The results implicate for the first time several transcription factors (e.g. Npas4, Esr2), proteases (Klk1b22), and the orphan 10-transmembrane transporter TMEM144 in the biology of GPR54/kisspeptin function in the hypothalamus. We show for the neuronal activity regulated transcription factor NPAS4, that distinct protein over-expression is seen in the hypothalamus and hippocampus in Gpr54 knockout mice. This links for the first time the hypothalamic-gonadal axis with this important regulator of inhibitory synapse formation. Similarly we confirm TMEM144 up-regulation in the hypothalamus by RNA in situ hybridization and western blot. CONCLUSIONS: Taken together, global transcriptional profiling shows that loss of GPR54 and kisspeptin are not fully equivalent in the mouse hypothalamus.

Loss of MLL5 results in pleiotropic hematopoietic defects, reduced neutrophil immune function, and extreme sensitivity to DNA demethylation.

MLL5 is a divergent member of the Drosophila Trithorax-related (SET) domain and plant homeodomain (PHD) domain-containing chromatin regulators that are involved in the regulation of transcriptional "memory" during differentiation. Human MLL5 is located on chromosome 7q22, which frequently is deleted in myeloid leukemias, suggesting a possible role in hemopoiesis. To address this question, we generated a loss-of-function allele (Mll5(tm1Apa)) in the murine Mll5 locus. Unlike other Mll genes, Mll5(tm1Apa) homozygous mice are viable but display defects in immunity and hematopoiesis. First, Mll5(tm1Apa) homozygous mice show increased susceptibility to spontaneous eye infections, associated with a cell-autonomous impairment of neutrophil function. Second, Mll5(tm1Apa/tm1Apa) mice exhibit a mild impairment of erythropoiesis. Third, Mll5(tm1Apa/tm1Apa) hematopoietic stem cells (HSCs) have impaired competitive repopulating capacity both under normal conditions and when subjected to self-renewal stimulation by NUP98-HOXA10. Fourth, Mll5(tm1Apa) homozygous HSCs show a dramatic sensitivity to DNA demethylation-induced differentiation (5-azadeoxycytidine). Taken together, our data show that MLL5 is involved in terminal myeloid differentiation and the regulation of HSC self-renewal by a mechanism that involves DNA methylation. These data warrant investigation of MLL5 expression levels as a predictive marker of demethylating-agent response in patients with myelodysplastic syndromes and leukemias and identify MLL5 as a key regulator of normal hematopoiesis.

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.

Novel function of the cyclin A binding site of E2F in regulating p53-induced apoptosis in response to DNA damage.

We demonstrate here that the E2F1 induced by DNA damage can bind to and promote the apoptotic function of p53 via the cyclin A binding site of E2F1. This function of E2F1 does not require its DP-1 binding, DNA binding, or transcriptional activity and is independent of mdm2. All the cyclin A binding E2F family members can interact and cooperate with p53 to induce apoptosis. This suggests a novel role for E2F in regulating apoptosis in response to DNA damage. Cyclin A, but not cyclin E, prevents E2F1 from interacting and cooperating with p53 to induce apoptosis. However, in response to DNA damage, cyclin A levels decrease, with a concomitant increase in E2F1-p53 complex formation. These results suggest that the binding of E2F1 to p53 can specifically stimulate the apoptotic function of p53 in response to DNA damage.

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.