Showing papers by "Gad Getz published in 2016"

Analysis of protein-coding genetic variation in 60,706 humans

Abstract: Large-scale reference data sets of human genetic variation are critical for the medical and functional interpretation of DNA sequence changes. Here we describe the aggregation and analysis of high-quality exome (protein-coding region) DNA sequence data for 60,706 individuals of diverse ancestries generated as part of the Exome Aggregation Consortium (ExAC). This catalogue of human genetic diversity contains an average of one variant every eight bases of the exome, and provides direct evidence for the presence of widespread mutational recurrence. We have used this catalogue to calculate objective metrics of pathogenicity for sequence variants, and to identify genes subject to strong selection against various classes of mutation; identifying 3,230 genes with near-complete depletion of predicted protein-truncating variants, with 72% of these genes having no currently established human disease phenotype. Finally, we demonstrate that these data can be used for the efficient filtering of candidate disease-causing variants, and for the discovery of human 'knockout' variants in protein-coding genes.

Distinct patterns of somatic genome alterations in lung adenocarcinomas and squamous cell carcinomas

Abstract: Matthew Meyerson, Ramaswamy Govindan and colleagues examine the exome sequences and copy number profiles of 660 lung adenocarcinoma and 484 lung squamous cell carcinoma tumors. They identify novel significantly mutated genes and amplification peaks and find that around half of the tumors have at least five predicted neoepitopes.

Molecular Mechanisms of Resistance to First- and Second-Generation ALK Inhibitors in ALK-Rearranged Lung Cancer

Abstract: Advanced, anaplastic lymphoma kinase (ALK)-positive lung cancer is currently treated with the first-generation ALK inhibitor crizotinib followed by more potent, second-generation ALK inhibitors (e.g., ceritinib, alectinib) upon progression. Second-generation inhibitors are generally effective even in the absence of crizotinib-resistant ALK mutations, likely reflecting incomplete inhibition of ALK by crizotinib in many cases. Herein, we analyzed 103 repeat biopsies from ALK-positive patients progressing on various ALK inhibitors. We find that each ALK inhibitor is associated with a distinct spectrum of ALK resistance mutations and that the frequency of one mutation - ALK G1202R - increases significantly after treatment with second-generation agents. To investigate strategies to overcome resistance to second-generation ALK inhibitors, we examine the activity of the third-generation ALK inhibitor lorlatinib in a series of ceritinib-resistant, patient-derived cell lines, and observe that the presence of ALK resistance mutations is highly predictive for sensitivity to lorlatinib, whereas those cell lines without ALK mutations are resistant.

Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma

Abstract: 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.

Tumor cells can follow distinct evolutionary paths to become resistant to epidermal growth factor receptor inhibition

Abstract: Although mechanisms of acquired resistance of epidermal growth factor receptor (EGFR)-mutant non-small-cell lung cancers to EGFR inhibitors have been identified, little is known about how resistant clones evolve during drug therapy. Here we observe that acquired resistance caused by the EGFR(T790M) gatekeeper mutation can occur either by selection of pre-existing EGFR(T790M)-positive clones or via genetic evolution of initially EGFR(T790M)-negative drug-tolerant cells. The path to resistance impacts the biology of the resistant clone, as those that evolved from drug-tolerant cells had a diminished apoptotic response to third-generation EGFR inhibitors that target EGFR(T790M); treatment with navitoclax, an inhibitor of the anti-apoptotic factors BCL-xL and BCL-2 restored sensitivity. We corroborated these findings using cultures derived directly from EGFR inhibitor-resistant patient tumors. These findings provide evidence that clinically relevant drug-resistant cancer cells can both pre-exist and evolve from drug-tolerant cells, and they point to therapeutic opportunities to prevent or overcome resistance in the clinic.

Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma

Abstract: 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.

Resensitization to Crizotinib by the Lorlatinib Alk Resistance Mutation L1198F.

Abstract: In a patient who had metastatic anaplastic lymphoma kinase (ALK)-rearranged lung cancer, resistance to crizotinib developed because of a mutation in the ALK kinase domain. This mutation is predicted to result in a substitution of cysteine by tyrosine at amino acid residue 1156 (C1156Y). Her tumor did not respond to a second-generation ALK inhibitor, but it did respond to lorlatinib (PF-06463922), a third-generation inhibitor. When her tumor relapsed, sequencing of the resistant tumor revealed an ALK L1198F mutation in addition to the C1156Y mutation. The L1198F substitution confers resistance to lorlatinib through steric interference with drug binding. However, L1198F paradoxically enhances binding to crizotinib, negating the effect of C1156Y and resensitizing resistant cancers to crizotinib. The patient received crizotinib again, and her cancer-related symptoms and liver failure resolved. (Funded by Pfizer and others; ClinicalTrials.gov number, NCT01970865.).

