Abstract:Abstract Cancer precision medicine is based on the idea that specific genomic features of a human tumor are associated with its drug response. Therefore, by assessing the genomic features of a tumor, its vulnerabilities could be predicted, and an optimal treatment could then be tailored. In order to achieve this goal, there is a burning need to develop a comprehensive “dependency map” that would enable to predict the “Achilles heel” of each tumor based on its genomic features. In an attempt to develop such cancer dependency map, the cancer community heavily relies on cancer models, such as genetically engineered mouse models (GEMMs), patient-derived xenografts (PDXs), and human cell lines (CLs). Importantly, xenografts and cell lines are used not just for cancer research, but also to predict the response of a specific patient to anti-cancer drugs (i.e., as tumor “avatars”). As cancer models are invaluable for cancer research, much effort has been recently invested in their genomic and phenotypic profiling. However, like any biological system, cancer models evolve. Whereas tumor heterogeneity and genomic instability have received much attention, the heterogeneity and instability of cancer models - and how the genomic evolution of these models affects their use in cancer research and cancer precision medicine - remain under-explored.In my postdoctoral work, I have studied the extent of genomic evolution in multiple cancer models, its biological origins and its functional implications. This work was comprised of three independent studies, which dissected various aspects of genomic evolution in three major cancer model systems: GEMMs, PDXs and CLs. In this talk, I will describe the main conclusions of my studies, and highlight important themes that emerge from their synthesis.In a first study (Ben-David et al. Nature Communications 2016), we generated the first comprehensive catalogue of copy-number alterations (CNAs) in cancer GEMMs. Mining this resource, we found that chromosomal aberrations accumulated late during breast tumorigenesis, and observed marked differences in CNA prevalence between mouse mammary tumours initiated with distinct drivers. Some aberrations were recurrent and unique to specific GEMMs, suggesting distinct driver-dependent routes to tumorigenesis. Using synteny-based oncogenomics of mouse and human data, we narrowed down the critical region of interest in one of the most recurrent chromosomal changes in breast cancer (loss of chromosome 1p), and identified a gene (Sfn) that cooperates with Erbb2 during breast cancer tumorigenesis. We experimentally validated that loss of Stratifin (SFN) promotes HER2-induced tumorigenesis. This study demonstrates how the natural genomic evolution during tumorigenesis in GEMMs can inform the identification of oncogenic networks that are important for the pathogenesis of human cancer.In a second study (Ben-David et al. Nature Genetics 2017), we followed the genomic evolution of PDXs throughout their derivation and in vivo propagation. We observed extensive evolution of the PDX copy-number landscapes, mostly driven by strong clonal dynamics leading to the expansion of pre-existing minor subclones. Clonal dynamics were strongest during derivation and early propagation of the models, and attenuated at later passages. Reproducible changes were observed across independent PDXs generated from the same primary tumor, indicating the involvement of selection, rather than mere genetic drift. Importantly, the rate of genomic evolution in PDXs was similar to that observed in patient-derived CLs. Tumor evolution in PDXs and in patients followed distinct trajectories. Specifically, aneuploidy events recurrently observed in primary tumors gradually disappeared in PDXs. The genomic stability of PDXs was associated with their response to chemotherapy and targeted drugs, suggesting that the genomic evolution of PDXs could affect the faithfulness of PDX-based prediction of patients' drug response.In a third study (Ben-David et al. Nature 2018), we compared genomic analyses of 106 CLs grown in two laboratories and revealed extensive genetic diversity. Follow-up comprehensive genomic characterization of 27 strains of the common breast cancer cell line MCF7 uncovered rapid genetic diversification. Similar results were obtained with multiple strains of 13 additional CLs. Genetic changes were associated with differential activation of gene expression programs and marked differences in cell morphology and proliferation. Barcoding experiments showed that CL evolution occurred as a result of positive clonal selection that is highly sensitive to culture conditions. Analyses of single cell-derived clones showed that ongoing instability quickly translated into CL heterogeneity. Functionally, testing of the 27 MCF7 strains against 321 anti-cancer compounds uncovered strikingly disparate drug response: over 80% of compounds that strongly inhibited some CL strains were completely inactive in others. These findings have broad implications for cell line use in basic cancer research and cancer precision medicine. Importantly, we suggest practical ways to mitigate the risks posed by CL genomic evolution and to facilitate maximally-reproducible cell line-based research. Moreover, we propose how to constructively build upon this phenomenon in future studies (e.g. to study of types of genetic variation that cannot be readily introduced experimentally, such as large chromosomal changes).In a fourth study (Ben-David et al. Submitted), we performed a deep genomic analysis of human cancer cell lines before and after the introduction of Cas9. We found upregulation of the p53 pathway upon Cas9 introduction, specifically in TP53-WT cell lines. This upregulation was confirmed at the mRNA and protein levels. Moreover, we found elevated levels of DNA repair transcriptional signatures in the Cas9 cell lines, and confirmed that the introduction of Cas9 induced DNA damage using immunofluorescence quantification of γH2AX. Genetic characterization of 42 cell line pairs showed that Cas9 introduction could lead to the emergence and to the expansion of p53-inactivating mutations. Competition experiments with isogenic TP53-WT/TP53-null cell lines confirmed that Cas9 introduction accelerated selection for p53-inactivating mutation. Lastly, we compared Cas9 activity across 719 human cell lines and found that Cas9 was significantly less active in TP53-WT cell lines, in line with functional p53 being a barrier to the efficient expression of Cas9 in human cells. These findings have broad implications for the proper use of CRISPR/Cas9 genome editing in basic research and in clinical applications. They also demonstrate the profound functional consequences of in vitro genomic evolution following a defined, common selection pressure.In summary, my postdoctoral studies revealed a previously under-appreciated genomic instability in multiple cancer models, and characterized the implications of this instability for the use of such models in cancer research. Combined, this body of work exposed extensive genetic evolution of cancer models, distinct from that seen in patients, which results in gene expression changes and disparate drug response. Our studies of the stability and faithfulness of cancer models improve our biological understanding of the advantages and limitations of the models, and can help guide their proper application. Importantly, the genomic instability of cancer models also presents a unique opportunity to use these models in novel and creative ways (Ben-David at al. Nature Reviews Cancer 2019; Ben-David and Amon Nature Reviews Genetics 2019). This work thus highlights both the perils and the opportunities of the natural evolution of cancer models. Citation Format: Uri Ben-David. Genomic evolution of cancer models: Perils and opportunities [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr NG02.

