Transcriptome Analysis of Recurrently Deregulated Genes Across Multiple

Total Page:16

File Type:pdf, Size:1020Kb

Transcriptome Analysis of Recurrently Deregulated Genes Across Multiple Published OnlineFirst November 9, 2015; DOI: 10.1158/0008-5472.CAN-15-0484 Cancer Integrated Systems and Technologies Research Transcriptome Analysis of Recurrently Deregulated Genes across Multiple Cancers Identifies New Pan-Cancer Biomarkers Bogumil Kaczkowski1, Yuji Tanaka1,2, Hideya Kawaji1,2,3, Albin Sandelin4, Robin Andersson4, Masayoshi Itoh1,3, Timo Lassmann1,5, the FANTOM5 consortium, Yoshihide Hayashizaki3, Piero Carninci1, and Alistair R.R. Forrest1,6 Abstract Genes that are commonly deregulated in cancer are clinically which are upregulated in cancer, defining promoters that overlap attractive as candidate pan-diagnostic markers and therapeutic with repetitive elements (especially SINE/Alu and LTR/ERV1 targets. To globally identify such targets, we compared Cap elements) that are often upregulated in cancer. Lastly, we docu- Analysis of Gene Expression profiles from 225 different cancer mented for the first time upregulation of multiple copies of cell lines and 339 corresponding primary cell samples to identify the REP522 interspersed repeat in cancer. Overall, our genome- transcripts that are deregulated recurrently in a broad range of wide expression profiling approach identified a comprehensive cancer types. Comparing RNA-seq data from 4,055 tumors and set of candidate biomarkers with pan-cancer potential, and 563 normal tissues profiled in the The Cancer Genome Atlas and extended the perspective and pathogenic significance of repetitive FANTOM5 datasets, we identified a core transcript set with ther- elements that are frequently activated during cancer progression. anostic potential. Our analyses also revealed enhancer RNAs, Cancer Res; 76(2); 216–26. Ó2015 AACR. Introduction time, cancers from different tissues can share some common features, for example, The Cancer Genome Atlas (TCGA) has Successful cancer treatment depends heavily on early detection found genes and pathways, DNA copy number alterations, muta- and diagnosis. Despite decades of research, relatively few bio- tions, methylation, and transcriptome changes that recur across markers are routinely used in clinics (e.g., CA-125 and PSA in 12 different primary tumor types (4). ovarian and prostate cancers, respectively; refs. 1, 2). There is a Here using Cap Analysis of Gene Expression (CAGE) data need for reliable and clinically applicable new cancer biomarkers collected for the Functional ANnoTation Of Mammalian for early detection. Cancers originating in the same tissue can be genome (FANTOM5) project (5), we identified mRNAs, very heterogeneous, often being derived from different cell types long-noncoding RNAs (lncRNA), enhancer RNAs (eRNA), and and having drastically different mutation profiles (3). At the same RNAs initiating from within repeat elements, which are recur- rently perturbed in cancer cell lines. To confirm that these transcripts are relevant to tumors, we compared their expres- 1RIKEN Center for Life Science Technologies, Division of Genomic Technologies, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa, sion in 4,055 primary tumors and 563 matching tissue sets Japan. 2RIKEN Advanced Center for Computing and Communication, RNA-seq profiled by the TCGA (6) and in a set of colorectal Preventive Medicine and Applied Genomics unit, 1-7-22 Suehiro-cho, tumor (7) samples profiled proteomically. Finally, for the most Tsurumi-ku, Yokohama, Japan. 3RIKEN Preventive Medicine & Diag- nosis Innovation Program, 2-1 Hirosawa, Wako, Saitama 351-0198, promising biomarker candidates we performed qRT-PCR vali- Japan. 4The Bioinformatics Centre, Department of Biology and Bio- dations in cancer cell lines and tumor cDNA panels. Taken tech Research and Innovation Centre (BRIC), University of Copenha- together, our analyses allowed for identification of a set of gen, Copenhagen, Denmark. 5Telethon Kids Institute, the University of Western Australia, Perth, Western Australia, Australia. 6Harry Perkins robust pan cancer biomarker candidates, which have the poten- Institute of Medical Research, QEII Medical Centre and Centre for tial for development as blood biomarkers for early detection Medical Research, the University of Western Australia, Nedlands, and for histological screening of biopsies. Western Australia, Australia. This work is part of the FANTOM5 project. Data download, Note: Supplementary data for this article are available at Cancer Research genomic tools, and copublished manuscripts have been summa- Online (http://cancerres.aacrjournals.org/). rized at the FANTOM5 website (8). Corresponding Authors: Bogumil Kaczkowski, RIKEN Center for Life Science Technologies (CLST), 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Materials and Methods Japan. Phone: 81-45-503-9222; Fax: 81-45-503-9216; E-mail: [email protected]; and Alistair R.R. Forrest, QEII Medical Centre and FANTOM5 data Centre for Medical Research, the University of Western Australia, 6 Verdun We used the cap analysis of gene expression (CAGE) data from Street, Nedlands, WA 6009, Australia. Phone: 61-8-6151-0780; Fax: 61-8-6151- the FANTOM5 project (libraries sequenced to a median depth of 4 0701; E-mail: [email protected] million mapped tags; ref. 5). We used 564 CAGE profiles: 225 doi: 10.1158/0008-5472.CAN-15-0484 cancer cell lines and 339 primary cells samples. We split the data Ó2015 American Association for Cancer Research. into three data sets: (i) matched solid, (ii) unmatched solid, and 216 Cancer Res; 76(2) January 15, 2016 Downloaded from cancerres.aacrjournals.org on September 30, 2021. © 2016 American Association for Cancer Research. Published OnlineFirst November 9, 2015; DOI: 10.1158/0008-5472.CAN-15-0484 Pan-Cancer Transcriptome A FANTOM5 DATA 225 cancer cell lines 339 primary cells 10 MATCHED origins UNMATCHED origins 2 MATCHED origins SOLID tumors SOLID cancers BLOOD cancers 72 cancer cell lines 102 cancer cell lines 51 cancer cell lines 65 primary cells 200 primary cells 74 primary cells edgeR and ON/OFF DEE pipelinei li DEE ppipeline ipelin DE pipeline analysis Solid cancer only differential expression Overlapping features Pan-cancer differential expression B CANCER Switching Expression shift NORMAL Figure 1. ON/OFF UP/DOWN p1@TERT p1@POLQ Summary of comparisons carried out to identify recurrently perturbed transcripts in the FANTOM5 cell line Upregulated dataset. A, differential expression (DE) pipeline applied to the FANTOM5 data. B, examples of differentially expressed promoters showing expression Melanocytes Mesothelium Melanocytes Mesothelium Lymphoid Myeloid Bone switching (ON and OFF) and expression Brain Kidney Liver Lung Breast Prostate Ovary Lymphoid Myeloid Bone Brain Kidney Liver Lung Breast Prostate Ovary shift (UP and DOWN). C, comparison p1@NAALADL1 p1@C13orf15 3 between promoter and gene level 2 differential expression (based on CAGE Downregulated data). Note: Although the majority of differentially expressed promoters reflect gene-wise differential expression, a significant fraction behave differently, for example, MPP2 or BCAT1. Lymphoid Brain Kidney Liver Lung Melanocytes Mesothelium Myeloid Bone Breast Prostate Ovary Melanocytes Mesothelium Lymphoid Myeloid Bone D, table summarizing the number of Brain Kidney Liver Lung Breast Prostate Ovary promoters and genes showing Blood Solid Blood Solid differential expression. Numbers in matched matched matched matched parentheses indicate numbers of unique C genes. FC 2 Promoter level log Promoter level Gene level log FC D 2 ′ Wi Non Gene- Promoter Gene- www.aacrjournals.org Cancer Res; 76(2) January 15, 2016 217 Downloaded from cancerres.aacrjournals.org on September 30, 2021. © 2016 American Association for Cancer Research. Published OnlineFirst November 9, 2015; DOI: 10.1158/0008-5472.CAN-15-0484 Kaczkowski et al. Table 1. The numbers of differentially expressed promoters (DPI clusters) and the type of genomic region they overlap, based on Gencode 19 Upregulated Downregulated Type of genomic region Pan cancer Solid only Pan cancer Solid only Total % Protein coding 434 455 92 354 1,335 63 LincRNA 45 38 2 7 92 4 Antisense 37 28 0 4 69 3 Pseudogene 12 9 2 4 27 1 Other ncRNAs 20 33 0 5 58 3 Unannotated 233 251 3 40 527 25 Total 781 814 99 414 2,108 100 (iii) matched blood (Supplementary Table S1A and S1C for list of upper quartile normalized RSEM count estimates with expression cancer types and sample annotation). The CAGE tag counts profiles of 20,531 genes in 4,618 samples. under 184,827 robust decomposition-based peak identification The counts were log2 transformed and used as an input (DPI) clusters (5) were used to represent a promoter-level expres- expression data to LIMMA. sion. For the enhancer activity, we used the CAGE tags counts The cancer versus normal comparison was performed using under 43,011 enhancer regions identified in ref. 9. equal weight for each solid cancer type, each type contributing equally to overall comparison. The P-values were adjusted for FANTOM5 differential expression analysis multiple testing by Benjamini–Hochberg method. The thresholds To identify up- and downregulated transcripts in cancer cell of fold change >2 and FDR <0.01 were used. lines versus normal primary cells, we used Genewise Negative Binomial Generalized Linear Models as implemented in edgeR Enrichment for cancer-related genes (10). The cancer versus normal comparison was performed We tested for the enrichment by applying a hypergeometric using glmLRT function. In matched solid comparison, we set test, using the significance threshold of P < 0.05. The list of equal weight for each solid cancer type,
Recommended publications
  • In-Depth Analysis of Genetic Variation Associated with Severe West Nile Viral Disease
    Article In-Depth Analysis of Genetic Variation Associated with Severe West Nile Viral Disease Megan E. Cahill 1, Mark Loeb 2, Andrew T. Dewan 1 and Ruth R. Montgomery 3,* 1 Center for Perinatal, Pediatric and Environmental Epidemiology, Department of Chronic Disease Epidemiology, Yale School of Public Health, 1 Church Street, New Haven, CT 06510, USA; [email protected] (M.E.C.); [email protected] (A.T.D.) 2 3208 Michael DeGroote Centre for Learning & Discovery, Division of Clinical Pathology, McMaster University, Hamilton, ON L8S 4L8, Canada; [email protected] 3 Department of Internal Medicine, Yale School of Medicine, 300 Cedar Street, New Haven, CT 06520, USA * Correspondence: [email protected] Received: 30 October 2020; Accepted: 3 December 2020; Published: 8 December 2020 Abstract: West Nile virus (WNV) is a mosquito-borne virus which causes symptomatic disease in a minority of infected humans. To identify novel genetic variants associated with severe disease, we utilized data from an existing case-control study of WNV and included population controls for an expanded analysis. We conducted imputation and gene-gene interaction analysis in the largest and most comprehensive genetic study conducted to date for West Nile neuroinvasive disease (WNND). Within the imputed West Nile virus dataset (severe cases n = 381 and asymptomatic/mild controls = 441), we found novel loci within the MCF.2 Cell Line Derived Transforming Sequence Like (MCF2L) gene (rs9549655 and rs2297192) through the individual loci analyses, although none reached statistical significance. Incorporating population controls from the Wisconsin Longitudinal Study on Aging (n = 9012) did not identify additional novel variants, a possible reflection of the cohort’s inclusion of individuals who could develop mild or severe WNV disease upon infection.
    [Show full text]
  • Core Transcriptional Regulatory Circuitries in Cancer
    Oncogene (2020) 39:6633–6646 https://doi.org/10.1038/s41388-020-01459-w REVIEW ARTICLE Core transcriptional regulatory circuitries in cancer 1 1,2,3 1 2 1,4,5 Ye Chen ● Liang Xu ● Ruby Yu-Tong Lin ● Markus Müschen ● H. Phillip Koeffler Received: 14 June 2020 / Revised: 30 August 2020 / Accepted: 4 September 2020 / Published online: 17 September 2020 © The Author(s) 2020. This article is published with open access Abstract Transcription factors (TFs) coordinate the on-and-off states of gene expression typically in a combinatorial fashion. Studies from embryonic stem cells and other cell types have revealed that a clique of self-regulated core TFs control cell identity and cell state. These core TFs form interconnected feed-forward transcriptional loops to establish and reinforce the cell-type- specific gene-expression program; the ensemble of core TFs and their regulatory loops constitutes core transcriptional regulatory circuitry (CRC). Here, we summarize recent progress in computational reconstitution and biologic exploration of CRCs across various human malignancies, and consolidate the strategy and methodology for CRC discovery. We also discuss the genetic basis and therapeutic vulnerability of CRC, and highlight new frontiers and future efforts for the study of CRC in cancer. Knowledge of CRC in cancer is fundamental to understanding cancer-specific transcriptional addiction, and should provide important insight to both pathobiology and therapeutics. 1234567890();,: 1234567890();,: Introduction genes. Till now, one critical goal in biology remains to understand the composition and hierarchy of transcriptional Transcriptional regulation is one of the fundamental mole- regulatory network in each specified cell type/lineage.
    [Show full text]
  • Genomic Correlates of Relationship QTL Involved in Fore- Versus Hind Limb Divergence in Mice
    Loyola University Chicago Loyola eCommons Biology: Faculty Publications and Other Works Faculty Publications 2013 Genomic Correlates of Relationship QTL Involved in Fore- Versus Hind Limb Divergence in Mice Mihaela Palicev Gunter P. Wagner James P. Noonan Benedikt Hallgrimsson James M. Cheverud Loyola University Chicago, [email protected] Follow this and additional works at: https://ecommons.luc.edu/biology_facpubs Part of the Biology Commons Recommended Citation Palicev, M, GP Wagner, JP Noonan, B Hallgrimsson, and JM Cheverud. "Genomic Correlates of Relationship QTL Involved in Fore- Versus Hind Limb Divergence in Mice." Genome Biology and Evolution 5(10), 2013. This Article is brought to you for free and open access by the Faculty Publications at Loyola eCommons. It has been accepted for inclusion in Biology: Faculty Publications and Other Works by an authorized administrator of Loyola eCommons. For more information, please contact [email protected]. This work is licensed under a Creative Commons Attribution-Noncommercial-No Derivative Works 3.0 License. © Palicev et al., 2013. GBE Genomic Correlates of Relationship QTL Involved in Fore- versus Hind Limb Divergence in Mice Mihaela Pavlicev1,2,*, Gu¨ nter P. Wagner3, James P. Noonan4, Benedikt Hallgrı´msson5,and James M. Cheverud6 1Konrad Lorenz Institute for Evolution and Cognition Research, Altenberg, Austria 2Department of Pediatrics, Cincinnati Children‘s Hospital Medical Center, Cincinnati, Ohio 3Yale Systems Biology Institute and Department of Ecology and Evolutionary Biology, Yale University 4Department of Genetics, Yale University School of Medicine 5Department of Cell Biology and Anatomy, The McCaig Institute for Bone and Joint Health and the Alberta Children’s Hospital Research Institute for Child and Maternal Health, University of Calgary, Calgary, Canada 6Department of Anatomy and Neurobiology, Washington University *Corresponding author: E-mail: [email protected].
    [Show full text]
  • Identification of the Binding Partners for Hspb2 and Cryab Reveals
    Brigham Young University BYU ScholarsArchive Theses and Dissertations 2013-12-12 Identification of the Binding arP tners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactions and Non- Redundant Roles for Small Heat Shock Proteins Kelsey Murphey Langston Brigham Young University - Provo Follow this and additional works at: https://scholarsarchive.byu.edu/etd Part of the Microbiology Commons BYU ScholarsArchive Citation Langston, Kelsey Murphey, "Identification of the Binding Partners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactions and Non-Redundant Roles for Small Heat Shock Proteins" (2013). Theses and Dissertations. 3822. https://scholarsarchive.byu.edu/etd/3822 This Thesis is brought to you for free and open access by BYU ScholarsArchive. It has been accepted for inclusion in Theses and Dissertations by an authorized administrator of BYU ScholarsArchive. For more information, please contact [email protected], [email protected]. Identification of the Binding Partners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactions and Non-Redundant Roles for Small Heat Shock Proteins Kelsey Langston A thesis submitted to the faculty of Brigham Young University in partial fulfillment of the requirements for the degree of Master of Science Julianne H. Grose, Chair William R. McCleary Brian Poole Department of Microbiology and Molecular Biology Brigham Young University December 2013 Copyright © 2013 Kelsey Langston All Rights Reserved ABSTRACT Identification of the Binding Partners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactors and Non-Redundant Roles for Small Heat Shock Proteins Kelsey Langston Department of Microbiology and Molecular Biology, BYU Master of Science Small Heat Shock Proteins (sHSP) are molecular chaperones that play protective roles in cell survival and have been shown to possess chaperone activity.
    [Show full text]
  • Molecular Dynamics and Evolutionary Aspects of the Transition from the Fully Grown Oocyte to Embryo
    Downloaded from genesdev.cshlp.org on September 26, 2021 - Published by Cold Spring Harbor Laboratory Press Cracking the egg: molecular dynamics and evolutionary aspects of the transition from the fully grown oocyte to embryo Alexei V. Evsikov,1,5 Joel H. Graber,1 J. Michael Brockman,1,2 Aleš Hampl,3 Andrea E. Holbrook,1 Priyam Singh,1,2 John J. Eppig,1 Davor Solter,1,4 and Barbara B. Knowles1 1The Jackson Laboratory, Bar Harbor, Maine 04609, USA; 2 Bioinformatics Program, Boston University, Boston, Massachusetts 02215, USA; 3Masaryk University Brno and Institute of Experimental Medicine, 625 00 Brno, Czech Republic; 4Max Planck Institute of Immunobiology, 79108 Freiburg, Germany Fully grown oocytes (FGOs) contain all the necessary transcripts to activate molecular pathways underlying the oocyte-to-embryo transition (OET). To elucidate this critical period of development, an extensive survey of the FGO transcriptome was performed by analyzing 19,000 expressed sequence tags of the Mus musculus FGO cDNA library. Expression of 5400 genes and transposable elements is reported. For a majority of genes expressed in mouse FGOs, homologs transcribed in eggs of Xenopus laevis or Ciona intestinalis were found, pinpointing evolutionary conservation of most regulatory cascades underlying the OET in chordates. A large proportion of identified genes belongs to several gene families with oocyte-restricted expression, a likely result of lineage-specific genomic duplications. Gene loss by mutation and expression in female germline of retrotransposed genes specific to M. musculus is documented. These findings indicate rapid diversification of genes involved in female reproduction. Comparison of the FGO and two-cell embryo transcriptomes demarcated the processes important for oogenesis from those involved in OET and identified novel motifs in maternal mRNAs associated with transcript stability.
    [Show full text]
  • A Computational Approach for Defining a Signature of Β-Cell Golgi Stress in Diabetes Mellitus
    Page 1 of 781 Diabetes A Computational Approach for Defining a Signature of β-Cell Golgi Stress in Diabetes Mellitus Robert N. Bone1,6,7, Olufunmilola Oyebamiji2, Sayali Talware2, Sharmila Selvaraj2, Preethi Krishnan3,6, Farooq Syed1,6,7, Huanmei Wu2, Carmella Evans-Molina 1,3,4,5,6,7,8* Departments of 1Pediatrics, 3Medicine, 4Anatomy, Cell Biology & Physiology, 5Biochemistry & Molecular Biology, the 6Center for Diabetes & Metabolic Diseases, and the 7Herman B. Wells Center for Pediatric Research, Indiana University School of Medicine, Indianapolis, IN 46202; 2Department of BioHealth Informatics, Indiana University-Purdue University Indianapolis, Indianapolis, IN, 46202; 8Roudebush VA Medical Center, Indianapolis, IN 46202. *Corresponding Author(s): Carmella Evans-Molina, MD, PhD ([email protected]) Indiana University School of Medicine, 635 Barnhill Drive, MS 2031A, Indianapolis, IN 46202, Telephone: (317) 274-4145, Fax (317) 274-4107 Running Title: Golgi Stress Response in Diabetes Word Count: 4358 Number of Figures: 6 Keywords: Golgi apparatus stress, Islets, β cell, Type 1 diabetes, Type 2 diabetes 1 Diabetes Publish Ahead of Print, published online August 20, 2020 Diabetes Page 2 of 781 ABSTRACT The Golgi apparatus (GA) is an important site of insulin processing and granule maturation, but whether GA organelle dysfunction and GA stress are present in the diabetic β-cell has not been tested. We utilized an informatics-based approach to develop a transcriptional signature of β-cell GA stress using existing RNA sequencing and microarray datasets generated using human islets from donors with diabetes and islets where type 1(T1D) and type 2 diabetes (T2D) had been modeled ex vivo. To narrow our results to GA-specific genes, we applied a filter set of 1,030 genes accepted as GA associated.
    [Show full text]
  • The Inactive X Chromosome Is Epigenetically Unstable and Transcriptionally Labile in Breast Cancer
    Supplemental Information The inactive X chromosome is epigenetically unstable and transcriptionally labile in breast cancer Ronan Chaligné1,2,3,8, Tatiana Popova1,4, Marco-Antonio Mendoza-Parra5, Mohamed-Ashick M. Saleem5 , David Gentien1,6, Kristen Ban1,2,3,8, Tristan Piolot1,7, Olivier Leroy1,7, Odette Mariani6, Hinrich Gronemeyer*5, Anne Vincent-Salomon*1,4,6,8, Marc-Henri Stern*1,4,6 and Edith Heard*1,2,3,8 Extended Experimental Procedures Cell Culture Human Mammary Epithelial Cells (HMEC, Invitrogen) were grown in serum-free medium (HuMEC, Invitrogen). WI- 38, ZR-75-1, SK-BR-3 and MDA-MB-436 cells were grown in Dulbecco’s modified Eagle’s medium (DMEM; Invitrogen) containing 10% fetal bovine serum (FBS). DNA Methylation analysis. We bisulfite-treated 2 µg of genomic DNA using Epitect bisulfite kit (Qiagen). Bisulfite converted DNA was amplified with bisulfite primers listed in Table S3. All primers incorporated a T7 promoter tag, and PCR conditions are available upon request. We analyzed PCR products by MALDI-TOF mass spectrometry after in vitro transcription and specific cleavage (EpiTYPER by Sequenom®). For each amplicon, we analyzed two independent DNA samples and several CG sites in the CpG Island. Design of primers and selection of best promoter region to assess (approx. 500 bp) were done by a combination of UCSC Genome Browser (http://genome.ucsc.edu) and MethPrimer (http://www.urogene.org). All the primers used are listed (Table S3). NB: MAGEC2 CpG analysis have been done with a combination of two CpG island identified in the gene core. Analysis of RNA allelic expression profiles (based on Human SNP Array 6.0) DNA and RNA hybridizations were normalized by Genotyping console.
    [Show full text]
  • Evaluation of Chromatin Accessibility in Prefrontal Cortex of Individuals with Schizophrenia
    ARTICLE DOI: 10.1038/s41467-018-05379-y OPEN Evaluation of chromatin accessibility in prefrontal cortex of individuals with schizophrenia Julien Bryois 1, Melanie E. Garrett2, Lingyun Song3, Alexias Safi3, Paola Giusti-Rodriguez 4, Graham D. Johnson 3, Annie W. Shieh13, Alfonso Buil5, John F. Fullard6, Panos Roussos 6,7,8, Pamela Sklar6, Schahram Akbarian 6, Vahram Haroutunian 6,9, Craig A. Stockmeier 10, Gregory A. Wray3,11, Kevin P. White12, Chunyu Liu13, Timothy E. Reddy 3,14, Allison Ashley-Koch2,15, Patrick F. Sullivan 1,4,16 & Gregory E. Crawford 3,17 1234567890():,; Schizophrenia genome-wide association studies have identified >150 regions of the genome associated with disease risk, yet there is little evidence that coding mutations contribute to this disorder. To explore the mechanism of non-coding regulatory elements in schizophrenia, we performed ATAC-seq on adult prefrontal cortex brain samples from 135 individuals with schizophrenia and 137 controls, and identified 118,152 ATAC-seq peaks. These accessible chromatin regions in the brain are highly enriched for schizophrenia SNP heritability. Accessible chromatin regions that overlap evolutionarily conserved regions exhibit an even higher heritability enrichment, indicating that sequence conservation can further refine functional risk variants. We identify few differences in chromatin accessibility between cases and controls, in contrast to thousands of age-related differential accessible chromatin regions. Altogether, we characterize chromatin accessibility in the human prefrontal cortex, the effect of schizophrenia and age on chromatin accessibility, and provide evidence that our dataset will allow for fine mapping of risk variants. 1 Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, SE-17177 Stockholm, Sweden.
    [Show full text]
  • Glioblastoma Stem Cells Induce Quiescence in Surrounding Neural Stem Cells Via Notch Signalling
    bioRxiv preprint doi: https://doi.org/10.1101/856062; this version posted November 29, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Glioblastoma stem cells induce quiescence in surrounding neural stem cells via Notch signalling. Katerina Lawlor1, Maria Angeles Marques-Torrejon2, Gopuraja Dharmalingham3, Yasmine El-Azhar1, Michael D. Schneider1, Steven M. Pollard2§ and Tristan A. Rodríguez1§ 1National Heart and Lung Institute, Imperial College London, Hammersmith Hospital Campus, Du Cane Road, London W12 0NN, United Kingdom. 2 MRC Centre for Regenerative Medicine & Edinburgh Cancer Research UK Centre, University of Edinburgh, Edinburgh, UK. 3MRC London Institute of Medical Sciences, Institute of Clinical Sciences, Imperial College London, UK §Authors for correspondence: [email protected] and [email protected] Running title: Glioblastoma stem cell competition Keyword: Neural stem cells, quiescence, glioblastoma, Notch, cell competition 1 bioRxiv preprint doi: https://doi.org/10.1101/856062; this version posted November 29, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. 1 Abstract 2 There is increasing evidence suggesting that adult neural stem cells (NSCs) are a cell of 3 origin of glioblastoma, the most aggressive form of malignant glioma. The earliest stages of 4 hyperplasia are not easy to explore, but likely involve a cross-talk between normal and 5 transformed NSCs. How normal cells respond to this cross-talk and if they expand or are 6 outcompeted is poorly understood.
    [Show full text]
  • Supplemental Materials ZNF281 Enhances Cardiac Reprogramming
    Supplemental Materials ZNF281 enhances cardiac reprogramming by modulating cardiac and inflammatory gene expression Huanyu Zhou, Maria Gabriela Morales, Hisayuki Hashimoto, Matthew E. Dickson, Kunhua Song, Wenduo Ye, Min S. Kim, Hanspeter Niederstrasser, Zhaoning Wang, Beibei Chen, Bruce A. Posner, Rhonda Bassel-Duby and Eric N. Olson Supplemental Table 1; related to Figure 1. Supplemental Table 2; related to Figure 1. Supplemental Table 3; related to the “quantitative mRNA measurement” in Materials and Methods section. Supplemental Table 4; related to the “ChIP-seq, gene ontology and pathway analysis” and “RNA-seq” and gene ontology analysis” in Materials and Methods section. Supplemental Figure S1; related to Figure 1. Supplemental Figure S2; related to Figure 2. Supplemental Figure S3; related to Figure 3. Supplemental Figure S4; related to Figure 4. Supplemental Figure S5; related to Figure 6. Supplemental Table S1. Genes included in human retroviral ORF cDNA library. Gene Gene Gene Gene Gene Gene Gene Gene Symbol Symbol Symbol Symbol Symbol Symbol Symbol Symbol AATF BMP8A CEBPE CTNNB1 ESR2 GDF3 HOXA5 IL17D ADIPOQ BRPF1 CEBPG CUX1 ESRRA GDF6 HOXA6 IL17F ADNP BRPF3 CERS1 CX3CL1 ETS1 GIN1 HOXA7 IL18 AEBP1 BUD31 CERS2 CXCL10 ETS2 GLIS3 HOXB1 IL19 AFF4 C17ORF77 CERS4 CXCL11 ETV3 GMEB1 HOXB13 IL1A AHR C1QTNF4 CFL2 CXCL12 ETV7 GPBP1 HOXB5 IL1B AIMP1 C21ORF66 CHIA CXCL13 FAM3B GPER HOXB6 IL1F3 ALS2CR8 CBFA2T2 CIR1 CXCL14 FAM3D GPI HOXB7 IL1F5 ALX1 CBFA2T3 CITED1 CXCL16 FASLG GREM1 HOXB9 IL1F6 ARGFX CBFB CITED2 CXCL3 FBLN1 GREM2 HOXC4 IL1F7
    [Show full text]
  • Reconstructing Cell Cycle Pseudo Time-Series Via Single-Cell Transcriptome Data—Supplement
    School of Natural Sciences and Mathematics Reconstructing Cell Cycle Pseudo Time-Series Via Single-Cell Transcriptome Data—Supplement UT Dallas Author(s): Michael Q. Zhang Rights: CC BY 4.0 (Attribution) ©2017 The Authors Citation: Liu, Zehua, Huazhe Lou, Kaikun Xie, Hao Wang, et al. 2017. "Reconstructing cell cycle pseudo time-series via single-cell transcriptome data." Nature Communications 8, doi:10.1038/s41467-017-00039-z This document is being made freely available by the Eugene McDermott Library of the University of Texas at Dallas with permission of the copyright owner. All rights are reserved under United States copyright law unless specified otherwise. File name: Supplementary Information Description: Supplementary figures, supplementary tables, supplementary notes, supplementary methods and supplementary references. CCNE1 CCNE1 CCNE1 CCNE1 36 40 32 34 32 35 30 32 28 30 30 28 28 26 24 25 Normalized Expression Normalized Expression Normalized Expression Normalized Expression 26 G1 S G2/M G1 S G2/M G1 S G2/M G1 S G2/M Cell Cycle Stage Cell Cycle Stage Cell Cycle Stage Cell Cycle Stage CCNE1 CCNE1 CCNE1 CCNE1 40 32 40 40 35 30 38 30 30 28 36 25 26 20 20 34 Normalized Expression Normalized Expression Normalized Expression 24 Normalized Expression G1 S G2/M G1 S G2/M G1 S G2/M G1 S G2/M Cell Cycle Stage Cell Cycle Stage Cell Cycle Stage Cell Cycle Stage Supplementary Figure 1 | High stochasticity of single-cell gene expression means, as demonstrated by relative expression levels of gene Ccne1 using the mESC-SMARTer data. For every panel, 20 sample cells were randomly selected for each of the three stages, followed by plotting the mean expression levels at each stage.
    [Show full text]
  • Supplementary Table 1: Adhesion Genes Data Set
    Supplementary Table 1: Adhesion genes data set PROBE Entrez Gene ID Celera Gene ID Gene_Symbol Gene_Name 160832 1 hCG201364.3 A1BG alpha-1-B glycoprotein 223658 1 hCG201364.3 A1BG alpha-1-B glycoprotein 212988 102 hCG40040.3 ADAM10 ADAM metallopeptidase domain 10 133411 4185 hCG28232.2 ADAM11 ADAM metallopeptidase domain 11 110695 8038 hCG40937.4 ADAM12 ADAM metallopeptidase domain 12 (meltrin alpha) 195222 8038 hCG40937.4 ADAM12 ADAM metallopeptidase domain 12 (meltrin alpha) 165344 8751 hCG20021.3 ADAM15 ADAM metallopeptidase domain 15 (metargidin) 189065 6868 null ADAM17 ADAM metallopeptidase domain 17 (tumor necrosis factor, alpha, converting enzyme) 108119 8728 hCG15398.4 ADAM19 ADAM metallopeptidase domain 19 (meltrin beta) 117763 8748 hCG20675.3 ADAM20 ADAM metallopeptidase domain 20 126448 8747 hCG1785634.2 ADAM21 ADAM metallopeptidase domain 21 208981 8747 hCG1785634.2|hCG2042897 ADAM21 ADAM metallopeptidase domain 21 180903 53616 hCG17212.4 ADAM22 ADAM metallopeptidase domain 22 177272 8745 hCG1811623.1 ADAM23 ADAM metallopeptidase domain 23 102384 10863 hCG1818505.1 ADAM28 ADAM metallopeptidase domain 28 119968 11086 hCG1786734.2 ADAM29 ADAM metallopeptidase domain 29 205542 11085 hCG1997196.1 ADAM30 ADAM metallopeptidase domain 30 148417 80332 hCG39255.4 ADAM33 ADAM metallopeptidase domain 33 140492 8756 hCG1789002.2 ADAM7 ADAM metallopeptidase domain 7 122603 101 hCG1816947.1 ADAM8 ADAM metallopeptidase domain 8 183965 8754 hCG1996391 ADAM9 ADAM metallopeptidase domain 9 (meltrin gamma) 129974 27299 hCG15447.3 ADAMDEC1 ADAM-like,
    [Show full text]