Single‐cell sequencing maps gene expression to mutational phylogenies in PDGF‐ and EGF‐driven gliomas

Glioblastoma multiforme (GBM) is the most common and aggressive type of primary brain tumor. Epidermal growth factor (EGF) and platelet‐derived growth factor (PDGF) receptors are frequently amplified and/or possess gain‐of‐function mutations in GBM. However, clinical trials of tyrosine‐kinase inhibitors have shown disappointing efficacy, in part due to intra‐tumor heterogeneity. To assess the effect of clonal heterogeneity on gene expression, we derived an approach to map single‐cell expression profiles to sequentially acquired mutations identified from exome sequencing. Using 288 single cells, we constructed high‐resolution phylogenies of EGF‐driven and PDGF‐driven GBMs, modeling transcriptional kinetics during tumor evolution. Descending the phylogenetic tree of a PDGF‐driven tumor corresponded to a progressive induction of an oligodendrocyte progenitor‐like cell type, expressing pro‐angiogenic factors. In contrast, phylogenetic analysis of an EGFR‐amplified tumor showed an up‐regulation of pro‐invasive genes. An in‐frame deletion in a specific dimerization domain of PDGF receptor correlates with an up‐regulation of growth pathways in a proneural GBM and enhances proliferation when ectopically expressed in glioma cell lines. In‐frame deletions in this domain are frequent in public GBM data.

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

[2]  Max A. Horlbeck,et al.  Single-cell analysis of long non-coding RNAs in the developing human neocortex , 2016, Genome Biology.

[3]  Hsien-Da Huang,et al.  miRTarBase 2016: updates to the experimentally validated miRNA-target interactions database , 2015, Nucleic Acids Res..

[4]  Gabor T. Marth,et al.  An integrated map of structural variation in 2,504 human genomes , 2015, Nature.

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

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

[7]  Michael C. Schatz,et al.  Interactive analysis and assessment of single-cell copy-number variations , 2015, Nature Methods.

[8]  V. Yong,et al.  ADAM-9 is a novel mediator of tenascin-C-stimulated invasiveness of brain tumor-initiating cells. , 2015, Neuro-oncology.

[9]  Timothy C Ryken,et al.  Toward precision medicine in glioblastoma: the promise and the challenges. , 2015, Neuro-oncology.

[10]  Sören Müller,et al.  Principles of miRNA–mRNA interactions: beyond sequence complementarity , 2015, Cellular and Molecular Life Sciences.

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

[12]  Parag Mallick,et al.  Neuronal Activity Promotes Glioma Growth through Neuroligin-3 Secretion , 2015, Cell.

[13]  Steven L Salzberg,et al.  HISAT: a fast spliced aligner with low memory requirements , 2015, Nature Methods.

[14]  William A. Flavahan,et al.  Periostin Secreted by Glioblastoma Stem Cells Recruits M2 Tumor-associated Macrophages and Promotes Malignant Growth , 2014, Nature Cell Biology.

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

[16]  Y. Pekarsky,et al.  Role of miR-15/16 in CLL , 2014, Cell Death and Differentiation.

[17]  Michael C. Schatz,et al.  Interactive analysis and quality assessment of single-cell copy-number variations , 2014, bioRxiv.

[18]  Mariano J. Alvarez,et al.  Identification of Causal Genetic Drivers of Human Disease through Systems-Level Analysis of Regulatory Networks , 2014, Cell.

[19]  Robin L. Anderson,et al.  BMP4 inhibits breast cancer metastasis by blocking myeloid-derived suppressor cell activity. , 2014, Cancer research.

[20]  Scott Chamberlain,et al.  Rphylip: an R interface for PHYLIP , 2014 .

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

[22]  A. Kriegstein,et al.  Control of outer radial glial stem cell mitosis in the human brain. , 2014, Cell Reports.

[23]  Harald Sontheimer,et al.  A neurocentric perspective on glioma invasion , 2014, Nature Reviews Neuroscience.

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

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

