Next-generation characterization of the Cancer Cell Line Encyclopedia.

Large panels of comprehensively characterized human cancer models, including the Cancer Cell Line Encyclopedia (CCLE), have provided a rigorous framework with which to study genetic variants, candidate targets, and small-molecule and biological therapeutics and to identify new marker-driven cancer dependencies. To improve our understanding of the molecular features that contribute to cancer phenotypes, including drug responses, here we have expanded the characterizations of cancer cell lines to include genetic, RNA splicing, DNA methylation, histone H3 modification, microRNA expression and reverse-phase protein array data for 1,072 cell lines from individuals of various lineages and ethnicities. Integration of these data with functional characterizations such as drug-sensitivity, short hairpin RNA knockdown and CRISPR-Cas9 knockout data reveals potential targets for cancer drugs and associated biomarkers. Together, this dataset and an accompanying public data portal provide a resource for the acceleration of cancer research using model cancer cell lines.

Contribution and clinical relevance of germline variation to the cancer transcriptome.

BACKGROUND: Somatic alterations in the cancer genome, some of which are associated with changes in gene expression, have been characterized in multiple studies across diverse cancer types. However, less is known about germline variants that influence tumor biology by shaping the cancer transcriptome. METHODS: We performed expression quantitative trait loci (eQTL) analyses using multi-dimensional data from The Cancer Genome Atlas to explore the role of germline variation in mediating the cancer transcriptome. After accounting for associations between somatic alterations and gene expression, we determined the contribution of inherited variants to the cancer transcriptome relative to that of somatic variants. Finally, we performed an interaction analysis using estimates of tumor cellularity to identify cell type-restricted eQTLs. RESULTS: The proportion of genes with at least one eQTL varied between cancer types, ranging between 0.8% in melanoma to 28.5% in thyroid cancer and was correlated more strongly with intratumor heterogeneity than with somatic alteration rates. Although contributions to variance in gene expression was low for most genes, some eQTLs accounted for more than 30% of expression of proximal genes. We identified cell type-restricted eQTLs in genes known to be cancer drivers including LPP and EZH2 that were associated with disease-specific mortality in TCGA but not associated with disease risk in published GWAS. Together, our results highlight the need to consider germline variation in interpreting cancer biology beyond risk prediction.

Correction of copy number induced false positives in CRISPR screens.

Cell autonomous cancer dependencies are now routinely identified using CRISPR loss-of-function viability screens. However, a bias exists that makes it difficult to assess the true essentiality of genes located in amplicons, since the entire amplified region can exhibit lethal scores. These false-positive hits can either be discarded from further analysis, which in cancer models can represent a significant number of hits, or methods can be developed to rescue the true-positives within amplified regions. We propose two methods to rescue true positive hits in amplified regions by correcting for this copy number artefact. The Local Drop Out (LDO) method uses the relative lethality scores within genomic regions to assess true essentiality and does not require additional orthogonal data (e.g. copy number value). LDO is meant to be used in screens covering a dense region of the genome (e.g. a whole chromosome or the whole genome). The General Additive Model (GAM) method models the screening data as a function of the known copy number values and removes the systematic effect from the measured lethality. GAM does not require the same density as LDO, but does require prior knowledge of the copy number values. Both methods have been developed with single sample experiments in mind so that the correction can be applied even in smaller screens. Here we demonstrate the efficacy of both methods at removing the copy number effect and rescuing hits from some of the amplified regions. We estimate a 70-80% decrease of false positive hits with either method in regions of high copy number compared to no correction.

Integrating common and rare genetic variation in diverse human populations.

Despite great progress in identifying genetic variants that influence human disease, most inherited risk remains unexplained. A more complete understanding requires genome-wide studies that fully examine less common alleles in populations with a wide range of ancestry. To inform the design and interpretation of such studies, we genotyped 1.6 million common single nucleotide polymorphisms (SNPs) in 1,184 reference individuals from 11 global populations, and sequenced ten 100-kilobase regions in 692 of these individuals. This integrated data set of common and rare alleles, called 'HapMap 3', includes both SNPs and copy number polymorphisms (CNPs). We characterized population-specific differences among low-frequency variants, measured the improvement in imputation accuracy afforded by the larger reference panel, especially in imputing SNPs with a minor allele frequency of

