Mapping the landscape of genetic dependencies in chordoma.

Identifying the spectrum of genes required for cancer cell survival can reveal essential cancer circuitry and therapeutic targets, but such a map remains incomplete for many cancer types. We apply genome-scale CRISPR-Cas9 loss-of-function screens to map the landscape of selectively essential genes in chordoma, a bone cancer with few validated targets. This approach confirms a known chordoma dependency, TBXT (T; brachyury), and identifies a range of additional dependencies, including PTPN11, ADAR, PRKRA, LUC7L2, SRRM2, SLC2A1, SLC7A5, FANCM, and THAP1. CDK6, SOX9, and EGFR, genes previously implicated in chordoma biology, are also recovered. We find genomic and transcriptomic features that predict specific dependencies, including interferon-stimulated gene expression, which correlates with ADAR dependence and is elevated in chordoma. Validating the therapeutic relevance of dependencies, small-molecule inhibitors of SHP2, encoded by PTPN11, have potent preclinical efficacy against chordoma. Our results generate an emerging map of chordoma dependencies to enable biological and therapeutic hypotheses.

Targeted brachyury degradation disrupts a highly specific autoregulatory program controlling chordoma cell identity.

Chordomas are rare spinal tumors addicted to expression of the developmental transcription factor brachyury. In chordomas, brachyury is super-enhancer associated and preferentially downregulated by pharmacologic transcriptional CDK inhibition, leading to cell death. To understand the underlying basis of this sensitivity, we dissect the brachyury transcription regulatory network and compare the consequences of brachyury degradation with transcriptional CDK inhibition. Brachyury defines the chordoma super-enhancer landscape and autoregulates through binding its super-enhancer, and its locus forms a transcriptional condensate. Transcriptional CDK inhibition and brachyury degradation disrupt brachyury autoregulation, leading to loss of its transcriptional condensate and transcriptional program. Compared with transcriptional CDK inhibition, which globally downregulates transcription, leading to cell death, brachyury degradation is much more selective, inducing senescence and sensitizing cells to anti-apoptotic inhibition. These data suggest that brachyury downregulation is a core tenet of transcriptional CDK inhibition and motivates developing strategies to target brachyury and its autoregulatory feedback loop.

Small-molecule targeting of brachyury transcription factor addiction in chordoma.

Chordoma is a primary bone cancer with no approved therapy1. The identification of therapeutic targets in this disease has been challenging due to the infrequent occurrence of clinically actionable somatic mutations in chordoma tumors2,3. Here we describe the discovery of therapeutically targetable chordoma dependencies via genome-scale CRISPR-Cas9 screening and focused small-molecule sensitivity profiling. These systematic approaches reveal that the developmental transcription factor T (brachyury; TBXT) is the top selectively essential gene in chordoma, and that transcriptional cyclin-dependent kinase (CDK) inhibitors targeting CDK7/12/13 and CDK9 potently suppress chordoma cell proliferation. In other cancer types, transcriptional CDK inhibitors have been observed to downregulate highly expressed, enhancer-associated oncogenic transcription factors4,5. In chordoma, we find that T is associated with a 1.5-Mb region containing 'super-enhancers' and is the most highly expressed super-enhancer-associated transcription factor. Notably, transcriptional CDK inhibition leads to preferential and concentration-dependent downregulation of cellular brachyury protein levels in all models tested. In vivo, CDK7/12/13-inhibitor treatment substantially reduces tumor growth. Together, these data demonstrate small-molecule targeting of brachyury transcription factor addiction in chordoma, identify a mechanism of T gene regulation that underlies this therapeutic strategy, and provide a blueprint for applying systematic genetic and chemical screening approaches to discover vulnerabilities in genomically quiet cancers.

Emerging Opportunities for Target Discovery in Rare Cancers.

