Fate mapping of human glioblastoma reveals an invariant stem cell hierarchy

Human glioblastomas harbour a subpopulation of glioblastoma stem cells that drive tumorigenesis. However, the origin of intratumoural functional heterogeneity between glioblastoma cells remains poorly understood. Here we study the clonal evolution of barcoded glioblastoma cells in an unbiased way following serial xenotransplantation to define their individual fate behaviours. Independent of an evolving mutational signature, we show that the growth of glioblastoma clones in vivo is consistent with a remarkably neutral process involving a conserved proliferative hierarchy rooted in glioblastoma stem cells. In this model, slow-cycling stem-like cells give rise to a more rapidly cycling progenitor population with extensive self-maintenance capacity, which in turn generates non-proliferative cells. We also identify rare ‘outlier’ clones that deviate from these dynamics, and further show that chemotherapy facilitates the expansion of pre-existing drug-resistant glioblastoma stem cells. Finally, we show that functionally distinct glioblastoma stem cells can be separately targeted using epigenetic compounds, suggesting new avenues for glioblastoma-targeted therapy.

[1]  Joachim M. Buhmann,et al.  TMARKER: A free software toolkit for histopathological cell counting and staining estimation , 2013, Journal of pathology informatics.

[2]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

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

[4]  M. Hirst,et al.  Barcoding reveals complex clonal dynamics of de novo transformed human mammary cells , 2015, Nature.

[5]  Joaquín Dopazo,et al.  Qualimap: evaluating next-generation sequencing alignment data , 2012, Bioinform..

[6]  R. Henkelman,et al.  Identification of human brain tumour initiating cells , 2004, Nature.

[7]  Bo Wen,et al.  Pharmacologic inhibition of the Menin-MLL interaction blocks progression of MLL leukemia in vivo. , 2015, Cancer cell.

[8]  Peter Rodgers,et al.  eulerAPE: Drawing Area-Proportional 3-Venn Diagrams Using Ellipses , 2014, PloS one.

[9]  Martin A Nowak,et al.  The age incidence of chronic myeloid leukemia can be explained by a one-mutation model , 2006, Proceedings of the National Academy of Sciences.

[10]  J. Dick,et al.  MLL5 Orchestrates a Cancer Self-Renewal State by Repressing the Histone Variant H3.3 and Globally Reorganizing Chromatin. , 2015, Cancer cell.

[11]  Deanna M. Church,et al.  ClinVar: public archive of relationships among sequence variation and human phenotype , 2013, Nucleic Acids Res..

[12]  Steffen Dickopf,et al.  A model of breast cancer heterogeneity reveals vascular mimicry as a driver of metastasis , 2015, Nature.

[13]  Z. Szallasi,et al.  Sequenza: allele-specific copy number and mutation profiles from tumor sequencing data , 2014, Annals of oncology : official journal of the European Society for Medical Oncology.

[14]  A. Sánchez-Danés,et al.  Defining the clonal dynamics leading to mouse skin tumour initiation , 2016, Nature.

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

[16]  B. Simons Deep sequencing as a probe of normal stem cell fate and preneoplasia in human epidermis , 2015, Proceedings of the National Academy of Sciences.

[17]  B. Simons,et al.  Tracing cellular dynamics in tissue development, maintenance and disease. , 2016, Current opinion in cell biology.

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

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

[20]  Aaron R. Quinlan,et al.  BIOINFORMATICS APPLICATIONS NOTE , 2022 .

[21]  Marc J. Williams,et al.  Identification of neutral tumor evolution across cancer types , 2016, Nature Genetics.

[22]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[23]  Alex A. Pollen,et al.  Molecular Identity of Human Outer Radial Glia during Cortical Development , 2015, Cell.

[24]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[25]  M. Greaves Cancer stem cells: back to Darwin? , 2010, Seminars in cancer biology.

[26]  Gary D Bader,et al.  SINGLE CELL DERIVED CLONAL ANALYSIS OF HUMAN GLIOBLASTOMA LINKS FUNCTIONAL AND GENOMIC HETEROGENEITY , 2014, Neuro-Oncology.

[27]  M. Stratton,et al.  High burden and pervasive positive selection of somatic mutations in normal human skin , 2015, Science.

[28]  J. Herman,et al.  Predicting lung cancer by detecting aberrant promoter methylation in sputum. , 2000, Cancer research.

[29]  Richard Durbin,et al.  Fast and accurate long-read alignment with Burrows–Wheeler transform , 2010, Bioinform..

