Analysis of Trans Esnps Infers Regulatory Network Architecture

Total Page:16

File Type:pdf, Size:1020Kb

Analysis of Trans Esnps Infers Regulatory Network Architecture Analysis of trans eSNPs infers regulatory network architecture Anat Kreimer Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in the Graduate School of Arts and Sciences COLUMBIA UNIVERSITY 2014 © 2014 Anat Kreimer All rights reserved ABSTRACT Analysis of trans eSNPs infers regulatory network architecture Anat Kreimer eSNPs are genetic variants associated with transcript expression levels. The characteristics of such variants highlight their importance and present a unique opportunity for studying gene regulation. eSNPs affect most genes and their cell type specificity can shed light on different processes that are activated in each cell. They can identify functional variants by connecting SNPs that are implicated in disease to a molecular mechanism. Examining eSNPs that are associated with distal genes can provide insights regarding the inference of regulatory networks but also presents challenges due to the high statistical burden of multiple testing. Such association studies allow: simultaneous investigation of many gene expression phenotypes without assuming any prior knowledge and identification of unknown regulators of gene expression while uncovering directionality. This thesis will focus on such distal eSNPs to map regulatory interactions between different loci and expose the architecture of the regulatory network defined by such interactions. We develop novel computational approaches and apply them to genetics-genomics data in human. We go beyond pairwise interactions to define network motifs, including regulatory modules and bi-fan structures, showing them to be prevalent in real data and exposing distinct attributes of such arrangements. We project eSNP associations onto a protein-protein interaction network to expose topological properties of eSNPs and their targets and highlight different modes of distal regulation. Overall, our work offers insights concerning the topological structure of human regulatory networks and the role genetics plays in shaping them. Tables of Contents List of Figures ................................................................................................................................. v List of Tables ............................................................................................................................... viii Acknowledgements ......................................................................................................................... x Chapter 1: Introduction ................................................................................................................... 1 Chapter 2: Inference of modules associated to eQTLs ................................................................... 7 2.1 Introduction ........................................................................................................................... 8 2.2 Results ................................................................................................................................. 13 2.2.1 Computational framework for detecting transcriptional modules ............................... 13 2.2.2 Modules’ topology ....................................................................................................... 16 2.2.3 Module’s score and filtering ........................................................................................ 18 2.2.4 Cis/trans-effects ........................................................................................................... 20 2.2.5 Independent cross validation by similar annotations from two sources of information and phenotypic analysis ........................................................................................................ 21 2.2.6 Comparison with standard approach to module construction ...................................... 23 2.2.7 Analysis of specific modules ....................................................................................... 24 2.3 Discussion ........................................................................................................................... 30 2.4 Materials and Methods ........................................................................................................ 32 i 2.4.1 Data details and processing .......................................................................................... 32 2.4.2 Step 1—nominal association testing ............................................................................ 32 2.4.3 Step 2—module construction, scoring and filtering .................................................... 34 2.4.4 Step 3—finding secondary SNPs ................................................................................. 36 2.4.5 Analysis of dependencies within modules ................................................................... 37 2.4.6 Module annotation ....................................................................................................... 39 2.4.7 Filtering modules using different criteria .................................................................... 39 2.4.8 Enrichment of cis-effects for main SNPs..................................................................... 40 Chapter 3: Co-regulated transcripts associated to cooperating eSNPs define bi-fan motifs in human gene networks ................................................................................................................... 41 3.1 Introduction ......................................................................................................................... 42 3.2 Results ................................................................................................................................. 45 3.2.1 Computational framework for associating pairs of SNPs with pairs of genes ............ 45 3.2.2 Distribution of genomic properties of eSNP sources and their gene targets ............... 52 3.2.3 Characterizing dependencies within cooperating quartets ........................................... 53 3.2.4 Identifying direction of effect between eSNP sources and gene targets ...................... 54 3.2.4 HLA quartet ................................................................................................................. 58 3.2.5 Functional enrichment of quartets ............................................................................... 60 3.2.6 Replication of quartet properties in a larger dataset .................................................... 60 3.3 Discussion ........................................................................................................................... 62 ii 3.4 Materials and Methods ........................................................................................................ 64 3.4.1 Data details and processing .......................................................................................... 64 3.4.2 Association testing ....................................................................................................... 64 3.4.3 Obtaining a random distribution of association test-statistics ..................................... 65 3.4.4 Creating and filtering quartets ..................................................................................... 65 3.4.5 Statistical challenges in comparing real vs. permuted quartets. .................................. 66 Chapter 4: Variants in exons and in transcription factors affect gene expression in trans ........... 67 4.1 Introduction ......................................................................................................................... 68 4.2 Results ................................................................................................................................. 72 4.2.1 Computational framework for mapping trans associations onto the PPI network ....... 72 4.2.2 Identifying topological properties of exonic eSNP interactions ............................. 74 4.2.3 Characterization of exon and transcription factor sources and targets ........................ 77 4.2.4 Co-expression of targets and cis-effects on the source gene ....................................... 82 4.2.5 Modular organization of eSNPs in TFs ....................................................................... 83 4.2.6 Support for eSNPs in TFs from different data sources ................................................ 85 4.2.7 Distribution of TF sources and targets in PPI functional clusters ............................... 86 4.2.8 Specific example of TF eSNP ...................................................................................... 86 4.2.9 Distribution of exonic sources and targets in PPI functional clusters ......................... 87 4.2.10 Specific example of exonic eSNP .............................................................................. 89 4.2.11 Mechanistic interpretation of exonic eSNPs .............................................................. 91 iii 4.3 Discussion ........................................................................................................................... 93 4.4 Materials and Methods ........................................................................................................ 99 4.3.1 Data details and processing .......................................................................................... 99 4.3.2 Association testing ....................................................................................................... 99 4.3.3 Obtaining a random
Recommended publications
  • Supplementary Information
    doi: 10.1038/nature08795 SUPPLEMENTARY INFORMATION Supplementary Discussion Population naming In some contexts, the indigenous hunter-gatherer and pastoralist peoples of southern Africa are referred to collectively as the Khoisan (Khoi-San) or more recently Khoesan (Khoe-San) people. This grouping is based on the unique linguistic use of click-consonants1. Many names, often country-specific, have been used by Bantu pastoralists and European settlers to describe the hunter-gatherers, including San, Saan, Sonqua, Soaqua, Souqua, Sanqua, Kwankhala, Basarwa, Batwa, Abathwa, Baroa, Bushmen, Bossiesmans, Bosjemans, or Bosquimanos. In addition, group-specific names such as !Kung and Khwe are often used for the broader population. The two most commonly used names, “San” and “Bushmen”, have both been associated with much controversy due to derogatory connotations2. “San” has become the more popular term used in Western literature, although “Bushmen” is arguably the more commonly recognized term within the communities. Since they have no collective name for themselves, the term Bushmen was selected for use in this paper as the term most familiar to the participants themselves. Regarding identification of individuals The five men identified in this study have all elected to have their identity made public knowledge. Thus we present two complete personal genomes (KB1 and ABT), a low-coverage personal genome (NB1), and personal exomes for all five men. On a scientific level, identification allows for current and future correlation of genetic data with demographic and medical histories. On a social level, identification allows for maximizing community benefit. For !Gubi, G/aq’o, D#kgao and !Aî, their name represents not only themselves, but importantly their extended family unit and a way of life severely under threat.
    [Show full text]
  • Inhibition of Breast and Brain Cancer Cell Growth by Bccipα, An
    Oncogene (2001) 20, 336 ± 345 ã 2001 Nature Publishing Group All rights reserved 0950 ± 9232/01 $15.00 www.nature.com/onc Inhibition of breast and brain cancer cell growth by BCCIPa,an evolutionarily conserved nuclear protein that interacts with BRCA2 Jingmei Liu1, Yuan Yuan1,2, Juan Huan2 and Zhiyuan Shen*,1 1Department of Molecular Genetics and Microbiology, University of New Mexico Health Sciences Center; 915 Camino de Salud, NE. Albuquerque, New Mexico, NM 87131, USA; 2Graduate Program of Molecular Genetics, College of Medicine, University of Illinois at Chicago, 900 S. Ashland Ave. Chicago, Illinois, IL 60607, USA BRCA2 is a tumor suppressor gene involved in mammary mouse BRCA2. It is expected that important functions tumorigenesis. Although important functions have been of BRCA2 reside in these conserved domains. Based on assigned to a few conserved domains of BRCA2, little is the functional analysis of the conserved BRCA2 known about the longest internal conserved domain domains, several models have been proposed for the encoded by exons 14 ± 24. We identi®ed a novel protein, role of BRCA2 in tumor suppression. designated BCCIPa, that interacts with part of the An N-terminus conserved domain in exon 3 (amino internal conserved region of human BRCA2. Human acids 48 ± 105) has been implicated in transcriptional BCCIP represents a family of proteins that are regulation of gene expression (Milner et al., 1997; evolutionarily conserved, and contain three distinct Nordling et al., 1998). Deletion of this region has been domains: an N-terminus acidic domain (NAD) of 30 ± identi®ed in breast cancers (Nordling et al., 1998).
    [Show full text]
  • Download.Cse.Ucsc.Edu/ Early Age of Onset (~2.5 Years) of This PRA Form in Goldenpath/Canfam2/Database/) Using Standard Settings
    Kropatsch et al. Canine Genetics and Epidemiology (2016) 3:7 DOI 10.1186/s40575-016-0037-x RESEARCH Open Access A large deletion in RPGR causes XLPRA in Weimaraner dogs Regina Kropatsch1*, Denis A. Akkad1, Matthias Frank2, Carsten Rosenhagen3, Janine Altmüller4,5, Peter Nürnberg4,6,7, Jörg T. Epplen1,8 and Gabriele Dekomien1 Abstract Background: Progressive retinal atrophy (PRA) belongs to a group of inherited retinal disorders associated with gradual vision impairment due to degeneration of retinal photoreceptors in various dog breeds. PRA is highly heterogeneous, with autosomal dominant, recessive or X-linked modes of inheritance. In this study we used exome sequencing to investigate the molecular genetic basis of a new type of PRA, which occurred spontaneously in a litter of German short-hair Weimaraner dogs. Results: Whole exome sequencing in two PRA-affected Weimaraner dogs identified a large deletion comprising the first four exons of the X-linked retinitis pigmentosa GTPase regulator (RPGR) gene known to be involved in human retinitis pigmentosa and canine PRA. Screening of 16 individuals in the corresponding pedigree of short-hair Weimaraners by qPCR, verified the deletion in hemizygous or heterozygous state in one male and six female dogs, respectively. The mutation was absent in 88 additional unrelated Weimaraners. The deletion was not detectable in the parents of one older female which transmitted the mutation to her offspring, indicating that the RPGR deletion represents a de novo mutation concerning only recent generations of the Weimaraner breed in Germany. Conclusion: Our results demonstrate the value of an existing DNA biobank combined with exome sequencing to identify the underlying genetic cause of a spontaneously occurring inherited disease.
    [Show full text]
  • Title a New Centrosomal Protein Regulates Neurogenesis By
    Title A new centrosomal protein regulates neurogenesis by microtubule organization Authors: Germán Camargo Ortega1-3†, Sven Falk1,2†, Pia A. Johansson1,2†, Elise Peyre4, Sanjeeb Kumar Sahu5, Loïc Broic4, Camino De Juan Romero6, Kalina Draganova1,2, Stanislav Vinopal7, Kaviya Chinnappa1‡, Anna Gavranovic1, Tugay Karakaya1, Juliane Merl-Pham8, Arie Geerlof9, Regina Feederle10,11, Wei Shao12,13, Song-Hai Shi12,13, Stefanie M. Hauck8, Frank Bradke7, Victor Borrell6, Vijay K. Tiwari§, Wieland B. Huttner14, Michaela Wilsch- Bräuninger14, Laurent Nguyen4 and Magdalena Götz1,2,11* Affiliations: 1. Institute of Stem Cell Research, Helmholtz Center Munich, German Research Center for Environmental Health, Munich, Germany. 2. Physiological Genomics, Biomedical Center, Ludwig-Maximilian University Munich, Germany. 3. Graduate School of Systemic Neurosciences, Biocenter, Ludwig-Maximilian University Munich, Germany. 4. GIGA-Neurosciences, Molecular regulation of neurogenesis, University of Liège, Belgium 5. Institute of Molecular Biology (IMB), Mainz, Germany. 6. Instituto de Neurociencias, Consejo Superior de Investigaciones Científicas and Universidad Miguel Hernández, Sant Joan d’Alacant, Spain. 7. Laboratory for Axon Growth and Regeneration, German Center for Neurodegenerative Diseases (DZNE), Bonn, Germany. 8. Research Unit Protein Science, Helmholtz Centre Munich, German Research Center for Environmental Health, Munich, Germany. 9. Protein Expression and Purification Facility, Institute of Structural Biology, Helmholtz Center Munich, German Research Center for Environmental Health, Munich, Germany. 10. Institute for Diabetes and Obesity, Monoclonal Antibody Core Facility, Helmholtz Center Munich, German Research Center for Environmental Health, Munich, Germany. 11. SYNERGY, Excellence Cluster of Systems Neurology, Biomedical Center, Ludwig- Maximilian University Munich, Germany. 12. Developmental Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, USA 13.
    [Show full text]
  • Whole-Genome Microarray Detects Deletions and Loss of Heterozygosity of Chromosome 3 Occurring Exclusively in Metastasizing Uveal Melanoma
    Anatomy and Pathology Whole-Genome Microarray Detects Deletions and Loss of Heterozygosity of Chromosome 3 Occurring Exclusively in Metastasizing Uveal Melanoma Sarah L. Lake,1 Sarah E. Coupland,1 Azzam F. G. Taktak,2 and Bertil E. Damato3 PURPOSE. To detect deletions and loss of heterozygosity of disease is fatal in 92% of patients within 2 years of diagnosis. chromosome 3 in a rare subset of fatal, disomy 3 uveal mela- Clinical and histopathologic risk factors for UM metastasis noma (UM), undetectable by fluorescence in situ hybridization include large basal tumor diameter (LBD), ciliary body involve- (FISH). ment, epithelioid cytomorphology, extracellular matrix peri- ϩ ETHODS odic acid-Schiff-positive (PAS ) loops, and high mitotic M . Multiplex ligation-dependent probe amplification 3,4 5 (MLPA) with the P027 UM assay was performed on formalin- count. Prescher et al. showed that a nonrandom genetic fixed, paraffin-embedded (FFPE) whole tumor sections from 19 change, monosomy 3, correlates strongly with metastatic death, and the correlation has since been confirmed by several disomy 3 metastasizing UMs. Whole-genome microarray analy- 3,6–10 ses using a single-nucleotide polymorphism microarray (aSNP) groups. Consequently, fluorescence in situ hybridization were performed on frozen tissue samples from four fatal dis- (FISH) detection of chromosome 3 using a centromeric probe omy 3 metastasizing UMs and three disomy 3 tumors with Ͼ5 became routine practice for UM prognostication; however, 5% years’ metastasis-free survival. to 20% of disomy 3 UM patients unexpectedly develop metas- tases.11 Attempts have therefore been made to identify the RESULTS. Two metastasizing UMs that had been classified as minimal region(s) of deletion on chromosome 3.12–15 Despite disomy 3 by FISH analysis of a small tumor sample were found these studies, little progress has been made in defining the key on MLPA analysis to show monosomy 3.
    [Show full text]
  • Dendritic Cell Differentiation Induced − Critical Role of AZI2 in GM-CSF
    Critical Role of AZI2 in GM-CSF−Induced Dendritic Cell Differentiation Masahiro Fukasaka, Daisuke Ori, Tatsukata Kawagoe, Satoshi Uematsu, Kenta Maruyama, Toshihiko Okazaki, This information is current as Tatsuya Kozaki, Tomoko Imamura, Sarang Tartey, Takashi of September 27, 2021. Mino, Takashi Satoh, Shizuo Akira and Osamu Takeuchi J Immunol published online 22 April 2013 http://www.jimmunol.org/content/early/2013/04/23/jimmun ol.1203155 Downloaded from Why The JI? Submit online. http://www.jimmunol.org/ • Rapid Reviews! 30 days* from submission to initial decision • No Triage! Every submission reviewed by practicing scientists • Fast Publication! 4 weeks from acceptance to publication *average by guest on September 27, 2021 Subscription Information about subscribing to The Journal of Immunology is online at: http://jimmunol.org/subscription Permissions Submit copyright permission requests at: http://www.aai.org/About/Publications/JI/copyright.html Email Alerts Receive free email-alerts when new articles cite this article. Sign up at: http://jimmunol.org/alerts The Journal of Immunology is published twice each month by The American Association of Immunologists, Inc., 1451 Rockville Pike, Suite 650, Rockville, MD 20852 Copyright © 2013 by The American Association of Immunologists, Inc. All rights reserved. Print ISSN: 0022-1767 Online ISSN: 1550-6606. Published April 24, 2013, doi:10.4049/jimmunol.1203155 The Journal of Immunology Critical Role of AZI2 in GM-CSF–Induced Dendritic Cell Differentiation Masahiro Fukasaka,*,† Daisuke Ori,*,‡,x Tatsukata Kawagoe,*,1 Satoshi Uematsu,*,2 Kenta Maruyama,* Toshihiko Okazaki,*,† Tatsuya Kozaki,* Tomoko Imamura,*,‡,x Sarang Tartey,*,‡,x Takashi Mino,‡,x Takashi Satoh,* Shizuo Akira,*,{ and Osamu Takeuchi*,‡,x TNFR-associated factor family member–associated NF-kB activator (TANK)–binding kinase 1 (TBK1) is critical for the activation of IFN regulatory factor 3 and type I IFN production upon virus infection.
    [Show full text]
  • The Molecular Genetic Architecture of Attention Deficit Hyperactivity Disorder
    Molecular Psychiatry (2015) 20, 289–297 © 2015 Macmillan Publishers Limited All rights reserved 1359-4184/15 www.nature.com/mp EXPERT REVIEW The molecular genetic architecture of attention deficit hyperactivity disorder Z Hawi1, TDR Cummins1, J Tong1, B Johnson1, R Lau1, W Samarrai2 and MA Bellgrove1 Attention deficit hyperactivity disorder (ADHD) is a common childhood behavioral condition which affects 2–10% of school age children worldwide. Although the underlying molecular mechanism for the disorder is poorly understood, familial, twin and adoption studies suggest a strong genetic component. Here we provide a state-of-the-art review of the molecular genetics of ADHD incorporating evidence from candidate gene and linkage designs, as well as genome-wide association (GWA) studies of common single-nucleotide polymorphisms (SNPs) and rare copy number variations (CNVs). Bioinformatic methods such as functional enrichment analysis and protein–protein network analysis are used to highlight biological processes of likely relevance to the aetiology of ADHD. Candidate gene associations of minor effect size have been replicated across a number of genes including SLC6A3, DRD5, DRD4, SLC6A4, LPHN3, SNAP-25, HTR1B, NOS1 and GIT1. Although case-control SNP-GWAS have had limited success in identifying common genetic variants for ADHD that surpass critical significance thresholds, quantitative trait designs suggest promising associations with Cadherin13 and glucose–fructose oxidoreductase domain 1 genes. Further, CNVs mapped to glutamate receptor genes (GRM1, GRM5, GRM7 and GRM8) have been implicated in the aetiology of the disorder and overlap with bioinformatic predictions based on ADHD GWAS SNP data regarding enriched pathways. Although increases in sample size across multi-center cohorts will likely yield important new results, we advocate that this must occur in parallel with a shift away from categorical case-control approaches that view ADHD as a unitary construct, towards dimensional approaches that incorporate endophenotypes and statistical classification methods.
    [Show full text]
  • Seq2pathway Vignette
    seq2pathway Vignette Bin Wang, Xinan Holly Yang, Arjun Kinstlick May 19, 2021 Contents 1 Abstract 1 2 Package Installation 2 3 runseq2pathway 2 4 Two main functions 3 4.1 seq2gene . .3 4.1.1 seq2gene flowchart . .3 4.1.2 runseq2gene inputs/parameters . .5 4.1.3 runseq2gene outputs . .8 4.2 gene2pathway . 10 4.2.1 gene2pathway flowchart . 11 4.2.2 gene2pathway test inputs/parameters . 11 4.2.3 gene2pathway test outputs . 12 5 Examples 13 5.1 ChIP-seq data analysis . 13 5.1.1 Map ChIP-seq enriched peaks to genes using runseq2gene .................... 13 5.1.2 Discover enriched GO terms using gene2pathway_test with gene scores . 15 5.1.3 Discover enriched GO terms using Fisher's Exact test without gene scores . 17 5.1.4 Add description for genes . 20 5.2 RNA-seq data analysis . 20 6 R environment session 23 1 Abstract Seq2pathway is a novel computational tool to analyze functional gene-sets (including signaling pathways) using variable next-generation sequencing data[1]. Integral to this tool are the \seq2gene" and \gene2pathway" components in series that infer a quantitative pathway-level profile for each sample. The seq2gene function assigns phenotype-associated significance of genomic regions to gene-level scores, where the significance could be p-values of SNPs or point mutations, protein-binding affinity, or transcriptional expression level. The seq2gene function has the feasibility to assign non-exon regions to a range of neighboring genes besides the nearest one, thus facilitating the study of functional non-coding elements[2]. Then the gene2pathway summarizes gene-level measurements to pathway-level scores, comparing the quantity of significance for gene members within a pathway with those outside a pathway.
    [Show full text]
  • CLASP2 Antibody Product Type
    PRODUCT INFORMATION Product name: CLASP2 antibody Product type: Primary antibodies Description: Rabbit polyclonal to CLASP2 Immunogen:3 synthetic peptides (human) conjugated to KLH Reacts with:Hu, Ms Tested applications:ELISA, WB and IF GENE INFORMATION Gene Symbol: CLASP2 Gene Name:cytoplasmic linker associated protein 2 Ensembl ID:ENSG00000163539 Entrez GeneID:23122 GenBank Accession number:AB014527 Swiss-Prot:O75122 Molecular weight of CLASP2: 165.9 & 108.6kDa Function:Microtubule plus-end tracking protein that promotes the stabilization of dynamic microtubules. Involved in the nucleation of noncentrosomal microtubules originating from the trans-Golgi network (TGN). Required for the polarization of the cytoplasmic microtubule arrays in migrating cells towards the leading edge of the cell. May act at the cell cortex to enhance the frequency of rescue of depolymerizing microtubules by attaching their plus- ends to cortical platforms composed of ERC1 and PHLDB2. This cortical microtubule stabilizing activity is regulated at least in part by phosphatidylinositol 3-kinase signaling. Also performs a similar stabilizing function at the kinetochore which is essential for the bipolar alignment of chromosomes on the mitotic spindle. Acts as a mediator of ERBB2- dependent stabilization of microtubules at the cell cortex. Expected subcellular localization:Cytoplasm › cytoskeleton. Cytoplasm › cytoskeleton › microtubule organizing center › centrosome. Chromosome › centromere › kinetochore. Cytoplasm › cytoskeleton › spindle. Golgi apparatus. Golgi apparatus › trans-Golgi network. Cell membrane. Cell projection › ruffle membrane. Note: Localizes to microtubule plus ends. Localizes to centrosomes, kinetochores and the mitotic spindle from prometaphase. Subsequently localizes to the spindle midzone from anaphase and to the midbody from telophase. In migrating cells localizes to the plus ends of microtubules within the cell body and to the entire microtubule lattice within the lamella.
    [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]
  • Is BOK Required for Apoptosis Induced by Endoplasmic Reticulum Stress?
    LETTER Is BOK required for apoptosis induced by endoplasmic reticulum stress? LETTER Yuniel Fernandez-Marreroa, Francine Keb,c, Nohemy Echeverrya,1, Philippe Bouilletb,c, Daniel Bachmanna, Andreas Strasserb,c, and Thomas Kaufmanna,2 The B-cell lymphoma 2 (BCL-2)-related ovarian killer comparable Chop levels). Given the critical role of BIM (BOK) shares sequence homology with the proapo- in ER stress-induced apoptosis, this reduction of Bim ptotic BCL-2 family members BAX and BAK. However, may fully account for the reported resistance to ER −/− Bok cells are not protected from classic apoptotic stress. It is unclear whether this reduction of Bim is par- triggers and evidence for a proapoptotic role of BOK ticular to these SV40 MEFs, which are prone to line-to- is derived mostly from overexpression studies (1). BOK line variations within the one genotype, or whether this localizes preferentially to the endoplasmic reticulum is also seen in primary cells (e.g., primary MEFs, which (ER) membrane, where it interacts with IP3-receptors were used for some experiments) from these mice. Im- −/− (2, 3). Using cells from their newly generated Bok portantly, we did not observe significant changes in Bim −/− mouse strain, Carpio et al. propose that BOK is a crit- levels in SV40 MEFs or tissues from our Bok mice (Fig. ical inducer of BAX/BAK-dependent apoptosis in re- 1A). Overall, our analysis of SV40 MEFs, primary MEFs, sponse to ER stress (4). This proposal is in contrast to myeloid progenitors, mast cells, and primary neutrophils our earlier report, in which we showed that loss of BOK did not support a proapoptotic role of BOK downstream did not confer resistance toward ER stress in several of ER stress (2) (Fig.
    [Show full text]
  • Statistical Methods for High-Dimensional Networked Data Analysis
    Statistical Methods for High-Dimensional Networked Data Analysis by Yan Zhou A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Biostatistics) in The University of Michigan 2015 Doctoral Committee: Professor Peter X. K. Song, Chair Professor Matthias Kretzler Assistant Professor Xiaoquan William Wen Professor Ji Zhu c Yan Zhou 2015 All Rights Reserved To my family ii ACKNOWLEDGEMENTS I would like to acknowledge many people for their help throughout my graduate study. I would first like to express my deepest gratitude to my advisor Dr. Peter X.- K. Song, for sharing his knowledge, encouraging me to work hard and constantly try to improve my work, and giving me the freedom to explore a diverse set of projects. He is the one who first sparked my interest in networked data analysis. I am not only respectful for his immense knowledge, motivation, and enthusiasm towards research, but also appreciate his invaluable guidance and enthusiastic encouragement during many difficult moments in my research. I could not have imagined having a better advisor and mentor for my Ph.D. study. I would also like to thank my other thesis committee members Dr. Ji Zhu, Dr. Xiaoquan Willian Wen and Dr. Matthias Kretzler for their help and advices that made me more productive than what I could achieve on my own. All enlightening discussions I had with them have lead to better contents of this dissertation. In addition, I would like to acknowledge my collaborators, Dr. Pei Wang, Dr. Betsy Lozoff and Dr. Fengji Geng, who offered me excellent opportunities and con- structive suggestions to apply my statistical knowledge to solve real-world problems, which greatly strengthened the scientific background of my dissertation research.
    [Show full text]