Mapping and sequencing of structural variation from eight human genomes.

Genetic variation among individual humans occurs on many different scales, ranging from gross alterations in the human karyotype to single nucleotide changes. Here we explore variation on an intermediate scale--particularly insertions, deletions and inversions affecting from a few thousand to a few million base pairs. We employed a clone-based method to interrogate this intermediate structural variation in eight individuals of diverse geographic ancestry. Our analysis provides a comprehensive overview of the normal pattern of structural variation present in these genomes, refining the location of 1,695 structural variants. We find that 50% were seen in more than one individual and that nearly half lay outside regions of the genome previously described as structurally variant. We discover 525 new insertion sequences that are not present in the human reference genome and show that many of these are variable in copy number between individuals. Complete sequencing of 261 structural variants reveals considerable locus complexity and provides insights into the different mutational processes that have shaped the human genome. These data provide the first high-resolution sequence map of human structural variation--a standard for genotyping platforms and a prelude to future individual genome sequencing projects.

Isocitrate dehydrogenase (IDH) mutations promote a reversible ZEB1/microRNA (miR)-200-dependent epithelial-mesenchymal transition (EMT).

Mutations in the genes encoding isocitrate dehydrogenase 1 and 2 (IDH1/2) occur in a variety of tumor types, resulting in production of the proposed oncometabolite, 2-hydroxyglutarate (2-HG). How mutant IDH and 2-HG alter signaling pathways to promote cancer, however, remains unclear. Additionally, there exist relatively few cell lines with IDH mutations. To examine the effect of endogenous IDH mutations and 2-HG, we created a panel of isogenic epithelial cell lines with either wild-type IDH1/2 or clinically relevant IDH1/2 mutations. Differences were noted in the ability of IDH mutations to cause robust 2-HG accumulation. IDH1/2 mutants that produce high levels of 2-HG cause an epithelial-mesenchymal transition (EMT)-like phenotype, characterized by changes in EMT-related gene expression and cellular morphology. 2-HG is sufficient to recapitulate aspects of this phenotype in the absence of an IDH mutation. In the cells types examined, mutant IDH-induced EMT is dependent on up-regulation of the transcription factor ZEB1 and down-regulation of the miR-200 family of microRNAs. Furthermore, sustained knockdown of IDH1 in IDH1 R132H mutant cells is sufficient to reverse many characteristics of EMT, demonstrating that continued expression of mutant IDH is required to maintain this phenotype. These results suggest mutant IDH proteins can reversibly deregulate discrete signaling pathways that contribute to tumorigenesis.

Efficiency and power as a function of sequence coverage, SNP array density, and imputation.

High coverage whole genome sequencing provides near complete information about genetic variation. However, other technologies can be more efficient in some settings by (a) reducing redundant coverage within samples and (b) exploiting patterns of genetic variation across samples. To characterize as many samples as possible, many genetic studies therefore employ lower coverage sequencing or SNP array genotyping coupled to statistical imputation. To compare these approaches individually and in conjunction, we developed a statistical framework to estimate genotypes jointly from sequence reads, array intensities, and imputation. In European samples, we find similar sensitivity (89%) and specificity (99.6%) from imputation with either 1× sequencing or 1 M SNP arrays. Sensitivity is increased, particularly for low-frequency polymorphisms (MAF < 5%), when low coverage sequence reads are added to dense genome-wide SNP arrays--the converse, however, is not true. At sites where sequence reads and array intensities produce different sample genotypes, joint analysis reduces genotype errors and identifies novel error modes. Our joint framework informs the use of next-generation sequencing in genome wide association studies and supports development of improved methods for genotype calling.

Accurately assessing the risk of schizophrenia conferred by rare copy-number variation affecting genes with brain function.