Rare cancers pose unique challenges to research due to their low incidence. Barriers include a scarcity of tissue and experimental models to enable basic research and insufficient patient accrual for clinical studies. Consequently, an understanding of the genetic and cellular features of many rare cancer types and their associated vulnerabilities has been lacking. However, new opportunities are emerging to facilitate discovery of therapeutic targets in rare cancers. Online platforms are allowing patients with rare cancers to organize on an unprecedented scale, tumor genome sequencing is now routinely performed in research and clinical settings, and the efficiency of patient-derived model generation has improved. New CRISPR/Cas9 and small-molecule libraries permit cancer dependency discovery in a rapid and systematic fashion. In parallel, large-scale studies of common cancers now provide reference datasets to help interpret rare cancer profiling data. Together, these advances motivate consideration of new research frameworks to accelerate rare cancer target discovery.

A functional landscape of resistance to ALK inhibition in lung cancer.

We conducted a large-scale functional genetic study to characterize mechanisms of resistance to ALK inhibition in ALK-dependent lung cancer cells. We identify members of known resistance pathways and additional putative resistance drivers. Among the latter were members of the P2Y purinergic receptor family of G-protein-coupled receptors (P2Y1, P2Y2, and P2Y6). P2Y receptors mediated resistance in part through a protein-kinase-C (PKC)-dependent mechanism. Moreover, PKC activation alone was sufficient to confer resistance to ALK inhibitors, whereas combined ALK and PKC inhibition restored sensitivity. We observed enrichment of gene signatures associated with several resistance drivers (including P2Y receptors) in crizotinib-resistant ALK-rearranged lung tumors compared to treatment-naive controls, supporting a role for these identified mechanisms in clinical ALK inhibitor resistance.

Genetic modifiers of EGFR dependence in non-small cell lung cancer.

Lung adenocarcinomas harboring activating mutations in the epidermal growth factor receptor (EGFR) represent a common molecular subset of non-small cell lung cancer (NSCLC) cases. EGFR mutations predict sensitivity to EGFR tyrosine kinase inhibitors (TKIs) and thus represent a dependency in NSCLCs harboring these alterations, but the genetic basis of EGFR dependence is not fully understood. Here, we applied an unbiased, ORF-based screen to identify genetic modifiers of EGFR dependence in EGFR-mutant NSCLC cells. This approach identified 18 kinase and kinase-related genes whose overexpression can substitute for EGFR in EGFR-dependent PC9 cells, and these genes include seven of nine Src family kinase genes, FGFR1, FGFR2, ITK, NTRK1, NTRK2, MOS, MST1R, and RAF1. A subset of these genes can complement loss of EGFR activity across multiple EGFR-dependent models. Unbiased gene-expression profiling of cells overexpressing EGFR bypass genes, together with targeted validation studies, reveals EGFR-independent activation of the MEK-ERK and phosphoinositide 3-kinase (PI3K)-AKT pathways. Combined inhibition of PI3K-mTOR and MEK restores EGFR dependence in cells expressing each of the 18 EGFR bypass genes. Together, these data uncover a broad spectrum of kinases capable of overcoming dependence on EGFR and underscore their convergence on the PI3K-AKT and MEK-ERK signaling axes in sustaining EGFR-independent survival.

A pan-cancer analysis of transcriptome changes associated with somatic mutations in U2AF1 reveals commonly altered splicing events.

Although recurrent somatic mutations in the splicing factor U2AF1 (also known as U2AF35) have been identified in multiple cancer types, the effects of these mutations on the cancer transcriptome have yet to be fully elucidated. Here, we identified splicing alterations associated with U2AF1 mutations across distinct cancers using DNA and RNA sequencing data from The Cancer Genome Atlas (TCGA). Using RNA-Seq data from 182 lung adenocarcinomas and 167 acute myeloid leukemias (AML), in which U2AF1 is somatically mutated in 3-4% of cases, we identified 131 and 369 splicing alterations, respectively, that were significantly associated with U2AF1 mutation. Of these, 30 splicing alterations were statistically significant in both lung adenocarcinoma and AML, including three genes in the Cancer Gene Census, CTNNB1, CHCHD7, and PICALM. Cell line experiments expressing U2AF1 S34F in HeLa cells and in 293T cells provide further support that these altered splicing events are caused by U2AF1 mutation. Consistent with the function of U2AF1 in 3' splice site recognition, we found that S34F/Y mutations cause preferences for CAG over UAG 3' splice site sequences. This report demonstrates consistent effects of U2AF1 mutation on splicing in distinct cancer cell types.

