Clonal evolution of glioblastoma under therapy

Glioblastoma (GBM) is the most common and aggressive primary brain tumor. To better understand how GBM evolves, we analyzed longitudinal genomic and transcriptomic data from 114 patients. The analysis shows a highly branched evolutionary pattern in which 63% of patients experience expression-based subtype changes. The branching pattern, together with estimates of evolutionary rate, suggests that relapse-associated clones typically existed years before diagnosis. Fifteen percent of tumors present hypermutation at relapse in highly expressed genes, with a clear mutational signature. We find that 11% of recurrence tumors harbor mutations in LTBP4, which encodes a protein binding to TGF-β. Silencing LTBP4 in GBM cells leads to suppression of TGF-β activity and decreased cell proliferation. In recurrent GBM with wild-type IDH1, high LTBP4 expression is associated with worse prognosis, highlighting the TGF-β pathway as a potential therapeutic target in GBM.

[1]  David T. W. Jones,et al.  Signatures of mutational processes in human cancer , 2013, Nature.

[2]  P. Armitage,et al.  The age distribution of cancer and a multi-stage theory of carcinogenesis , 1954, British Journal of Cancer.

[3]  C. Nordling A New Theory on the Cancer-inducing Mechanism , 1953, British Journal of Cancer.

[4]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[5]  Trevor Bedford,et al.  Overdispersion of the molecular clock: temporal variation of gene-specific substitution rates in Drosophila. , 2008, Molecular biology and evolution.

[6]  J. Baselga,et al.  TGF-β Receptor Inhibitors Target the CD44(high)/Id1(high) Glioma-Initiating Cell Population in Human Glioblastoma. , 2010, Cancer cell.

[7]  Luca Laurenti,et al.  Tumor evolutionary directed graphs and the history of chronic lymphocytic leukemia , 2014, eLife.

[8]  D. Mercola,et al.  Eradication of established intracranial rat gliomas by transforming growth factor beta antisense gene therapy. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[9]  Tracy T Batchelor,et al.  A hypermutation phenotype and somatic MSH6 mutations in recurrent human malignant gliomas after alkylator chemotherapy. , 2006, Cancer research.

[10]  T. Golub,et al.  The Genomic Landscape of Juvenile Myelomonocytic Leukemia , 2015, Nature Genetics.

[11]  D. Haussler,et al.  The Somatic Genomic Landscape of Glioblastoma , 2013, Cell.

[12]  A. Iavarone,et al.  Detection, Characterization, and Inhibition of FGFR–TACC Fusions in IDH Wild-type Glioma , 2015, Clinical Cancer Research.

[13]  Raul Rabadan,et al.  MutComFocal: an integrative approach to identifying recurrent and focal genomic alterations in tumor samples , 2013, BMC Systems Biology.

[14]  D. Bigner,et al.  EGF mutant receptor vIII as a molecular target in cancer therapy. , 2001, Endocrine-related cancer.

[15]  Tom H. Pringle,et al.  The human genome browser at UCSC. , 2002, Genome research.

[16]  J. Otte,et al.  Disruption of the gene encoding the latent transforming growth factor-beta binding protein 4 (LTBP-4) causes abnormal lung development, cardiomyopathy, and colorectal cancer. , 2002, Genes & development.

[17]  Shawn M. Gillespie,et al.  Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma , 2014, Science.

[18]  Raul Rabadan,et al.  SAVI: a statistical algorithm for variant frequency identification , 2013, BMC Systems Biology.

[19]  K. Miyazono,et al.  A role of the latent TGF‐beta 1‐binding protein in the assembly and secretion of TGF‐beta 1. , 1991, The EMBO journal.

[20]  A. Iavarone,et al.  The HECT-domain ubiquitin ligase Huwe1 controls neural differentiation and proliferation by destabilizing the N-Myc oncoprotein , 2008, Nature Cell Biology.

[21]  Jill S Barnholtz-Sloan,et al.  Whole-genome and multisector exome sequencing of primary and post-treatment glioblastoma reveals patterns of tumor evolution , 2015, Genome research.

[22]  G. Getz,et al.  GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers , 2011, Genome Biology.

[23]  A. McKenna,et al.  Absolute quantification of somatic DNA alterations in human cancer , 2012, Nature Biotechnology.

[24]  P. Nowell The clonal evolution of tumor cell populations. , 1976, Science.

[25]  V. Balasubramaniyan,et al.  TGF-β as a therapeutic target in high grade gliomas - promises and challenges. , 2013, Biochemical pharmacology.

[26]  Raul Rabadan,et al.  An information theoretic method to identify combinations of genomic alterations that promote glioblastoma. , 2015, Journal of molecular cell biology.

[27]  A. Verma MGMT Gene Silencing and Benefit From Temozolomide in Glioblastoma , 2006 .