Somatic ERCC2 mutations are associated with a distinct genomic signature in urothelial tumors.

Abstract: Alterations in DNA repair pathways are common in tumors and can result in characteristic mutational signatures; however, a specific mutational signature associated with somatic alterations in the nucleotide- excision repair (NER) pathway has not yet been identified. Here we examine the mutational processes operating in urothelial cancer, a tumor type in which the core NER gene ERCC2 is significantly mutated. Analysis of three independent urothelial tumor cohorts demonstrates a strong association between somatic ERCC2 mutations and the activity of a mutational signature characterized by a broad spectrum of base changes. In addition, we note an association between the activity of this signature and smoking that is independent of ERCC2 mutation status, providing genomic evidence of tobacco-related mutagenesis in urothelial cancer. Together, these analyses identify an NER-related mutational signature and highlight the related roles of DNA damage and subsequent DNA repair in shaping tumor mutational landscape.

Clonal evolution in patients with chronic lymphocytic leukaemia developing resistance to BTK inhibition

Abstract: Resistance to the Bruton's tyrosine kinase (BTK) inhibitor ibrutinib has been attributed solely to mutations in BTK and related pathway molecules Using whole-exome and deep-targeted sequencing, we dissect evolution of ibrutinib resistance in serial samples from five chronic lymphocytic leukaemia patients In two patients, we detect BTK-C481S mutation or multiple PLCG2 mutations The other three patients exhibit an expansion of clones harbouring del(8p) with additional driver mutations (EP300, MLL2 and EIF2A), with one patient developing trans-differentiation into CD19-negative histiocytic sarcoma Using droplet-microfluidic technology and growth kinetic analyses, we demonstrate the presence of ibrutinib-resistant subclones and estimate subclone size before treatment initiation Haploinsufficiency of TRAIL-R, a consequence of del(8p), results in TRAIL insensitivity, which may contribute to ibrutinib resistance These findings demonstrate that the ibrutinib therapy favours selection and expansion of rare subclones already present before ibrutinib treatment, and provide insight into the heterogeneity of genetic changes associated with ibrutinib resistance

The genomic landscape and evolution of endometrial carcinoma progression and abdominopelvic metastasis.

Abstract: Recent studies have detailed the genomic landscape of primary endometrial cancers, but the evolution of these cancers into metastases has not been characterized. We performed whole-exome sequencing of 98 tumor biopsies including complex atypical hyperplasias, primary tumors and paired abdominopelvic metastases to survey the evolutionary landscape of endometrial cancer. We expanded and reanalyzed The Cancer Genome Atlas (TCGA) data, identifying new recurrent alterations in primary tumors, including mutations in the estrogen receptor cofactor gene NRIP1 in 12% of patients. We found that likely driver events were present in both primary and metastatic tissue samples, with notable exceptions such as ARID1A mutations. Phylogenetic analyses indicated that the sampled metastases typically arose from a common ancestral subclone that was not detected in the primary tumor biopsy. These data demonstrate extensive genetic heterogeneity in endometrial cancers and relative homogeneity across metastatic sites.

Dramatic Response of BRAF V600E Mutant Papillary Craniopharyngioma to Targeted Therapy

Abstract: We recently reported that BRAF V600E is the principal oncogenic driver of papillary craniopharyngioma, a highly morbid intracranial tumor commonly refractory to treatment. Here, we describe our treatment of a man age 39 years with multiply recurrent BRAF V600E craniopharyngioma using dabrafenib (150mg, orally twice daily) and trametinib (2mg, orally twice daily). After 35 days of treatment, tumor volume was reduced by 85%. Mutations that commonly mediate resistance to MAPK pathway inhibition were not detected in a post-treatment sample by whole exome sequencing. A blood-based BRAF V600E assay detected circulating BRAF V600E in the patient's blood. Re-evaluation of the existing management paradigms for craniopharyngioma is warranted, as patient morbidity might be reduced by noninvasive mutation testing and neoadjuvant-targeted treatment.