Mapping the hallmarks of lung adenocarcinoma with massively parallel sequencing.

Lung adenocarcinoma, the most common subtype of non-small cell lung cancer, is responsible for more than 500,000 deaths per year worldwide. Here, we report exome and genome sequences of 183 lung adenocarcinoma tumor/normal DNA pairs. These analyses revealed a mean exonic somatic mutation rate of 12.0 events/megabase and identified the majority of genes previously reported as significantly mutated in lung adenocarcinoma. In addition, we identified statistically recurrent somatic mutations in the splicing factor gene U2AF1 and truncating mutations affecting RBM10 and ARID1A. Analysis of nucleotide context-specific mutation signatures grouped the sample set into distinct clusters that correlated with smoking history and alterations of reported lung adenocarcinoma genes. Whole-genome sequence analysis revealed frequent structural rearrangements, including in-frame exonic alterations within EGFR and SIK2 kinases. The candidate genes identified in this study are attractive targets for biological characterization and therapeutic targeting of lung adenocarcinoma.

Inhibitor-sensitive FGFR1 amplification in human non-small cell lung cancer.

BACKGROUND: Squamous cell lung carcinomas account for approximately 25% of new lung carcinoma cases and 40,000 deaths per year in the United States. Although there are multiple genomically targeted therapies for lung adenocarcinoma, none has yet been reported in squamous cell lung carcinoma. METHODOLOGY/PRINCIPAL FINDINGS: Using SNP array analysis, we found that a region of chromosome segment 8p11-12 containing three genes-WHSC1L1, LETM2, and FGFR1-is amplified in 3% of lung adenocarcinomas and 21% of squamous cell lung carcinomas. Furthermore, we demonstrated that a non-small cell lung carcinoma cell line harboring focal amplification of FGFR1 is dependent on FGFR1 activity for cell growth, as treatment of this cell line either with FGFR1-specific shRNAs or with FGFR small molecule enzymatic inhibitors leads to cell growth inhibition. CONCLUSIONS/SIGNIFICANCE: These studies show that FGFR1 amplification is common in squamous cell lung cancer, and that FGFR1 may represent a promising therapeutic target in non-small cell lung cancer.

Small-molecule inducers of insulin expression in pancreatic alpha-cells.

High-content screening for small-molecule inducers of insulin expression identified the compound BRD7389, which caused alpha-cells to adopt several morphological and gene expression features of a beta-cell state. Assay-performance profile analysis suggests kinase inhibition as a mechanism of action, and we show that biochemical and cellular inhibition of the RSK kinase family by BRD7389 is likely related to its ability induce a beta-cell-like state. BRD7389 also increases the endocrine cell content and function of donor human pancreatic islets in culture.

Highly parallel identification of essential genes in cancer cells.

More complete knowledge of the molecular mechanisms underlying cancer will improve prevention, diagnosis and treatment. Efforts such as The Cancer Genome Atlas are systematically characterizing the structural basis of cancer, by identifying the genomic mutations associated with each cancer type. A powerful complementary approach is to systematically characterize the functional basis of cancer, by identifying the genes essential for growth and related phenotypes in different cancer cells. Such information would be particularly valuable for identifying potential drug targets. Here, we report the development of an efficient, robust approach to perform genome-scale pooled shRNA screens for both positive and negative selection and its application to systematically identify cell essential genes in 12 cancer cell lines. By integrating these functional data with comprehensive genetic analyses of primary human tumors, we identified known and putative oncogenes such as EGFR, KRAS, MYC, BCR-ABL, MYB, CRKL, and CDK4 that are essential for cancer cell proliferation and also altered in human cancers. We further used this approach to identify genes involved in the response of cancer cells to tumoricidal agents and found 4 genes required for the response of CML cells to imatinib treatment: PTPN1, NF1, SMARCB1, and SMARCE1, and 5 regulators of the response to FAS activation, FAS, FADD, CASP8, ARID1A and CBX1. Broad application of this highly parallel genetic screening strategy will not only facilitate the rapid identification of genes that drive the malignant state and its response to therapeutics but will also enable the discovery of genes that participate in any biological process.