A Tumorigenic Index for Quantitative Analysis of Liver Cancer Initiation and Progression

Total Page:16

File Type:pdf, Size:1020Kb

A Tumorigenic Index for Quantitative Analysis of Liver Cancer Initiation and Progression A tumorigenic index for quantitative analysis of liver cancer initiation and progression Gaowei Wanga, Xiaolin Luoa, Yan Lianga, Kota Kanekoa, Hairi Lib, Xiang-Dong Fub, and Gen-Sheng Fenga,1 aDepartment of Pathology, Division of Biological Sciences, Moores Cancer Center, University of California San Diego, La Jolla, CA 92093; and bDepartment of Cellular and Molecular Medicine, University of California San Diego, La Jolla, CA 92093 Edited by Roger J. Davis, University of Massachusetts Medical School, Worcester, MA, and approved November 13, 2019 (received for review June 29, 2019) Primary liver cancer develops from multifactorial etiologies, result- prognosis of liver cancer patients with diverse etiologies and also in ing in extensive genomic heterogeneity. To probe the common diagnosis of precancer patients with steatosis, fibrosis, or cirrhosis. mechanism of hepatocarcinogenesis, we interrogated temporal Upon further development and optimization, this TI platform may gene expression profiles in a group of mouse models with hepatic become a quantitative means for early detection of liver tumor steatosis, fibrosis, inflammation, and, consequently, tumorigenesis. initiation or risk assessment of chronic disease patients, based on Instead of anticipated progressive changes, we observed a sudden comprehensive analysis of the whole transcriptome rather than a molecular switch at a critical precancer stage, by developing analyt- few biomarker molecules. We believe that the principle and ra- ical platform that focuses on transcription factor (TF) clusters. tionale of this TI derivation method can be extended from HCC to Coarse-grained network modeling demonstrated that an abrupt other types of cancer, to develop a quantitative analytical tool for transcriptomic transition occurred once changes were accumulated cancer prediction and prognosis in a given organ or tissue. to reach a threshold. Based on the experimental and bioinformatic data analyses as well as mathematical modeling, we derived a tu- Results morigenic index (TI) to quantify tumorigenic signal strengths. The TI Temporal Gene Expression Patterns in Liver Tumorigenesis. Pten is a is powerful in predicting the disease status of patients with meta- tumor suppressor that counteracts the PI3K/Akt signaling path- bolic disorders and also the tumor stages and prognosis of liver way, and Pten deficiency has often been detected in liver can- cancer patients with diverse backgrounds. This work establishes a cer patients (5, 6). Targeted deletion of Pten in hepatocytes quantitative tool for triage of liver cancer patients and also for can- induced NAFLD and, subsequently, NASH, followed by tumor cer risk assessment of chronic liver disease patients. development in mice (7, 8). The NASH-driven pathogenic pro- MEDICAL SCIENCES cess in Pten-deficient liver was accelerated and aggravated by liver cancer | tumorigenic index | quantitative analysis | additional deletion of Shp2/Ptpn11 (9), a newly identified liver transcription factor clusters tumor suppressor (10, 11). These mutant mouse lines with de- fined genetic background and characterized tumor phenotype he incidences and mortality of liver cancer, mainly hepato- constitute an ideal group of animal models to dissect stepwise Tcellular carcinoma (HCC), are increasing rapidly worldwide mechanisms underlying NASH-HCC development. Toward this (1). Diverse risk factors for primary liver cancer have been goal, we isolated liver samples at multiple time points from mice identified, including infection of hepatitis B virus and hepatitis C with hepatocyte-specific deletion of Pten (PKO), Shp2 (SKO), virus, alcohol abuse, nonalcoholic fatty liver disease (NAFLD) or nonalcoholic steatohepatitis (NASH), and intake of aflatoxin B1. Significance Consistent with the complex and multifactorial etiologies, multiomics analyses of human liver tumor samples have identified The mechanisms of hepatocarcinogenesis are poorly understood. vast genomic heterogeneity, molecular and cellular defects, In this study, we interrogated the temporal gene expression metabolic reprogramming, and subtypes of tumors as well as profiles in mouse livers during progression of chronic fatty liver altered tumor microenvironment in the liver (2–4). diseases to carcinogenesis. By establishing an analytical system However, it remains to be determined if any common molecular that focuses on transcription factor (TF) clusters, we identified a signatures in the transcriptomes exist for liver cancer, despite their sudden switch from healthy liver to tumor tissues in the tran- considerable genomic heterogeneity. Furthermore, little is known scriptomes at precancer stage, prior to detection of any tumor about the kinetics and fashions, either gradual accumulation or nodule. We further developed a platform to calculate tumorigenic dramatic transition, in generation of cell-intrinsic and cell-extrinsic index (TI) based on the transcriptome, especially the TF clusters, to signals that are intertwined to drive malignant transformation of measure tumorigenic signal strengths. This quantitative analytical hepatocytes and tumor initiation. To dissect the stepwise onco- tool of TI is powerful in assessing cancer risk of chronic liver dis- genic signals and mechanisms at the precancer stages, we chose ease patients and in predicting tumor stages and prognosis of to work on mouse models that recapitulate key pathogenic fea- liver cancer patients. tures in human liver cancer. By pathological examination, RNA- sequencing (RNA-seq), and bioinformatic data analysis, we iden- Author contributions: G.-S.F. designed research; G.W., X.L., and G.-S.F. performed re- search; G.W., X.L., K.K., H.L., X.-D.F., and G.-S.F. contributed new reagents/analytic tools; tified a sudden switch in transcription factor (TF) clusters at a G.W., X.L., Y.L., and G.-S.F. analyzed data; and G.W. and G.-S.F. wrote the paper. precancer stage of hepatocarcinogenesis. Based on a multilayer The authors declare no competing interest. analysis of the TF clusters and mathematic modeling, we de- This article is a PNAS Direct Submission. veloped a tumorigenic index (TI) calculation system for quantita- tive measurement of liver tumorigenesis. Although this platform Published under the PNAS license. was established based on the transcriptomic data derived from Data deposition: The RNA-seq data have been deposited in the National Center for Bio- technology Information Gene Expression Omnibus (NCBI GEO) database under ID code genetically modified mouse tumor models, we have found appli- GEO: GSE123427. Codes have been deposited in GitHub (https://github.com/wanyewang1/ cable effectiveness of the derived TI values in determination of Index_model). disease status in other mouse models of chronic or malignant liver 1To whom correspondence may be addressed. Email: [email protected]. diseases with diverse backgrounds. Furthermore, using the TI This article contains supporting information online at https://www.pnas.org/lookup/suppl/ method developed from animal models to interrogate the hu- doi:10.1073/pnas.1911193116/-/DCSupplemental. man patients’ data, we demonstrated the power of the TI tool in www.pnas.org/cgi/doi/10.1073/pnas.1911193116 PNAS Latest Articles | 1of8 Downloaded by guest on September 24, 2021 Pten and Shp2 (DKO), and wild-type (WT) control (Fig. 1A). point, we identified genes that were significantly changed in Based on the developmental or pathological features, these liver mutant livers at different stages (Dataset S1). The resulting samples were divided into 4 groups: 1) youth, 2) adult, 3) pre- heatmap depicts the progressive changes in liver transcriptomes cancer, and 4) cancer (Fig. 1A). Representative histological or of the 3 mutant mouse lines (Fig. 1B), which correlated well with pathological properties of the specimens were shown (SI Ap- the kinetics and severity of tumor progression in SKO, PKO, and pendix, Fig. S1) and were also described in detail previously (9). DKO livers (SI Appendix, Fig. S1). Highlighted in Fig. 1B are RNA-seq was performed for all liver samples, and the data some significantly changed genes related to the MAPK pathway, quality analysis showed high levels of correlation for biological fatty acid, glucose or bile acid metabolism, extracellular matrix, replicates in each group (SI Appendix, Fig. S2A), with reduced epigenetic machinery, and inflammatory response. expression of Shp2 and/or Pten in corresponding mutants as By comparing the transcriptomes between PKO and WT livers expected (SI Appendix, Fig. S2 B and C). By comparing expres- at multiple time points, we identified significantly changed bi- sion levels in each mutant with WT control at the same time ological processes in the mutant liver using gene set enrichment A B 1 Pkm Jun Time (month) 0.5 0 Mapk3 1 2 3 4 5 7 912 16 −0.5 Mmp2 WT −1 Abhd5 PKO Tnf SKO 1M 2M 3M 4M 5M 7M_T 12M_T 16M_T 16M_T 1M 2M 7M 12M 1M 2M 7M_T 12M_T T T T DKO Kmt2c Kat6a Slc27a2 Youth Adult Pre-cancer Cancer Cyp7b1 PKO SKO DKO CDE F 1 1 1 Extracellular JAG2_NOTCH1 1 TRIGLYCERIDE_BIOSYNTHESIS BMP4_ACVR2B BMP7_ACVR2B Hdac5 EXTRACELLULAR_MATRIX_ORGANIZATION Structure CSF2_CSF2RA 0.5 EFNB3_EPHB2 INTEGRIN1_PATHWAY PDGFB_PDGFRB 0.5 0.5 0.5 ECM_RECEPTOR_INTERACTION Fatty acid FGF9_FGFR3 EFNA5_EPHA8 Mta3 ECM_AFFILIATED 0 TGFB3_TGFBR2 metabolism CSF1_CSF1R AVB3_INTEGRIN_PATHWAY IGF2_IGF2R Top2a JAG1_NOTCH4 −0.5 0 0 0 COLLAGEN_FORMATION Inflammatory MDK_PTPRB TNFSF12_TNFRSF12 CORE_MATRISOME NTF5_NTRK2 Dnmt1 response ECM_GLYCOPROTEINS −1 IL5_CSF2RB IFNG_IFNGR1
Recommended publications
  • THE ROLE of HOMEODOMAIN TRANSCRIPTION FACTOR IRX5 in CARDIAC CONTRACTILITY and HYPERTROPHIC RESPONSE by © COPYRIGHT by KYOUNG H
    THE ROLE OF HOMEODOMAIN TRANSCRIPTION FACTOR IRX5 IN CARDIAC CONTRACTILITY AND HYPERTROPHIC RESPONSE By KYOUNG HAN KIM A THESIS SUBMITTED IN CONFORMITY WITH THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY GRADUATE DEPARTMENT OF PHYSIOLOGY UNIVERSITY OF TORONTO © COPYRIGHT BY KYOUNG HAN KIM (2011) THE ROLE OF HOMEODOMAIN TRANSCRIPTION FACTOR IRX5 IN CARDIAC CONTRACTILITY AND HYPERTROPHIC RESPONSE KYOUNG HAN KIM DOCTOR OF PHILOSOPHY GRADUATE DEPARTMENT OF PHYSIOLOGY UNIVERSITY OF TORONTO 2011 ABSTRACT Irx5 is a homeodomain transcription factor that negatively regulates cardiac fast transient + outward K currents (Ito,f) via the KV4.2 gene and is thereby a major determinant of the transmural repolarization gradient. While Ito,f is invariably reduced in heart disease and changes in Ito,f can modulate both cardiac contractility and hypertrophy, less is known about a functional role of Irx5, and its relationship with Ito,f, in the normal and diseased heart. Here I show that Irx5 plays crucial roles in the regulation of cardiac contractility and proper adaptive hypertrophy. Specifically, Irx5-deficient (Irx5-/-) hearts had reduced cardiac contractility and lacked the normal regional difference in excitation-contraction with decreased action potential duration, Ca2+ transients and myocyte shortening in sub-endocardial, but not sub-epicardial, myocytes. In addition, Irx5-/- mice showed less cardiac hypertrophy, but increased interstitial fibrosis and greater contractility impairment following pressure overload. A defect in hypertrophic responses in Irx5-/- myocardium was confirmed in cultured neonatal mouse ventricular myocytes, exposed to norepinephrine while being restored with Irx5 replacement. Interestingly, studies using mice ii -/- virtually lacking Ito,f (i.e. KV4.2-deficient) showed that reduced contractility in Irx5 mice was completely restored by loss of KV4.2, whereas hypertrophic responses to pressure-overload in hearts remained impaired when both Irx5 and Ito,f were absent.
    [Show full text]
  • Relb Deficiency in Dendritic Cells Protects from Autoimmune
    RelB Deficiency in Dendritic Cells Protects from Autoimmune Inflammation Due to Spontaneous Accumulation of Tissue T Regulatory Cells This information is current as of September 24, 2021. Nico Andreas, Maria Potthast, Anna-Lena Geiselhöringer, Garima Garg, Renske de Jong, Julia Riewaldt, Dennis Russkamp, Marc Riemann, Jean-Philippe Girard, Simon Blank, Karsten Kretschmer, Carsten Schmidt-Weber, Thomas Korn, Falk Weih and Caspar Ohnmacht Downloaded from J Immunol 2019; 203:2602-2613; Prepublished online 2 October 2019; doi: 10.4049/jimmunol.1801530 http://www.jimmunol.org/content/203/10/2602 http://www.jimmunol.org/ Supplementary http://www.jimmunol.org/content/suppl/2019/10/01/jimmunol.180153 Material 0.DCSupplemental References This article cites 74 articles, 23 of which you can access for free at: http://www.jimmunol.org/content/203/10/2602.full#ref-list-1 by guest on September 24, 2021 Why The JI? Submit online. • 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 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 Author Choice Freely available online through The Journal of Immunology Author Choice option 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 © 2019 by The American Association of Immunologists, Inc.
    [Show full text]
  • Identifying and Mapping Cell-Type-Specific Chromatin PNAS PLUS Programming of Gene Expression
    Identifying and mapping cell-type-specific chromatin PNAS PLUS programming of gene expression Troels T. Marstranda and John D. Storeya,b,1 aLewis-Sigler Institute for Integrative Genomics, and bDepartment of Molecular Biology, Princeton University, Princeton, NJ 08544 Edited by Wing Hung Wong, Stanford University, Stanford, CA, and approved January 2, 2014 (received for review July 2, 2013) A problem of substantial interest is to systematically map variation Relating DHS to gene-expression levels across multiple cell in chromatin structure to gene-expression regulation across con- types is challenging because the DHS represents a continuous ditions, environments, or differentiated cell types. We developed variable along the genome not bound to any specific region, and and applied a quantitative framework for determining the exis- the relationship between DHS and gene expression is largely tence, strength, and type of relationship between high-resolution uncharacterized. To exploit variation across cell types and test chromatin structure in terms of DNaseI hypersensitivity and genome- for cell-type-specific relationships between DHS and gene expres- wide gene-expression levels in 20 diverse human cell types. We sion, the measurement units must be placed on a common scale, show that ∼25% of genes show cell-type-specific expression ex- the continuous DHS measure associated to each gene in a well- plained by alterations in chromatin structure. We find that distal defined manner, and all measurements considered simultaneously. regions of chromatin structure (e.g., ±200 kb) capture more genes Moreover, the chromatin and gene-expression relationship may with this relationship than local regions (e.g., ±2.5 kb), yet the local only manifest in a single cell type, making standard measures of regions show a more pronounced effect.
    [Show full text]
  • Molecular Profile of Tumor-Specific CD8+ T Cell Hypofunction in a Transplantable Murine Cancer Model
    Downloaded from http://www.jimmunol.org/ by guest on September 25, 2021 T + is online at: average * The Journal of Immunology , 34 of which you can access for free at: 2016; 197:1477-1488; Prepublished online 1 July from submission to initial decision 4 weeks from acceptance to publication 2016; doi: 10.4049/jimmunol.1600589 http://www.jimmunol.org/content/197/4/1477 Molecular Profile of Tumor-Specific CD8 Cell Hypofunction in a Transplantable Murine Cancer Model Katherine A. Waugh, Sonia M. Leach, Brandon L. Moore, Tullia C. Bruno, Jonathan D. Buhrman and Jill E. Slansky J Immunol cites 95 articles Submit online. Every submission reviewed by practicing scientists ? is published twice each month by Receive free email-alerts when new articles cite this article. Sign up at: http://jimmunol.org/alerts http://jimmunol.org/subscription Submit copyright permission requests at: http://www.aai.org/About/Publications/JI/copyright.html http://www.jimmunol.org/content/suppl/2016/07/01/jimmunol.160058 9.DCSupplemental This article http://www.jimmunol.org/content/197/4/1477.full#ref-list-1 Information about subscribing to The JI No Triage! Fast Publication! Rapid Reviews! 30 days* Why • • • Material References Permissions Email Alerts Subscription Supplementary The Journal of Immunology The American Association of Immunologists, Inc., 1451 Rockville Pike, Suite 650, Rockville, MD 20852 Copyright © 2016 by The American Association of Immunologists, Inc. All rights reserved. Print ISSN: 0022-1767 Online ISSN: 1550-6606. This information is current as of September 25, 2021. The Journal of Immunology Molecular Profile of Tumor-Specific CD8+ T Cell Hypofunction in a Transplantable Murine Cancer Model Katherine A.
    [Show full text]
  • Feedback Regulation Between Initiation and Maturation Networks Orchestrates the Chromatin Dynamics of Epidermal Lineage
    bioRxiv preprint doi: https://doi.org/10.1101/349308; this version posted June 18, 2018. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Li et al., p. 1 Feedback Regulation between Initiation and Maturation Networks Orchestrates the Chromatin Dynamics of Epidermal Lineage Commitment Lingjie Li1,3,4, Yong Wang2,4,7,8*, Jessica L. Torkelson1,3*, Gautam Shankar1, Jillian M. Pattison1,3, Hanson H. Zhen1,3, Zhana Duren2,4,7, Fengqin Fang5, Sandra P. Melo1, Samantha N. Piekos1,3, Jiang Li1, Eric J. Liaw1, Lang Chen7, Rui Li1,4, Marius Wernig6, Wing H. Wong2,4, Howard Y. Chang1,4, Anthony E. Oro1,3,9 1 Program in Epithelial Biology and Department of Dermatology 2 Department of Statistics and Biomedical Data Science 3 Center for Definitive and Curative Medicine 4 Center for Personal Dynamic Regulome 5 Division of Immunology and Rheumatology, Department of Medicine, 6 Institute for Stem Cell Biology and Regenerative Medicine, Department of Pathology, Stanford University School of Medicine, Stanford, CA 94305, USA. 7 CEMS, NCMIS, MDIS, Academy of Mathematics & Systems Science, Chinese Academy of Sciences, Beijing,100080, China 8 Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, 650223, China *These authors made equal and independent contributions. 9 Correspondence to Lead Contact: Anthony E. Oro at [email protected] bioRxiv preprint doi: https://doi.org/10.1101/349308; this version posted June 18, 2018.
    [Show full text]
  • Modes of Interaction of KMT2 Histone H3 Lysine 4 Methyltransferase/COMPASS Complexes with Chromatin
    cells Review Modes of Interaction of KMT2 Histone H3 Lysine 4 Methyltransferase/COMPASS Complexes with Chromatin Agnieszka Bochy ´nska,Juliane Lüscher-Firzlaff and Bernhard Lüscher * ID Institute of Biochemistry and Molecular Biology, Medical School, RWTH Aachen University, Pauwelsstrasse 30, 52057 Aachen, Germany; [email protected] (A.B.); jluescher-fi[email protected] (J.L.-F.) * Correspondence: [email protected]; Tel.: +49-241-8088850; Fax: +49-241-8082427 Received: 18 January 2018; Accepted: 27 February 2018; Published: 2 March 2018 Abstract: Regulation of gene expression is achieved by sequence-specific transcriptional regulators, which convey the information that is contained in the sequence of DNA into RNA polymerase activity. This is achieved by the recruitment of transcriptional co-factors. One of the consequences of co-factor recruitment is the control of specific properties of nucleosomes, the basic units of chromatin, and their protein components, the core histones. The main principles are to regulate the position and the characteristics of nucleosomes. The latter includes modulating the composition of core histones and their variants that are integrated into nucleosomes, and the post-translational modification of these histones referred to as histone marks. One of these marks is the methylation of lysine 4 of the core histone H3 (H3K4). While mono-methylation of H3K4 (H3K4me1) is located preferentially at active enhancers, tri-methylation (H3K4me3) is a mark found at open and potentially active promoters. Thus, H3K4 methylation is typically associated with gene transcription. The class 2 lysine methyltransferases (KMTs) are the main enzymes that methylate H3K4. KMT2 enzymes function in complexes that contain a necessary core complex composed of WDR5, RBBP5, ASH2L, and DPY30, the so-called WRAD complex.
    [Show full text]
  • I LITERATURE-BASED DISCOVERY of KNOWN and POTENTIAL NEW
    LITERATURE-BASED DISCOVERY OF KNOWN AND POTENTIAL NEW MECHANISMS FOR RELATING THE STATUS OF CHOLESTEROL TO THE PROGRESSION OF BREAST CANCER BY YU WANG THESIS Submitted in partial fulfillment of the requirements for the degree of Master of Science in Bioinformatics with a concentration in Library and Information Science in the Graduate College of the University of Illinois at Urbana-Champaign, 2019 Urbana, Illinois Adviser: Professor Vetle I. Torvik Professor Erik Russell Nelson i ABSTRACT Breast cancer has been studied for a long period of time and from a variety of perspectives in order to understand its pathogeny. The pathogeny of breast cancer can be classified into two groups: hereditary and spontaneous. Although cancer in general is considered a genetic disease, spontaneous factors are responsible for most of the pathogeny of breast cancer. In other words, breast cancer is more likely to be caused and deteriorated by the dysfunction of a physical molecule than be caused by germline mutation directly. Interestingly, cholesterol, as one of those molecules, has been discovered to correlate with breast cancer risk. However, the mechanisms of how cholesterol helps breast cancer progression are not thoroughly understood. As a result, this study aims to study known and discover potential new mechanisms regarding to the correlation of cholesterol and breast cancer progression using literature review and literature-based discovery. The known mechanisms are further classified into four groups: cholesterol membrane content, transport of cholesterol, cholesterol metabolites, and other. The potential mechanisms, which are intended to provide potential new treatments, have been identified and checked for feasibility by an expert.
    [Show full text]
  • Watsonjn2018.Pdf (1.780Mb)
    UNIVERSITY OF CENTRAL OKLAHOMA Edmond, Oklahoma Department of Biology Investigating Differential Gene Expression in vivo of Cardiac Birth Defects in an Avian Model of Maternal Phenylketonuria A THESIS SUBMITTED TO THE GRADUATE FACULTY In partial fulfillment of the requirements For the degree of MASTER OF SCIENCE IN BIOLOGY By Jamie N. Watson Edmond, OK June 5, 2018 J. Watson/Dr. Nikki Seagraves ii J. Watson/Dr. Nikki Seagraves Acknowledgements It is difficult to articulate the amount of gratitude I have for the support and encouragement I have received throughout my master’s thesis. Many people have added value and support to my life during this time. I am thankful for the education, experience, and friendships I have gained at the University of Central Oklahoma. First, I would like to thank Dr. Nikki Seagraves for her mentorship and friendship. I lucked out when I met her. I have enjoyed working on this project and I am very thankful for her support. I would like thank Thomas Crane for his support and patience throughout my master’s degree. I would like to thank Dr. Shannon Conley for her continued mentorship and support. I would like to thank Liz Bullen and Dr. Eric Howard for their training and help on this project. I would like to thank Kristy Meyer for her friendship and help throughout graduate school. I would like to thank my committee members Dr. Robert Brennan and Dr. Lilian Chooback for their advisement on this project. Also, I would like to thank the biology faculty and staff. I would like to thank the Seagraves lab members: Jailene Canales, Kayley Pate, Mckayla Muse, Grace Thetford, Kody Harvey, Jordan Guffey, and Kayle Patatanian for their hard work and support.
    [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]
  • 4-6 Weeks Old Female C57BL/6 Mice Obtained from Jackson Labs Were Used for Cell Isolation
    Methods Mice: 4-6 weeks old female C57BL/6 mice obtained from Jackson labs were used for cell isolation. Female Foxp3-IRES-GFP reporter mice (1), backcrossed to B6/C57 background for 10 generations, were used for the isolation of naïve CD4 and naïve CD8 cells for the RNAseq experiments. The mice were housed in pathogen-free animal facility in the La Jolla Institute for Allergy and Immunology and were used according to protocols approved by the Institutional Animal Care and use Committee. Preparation of cells: Subsets of thymocytes were isolated by cell sorting as previously described (2), after cell surface staining using CD4 (GK1.5), CD8 (53-6.7), CD3ε (145- 2C11), CD24 (M1/69) (all from Biolegend). DP cells: CD4+CD8 int/hi; CD4 SP cells: CD4CD3 hi, CD24 int/lo; CD8 SP cells: CD8 int/hi CD4 CD3 hi, CD24 int/lo (Fig S2). Peripheral subsets were isolated after pooling spleen and lymph nodes. T cells were enriched by negative isolation using Dynabeads (Dynabeads untouched mouse T cells, 11413D, Invitrogen). After surface staining for CD4 (GK1.5), CD8 (53-6.7), CD62L (MEL-14), CD25 (PC61) and CD44 (IM7), naïve CD4+CD62L hiCD25-CD44lo and naïve CD8+CD62L hiCD25-CD44lo were obtained by sorting (BD FACS Aria). Additionally, for the RNAseq experiments, CD4 and CD8 naïve cells were isolated by sorting T cells from the Foxp3- IRES-GFP mice: CD4+CD62LhiCD25–CD44lo GFP(FOXP3)– and CD8+CD62LhiCD25– CD44lo GFP(FOXP3)– (antibodies were from Biolegend). In some cases, naïve CD4 cells were cultured in vitro under Th1 or Th2 polarizing conditions (3, 4).
    [Show full text]
  • Supplemental Materials ZNF281 Enhances Cardiac Reprogramming
    Supplemental Materials ZNF281 enhances cardiac reprogramming by modulating cardiac and inflammatory gene expression Huanyu Zhou, Maria Gabriela Morales, Hisayuki Hashimoto, Matthew E. Dickson, Kunhua Song, Wenduo Ye, Min S. Kim, Hanspeter Niederstrasser, Zhaoning Wang, Beibei Chen, Bruce A. Posner, Rhonda Bassel-Duby and Eric N. Olson Supplemental Table 1; related to Figure 1. Supplemental Table 2; related to Figure 1. Supplemental Table 3; related to the “quantitative mRNA measurement” in Materials and Methods section. Supplemental Table 4; related to the “ChIP-seq, gene ontology and pathway analysis” and “RNA-seq” and gene ontology analysis” in Materials and Methods section. Supplemental Figure S1; related to Figure 1. Supplemental Figure S2; related to Figure 2. Supplemental Figure S3; related to Figure 3. Supplemental Figure S4; related to Figure 4. Supplemental Figure S5; related to Figure 6. Supplemental Table S1. Genes included in human retroviral ORF cDNA library. Gene Gene Gene Gene Gene Gene Gene Gene Symbol Symbol Symbol Symbol Symbol Symbol Symbol Symbol AATF BMP8A CEBPE CTNNB1 ESR2 GDF3 HOXA5 IL17D ADIPOQ BRPF1 CEBPG CUX1 ESRRA GDF6 HOXA6 IL17F ADNP BRPF3 CERS1 CX3CL1 ETS1 GIN1 HOXA7 IL18 AEBP1 BUD31 CERS2 CXCL10 ETS2 GLIS3 HOXB1 IL19 AFF4 C17ORF77 CERS4 CXCL11 ETV3 GMEB1 HOXB13 IL1A AHR C1QTNF4 CFL2 CXCL12 ETV7 GPBP1 HOXB5 IL1B AIMP1 C21ORF66 CHIA CXCL13 FAM3B GPER HOXB6 IL1F3 ALS2CR8 CBFA2T2 CIR1 CXCL14 FAM3D GPI HOXB7 IL1F5 ALX1 CBFA2T3 CITED1 CXCL16 FASLG GREM1 HOXB9 IL1F6 ARGFX CBFB CITED2 CXCL3 FBLN1 GREM2 HOXC4 IL1F7
    [Show full text]
  • Olig1 and Sox10 Interact Synergistically to Drivemyelin Basic
    The Journal of Neuroscience, December 26, 2007 • 27(52):14375–14382 • 14375 Cellular/Molecular Olig1 and Sox10 Interact Synergistically to Drive Myelin Basic Protein Transcription in Oligodendrocytes Huiliang Li,1 Yan Lu,2 Hazel K. Smith,1 and William D. Richardson1 1Wolfson Institute for Biomedical Research and Department of Biology, University College London, London WC1E 6BT, United Kingdom, and 2Medical Research Council, Clinical Sciences Centre, Imperial College London, London W12 0NN, United Kingdom The oligodendrocyte lineage genes (Olig1/2), encoding basic helix-loop-helix transcription factors, were first identified in screens for master regulators of oligodendrocyte development. OLIG1 is important for differentiation of oligodendrocyte precursors into myelin- forming oligodendrocytes during development and is thought to play a crucial role in remyelination during multiple sclerosis. However, itisstillunclearhowOLIG1interactswithitstranscriptionalcofactorsandDNAtargets.OLIG1wasreportedlyrestrictedtomammals,but we demonstrate here that zebrafish and other teleosts also possess an OLIG1 homolog. In zebrafish, as in mammals, Olig1 is expressed in the oligodendrocyte lineage. Olig1 associates physically with another myelin-associated transcription factor, Sox10, and the Olig1/Sox10 complex activates mbp (myelin basic protein) transcription via conserved DNA sequence motifs in the mbp promoter region. In contrast, Olig2 does not bind to Sox10 in zebrafish, although both OLIG1 and OLIG2 bind SOX10 in mouse. Key words: Olig1; Olig2; Sox10; Mbp; oligodendrocyte; myelin; zebrafish; mouse; evolution; development Introduction directly regulates Mbp transcription (Stolt et al., 2002), and over- Myelin, the multilayered glial sheath around axons, is one of the expression of SOX10 alone is sufficient to induce myelin gene defining features of jawed vertebrates (gnathostomes). It is expression in embryonic chick spinal cord (Liu et al., 2007).
    [Show full text]