Transcriptomic and Proteomic Profiling Provides Insight Into

Total Page:16

File Type:pdf, Size:1020Kb

Transcriptomic and Proteomic Profiling Provides Insight Into BASIC RESEARCH www.jasn.org Transcriptomic and Proteomic Profiling Provides Insight into Mesangial Cell Function in IgA Nephropathy † † ‡ Peidi Liu,* Emelie Lassén,* Viji Nair, Celine C. Berthier, Miyuki Suguro, Carina Sihlbom,§ † | † Matthias Kretzler, Christer Betsholtz, ¶ Börje Haraldsson,* Wenjun Ju, Kerstin Ebefors,* and Jenny Nyström* *Department of Physiology, Institute of Neuroscience and Physiology, §Proteomics Core Facility at University of Gothenburg, University of Gothenburg, Gothenburg, Sweden; †Division of Nephrology, Department of Internal Medicine and Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan; ‡Division of Molecular Medicine, Aichi Cancer Center Research Institute, Nagoya, Japan; |Department of Immunology, Genetics and Pathology, Uppsala University, Uppsala, Sweden; and ¶Integrated Cardio Metabolic Centre, Karolinska Institutet Novum, Huddinge, Sweden ABSTRACT IgA nephropathy (IgAN), the most common GN worldwide, is characterized by circulating galactose-deficient IgA (gd-IgA) that forms immune complexes. The immune complexes are deposited in the glomerular mesangium, leading to inflammation and loss of renal function, but the complete pathophysiology of the disease is not understood. Using an integrated global transcriptomic and proteomic profiling approach, we investigated the role of the mesangium in the onset and progression of IgAN. Global gene expression was investigated by microarray analysis of the glomerular compartment of renal biopsy specimens from patients with IgAN (n=19) and controls (n=22). Using curated glomerular cell type–specific genes from the published literature, we found differential expression of a much higher percentage of mesangial cell–positive standard genes than podocyte-positive standard genes in IgAN. Principal coordinate analysis of expression data revealed clear separation of patient and control samples on the basis of mesangial but not podocyte cell–positive standard genes. Additionally, patient clinical parameters (serum creatinine values and eGFRs) significantly correlated with Z scores derived from the expression profile of mesangial cell–positive standard genes. Among patients groupedaccordingtoOxfordMESTscore,patientswith segmental glomerulosclerosis had a significantly higher mesangial cell–positive standard gene Z score than patients without segmental glomerulosclerosis. By investi- gating mesangial cell proteomics and glomerular transcriptomics, we identified 22 common pathways induced in mesangial cells by gd-IgA, most of which mediate inflammation. The genes, proteins, and corresponding pathways identified provide novel insights into the pathophysiologic mechanisms leading to IgAN. J Am Soc Nephrol 28: ccc–ccc, 2017. doi: https://doi.org/10.1681/ASN.2016101103 The glomerular mesangium is a key element for a The diagnosis relies on histopathologic evaluation number of human kidney diseases, including dia- with findings of deposits of galactose-deficient IgA betic nephropathy, SLE, and IgA nephropathy (IgAN).1 However, knowledge of the role of Received October 18, 2016. Accepted May 5, 2017. mesangial cells in these diseases is not complete, Published online ahead of print. Publication date available at partly due to difficulties in developing appropriate www.jasn.org. 2 animal models and also because of the lack of Correspondence: Dr. Jenny Nyström, Institute of Neuroscience mesangial cell–specificmarkers.IgANisthemostcom- and Physiology, University of Gothenburg, Box 430, SE-405 30 monGNintheworld,and30%–40% of the patients Gothenburg, Sweden. Email: [email protected] develop ESRD within 20–30 years after diagnosis.3 Copyright © 2017 by the American Society of Nephrology J Am Soc Nephrol 28: ccc–ccc, 2017 ISSN : 1046-6673/2810-ccc 1 2 BASIC RESEARCH Journal of the American Society of Nephrology www.jasn.org Table 1. Clinical data for patients with IgAN Immuno Oxford Classification Urine Albumin-to-Creatinine Creatinine, BP Sex Age, yr eGFRa CKD Stage BP Steroids suppressive Progressionb Ratio mmol/L Medication Therapy MEST Woman 54 229 132 39.4 3 130/80 Yes No No 0 0 1 1 24.55 Man 46 323 155 45.6 3 150/90 Yes No No 1 1 1 1 0.60 Woman 45 155 132 40.5 3 140/80 Yes No No 1 0 1 1 23.05 Woman 41 189 165 33.0 3 130/80 Yes No No 0 1 1 1 25.37 Man 49 37 124 58.4 3 125/80 No No No 0 1 0 1 23.86 Man 48 65 139 51.3 3 125/85 No No No 0 0 1 1 218.99 Man 58 70 107 65.6 2 140/80 Yes No No 0 0 1 0 28.49 Man 35 14 103 80.7 2 120/80 Yes No No 0 0 1 0 2.40 Man 58 61 117 67.9 2 130/90 Yes No No 0 0 1 2 20.84 Man 34 5 124 64.9 2 126/85 Yes No No 0 0 1 1 1.18 Man 53 74 96 77.4 2 110/70 No Yesc No 0 1 0 0 227.74 Man 61 560d 84 86.0 2 150/95 Yes No No 0 0 1 0 23.94 Man 65 490d 93 74.0 2 155/95 Yes No No 0 0 1 1 21.88 Woman 17 15 71 108.2 1 130/80 No No No 1 1 1 0 22.24 Man 30 75 91 97.1 1 130/90 Yes No No 0 0 1 0 4.49 Man 20 1.2 76 124.8 1 100/60 No No No 0 0 0 0 26.94 Woman 23 0.6 68 109.3 1 110/65 No No No 0 0 0 0 3.92 Man 24 12 77 125.6 1 120/70 No No No 0 0 0 0 217.45 Woman 23 13 66 113.3 1 120/70 No No No 0 1 1 0 3.77 M, mesangial hypercellularity; E, endocapillary hypercellularity; S, segmental glomerulosclerosis; T, tubular atrophy/interstitial. aeGFR has been calculated using the CKD-EPI creatinine formula.bProgression is change in eGFR over time. c J Am Soc Nephrol Before the biopsy. dCalculated from total urinary protein. 28: ccc – ccc ,2017 www.jasn.org BASIC RESEARCH addition, by integration of the profiling information from both microarray and MS analyses, we have found several com- mon significantly differentially expressed pathways at both the transcriptomic and proteomic levels. Most of these are inflam- matory pathways that might be critical in the understanding of the underlying mechanisms of IgAN. RESULTS Statistics Reveal Important Transcriptomic Information about IgAN In Vivo After quality check, normalized microarray data from 19 pa- tients diagnosed with IgAN (one patient sample was run in duplicate in two different batches as an internal control) and 22 living donor controls were further analyzed (Table 1). A hierarchical clustering separated patient samples from the controls (Figure 1), indicating a distinct transcriptomic dif- ference between patients and controls. Using the Significant Figure 1. Hierarchical clustering of patient and control samples. Analysis of Microarray (SAM)11 method, we found 736 genes Clustering of the glomerular microarray data on the basis of whole- to be significantly differentially expressed between the groups genome expression values using Euclidean distance with the Ward (Supplemental Table 1). The significant genes were analyzed averaging method. Using whole-genome expression profiles, using the Ingenuity Pathway Analysis tool to investigate the samples from patients with IgAN were clearly distinguishable fi from living donor controls. corresponding signi cant pathways. 113 pathways were found to be significant with Fisher exact P,0.05, and most of the pathways were related to inflammation, cytokines, and growth factors. The (gd-IgA) in the glomerular mesangial region in addition to top 25 differentially expressed pathways are listed in Table 2. mesangial cell proliferation and mesangial matrix expansion.4 In an earlier study, it was shown that extracellular matrix Moreover, the amount of gd-IgA deposited is not correlated genes play an important role in IgAN in vitro.9 Toconfirm these with mesangial cell damage or kidney injury.5 However, data in patient material, we investigated Gene Ontology (GO) mesangialcellproliferationandinflammation induced by terms using the significantly differentially expressed genes. As gd-IgA deposits are important for disease development.6–8 seen in Table 3, the top three enriched GO terms for biologic In our previous study,9 we showed that mesangial cells de- processes are extracellular matrix disassembly (GO:0022617), rived from patients with IgAN have increased susceptibility to extracellular structure organization (GO:0043062), and extra- gd-IgA compared with mesangial cells from control patients. cellular matrix organization (GO:0030198). Additionally, Toincrease the understanding of the role of the mesangium in considering GO terms of cellular components as seen in Table IgAN, we investigated IgAN from a broader perspective using 3, extracellular matrix component (GO:0044420) and extra- both a transcriptomic and a proteomic approach. To achieve cellular matrix (GO:0031012) were also found to be highly this, we conducted a human genome microarray analysis of enriched along with basement membrane (GO:0005604) glomeruli from patients with IgAN and healthy controls. The and cell-cell adherens junction (GO:0005913). results were combined with proteomic data obtained from mass spectrometry (MS) analysis of mesangial cells treated Mesangial Cell–Positive Standard Genes Have a Strong with patient-derived gd-IgA. Effect in Patients with IgAN In 2013, Ju et al.10 developed an iterative algorithm “in silico On the basis of the algorithm presented in the study by Ju nanodissection” that predicts cell type–specifictranscriptsus- et al.,10 we investigated the expression profile of cell type– ing large-scale datasets. On the basis of literature curation, positive standard genes in the microarray data. In total, 35 they defined mesangial cell– or podocyte-positive standard positive standard genes were identified in mesangial cells, and genes as genes that are specifically expressed in mesangial cells 50 were identified in podocytes (for details of specific stan- or podocytes in the kidney. In this study, we have used these dard genes see Supplemental Table 2).
Recommended publications
  • Supplementary Data
    Figure 2S 4 7 A - C 080125 CSCs 080418 CSCs - + IFN-a 48 h + IFN-a 48 h + IFN-a 72 h 6 + IFN-a 72 h 3 5 MRFI 4 2 3 2 1 1 0 0 MHC I MHC II MICA MICB ULBP-1 ULBP-2 ULBP-3 ULBP-4 MHC I MHC II MICA MICB ULBP-1 ULBP-2 ULBP-3 ULBP-4 7 B 13 080125 FBS - D 080418 FBS - + IFN-a 48 h 12 + IFN-a 48 h + IFN-a 72 h + IFN-a 72 h 6 080125 FBS 11 10 5 9 8 4 7 6 3 MRFI 5 4 2 3 2 1 1 0 0 MHC I MHC II MICA MICB ULBP-1 ULBP-2 ULBP-3 ULBP-4 MHC I MHC II MICA MICB ULBP-1 ULBP-2 ULBP-3 ULBP-4 Molecule Molecule FIGURE 4S FIGURE 5S Panel A Panel B FIGURE 6S A B C D Supplemental Results Table 1S. Modulation by IFN-α of APM in GBM CSC and FBS tumor cell lines. Molecule * Cell line IFN-α‡ HLA β2-m# HLA LMP TAP1 TAP2 class II A A HC§ 2 7 10 080125 CSCs - 1∞ (1) 3 (65) 2 (91) 1 (2) 6 (47) 2 (61) 1 (3) 1 (2) 1 (3) + 2 (81) 11 (80) 13 (99) 1 (3) 8 (88) 4 (91) 1 (2) 1 (3) 2 (68) 080125 FBS - 2 (81) 4 (63) 4 (83) 1 (3) 6 (80) 3 (67) 2 (86) 1 (3) 2 (75) + 2 (99) 14 (90) 7 (97) 5 (75) 7 (100) 6 (98) 2 (90) 1 (4) 3 (87) 080418 CSCs - 2 (51) 1 (1) 1 (3) 2 (47) 2 (83) 2 (54) 1 (4) 1 (2) 1 (3) + 2 (81) 3 (76) 5 (75) 2 (50) 2 (83) 3 (71) 1 (3) 2 (87) 1 (2) 080418 FBS - 1 (3) 3 (70) 2 (88) 1 (4) 3 (87) 2 (76) 1 (3) 1 (3) 1 (2) + 2 (78) 7 (98) 5 (99) 2 (94) 5 (100) 3 (100) 1 (4) 2 (100) 1 (2) 070104 CSCs - 1 (2) 1 (3) 1 (3) 2 (78) 1 (3) 1 (2) 1 (3) 1 (3) 1 (2) + 2 (98) 8 (100) 10 (88) 4 (89) 3 (98) 3 (94) 1 (4) 2 (86) 2 (79) * expression of APM molecules was evaluated by intracellular staining and cytofluorimetric analysis; ‡ cells were treatead or not (+/-) for 72 h with 1000 IU/ml of IFN-α; # β-2 microglobulin; § β-2 microglobulin-free HLA-A heavy chain; ∞ values are indicated as ratio between the mean of fluorescence intensity of cells stained with the selected mAb and that of the negative control; bold values indicate significant MRFI (≥ 2).
    [Show full text]
  • Genome-Wide Association Study of Increasing Suicidal Ideation During Antidepressant Treatment in the GENDEP Project
    The Pharmacogenomics Journal (2012) 12, 68–77 & 2012 Macmillan Publishers Limited. All rights reserved 1470-269X/12 www.nature.com/tpj ORIGINAL ARTICLE Genome-wide association study of increasing suicidal ideation during antidepressant treatment in the GENDEP project NPerroud1,2,3,RUher1, Suicidal thoughts during antidepressant treatment have been the focus of 1 2,3,4 several candidate gene association studies. The aim of the present genome- MYM Ng , M Guipponi , wide association study was to identify additional genetic variants involved in 5 6 JHauser, N Henigsberg , increasing suicidal ideation during escitalopram and nortriptyline treatment. WMaier7,OMors8, A total of 706 adult participants of European ancestry, treated for major M Gennarelli9, M Rietschel10, depression with escitalopram or nortriptyline over 12 weeks in the Genome- DSouery11, MZ Dernovsek12, Based Therapeutic Drugs for Depression (GENDEP) study were genotyped 13 14 with Illumina Human 610-Quad Beadchips (Illumina, San Diego, CA, USA). AS Stamp ,MLathrop , A total of 244 subjects experienced an increase in suicidal ideation AFarmer1, G Breen1,KJAitchison1, during follow-up. The genetic marker most significantly associated with CM Lewis1,IWCraig1 increasing suicidality (8.28 Â 10À7) was a single-nucleotide polymorphism and P McGuffin1 (SNP; rs11143230) located 30 kb downstream of a gene encoding guanine deaminase (GDA) on chromosome 9q21.13. Two suggestive drug-specific 1MRC Social, Genetic and Developmental Psychiatry associations within KCNIP4 (Kv channel-interacting protein 4; chromosome Centre, Institute of Psychiatry, King’s College 4p15.31) and near ELP3 (elongation protein 3 homolog; chromosome London, London, UK; 2Department of Psychiatry, 8p21.1) were found in subjects treated with escitalopram.
    [Show full text]
  • Screening and Identification of Key Biomarkers in Clear Cell Renal Cell Carcinoma Based on Bioinformatics Analysis
    bioRxiv preprint doi: https://doi.org/10.1101/2020.12.21.423889; this version posted December 23, 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. Screening and identification of key biomarkers in clear cell renal cell carcinoma based on bioinformatics analysis Basavaraj Vastrad1, Chanabasayya Vastrad*2 , Iranna Kotturshetti 1. Department of Biochemistry, Basaveshwar College of Pharmacy, Gadag, Karnataka 582103, India. 2. Biostatistics and Bioinformatics, Chanabasava Nilaya, Bharthinagar, Dharwad 580001, Karanataka, India. 3. Department of Ayurveda, Rajiv Gandhi Education Society`s Ayurvedic Medical College, Ron, Karnataka 562209, India. * Chanabasayya Vastrad [email protected] Ph: +919480073398 Chanabasava Nilaya, Bharthinagar, Dharwad 580001 , Karanataka, India bioRxiv preprint doi: https://doi.org/10.1101/2020.12.21.423889; this version posted December 23, 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. Abstract Clear cell renal cell carcinoma (ccRCC) is one of the most common types of malignancy of the urinary system. The pathogenesis and effective diagnosis of ccRCC have become popular topics for research in the previous decade. In the current study, an integrated bioinformatics analysis was performed to identify core genes associated in ccRCC. An expression dataset (GSE105261) was downloaded from the Gene Expression Omnibus database, and included 26 ccRCC and 9 normal kideny samples. Assessment of the microarray dataset led to the recognition of differentially expressed genes (DEGs), which was subsequently used for pathway and gene ontology (GO) enrichment analysis.
    [Show full text]
  • A Review on Agmatinase Inhibitors
    www.ijcrt.org © 2020 IJCRT | Volume 8, Issue 12 December 2020 | ISSN: 2320-2882 A review on Agmatinase inhibitors 1Sunnica Biswas, 2Mr. R. T. Lohiya, 3Dr. Milind Umekar *1&2 Department Of Pharmaceutical Chemistry Smt. Kishoriti Bhoyar College of Pharmacy, Kamptee, Dst. Nagpur, 441002 Abstract : Agmatine is the product of arginine decarboxylation and can be hydrolyzed by agmatinase to putrescine, the precursor for biosynthesis of higher polyamines, spermidine, and spermine. Besides being an intermediate in polyamine metabolism, recent findings indicate that agmatine may play important regulatory roles in mammals. Agmatine, 4-aminobutyl guanidine, has recently been found in various mammalian organs and is thought to act as a neurotransmitter or neuromodulatory agent. The present study is to do a review on agmatine and its synthesized analogues till now for agmatinase inhibitory action. Agmatinase is a binuclear manganese metalloenzyme and belongs to the ureohydrolase superfamily that includes arginase, formiminoglutamase, and proclavaminate amidinohydrolase. Compared with a wealth of structural information available for arginases, no three dimensional structure of agmatinase has been reported. Agmatinase is an enzyme which blocks the mammalian agmatine which is ultimately responsible for the agmatine degradation in the body. Agmatinase is an enzyme which regulates the half life of agmatine in the brain. Hence a selective inhibitor of brain agmatinase is required. Several derivatives of agmatine are synthesized previously for agmatinase inhibitory activity but none of them showed selective inhibition. PZC (Piperazinecarboxamidine) is a derivative of agmatine or guanidine is expected to show selective inhibition of human agmatinase. A detailed review is carried out in order to understand the agmatinase inhibitor.
    [Show full text]
  • PARSANA-DISSERTATION-2020.Pdf
    DECIPHERING TRANSCRIPTIONAL PATTERNS OF GENE REGULATION: A COMPUTATIONAL APPROACH by Princy Parsana A dissertation submitted to The Johns Hopkins University in conformity with the requirements for the degree of Doctor of Philosophy Baltimore, Maryland July, 2020 © 2020 Princy Parsana All rights reserved Abstract With rapid advancements in sequencing technology, we now have the ability to sequence the entire human genome, and to quantify expression of tens of thousands of genes from hundreds of individuals. This provides an extraordinary opportunity to learn phenotype relevant genomic patterns that can improve our understanding of molecular and cellular processes underlying a trait. The high dimensional nature of genomic data presents a range of computational and statistical challenges. This dissertation presents a compilation of projects that were driven by the motivation to efficiently capture gene regulatory patterns in the human transcriptome, while addressing statistical and computational challenges that accompany this data. We attempt to address two major difficulties in this domain: a) artifacts and noise in transcriptomic data, andb) limited statistical power. First, we present our work on investigating the effect of artifactual variation in gene expression data and its impact on trans-eQTL discovery. Here we performed an in-depth analysis of diverse pre-recorded covariates and latent confounders to understand their contribution to heterogeneity in gene expression measurements. Next, we discovered 673 trans-eQTLs across 16 human tissues using v6 data from the Genotype Tissue Expression (GTEx) project. Finally, we characterized two trait-associated trans-eQTLs; one in Skeletal Muscle and another in Thyroid. Second, we present a principal component based residualization method to correct gene expression measurements prior to reconstruction of co-expression networks.
    [Show full text]
  • Identification of the Binding Partners for Hspb2 and Cryab Reveals
    Brigham Young University BYU ScholarsArchive Theses and Dissertations 2013-12-12 Identification of the Binding arP tners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactions and Non- Redundant Roles for Small Heat Shock Proteins Kelsey Murphey Langston Brigham Young University - Provo Follow this and additional works at: https://scholarsarchive.byu.edu/etd Part of the Microbiology Commons BYU ScholarsArchive Citation Langston, Kelsey Murphey, "Identification of the Binding Partners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactions and Non-Redundant Roles for Small Heat Shock Proteins" (2013). Theses and Dissertations. 3822. https://scholarsarchive.byu.edu/etd/3822 This Thesis is brought to you for free and open access by BYU ScholarsArchive. It has been accepted for inclusion in Theses and Dissertations by an authorized administrator of BYU ScholarsArchive. For more information, please contact [email protected], [email protected]. Identification of the Binding Partners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactions and Non-Redundant Roles for Small Heat Shock Proteins Kelsey Langston A thesis submitted to the faculty of Brigham Young University in partial fulfillment of the requirements for the degree of Master of Science Julianne H. Grose, Chair William R. McCleary Brian Poole Department of Microbiology and Molecular Biology Brigham Young University December 2013 Copyright © 2013 Kelsey Langston All Rights Reserved ABSTRACT Identification of the Binding Partners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactors and Non-Redundant Roles for Small Heat Shock Proteins Kelsey Langston Department of Microbiology and Molecular Biology, BYU Master of Science Small Heat Shock Proteins (sHSP) are molecular chaperones that play protective roles in cell survival and have been shown to possess chaperone activity.
    [Show full text]
  • Supplemental Table S1
    Entrez Gene Symbol Gene Name Affymetrix EST Glomchip SAGE Stanford Literature HPA confirmed Gene ID Profiling profiling Profiling Profiling array profiling confirmed 1 2 A2M alpha-2-macroglobulin 0 0 0 1 0 2 10347 ABCA7 ATP-binding cassette, sub-family A (ABC1), member 7 1 0 0 0 0 3 10350 ABCA9 ATP-binding cassette, sub-family A (ABC1), member 9 1 0 0 0 0 4 10057 ABCC5 ATP-binding cassette, sub-family C (CFTR/MRP), member 5 1 0 0 0 0 5 10060 ABCC9 ATP-binding cassette, sub-family C (CFTR/MRP), member 9 1 0 0 0 0 6 79575 ABHD8 abhydrolase domain containing 8 1 0 0 0 0 7 51225 ABI3 ABI gene family, member 3 1 0 1 0 0 8 29 ABR active BCR-related gene 1 0 0 0 0 9 25841 ABTB2 ankyrin repeat and BTB (POZ) domain containing 2 1 0 1 0 0 10 30 ACAA1 acetyl-Coenzyme A acyltransferase 1 (peroxisomal 3-oxoacyl-Coenzyme A thiol 0 1 0 0 0 11 43 ACHE acetylcholinesterase (Yt blood group) 1 0 0 0 0 12 58 ACTA1 actin, alpha 1, skeletal muscle 0 1 0 0 0 13 60 ACTB actin, beta 01000 1 14 71 ACTG1 actin, gamma 1 0 1 0 0 0 15 81 ACTN4 actinin, alpha 4 0 0 1 1 1 10700177 16 10096 ACTR3 ARP3 actin-related protein 3 homolog (yeast) 0 1 0 0 0 17 94 ACVRL1 activin A receptor type II-like 1 1 0 1 0 0 18 8038 ADAM12 ADAM metallopeptidase domain 12 (meltrin alpha) 1 0 0 0 0 19 8751 ADAM15 ADAM metallopeptidase domain 15 (metargidin) 1 0 0 0 0 20 8728 ADAM19 ADAM metallopeptidase domain 19 (meltrin beta) 1 0 0 0 0 21 81792 ADAMTS12 ADAM metallopeptidase with thrombospondin type 1 motif, 12 1 0 0 0 0 22 9507 ADAMTS4 ADAM metallopeptidase with thrombospondin type 1
    [Show full text]
  • A Computational Approach for Defining a Signature of Β-Cell Golgi Stress in Diabetes Mellitus
    Page 1 of 781 Diabetes A Computational Approach for Defining a Signature of β-Cell Golgi Stress in Diabetes Mellitus Robert N. Bone1,6,7, Olufunmilola Oyebamiji2, Sayali Talware2, Sharmila Selvaraj2, Preethi Krishnan3,6, Farooq Syed1,6,7, Huanmei Wu2, Carmella Evans-Molina 1,3,4,5,6,7,8* Departments of 1Pediatrics, 3Medicine, 4Anatomy, Cell Biology & Physiology, 5Biochemistry & Molecular Biology, the 6Center for Diabetes & Metabolic Diseases, and the 7Herman B. Wells Center for Pediatric Research, Indiana University School of Medicine, Indianapolis, IN 46202; 2Department of BioHealth Informatics, Indiana University-Purdue University Indianapolis, Indianapolis, IN, 46202; 8Roudebush VA Medical Center, Indianapolis, IN 46202. *Corresponding Author(s): Carmella Evans-Molina, MD, PhD ([email protected]) Indiana University School of Medicine, 635 Barnhill Drive, MS 2031A, Indianapolis, IN 46202, Telephone: (317) 274-4145, Fax (317) 274-4107 Running Title: Golgi Stress Response in Diabetes Word Count: 4358 Number of Figures: 6 Keywords: Golgi apparatus stress, Islets, β cell, Type 1 diabetes, Type 2 diabetes 1 Diabetes Publish Ahead of Print, published online August 20, 2020 Diabetes Page 2 of 781 ABSTRACT The Golgi apparatus (GA) is an important site of insulin processing and granule maturation, but whether GA organelle dysfunction and GA stress are present in the diabetic β-cell has not been tested. We utilized an informatics-based approach to develop a transcriptional signature of β-cell GA stress using existing RNA sequencing and microarray datasets generated using human islets from donors with diabetes and islets where type 1(T1D) and type 2 diabetes (T2D) had been modeled ex vivo. To narrow our results to GA-specific genes, we applied a filter set of 1,030 genes accepted as GA associated.
    [Show full text]
  • Supplementary Table 3 Complete List of RNA-Sequencing Analysis of Gene Expression Changed by ≥ Tenfold Between Xenograft and Cells Cultured in 10%O2
    Supplementary Table 3 Complete list of RNA-Sequencing analysis of gene expression changed by ≥ tenfold between xenograft and cells cultured in 10%O2 Expr Log2 Ratio Symbol Entrez Gene Name (culture/xenograft) -7.182 PGM5 phosphoglucomutase 5 -6.883 GPBAR1 G protein-coupled bile acid receptor 1 -6.683 CPVL carboxypeptidase, vitellogenic like -6.398 MTMR9LP myotubularin related protein 9-like, pseudogene -6.131 SCN7A sodium voltage-gated channel alpha subunit 7 -6.115 POPDC2 popeye domain containing 2 -6.014 LGI1 leucine rich glioma inactivated 1 -5.86 SCN1A sodium voltage-gated channel alpha subunit 1 -5.713 C6 complement C6 -5.365 ANGPTL1 angiopoietin like 1 -5.327 TNN tenascin N -5.228 DHRS2 dehydrogenase/reductase 2 leucine rich repeat and fibronectin type III domain -5.115 LRFN2 containing 2 -5.076 FOXO6 forkhead box O6 -5.035 ETNPPL ethanolamine-phosphate phospho-lyase -4.993 MYO15A myosin XVA -4.972 IGF1 insulin like growth factor 1 -4.956 DLG2 discs large MAGUK scaffold protein 2 -4.86 SCML4 sex comb on midleg like 4 (Drosophila) Src homology 2 domain containing transforming -4.816 SHD protein D -4.764 PLP1 proteolipid protein 1 -4.764 TSPAN32 tetraspanin 32 -4.713 N4BP3 NEDD4 binding protein 3 -4.705 MYOC myocilin -4.646 CLEC3B C-type lectin domain family 3 member B -4.646 C7 complement C7 -4.62 TGM2 transglutaminase 2 -4.562 COL9A1 collagen type IX alpha 1 chain -4.55 SOSTDC1 sclerostin domain containing 1 -4.55 OGN osteoglycin -4.505 DAPL1 death associated protein like 1 -4.491 C10orf105 chromosome 10 open reading frame 105 -4.491
    [Show full text]
  • 1 Metabolic Dysfunction Is Restricted to the Sciatic Nerve in Experimental
    Page 1 of 255 Diabetes Metabolic dysfunction is restricted to the sciatic nerve in experimental diabetic neuropathy Oliver J. Freeman1,2, Richard D. Unwin2,3, Andrew W. Dowsey2,3, Paul Begley2,3, Sumia Ali1, Katherine A. Hollywood2,3, Nitin Rustogi2,3, Rasmus S. Petersen1, Warwick B. Dunn2,3†, Garth J.S. Cooper2,3,4,5* & Natalie J. Gardiner1* 1 Faculty of Life Sciences, University of Manchester, UK 2 Centre for Advanced Discovery and Experimental Therapeutics (CADET), Central Manchester University Hospitals NHS Foundation Trust, Manchester Academic Health Sciences Centre, Manchester, UK 3 Centre for Endocrinology and Diabetes, Institute of Human Development, Faculty of Medical and Human Sciences, University of Manchester, UK 4 School of Biological Sciences, University of Auckland, New Zealand 5 Department of Pharmacology, Medical Sciences Division, University of Oxford, UK † Present address: School of Biosciences, University of Birmingham, UK *Joint corresponding authors: Natalie J. Gardiner and Garth J.S. Cooper Email: [email protected]; [email protected] Address: University of Manchester, AV Hill Building, Oxford Road, Manchester, M13 9PT, United Kingdom Telephone: +44 161 275 5768; +44 161 701 0240 Word count: 4,490 Number of tables: 1, Number of figures: 6 Running title: Metabolic dysfunction in diabetic neuropathy 1 Diabetes Publish Ahead of Print, published online October 15, 2015 Diabetes Page 2 of 255 Abstract High glucose levels in the peripheral nervous system (PNS) have been implicated in the pathogenesis of diabetic neuropathy (DN). However our understanding of the molecular mechanisms which cause the marked distal pathology is incomplete. Here we performed a comprehensive, system-wide analysis of the PNS of a rodent model of DN.
    [Show full text]
  • A SARS-Cov-2-Human Protein-Protein Interaction Map Reveals Drug Targets and Potential Drug-Repurposing
    A SARS-CoV-2-Human Protein-Protein Interaction Map Reveals Drug Targets and Potential Drug-Repurposing Supplementary Information Supplementary Discussion All SARS-CoV-2 protein and gene functions described in the subnetwork appendices, including the text below and the text found in the individual bait subnetworks, are based on the functions of homologous genes from other coronavirus species. These are mainly from SARS-CoV and MERS-CoV, but when available and applicable other related viruses were used to provide insight into function. The SARS-CoV-2 proteins and genes listed here were designed and researched based on the gene alignments provided by Chan et. al. 1 2020 .​ Though we are reasonably sure the genes here are well annotated, we want to note that not every ​ protein has been verified to be expressed or functional during SARS-CoV-2 infections, either in vitro or in vivo. ​ ​ In an effort to be as comprehensive and transparent as possible, we are reporting the sub-networks of these functionally unverified proteins along with the other SARS-CoV-2 proteins. In such cases, we have made notes within the text below, and on the corresponding subnetwork figures, and would advise that more caution be taken when examining these proteins and their molecular interactions. Due to practical limits in our sample preparation and data collection process, we were unable to generate data for proteins corresponding to Nsp3, Orf7b, and Nsp16. Therefore these three genes have been left out of the following literature review of the SARS-CoV-2 proteins and the protein-protein interactions (PPIs) identified in this study.
    [Show full text]
  • Supplementary Table 1: Adhesion Genes Data Set
    Supplementary Table 1: Adhesion genes data set PROBE Entrez Gene ID Celera Gene ID Gene_Symbol Gene_Name 160832 1 hCG201364.3 A1BG alpha-1-B glycoprotein 223658 1 hCG201364.3 A1BG alpha-1-B glycoprotein 212988 102 hCG40040.3 ADAM10 ADAM metallopeptidase domain 10 133411 4185 hCG28232.2 ADAM11 ADAM metallopeptidase domain 11 110695 8038 hCG40937.4 ADAM12 ADAM metallopeptidase domain 12 (meltrin alpha) 195222 8038 hCG40937.4 ADAM12 ADAM metallopeptidase domain 12 (meltrin alpha) 165344 8751 hCG20021.3 ADAM15 ADAM metallopeptidase domain 15 (metargidin) 189065 6868 null ADAM17 ADAM metallopeptidase domain 17 (tumor necrosis factor, alpha, converting enzyme) 108119 8728 hCG15398.4 ADAM19 ADAM metallopeptidase domain 19 (meltrin beta) 117763 8748 hCG20675.3 ADAM20 ADAM metallopeptidase domain 20 126448 8747 hCG1785634.2 ADAM21 ADAM metallopeptidase domain 21 208981 8747 hCG1785634.2|hCG2042897 ADAM21 ADAM metallopeptidase domain 21 180903 53616 hCG17212.4 ADAM22 ADAM metallopeptidase domain 22 177272 8745 hCG1811623.1 ADAM23 ADAM metallopeptidase domain 23 102384 10863 hCG1818505.1 ADAM28 ADAM metallopeptidase domain 28 119968 11086 hCG1786734.2 ADAM29 ADAM metallopeptidase domain 29 205542 11085 hCG1997196.1 ADAM30 ADAM metallopeptidase domain 30 148417 80332 hCG39255.4 ADAM33 ADAM metallopeptidase domain 33 140492 8756 hCG1789002.2 ADAM7 ADAM metallopeptidase domain 7 122603 101 hCG1816947.1 ADAM8 ADAM metallopeptidase domain 8 183965 8754 hCG1996391 ADAM9 ADAM metallopeptidase domain 9 (meltrin gamma) 129974 27299 hCG15447.3 ADAMDEC1 ADAM-like,
    [Show full text]