Analysis Combining Correlated Glaucoma Traits Identifies Five New Risk Loci for Open-Angle Glaucoma

Total Page:16

File Type:pdf, Size:1020Kb

Analysis Combining Correlated Glaucoma Traits Identifies Five New Risk Loci for Open-Angle Glaucoma Analysis combining correlated glaucoma traits identifies five new risk loci for open-angle glaucoma The Harvard community has made this article openly available. Please share how this access benefits you. Your story matters Citation Gharahkhani, P., K. P. Burdon, J. N. Cooke Bailey, A. W. Hewitt, M. H. Law, L. R. Pasquale, J. H. Kang, et al. 2018. “Analysis combining correlated glaucoma traits identifies five new risk loci for open- angle glaucoma.” Scientific Reports 8 (1): 3124. doi:10.1038/ s41598-018-20435-9. http://dx.doi.org/10.1038/s41598-018-20435-9. Published Version doi:10.1038/s41598-018-20435-9 Citable link http://nrs.harvard.edu/urn-3:HUL.InstRepos:35014999 Terms of Use This article was downloaded from Harvard University’s DASH repository, and is made available under the terms and conditions applicable to Other Posted Material, as set forth at http:// nrs.harvard.edu/urn-3:HUL.InstRepos:dash.current.terms-of- use#LAA www.nature.com/scientificreports OPEN Analysis combining correlated glaucoma traits identifes fve new risk loci for open-angle glaucoma Received: 23 August 2017 Puya Gharahkhani 1, Kathryn P. Burdon 2, Jessica N. Cooke Bailey 3, Alex W. Hewitt 2, Accepted: 18 January 2018 Matthew H. Law 1, Louis R. Pasquale 4,5, Jae H. Kang 5, Jonathan L. Haines3, Published: xx xx xxxx Emmanuelle Souzeau 6, Tiger Zhou6, Owen M. Siggs 6, John Landers6, Mona Awadalla6, Shiwani Sharma6, Richard A. Mills6, Bronwyn Ridge6, David Lynn7, Robert Casson8, Stuart L. Graham9, Ivan Goldberg10, Andrew White10,11, Paul R. Healey10,11, John Grigg10, Mitchell Lawlor10, Paul Mitchell11, Jonathan Ruddle12, Michael Coote12, Mark Walland12, Stephen Best13, Andrea Vincent13, Jesse Gale14, Graham RadfordSmith1,15, David C. Whiteman1, Grant W. Montgomery1,16, Nicholas G. Martin1, David A Mackey2,17, Janey L. Wiggs4, Stuart MacGregor1, Jamie E. Craig6 & The NEIGHBORHOOD consortium* Open-angle glaucoma (OAG) is a major cause of blindness worldwide. To identify new risk loci for OAG, we performed a genome-wide association study in 3,071 OAG cases and 6,750 unscreened controls, and meta-analysed the results with GWAS data for intraocular pressure (IOP) and optic disc parameters (the overall meta-analysis sample size varying between 32,000 to 48,000 participants), which are glaucoma-related traits. We identifed and independently validated four novel genome-wide signifcant associations within or near MYOF and CYP26A1, LINC02052 and CRYGS, LMX1B, and LMO7 using single variant tests, one additional locus (C9) using gene-based tests, and two genetic pathways - “response to fuid shear stress” and “abnormal retina morphology” - in pathway-based tests. Interestingly, some of the new risk loci contribute to risk of other genetically-correlated eye diseases including myopia and age-related macular degeneration. To our knowledge, this study is the frst integrative study to combine genetic data from OAG and its correlated traits to identify new risk variants and genetic pathways, highlighting the future potential of combining genetic data from genetically-correlated eye traits for the purpose of gene discovery and mapping. OAG is characterized by optic nerve damage and progressive loss of peripheral vision, with many patients remain- ing undiagnosed until severe irreversible vision loss has occurred1,2. OAG has a signifcant genetic component 1QIMR Berghofer Medical Research Institute, Brisbane, Queensland, Australia. 2University of Tasmania, Hobart, Tasmania, Australia. 3Population and Quantitative Health Sciences, Institute for Computational Biology, Case Western Reserve University School of Medicine, Cleveland, OH, USA. 4Department of Ophthalmology, Harvard Medical School, Massachusetts Eye and Ear Infirmary, Boston, MA, USA. 5Channing Division of Network Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, USA. 6Department of Ophthalmology, Flinders University, Adelaide, South Australia, Australia. 7South Australian Health & Medical Research Institute, School of Medicine, Flinders University, Adelaide, South Australia, Australia. 8South Australian Institute of Ophthalmology, University of Adelaide, Adelaide, South Australia, Australia. 9Ophthalmology and Vision Science, Macquarie University, Sydney, New South Wales, Australia. 10Department of Ophthalmology, University of Sydney, Sydney, Australia. 11Centre for Vision Research, The Westmead Institute for Medical Research, University of Sydney, Westmead, NSW, Australia. 12Centre for Eye Research Australia (CERA), University of Melbourne, Royal Victorian Eye and Ear Hospital, Melbourne, Victoria, Australia. 13Department of Ophthalmology, University of Auckland, Auckland, New Zealand. 14Department of Ophthalmology, University of Otago, Dunedin, Otago, New Zealand. 15School of Medicine, University of Queensland, Herston Campus, Brisbane, QLD, Australia. 16Institute for Molecular Bioscience, The University of Queensland, Brisbane, Queensland, Australia. 17Centre for Ophthalmology and Visual Science, Lions Eye Institute, University of Western Australia, Perth, Australia. *A comprehensive list of consortium members appears at the end of the paper. Correspondence and requests for materials should be addressed to P.G. (email: [email protected]) or J.E.C. (email: [email protected]) SCIENTIFIC REPORTS | (2018) 8:3124 | DOI:10.1038/s41598-018-20435-9 1 www.nature.com/scientificreports/ Meta-analysis Chr SNP Risk allele P-value Analysis heterogeneity P Nearest Genes 10 rs72815193 G 6.10 × 10−10 OAG + VCDR 0.31 MYOF and XRCC6P1 3 rs56962872 G 2.81 × 10−8 OAG + VCDR 0.51 LINC02052 and CRYGS 9 rs6478746 G 4.54 × 10−8 OAG + CA 0.44 LOC105376277 and LMX1B 1 rs148639588 T 3.53 × 10−8 OAG + CA 0.71 COL11A1 Table 1. Association results for the best SNPs within the genome-wide signifcant regions in meta-analyses of ANZRAG OAG and the endophenotypes. Efect sizes of these SNPs on OAG are presented in Table 2. OAG, open-angle glaucoma; CA, cup area; VCDR, vertical cup to disk ratio. ANZRAG NEIGHBORHOOD combined Efect Other Chr SNP allele allele OR SE P-value OR SE P-value OR SE P-value 10 rs4918865^ C G 1.149 0.03 1.89 × 10−5 1.086 0.04 0.01958 1.119 0.02 2.31 × 10−6 3 rs56962872 A G 0.862 0.04 2.91 × 10−5 0.892 0.04 0.002324 0.876 0.03 3.03 × 10−7 9 rs6478746 A G 0.853 0.04 1.09 × 10−5 0.909 0.04 0.01231 0.879 0.03 9.10 × 10−7 13 rs9530458 T C 1.158 0.03 4.50 × 10−6 1.138 0.03 0.00018 1.148 0.02 3.45 × 10−9 Table 2. GWAS statistics for the new OAG loci in ANZRAG (the discovery OAG set), NEIGHBORHOOD (the replication OAG set), and combined (fxed-efect meta-analysis). ^rs4918865 in high LD with rs72815193, LD r2 = 0.93; OR, odds ratio; SE, standard error of regression coefcent. with a relative risk of over 9 in frst-degree relatives of afected individuals compared to relatives of unafected people3. Our previous genome-wide association studies (GWAS) have reproducibly identifed several risk loci for OAG including TMCO1, CDKN2B-AS1, SIX6, CAV1, CAV2, ABCA1, AFAP1, GMDS, ARHGEF12, TXNRD2, ATXN2, and FOXC14–9. However, the majority of the genetic variance contributing to OAG remains unexplained, emphasizing that further studies to identify additional risk loci for OAG are required in order to make genetic risk prediction more clinically useful. Optic disk parameters including cup area (CA; the central area), disc area (DA; the total area of optic disc including cup area and the surrounding area containing axons of the retinal ganglion cells), and vertical cup-disc ratio (VCDR; the ratio of the vertical diameter of cup area to the vertical diameter of the optic disc) are key measurements used to assess OAG diagnosis and progression10. Elevated intraocular pressure (IOP) is the major known risk factor for OAG2. We refer to CA, DA, VCDR, and IOP as OAG endophenotypes or quantitative traits. Tere are high genetic correlations between these quantitative traits and OAG, with several of the already known risk loci for these traits overlapping with each other and with OAG loci, demonstrating their utility as endo- phenotypes11–14. Tese fndings suggest that combining genetic data from OAG and its endophenotypes has the potential to increase the probability of identifying genetic variants that are common between traits, thus enabling the extraction of greater genetic power from valuable disease cohorts. In this study, we sought to identify additional risk loci contributing to OAG susceptibility by (1) increasing the sample size for OAG, (2) combining GWAS data from OAG and its endophenotypes in order to increase our statistical power to identify new risk loci for OAG, and (3) applying gene and pathway based approaches. Results In total, 3,071 OAG cases from the Australian & New Zealand Registry of Advanced Glaucoma (ANZRAG) obtained in three phases of data collection, and 6,750 unscreened controls of European descent were used as the GWAS discovery dataset in this study (Supplementary Table 1). Five loci were associated with OAG at genome-wide signifcance level in the meta-analysis of GWAS results between the three phases of ANZRAG data (P < 5 × 10−8), including regions near or within CDKN2B-AS1, ABCA1, C14orf39 and SIX6, TMCO1, and ARHGEF12, all of which are now well established risk loci for OAG5–7,9. Manhattan and Q-Q plots are shown in Supplementary Figure 1. Genomic infation factor lambda was 1.006 for this analysis. Next, to increase the power of this study to identify new risk loci for OAG, we performed GWAS meta-analyses of ANZRAG OAG and each of the endophenotypes (CA, DA, VCDR, and IOP) that we obtained from our previous study11 (Supplementary Table 1). Before performing the meta-analyses, we confrmed validity of the endophenotypes for OAG by showing that there were signifcant genetic correlations between OAG and the endo- phenotypes (ranging between 20% and 47%, Supplementary Table 2; p ≤ 0.018) using the cross-trait bivariate LD score regression approach15.
Recommended publications
  • Entrez Symbols Name Termid Termdesc 117553 Uba3,Ube1c
    Entrez Symbols Name TermID TermDesc 117553 Uba3,Ube1c ubiquitin-like modifier activating enzyme 3 GO:0016881 acid-amino acid ligase activity 299002 G2e3,RGD1310263 G2/M-phase specific E3 ubiquitin ligase GO:0016881 acid-amino acid ligase activity 303614 RGD1310067,Smurf2 SMAD specific E3 ubiquitin protein ligase 2 GO:0016881 acid-amino acid ligase activity 308669 Herc2 hect domain and RLD 2 GO:0016881 acid-amino acid ligase activity 309331 Uhrf2 ubiquitin-like with PHD and ring finger domains 2 GO:0016881 acid-amino acid ligase activity 316395 Hecw2 HECT, C2 and WW domain containing E3 ubiquitin protein ligase 2 GO:0016881 acid-amino acid ligase activity 361866 Hace1 HECT domain and ankyrin repeat containing, E3 ubiquitin protein ligase 1 GO:0016881 acid-amino acid ligase activity 117029 Ccr5,Ckr5,Cmkbr5 chemokine (C-C motif) receptor 5 GO:0003779 actin binding 117538 Waspip,Wip,Wipf1 WAS/WASL interacting protein family, member 1 GO:0003779 actin binding 117557 TM30nm,Tpm3,Tpm5 tropomyosin 3, gamma GO:0003779 actin binding 24779 MGC93554,Slc4a1 solute carrier family 4 (anion exchanger), member 1 GO:0003779 actin binding 24851 Alpha-tm,Tma2,Tmsa,Tpm1 tropomyosin 1, alpha GO:0003779 actin binding 25132 Myo5b,Myr6 myosin Vb GO:0003779 actin binding 25152 Map1a,Mtap1a microtubule-associated protein 1A GO:0003779 actin binding 25230 Add3 adducin 3 (gamma) GO:0003779 actin binding 25386 AQP-2,Aqp2,MGC156502,aquaporin-2aquaporin 2 (collecting duct) GO:0003779 actin binding 25484 MYR5,Myo1e,Myr3 myosin IE GO:0003779 actin binding 25576 14-3-3e1,MGC93547,Ywhah
    [Show full text]
  • Open Dogan Phdthesis Final.Pdf
    The Pennsylvania State University The Graduate School Eberly College of Science ELUCIDATING BIOLOGICAL FUNCTION OF GENOMIC DNA WITH ROBUST SIGNALS OF BIOCHEMICAL ACTIVITY: INTEGRATIVE GENOME-WIDE STUDIES OF ENHANCERS A Dissertation in Biochemistry, Microbiology and Molecular Biology by Nergiz Dogan © 2014 Nergiz Dogan Submitted in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy August 2014 ii The dissertation of Nergiz Dogan was reviewed and approved* by the following: Ross C. Hardison T. Ming Chu Professor of Biochemistry and Molecular Biology Dissertation Advisor Chair of Committee David S. Gilmour Professor of Molecular and Cell Biology Anton Nekrutenko Professor of Biochemistry and Molecular Biology Robert F. Paulson Professor of Veterinary and Biomedical Sciences Philip Reno Assistant Professor of Antropology Scott B. Selleck Professor and Head of the Department of Biochemistry and Molecular Biology *Signatures are on file in the Graduate School iii ABSTRACT Genome-wide measurements of epigenetic features such as histone modifications, occupancy by transcription factors and coactivators provide the opportunity to understand more globally how genes are regulated. While much effort is being put into integrating the marks from various combinations of features, the contribution of each feature to accuracy of enhancer prediction is not known. We began with predictions of 4,915 candidate erythroid enhancers based on genomic occupancy by TAL1, a key hematopoietic transcription factor that is strongly associated with gene induction in erythroid cells. Seventy of these DNA segments occupied by TAL1 (TAL1 OSs) were tested by transient transfections of cultured hematopoietic cells, and 56% of these were active as enhancers. Sixty-six TAL1 OSs were evaluated in transgenic mouse embryos, and 65% of these were active enhancers in various tissues.
    [Show full text]
  • Identification and Characterization of TPRKB Dependency in TP53 Deficient Cancers
    Identification and Characterization of TPRKB Dependency in TP53 Deficient Cancers. by Kelly Kennaley A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Molecular and Cellular Pathology) in the University of Michigan 2019 Doctoral Committee: Associate Professor Zaneta Nikolovska-Coleska, Co-Chair Adjunct Associate Professor Scott A. Tomlins, Co-Chair Associate Professor Eric R. Fearon Associate Professor Alexey I. Nesvizhskii Kelly R. Kennaley [email protected] ORCID iD: 0000-0003-2439-9020 © Kelly R. Kennaley 2019 Acknowledgements I have immeasurable gratitude for the unwavering support and guidance I received throughout my dissertation. First and foremost, I would like to thank my thesis advisor and mentor Dr. Scott Tomlins for entrusting me with a challenging, interesting, and impactful project. He taught me how to drive a project forward through set-backs, ask the important questions, and always consider the impact of my work. I’m truly appreciative for his commitment to ensuring that I would get the most from my graduate education. I am also grateful to the many members of the Tomlins lab that made it the supportive, collaborative, and educational environment that it was. I would like to give special thanks to those I’ve worked closely with on this project, particularly Dr. Moloy Goswami for his mentorship, Lei Lucy Wang, Dr. Sumin Han, and undergraduate students Bhavneet Singh, Travis Weiss, and Myles Barlow. I am also grateful for the support of my thesis committee, Dr. Eric Fearon, Dr. Alexey Nesvizhskii, and my co-mentor Dr. Zaneta Nikolovska-Coleska, who have offered guidance and critical evaluation since project inception.
    [Show full text]
  • Supplement 1 Microarray Studies
    EASE Categories Significantly Enriched in vs MG vs vs MGC4-2 Pt1-C vs C4-2 Pt1-C UP-Regulated Genes MG System Gene Category EASE Global MGRWV Pt1-N RWV Pt1-N Score FDR GO Molecular Extracellular matrix cellular construction 0.0008 0 110 genes up- Function Interpro EGF-like domain 0.0009 0 regulated GO Molecular Oxidoreductase activity\ acting on single dono 0.0015 0 Function GO Molecular Calcium ion binding 0.0018 0 Function Interpro Laminin-G domain 0.0025 0 GO Biological Process Cell Adhesion 0.0045 0 Interpro Collagen Triple helix repeat 0.0047 0 KEGG pathway Complement and coagulation cascades 0.0053 0 KEGG pathway Immune System – Homo sapiens 0.0053 0 Interpro Fibrillar collagen C-terminal domain 0.0062 0 Interpro Calcium-binding EGF-like domain 0.0077 0 GO Molecular Cell adhesion molecule activity 0.0105 0 Function EASE Categories Significantly Enriched in Down-Regulated Genes System Gene Category EASE Global Score FDR GO Biological Process Copper ion homeostasis 2.5E-09 0 Interpro Metallothionein 6.1E-08 0 Interpro Vertebrate metallothionein, Family 1 6.1E-08 0 GO Biological Process Transition metal ion homeostasis 8.5E-08 0 GO Biological Process Heavy metal sensitivity/resistance 1.9E-07 0 GO Biological Process Di-, tri-valent inorganic cation homeostasis 6.3E-07 0 GO Biological Process Metal ion homeostasis 6.3E-07 0 GO Biological Process Cation homeostasis 2.1E-06 0 GO Biological Process Cell ion homeostasis 2.1E-06 0 GO Biological Process Ion homeostasis 2.1E-06 0 GO Molecular Helicase activity 2.3E-06 0 Function GO Biological
    [Show full text]
  • 1 Supporting Information for a Microrna Network Regulates
    Supporting Information for A microRNA Network Regulates Expression and Biosynthesis of CFTR and CFTR-ΔF508 Shyam Ramachandrana,b, Philip H. Karpc, Peng Jiangc, Lynda S. Ostedgaardc, Amy E. Walza, John T. Fishere, Shaf Keshavjeeh, Kim A. Lennoxi, Ashley M. Jacobii, Scott D. Rosei, Mark A. Behlkei, Michael J. Welshb,c,d,g, Yi Xingb,c,f, Paul B. McCray Jr.a,b,c Author Affiliations: Department of Pediatricsa, Interdisciplinary Program in Geneticsb, Departments of Internal Medicinec, Molecular Physiology and Biophysicsd, Anatomy and Cell Biologye, Biomedical Engineeringf, Howard Hughes Medical Instituteg, Carver College of Medicine, University of Iowa, Iowa City, IA-52242 Division of Thoracic Surgeryh, Toronto General Hospital, University Health Network, University of Toronto, Toronto, Canada-M5G 2C4 Integrated DNA Technologiesi, Coralville, IA-52241 To whom correspondence should be addressed: Email: [email protected] (M.J.W.); yi- [email protected] (Y.X.); Email: [email protected] (P.B.M.) This PDF file includes: Materials and Methods References Fig. S1. miR-138 regulates SIN3A in a dose-dependent and site-specific manner. Fig. S2. miR-138 regulates endogenous SIN3A protein expression. Fig. S3. miR-138 regulates endogenous CFTR protein expression in Calu-3 cells. Fig. S4. miR-138 regulates endogenous CFTR protein expression in primary human airway epithelia. Fig. S5. miR-138 regulates CFTR expression in HeLa cells. Fig. S6. miR-138 regulates CFTR expression in HEK293T cells. Fig. S7. HeLa cells exhibit CFTR channel activity. Fig. S8. miR-138 improves CFTR processing. Fig. S9. miR-138 improves CFTR-ΔF508 processing. Fig. S10. SIN3A inhibition yields partial rescue of Cl- transport in CF epithelia.
    [Show full text]
  • Enzyme DHRS7
    Toward the identification of a function of the “orphan” enzyme DHRS7 Inauguraldissertation zur Erlangung der Würde eines Doktors der Philosophie vorgelegt der Philosophisch-Naturwissenschaftlichen Fakultät der Universität Basel von Selene Araya, aus Lugano, Tessin Basel, 2018 Originaldokument gespeichert auf dem Dokumentenserver der Universität Basel edoc.unibas.ch Genehmigt von der Philosophisch-Naturwissenschaftlichen Fakultät auf Antrag von Prof. Dr. Alex Odermatt (Fakultätsverantwortlicher) und Prof. Dr. Michael Arand (Korreferent) Basel, den 26.6.2018 ________________________ Dekan Prof. Dr. Martin Spiess I. List of Abbreviations 3α/βAdiol 3α/β-Androstanediol (5α-Androstane-3α/β,17β-diol) 3α/βHSD 3α/β-hydroxysteroid dehydrogenase 17β-HSD 17β-Hydroxysteroid Dehydrogenase 17αOHProg 17α-Hydroxyprogesterone 20α/βOHProg 20α/β-Hydroxyprogesterone 17α,20α/βdiOHProg 20α/βdihydroxyprogesterone ADT Androgen deprivation therapy ANOVA Analysis of variance AR Androgen Receptor AKR Aldo-Keto Reductase ATCC American Type Culture Collection CAM Cell Adhesion Molecule CYP Cytochrome P450 CBR1 Carbonyl reductase 1 CRPC Castration resistant prostate cancer Ct-value Cycle threshold-value DHRS7 (B/C) Dehydrogenase/Reductase Short Chain Dehydrogenase Family Member 7 (B/C) DHEA Dehydroepiandrosterone DHP Dehydroprogesterone DHT 5α-Dihydrotestosterone DMEM Dulbecco's Modified Eagle's Medium DMSO Dimethyl Sulfoxide DTT Dithiothreitol E1 Estrone E2 Estradiol ECM Extracellular Membrane EDTA Ethylenediaminetetraacetic acid EMT Epithelial-mesenchymal transition ER Endoplasmic Reticulum ERα/β Estrogen Receptor α/β FBS Fetal Bovine Serum 3 FDR False discovery rate FGF Fibroblast growth factor HEPES 4-(2-Hydroxyethyl)-1-Piperazineethanesulfonic Acid HMDB Human Metabolome Database HPLC High Performance Liquid Chromatography HSD Hydroxysteroid Dehydrogenase IC50 Half-Maximal Inhibitory Concentration LNCaP Lymph node carcinoma of the prostate mRNA Messenger Ribonucleic Acid n.d.
    [Show full text]
  • Spatiotemporal Regulation of Rap Guanine Nucleotide Exchange Factors
    Spatiotemporal Regulation of Rap Guanine Nucleotide Exchange Factors Sarah Valeria Consonni ISBN: 9789039361290 © S. V. Consonni, 2014 No part of this thesis may be reproduced in any form without prior written permission of the author. Printed by Gildeprint, The Netherlands Cover picture: “Enough already!” by Aaron de la Cruz, part of “In the family” project aimed at raising awereness about breast cancer. Patients diagnosed with breast cancer filled in their “cell” within the piece. Invitation picture by Jacob, 14 months old. The printing of this thesis was financially supported by the UMC Utrecht and Universiteit Utrecht Spatiotemporal Regulation of Rap Guanine Nucleotide Exchange Factors Regulatie van Rap Guanine Nucleotide Exchange Factors in Plaats en Tijd (met een samenvatting in het Nederlands) Regolazione dei Fattori di Scambio dei Nucleotidi Guaninici per Rap (con un riassunto in Italiano) Proefschrift ter verkrijging van de graad van doctor aan de Universiteit Utrecht op gezag van de rector magnificus, Prof. Dr. G.J. van der Zwaan, ingevolge het besluit van het college voor promoties in het openbaar te verdedigen op donderdag 15 mei 2014 des middags te 2.30 uur door Sarah Valeria Consonni geboren op 12 mei 1986 te Mariano Comense (CO), Italië Promotor: Prof. dr. J.L. Bos What would life be if we had no courage to attemp anything? Vincent van Gogh Table of contents Chapter 1 General Introduction 9 Chapter 2 cAMP regulates DEP domain-mediated binding of Epac1 to phosphatidic acid at the plasma membrane 19 Addendum Critical differences
    [Show full text]
  • Pan-Cancer Analysis of Homozygous Deletions in Primary Tumours Uncovers Rare Tumour Suppressors
    Corrected: Author correction; Corrected: Author correction ARTICLE DOI: 10.1038/s41467-017-01355-0 OPEN Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors Jiqiu Cheng1,2, Jonas Demeulemeester 3,4, David C. Wedge5,6, Hans Kristian M. Vollan2,3, Jason J. Pitt7,8, Hege G. Russnes2,9, Bina P. Pandey1, Gro Nilsen10, Silje Nord2, Graham R. Bignell5, Kevin P. White7,11,12,13, Anne-Lise Børresen-Dale2, Peter J. Campbell5, Vessela N. Kristensen2, Michael R. Stratton5, Ole Christian Lingjærde 10, Yves Moreau1 & Peter Van Loo 3,4 1234567890 Homozygous deletions are rare in cancers and often target tumour suppressor genes. Here, we build a compendium of 2218 primary tumours across 12 human cancer types and sys- tematically screen for homozygous deletions, aiming to identify rare tumour suppressors. Our analysis defines 96 genomic regions recurrently targeted by homozygous deletions. These recurrent homozygous deletions occur either over tumour suppressors or over fragile sites, regions of increased genomic instability. We construct a statistical model that separates fragile sites from regions showing signatures of positive selection for homozygous deletions and identify candidate tumour suppressors within those regions. We find 16 established tumour suppressors and propose 27 candidate tumour suppressors. Several of these genes (including MGMT, RAD17, and USP44) show prior evidence of a tumour suppressive function. Other candidate tumour suppressors, such as MAFTRR, KIAA1551, and IGF2BP2, are novel. Our study demonstrates how rare tumour suppressors can be identified through copy number meta-analysis. 1 Department of Electrical Engineering (ESAT) and iMinds Future Health Department, University of Leuven, Kasteelpark Arenberg 10, B-3001 Leuven, Belgium.
    [Show full text]
  • Human Lectins, Their Carbohydrate Affinities and Where to Find Them
    biomolecules Review Human Lectins, Their Carbohydrate Affinities and Where to Review HumanFind Them Lectins, Their Carbohydrate Affinities and Where to FindCláudia ThemD. Raposo 1,*, André B. Canelas 2 and M. Teresa Barros 1 1, 2 1 Cláudia D. Raposo * , Andr1 é LAQVB. Canelas‐Requimte,and Department M. Teresa of Chemistry, Barros NOVA School of Science and Technology, Universidade NOVA de Lisboa, 2829‐516 Caparica, Portugal; [email protected] 12 GlanbiaLAQV-Requimte,‐AgriChemWhey, Department Lisheen of Chemistry, Mine, Killoran, NOVA Moyne, School E41 of ScienceR622 Co. and Tipperary, Technology, Ireland; canelas‐ [email protected] NOVA de Lisboa, 2829-516 Caparica, Portugal; [email protected] 2* Correspondence:Glanbia-AgriChemWhey, [email protected]; Lisheen Mine, Tel.: Killoran, +351‐212948550 Moyne, E41 R622 Tipperary, Ireland; [email protected] * Correspondence: [email protected]; Tel.: +351-212948550 Abstract: Lectins are a class of proteins responsible for several biological roles such as cell‐cell in‐ Abstract:teractions,Lectins signaling are pathways, a class of and proteins several responsible innate immune for several responses biological against roles pathogens. such as Since cell-cell lec‐ interactions,tins are able signalingto bind to pathways, carbohydrates, and several they can innate be a immuneviable target responses for targeted against drug pathogens. delivery Since sys‐ lectinstems. In are fact, able several to bind lectins to carbohydrates, were approved they by canFood be and a viable Drug targetAdministration for targeted for drugthat purpose. delivery systems.Information In fact, about several specific lectins carbohydrate were approved recognition by Food by andlectin Drug receptors Administration was gathered for that herein, purpose. plus Informationthe specific organs about specific where those carbohydrate lectins can recognition be found by within lectin the receptors human was body.
    [Show full text]
  • Appendix 2. Significantly Differentially Regulated Genes in Term Compared with Second Trimester Amniotic Fluid Supernatant
    Appendix 2. Significantly Differentially Regulated Genes in Term Compared With Second Trimester Amniotic Fluid Supernatant Fold Change in term vs second trimester Amniotic Affymetrix Duplicate Fluid Probe ID probes Symbol Entrez Gene Name 1019.9 217059_at D MUC7 mucin 7, secreted 424.5 211735_x_at D SFTPC surfactant protein C 416.2 206835_at STATH statherin 363.4 214387_x_at D SFTPC surfactant protein C 295.5 205982_x_at D SFTPC surfactant protein C 288.7 1553454_at RPTN repetin solute carrier family 34 (sodium 251.3 204124_at SLC34A2 phosphate), member 2 238.9 206786_at HTN3 histatin 3 161.5 220191_at GKN1 gastrokine 1 152.7 223678_s_at D SFTPA2 surfactant protein A2 130.9 207430_s_at D MSMB microseminoprotein, beta- 99.0 214199_at SFTPD surfactant protein D major histocompatibility complex, class II, 96.5 210982_s_at D HLA-DRA DR alpha 96.5 221133_s_at D CLDN18 claudin 18 94.4 238222_at GKN2 gastrokine 2 93.7 1557961_s_at D LOC100127983 uncharacterized LOC100127983 93.1 229584_at LRRK2 leucine-rich repeat kinase 2 HOXD cluster antisense RNA 1 (non- 88.6 242042_s_at D HOXD-AS1 protein coding) 86.0 205569_at LAMP3 lysosomal-associated membrane protein 3 85.4 232698_at BPIFB2 BPI fold containing family B, member 2 84.4 205979_at SCGB2A1 secretoglobin, family 2A, member 1 84.3 230469_at RTKN2 rhotekin 2 82.2 204130_at HSD11B2 hydroxysteroid (11-beta) dehydrogenase 2 81.9 222242_s_at KLK5 kallikrein-related peptidase 5 77.0 237281_at AKAP14 A kinase (PRKA) anchor protein 14 76.7 1553602_at MUCL1 mucin-like 1 76.3 216359_at D MUC7 mucin 7,
    [Show full text]
  • Network-Based Genetic Profiling, and Therapeutic Target Identification Of
    bioRxiv preprint doi: https://doi.org/10.1101/480632; this version posted November 29, 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. Network-based genetic profiling, and therapeutic target identification of Thyroid Cancer 1st Md. Ali Hossain 2nd Tania Akter Asa 3rd Julian Quinn 4rdMd. Mijanur Rahman Dept. of CSE Dept. of EEE Bone biology divisions Dept. of CSE MIU, Jahangirnagar University Islamic University Garvan Institute of Medical Research JKKNIU Dhaka, Bangladesh Kushtia, Bangladesh Sydney, Australia Bangladesh ali:cse:bd@gmail:com tania:eee:iu@gmail:com j:quinn@garvan:org:au mijan cse@yahoo:com 5th Fazlul Huq 6th Mohammad Ali Moni Biomedical Sciences, Faculty of Medicine and Health Biomedical Sciences, Faculty of Medicine and Health The University of Sydney The University of Sydney Sydney, Australia Sydney, Australia fazlul:huq@sydney:edu:au mohammad:moni@sydney:edu:au Abstract—Molecular mechanisms that underlie the pathogene- I. INTRODUCTION sis and progression of malignant thyroid cancer (TC) are poorly understood. In this study, we employed network-based integrative Thyroid cancer (TC) is the most common endocrine ma- analyses of TC lesions to identify key molecules that are possible lignant disease [24], and although it has a relatively low hub genes and proteins in molecular pathways active in TC. We thus studied a microarray gene expression dataset (GSE82208, mortality rate compared to most other common metastatic n=52) that compared TC to benign thyroid tumours (follicular diseases its incidence is rising and is now the fifth most thyroid adenomas) in order to identify differentially expressed common cancer disease in women, notably affecting those genes (DEGs) in TC.
    [Show full text]
  • LMO7 Deficiency Reveals the Significance of the Cuticular Plate For
    ARTICLE https://doi.org/10.1038/s41467-019-09074-4 OPEN LMO7 deficiency reveals the significance of the cuticular plate for hearing function Ting-Ting Du1, James B. Dewey2, Elizabeth L. Wagner1, Runjia Cui3, Jinho Heo4, Jeong-Jin Park5, Shimon P. Francis1, Edward Perez-Reyes6, Stacey J. Guillot7, Nicholas E. Sherman5, Wenhao Xu8, John S Oghalai2, Bechara Kachar3 & Jung-Bum Shin1 Sensory hair cells, the mechanoreceptors of the auditory and vestibular systems, harbor 1234567890():,; two specialized elaborations of the apical surface, the hair bundle and the cuticular plate. In contrast to the extensively studied mechanosensory hair bundle, the cuticular plate is not as well understood. It is believed to provide a rigid foundation for stereocilia motion, but specifics about its function, especially the significance of its integrity for long-term maintenance of hair cell mechanotransduction, are not known. We discovered that a hair cell protein called LIM only protein 7 (LMO7) is specifically localized in the cuticular plate and the cell junction. Lmo7 KO mice suffer multiple cuticular plate deficiencies, including reduced filamentous actin density and abnormal stereociliar rootlets. In addition to the cuticular plate defects, older Lmo7 KO mice develop abnormalities in inner hair cell stereocilia. Together, these defects affect cochlear tuning and sensitivity and give rise to late-onset progressive hearing loss. 1 Department of Neuroscience, University of Virginia, Charlottesville, VA 22908, USA. 2 Caruso Department of Otolaryngology-Head and Neck Surgery, University of Southern California, Los Angeles, CA 90033, USA. 3 National Institute for Deafness and Communications Disorders, National Institute of Health, Bethesda, MD 20892, USA. 4 Center for Cell Signaling and Department of Microbiology, Immunology and Cancer Biology, University of Virginia, Charlottesville, VA 22908, USA.
    [Show full text]