[28]  Stephen Yip,et al.  MSH6 Mutations Arise in Glioblastomas during Temozolomide Therapy and Mediate Temozolomide Resistance , 2009, Clinical Cancer Research.

[29]  J. Uhm Comprehensive genomic characterization defines human glioblastoma genes and core pathways , 2009 .

[30]  Cole Trapnell,et al.  TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions , 2013, Genome Biology.

[31]  Robert J Dempsey,et al.  Science Times , 2022 .

[32]  Chibo Hong,et al.  DNA Methylation and Somatic Mutations Converge on the Cell Cycle and Define Similar Evolutionary Histories in Brain Tumors. , 2015, Cancer cell.

[33]  Raul Rabadan,et al.  The integrated landscape of driver genomic alterations in glioblastoma , 2013, Nature Genetics.

[34]  Chris Wiggins,et al.  Pegasus: a comprehensive annotation and prediction tool for detection of driver gene fusions in cancer , 2014, BMC Systems Biology.

[35]  T. D. Schneider,et al.  Sequence logos: a new way to display consensus sequences. , 1990, Nucleic acids research.

[36]  Martin J. van den Bent,et al.  Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. , 2005, The New England journal of medicine.

[37]  Satoru Miyano,et al.  Mutational landscape and clonal architecture in grade II and III gliomas , 2015, Nature Genetics.

[38]  B. Kamińska,et al.  TGF beta signaling and its role in glioma pathogenesis. , 2013, Advances in experimental medicine and biology.

[39]  K. Aldape,et al.  Mesenchymal high-grade glioma is maintained by the ID-RAP1 axis. , 2013, The Journal of clinical investigation.

[40]  Zoltan Szallasi,et al.  In silico prediction of tumor antigens derived from functional missense mutations of the cancer gene census , 2012, Oncoimmunology.

[41]  Peter Canoll,et al.  MRI-localized biopsies reveal subtype-specific differences in molecular and cellular composition at the margins of glioblastoma , 2014, Proceedings of the National Academy of Sciences.

[42]  Rebecca A Betensky,et al.  Loss of the mismatch repair protein MSH6 in human glioblastomas is associated with tumor progression during temozolomide treatment. , 2007, Clinical cancer research : an official journal of the American Association for Cancer Research.

[43]  Benjamin E. Gross,et al.  Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal , 2013, Science Signaling.

[44]  J. Massagué,et al.  TGFβ in Cancer , 2008, Cell.

[45]  Peter J. Campbell,et al.  Evolution of the cancer genome , 2012, Nature Reviews Genetics.

[46]  K. Hoang-Xuan,et al.  Primary brain tumours in adults , 2003, The Lancet.

[47]  M. Wigler,et al.  Circular binary segmentation for the analysis of array-based DNA copy number data. , 2004, Biostatistics.

[48]  Nordling Co A New Theory on the Cancer-inducing Mechanism , 1953 .

[49]  S. Gabriel,et al.  Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. , 2010, Cancer cell.

[50]  Trevor Bedford,et al.  Overdispersion of the Molecular Clock Varies Between Yeast, Drosophila and Mammals , 2008, Genetics.

[51]  J. Uhm,et al.  The transcriptional network for mesenchymal transformation of brain tumours , 2010 .

[52]  V. P. Collins,et al.  Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics , 2013, Proceedings of the National Academy of Sciences.

[53]  Andrew J. Blumberg,et al.  Moduli Spaces of Phylogenetic Trees Describing Tumor Evolutionary Patterns , 2014, Brain Informatics and Health.

[54]  Christopher A. Maher,et al.  ChimeraScan: a tool for identifying chimeric transcription in sequencing data , 2011, Bioinform..

[55]  Qi-En Wang,et al.  Review Article TGF-β signaling and its targeting for glioma treatment , 2015 .

[56]  C. Sander,et al.  Mutual exclusivity analysis identifies oncogenic network modules. , 2012, Genome research.

[57]  Cole Trapnell,et al.  Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. , 2010, Nature biotechnology.

[58]  Benjamin E. Gross,et al.  The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. , 2012, Cancer discovery.

[59]  E. E. Gresch Genetic Alterations During Colorectal-Tumor Development , 1989 .

[60]  In-Hee Lee,et al.  Spatiotemporal Evolution of the Primary Glioblastoma Genome. , 2015, Cancer cell.

[61]  A. Bouchard-Côté,et al.  PyClone: statistical inference of clonal population structure in cancer , 2014, Nature Methods.

[62]  B. Giusti,et al.  EXCAVATOR: detecting copy number variants from whole-exome sequencing data , 2013, Genome Biology.

[63]  John N. Weinstein,et al.  PRADA: pipeline for RNA sequencing data analysis , 2014, Bioinform..

[64]  D. Brat,et al.  Transforming Fusions of FGFR and TACC Genes in Human Glioblastoma , 2012, Science.