[26]  P. Kharchenko,et al.  Bayesian approach to single-cell differential expression analysis , 2014, Nature Methods.

[27]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

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

[29]  J. Huse,et al.  Osteopontin-CD44 signaling in the glioma perivascular niche enhances cancer stem cell phenotypes and promotes aggressive tumor growth. , 2014, Cell stem cell.

[30]  Charlotte Soneson,et al.  A glioma classification scheme based on coexpression modules of EGFR and PDGFRA , 2014, Proceedings of the National Academy of Sciences.

[31]  P. Verrelle,et al.  Role of Akt in human malignant glioma: from oncogenesis to tumor aggressiveness , 2014, Journal of Neuro-Oncology.

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

[33]  Natalio Krasnogor,et al.  JEPETTO: a Cytoscape plugin for gene set enrichment and topological analysis based on interaction networks , 2013, Bioinform..

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

[35]  Barbara S. Paugh,et al.  Novel oncogenic PDGFRA mutations in pediatric high-grade gliomas. , 2013, Cancer research.

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

[37]  Xiaoyan Chen,et al.  Platelet-derived growth factors and their receptors: structural and functional perspectives. , 2013, Biochimica et biophysica acta.

[38]  F. Guillemot,et al.  Ascl1/Mash1 Promotes Brain Oligodendrogenesis during Myelination and Remyelination , 2013, The Journal of Neuroscience.

[39]  Jing Wang,et al.  WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013 , 2013, Nucleic Acids Res..

[40]  K. Aldape,et al.  αvβ8 integrin interacts with RhoGDI1 to regulate Rac1 and Cdc42 activation and drive glioblastoma cell invasion , 2013, Molecular biology of the cell.

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

[42]  Michael Weller,et al.  Standards of care for treatment of recurrent glioblastoma--are we there yet? , 2013, Neuro-oncology.

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

[44]  Y. Suh,et al.  c-Jun N-terminal kinase has a pivotal role in the maintenance of self-renewal and tumorigenicity in glioma stem-like cells , 2012, Oncogene.

[45]  B. Vought,et al.  LINGO-1, a Transmembrane Signaling Protein, Inhibits Oligodendrocyte Differentiation and Myelination through Intercellular Self-interactions , 2012, The Journal of Biological Chemistry.

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

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

[48]  W. Fellows-Mayle,et al.  COX-2 blockade suppresses gliomagenesis by inhibiting myeloid-derived suppressor cells. , 2011, Cancer research.

[49]  Gary D. Bader,et al.  Pathway Commons, a web resource for biological pathway data , 2010, Nucleic Acids Res..

[50]  M. Wegner,et al.  Olig2 regulates Sox10 expression in oligodendrocyte precursors through an evolutionary conserved distal enhancer , 2010, Nucleic acids research.

[51]  Gary D. Bader,et al.  GeneMANIA Cytoscape plugin: fast gene function predictions on the desktop , 2010, Bioinform..

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

[53]  A. Epstein,et al.  Characterization of Cytokine-Induced Myeloid-Derived Suppressor Cells from Normal Human Peripheral Blood Mononuclear Cells , 2010, The Journal of Immunology.

[54]  Li Guo,et al.  Preparation of DNA Ladder Based on Multiplex PCR Technique , 2010, Journal of nucleic acids.

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

[56]  C. Croce,et al.  miR-15a and miR-16-1 in cancer: discovery, function and future perspectives , 2010, Cell Death and Differentiation.

[57]  Aaron R. Quinlan,et al.  Bioinformatics Applications Note Genome Analysis Bedtools: a Flexible Suite of Utilities for Comparing Genomic Features , 2022 .

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

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

[60]  F. Lazeyras,et al.  Cooperative expression of junctional adhesion molecule‐C and ‐B supports growth and invasion of glioma , 2009, Glia.

[61]  P. C. de Witt Hamer Small molecule kinase inhibitors in glioblastoma: a systematic review of clinical studies. , 2010, Neuro-oncology.

