Spatio-Temporal Gene Discovery for Autism Spectrum Disorder

Total Page:16

File Type:pdf, Size:1020Kb

Spatio-Temporal Gene Discovery for Autism Spectrum Disorder bioRxiv preprint doi: https://doi.org/10.1101/256693; this version posted March 8, 2018. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Spatio-Temporal Gene Discovery for Autism Spectrum Disorder Utku Norman∗ and A. Ercument Cicek∗;y ∗Computer Engineering Department, Bilkent University, Ankara, 06800, Turkey yComputational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA, 15213, USA Correspondence: [email protected] Abstract—Whole exome sequencing (WES) studies for Autism full picture of the genetic architecture will take a long time Spectrum Disorder (ASD) could identify only around six dozen and will be financially costly. Guilt-by-association based gene risk genes to date because the genetic architecture of the discovery techniques came handy to predict ASD risk genes. disorder is highly complex. To speed the gene discovery process up, a few network-based ASD gene discovery algorithms were They assume these genes are working as a functional cluster. proposed. Although these methods use static gene interaction Starting from already known risk genes, these techniques networks, functional clustering of genes is bound to evolve predict a cluster of closely interacting genes. There are only a during neurodevelopment and disruptions are likely to have a few network-based ASD-tailored gene discovery algorithms cascading effect on the future associations. Thus, approaches that in the literature (Gilman et al., 2011; Liu et al., 2014; disregard the dynamic nature of neurodevelopment are limited in power. Here, we present a spatio-temporal gene discovery Hormozdiari et al., 2014) that are later improved (Gilman algorithm for ASD, which leverages information from evolving et al., 2012; Liu et al., 2015) and these algorithms are de- gene coexpression networks of neurodevelopment. The algorithm scribed in the Background section. Despite their fundamental solves a variant of prize-collecting Steiner forest-based problem differences, all of these methods have one point in common: on coexpression networks to model neurodevelopment and trans- The biological gene interaction networks that they use are fer information from precursor neurodevelopmental windows. The decisions made by the algorithm can be traced back, static. None considers the fact that gene interactions (coex- adding interpretability to the results. We apply the algorithm pressions) evolve over time: It was clearly demonstrated that on WES data of 3,871 samples and identify risk clusters using different neurodevelopmental spatio-temporal windows have BrainSpan coexpression networks of early- and mid-fetal periods. different topologies and consequently, the clustering of ASD On an independent dataset, we show that incorporation of the susceptibility genes drastically change (Willsey et al., 2013). temporal dimension increases the prediction power: Predicted clusters are hit more and show higher enrichment in ASD-related Moreover, dysregulation of pathways in earlier periods have functions compared to the state-of-the-art. Code is available at cascading effects on the circuitry of the future time periods. http://ciceklab.cs.bilkent.edu.tr/ST-Steiner/. For instance, it is shown that ASD risk increases 36-fold with a Index Terms—Gene Discovery, Spatio-Temporal Networks, damage in early cerebellar circuitry, but a damage in adulthood Autism Spectrum Disorder, prize-collecting Steiner tree problem. does not cause dysfunction in social interactions (Wang et al., 2014). Furthermore, early cerebellar damage has been shown to result in a cascade of long-term deficits in cerebellum and I. INTRODUCTION cause ASD (de la Torre-Ubieta et al., 2016). Thus, we argue Autism Spectrum Disorder (ASD) is a common neurode- that the state-of-the-art methods are limited in their prediction velopmental disorder that affects ∼1.5% of the children in the power since static networks they use would fail to capture the US (Autism and Investigators, 2014). Recent whole exome se- dynamic nature of neurodevelopment. quencing (WES) efforts have paved the way for identification In this article, we propose a novel ASD gene discovery of dozens of ASD risk genes (O’Roak et al., 2012; Iossifov algorithm termed ST-Steiner. The algorithm modifies the prize- et al., 2012; Neale et al., 2012; Sanders et al., 2012; De Rubeis collecting Steiner forest (PCSF) problem and extends it to et al., 2014; Iossifov et al., 2014; Sanders et al., 2015). Unfor- spatio-temporal networks in order to mimic neurodevelopment. tunately, this number corresponds to only a small portion of Instead of performing gene discovery on a single network the large genetic puzzle, which is expected to contain around or any number of networks separately, ST-Steiner solves an a thousand genes (He et al., 2013). Detection of de novo optimization problem progressively over a cascade of spatio- loss-of-function mutations (dnLoF) has been the key for gene temporal networks while leveraging information coming from discovery due to their high signal-to-noise ratio. However, earlier neurodevelopmental periods. The algorithm has three such mutations are rare and they affect a diverse set of genes. novel aspects: (i) For the first time, the problem is solved Thus, for most of the genes, the rarity and diversity of variants on spatio-temporal coexpression networks and the dynamic prevent statistically significant signals from being observed. nature of neurodevelopment is taken into account; (ii) The Therefore, our yield of discovered genes is still low, despite results are more interpretable compared to other methods analyzing thousands of trios. The journey towards getting the in the literature, since the decisions made by the algorithm bioRxiv preprint doi: https://doi.org/10.1101/256693; this version posted March 8, 2018. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. can be traced back on the spatio-temporal cascade; (iii) As A more recent method is DAWN, a hidden Markov random the problem is formulated as a PCSF prediction problem, field based approach. It predicts ASD genes by assigning the algorithm predicts only the genes that are essential for posterior risk scores to each gene, based on their interactions connectivity of known risk genes. Note that, this is in contrast with other genes in a gene coexpression network (Liu et al., to the other approaches which can return redundant paths 2014, 2015). Inputs to DAWN are prior TADA risk scores per between known risk genes. We apply the algorithm on exome gene (He et al., 2013), and a partial coexpression network. A sequencing data from (De Rubeis et al., 2014), and identify gene is assigned a higher risk score if the gene has a high gene clusters using 2 gene coexpression network cascades prior risk, or is highly interacting with other risk genes. of early-fetal and mid-fetal periods. For the first time, we Latest is the MAGI algorithm (Hormozdiari et al., 2014). benchmark the performances of the network-based ASD gene MAGI also uses a seed-and-extend based approach like NET- discovery algorithms on common training data. We validate BAG. Using a combination of PPI networks and a coexpres- the predicted clusters with independent exome sequencing sion network, it generates seed pathways which are enriched data from (Iossifov et al., 2014) and find that ST-Steiner’s with dnLoFs in cases, compared to controls. Then, it merges predictions (i) identify more risk genes with less predictions the pathways as long as the cluster score is improved. Finally, compared to the state-of-the-art methods, (ii) overlap more some genes are removed, added or swapped in order to with targets of known ASD-related transcription factors and reach a local maximum score. Using a permutation test, the (iii) enriched more in ASD-related KEGG pathways. We significance of the detected modules are assessed. In addition, also show in various controlled settings that using temporal MAGI finds an extended module which is the union of high information boosts the prediction power, which supports our scoring suboptimal clusters that overlap with the best module. claim that the clustering of risk genes is spatio-temporal. Besides, in order to obtain additional modules, genes of the We demonstrate that prediction of genes like NOTCH2 and best cluster are removed and MAGI is run again to obtain a NOTCH3 can be traced back to their role in neuron differ- new cluster. entiation in embryonic period. Finally, we predict new genes A final method takes a different path: Even though it uses a with the incorporation of time dimension which point to a brain specific functional gene interaction network, in contrast disruption of organization of microtubules in the embryonic to ST-Steiner and the three methods mentioned above, it does period. not aim to predict a cluster of genes (Krishnan et al., 2016). Instead, it uses the known ASD genes and their connectivity II. BACKGROUND patterns as features to train an SVM classifier, which is used to assign every gene a probability of being associated to ASD. 1) Gene discovery in the context of neurodevelopmental 2) No consensus on the network choice: Module detection disorders: Network-based gene discovery and disease gene methods in the literature assume that ASD risk genes must prioritization have been very popular research areas in the last be clustered in a gene interaction network. Yet, the nature of decade (Bebek et al., 2012; Bromberg, 2013; Cowen et al., the networks they prefer differs. NETBAG constructs its own 2017). While these methods proved their use mainly in the network through a comprehensive analysis of many annotation context of cancer, a parallel track has opened up for network- resources (e.g.
Recommended publications
  • Condensation of Prometaphase Chromo-Somes
    Condensation of Prometaphase Chromo- somes D'Eustachio, P., Matthews, L., Orlic-Milacic, M. European Bioinformatics Institute, New York University Langone Medical Center, Ontario Institute for Cancer Research, Oregon Health and Science University. The contents of this document may be freely copied and distributed in any media, provided the authors, plus the institutions, are credited, as stated under the terms of Creative Commons Attribution 4.0 Inter- national (CC BY 4.0) License. For more information see our license. 09/03/2019 Introduction Reactome is open-source, open access, manually curated and peer-reviewed pathway database. Pathway annotations are authored by expert biologists, in collaboration with Reactome editorial staff and cross- referenced to many bioinformatics databases. A system of evidence tracking ensures that all assertions are backed up by the primary literature. Reactome is used by clinicians, geneticists, genomics research- ers, and molecular biologists to interpret the results of high-throughput experimental studies, by bioin- formaticians seeking to develop novel algorithms for mining knowledge from genomic studies, and by systems biologists building predictive models of normal and disease variant pathways. The development of Reactome is supported by grants from the US National Institutes of Health (P41 HG003751), University of Toronto (CFREF Medicine by Design), European Union (EU STRP, EMI-CD), and the European Molecular Biology Laboratory (EBI Industry program). Literature references Fabregat, A., Sidiropoulos, K., Viteri, G., Forner, O., Marin-Garcia, P., Arnau, V. et al. (2017). Reactome pathway ana- lysis: a high-performance in-memory approach. BMC bioinformatics, 18, 142. ↗ Sidiropoulos, K., Viteri, G., Sevilla, C., Jupe, S., Webber, M., Orlic-Milacic, M.
    [Show full text]
  • Condensins, Chromatin Remodeling and Gene Transcription
    Chapter 4 Condensins, Chromatin Remodeling and Gene Transcription Laurence O. W. Wilson and Aude M. Fahrer Additional information is available at the end of the chapter http://dx.doi.org/10.5772/55732 1. Introduction Condensin complexes condense chromosomes during mitosis, turning the diffuse interphase chromosomes into the familiar X-shaped compact chromosomes that segregate during cell division. More recently a second role for condensins has emerged, in the epigenetic regulation of interphase gene transcription. This second fascinating role is very difficult to study since defects in condensin will usually interfere with mitosis and result in cell death. While several excellent reviews of condensin function have recently been published (see [1-3]), in this article we concentrate on the epigenetic functions of condensins and how they can be studied. We also provide the first summary of condensin protein splice forms, a largely overlooked contributor to condensin variation. The DNA within a cell is too large to fit inside if left in its unwound state. In order to accom‐ modate the genetic material, a cell packages this DNA as chromati:, a combination of DNA bound to protein. The DNA is wound around protein complexes known as histones to form nucleosomes, which form a “beads-on-a-string” structure. These nucleosomes can be com‐ pacted further to produce a highly condensed structure that fits inside the nucleus of a cell. The regulation of this structure is vital for cell growth and survival. During mitosis, chromatin must be unwound so that it can be properly replicated and then repackaged into sister chromosomes that must then be segregated into the dividing cells.
    [Show full text]
  • Sleeping Beauty Transposon Mutagenesis Identifies Genes That
    Sleeping Beauty transposon mutagenesis identifies PNAS PLUS genes that cooperate with mutant Smad4 in gastric cancer development Haruna Takedaa,b, Alistair G. Rustc,d, Jerrold M. Warda, Christopher Chin Kuan Yewa, Nancy A. Jenkinsa,e, and Neal G. Copelanda,e,1 aDivision of Genomics and Genetics, Institute of Molecular and Cell Biology, Agency for Science, Technology and Research, Singapore 138673; bDepartment of Pathology, School of Medicine, Kanazawa Medical University, Ishikawa 920-0293, Japan; cExperimental Cancer Genetics, Wellcome Trust Sanger Institute, Cambridge CB10 1HH, United Kingdom; dTumour Profiling Unit, The Institute of Cancer Research, Chester Beatty Laboratories, London SW3 6JB, United Kingdom; and eCancer Research Program, Houston Methodist Research Institute, Houston, TX 77030 Contributed by Neal G. Copeland, February 27, 2016 (sent for review October 15, 2015; reviewed by Yoshiaki Ito and David A. Largaespada) Mutations in SMAD4 predispose to the development of gastroin- animal models that mimic human GC, researchers have infected testinal cancer, which is the third leading cause of cancer-related mice with H. pylori and then, treated them with carcinogens. They deaths. To identify genes driving gastric cancer (GC) development, have also used genetic engineering to develop a variety of trans- we performed a Sleeping Beauty (SB) transposon mutagenesis genic and KO mouse models of GC (10). Smad4 KO mice are one + − screen in the stomach of Smad4 / mutant mice. This screen iden- GC model that has been of particular interest to us (11, 12). tified 59 candidate GC trunk drivers and a much larger number of Heterozygous Smad4 KO mice develop polyps in the pyloric re- candidate GC progression genes.
    [Show full text]
  • Blocking the Short Isoform of Augmenter of Liver Regeneration Inhibits Proliferation of Human Multiple Myeloma U266 Cells Via T
    ONCOLOGY LETTERS 21: 197, 2021 Blocking the short isoform of augmenter of liver regeneration inhibits proliferation of human multiple myeloma U266 cells via the MAPK/STAT3/cell cycle signaling pathway WENQI HUANG1,2, HANG SUN1, TING HU1, DONGJU ZHU3,4, XIANLI LONG1, HUI GUO1 and QI LIU1 1Key Laboratory of Molecular Biology for Infectious Diseases (Ministry of Education), Department of Infectious Diseases, Institute for Viral Hepatitis, 2Department of Intensive Care Medicine, The Second Affiliated Hospital of Chongqing Medical University, Chongqing 400010; 3Department of Nephrology, The Affiliated Hospital of Panzhihua University, Panzhihua, Sichuan 617000; 4Department of Nephrology, The Second Affiliated Hospital of Chongqing Medical University, Chongqing 400010, P.R. China Received April 29, 2020; Accepted November 20, 2020 DOI: 10.3892/ol.2021.12458 Abstract. Multiple myeloma (MM) is the second most anti‑ALR McAb decreased cell proliferation via the MAPK, common haematological malignancy and remains an incur‑ STAT3 and cell cycle signaling pathways without increasing able disease, with most patients relapsing and requiring apoptosis. Thus, 15‑kDa‑ALR may be a new therapeutic further treatment. Augmenter of liver regeneration (ALR) is a target for myeloma. vital protein affecting fundamental processes such as energy transduction, cell survival and regeneration. Silencing ALR Introduction inhibits cell proliferation and triggers apoptosis in human MM U266 cells. However, little is known about the role Multiple myeloma (MM) accounts for ~13% of hematological of 15‑kDa‑ALR on MM. In the present study, the role of cancers, with an estimated 24,280 to 30,330 new cases and 15‑kDa‑ALR in human MM cells was investigated.
    [Show full text]
  • Human Brain Organoids Reveal Accelerated Development of Cortical Neuron Classes As a Shared Feature of Autism Risk Genes
    bioRxiv preprint doi: https://doi.org/10.1101/2020.11.10.376509; this version posted November 12, 2020. 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. Human brain organoids reveal accelerated development of cortical neuron classes as a shared feature of autism risk genes Bruna Paulsen1,2,†, Silvia Velasco1,2,†,#, Amanda J. Kedaigle1,2,3,†, Martina Pigoni1,2,†, Giorgia Quadrato4,5 Anthony Deo2,6,7,8, Xian Adiconis2,3, Ana Uzquiano1,2, Kwanho Kim1,2,3, Sean K. Simmons2,3, Kalliopi Tsafou2, Alex Albanese9, Rafaela Sartore1,2, Catherine Abbate1,2, Ashley Tucewicz1,2, Samantha Smith1,2, Kwanghun Chung9,10,11,12, Kasper Lage2,13, Aviv Regev3,14, Joshua Z. Levin2,3, Paola Arlotta1,2,# † These authors contributed equally to the work # Correspondence should be addressed to [email protected] and [email protected] 1 Department of Stem Cell and Regenerative Biology, Harvard University, Cambridge, MA 02138, USA 2 Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA 3 Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA 4 Department of Stem Cell Biology and Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, CA 90033, USA; 5 Eli and Edythe Broad CIRM Center for Regenerative Medicine and Stem Cell Research at the University of Southern California, Los Angeles, CA 90033, USA. 6 Department of Psychiatry,
    [Show full text]
  • Supplementary Methods
    Supplementary methods Human lung tissues and tissue microarray (TMA) All human tissues were obtained from the Lung Cancer Specialized Program of Research Excellence (SPORE) Tissue Bank at the M.D. Anderson Cancer Center (Houston, TX). A collection of 26 lung adenocarcinomas and 24 non-tumoral paired tissues were snap-frozen and preserved in liquid nitrogen for total RNA extraction. For each tissue sample, the percentage of malignant tissue was calculated and the cellular composition of specimens was determined by histological examination (I.I.W.) following Hematoxylin-Eosin (H&E) staining. All malignant samples retained contained more than 50% tumor cells. Specimens resected from NSCLC stages I-IV patients who had no prior chemotherapy or radiotherapy were used for TMA analysis by immunohistochemistry. Patients who had smoked at least 100 cigarettes in their lifetime were defined as smokers. Samples were fixed in formalin, embedded in paraffin, stained with H&E, and reviewed by an experienced pathologist (I.I.W.). The 413 tissue specimens collected from 283 patients included 62 normal bronchial epithelia, 61 bronchial hyperplasias (Hyp), 15 squamous metaplasias (SqM), 9 squamous dysplasias (Dys), 26 carcinomas in situ (CIS), as well as 98 squamous cell carcinomas (SCC) and 141 adenocarcinomas. Normal bronchial epithelia, hyperplasia, squamous metaplasia, dysplasia, CIS, and SCC were considered to represent different steps in the development of SCCs. All tumors and lesions were classified according to the World Health Organization (WHO) 2004 criteria. The TMAs were prepared with a manual tissue arrayer (Advanced Tissue Arrayer ATA100, Chemicon International, Temecula, CA) using 1-mm-diameter cores in triplicate for tumors and 1.5 to 2-mm cores for normal epithelial and premalignant lesions.
    [Show full text]
  • A Yeast Phenomic Model for the Influence of Warburg Metabolism on Genetic Buffering of Doxorubicin Sean M
    Santos and Hartman Cancer & Metabolism (2019) 7:9 https://doi.org/10.1186/s40170-019-0201-3 RESEARCH Open Access A yeast phenomic model for the influence of Warburg metabolism on genetic buffering of doxorubicin Sean M. Santos and John L. Hartman IV* Abstract Background: The influence of the Warburg phenomenon on chemotherapy response is unknown. Saccharomyces cerevisiae mimics the Warburg effect, repressing respiration in the presence of adequate glucose. Yeast phenomic experiments were conducted to assess potential influences of Warburg metabolism on gene-drug interaction underlying the cellular response to doxorubicin. Homologous genes from yeast phenomic and cancer pharmacogenomics data were analyzed to infer evolutionary conservation of gene-drug interaction and predict therapeutic relevance. Methods: Cell proliferation phenotypes (CPPs) of the yeast gene knockout/knockdown library were measured by quantitative high-throughput cell array phenotyping (Q-HTCP), treating with escalating doxorubicin concentrations under conditions of respiratory or glycolytic metabolism. Doxorubicin-gene interaction was quantified by departure of CPPs observed for the doxorubicin-treated mutant strain from that expected based on an interaction model. Recursive expectation-maximization clustering (REMc) and Gene Ontology (GO)-based analyses of interactions identified functional biological modules that differentially buffer or promote doxorubicin cytotoxicity with respect to Warburg metabolism. Yeast phenomic and cancer pharmacogenomics data were integrated to predict differential gene expression causally influencing doxorubicin anti-tumor efficacy. Results: Yeast compromised for genes functioning in chromatin organization, and several other cellular processes are more resistant to doxorubicin under glycolytic conditions. Thus, the Warburg transition appears to alleviate requirements for cellular functions that buffer doxorubicin cytotoxicity in a respiratory context.
    [Show full text]
  • Proteomics Analysis of Cellular Proteins Co-Immunoprecipitated with Nucleoprotein of Influenza a Virus (H7N9)
    Article Proteomics Analysis of Cellular Proteins Co-Immunoprecipitated with Nucleoprotein of Influenza A Virus (H7N9) Ningning Sun 1,:, Wanchun Sun 2,:, Shuiming Li 3, Jingbo Yang 1, Longfei Yang 1, Guihua Quan 1, Xiang Gao 1, Zijian Wang 1, Xin Cheng 1, Zehui Li 1, Qisheng Peng 2,* and Ning Liu 1,* Received: 26 August 2015 ; Accepted: 22 October 2015 ; Published: 30 October 2015 Academic Editor: David Sheehan 1 Central Laboratory, Jilin University Second Hospital, Changchun 130041, China; [email protected] (N.S.); [email protected] (J.Y.); [email protected] (L.Y.); [email protected] (G.Q.); [email protected] (X.G.); [email protected] (Z.W.); [email protected] (X.C.); [email protected] (Z.L.) 2 Key Laboratory of Zoonosis Research, Ministry of Education, Institute of Zoonosis, Jilin University, Changchun 130062, China; [email protected] 3 College of Life Sciences, Shenzhen University, Shenzhen 518057, China; [email protected] * Correspondence: [email protected] (Q.P.); [email protected] (N.L.); Tel./Fax: +86-431-8879-6510 (Q.P. & N.L.) : These authors contributed equally to this work. Abstract: Avian influenza A viruses are serious veterinary pathogens that normally circulate among avian populations, causing substantial economic impacts. Some strains of avian influenza A viruses, such as H5N1, H9N2, and recently reported H7N9, have been occasionally found to adapt to humans from other species. In order to replicate efficiently in the new host, influenza viruses have to interact with a variety of host factors. In the present study, H7N9 nucleoprotein was transfected into human HEK293T cells, followed by immunoprecipitated and analyzed by proteomics approaches.
    [Show full text]
  • Comparative Transcriptomics Reveals Similarities and Differences
    Seifert et al. BMC Cancer (2015) 15:952 DOI 10.1186/s12885-015-1939-9 RESEARCH ARTICLE Open Access Comparative transcriptomics reveals similarities and differences between astrocytoma grades Michael Seifert1,2,5*, Martin Garbe1, Betty Friedrich1,3, Michel Mittelbronn4 and Barbara Klink5,6,7 Abstract Background: Astrocytomas are the most common primary brain tumors distinguished into four histological grades. Molecular analyses of individual astrocytoma grades have revealed detailed insights into genetic, transcriptomic and epigenetic alterations. This provides an excellent basis to identify similarities and differences between astrocytoma grades. Methods: We utilized public omics data of all four astrocytoma grades focusing on pilocytic astrocytomas (PA I), diffuse astrocytomas (AS II), anaplastic astrocytomas (AS III) and glioblastomas (GBM IV) to identify similarities and differences using well-established bioinformatics and systems biology approaches. We further validated the expression and localization of Ang2 involved in angiogenesis using immunohistochemistry. Results: Our analyses show similarities and differences between astrocytoma grades at the level of individual genes, signaling pathways and regulatory networks. We identified many differentially expressed genes that were either exclusively observed in a specific astrocytoma grade or commonly affected in specific subsets of astrocytoma grades in comparison to normal brain. Further, the number of differentially expressed genes generally increased with the astrocytoma grade with one major exception. The cytokine receptor pathway showed nearly the same number of differentially expressed genes in PA I and GBM IV and was further characterized by a significant overlap of commonly altered genes and an exclusive enrichment of overexpressed cancer genes in GBM IV. Additional analyses revealed a strong exclusive overexpression of CX3CL1 (fractalkine) and its receptor CX3CR1 in PA I possibly contributing to the absence of invasive growth.
    [Show full text]
  • Genome-Wide Analysis of HPV Integration in Human Cancers Reveals Recurrent, Focal Genomic Instability
    Downloaded from genome.cshlp.org on October 2, 2021 - Published by Cold Spring Harbor Laboratory Press Genome-wide analysis of HPV integration in human cancers reveals recurrent, focal genomic instability Keiko Akagi*a,b,c, Jingfeng Li*a,b,c, Tatevik R. Broutianb,d, Hesed Padilla-Nashf, Weihong Xiaob,d, Bo Jiangb,d, James W. Roccog,h, Theodoros N. Teknosi, Bhavna Kumari, Danny Wangsaf, Dandan Hea,b,c, Thomas Riedf, David E. Symer** ŧ a,b,c,d,e, Maura L. Gillison** ŧ b,d aHuman Cancer Genetics Program and bViral Oncology Program, Departments of cMolecular Virology, Immunology and Medical Genetics, dInternal Medicine and eBioinformatics, The Ohio State University Comprehensive Cancer Center, Columbus OH; fCancer Genomics Section, Center for Cancer Research, National Cancer Institute, Bethesda, MD; gCenter for Cancer Research and Department of Surgery, Massachusetts General Hospital, Boston, MA; hDepartment of Otolaryngology, Massachusetts Eye and Ear Infirmary, Harvard Medical School, Boston, MA; iDepartment of Otolaryngology-Head and Neck Surgery, The Ohio State University Medical Center *These authors contributed equally to this work **These authors contributed equally to this work ŧ Corresponding authors David E. Symer, M.D., Ph.D. [email protected] Maura L. Gillison M.D., Ph.D. [email protected] Downloaded from genome.cshlp.org on October 2, 2021 - Published by Cold Spring Harbor Laboratory Press Akagi and Li et al SUMMARY Genomic instability is a hallmark of human cancers, including the 5% caused by human papillomavirus (HPV). Here we report a striking association between HPV integration and adjacent host genomic structural variation in human cancer cell lines and primary tumors.
    [Show full text]
  • Histone 4 Lysine 20 Methylation: a Case for Neurodevelopmental Disease
    biology Review Histone 4 Lysine 20 Methylation: A Case for Neurodevelopmental Disease Rochelle N. Wickramasekara and Holly A. F. Stessman * Department of Pharmacology, School of Medicine, Creighton University, Omaha, NE 68178, USA; [email protected] * Correspondence: [email protected]; Tel.: +1-402-280-2255 Received: 28 December 2018; Accepted: 26 February 2019; Published: 3 March 2019 Abstract: Neurogenesis is an elegantly coordinated developmental process that must maintain a careful balance of proliferation and differentiation programs to be compatible with life. Due to the fine-tuning required for these processes, epigenetic mechanisms (e.g., DNA methylation and histone modifications) are employed, in addition to changes in mRNA transcription, to regulate gene expression. The purpose of this review is to highlight what we currently know about histone 4 lysine 20 (H4K20) methylation and its role in the developing brain. Utilizing publicly-available RNA-Sequencing data and published literature, we highlight the versatility of H4K20 methyl modifications in mediating diverse cellular events from gene silencing/chromatin compaction to DNA double-stranded break repair. From large-scale human DNA sequencing studies, we further propose that the lysine methyltransferase gene, KMT5B (OMIM: 610881), may fit into a category of epigenetic modifier genes that are critical for typical neurodevelopment, such as EHMT1 and ARID1B, which are associated with Kleefstra syndrome (OMIM: 610253) and Coffin-Siris syndrome (OMIM: 135900), respectively. Based on our current knowledge of the H4K20 methyl modification, we discuss emerging themes and interesting questions on how this histone modification, and particularly KMT5B expression, might impact neurodevelopment along with current challenges and potential avenues for future research.
    [Show full text]
  • The Genetic Program of Pancreatic Beta-Cell Replication in Vivo
    Page 1 of 65 Diabetes The genetic program of pancreatic beta-cell replication in vivo Agnes Klochendler1, Inbal Caspi2, Noa Corem1, Maya Moran3, Oriel Friedlich1, Sharona Elgavish4, Yuval Nevo4, Aharon Helman1, Benjamin Glaser5, Amir Eden3, Shalev Itzkovitz2, Yuval Dor1,* 1Department of Developmental Biology and Cancer Research, The Institute for Medical Research Israel-Canada, The Hebrew University-Hadassah Medical School, Jerusalem 91120, Israel 2Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, Israel. 3Department of Cell and Developmental Biology, The Silberman Institute of Life Sciences, The Hebrew University of Jerusalem, Jerusalem 91904, Israel 4Info-CORE, Bioinformatics Unit of the I-CORE Computation Center, The Hebrew University and Hadassah, The Institute for Medical Research Israel- Canada, The Hebrew University-Hadassah Medical School, Jerusalem 91120, Israel 5Endocrinology and Metabolism Service, Department of Internal Medicine, Hadassah-Hebrew University Medical Center, Jerusalem 91120, Israel *Correspondence: [email protected] Running title: The genetic program of pancreatic β-cell replication 1 Diabetes Publish Ahead of Print, published online March 18, 2016 Diabetes Page 2 of 65 Abstract The molecular program underlying infrequent replication of pancreatic beta- cells remains largely inaccessible. Using transgenic mice expressing GFP in cycling cells we sorted live, replicating beta-cells and determined their transcriptome. Replicating beta-cells upregulate hundreds of proliferation- related genes, along with many novel putative cell cycle components. Strikingly, genes involved in beta-cell functions, namely glucose sensing and insulin secretion were repressed. Further studies using single molecule RNA in situ hybridization revealed that in fact, replicating beta-cells double the amount of RNA for most genes, but this upregulation excludes genes involved in beta-cell function.
    [Show full text]