Investigators have linked rare copy number variation (CNVs) to neuropsychiatric diseases, such as schizophrenia. One hypothesis is that CNV events cause disease by affecting genes with specific brain functions. Under these circumstances, we expect that CNV events in cases should impact brain-function genes more frequently than those events in controls. Previous publications have applied "pathway" analyses to genes within neuropsychiatric case CNVs to show enrichment for brain-functions. While such analyses have been suggestive, they often have not rigorously compared the rates of CNVs impacting genes with brain function in cases to controls, and therefore do not address important confounders such as the large size of brain genes and overall differences in rates and sizes of CNVs. To demonstrate the potential impact of confounders, we genotyped rare CNV events in 2,415 unaffected controls with Affymetrix 6.0; we then applied standard pathway analyses using four sets of brain-function genes and observed an apparently highly significant enrichment for each set. The enrichment is simply driven by the large size of brain-function genes. Instead, we propose a case-control statistical test, cnv-enrichment-test, to compare the rate of CNVs impacting specific gene sets in cases versus controls. With simulations, we demonstrate that cnv-enrichment-test is robust to case-control differences in CNV size, CNV rate, and systematic differences in gene size. Finally, we apply cnv-enrichment-test to rare CNV events published by the International Schizophrenia Consortium (ISC). This approach reveals nominal evidence of case-association in neuronal-activity and the learning gene sets, but not the other two examined gene sets. The neuronal-activity genes have been associated in a separate set of schizophrenia cases and controls; however, testing in independent samples is necessary to definitively confirm this association. Our method is implemented in the PLINK software package.

Rapid assessment of genetic ancestry in populations of unknown origin by genome-wide genotyping of pooled samples.

As we move forward from the current generation of genome-wide association (GWA) studies, additional cohorts of different ancestries will be studied to increase power, fine map association signals, and generalize association results to additional populations. Knowledge of genetic ancestry as well as population substructure will become increasingly important for GWA studies in populations of unknown ancestry. Here we propose genotyping pooled DNA samples using genome-wide SNP arrays as a viable option to efficiently and inexpensively estimate admixture proportion and identify ancestry informative markers (AIMs) in populations of unknown origin. We constructed DNA pools from African American, Native Hawaiian, Latina, and Jamaican samples and genotyped them using the Affymetrix 6.0 array. Aided by individual genotype data from the African American cohort, we established quality control filters to remove poorly performing SNPs and estimated allele frequencies for the remaining SNPs in each panel. We then applied a regression-based method to estimate the proportion of admixture in each cohort using the allele frequencies estimated from pooling and populations from the International HapMap Consortium as reference panels, and identified AIMs unique to each population. In this study, we demonstrated that genotyping pooled DNA samples yields estimates of admixture proportion that are both consistent with our knowledge of population history and similar to those obtained by genotyping known AIMs. Furthermore, through validation by individual genotyping, we demonstrated that pooling is quite effective for identifying SNPs with large allele frequency differences (i.e., AIMs) and that these AIMs are able to differentiate two closely related populations (HapMap JPT and CHB).

Common genetic variants and modification of penetrance of BRCA2-associated breast cancer.

The considerable uncertainty regarding cancer risks associated with inherited mutations of BRCA2 is due to unknown factors. To investigate whether common genetic variants modify penetrance for BRCA2 mutation carriers, we undertook a two-staged genome-wide association study in BRCA2 mutation carriers. In stage 1 using the Affymetrix 6.0 platform, 592,163 filtered SNPs genotyped were available on 899 young (<40 years) affected and 804 unaffected carriers of European ancestry. Associations were evaluated using a survival-based score test adjusted for familial correlations and stratified by country of the study and BRCA2*6174delT mutation status. The genomic inflation factor (λ) was 1.011. The stage 1 association analysis revealed multiple variants associated with breast cancer risk: 3 SNPs had p-values<10(-5) and 39 SNPs had p-values<10(-4). These variants included several previously associated with sporadic breast cancer risk and two novel loci on chromosome 20 (rs311499) and chromosome 10 (rs16917302). The chromosome 10 locus was in ZNF365, which contains another variant that has recently been associated with breast cancer in an independent study of unselected cases. In stage 2, the top 85 loci from stage 1 were genotyped in 1,264 cases and 1,222 controls. Hazard ratios (HR) and 95% confidence intervals (CI) for stage 1 and 2 were combined and estimated using a retrospective likelihood approach, stratified by country of residence and the most common mutation, BRCA2*6174delT. The combined per allele HR of the minor allele for the novel loci rs16917302 was 0.75 (95% CI 0.66-0.86, ) and for rs311499 was 0.72 (95% CI 0.61-0.85, ). FGFR2 rs2981575 had the strongest association with breast cancer risk (per allele HR = 1.28, 95% CI 1.18-1.39, ). These results indicate that SNPs that modify BRCA2 penetrance identified by an agnostic approach thus far are limited to variants that also modify risk of sporadic BRCA2 wild-type breast cancer.

