Supplement for TWO-SIGMA-G S1 Additional Type-I Error and Power Results

Total Page:16

File Type:pdf, Size:1020Kb

Supplement for TWO-SIGMA-G S1 Additional Type-I Error and Power Results Supplement for TWO-SIGMA-G Eric Van Buren, Ming Hu, Liang Chen, John Wrobel, Kirk Wilhelmsen, Lishan Su, Yun Li, Di Wu S1 Additional Type-I Error and Power Results Independent Genes (No IGC) Genes Simulated with IGC (A) No Gene-Level Random Effects (B) No Gene-Level Random Effects Test Size = 30, Ref. Size = 30 Test Size = 30, Ref. Size = 30 0.03 0.03 0.02 0.02 0.01 0.01 Type-I Error Type-I Error Observed Set-Level Observed Set-Level 0.00 0.00 0.0000 0.0025 0.0050 0.0075 0.0100 0.0000 0.0025 0.0050 0.0075 0.0100 Significance Threshold Significance Threshold CAMERA MAST TWO-SIGMA-G CAMERA MAST TWO-SIGMA-G Genes Simulated with IGC Genes Simulated with IGC (C) Gene-Level Random Effects Present (D) Gene-Level Random Effects Test Size = 30, Ref. Size = 30 Incorrectly Absent 0.03 Test Size = 30, Ref. Size = 30 0.03 0.02 0.02 0.01 0.01 Type-I Error Type-I Error Observed Set-Level 0.00 Observed Set-Level 0.00 0.0000 0.0025 0.0050 0.0075 0.0100 0.0000 0.0025 0.0050 0.0075 0.0100 Significance Threshold Significance Threshold CAMERA MAST TWO-SIGMA-G CAMERA MAST TWO-SIGMA-G Supplementary Figure S1: Type-I error performance of CAMERA, MAST, and TWO-SIGMA-G as significance threshold varies from 0 to .01. Each panel varies the existence of IGC between genes in the test set and the presence of gene-level random effect terms in gene-level model (CAMERA never includes gene-level random effect terms). Each plot combines six different settings, and 10 replicates per setting, which vary both the magnitude of the average inter-gene correlation (where applicable) in the test set and the nature of the correlation structure via the introduction of other individual-level covariates. See the Methods section of the main text for more details regarding the simulation procedure. 1 Independent Genes (No IGC) Genes Simulated with IGC (A) No Gene-Level Random Effects (B) No Gene-Level Random Effects Test Size = 30, Ref. Size = 30 Test Size = 30, Ref. Size = 30 0.100 0.100 0.075 0.075 0.050 0.050 0.025 0.025 Type-I Error Type-I Error Observed Set-Level Observed Set-Level 0.000 0.000 0.00 0.01 0.02 0.03 0.04 0.05 0.00 0.01 0.02 0.03 0.04 0.05 Significance Threshold Significance Threshold CAMERA MAST TWO-SIGMA-G CAMERA MAST TWO-SIGMA-G Genes Simulated with IGC Genes Simulated with IGC (C) Gene-Level Random Effects Present (D) Gene-Level Random Effects Test Size = 30, Ref. Size = 30 Incorrectly Absent 0.100 Test Size = 30, Ref. Size = 30 0.100 0.075 0.075 0.050 0.050 0.025 Type-I Error 0.025 Type-I Error Observed Set-Level 0.000 Observed Set-Level 0.000 0.00 0.01 0.02 0.03 0.04 0.05 0.00 0.01 0.02 0.03 0.04 0.05 Significance Threshold Significance Threshold CAMERA MAST TWO-SIGMA-G CAMERA MAST TWO-SIGMA-G Supplementary Figure S2: Type-I error performance of CAMERA, MAST, and TWO-SIGMA-G as significance threshold varies from 0 to .05. Each panel varies the existence of IGC between genes in the test set and the presence of gene-level random effect terms in gene-level model (CAMERA never includes gene-level random effect terms). Each plot combines six different settings, and 10 replicates per setting, which vary both the magnitude of the average inter-gene correlation (where applicable) in the test set and the nature of the correlation structure via the introduction of other individual-level covariates. See the Methods section of the main text for more details regarding the simulation procedure. 2 Independent Genes (No IGC) Genes Simulated with IGC (A) No Gene-Level Random Effects (B) No Gene-Level Random Effects Test Size = 30, Ref. Size = 100 Test Size = 30, Ref. Size = 100 0.4 0.4 0.3 0.3 0.2 0.2 0.1 0.1 0.0 0.0 Set-Level Type-I Error Unadjusted Adjusted Set-Level Type-I Error Unadjusted Adjusted p-values p-values CAMERA MAST TWO-SIGMA-G CAMERA MAST TWO-SIGMA-G Genes Simulated with IGC Genes Simulated with IGC (C) Gene-Level Random Effects Present (D) Gene-Level Random Effects Test Size = 30, Ref. Size = 100 Incorrectly Absent Test Size = 30, Ref. Size = 100 0.4 0.4 0.3 0.3 0.2 0.2 0.1 0.1 0.0 0.0 Set-Level Type-I Error Unadjusted Adjusted Set-Level Type-I Error Unadjusted Adjusted p-values p-values CAMERA MAST TWO-SIGMA-G CAMERA MAST TWO-SIGMA-G Supplementary Figure S3: Type-I error performance of CAMERA, MAST, and TWO-SIGMA-G using a reference set size of 100 genes. Each panel varies the existence of IGC between genes in the test set and the presence of gene-level random effect terms in gene-level model (CAMERA never includes gene-level random effect terms). Within each panel, both unadjusted and adjusted set-level p-values are plotted (unadjusted p-values are unavailable for MAST). Each boxplot aggregates six different settings which vary both the magnitude of the average inter-gene correlation (where applicable) in the test set and the nature of the correlation structure via the introduction of other individual-level covariates. Such settings are intended to represent the diversity seen in real data sets to paint an accurate picture of testing properties over a wide range of gene sets. Each of the six settings is further composed of 10 replicates which vary only random seed to mimic the impact of a different starting pool of cells from which genes were simulated. See the Methods section of the main text for more details regarding the simulation procedure. 3 Genes Simulated with IGC Genes Simulated with IGC (A) No Gene-Level Random Effects (B) Gene-Level Random Effects Test Size = 30, Ref. Size = 30 Incorrectly Absent 0.3 Test Size = 30, Ref. Size = 30 0.3 0.2 0.2 0.1 0.1 Set-Level0.0 Type-I Error 0.0 T0, T20, T50, Set-Level Type-I Error T0, T20, T50, R0 R20 R50 R0 R20 R50 Configuration (% of Genes DE in Configuration (% of Genes DE in Test/Reference Sets) Test/Reference Sets) CAMERA MAST TWO-SIGMA-G CAMERA MAST TWO-SIGMA-G Genes Simulated with IGC Genes Simulated with IGC (C) No Gene-Level Random Effects (D) Gene-Level Random Effects Test Size = 30, Ref. Size = 100 Incorrectly Absent 0.3 Test Size = 30, Ref. Size = 100 0.3 0.2 0.2 0.1 0.1 Set-Level0.0 Type-I Error 0.0 T0, T20, T50, Set-Level Type-I Error T0, T20, T50, R0 R20 R50 R0 R20 R50 Configuration (% of Genes DE in Configuration (% of Genes DE in Test/Reference Sets) Test/Reference Sets) CAMERA MAST TWO-SIGMA-G CAMERA MAST TWO-SIGMA-G Supplementary Figure S4: Type-I error performance of TWO-SIGMA-G, CAM- ERA, and MAST for various set-level null hypotheses. Each panel varies the refer- ence set size and the presence of gene-level random effect terms in gene-level model. Scenarios along the x-axis of each panel vary the percentage of genes that are differentially expressed (with the same effect size) in the test and reference sets. Each boxplot aggregates six differ- ent settings which vary both the magnitude of the average inter-gene correlation in the test set and the nature of the correlation structure via the introduction of other individual-level covariates. Such settings are intended to represent the diversity seen in real data sets to paint an accurate picture of testing properties over a wide range of gene sets. Each of the six settings is further composed of 10 replicates which vary only random seed to mimic the impact of a different starting pool of cells from which genes were simulated. See the Methods section of the main text for more details regarding the simulation procedure. 4 Genes Simulated with IGC Genes Simulated with IGC (A) No Gene-Level Random Effects (B) No Gene-Level Random Effects Test Size = 30, Ref. Size = 100 Test Size = 30, Ref. Size = 30 0.100 0.100 0.075 0.075 0.050 0.050 0.025 0.025 Set-Level0.000 Type-I Error Set-Level0.000 Type-I Error T0, T20, T50, T0, T20, T50, R0 R20 R50 R0 R20 R50 Configuration (% of Genes DE in Configuration (% of Genes DE in Test/Reference Sets) Test/Reference Sets) iDEA PAGE TWO-SIGMA-G iDEA PAGE TWO-SIGMA-G Supplementary Figure S5: Type-I error performance of iDEA, PAGE, and TWO- SIGMA for various set-level null hypotheses for genes simulated with IGC. Refer- ence set sizes of 100 and 30 are shown. Each boxplot aggregates six different settings which vary both the magnitude of the average inter-gene correlation in the test set and the nature of the correlation structure via the introduction of other individual-level covariates. Such settings are intended to represent the diversity seen in real data sets to paint an accurate picture of testing properties over a wide range of gene sets. Each of the six settings is further composed of 10 replicates which vary only random seed to mimic the impact of a different starting pool of cells from which genes were simulated.
Recommended publications
  • Low Abundance of the Matrix Arm of Complex I in Mitochondria Predicts Longevity in Mice
    ARTICLE Received 24 Jan 2014 | Accepted 9 Apr 2014 | Published 12 May 2014 DOI: 10.1038/ncomms4837 OPEN Low abundance of the matrix arm of complex I in mitochondria predicts longevity in mice Satomi Miwa1, Howsun Jow2, Karen Baty3, Amy Johnson1, Rafal Czapiewski1, Gabriele Saretzki1, Achim Treumann3 & Thomas von Zglinicki1 Mitochondrial function is an important determinant of the ageing process; however, the mitochondrial properties that enable longevity are not well understood. Here we show that optimal assembly of mitochondrial complex I predicts longevity in mice. Using an unbiased high-coverage high-confidence approach, we demonstrate that electron transport chain proteins, especially the matrix arm subunits of complex I, are decreased in young long-living mice, which is associated with improved complex I assembly, higher complex I-linked state 3 oxygen consumption rates and decreased superoxide production, whereas the opposite is seen in old mice. Disruption of complex I assembly reduces oxidative metabolism with concomitant increase in mitochondrial superoxide production. This is rescued by knockdown of the mitochondrial chaperone, prohibitin. Disrupted complex I assembly causes premature senescence in primary cells. We propose that lower abundance of free catalytic complex I components supports complex I assembly, efficacy of substrate utilization and minimal ROS production, enabling enhanced longevity. 1 Institute for Ageing and Health, Newcastle University, Newcastle upon Tyne NE4 5PL, UK. 2 Centre for Integrated Systems Biology of Ageing and Nutrition, Newcastle University, Newcastle upon Tyne NE4 5PL, UK. 3 Newcastle University Protein and Proteome Analysis, Devonshire Building, Devonshire Terrace, Newcastle upon Tyne NE1 7RU, UK. Correspondence and requests for materials should be addressed to T.v.Z.
    [Show full text]
  • UQCRFS1N Assembles Mitochondrial Respiratory Complex-III Into an Asymmetric 21-Subunit Dimer
    Protein Cell 2018, 9(6):586–591 https://doi.org/10.1007/s13238-018-0515-x Protein & Cell LETTER TO THE EDITOR UQCRFS1N assembles mitochondrial respiratory complex-III into an asymmetric 21-subunit dimer Dear Editor, Although several structures have been solved with very high resolution, the full length N-terminal processed peptide (1– Mitochondrial respiratory chain consists of four multimeric 78 amino acids, UQCRFS1N) of the iron-sulfur Rieske pro- protein complexes, Complex I-IV (CI, NADH dehydroge- tein (UQCRFS1) subunit has not been assigned in all of nase; CII, succinate:ubiquinone oxidoreductase; CIII, cyto- these structures (Table 1). UQCRFS1N is the N-terminal chrome bc1 complex; and CIV, cytochrome c oxidase). Cell mitochondrial targeting sequence of UQCRFS1, and after its These four complexes transfer electrons from NADH or cleavage from the precursor, this small peptide remains & FADH to oxygen and pump protons from mitochondrial 2 bound to CIII with unknown functions. In this letter, we show matrix to intermembrane space, generating electrochemical that one UQCRFS1N links the two 10-subunit CIII protomers gradient across the inner membrane which is harnessed by together to form the intact CIII, which resultantly contains complex V to synthesize ATP, providing the majority of only 21 subunits rather than previously assumed 22 subunits energy acquired by living organisms. Respiratory chain Protein (Fig. 1A and 1B). complexes were reported to interact with each other to form Firstly, we rebuilt the high-resolution crystal structures of supercomplexes, even megacomplex (Guo et al., 2017). bovine CIII (PDB: 2A06) (Huang et al., 2005) and chicken However, despite decades of intensive research, many CIII (PDB:3TGU) (Hao et al., 2012).
    [Show full text]
  • Differential Expression of Multiple Disease-Related Protein Groups
    brain sciences Article Differential Expression of Multiple Disease-Related Protein Groups Induced by Valproic Acid in Human SH-SY5Y Neuroblastoma Cells 1,2, 1, 1 1 Tsung-Ming Hu y, Hsiang-Sheng Chung y, Lieh-Yung Ping , Shih-Hsin Hsu , Hsin-Yao Tsai 1, Shaw-Ji Chen 3,4 and Min-Chih Cheng 1,* 1 Department of Psychiatry, Yuli Branch, Taipei Veterans General Hospital, Hualien 98142, Taiwan; [email protected] (T.-M.H.); [email protected] (H.-S.C.); [email protected] (L.-Y.P.); fi[email protected] (S.-H.H.); [email protected] (H.-Y.T.) 2 Department of Future Studies and LOHAS Industry, Fo Guang University, Jiaosi, Yilan County 26247, Taiwan 3 Department of Psychiatry, Mackay Medical College, New Taipei City 25245, Taiwan; [email protected] 4 Department of Psychiatry, Taitung Mackay Memorial Hospital, Taitung County 95064, Taiwan * Correspondence: [email protected]; Tel.: +886-3888-3141 (ext. 475) These authors contributed equally to this work. y Received: 10 July 2020; Accepted: 8 August 2020; Published: 12 August 2020 Abstract: Valproic acid (VPA) is a multifunctional medication used for the treatment of epilepsy, mania associated with bipolar disorder, and migraine. The pharmacological effects of VPA involve a variety of neurotransmitter and cell signaling systems, but the molecular mechanisms underlying its clinical efficacy is to date largely unknown. In this study, we used the isobaric tags for relative and absolute quantitation shotgun proteomic analysis to screen differentially expressed proteins in VPA-treated SH-SY5Y cells. We identified changes in the expression levels of multiple proteins involved in Alzheimer’s disease, Parkinson’s disease, chromatin remodeling, controlling gene expression via the vitamin D receptor, ribosome biogenesis, ubiquitin-mediated proteolysis, and the mitochondrial oxidative phosphorylation and electron transport chain.
    [Show full text]
  • Human Mitochondrial Pathologies of the Respiratory Chain and ATP Synthase: Contributions from Studies of Saccharomyces Cerevisiae
    life Review Human Mitochondrial Pathologies of the Respiratory Chain and ATP Synthase: Contributions from Studies of Saccharomyces cerevisiae Leticia V. R. Franco 1,2,* , Luca Bremner 1 and Mario H. Barros 2 1 Department of Biological Sciences, Columbia University, New York, NY 10027, USA; [email protected] 2 Department of Microbiology,Institute of Biomedical Sciences, Universidade de Sao Paulo, Sao Paulo 05508-900, Brazil; [email protected] * Correspondence: [email protected] Received: 27 October 2020; Accepted: 19 November 2020; Published: 23 November 2020 Abstract: The ease with which the unicellular yeast Saccharomyces cerevisiae can be manipulated genetically and biochemically has established this organism as a good model for the study of human mitochondrial diseases. The combined use of biochemical and molecular genetic tools has been instrumental in elucidating the functions of numerous yeast nuclear gene products with human homologs that affect a large number of metabolic and biological processes, including those housed in mitochondria. These include structural and catalytic subunits of enzymes and protein factors that impinge on the biogenesis of the respiratory chain. This article will review what is currently known about the genetics and clinical phenotypes of mitochondrial diseases of the respiratory chain and ATP synthase, with special emphasis on the contribution of information gained from pet mutants with mutations in nuclear genes that impair mitochondrial respiration. Our intent is to provide the yeast mitochondrial specialist with basic knowledge of human mitochondrial pathologies and the human specialist with information on how genes that directly and indirectly affect respiration were identified and characterized in yeast. Keywords: mitochondrial diseases; respiratory chain; yeast; Saccharomyces cerevisiae; pet mutants 1.
    [Show full text]
  • Single Cell Transcriptomics in Schizophrenia Postmortem Brain: Moving Beyond Bulk Lysate
    SINGLE CELL TRANSCRIPTOMICS IN SCHIZOPHRENIA POSTMORTEM BRAIN: MOVING BEYOND BULK LYSATE Richard Crist June 8th, 2020 Schizophrenia Transcriptomics ■ Microarray and RNA sequencing ■ Differentially expressed genes across many cortical and sub-cortical regions – Dorsolateral prefrontal cortex (dlPFC) (Fillman et al, 2013) – Anterior Cingulate Cortex (Zhao et al, 2015; Hong et al, 2013) – Superior temporal gyrus (Wu et al, 2012) – Hippocampus (Hwang et al, 2013; Kohen et al, 2014) – Amygdala (Chang et al, 2017) ■ Enrichment of pathways and gene networks – Neural development – Axon guidance – Inflammation and immune-related proteins CommonMind Consortium ■ Largest transcriptomic analysis of schizophrenia – 258 cases/279 controls – RNAseq in dlPFC ■ 693 differentially expressed genes Fromer et al, 2016 Cell Diversity in Postmortem Brain ■ Brain, like all tissues, consists of many cell types – Major cell populations (e.g. astrocytes) – Distinct sub-populations (e.g. PVALB+ interneurons) ■ Problems in assessing differential expression in bulk lysate – Inability to identify which cells are affected – Missed expression changes in less common cell types Penney et al, 2019 Schizophrenia Single Cell Transcriptomics ■ Immunofluorescence and laser capture microdissection to collect individual populations of cells ■ Layer III/V pyramidal neurons (Arion et al, 2017) – 72 PFC samples – 36 cases/36 controls – 100 cells per layer for each sample – Expression assessed by microarray – 1,783 differentially expressed probe sets corresponding to 1,420 genes
    [Show full text]
  • Alzheimer's Disease Is Associated With
    Alzheimer’s disease is associated with reduced expression of energy metabolism genes in posterior cingulate neurons Winnie S. Liang*†, Eric M. Reiman*†‡§, Jon Valla†¶, Travis Dunckley*†, Thomas G. Beach†ʈ, Andrew Grover†ʈ, Tracey L. Niedzielko†¶, Lonnie E. Schneider†¶, Diego Mastroeni†ʈ, Richard Caselli†**, Walter Kukull††, John C. Morris‡‡, Christine M. Hulette§§, Donald Schmechel§§, Joseph Rogers†ʈ, and Dietrich A. Stephan*†¶¶ *Neurogenomics Division, Translational Genomics Research Institute, 445 North Fifth Street, Phoenix, AZ 85004; †Arizona Alzheimer’s Consortium, 901 East Willetta Street, Phoenix, AZ 85006; ‡Banner Alzheimer’s Institute, 901 East Willetta Street Phoenix, AZ 85006; §Department of Psychiatry and Evelyn F. McKnight Brain Institute, University of Arizona, 1501 North Campbell Avenue, Tucson, AZ 85724; ¶Barrow Neurological Institute, 350 West Thomas Road, Phoenix, AZ 85013; ʈSun Health Research Institute, 10515 West Santa Fe Drive, Sun City, AZ 85351; **Department of Neurology, Mayo Clinic, 13400 East Shea Boulevard, Scottsdale, AZ 85259; ††National Alzheimer’s Coordinating Center, 4311 11th Avenue NE, No. 300, Seattle, WA 98105; ‡‡Washington University Alzheimer’s Disease Research Center, Washington University School of Medicine, 4488 Forest Park Avenue, Suite 101, St. Louis, MO 63108; and §§Bryan Alzheimer’s Disease Research Center, Duke University Medical Center, 2200 West Main Street, Suite A200, Durham, NC 27705 Edited by Marcus E. Raichle, Washington University School of Medicine, St. Louis, MO, and approved January 15, 2008 (received for review September 28, 2007) Alzheimer’s disease (AD) is associated with regional reductions in of genetic risk for AD) and were progressive in late-middle-aged fluorodeoxyglucose positron emission tomography (FDG PET) mea- persons (19).
    [Show full text]
  • Genomic Annotation of the Meningioma Tumor Suppressor Locus on Chromosome 1P34
    Oncogene (2004) 23, 1014–1020 & 2004 Nature Publishing Group All rights reserved 0950-9232/04 $25.00 www.nature.com/onc Genomic annotation of the meningioma tumor suppressor locus on chromosome 1p34 Erik P Sulman1,2, Peter S White2,3 and Garrett M Brodeur*,2,3 1The Fels Institute for Molecular Biology and Cancer Research, Temple University School of Medicine, Philadelphia, PA 19140, USA; 2Division of Oncology, Children’s Hospital of Philadelphia, University of Pennsylvania School of Medicine, Philadelphia, PA 19104, USA; 3Department of Pediatrics, University of Pennsylvania School of Medicine, Philadelphia, PA 19104, USA Meningioma is a frequently occurring tumor of the 13–20% of adult brain tumors, with an overall incidence meninges surrounding the central nervous system. Loss of approximately 7–8/100 000 per year, including of the short arm of chromosome 1 (1p) is the second most asymptomatic cases discovered at autopsy (Nakasu frequent chromosomal abnormality observed in these et al., 1987; Bondy and Ligon, 1996). Although usually tumors. Previously,we identified a 3.7 megabase (Mb) benign, 5–10% of tumors are malignant and invasive region of consistent deletion on 1p33–p34 in a panel of 157 and 20–40% of patients suffer from recurrent disease tumors. Loss of this region was associated with advanced (Zu¨ lch, 1979; Mathiesen et al., 1996). Meningiomas are disease and predictive for tumor relapse. In this report,a also associated with Neurofibromatosis type 2, occur- high-resolution integrated map of the region was con- ring nearly as frequently in that disease as vestibular structed (CompView) to identify all markers in the schwannomas (Evans et al., 1992).
    [Show full text]
  • Produktinformation
    Produktinformation Diagnostik & molekulare Diagnostik Laborgeräte & Service Zellkultur & Verbrauchsmaterial Forschungsprodukte & Biochemikalien Weitere Information auf den folgenden Seiten! See the following pages for more information! Lieferung & Zahlungsart Lieferung: frei Haus Bestellung auf Rechnung SZABO-SCANDIC Lieferung: € 10,- HandelsgmbH & Co KG Erstbestellung Vorauskassa Quellenstraße 110, A-1100 Wien T. +43(0)1 489 3961-0 Zuschläge F. +43(0)1 489 3961-7 [email protected] • Mindermengenzuschlag www.szabo-scandic.com • Trockeneiszuschlag • Gefahrgutzuschlag linkedin.com/company/szaboscandic • Expressversand facebook.com/szaboscandic SANTA CRUZ BIOTECHNOLOGY, INC. UQCRH siRNA (m): sc-106675 BACKGROUND PRODUCT UQCRH (ubiquinol-cytochrome c reductase hinge protein), also known as UQCRH siRNA (m) is a pool of 3 target-specific 19-25 nt siRNAs designed QCR6, is a 91 amino acid protein that localizes to the inner mitochondrial to knock down gene expression. Each vial contains 3.3 nmol of lyophilized membrane and exists as a component of the ubiquinol-cytochrome c reduc- siRNA, sufficient for a 10 µM solution once resuspended using protocol tase complex (known as complex III or as the cytochrome b-c1 complex). below. Suitable for 50-100 transfections. Also see UQCRH shRNA Functioning in conjunction with a variety of other proteins, UQCRH plays a role Plasmid (m): sc-106675-SH and UQCRH shRNA (m) Lentiviral Particles: in the mitochondrial respiratory chain and it thought to mediate cytochrome sc-106675-V as alternate gene silencing products. complex formation. The gene encoding UQCRH maps to human chromosome 1, For independent verification of UQCRH (m) gene silencing results, we also which spans 260 million base pairs, contains over 3,000 genes and comprises provide the individual siRNA duplex components.
    [Show full text]
  • UQCRH Downregulation Promotes Warburg Effect in Renal Cell
    www.nature.com/scientificreports OPEN UQCRH downregulation promotes Warburg efect in renal cell carcinoma cells Yanting Luo1,2,6, Louise Medina Bengtsson1,2,6, Xuechun Wang1,2,6, Tianhe Huang1,2, Guoqiang Liu1,2,3, Sean Murphy1,2, Caiqin Wang1,2, John Koren III4, Zachary Schafer1,2 & Xin Lu1,2,3,5* Ubiquinol-cytochrome c reductase hinge protein (UQCRH) is the hinge protein for the multi-subunit complex III of the mitochondrial electron transport chain and is involved in the electron transfer reaction between cytochrome c1 and c. Recent genome-wide transcriptomic and epigenomic profling of clear cell renal cell carcinoma (ccRCC) by The Cancer Genome Atlas (TCGA) identifed UQCRH as the top-ranked gene showing inverse correlation between DNA hypermethylation and mRNA downregulation. The function and underlying mechanism of UQCRH in the Warburg efect metabolism of ccRCC have not been characterized. Here, we verifed the clinical association of low UQCRH expression and shorter survival of ccRCC patients through in silico analysis and identifed KMRC2 as a highly relevant ccRCC cell line that displays hypermethylation-induced UQCRH extinction. Ectopic overexpression of UQCRH in KMRC2 restored mitochondrial membrane potential, increased oxygen consumption, and attenuated the Warburg efect at the cellular level. UQCRH overexpression in KMRC2 induced higher apoptosis and slowed down in vitro and in vivo tumor growth. UQCRH knockout by CRISPR/Cas9 had little impact on the metabolism and proliferation of 786O ccRCC cell line, suggesting the dispensable role of UQCRH in cells that have entered a Warburg-like state through other mechanisms. Together, our study suggests that loss of UQCRH expression by hypermethylation may promote kidney carcinogenesis through exacerbating the functional decline of mitochondria thus reinforcing the Warburg efect.
    [Show full text]
  • Mutations in CYC1, Encoding Cytochrome C1 Subunit Of
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector REPORT CYC1 c Mutations in , Encoding Cytochrome 1 Subunit of Respiratory Chain Complex III, Cause Insulin-Responsive Hyperglycemia Pauline Gaignard,1,21 Minal Menezes,2,21 Manuel Schiff,3,4,5 Aure´lien Bayot,3,4 Malgorzata Rak,3,4 He´le`ne Ogier de Baulny,5 Chen-Hsien Su,6 Mylene Gilleron,7,8 Anne Lombes,7,8 Heni Abida,6 Alexander Tzagoloff,6 Lisa Riley,9 Sandra T. Cooper,2,10 Kym Mina,11,12 Padma Sivadorai,13 Mark R. Davis,13 Richard J.N. Allcock,14,15 Nina Kresoje,14 Nigel G. Laing,16,17 David R. Thorburn,18,19 Abdelhamid Slama,1 John Christodoulou,2,9,20 and Pierre Rustin3,4,* Many individuals with abnormalities of mitochondrial respiratory chain complex III remain genetically undefined. Here, we report > > CYC1 c mutations (c.288G T [p.Trp96Cys] and c.643C T [p.Leu215Phe]) in , encoding the cytochrome 1 subunit of complex III, in two unrelated children presenting with recurrent episodes of ketoacidosis and insulin-responsive hyperglycemia. Cytochrome c1, the heme-containing component of complex III, mediates the transfer of electrons from the Rieske iron-sulfur protein to cytochrome c c . Cytochrome 1 is present at reduced levels in the skeletal muscle and skin fibroblasts of affected individuals. Moreover, studies on yeast mutants and affected individuals’ fibroblasts have shown that exogenous expression of wild-type CYC1 rescues complex III activity, demonstrating the deleterious effect of each mutation on cytochrome c1 stability and complex III activity.
    [Show full text]
  • Chemico-Biological Interactions 196 (2012) 89–95
    Chemico-Biological Interactions 196 (2012) 89–95 Contents lists available at ScienceDirect Chemico-Biological Interactions journal homepage: www.elsevier.com/locate/chembioint Exposure to sodium tungstate and Respiratory Syncytial Virus results in hematological/immunological disease in C57BL/6J mice ⇑ Cynthia D. Fastje a, , Kevin Harper a, Chad Terry a, Paul R. Sheppard b, Mark L. Witten a,1 a Steele Children’s Research Center, PO Box 245073, University of Arizona, Tucson, AZ 85724-5073, USA b Laboratory of Tree-Ring Research, PO Box 210058, University of Arizona, Tucson, AZ 85721-0058, USA article info abstract Article history: The etiology of childhood leukemia is not known. Strong evidence indicates that precursor B-cell Acute Available online 1 May 2011 Lymphoblastic Leukemia (Pre-B ALL) is a genetic disease originating in utero. Environmental exposures in two concurrent, childhood leukemia clusters have been profiled and compared with geographically Keywords: similar control communities. The unique exposures, shared in common by the leukemia clusters, have Tungsten been modeled in C57BL/6 mice utilizing prenatal exposures. This previous investigation has suggested Respiratory Syncytial Virus in utero exposure to sodium tungstate (Na2WO4) may result in hematological/immunological disease Childhood leukemia through genes associated with viral defense. The working hypothesis is (1) in addition to spontaneously and/or chemically generated genetic lesions forming pre-leukemic clones, in utero exposure to Na2WO4 increases genetic susceptibility to viral influence(s); (2) postnatal exposure to a virus possessing the 1FXXKXFXXA/V9 peptide motif will cause an unnatural immune response encouraging proliferation in the B-cell precursor compartment. This study reports the results of exposing C57BL/6J mice to Na2WO4 3 in utero via water (15 ppm, ad libetum) and inhalation (mean concentration PM5 3.33 mg/m ) and to Respiratory Syncytial Virus (RSV) within 2 weeks of weaning.
    [Show full text]
  • Mitochondrial Complex III Deficiency with Ketoacidosis and Hyperglycemia Mimicking Neonatal Diabetes
    JIMD Reports DOI 10.1007/8904_2016_557 CASE REPORT Mitochondrial Complex III Deficiency with Ketoacidosis and Hyperglycemia Mimicking Neonatal Diabetes Natascia Anastasio • Maja Tarailo-Graovac • Reem Al-Khalifah • Laurent Legault • Britt Drogemoller • Colin J.D. Ross • Wyeth W. Wasserman • Clara van Karnebeek • Daniela Buhas Received: 20 January 2016 /Revised: 10 March 2016 /Accepted: 15 March 2016 /Published online: 14 April 2016 # SSIEM and Springer-Verlag Berlin Heidelberg 2016 Abstract Hyperglycemia is a rare presenting symptom of tation, whole-exome sequencing (WES) was performed on mitochondrial disorders. We report a case of a young girl the parent–offspring trio. The patient was homozygous for who presented shortly after birth with ketoacidosis, hyper- the c.643C>T (p.Leu215Phe) variant in CYC1, a nuclear lactatemia, hyperammonemia, and insulin-responsive gene which encodes cytochrome c1, a subunit of respiratory hyperglycemia. Initial metabolic work-up suggested mito- chain complex III. Variants in this gene have only been chondrial dysfunction. Given our patient’s unusual presen- previously reported in two patients with similar presenta- tion, one of whom carries the same variant as our patient who is also of Sri Lankan origin. Communicated by: Garry Brown Primary complex III deficiencies are rare and its Competing interests: None declared phenotypes can vary significantly, even among patients N. Anastasio (*) : D. Buhas Department of Medical Genetics, McGill University, 1001 Boulevard with the same genotype. Décarie, Montréal, QC, Canada H4A 3J1 e-mail: [email protected] M. Tarailo-Graovac Centre for Molecular Medicine and Therapeutics, University of British Introduction Columbia, 950 West 28th Avenue, Vancouver, BC, Canada V5Z 4H4 : R.
    [Show full text]