Abstract 4307: Quantifying the cell-intrinsic and cell-extrinsic growth effects of cancer-causing mutations in vivo

Abstract 11: Analysis of 30 million tumor time course reveals that inactivating tumor suppressors produces complex gene-specific, time-dependent growth dynamics

In Vivo Functional Genomics Identifies Essentiality of Potassium Homeostasis in Medulloblastoma

Abstract A014: Models of neighbors and space: frequency and density dependent dynamics of tumor evolution

Abstract 4676: Deconvolution of extrinsic and intrinsic factors influencing tumor growth in solid cancer PDX models to support design, analysis and application of PDX-based pharmacology studies

Abstract IA014: Detecting epistatic effects among somatic cancer mutations across oncogenesis

Abstract IA010: Unraveling mutation and selection to better understand early cancer evolution

Functional genomic landscape of T-cell mediated cytotoxicity.

Abstract 859: Simulations of tumor heterogeneity impact on treatment response using a mechanistic model of EGFR-mutant lung adenocarcinoma

Abstract 2286: Clonal sweep dynamics: A marker of tumor fitness and predictor of clinical outcomes

Abstract C030: Evolutionary and epistatic analyses reveal genic interactions with KRAS during malignant progression of pancreatic ductal adenocarcinoma

Abstract IA024: Unraveling potential genetic drivers of intratumor transcriptional heterogeneity and phenotypic plasticity in malignant cells

Abstract A027: Detecting and modeling frequency-dependent evolution in serially passaged primary triple negative breast cancer patient derived xenografts

Abstract PR013: An agent-based modelling framework to study cell plasticity in non-small cell lung cancer

Accumulation of driver and passenger mutations during tumor progression

Abstract 6930: Spatial analysis of multi-ancestral tumors to understand how intercellular signaling impacts tumor invasion

Abstract A005: Quantification of cancer effect sizes of somatic variants reveals stage-dependent evolution in primary and metastatic skin cutaneous melanoma

Abstract NG02: Genomic evolution of cancer models: Perils and opportunities

Abstract C057: Delineating tissue-specific oncogene selection with autochthonous CRISPR activation screening

An evolutionary model of tumor cell kinetics and the emergence of molecular heterogeneity driving Gompertzian growth

Abstract 2928: Solving the cancer mutation conundrum: A single cell, massively parallel approach for cancer mutation discovery, genome modelling and functional characterization