Generating a Murine PTEN Null Cell Line to Discover the Key Role of p110ß-PAK1 in Castration-Resistant Prostate Cancer Invasion.

Although androgen deprivation treatment often effectively decreases prostate cancer, incurable metastatic castration-resistant prostate cancer (CRPC) eventually occurs. It is important to understand how CRPC metastasis progresses, which is not clearly defined. The loss of PTEN, a phosphatase to dephosphorylate phosphatidylinositol 3,4,5-trisphosphate in the PI3K pathway, occurs in up to 70% to 80% of CRPC. We generated a mouse androgen-independent prostate cancer cell line (PKO) from PTEN null and Hi-Myc transgenic mice in C57BL/6 background. We confirmed that this PKO cell line has an activated PI3K pathway and can metastasize into the femur and tibia of immunodeficient nude and immunocompetent C57BL/6 mice. In vitro, we found that androgen deprivation significantly enhanced PKO cell migration/invasion via the p110ß isoform-depended PAK1-MAPK activation. Inhibition of the p110ß-PAK1 axis significantly decreased prostate cancer cell migration/invasion. Of note, our analysis using clinical samples showed that PAK1 is more activated in CRPC than in advanced prostate cancer; high PAK1/phosphorylated-PAK1 levels are associated with decreased survival rates in patients with CRPC. All the information suggests that this cell line reflects the characteristics of CRPC cells and can be applied to dissect the mechanism of CRPC initiation and progression. This study also shows that PAK1 is a potential target for CRPC treatment. IMPLICATIONS: This study uses a newly generated PTEN null prostate cancer cell line to define a critical functional role of p110ß-PAK1 in CRPC migration/invasion. This study also shows that the p110ß-PAK1 axis can potentially be a therapeutic target in CRPC metastasis.

Expressed Barcoding Enables High-Resolution Tracking of the Evolution of Drug Tolerance.

For a majority of patients with non-small cell lung cancer with EGFR mutations, treatment with EGFR inhibitors (EGFRi) induces a clinical response. Despite this initial reduction in tumor size, residual disease persists that leads to disease relapse. Elucidating the preexisting biological differences between sensitive cells and surviving drug-tolerant persister cells and deciphering how drug-tolerant cells evolve in response to treatment could help identify strategies to improve the efficacy of EGFRi. In this study, we tracked the origins and clonal evolution of drug-tolerant cells at a high resolution by using an expressed barcoding system coupled with single-cell RNA sequencing. This platform enabled longitudinal profiling of gene expression and drug sensitivity in response to EGFRi across a large number of clones. Drug-tolerant cells had higher expression of key survival pathways such as YAP and EMT at baseline and could also differentially adapt their gene expression following EGFRi treatment compared with sensitive cells. In addition, drug combinations targeting common downstream components (MAPK) or orthogonal factors (chemotherapy) showed greater efficacy than EGFRi alone, which is attributable to broader targeting of the heterogeneous EGFRi-tolerance mechanisms present in tumors. Overall, this approach facilitates thorough examination of clonal evolution in response to therapy that could inform the development of improved diagnostic approaches and treatment strategies for targeting drug-tolerant cells. SIGNIFICANCE: The evolution and heterogeneity of EGFR inhibitor tolerance are identified in a large number of clones at enhanced cellular and temporal resolution using an expressed barcode technology coupled with single-cell RNA sequencing.