Systematic functional interrogation of rare cancer variants identifies oncogenic alleles

Abstract: 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. ©2016 AACR. See related commentary by Cho and Collisson, [p. 694][1] . This article is highlighted in the In This Issue feature, [p. 681][2] [1]: /lookup/volpage/6/694?iss=7 [2]: /lookup/volpage/6/681?iss=7

Genomic evolution and chemoresistance in germ-cell tumours

Abstract: Germ-cell tumours (GCTs) are derived from germ cells and occur most frequently in the testes. GCTs are histologically heterogeneous and distinctly curable with chemotherapy. Gains of chromosome arm 12p and aneuploidy are nearly universal in GCTs, but specific somatic genomic features driving tumour initiation, chemosensitivity and progression are incompletely characterized. Here, using clinical whole-exome and transcriptome sequencing of precursor, primary (testicular and mediastinal) and chemoresistant metastatic human GCTs, we show that the primary somatic feature of GCTs is highly recurrent chromosome arm level amplifications and reciprocal deletions (reciprocal loss of heterozygosity), variations that are significantly enriched in GCTs compared to 19 other cancer types. These tumours also acquire KRAS mutations during the development from precursor to primary disease, and primary testicular GCTs (TGCTs) are uniformly wild type for TP53. In addition, by functional measurement of apoptotic signalling (BH3 profiling) of fresh tumour and adjacent tissue, we find that primary TGCTs have high mitochondrial priming that facilitates chemotherapy-induced apoptosis. Finally, by phylogenetic analysis of serial TGCTs that emerge with chemotherapy resistance, we show how TGCTs gain additional reciprocal loss of heterozygosity and that this is associated with loss of pluripotency markers (NANOG and POU5F1) in chemoresistant teratomas or transformed carcinomas. Our results demonstrate the distinct genomic features underlying the origins of this disease and associated with the chemosensitivity phenotype, as well as the rare progression to chemoresistance. These results identify the convergence of cancer genomics, mitochondrial priming and GCT evolution, and may provide insights into chemosensitivity and resistance in other cancers.

Whole Exome Sequencing Identifies TSC1/TSC2 Biallelic Loss as the Primary and Sufficient Driver Event for Renal Angiomyolipoma Development

Abstract: Renal angiomyolipoma is a kidney tumor in the perivascular epithelioid (PEComa) family that is common in patients with Tuberous Sclerosis Complex (TSC) and Lymphangioleiomyomatosis (LAM) but occurs rarely sporadically. Though histologically benign, renal angiomyolipoma can cause life-threatening hemorrhage and kidney failure. Both angiomyolipoma and LAM have mutations in TSC2 or TSC1. However, the frequency and contribution of other somatic events in tumor development is unknown. We performed whole exome sequencing in 32 resected tumor samples (n = 30 angiomyolipoma, n = 2 LAM) from 15 subjects, including three with TSC. Two germline and 22 somatic inactivating mutations in TSC2 were identified, and one germline TSC1 mutation. Twenty of 32 (62%) samples showed copy neutral LOH (CN-LOH) in TSC2 or TSC1 with at least 8 different LOH regions, and 30 of 32 (94%) had biallelic loss of either TSC2 or TSC1. Whole exome sequencing identified a median of 4 somatic non-synonymous coding region mutations (other than in TSC2/TSC1), a mutation rate lower than nearly all other cancer types. Three genes with mutations were known cancer associated genes (BAP1, ARHGAP35 and SPEN), but they were mutated in a single sample each, and were missense variants with uncertain functional effects. Analysis of sixteen angiomyolipomas from a TSC subject showed both second hit point mutations and CN-LOH in TSC2, many of which were distinct, indicating that they were of independent clonal origin. However, three tumors had two shared mutations in addition to private somatic mutations, suggesting a branching evolutionary pattern of tumor development following initiating loss of TSC2. Our results indicate that TSC2 and less commonly TSC1 alterations are the primary essential driver event in angiomyolipoma/LAM, whereas other somatic mutations are rare and likely do not contribute to tumor development.

Local genetic effects on gene expression across 44 human tissues

