Publications

2016

Tirosh I, Venteicher AS, Hebert C, Escalante LE, Patel AP, Yizhak K, Fisher JM, Rodman C, Mount C, Filbin MG, et al. Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma. Nature. 2016;539:309–313.
Although human tumours are shaped by the genetic evolution of cancer cells, evidence also suggests that they display hierarchies related to developmental pathways and epigenetic programs in which cancer stem cells (CSCs) can drive tumour growth and give rise to differentiated progeny. Yet, unbiased evidence for CSCs in solid human malignancies remains elusive. Here we profile 4,347 single cells from six IDH1 or IDH2 mutant human oligodendrogliomas by RNA sequencing (RNA-seq) and reconstruct their developmental programs from genome-wide expression signatures. We infer that most cancer cells are differentiated along two specialized glial programs, whereas a rare subpopulation of cells is undifferentiated and associated with a neural stem cell expression program. Cells with expression signatures for proliferation are highly enriched in this rare subpopulation, consistent with a model in which CSCs are primarily responsible for fuelling the growth of oligodendroglioma in humans. Analysis of copy number variation (CNV) shows that distinct CNV sub-clones within tumours display similar cellular hierarchies, suggesting that the architecture of oligodendroglioma is primarily dictated by developmental programs. Subclonal point mutation analysis supports a similar model, although a full phylogenetic tree would be required to definitively determine the effect of genetic evolution on the inferred hierarchies. Our single-cell analyses provide insight into the cellular architecture of oligodendrogliomas at single-cell resolution and support the cancer stem cell model, with substantial implications for disease management.
Aguirre AJ, Meyers RM, Weir BA, Vazquez F, Zhang CZ, Ben-David U, Cook A, Ha G, Harrington WF, Doshi MB, et al. Genomic Copy Number Dictates a Gene-Independent Cell Response to CRISPR/Cas9 Targeting. Cancer Discov. 2016;6:914–29.
UNLABELLED: The CRISPR/Cas9 system enables genome editing and somatic cell genetic screens in mammalian cells. We performed genome-scale loss-of-function screens in 33 cancer cell lines to identify genes essential for proliferation/survival and found a strong correlation between increased gene copy number and decreased cell viability after genome editing. Within regions of copy-number gain, CRISPR/Cas9 targeting of both expressed and unexpressed genes, as well as intergenic loci, led to significantly decreased cell proliferation through induction of a G2 cell-cycle arrest. By examining single-guide RNAs that map to multiple genomic sites, we found that this cell response to CRISPR/Cas9 editing correlated strongly with the number of target loci. These observations indicate that genome targeting by CRISPR/Cas9 elicits a gene-independent antiproliferative cell response. This effect has important practical implications for the interpretation of CRISPR/Cas9 screening data and confounds the use of this technology for the identification of essential genes in amplified regions. SIGNIFICANCE: We found that the number of CRISPR/Cas9-induced DNA breaks dictates a gene-independent antiproliferative response in cells. These observations have practical implications for using CRISPR/Cas9 to interrogate cancer gene function and illustrate that cancer cells are highly sensitive to site-specific DNA damage, which may provide a path to novel therapeutic strategies. Cancer Discov; 6(8); 914-29. (c)2016 AACR.See related commentary by Sheel and Xue, p. 824See related article by Munoz et al., p. 900This article is highlighted in the In This Issue feature, p. 803.
Ben-David U, Ha G, Khadka P, Jin X, Wong B, Franke L, Golub TR. The landscape of chromosomal aberrations in breast cancer mouse models reveals driver-specific routes to tumorigenesis. Nat Commun. 2016;7:12160.
Aneuploidy and copy-number alterations (CNAs) are a hallmark of human cancer. Although genetically engineered mouse models (GEMMs) are commonly used to model human cancer, their chromosomal landscapes remain underexplored. Here we use gene expression profiles to infer CNAs in 3,108 samples from 45 mouse models, providing the first comprehensive catalogue of chromosomal aberrations in cancer GEMMs. Mining this resource, we find that most chromosomal aberrations accumulate late during breast tumorigenesis, and observe marked differences in CNA prevalence between mouse mammary tumours initiated with distinct drivers. Some aberrations are recurrent and unique to specific GEMMs, suggesting distinct driver-dependent routes to tumorigenesis. Synteny-based comparison of mouse and human tumours narrows critical regions in CNAs, thereby identifying candidate driver genes. We experimentally validate that loss of Stratifin (SFN) promotes HER2-induced tumorigenesis in human cells. These results demonstrate the power of GEMM CNA analysis to inform the pathogenesis of human cancer.
Kim E, Ilic N, Shrestha Y, Zou L, Kamburov A, Zhu C, Yang X, Lubonja R, Tran N, Nguyen C, et al. Systematic Functional Interrogation of Rare Cancer Variants Identifies Oncogenic Alleles. Cancer Discov. 2016;6:714–26.
UNLABELLED: Cancer genome characterization efforts now provide an initial view of the somatic alterations in primary tumors. However, most point mutations occur at low frequency, and the function of these alleles remains undefined. We have developed a scalable systematic approach to interrogate the function of cancer-associated gene variants. We subjected 474 mutant alleles curated from 5,338 tumors to pooled in vivo tumor formation assays and gene expression profiling. We identified 12 transforming alleles, including two in genes (PIK3CB, POT1) that have not been shown to be tumorigenic. One rare KRAS allele, D33E, displayed tumorigenicity and constitutive activation of known RAS effector pathways. By comparing gene expression changes induced upon expression of wild-type and mutant alleles, we inferred the activity of specific alleles. Because alleles found to be mutated only once in 5,338 tumors rendered cells tumorigenic, these observations underscore the value of integrating genomic information with functional studies. SIGNIFICANCE: Experimentally inferring the functional status of cancer-associated mutations facilitates the interpretation of genomic information in cancer. Pooled in vivo screen and gene expression profiling identified functional variants and demonstrated that expression of rare variants induced tumorigenesis. Variant phenotyping through functional studies will facilitate defining key somatic events in cancer. Cancer Discov; 6(7); 714-26. (c)2016 AACR.See related commentary by Cho and Collisson, p. 694This article is highlighted in the In This Issue feature, p. 681.
Berger AH, Brooks AN, Wu X, Shrestha Y, Chouinard C, Piccioni F, Bagul M, Kamburov A, Imielinski M, Hogstrom L, et al. High-throughput Phenotyping of Lung Cancer Somatic Mutations. Cancer Cell. 2016;30:214–228.
Recent genome sequencing efforts have identified millions of somatic mutations in cancer. However, the functional impact of most variants is poorly understood. Here we characterize 194 somatic mutations identified in primary lung adenocarcinomas. We present an expression-based variant-impact phenotyping (eVIP) method that uses gene expression changes to distinguish impactful from neutral somatic mutations. eVIP identified 69% of mutations analyzed as impactful and 31% as functionally neutral. A subset of the impactful mutations induces xenograft tumor formation in mice and/or confers resistance to cellular EGFR inhibition. Among these impactful variants are rare somatic, clinically actionable variants including EGFR S645C, ARAF S214C and S214F, ERBB2 S418T, and multiple BRAF variants, demonstrating that rare mutations can be functionally important in cancer.
Chapuy B, Roemer MG, Stewart C, Tan Y, Abo RP, Zhang L, Dunford AJ, Meredith DM, Thorner AR, Jordanova ES, et al. Targetable genetic features of primary testicular and primary central nervous system lymphomas. Blood. 2016;127:869–81.
Primary central nervous system lymphomas (PCNSLs) and primary testicular lymphomas (PTLs) are extranodal large B-cell lymphomas (LBCLs) with inferior responses to current empiric treatment regimens. To identify targetable genetic features of PCNSL and PTL, we characterized their recurrent somatic mutations, chromosomal rearrangements, copy number alterations (CNAs), and associated driver genes, and compared these comprehensive genetic signatures to those of diffuse LBCL and primary mediastinal large B-cell lymphoma (PMBL). These studies identify unique combinations of genetic alterations in discrete LBCL subtypes and subtype-selective bases for targeted therapy. PCNSLs and PTLs frequently exhibit genomic instability, and near-uniform, often biallelic, CDKN2A loss with rare TP53 mutations. PCNSLs and PTLs also use multiple genetic mechanisms to target key genes and pathways and exhibit near-uniform oncogenic Toll-like receptor signaling as a result of MYD88 mutation and/or NFKBIZ amplification, frequent concurrent B-cell receptor pathway activation, and deregulation of BCL6. Of great interest, PCNSLs and PTLs also have frequent 9p24.1/PD-L1/PD-L2 CNAs and additional translocations of these loci, structural bases of immune evasion that are shared with PMBL.
Lohr JG, Kim S, Gould J, Knoechel B, Drier Y, Cotton MJ, Gray D, Birrer N, Wong B, Ha G, et al. Genetic interrogation of circulating multiple myeloma cells at single-cell resolution. Sci Transl Med. 2016;8:363ra147.
Multiple myeloma (MM) remains an incurable disease, with a treatment-refractory state eventually developing in all patients. Constant clonal evolution and genetic heterogeneity of MM are a likely explanation for the emergence of drug-resistant disease. Monitoring of MM genomic evolution on therapy by serial bone marrow biopsy is unfortunately impractical because it involves an invasive and painful procedure. We describe how noninvasive and highly sensitive isolation and characterization of circulating tumor cells (CTCs) from peripheral blood at single-cell resolution recapitulate MM in the bone marrow. We demonstrate that CTCs provide the same genetic information as bone marrow MM cells and even reveal mutations with greater sensitivity than bone marrow biopsies in some cases. Single CTC RNA sequencing enables classification of MM and quantitative assessment of genes that are relevant for prognosis. We propose that the genomic characterization of CTCs should be included in clinical trials to follow the emergence of resistant subclones after MM therapy.
Kryukov GV, Wilson FH, Ruth JR, Paulk J, Tsherniak A, Marlow SE, Vazquez F, Weir BA, Fitzgerald ME, Tanaka M, et al. MTAP deletion confers enhanced dependency on the PRMT5 arginine methyltransferase in cancer cells. Science. 2016;351:1214–8.
The discovery of cancer dependencies has the potential to inform therapeutic strategies and to identify putative drug targets. Integrating data from comprehensive genomic profiling of cancer cell lines and from functional characterization of cancer cell dependencies, we discovered that loss of the enzyme methylthioadenosine phosphorylase (MTAP) confers a selective dependence on protein arginine methyltransferase 5 (PRMT5) and its binding partner WDR77. MTAP is frequently lost due to its proximity to the commonly deleted tumor suppressor gene, CDKN2A. We observed increased intracellular concentrations of methylthioadenosine (MTA, the metabolite cleaved by MTAP) in cells harboring MTAP deletions. Furthermore, MTA specifically inhibited PRMT5 enzymatic activity. Administration of either MTA or a small-molecule PRMT5 inhibitor showed a modest preferential impairment of cell viability for MTAP-null cancer cell lines compared with isogenic MTAP-expressing counterparts. Together, our findings reveal PRMT5 as a potential vulnerability across multiple cancer lineages augmented by a common "passenger" genomic alteration.