Blocking PI3K p110ß Attenuates Development of PTEN-Deficient Castration-Resistant Prostate Cancer.

A common outcome of androgen deprivation in prostate cancer therapy is disease relapse and progression to castration-resistant prostate cancer (CRPC) via multiple mechanisms. To gain insight into the recent clinical findings that highlighted genomic alterations leading to hyperactivation of PI3K, we examined the roles of the commonly expressed p110 catalytic isoforms of PI3K in a murine model of Pten-null invasive CRPC. While blocking p110α had negligible effects in the development of Pten-null invasive CRPC, either genetic or pharmacologic perturbation of p110ß dramatically slowed CRPC initiation and progression. Once fully established, CRPC tumors became partially resistant to p110ß inhibition, indicating the acquisition of new dependencies. Driven by our genomic analyses highlighting potential roles for the p110ß/RAC/PAK1 and ß-catenin pathways in CRPC, we found that combining p110ß with RAC/PAK1 or tankyrase inhibitors significantly reduced the growth of murine and human CRPC organoids in vitro and in vivo. Because p110ß activity is dispensable for most physiologic processes, our studies support novel therapeutic strategies both for preventing disease progression into CRPC and for treating CRPC. IMPLICATIONS: This work establishes p110ß as a promising target for preventing the progression of primary PTEN-deficient prostate tumors to CRPC, and for treating established CRPC in combination with RAC/PAK1 or tankyrase inhibitors.

The efficacy of detecting variants with small effects on the Affymetrix 6.0 platform using pooled DNA.

Genome-wide genotyping of a cohort using pools rather than individual samples has long been proposed as a cost-saving alternative for performing genome-wide association (GWA) studies. However, successful disease gene mapping using pooled genotyping has thus far been limited to detecting common variants with large effect sizes, which tend not to exist for many complex common diseases or traits. Therefore, for DNA pooling to be a viable strategy for conducting GWA studies, it is important to determine whether commonly used genome-wide SNP array platforms such as the Affymetrix 6.0 array can reliably detect common variants of small effect sizes using pooled DNA. Taking obesity and age at menarche as examples of human complex traits, we assessed the feasibility of genome-wide genotyping of pooled DNA as a single-stage design for phenotype association. By individually genotyping the top associations identified by pooling, we obtained a 14- to 16-fold enrichment of SNPs nominally associated with the phenotype, but we likely missed the top true associations. In addition, we assessed whether genotyping pooled DNA can serve as an inexpensive screen as the second stage of a multi-stage design with a large number of samples by comparing the most cost-effective 3-stage designs with 80% power to detect common variants with genotypic relative risk of 1.1, with and without pooling. Given the current state of the specific technology we employed and the associated genotyping costs, we showed through simulation that a design involving pooling would be 1.07 times more expensive than a design without pooling. Thus, while a significant amount of information exists within the data from pooled DNA, our analysis does not support genotyping pooled DNA as a means to efficiently identify common variants contributing small effects to phenotypes of interest. While our conclusions were based on the specific technology and study design we employed, the approach presented here will be useful for evaluating the utility of other or future genome-wide genotyping platforms in pooled DNA studies.

Association between microdeletion and microduplication at 16p11.2 and autism.

BACKGROUND: Autism spectrum disorder is a heritable developmental disorder in which chromosomal abnormalities are thought to play a role. METHODS: As a first component of a genomewide association study of families from the Autism Genetic Resource Exchange (AGRE), we used two novel algorithms to search for recurrent copy-number variations in genotype data from 751 multiplex families with autism. Specific recurrent de novo events were further evaluated in clinical-testing data from Children's Hospital Boston and in a large population study in Iceland. RESULTS: Among the AGRE families, we observed five instances of a de novo deletion of 593 kb on chromosome 16p11.2. Using comparative genomic hybridization, we observed the identical deletion in 5 of 512 children referred to Children's Hospital Boston for developmental delay, mental retardation, or suspected autism spectrum disorder, as well as in 3 of 299 persons with autism in an Icelandic population; the deletion was also carried by 2 of 18,834 unscreened Icelandic control subjects. The reciprocal duplication of this region occurred in 7 affected persons in AGRE families and 4 of the 512 children from Children's Hospital Boston. The duplication also appeared to be a high-penetrance risk factor. CONCLUSIONS: We have identified a novel, recurrent microdeletion and a reciprocal microduplication that carry substantial susceptibility to autism and appear to account for approximately 1% of cases. We did not identify other regions with similar aggregations of large de novo mutations.

