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

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.

[1]  C. Ball,et al.  Identification of genes periodically expressed in the human cell cycle and their expression in tumors. , 2002, Molecular biology of the cell.

[2]  Rebecca A Betensky,et al.  Glioma test array for use with formalin-fixed, paraffin-embedded tissue: array comparative genomic hybridization correlates with loss of heterozygosity and fluorescence in situ hybridization. , 2006, The Journal of molecular diagnostics : JMD.

[3]  A. Iafrate,et al.  Expression of Oligodendroglial and Astrocytic Lineage Markers in Diffuse Gliomas: Use of YKL-40, ApoE, ASCL1, and NKX2-2 , 2006, Journal of neuropathology and experimental neurology.

[4]  D. Louis WHO classification of tumours of the central nervous system , 2007 .

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

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

[7]  K. Miyazono,et al.  Autocrine TGF- Signaling Maintains Tumorigenicity of Glioma-Initiating Cells through Sry-Related HM , 2009 .

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

[9]  Federico Calegari,et al.  Cdks and cyclins link G1 length and differentiation of embryonic, neural and hematopoietic stem cells , 2010, Cell cycle.

[10]  Kristian Cibulskis,et al.  ContEst: estimating cross-contamination of human samples in next-generation sequencing data , 2011, Bioinform..

[11]  L. Luo,et al.  Mosaic Analysis with Double Markers Reveals Tumor Cell of Origin in Glioma , 2011, Cell.

[12]  Colin N. Dewey,et al.  RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome , 2011, BMC Bioinformatics.

[13]  Noemi Andor,et al.  Asymmetry-defective oligodendrocyte progenitors are glioma precursors. , 2011, Cancer cell.

[14]  J. Doudna,et al.  Mutations in CIC and FUBP1 Contribute to Human Oligodendroglioma , 2011 .

[15]  Eric A Bushong,et al.  Dedifferentiation of Neurons and Astrocytes by Oncogenes Can Induce Gliomas in Mice , 2012, Science.

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

[17]  Trevor J Pugh,et al.  Discovery and characterization of artifactual mutations in deep coverage targeted capture sequencing data due to oxidative DNA damage during sample preparation , 2013, Nucleic acids research.

[18]  H Aburatani,et al.  The critical role of cyclin D2 in cell cycle progression and tumorigenicity of glioblastoma stem cells , 2013, Oncogene.

[19]  B. Bernstein,et al.  Epigenetic Reprogramming in Cancer , 2013, Science.

[20]  M. A. Khan,et al.  The chromatin remodeler CHD7 regulates adult neurogenesis via activation of SoxC transcription factors. , 2013, Cell stem cell.

[21]  Simon Kasif,et al.  An aberrant transcription factor network essential for Wnt signaling and stem cell maintenance in glioblastoma. , 2013, Cell reports.

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

[23]  S. Gygi,et al.  Identification of a Unique TGF-β Dependent Molecular and Functional Signature in Microglia , 2013, Nature Neuroscience.

[24]  Simon Kasif,et al.  Reconstructing and Reprogramming the Tumor-Propagating Potential of Glioblastoma Stem-like Cells , 2014, Cell.

[25]  A. Krešo,et al.  Evolution of the cancer stem cell model. , 2014, Cell stem cell.

[26]  Allan R. Jones,et al.  Transcriptional Landscape of the Prenatal Human Brain , 2014, Nature.

[27]  Åsa K. Björklund,et al.  Full-length RNA-seq from single cells using Smart-seq2 , 2014, Nature Protocols.

[28]  Rona S. Gertner,et al.  Single cell RNA Seq reveals dynamic paracrine control of cellular variation , 2014, Nature.

[29]  T. Maniatis,et al.  An RNA-Sequencing Transcriptome and Splicing Database of Glia, Neurons, and Vascular Cells of the Cerebral Cortex , 2014, The Journal of Neuroscience.

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

[31]  J. Rich,et al.  Cancer stem cells in glioblastoma , 2015, Genes & development.

[32]  A. Moiyadi,et al.  ETV/Pea3 family transcription factor‐encoding genes are overexpressed in CIC‐mutant oligodendrogliomas , 2015, Genes, chromosomes & cancer.

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

[34]  Evan Z. Macosko,et al.  Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets , 2015, Cell.

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

[36]  David W. Nauen,et al.  Single-Cell RNA-Seq with Waterfall Reveals Molecular Cascades underlying Adult Neurogenesis. , 2015, Cell stem cell.

[37]  Monika S. Kowalczyk,et al.  Single-cell RNA-seq reveals changes in cell cycle and differentiation programs upon aging of hematopoietic stem cells , 2015, Genome research.

[38]  A. Regev,et al.  Spatial reconstruction of single-cell gene expression data , 2015 .

[39]  Steven J. M. Jones,et al.  Comprehensive, Integrative Genomic Analysis of Diffuse Lower-Grade Gliomas. , 2015, The New England journal of medicine.

[40]  G. Reifenberger,et al.  The 2016 World Health Organization Classification of Tumors of the Central Nervous System: a summary , 2016, Acta Neuropathologica.

[41]  E. Chang,et al.  Purification and Characterization of Progenitor and Mature Human Astrocytes Reveals Transcriptional and Functional Differences with Mouse , 2016, Neuron.