[30]  Sohrab P. Shah,et al.  Dynamics of genomic clones in breast cancer patient xenografts at single-cell resolution , 2014, Nature.

[31]  Trevor J Pugh,et al.  Oncotator: Cancer Variant Annotation Tool , 2015, Human mutation.

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

[33]  Andrew Mugler,et al.  Analytic methods for modeling stochastic regulatory networks. , 2010, Methods in molecular biology.

[34]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[35]  C. Gardiner Stochastic Methods: A Handbook for the Natural and Social Sciences , 2009 .

[36]  Mark Bernstein,et al.  Glioma stem cell lines expanded in adherent culture have tumor-specific phenotypes and are suitable for chemical and genetic screens. , 2009, Cell stem cell.

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

[38]  Tzong-Shiue Yu,et al.  A restricted cell population propagates glioblastoma growth after chemotherapy , 2012 .

[39]  Jos Jonkers,et al.  A high-throughput splinkerette-PCR method for the isolation and sequencing of retroviral insertion sites , 2009, Nature Protocols.

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

[41]  Gabor T. Marth,et al.  A global reference for human genetic variation , 2015, Nature.

[42]  Renata Walewska,et al.  Chromatin accessibility maps of chronic lymphocytic leukaemia identify subtype-specific epigenome signatures and transcription regulatory networks , 2016, Nature Communications.

[43]  Martin Hirst,et al.  DNA barcoding reveals diverse growth kinetics of human breast tumour subclones in serially passaged xenografts , 2014, Nature Communications.

[44]  Clifford A. Meyer,et al.  Model-based Analysis of ChIP-Seq (MACS) , 2008, Genome Biology.

[45]  Howard Y. Chang,et al.  Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position , 2013, Nature Methods.

[46]  Elizabeth M. Smigielski,et al.  dbSNP: the NCBI database of genetic variation , 2001, Nucleic Acids Res..

[47]  Joshua M. Korn,et al.  Studying clonal dynamics in response to cancer therapy using high-complexity barcoding , 2015, Nature Medicine.

[48]  Giovanni Coppola,et al.  Identification of an Efficient Gene Expression Panel for Glioblastoma Classification , 2016, PloS one.

[49]  C. Allis,et al.  Use of human embryonic stem cells to model pediatric gliomas with H3.3K27M histone mutation , 2014, Science.

[50]  Kenny Q. Ye,et al.  An integrated map of genetic variation from 1,092 human genomes , 2012, Nature.

[51]  D. E. Barton,et al.  The Elements of Stochastic Processes with Applications to the Natural Sciences , 1964 .

[52]  Nicolò Riggi,et al.  EZH2 is essential for glioblastoma cancer stem cell maintenance. , 2009, Cancer research.

[53]  Mingming Jia,et al.  COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer , 2010, Nucleic Acids Res..

[54]  Christopher A. Miller,et al.  VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. , 2012, Genome research.

[55]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[56]  Benjamin D. Simons,et al.  Defining the mode of tumour growth by clonal analysis , 2012, Nature.

[57]  C. Curtis,et al.  A Big Bang model of human colorectal tumor growth , 2015, Nature Genetics.

[58]  Mariella G. Filbin,et al.  Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma , 2016, Nature.

[59]  Steven J. M. Jones,et al.  Quiescent sox2(+) cells drive hierarchical growth and relapse in sonic hedgehog subgroup medulloblastoma. , 2014, Cancer cell.

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

[61]  M. Hirst,et al.  Clonal analysis via barcoding reveals diverse growth and differentiation of transplanted mouse and human mammary stem cells. , 2014, Cell stem cell.

[62]  I. Glauche,et al.  Limitations and challenges of genetic barcode quantification , 2017, Scientific Reports.

[63]  Martin A. Nowak,et al.  A spatial model predicts that dispersal and cell turnover limit intratumour heterogeneity , 2015, Nature.

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

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

[66]  N. Potter,et al.  Genetic and Functional Diversity of Propagating Cells in Glioblastoma , 2014, Stem cell reports.

[67]  G. Smyth,et al.  ELDA: extreme limiting dilution analysis for comparing depleted and enriched populations in stem cell and other assays. , 2009, Journal of immunological methods.

[68]  Gary D Bader,et al.  Single cell-derived clonal analysis of human glioblastoma links functional and genomic heterogeneity , 2015, Proceedings of the National Academy of Sciences.

[69]  James Y. Zou Analysis of protein-coding genetic variation in 60,706 humans , 2015, Nature.

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

[71]  Norman T. J. Bailey The Elements of Stochastic Processes with Applications to the Natural Sciences , 1964 .