Patient-derived glioblastoma cultures as a tool for small-molecule drug discovery.

There is a compelling need for new therapeutic strategies for glioblastoma multiforme (GBM). Preclinical target and therapeutic discovery for GBMs is primarily conducted using cell lines grown in serum-containing media, such as U-87 MG, which do not reflect the gene expression profiles of tumors found in GBM patients. To address this lack of representative models, we sought to develop a panel of patient-derived GBM models and characterize their genomic features, using RNA sequencing (RNA-seq) and growth characteristics, both when grown as neurospheres in culture, and grown orthotopically as xenografts in mice. When we compared these with commonly used GBM cell lines in the Cancer Cell Line Encyclopedia (CCLE), we found these patient-derived models to have greater diversity in gene expression and to better correspond to GBMs directly sequenced from patient tumor samples. We also evaluated the potential of these models for targeted therapy, by using the genomic characterization to identify small molecules that inhibit the growth of distinct subsets of GBMs, paving the way for precision medicines for GBM.

A Proof of Concept for Biomarker-Guided Targeted Therapy against Ovarian Cancer Based on Patient-Derived Tumor Xenografts.

Advanced ovarian cancers are a leading cause of cancer-related death in women and are currently treated with surgery and chemotherapy. This standard of care is often temporarily successful but exhibits a high rate of relapse, after which, treatment options are few. Here we investigate whether biomarker-guided use of multiple targeted therapies, including small molecules and antibody-drug conjugates, is a viable alternative. A panel of patient-derived ovarian cancer xenografts (PDX), similar in genetics and chemotherapy responsiveness to human tumors, was exposed to 21 monotherapies and combination therapies. Three monotherapies and one combination were found to be active in different subsets of PDX. Analysis of gene expression data identified biomarkers associated with responsiveness to each of the three targeted therapies, none of which directly inhibits an oncogenic driver. While no single treatment had as high a response rate as chemotherapy, nearly 90% of PDXs were eligible for and responded to at least one biomarker-guided treatment, including tumors resistant to standard chemotherapy. The distribution of biomarker positivity in The Cancer Genome Atlas data suggests the potential for a similar precision approach in human patients. SIGNIFICANCE: This study exploits a panel of patient-derived xenografts to demonstrate that most ovarian tumors can be matched to effective biomarker-guided treatments.

Tumor-derived IFN triggers chronic pathway agonism and sensitivity to ADAR loss.

Interferons (IFNs) are cytokines that play a critical role in limiting infectious and malignant diseases 1-4 . Emerging data suggest that the strength and duration of IFN signaling can differentially impact cancer therapies, including immune checkpoint blockade 5-7 . Here, we characterize the output of IFN signaling, specifically IFN-stimulated gene (ISG) signatures, in primary tumors from The Cancer Genome Atlas. While immune infiltration correlates with the ISG signature in some primary tumors, the existence of ISG signature-positive tumors without evident infiltration of IFN-producing immune cells suggests that cancer cells per se can be a source of IFN production. Consistent with this hypothesis, analysis of patient-derived tumor xenografts propagated in immune-deficient mice shows evidence of ISG-positive tumors that correlates with expression of human type I and III IFNs derived from the cancer cells. Mechanistic studies using cell line models from the Cancer Cell Line Encyclopedia that harbor ISG signatures demonstrate that this is a by-product of a STING-dependent pathway resulting in chronic tumor-derived IFN production. This imposes a transcriptional state on the tumor, poising it to respond to the aberrant accumulation of double-stranded RNA (dsRNA) due to increased sensor levels (MDA5, RIG-I and PKR). By interrogating our functional short-hairpin RNA screen dataset across 398 cancer cell lines, we show that this ISG transcriptional state creates a novel genetic vulnerability. ISG signature-positive cancer cells are sensitive to the loss of ADAR, a dsRNA-editing enzyme that is also an ISG. A genome-wide CRISPR genetic suppressor screen reveals that the entire type I IFN pathway and the dsRNA-activated kinase, PKR, are required for the lethality induced by ADAR depletion. Therefore, tumor-derived IFN resulting in chronic signaling creates a cellular state primed to respond to dsRNA accumulation, rendering ISG-positive tumors susceptible to ADAR loss.