[62]  Kai Ye,et al.  Pindel: a pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads , 2009, Bioinform..

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

[64]  Mauro Biffoni,et al.  The miR-15a–miR-16-1 cluster controls prostate cancer by targeting multiple oncogenic activities , 2008, Nature Medicine.

[65]  John D. Lambris,et al.  Modulation of the anti-tumor immune response by complement , 2008, Nature Immunology.

[66]  Christopher Beadle,et al.  The role of myosin II in glioma invasion of the brain. , 2008, Molecular biology of the cell.

[67]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[68]  Keith L. Ligon,et al.  Coactivation of Receptor Tyrosine Kinases Affects the Response of Tumor Cells to Targeted Therapies , 2007, Science.

[69]  J. Rubenstein,et al.  Dlx1 and Dlx2 Control Neuronal versus Oligodendroglial Cell Fate Acquisition in the Developing Forebrain , 2007, Neuron.

[70]  Shannan J. Ho Sui,et al.  oPOSSUM: integrated tools for analysis of regulatory motif over-representation , 2007, Nucleic Acids Res..

[71]  Oscar Gonzalez-Perez,et al.  Origin of Oligodendrocytes in the Subventricular Zone of the Adult Brain , 2006, The Journal of Neuroscience.

[72]  Palma Iannarelli,et al.  Competing waves of oligodendrocytes in the forebrain and postnatal elimination of an embryonic lineage , 2006, Nature Neuroscience.

[73]  J. Stewart,et al.  CD155/PVR enhances glioma cell dispersal by regulating adhesion signaling and focal adhesion dynamics. , 2005, Cancer research.

[74]  J. Baudier,et al.  Nuclear expression of S100B in oligodendrocyte progenitor cells correlates with differentiation toward the oligodendroglial lineage and modulates oligodendrocytes maturation , 2004, Molecular and Cellular Neuroscience.

[75]  W. Mandemakers,et al.  An Oligodendrocyte Lineage-Specific Semaphorin, Sema5A, Inhibits Axon Growth by Retinal Ganglion Cells , 2004, The Journal of Neuroscience.

[76]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[77]  A. Barclay,et al.  Membrane proteins with immunoglobulin-like domains--a master superfamily of interaction molecules. , 2003, Seminars in immunology.

[78]  P. Dirks,et al.  A human brain tumor-derived PDGFR-α deletion mutant is transforming , 2003, Oncogene.

[79]  H. Katayama,et al.  Increased mitotic phosphorylation of histone H3 attributable to AIM-1/Aurora-B overexpression contributes to chromosome number instability. , 2002, Cancer research.

[80]  David J. Anderson,et al.  The bHLH Transcription Factors OLIG2 and OLIG1 Couple Neuronal and Glial Subtype Specification , 2002, Cell.

[81]  M. Wegner,et al.  Terminal differentiation of myelin-forming oligodendrocytes depends on the transcription factor Sox10. , 2002, Genes & development.

[82]  P. Slocombe,et al.  Meltrin gamma(ADAM-9) mediates cellular adhesion through alpha(6)beta(1 )integrin, leading to a marked induction of fibroblast cell motility. , 2000, Journal of cell science.

[83]  P. Davies,et al.  Novel microtubule‐associated protein‐2 isoform is expressed early in human oligodendrocyte maturation , 2000, Glia.

[84]  H. Katayama,et al.  Multinuclearity and increased ploidy caused by overexpression of the aurora- and Ipl1-like midbody-associated protein mitotic kinase in human cancer cells. , 1998, Cancer research.

[85]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[86]  T. Yoshimoto,et al.  Amplification of alpha-platelet-derived growth factor receptor gene lacking an exon coding for a portion of the extracellular region in a primary brain tumor of glial origin. , 1992, Oncogene.

[87]  W. Fitch,et al.  Construction of phylogenetic trees. , 1967, Science.

[88]  the original work is properly cited. , 2022 .