Abstract: Expression quantitative trait locus (eQTL) mapping provides a powerful means to identify functional variants influencing gene expression and disease pathogenesis. We report the identification of cis-eQTLs from 7,051 post-mortem samples representing 44 tissues and 449 individuals as part of the Genotype-Tissue Expression (GTEx) project. We find a cis-eQTL for 88% of all annotated protein-coding genes, with one-third having multiple independent effects. We identify numerous tissue-specific cis-eQTLs, highlighting the unique functional impact of regulatory variation in diverse tissues. By integrating large-scale functional genomics data and state-of-the-art fine-mapping algorithms, we identify multiple features predictive of tissue-specific and shared regulatory effects. We improve estimates of cis-eQTL sharing and effect sizes using allele specific expression across tissues. Finally, we demonstrate the utility of this large compendium of cis-eQTLs for understanding the tissue-specific etiology of complex traits, including coronary artery disease. The GTEx project provides an exceptional resource that has improved our understanding of gene regulation across tiss

Characterizing genomic alterations in cancer by complementary functional associations

Abstract: Systematic efforts to sequence the cancer genome have identified large numbers of mutations and copy number alterations in human cancers. However, elucidating the functional consequences of these variants, and their interactions to drive or maintain oncogenic states, remains a challenge in cancer research. We developed REVEALER, a computational method that identifies combinations of mutually exclusive genomic alterations correlated with functional phenotypes, such as the activation or gene dependency of oncogenic pathways or sensitivity to a drug treatment. We used REVEALER to uncover complementary genomic alterations associated with the transcriptional activation of β-catenin and NRF2, MEK-inhibitor sensitivity, and KRAS dependency. REVEALER successfully identified both known and new associations, demonstrating the power of combining functional profiles with extensive characterization of genomic alterations in cancer genomes.

Rapid and efficient analysis of 20,000 RNA-seq samples with Toil

Abstract: Toil is portable, open-source workflow software that supports contemporary workflow definition languages and can be used to securely and reproducibly run scientific workflows efficiently at large-scale. To demonstrate Toil, we processed over 20,000 RNA-seq samples to create a consistent meta-analysis of five datasets free of computational batch effects that we make freely available. Nearly all the samples were analysed in under four days using a commercial cloud cluster of 32,000 preemptable cores.

Abstract 4368: High-throughput phenotyping of lung cancer somatic mutations

Abstract: Recent cancer genome sequencing and analysis has identified millions of somatic mutations in cancer. However, the functional impact of most variants is poorly understood, limiting the use of this genetic knowledge for clinical decision-making. Here we describe a new high-throughput approach, expression-based variant impact phenotyping (eVIP), which uses gene expression changes to infer somatic mutation impact. We generated a lentiviral expression library representing 53 genes and 194 somatic mutations identified in primary lung adenocarcinomas. Next, we introduced this library into A549 lung adenocarcinoma cells and 96 hours later performed gene expression profiling using Luminex-based L1000 profiling. We built a computational pipeline, eVIP, to compare mutant and wild-type expression signatures to infer whether variants were gain-of-function, change-of-function, loss-of-function, or neutral. Overall, eVIP identified 69% of mutations as impactful whereas 31% appeared functionally neutral. A very high rate, 92%, of missense mutations in the KEAP1 and STK11 tumor suppressor genes were found to inactivate or diminish protein function. As a complementary approach, we assessed which mutations are epistatic to EGFR or capable of initiating xenograft tumor formation in vivo. A subset of the impactful mutations identified by eVIP could induce xenograft tumor formation in mice and/or confer resistance to cellular EGFR inhibition. Among these mutations were 20 rare or non-canonical somatic variants in clinically-actionable or -relevant oncogenes including EGFR S645C, ARAF S214C and S214F, ERBB2 S418T, and PIK3CA E600K. eVIP can, in principle, characterize any genetic variant, independent of prior knowledge of gene function. Further application of eVIP should significantly advance the pace of functional characterization of mutations identified from genome sequencing. Citation Format: Alice H. Berger, Angela N. Brooks, Xiaoyun Wu, Yashaswi Shrestha, Candace Chouinard, Federica Piccioni, Mukta Bagul, Atanas Kamburov, Marcin Imielinski, Larson Hogstrom, Cong Zhu, Xiaoping Yang, Sasha Pantel, Ryo Sakai, Nathan Kaplan, David Root, Rajiv Narayan, Ted Natoli, David Lahr, Itay Tirosh, Pablo Tamayo, Gad Getz, Bang Wong, John Doench, Aravind Subramanian, Todd R. Golub, Matthew Meyerson, Jesse S. Boehm. High-throughput phenotyping of lung cancer somatic mutations. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 4368.