The phenotypes of proliferating glioblastoma cells reside on a single axis of variation.

Although tumor-propagating cells can be derived from glioblastomas (GBMs) of the proneural and mesenchymal subtypes, a glioma stem-like cell (GSC) of the classical subtype has not been identified. It is unclear if mesenchymal GSCs (mGSCs) and/or proneural GSCs (pGSCs) alone are sufficient to generate the heterogeneity observed in GBM. We performed single-cell/nuclei RNA-sequencing of 28 gliomas, and single-cell ATAC-sequencing for 8 cases. We find that GBM GSCs reside on a single axis of variation, ranging from proneural to mesenchymal. In silico lineage tracing using both transcriptomics and genetics supports mGSCs as the progenitors of pGSCs. Dual inhibition of pGSC-enriched and mGSC-enriched growth and survival pathways provides a more complete treatment than combinations targeting one GSC phenotype alone. This study sheds light on a long-standing debate regarding lineage relationships among GSCs and presents a paradigm by which personalized combination therapies can be derived from single-cell RNA signatures, to overcome intra-tumor heterogeneity.

[1]  Wei Shi,et al.  featureCounts: an efficient general purpose program for assigning sequence reads to genomic features , 2013, Bioinform..

[2]  Ming Tang,et al.  TumorFusions: an integrative resource for cancer-associated transcript fusions , 2017, Nucleic Acids Res..

[3]  B. Ren,et al.  Fast and Accurate Clustering of Single Cell Epigenomes Reveals Cis-Regulatory Elements in Rare Cell Types , 2019 .

[4]  S. Quake,et al.  A survey of human brain transcriptome diversity at the single cell level , 2015, Proceedings of the National Academy of Sciences.

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

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

[7]  A. Falini,et al.  The proneural gene ASCL1 governs the transcriptional subgroup affiliation in glioblastoma stem cells by directly repressing the mesenchymal gene NDRG1 , 2018, Cell Death & Differentiation.

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

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

[10]  Emanuel F Petricoin,et al.  Discordant inheritance of chromosomal and extrachromosomal DNA elements contributes to dynamic disease evolution in glioblastoma , 2018, Nature Genetics.

[11]  W. Krzyzosiak,et al.  Sequence-non-specific effects of RNA interference triggers and microRNA regulators , 2009, Nucleic acids research.

[12]  Martin J. Aryee,et al.  Lineage Tracing in Humans Enabled by Mitochondrial Mutations and Single-Cell Genomics , 2019, Cell.

[13]  Paul Hoffman,et al.  Integrating single-cell transcriptomic data across different conditions, technologies, and species , 2018, Nature Biotechnology.

[14]  Sören Müller,et al.  CONICS integrates scRNA-seq with DNA sequencing to map gene expression to tumor sub-clones , 2018, Bioinform..

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

[16]  Mariella G. Filbin,et al.  An Integrative Model of Cellular States, Plasticity, and Genetics for Glioblastoma , 2019, Cell.

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

[18]  Alex A. Pollen,et al.  Spatiotemporal gene expression trajectories reveal developmental hierarchies of the human cortex , 2017, Science.

[19]  Lisa C. Wallace,et al.  Targeting Glioma Stem Cells through Combined BMI1 and EZH2 Inhibition , 2017, Nature Medicine.

[20]  S. Halgamuge,et al.  Inferring copy number and genotype in tumour exome data , 2014, BMC Genomics.

[21]  Erik Sundström,et al.  RNA velocity of single cells , 2018, Nature.

[22]  A. Kriegstein,et al.  Single-cell profiling of human gliomas reveals macrophage ontogeny as a basis for regional differences in macrophage activation in the tumor microenvironment , 2017, Genome Biology.

[23]  R. Verhaak,et al.  GlioVis data portal for visualization and analysis of brain tumor expression datasets. , 2017, Neuro-oncology.

[24]  Aaron Diaz,et al.  SCell: integrated analysis of single-cell RNA-seq data , 2016, Bioinform..

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

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

[27]  Ulrike Groemping,et al.  Relative Importance for Linear Regression in R: The Package relaimpo , 2006 .

[28]  Lindy E. Barrett,et al.  Self-renewal does not predict tumor growth potential in mouse models of high-grade glioma. , 2012, Cancer cell.

[29]  Jeffrey T Leek,et al.  Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown , 2016, Nature Protocols.

[30]  Mariella G. Filbin,et al.  Decoupling genetics, lineages, and microenvironment in IDH-mutant gliomas by single-cell RNA-seq , 2017, Science.

[31]  Sören Müller,et al.  Single-cell Co-expression Subnetwork Analysis , 2017, Scientific Reports.

[32]  Alex A. Pollen,et al.  Single‐cell sequencing maps gene expression to mutational phylogenies in PDGF‐ and EGF‐driven gliomas , 2016, Molecular Systems Biology.