Targeting FGFR overcomes EMT-mediated resistance in EGFR mutant non-small cell lung cancer.

Evolved resistance to tyrosine kinase inhibitor (TKI)-targeted therapies remains a major clinical challenge. In epidermal growth factor receptor (EGFR) mutant non-small-cell lung cancer (NSCLC), failure of EGFR TKIs can result from both genetic and epigenetic mechanisms of acquired drug resistance. Widespread reports of histologic and gene expression changes consistent with an epithelial-to-mesenchymal transition (EMT) have been associated with initially surviving drug-tolerant persister cells, which can seed bona fide genetic mechanisms of resistance to EGFR TKIs. While therapeutic approaches targeting fully resistant cells, such as those harboring an EGFRT790M mutation, have been developed, a clinical strategy for preventing the emergence of persister cells remains elusive. Using mesenchymal cell lines derived from biopsies of patients who progressed on EGFR TKI as surrogates for persister populations, we performed whole-genome CRISPR screening and identified fibroblast growth factor receptor 1 (FGFR1) as the top target promoting survival of mesenchymal EGFR mutant cancers. Although numerous previous reports of FGFR signaling contributing to EGFR TKI resistance in vitro exist, the data have not yet been sufficiently compelling to instigate a clinical trial testing this hypothesis, nor has the role of FGFR in promoting the survival of persister cells been elucidated. In this study, we find that combining EGFR and FGFR inhibitors inhibited the survival and expansion of EGFR mutant drug-tolerant cells over long time periods, preventing the development of fully resistant cancers in multiple vitro models and in vivo. These results suggest that dual EGFR and FGFR blockade may be a promising clinical strategy for both preventing and overcoming EMT-associated acquired drug resistance and provide motivation for the clinical study of combined EGFR and FGFR inhibition in EGFR-mutated NSCLCs.

Antitumor Properties of RAF709, a Highly Selective and Potent Inhibitor of RAF Kinase Dimers, in Tumors Driven by Mutant RAS or BRAF.

Resistance to the RAF inhibitor vemurafenib arises commonly in melanomas driven by the activated BRAF oncogene. Here, we report antitumor properties of RAF709, a novel ATP-competitive kinase inhibitor with high potency and selectivity against RAF kinases. RAF709 exhibited a mode of RAF inhibition distinct from RAF monomer inhibitors such as vemurafenib, showing equal activity against both RAF monomers and dimers. As a result, RAF709 inhibited MAPK signaling activity in tumor models harboring either BRAFV600 alterations or mutant N- and KRAS-driven signaling, with minimal paradoxical activation of wild-type RAF. In cell lines and murine xenograft models, RAF709 demonstrated selective antitumor activity in tumor cells harboring BRAF or RAS mutations compared with cells with wild-type BRAF and RAS genes. RAF709 demonstrated a direct pharmacokinetic/pharmacodynamic relationship in in vivo tumor models harboring KRAS mutation. Furthermore, RAF709 elicited regression of primary human tumor-derived xenograft models with BRAF, NRAS, or KRAS mutations with excellent tolerability. Our results support further development of inhibitors like RAF709, which represents a next-generation RAF inhibitor with unique biochemical and cellular properties that enables antitumor activities in RAS-mutant tumors.Significance: In an effort to develop RAF inhibitors with the appropriate pharmacological properties to treat RAS mutant tumors, RAF709, a compound with potency, selectivity, and in vivo properties, was developed that will allow preclinical therapeutic hypothesis testing, but also provide an excellent probe to further unravel the complexities of RAF kinase signaling. Cancer Res; 78(6); 1537-48. ©2018 AACR.