Whole-genome and multisector exome sequencing of primary and post-treatment glioblastoma reveals patterns of tumor evolution

Glioblastoma (GBM) is a prototypical heterogeneous brain tumor refractory to conventional therapy. A small residual population of cells escapes surgery and chemoradiation, resulting in a typically fatal tumor recurrence ∼ 7 mo after diagnosis. Understanding the molecular architecture of this residual population is critical for the development of successful therapies. We used whole-genome sequencing and whole-exome sequencing of multiple sectors from primary and paired recurrent GBM tumors to reconstruct the genomic profile of residual, therapy resistant tumor initiating cells. We found that genetic alteration of the p53 pathway is a primary molecular event predictive of a high number of subclonal mutations in glioblastoma. The genomic road leading to recurrence is highly idiosyncratic but can be broadly classified into linear recurrences that share extensive genetic similarity with the primary tumor and can be directly traced to one of its specific sectors, and divergent recurrences that share few genetic alterations with the primary tumor and originate from cells that branched off early during tumorigenesis. Our study provides mechanistic insights into how genetic alterations in primary tumors impact the ensuing evolution of tumor cells and the emergence of subclonal heterogeneity.

[1]  Obi L. Griffith,et al.  SciClone: Inferring Clonal Architecture and Tracking the Spatial and Temporal Patterns of Tumor Evolution , 2014, PLoS Comput. Biol..

[2]  J. C. Love,et al.  EGFR variant heterogeneity in glioblastoma resolved through single-nucleus sequencing. , 2014, Cancer discovery.

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

[4]  P. A. Futreal,et al.  Genomic architecture and evolution of clear cell renal cell carcinomas defined by multiregion sequencing , 2014, Nature Genetics.

[5]  Steven J. M. Jones,et al.  Mutational Analysis Reveals the Origin and Therapy-Driven Evolution of Recurrent Glioma , 2014, Science.

[6]  Noemi Andor,et al.  EXPANDS: expanding ploidy and allele frequency on nested subpopulations , 2013, Bioinform..

[7]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

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

[9]  V. Pant,et al.  The p53–Mdm2 feedback loop protects against DNA damage by inhibiting p53 activity but is dispensable for p53 stability, development, and longevity , 2013, Genes & development.

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

[11]  Ken Chen,et al.  A survey of intragenic breakpoints in glioblastoma identifies a distinct subset associated with poor survival. , 2013, Genes & development.

[12]  Krishnendu Chatterjee,et al.  Evolutionary dynamics of cancer in response to targeted combination therapy , 2013, eLife.

[13]  Carlos Caldas,et al.  The implications of clonal genome evolution for cancer medicine. , 2013, The New England journal of medicine.

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

[15]  A. Sivachenko,et al.  Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples , 2013, Nature Biotechnology.

[16]  Giovanni Parmigiani,et al.  Half or more of the somatic mutations in cancers of self-renewing tissues originate prior to tumor initiation , 2013, Proceedings of the National Academy of Sciences.

[17]  A. McKenna,et al.  Evolution and Impact of Subclonal Mutations in Chronic Lymphocytic Leukemia , 2012, Cell.

[18]  Steven A. Roberts,et al.  Mutational heterogeneity in cancer and the search for new cancer-associated genes , 2013 .

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

[20]  Irmtraud M. Meyer,et al.  The clonal and mutational evolution spectrum of primary triple-negative breast cancers , 2012, Nature.

[21]  A. Børresen-Dale,et al.  The Life History of 21 Breast Cancers , 2012, Cell.

[22]  M. Adams,et al.  Characterizing Mutational Heterogeneity in a Glioblastoma Patient with Double Recurrence , 2012, PloS one.

[23]  Lynda Chin,et al.  Emerging insights into the molecular and cellular basis of glioblastoma. , 2012, Genes & development.

[24]  P. A. Futreal,et al.  Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. , 2012, The New England journal of medicine.

