Heterochromatin Protein 1A Mediates Development and Aggressiveness of Neuroendocrine Prostate Cancer Xinpei Ci1,2, Jun Hao1,2, Xin Dong2, Stephen Y

Total Page:16

File Type:pdf, Size:1020Kb

Heterochromatin Protein 1A Mediates Development and Aggressiveness of Neuroendocrine Prostate Cancer Xinpei Ci1,2, Jun Hao1,2, Xin Dong2, Stephen Y Published OnlineFirst February 27, 2018; DOI: 10.1158/0008-5472.CAN-17-3677 Cancer Tumor Biology and Immunology Research Heterochromatin Protein 1a Mediates Development and Aggressiveness of Neuroendocrine Prostate Cancer Xinpei Ci1,2, Jun Hao1,2, Xin Dong2, Stephen Y. Choi1,2, Hui Xue2, Rebecca Wu2, Sifeng Qu1,2, Peter W. Gout2, Fang Zhang2, Anne M. Haegert1, Ladan Fazli1, Francesco Crea3, Christopher J. Ong1, Amina Zoubeidi1, Housheng H. He4,5, Martin E. Gleave1, Colin C. Collins1, Dong Lin1,2, and Yuzhuo Wang1,2 Abstract Neuroendocrine prostate cancer (NEPC) is a lethal subtype NEPC samples. HP1a knockdown in the NCI-H660 NEPC cell of prostate cancer arising mostly from adenocarcinoma via line inhibited proliferation, ablated colony formation, and neuroendocrine transdifferentiation following androgen dep- induced apoptotic cell death, ultimately leading to tumor rivation therapy. Mechanisms contributing to both NEPC growth arrest. Its ectopic expression significantly promoted development and its aggressiveness remain elusive. In light of NE transdifferentiation in adenocarcinoma cells subjected to the fact that hyperchromatic nuclei are a distinguishing histo- androgen deprivation treatment. Mechanistically, HP1a pathologic feature of NEPC, we utilized transcriptomic anal- reduced expression of androgen receptor and RE1 silencing yses of our patient-derived xenograft (PDX) models, multiple transcription factor and enriched the repressive trimethylated clinical cohorts, and genetically engineered mouse models to histone H3 at Lys9 mark on their respective gene promoters. identify 36 heterochromatin-related genes that are significantly These observations indicate a novel mechanism underlying enriched in NEPC. Longitudinal analysis using our unique, NEPC development mediated by abnormally expressed hetero- first-in-field PDX model of adenocarcinoma-to-NEPC transdif- chromatin genes, with HP1a as an early functional mediator ferentiation revealed that, among those 36 heterochromatin- and a potential therapeutic target for NEPC prevention and related genes, heterochromatin protein 1a (HP1a) expression management. increased early and steadily during NEPC development and Significance: Heterochromatin proteins play a fundamental remained elevated in the developed NEPC tumor. Its elevated role in NEPC, illuminating new therapeutic targets for this aggres- expression was further confirmed in multiple PDX and clinical sive disease. Cancer Res; 78(10); 2691–704. Ó2018 AACR. Introduction androgen deprivation treatment (ADT), resisting dependence on AR signaling as an adaptive response. As next-generation AR Neuroendocrine prostate cancer (NEPC) has become a clinical pathway inhibitors (ARPI) such as enzalutamide (EnZ) and challenge in the management of castration-resistant prostate abiraterone have made substantial improvements in managing cancer (CRPC). Although de novo cases are rare, NEPC as a special CRPC adenocarcinomas in recent years (1), it is expected that the subtype of CRPC (10%–20%) is thought to occur via NE incidence of NEPC will further increase (2). Unfortunately, the transdifferentiation of prostate adenocarcinomas in response to overall median survival of NEPC with small cell feature is less than 1 year, primarily due to its aggressiveness and limited available treatment options (3). As such, a better understanding of the 1 Vancouver Prostate Centre, Department of Urologic Sciences, University of mechanisms underlying NEPC development remains much need- 2 British Columbia, Vancouver, British Columbia, Canada. Department of Exper- ed in order to develop more effective therapeutics. imental Therapeutics, BC Cancer, Vancouver, British Columbia, Canada. 3School of Life, Health and Chemical Sciences, the Open University, Walton Hall, Milton One distinct histopathologic feature of NEPC cells is the fi Keynes, United Kingdom. 4Princess Margaret Cancer Center, University Health frequent manifestation of hyperchromatic nuclei with nely Network, Toronto, Ontario, Canada. 5Department of Medical Biophysics, Uni- dispersed chromatin and inconspicuous nucleoli, a phenomenon versity of Toronto, Toronto, Ontario, Canada. known as "salt and pepper" chromatin (3, 4). This is in contrast to Note: Supplementary data for this article are available at Cancer Research prostatic adenocarcinoma cells, which tend to have enlarged Online (http://cancerres.aacrjournals.org/). nuclei with prominent nucleoli (3, 4). This special hematoxy- Corresponding Authors: Yuzhuo Wang, BC Cancer/UBC, 675 W. 10th Ave., lin-staining characteristic suggests a distinct NEPC heterochro- Vancouver, BC V5Z 1L3, Canada. Phone: 604-675-8013; Fax: 604-675-8019; matin pattern (5). Heterochromatin is a condensed and transcrip- E-mail: [email protected]; and Dong Lin, Vancouver Prostate Centre, Department tionally inert chromosome conformation, regulated epigenetical- of Urologic Sciences, University of British Columbia, 2660 Oak Street, Vancou- ly by precise and dedicated machineries (6, 7). Multiple studies ver, BC V6H3Z6, Canada. Phone: 604-675-7013; E-mail: have demonstrated that epigenetic regulation is one major mech- [email protected] anism underlying NEPC development (8–11). An NEPC-specific doi: 10.1158/0008-5472.CAN-17-3677 heterochromatin gene signature can thus shed light to help better Ó2018 American Association for Cancer Research. understand the disease and identify novel therapeutic targets. www.aacrjournals.org 2691 Downloaded from cancerres.aacrjournals.org on October 2, 2021. © 2018 American Association for Cancer Research. Published OnlineFirst February 27, 2018; DOI: 10.1158/0008-5472.CAN-17-3677 Ci et al. Another distinct feature of NEPC is the lineage alteration associ- the repressive histone mark H3K9me3 on their respective gene ated with a decrease or loss of crucial adenocarcinoma lineage– promoters. specific transcription factors, AR, FOXA1, and RE1 silencing transcription factor (REST; refs. 12–15). This leads to the repres- Materials and Methods sion of AR-regulated genes (e.g., prostate-specific antigen) and PDXs and clinical datasets gain of neuroendocrine markers [e.g., neural cell adhesion mol- All PDX tumor lines were grafted in NSG mice as previously ecule 1 (NCAM1/CD56), neuronal-specific enolase (NSE), chro- described (21). This study followed the ethical guidelines stated in mogranin A (CHGA), and synaptophysin (SYP); refs. 12–14]. the Declaration of Helsinki, specimens were obtained from Although recent studies have reported that EZH2 inhibits AR patients with their informed written consent form following a expression and SRRM4 mediates REST splicing (8, 16), the protocol (#H09-01628) approved by the Institutional Review mechanisms underlying the constant transcriptional suppression Board of the University of British Columbia (UBC). Animal of AR and REST are still not well understood. studies under the protocol #A17-0165 were approved by UBC's One of the major hurdles in studying NEPC is the lack of Animal Care and Use Committee. The LTL331-castrated tissues clinically relevant models, but substantial progress has been and the NEPC-relapsed LTL331R tissues were harvested at differ- achieved recently in modeling NEPC development. Employing ent time points after host castration (22). Transcriptomic analysis genetically engineered mouse (GEM) models, ectopic gain of for all PDX tumors, with the exception of the LTL331-331R N-myc and concomitant loss of Rb1 and Tp53 have both been castration time-series samples, was performed using GE 8 Â demonstrated to induce de novo NEPC-like tumors (8, 11). 60K microarray, as previously described (21). Transcriptomic Potent ADT using abiraterone in Tp53/Pten double-deficient analysis of the LTL331-331R time-series was done using RNA- mice also promoted overt neuroendocrine (NE) transdifferen- sequencing data. tiation of luminal adenocarcinoma cells (17). Employing engi- The clinical cohorts used in this study are as follows: RNA-seq neered human primary cells or cell lines, N-myc, SOX2, BRN2, data for the Beltran 2011 cohort (30 adenocarcinomas, 6 NEPCs) and SRRM4 overexpressions have all been shown to promote were from Weill Medical College of Cornell University (23); RNA- NE transdifferentiation (16, 18–20). Our laboratory has estab- Seq data for the Beltran 2016 (34 CRPC-adenocarcinomas, 15 lished over 45 high-fidelity patient-derived xenograft (PDX) NEPCs), the Grasso 2012 (31 CRPC-adenocarcinomas, 4 NEPCs), models of prostate cancer including 5 NEPCs (www.livingtu and the Kumar 2016 (149 CRPC-adenocarcinomas) cohorts were morlab.ca; ref. 21). Among them, LTL331/331R is the first-in- accessed through the cBioPortal (10, 24–26); microarray data for field and unique PDX model of adenocarcinoma-to-NEPC the Varambally 2015 cohort (6 CRPC-adenocarcinomas) were transdifferentiation. Upon host castration, the primary adeno- accessed from the Gene Expression Omnibus (GEO) database carcinoma (LTL331) initially regresses but relapses within few (GSE3325; ref. 27). months as typical NEPC (LTL331R; ref. 21). Importantly, the whole transdifferentiation process observed in the LTL331/ 331R model is predictive of disease progression and is fully Human prostate cancer specimens recapitulated in the donor patient (22), suggesting a strong Prostate cancer specimens were obtained from the Vancouver clinical relevance. Because an overwhelming majority of Prostate Centre Tissue Bank (103 hormone-na€ve primary pros- expressed genes accounting for the NE phenotype in terminal tatic adenocarcinomas, 120 CRPC-adenocarcinomas,
Recommended publications
  • 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]
  • Harnessing Gene Expression Profiles for the Identification of Ex Vivo Drug
    cancers Article Harnessing Gene Expression Profiles for the Identification of Ex Vivo Drug Response Genes in Pediatric Acute Myeloid Leukemia David G.J. Cucchi 1 , Costa Bachas 1 , Marry M. van den Heuvel-Eibrink 2,3, Susan T.C.J.M. Arentsen-Peters 3, Zinia J. Kwidama 1, Gerrit J. Schuurhuis 1, Yehuda G. Assaraf 4, Valérie de Haas 3 , Gertjan J.L. Kaspers 3,5 and Jacqueline Cloos 1,* 1 Hematology, Cancer Center Amsterdam, Amsterdam UMC, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands; [email protected] (D.G.J.C.); [email protected] (C.B.); [email protected] (Z.J.K.); [email protected] (G.J.S.) 2 Department of Pediatric Oncology/Hematology, Erasmus MC–Sophia Children’s Hospital, 3015 CN Rotterdam, The Netherlands; [email protected] 3 Princess Máxima Center for Pediatric Oncology, 3584 CS Utrecht, The Netherlands; [email protected] (S.T.C.J.M.A.-P.); [email protected] (V.d.H.); [email protected] (G.J.L.K.) 4 The Fred Wyszkowski Cancer Research, Laboratory, Department of Biology, Technion-Israel Institute of Technology, 3200003 Haifa, Israel; [email protected] 5 Emma’s Children’s Hospital, Amsterdam UMC, Vrije Universiteit Amsterdam, Pediatric Oncology, 1081 HV Amsterdam, The Netherlands * Correspondence: [email protected] Received: 21 April 2020; Accepted: 12 May 2020; Published: 15 May 2020 Abstract: Novel treatment strategies are of paramount importance to improve clinical outcomes in pediatric AML. Since chemotherapy is likely to remain the cornerstone of curative treatment of AML, insights in the molecular mechanisms that determine its cytotoxic effects could aid further treatment optimization.
    [Show full text]
  • The Roles of Histone Deacetylase 5 and the Histone Methyltransferase Adaptor WDR5 in Myc Oncogenesis
    The Roles of Histone Deacetylase 5 and the Histone Methyltransferase Adaptor WDR5 in Myc oncogenesis By Yuting Sun This thesis is submitted in fulfilment of the requirements for the degree of Doctor of Philosophy at the University of New South Wales Children’s Cancer Institute Australia for Medical Research School of Women’s and Children’s Health, Faculty of Medicine University of New South Wales Australia August 2014 PLEASE TYPE THE UNIVERSITY OF NEW SOUTH WALES Thesis/Dissertation Sheet Surname or Family name: Sun First name: Yuting Other name/s: Abbreviation for degree as given in the University calendar: PhD School : School of·Women's and Children's Health Faculty: Faculty of Medicine Title: The Roles of Histone Deacetylase 5 and the Histone Methyltransferase Adaptor WDR5 in Myc oncogenesis. Abstract 350 words maximum: (PLEASE TYPE) N-Myc Induces neuroblastoma by regulating the expression of target genes and proteins, and N-Myc protein is degraded by Fbxw7 and NEDD4 and stabilized by Aurora A. The class lla histone deacetylase HDAC5 suppresses gene transcription, and blocks myoblast and leukaemia cell differentiation. While histone H3 lysine 4 (H3K4) trimethylation at target gene promoters is a pre-requisite for Myc· induced transcriptional activation, WDRS, as a histone H3K4 methyltransferase presenter, is required for H3K4 methylation and transcriptional activation mediated by a histone H3K4 methyltransferase complex. Here, I investigated the roles of HDAC5 and WDR5 in N-Myc overexpressing neuroblastoma. I have found that N-Myc upregulates HDAC5 protein expression, and that HDAC5 represses NEDD4 gene expression, increases Aurora A gene expression and consequently upregulates N-Myc protein expression in neuroblastoma cells.
    [Show full text]
  • (RRP1B) Modulates Metastasis Through Regulation of Histone Methylation
    Published OnlineFirst August 4, 2014; DOI: 10.1158/1541-7786.MCR-14-0167 Molecular Cancer Oncogenes and Tumor Suppressors Research Metastasis-Associated Protein Ribosomal RNA Processing 1 Homolog B (RRP1B) Modulates Metastasis through Regulation of Histone Methylation Minnkyong Lee1, Amy M. Dworkin1, Jens Lichtenberg1, Shashank J. Patel1, Niraj S. Trivedi2, Derek Gildea2, David M. Bodine1, and Nigel P.S. Crawford1 Abstract Overexpression of ribosomal RNA processing 1 homolog B (RRP1B) induces a transcriptional profile that accurately predicts patient outcome in breast cancer. However, the mechanism by which RRP1B modulates transcription is unclear. Here, the chromatin-binding properties of RRP1B were examined to define how it regulates metastasis-associated transcription. To identify genome-wide RRP1B-binding sites, high-throughput ChIP-seq was performed in the human breast cancer cell line MDA-MB-231 and HeLa cells using antibodies against endogenous RRP1B. Global changes in repressive marks such as histone H3 lysine 9 trimethylation (H3K9me3) were also examined by ChIP-seq. Analysis of these samples identified 339 binding regions in MDA- MB-231 cells and 689 RRP1B-binding regions in HeLa cells. Among these, 136 regions were common to both cell lines. Gene expression analyses of these RRP1B-binding regions revealed that transcriptional repression is the primary result of RRP1B binding to chromatin. ChIP-reChIP assays demonstrated that RRP1B co-occupies loci with decreased gene expression with the heterochromatin-associated proteins, tripartite motif-containing protein 28 (TRIM28/KAP1), and heterochromatin protein 1-a (CBX5/HP1a). RRP1B occupancy at these loci was also associated with higher H3K9me3 levels, indicative of heterochromatinization mediated by the TRIM28/HP1a complex.
    [Show full text]
  • UNIVERSITY of CALIFORNIA, IRVINE Combinatorial Regulation By
    UNIVERSITY OF CALIFORNIA, IRVINE Combinatorial regulation by maternal transcription factors during activation of the endoderm gene regulatory network DISSERTATION submitted in partial satisfaction of the requirements for the degree of DOCTOR OF PHILOSOPHY in Biological Sciences by Kitt D. Paraiso Dissertation Committee: Professor Ken W.Y. Cho, Chair Associate Professor Olivier Cinquin Professor Thomas Schilling 2018 Chapter 4 © 2017 Elsevier Ltd. © 2018 Kitt D. Paraiso DEDICATION To the incredibly intelligent and talented people, who in one way or another, helped complete this thesis. ii TABLE OF CONTENTS Page LIST OF FIGURES vii LIST OF TABLES ix LIST OF ABBREVIATIONS X ACKNOWLEDGEMENTS xi CURRICULUM VITAE xii ABSTRACT OF THE DISSERTATION xiv CHAPTER 1: Maternal transcription factors during early endoderm formation in 1 Xenopus Transcription factors co-regulate in a cell type-specific manner 2 Otx1 is expressed in a variety of cell lineages 4 Maternal otx1 in the endodermal conteXt 5 Establishment of enhancers by maternal transcription factors 9 Uncovering the endodermal gene regulatory network 12 Zygotic genome activation and temporal control of gene eXpression 14 The role of maternal transcription factors in early development 18 References 19 CHAPTER 2: Assembly of maternal transcription factors initiates the emergence 26 of tissue-specific zygotic cis-regulatory regions Introduction 28 Identification of maternal vegetally-localized transcription factors 31 Vegt and OtX1 combinatorially regulate the endodermal 33 transcriptome iii
    [Show full text]
  • Integration of 198 Chip-Seq Datasets Reveals Human Cis-Regulatory Regions
    Integration of 198 ChIP-seq datasets reveals human cis-regulatory regions Hamid Bolouri & Larry Ruzzo Thanks to Stephen Tapscott, Steven Henikoff & Zizhen Yao Slides will be available from: http://labs.fhcrc.org/bolouri/ Email [email protected] for manuscript (Bolouri & Ruzzo, submitted) Kleinjan & van Heyningen, Am. J. Hum. Genet., 2005, (76)8–32 Epstein, Briefings in Func. Genom. & Protemoics, 2009, 8(4)310-16 Regulation of SPi1 (Sfpi1, PU.1 protein) expression – part 1 miR155*, miR569# ~750nt promoter ~250nt promoter The antisense RNA • causes translational stalling • has its own promoter • requires distal SPI1 enhancer • is transcribed with/without SPI1. # Hikami et al, Arthritis & Rheumatism, 2011, 63(3):755–763 * Vigorito et al, 2007, Immunity 27, 847–859 Ebralidze et al, Genes & Development, 2008, 22: 2085-2092. Regulation of SPi1 expression – part 2 (mouse coordinates) Bidirectional ncRNA transcription proportional to PU.1 expression PU.1/ELF1/FLI1/GLI1 GATA1 GATA1 Sox4/TCF/LEF PU.1 RUNX1 SP1 RUNX1 RUNX1 SP1 ELF1 NF-kB SATB1 IKAROS PU.1 cJun/CEBP OCT1 cJun/CEBP 500b 500b 500b 500b 500b 750b 500b -18Kb -14Kb -12Kb -10Kb -9Kb Chou et al, Blood, 2009, 114: 983-994 Hoogenkamp et al, Molecular & Cellular Biology, 2007, 27(21):7425-7438 Zarnegar & Rothenberg, 2010, Mol. & cell Biol. 4922-4939 An NF-kB binding-site variant in the SPI1 URE reduces PU.1 expression & is GGGCCTCCCC correlated with AML GGGTCTTCCC Bonadies et al, Oncogene, 2009, 29(7):1062-72. SATB1 binding site A distal single nucleotide polymorphism alters long- range regulation of the PU.1 gene in acute myeloid leukemia Steidl et al, J Clin Invest.
    [Show full text]
  • Lamina-Associated Domains: Peripheral Matters and Internal Affairs Nolwenn Briand1,2 and Philippe Collas1,2*
    Briand and Collas Genome Biology (2020) 21:85 https://doi.org/10.1186/s13059-020-02003-5 REVIEW Open Access Lamina-associated domains: peripheral matters and internal affairs Nolwenn Briand1,2 and Philippe Collas1,2* * Correspondence: philippe.collas@ medisin.uio.no Abstract 1Department of Molecular Medicine, Faculty of Medicine, University of At the nuclear periphery, associations of chromatin with the nuclear lamina through Oslo, PO Box 1112 Blindern, 0317 lamina-associated domains (LADs) aid functional organization of the genome. We Oslo, Norway review the organization of LADs and provide evidence of LAD heterogeneity from 2Department of Immunology and Transfusion Medicine, Oslo cell ensemble and single-cell data. LADs are typically repressive environments in the University Hospital, 0424 Oslo, genome; nonetheless, we discuss findings of lamin interactions with regulatory Norway elements of active genes, and the role lamins may play in genome regulation. We address the relationship between LADs and other genome organizers, and the involvement of LADs in laminopathies. The current data lay the basis for future studies on the significance of lamin-chromatin interactions in health and disease. Keywords: 3D genome, Chromatin, LAD, Lamin A mutation, Nuclear envelope, Nuclear lamin, Radial positioning Introduction A snapshot of the genome in three dimensions In the eukaryotic nucleus, individual chromosomes occupy non-random territories and a radial (i.e., center-to-periphery) disposition enabled by sequestering specific chromo- somal domains to the nuclear envelope. Current views of the three-dimensional (3D) genome depict a hierarchical architecture of chromatin [1, 2]. Chromosome conform- ation capture combined with high-throughput sequencing (a method called Hi-C [3]) makes it possible to map interactions within and between chromosomes across ge- nomes [4].
    [Show full text]
  • Prognostic Values of Distinct CBX Family Members in Breast Cancer Liang, Yuan-Ke; Lin, Hao-Yu; Chen, Chun-Fa; Zeng, De
    University of Groningen Prognostic values of distinct CBX family members in breast cancer Liang, Yuan-Ke; Lin, Hao-Yu; Chen, Chun-Fa; Zeng, De Published in: Oncotarget DOI: 10.18632/oncotarget.21325 IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check the document version below. Document Version Publisher's PDF, also known as Version of record Publication date: 2017 Link to publication in University of Groningen/UMCG research database Citation for published version (APA): Liang, Y-K., Lin, H-Y., Chen, C-F., & Zeng, D. (2017). Prognostic values of distinct CBX family members in breast cancer. Oncotarget, 8(54), 92375-92387. https://doi.org/10.18632/oncotarget.21325 Copyright Other than for strictly personal use, it is not permitted to download or to forward/distribute the text or part of it without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license (like Creative Commons). The publication may also be distributed here under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license. More information can be found on the University of Groningen website: https://www.rug.nl/library/open-access/self-archiving-pure/taverne- amendment. Take-down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. Downloaded from the University of Groningen/UMCG research database (Pure): http://www.rug.nl/research/portal.
    [Show full text]
  • Virtual Chip-Seq: Predicting Transcription Factor Binding
    bioRxiv preprint doi: https://doi.org/10.1101/168419; this version posted March 12, 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 Virtual ChIP-seq: predicting transcription factor binding 2 by learning from the transcriptome 1,2,3 1,2,3,4,5 3 Mehran Karimzadeh and Michael M. Hoffman 1 4 Department of Medical Biophysics, University of Toronto, Toronto, ON, Canada 2 5 Princess Margaret Cancer Centre, Toronto, ON, Canada 3 6 Vector Institute, Toronto, ON, Canada 4 7 Department of Computer Science, University of Toronto, Toronto, ON, Canada 5 8 Lead contact: michael.hoff[email protected] 9 March 8, 2019 10 Abstract 11 Motivation: 12 Identifying transcription factor binding sites is the first step in pinpointing non-coding mutations 13 that disrupt the regulatory function of transcription factors and promote disease. ChIP-seq is 14 the most common method for identifying binding sites, but performing it on patient samples is 15 hampered by the amount of available biological material and the cost of the experiment. Existing 16 methods for computational prediction of regulatory elements primarily predict binding in genomic 17 regions with sequence similarity to known transcription factor sequence preferences. This has limited 18 efficacy since most binding sites do not resemble known transcription factor sequence motifs, and 19 many transcription factors are not even sequence-specific. 20 Results: 21 We developed Virtual ChIP-seq, which predicts binding of individual transcription factors in new 22 cell types using an artificial neural network that integrates ChIP-seq results from other cell types 23 and chromatin accessibility data in the new cell type.
    [Show full text]
  • Role of Estrogen Receptor in Breast Cancer Cell Gene Expression
    4046 MOLECULAR MEDICINE REPORTS 13: 4046-4050, 2016 Role of estrogen receptor in breast cancer cell gene expression YABING ZHENG1, XIYING SHAO1, YUAN HUANG1, LEI SHI1, BO CHEN2, XIAOJIA WANG1, HONGJIAN YANG3, ZHANHONG CHEN1 and XIPING ZHANG3 Departments of 1Medical Oncology (Breast), 2Pathology and 3Breast Surgery, Zhejiang Cancer Hospital, Hangzhou, Zhejiang 310022, P.R. China Received April 28, 2015; Accepted February 23, 2016 DOI: 10.3892/mmr.2016.5018 Abstract. The aim of the present study was to establish the Europe in 2012, and the number of breast cancer-associated underlying regulatory mechanism of estrogen receptor (ER) mortalities is 131,000 (6). Furthermore, breast cancer is in breast cancer cell gene expression. A gene expression the most common cause of cancer-associated mortality in profile accession( no. GSE11324) was downloaded from the females. Therefore, it is essential to understand its molecular Gene Expression Omnibus (GEO) database. Differentially mechanism and develop more effective therapeutic methods expressed genes (DEGs) from an estrogen treatment group and for breast cancer treatment. a control group were identified. Chromatin immunoprecipita- The estrogen receptor (ER) is critical in determining the tion with high-throughput sequencing data (series GSE25710) phenotype of human breast cancers and is one of the most was obtained from the GEO for the ER binding sites, and important therapeutic targets (7). Furthermore, certain studies binding and expression target analysis was performed. A total have suggested that activation of ER is responsible for various of 3,122 DEGs were obtained and ER was demonstrated to biological processes, including cell growth and differentia- exhibit inhibition and activation roles during the regulation tion, and programmed cell death (8,9).
    [Show full text]
  • Lamin B1 Is Required for Mature Neuron-Specific Gene Expression
    ARTICLE Received 4 Jul 2016 | Accepted 28 Feb 2017 | Published 20 Apr 2017 DOI: 10.1038/ncomms15098 OPEN Lamin B1 is required for mature neuron-specific gene expression during olfactory sensory neuron differentiation Crystal M. Gigante1,2, Michele Dibattista3,4, Frederick N. Dong1, Xiaobin Zheng2, Sibiao Yue2, Stephen G. Young5, Johannes Reisert3, Yixian Zheng2 & Haiqing Zhao1 B-type lamins are major constituents of the nuclear lamina in all metazoan cells, yet have specific roles in the development of certain cell types. Although they are speculated to regulate gene expression in developmental contexts, a direct link between B-type lamins and developmental gene expression in an in vivo system is currently lacking. Here, we identify lamin B1 as a key regulator of gene expression required for the formation of functional olfactory sensory neurons. By using targeted knockout in olfactory epithelial stem cells in adult mice, we show that lamin B1 deficient neurons exhibit attenuated response to odour stimulation. This deficit can be explained by decreased expression of genes involved in mature neuron function, along with increased expression of genes atypical of the olfactory lineage. These results support that the broadly expressed lamin B1 regulates expression of a subset of genes involved in the differentiation of a specific cell type. 1 Department of Biology, The Johns Hopkins University, Baltimore, Maryland 21218, USA. 2 Department of Embryology, Carnegie Institution for Science, Baltimore, Maryland 21218, USA. 3 Monell Chemical Senses Center, Philadelphia, Pennsylvania 19104, USA. 4 Department of Basic Medical Sciences, Neuroscience and Sensory Organs, University of Bari ‘A. Moro’, Bari 70121, Italy. 5 Department of Medicine, Molecular Biology Institute and Department of Human Genetics, University of California, Los Angeles, California 90095, USA.
    [Show full text]
  • Integrated Computational Approach to the Analysis of RNA-Seq Data Reveals New Transcriptional Regulators of Psoriasis
    OPEN Experimental & Molecular Medicine (2016) 48, e268; doi:10.1038/emm.2016.97 & 2016 KSBMB. All rights reserved 2092-6413/16 www.nature.com/emm ORIGINAL ARTICLE Integrated computational approach to the analysis of RNA-seq data reveals new transcriptional regulators of psoriasis Alena Zolotarenko1, Evgeny Chekalin1, Alexandre Mesentsev1, Ludmila Kiseleva2, Elena Gribanova2, Rohini Mehta3, Ancha Baranova3,4,5,6, Tatiana V Tatarinova6,7,8, Eleonora S Piruzian1 and Sergey Bruskin1,5 Psoriasis is a common inflammatory skin disease with complex etiology and chronic progression. To provide novel insights into the regulatory molecular mechanisms of the disease, we performed RNA sequencing analysis of 14 pairs of skin samples collected from patients with psoriasis. Subsequent pathway analysis and extraction of the transcriptional regulators governing psoriasis-associated pathways was executed using a combination of the MetaCore Interactome enrichment tool and the cisExpress algorithm, followed by comparison to a set of previously described psoriasis response elements. A comparative approach allowed us to identify 42 core transcriptional regulators of the disease associated with inflammation (NFκB, IRF9, JUN, FOS, SRF), the activity of T cells in psoriatic lesions (STAT6, FOXP3, NFATC2, GATA3, TCF7, RUNX1), the hyper- proliferation and migration of keratinocytes (JUN, FOS, NFIB, TFAP2A, TFAP2C) and lipid metabolism (TFAP2, RARA, VDR). In addition to the core regulators, we identified 38 transcription factors previously not associated with the disease that can clarify the pathogenesis of psoriasis. To illustrate these findings, we analyzed the regulatory role of one of the identified transcription factors (TFs), FOXA1. Using ChIP-seq and RNA-seq data, we concluded that the atypical expression of the FOXA1 TF is an important player in the disease as it inhibits the maturation of naive T cells into the (CD4+FOXA1+CD47+CD69+PD-L1(hi) FOXP3 − ) regulatory T cell subpopulation, therefore contributing to the development of psoriatic skin lesions.
    [Show full text]