[25]  Debyani Chakravarty,et al.  Intratumoral heterogeneity of receptor tyrosine kinases EGFR and PDGFRA amplification in glioblastoma defines subpopulations with distinct growth factor response , 2012, Proceedings of the National Academy of Sciences.

[26]  David T. W. Jones,et al.  Genome Sequencing of Pediatric Medulloblastoma Links Catastrophic DNA Rearrangements with TP53 Mutations , 2012, Cell.

[27]  Joshua F. McMichael,et al.  Clonal evolution in relapsed acute myeloid leukemia revealed by whole genome sequencing , 2011, Nature.

[28]  Rebecca A Betensky,et al.  Mosaic amplification of multiple receptor tyrosine kinase genes in glioblastoma. , 2011, Cancer cell.

[29]  Paul T. Spellman,et al.  Parent-specific copy number in paired tumor-normal studies using circular binary segmentation , 2011, Bioinform..

[30]  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.

[31]  J. Troge,et al.  Tumour evolution inferred by single-cell sequencing , 2011, Nature.

[32]  Raymond Sawaya,et al.  Immunologic escape after prolonged progression-free survival with epidermal growth factor receptor variant III peptide vaccination in patients with newly diagnosed glioblastoma. , 2010, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[33]  C. Perou,et al.  Allele-specific copy number analysis of tumors , 2010, Proceedings of the National Academy of Sciences.

[34]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[35]  H. Hakonarson,et al.  ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data , 2010, Nucleic acids research.

[36]  R. Wilson,et al.  Identification of a CpG island methylator phenotype that defines a distinct subgroup of glioma. , 2010, Cancer cell.

[37]  Erwin G. Van Meir,et al.  Exciting New Advances in Neuro‐Oncology: The Avenue to a Cure for Malignant Glioma , 2010, CA: a cancer journal for clinicians.

[38]  J. Salk Clonal evolution in cancer , 2010 .

[39]  Yuan Qi,et al.  Integrated Genomic Analysis Identifies Clinically Relevant Subtypes of Glioblastoma Characterized by Abnormalities in PDGFRA , IDH 1 , EGFR , and NF 1 Citation Verhaak , 2010 .

[40]  Larry Norton,et al.  Tumor Self-Seeding by Circulating Cancer Cells , 2009, Cell.

[41]  Terence P. Speed,et al.  A single-array preprocessing method for estimating full-resolution raw copy numbers from all Affymetrix genotyping arrays including GenomeWideSNP 5 & 6 , 2009, Bioinform..

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

[43]  Santosh Kesari,et al.  Malignant gliomas in adults. , 2008, The New England journal of medicine.

[44]  R. Gilbertson,et al.  Making a tumour's bed: glioblastoma stem cells and the vascular niche , 2007, Nature Reviews Cancer.

[45]  Mark W. Dewhirst,et al.  Glioma stem cells promote radioresistance by preferential activation of the DNA damage response , 2006, Nature.

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

[47]  Marina Konopleva,et al.  MDM2 antagonists induce p53-dependent apoptosis in AML: implications for leukemia therapy. , 2005, Blood.

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

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

[50]  M. Nowak,et al.  Evolutionary Dynamics of Biological Games , 2004, Science.

[51]  R. Solé,et al.  Metapopulation dynamics and spatial heterogeneity in cancer , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[52]  V. Rotter,et al.  The onset of p53-dependent DNA repair or apoptosis is determined by the level of accumulated damaged DNA. , 2002, Carcinogenesis.

[53]  Erwin G. Van Meir,et al.  Initiation of human astrocytoma by clonal evolution of cells with progressive loss of p53 functions in a patient with a 283H TP53 germ-line mutation: evidence for a precursor lesion. , 2002, Cancer research.

[54]  Erwin G. Van Meir,et al.  Cells with TP53 mutations in low grade astrocytic tumors evolve clonally to malignancy and are an unfavorable prognostic factor , 1999, Oncogene.

[55]  F. Cleton Evolution of Cancer , 1991, British Journal of Cancer.