ZNF384-Related Fusion Genes Define a Subgroup of Childhood B-Cell

Total Page:16

File Type:pdf, Size:1020Kb

ZNF384-Related Fusion Genes Define a Subgroup of Childhood B-Cell Acute Lymphoblastic Leukemia SUPPLEMENTARY APPENDIX ZNF384 -related fusion genes define a subgroup of childhood B-cell precursor acute lymphoblastic leukemia with a characteristic immunotype Shinsuke Hirabayashi, 1,2 Kentaro Ohki, 1 Kazuhiko Nakabayashi, 3 Hitoshi Ichikawa, 4 Yukihide Momozawa, 5 Kohji Oka - mura, 6 Akinori Yaguchi, 1,7 Kazuki Terada, 1 Yuya Saito, 1,8 Ai Yoshimi, 1,9 Hiroko Ogata-Kawata, 3 Hiromi Sakamoto, 4 Moto - hiro Kato, 1,10 Junya Fujimura, 7 Moeko Hino, 11 Akitoshi Kinoshita, 12 Harumi Kakuda, 13 Hidemitsu Kurosawa, 14 Keisuke Kato, 9 Ryosuke Kajiwara, 15 Koichi Moriwaki, 16 Tsuyoshi Morimoto, 17 Kozue Nakamura, 18 Yasushi Noguchi, 19 Tomoo Osumi, 1,20 Kazuo Sakashita, 21 Junko Takita, 22 Yuki Yuza, 8 Koich Matsuda, 23 Teruhiko Yoshida, 4 Kenji Matsumoto, 24 Kenichiro Hata, 3 Michiaki Kubo, 5 Yoichi Matsubara, 25 Takashi Fukushima, 26 Katsuyoshi Koh, 27 Atsushi Manabe, 2 Akira Ohara 28 and Nobutaka Kiyokawa 1 for the Tokyo Children’s Cancer Study Group (TCCSG) 1Department of Pediatric Hematology and Oncology Research, National Research Institute for Child Health and Development, Seta - gaya-ku, Tokyo; 2Department of Pediatrics, St. Luke's International Hospital, Chuo-ku, Tokyo; 3Department of Maternal-Fetal Biology, National Research Institute for Child Health and Development, Setagaya-ku, Tokyo; 4Division of Genetics, National Cancer Center Re - search Institute, Chuo-ku, Tokyo; 5Laboratory for Genotyping Development, Center for Integrative Medical Sciences (IMS), RIKEN, Yoko - hama-shi, Kanagawa; 6Department of Systems BioMedicine, National Research Institute for Child Health and Development, Setagaya-ku, Tokyo; 7Department of Pediatrics and Adolescent Medicine, Juntendo University Graduate School of Medicine, Bunkyo-ku, Tokyo; 8Department of Hematology/Oncology, Tokyo Metropolitan Children’s Medical Center, Fuchu-shi, Tokyo; 9Division of Pediatric Hematology and Oncology, Ibaraki Children’s Hospital, Mito-shi, Ibaraki; 10 Division of Stem Cell Transplant and Cellular Therapy, Chil - dren’s Cancer Center, National Center for Child Health and Development, Setagaya-ku, Tokyo; 11 Department of Pediatrics, Chiba Uni - versity Graduate School of Medicine, Chiba-shi, Chiba; 12 Department of Pediatrics, St. Marianna University School of Medicine, Kawasaki-shi, Kanagawa; 13 Department of Haematology/Oncology, Chiba Children’s Hospital, Chiba-shi, Chiba; 14 Department of Pedi - atrics, Dokkyo Medical University, Mibu, Tochigi; 15 Department of Pediatrics, Yokohama City University Hospital, Yokohama-shi, Kana - gawa; 16 Department of Pediatrics, Saitama Medical Center, Saitama Medical University, Kawagoe-shi, Saitama; 17 Department of Pediatrics, Tokai University School of Medicine, Isehara-shi, Kanagawa; 18 Department of Pediatrics, Teikyo University School of Medi - cine, Itabashi-ku, Tokyo; 19 Department of Pediatrics, Japanese Red Cross Narita Hospital, Narita-shi, Chiba; 20 Division of Leukemia and Lymphoma, Children’s Cancer Center, National Center for Child Health and Development, Setagaya-ku, Tokyo; 21 Department of Hema - tology/Oncology, Nagano Children's Hospital, Azumino-shi, Nagano; 22 Department of Pediatrics, Graduate School of Medicine, Univer - sity of Tokyo, Bunkyo-ku, Tokyo; 23 Laboratory of Clinical Sequence,Department of Computational biology and medical Sciences, Graduate school of Frontier Sciences, The University of Tokyo, Minato-ku, Tokyo; 24 Department of Allergy and Clinical Immunology, Na - tional Research Institute for Child Health and Development, Setagaya-ku, Tokyo; 25 National Research Institute for Child Health and De - velopment, Setagaya-ku, Tokyo; 26 Department of Child Health, Faculty of Medicine, University of Tsukuba, Tsukuba-shi, Ibaraki; 27 Department of Hematology/Oncology, Saitama Children’s Medical Center, Saitama-shi, Saitama and 28 Department of Pediatrics, Toho University Omori Medical Center, Ohta-ku, Tokyo, Japan ©2017 Ferrata Storti Foundation. This is an open-access paper. doi:10.3324/haematol. 2016.151035 Received: June 14, 2016. Accepted: September 14, 2016. Pre-published: September 15, 2016. Correspondence: [email protected] or [email protected] Supplementary Methods Cases As a result of whole transcriptome sequencing performed after the publication of our previous report1 on the samples obtained from the cases diagnosed initially as "B-others" enrolled in the L0416/0616 study, we identified one of each minor BCR-ABL, ETV6-RUNX1, and TCF3-PBX1-positive cases and, thus, these 3 cases were removed from the B-others group and re-classified to each category of genetic alteration. As a consequence, 333 ALL patients enrolled in the TCCSG L0416/0616 study included 130 B-others patients and 16 BCR-ABL1+, 18 E2A-PBX1+, 54 ETV6-RUNX1+, 3 MLL-AF4+, 2 MLL-AF9+, 65 hyperdiploid, and 3 hypodiploid (291 BCP-ALL patients) as well as 41 T-ALL and 1 unclassified ALL (Supplementary Table 1). In addition to the 15 TCF3-ZNF384-positive patients presented in the manuscript, we identified 1 more patient harboring TCF3-ZNF384 (Case26). Case-26 was female, being 5 years old at presentation, and TCF3-ZNF384 (fusion of TCF3 Ex13 and ZNF384 Ex2) was identified in the BM specimen obtained at the 3rd relapse. As additional genetic mutation, kinase activating point mutation in FLT3 and CDKN2A/B deletion were identified by whole exome sequencing and MLPA, respectively. This patient received SCT as a salvage therapy but died. Since other clinical information of on this patient was not available, she was excluded from our study. Whole transcriptome sequencing and detection of fusion genes Total RNAs were extracted using the miRNeasy Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer's instruction. After qualification using the Agilent RNA 6000 Nano Kit (Agilent Technologies, Santa Clara, CA, USA), cDNA libraries were prepared from 1 µg of total RNA using the TruSeq RNA sample preparation kit v2 (Illumina, Inc., San Diego, CA, USA, catalog # RS-122-2001). The resultant libraries were quantified using the KAPA Library Quantification Kit (KAPA Biosystems, Inc., Woburn, MA, USA, catalog # KK4835) and checked for quality and size using the Agilent High Sensitivity DNA Kit (Agilent, catalog # 5067-4626). The samples were loaded onto cBot (Illumina) for clustering on a flow cell, and the flow cell was then sequenced using HiSeq2500 (Illumina) according to the manufacturer’s instructions. A paired-end (2 x 101) run was performed using the SBS Kit v4-HS (Illumina, catalog # FC-401-3001). More than 9 Gb of sequence was obtained for each library by paired-end sequencing (126 bp x2). To avoid multiple counting of each fusion transcript, RNA sequencing data were used after the removal of paired-end reads with the identical nucleotide sequence, which had probably been derived from PCR duplicates during library preparation. For the prediction of fusion genes, the deFuse program2 was used as described previously.3 After applying default filtering of this program, potential alternative splicing and read-through products that the program predicted were eliminated, and candidates that had exon boundary junctions were selected. Microarray and data analysis The data were normalized by the MAS5 algorithm and the baseline was transformed to the median of all samples. The data were filtered with the following steps. (1) Genes with a signal intensity lower than 1.0 in each sample were eliminated (Filtered on expression). (2) Genes that were scored as absent in all samples were eliminated (Filtered on Flags [Present and Marginal]). (3) For the multi-group comparison between ZNF384-related fusion genes positive ALL and TCF3-PBX1-positive or ZNF384 wild-type ALL, one-way analysis of variance (ANOVA) was conducted selecting probe sets with a p-value of <0.01. For the direct comparison between TCF3-ZNF384-positive and TCF3-PBX1-positive or ZNF384 wild-type ALL, a moderated t-test was conducted selecting probe sets with a p-value of <0.01. Performing fold-change analysis using filtered genes, genes were selected that exhibited 2.0-fold or more higher or lower expressions in TCF3-ZNF384-positive ALL in comparison with TCF3-PBX1-positive (Figure 4A, Supplementary Figure 2C, Supplementary Table 6) or ZNF384 wild-type ALL (Supplementary Figure 2F, Supplementary Tables 8). For direct multigroup comparison including conventional genetic subtypes and ZNF384-related fusion gene positive ALL or B-others, hierarchical clustering analysis was performed on the filtered microarray data using the selected classifying gene probe sets reported by Roberts et al. 4 and Den Boer et al.5-6 that can separate the major subtypes of pediatric ALL in distinct clusters. The functional analyses of genes characteristic of TCF3-ZNF384-positive ALLs in comparison with TCF3-PBX1-positive or ZNF384 wild-type ALLs were done using gene set enrichment analysis (GSEA) according to the developers instructions [http://software.broadinstitute.org/gsea/index.jsp] as well as previous reports.7-9 Ranked gene lists of the up- and down-regulated genes in the signatures of TCF3-ZNF384-positive ALL (10 cases) in comparison with TCF3-PBX1-positive (19 cases) or ZNF384 wild-type ALL (83 out of 104 cases that located in distinct cluster from ZNF384-related fusion genes positive cases by hierarchical clustering) were created. The analysis included gene sets for hematopoietic stage6 and human immunologic signatures (C7) from MsigDB (http://software.broadinstitute.org/gsea/index.jsp). GEO accession number, GSE79533. Exome data analyses: Mapping
Recommended publications
  • Genome-Wide Characterization of PRE-1 Reveals a Hidden Evolutionary Relationship Between Suidae and Primates
    bioRxiv preprint doi: https://doi.org/10.1101/025791; this version posted August 31, 2015. 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 4.0 International license. Genome-wide characterization of PRE-1 reveals a hidden evolutionary relationship between suidae and primates Hao Yu1,, Qingyan Wu1, Jing Zhag1, Ying zhang1, Chao Lu1, Yunyun Cheng1, Zhihui Zhao1, Andreas Windemuth3, Di Liu2,, Linlin Hao1 1 College of Animal Science, Jilin University, Changchun 130062, China 2 Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China 3 Abcam, Firefly BioWorks Inc, United States Corresponding author: Linlin Hao ([email protected]); Di liu ([email protected]); Andreas Windemuth ([email protected]) Abstract We identified and characterized a free PRE-1 element inserted into the promoter region of the porcine IGFBP7 gene whose integration mechanisms into the genome, including copy number, distribution preferences, capacity to exonize and phyloclustering pattern are similar to that of the primate Alu element. 98% of these PRE-1 elements also contain two conserved internal AluI restriction enzyme recognition sites, and the RNA structure of PRE1 can be folded into a two arms model like the Alu RNA structure. It is more surprising that the length of the PRE-1 fragments is nearly the same in 20 chromosomes and positively correlated to its fracture site frequency. All of these fracture sites are close to the mutation hot spots of PRE-1 families, and most of these hot spots are located in the non-complementary fragile regions of the PRE-1 RNA structure.
    [Show full text]
  • Autophagy: from Basic Science to Clinical Application
    nature publishing group REVIEW See COMMENTARY page XX Autophagy: from basic science to clinical application J Va n L i m b e r g e n 1 , 2 , 3 , C S t e v e n s 4 , E R N i m m o 1 , D C W i l s o n 2 , 3 a n d J S a t s a n g i 1 Autophagy is a cellular pathway involved in protein and organelle degradation, which is likely to represent an innate adaptation to starvation. In times of nutrient deficiency, the cell can self-digest and recycle some nonessential components through nonselective autophagy, thus sustaining minimal growth requirements until a food source becomes available. Over recent years, autophagy has been implicated in an increasing number of clinical scenarios, notably infectious diseases, cancer, neurodegenerative diseases, and autoimmunity. The recent identification of the importance of autophagy genes in the genetic susceptibility to Crohn ’ s disease suggests that a selective autophagic response may play a crucial role in the pathogenesis of common complex immune-mediated diseases. In this review, we discuss the autophagic mechanisms, their molecular regulation, and summarize their clinical relevance. This progress has led to great interest in the therapeutic potential of manipulation of both selective and nonselective autophagy in established disease. INTRODUCTION The ability to adapt to environmental change is essential for sur- Autophagy encompasses several distinct processes involving vival. This is true for the organism as a whole and for individual the delivery of portions of the cytoplasm to the lysosome for cells alike.
    [Show full text]
  • F2RL2 Antibody Cat
    F2RL2 Antibody Cat. No.: 56-323 F2RL2 Antibody F2RL2 Antibody immunohistochemistry analysis in formalin fixed and paraffin embedded human heart tissue followed by peroxidase conjugation of the secondary antibody and DAB staining. Specifications HOST SPECIES: Rabbit SPECIES REACTIVITY: Human This F2RL2 antibody is generated from rabbits immunized with a KLH conjugated IMMUNOGEN: synthetic peptide between 21-50 amino acids from the N-terminal region of human F2RL2. TESTED APPLICATIONS: IHC-P, WB For WB starting dilution is: 1:1000 APPLICATIONS: For IHC-P starting dilution is: 1:10~50 PREDICTED MOLECULAR 43 kDa WEIGHT: September 25, 2021 1 https://www.prosci-inc.com/f2rl2-antibody-56-323.html Properties This antibody is purified through a protein A column, followed by peptide affinity PURIFICATION: purification. CLONALITY: Polyclonal ISOTYPE: Rabbit Ig CONJUGATE: Unconjugated PHYSICAL STATE: Liquid BUFFER: Supplied in PBS with 0.09% (W/V) sodium azide. CONCENTRATION: batch dependent Store at 4˚C for three months and -20˚C, stable for up to one year. As with all antibodies STORAGE CONDITIONS: care should be taken to avoid repeated freeze thaw cycles. Antibodies should not be exposed to prolonged high temperatures. Additional Info OFFICIAL SYMBOL: F2RL2 Proteinase-activated receptor 3, PAR-3, Coagulation factor II receptor-like 2, Thrombin ALTERNATE NAMES: receptor-like 2, F2RL2, PAR3 ACCESSION NO.: O00254 GENE ID: 2151 USER NOTE: Optimal dilutions for each application to be determined by the researcher. Background and References Coagulation factor II (thrombin) receptor-like 2 (F2RL2) is a member of the large family of 7-transmembrane-region receptors that couple to guanosine-nucleotide-binding proteins.
    [Show full text]
  • The Title of the Article
    Mechanism-Anchored Profiling Derived from Epigenetic Networks Predicts Outcome in Acute Lymphoblastic Leukemia Xinan Yang, PhD1, Yong Huang, MD1, James L Chen, MD1, Jianming Xie, MSc2, Xiao Sun, PhD2, Yves A Lussier, MD1,3,4§ 1Center for Biomedical Informatics and Section of Genetic Medicine, Department of Medicine, The University of Chicago, Chicago, IL 60637 USA 2State Key Laboratory of Bioelectronics, Southeast University, 210096 Nanjing, P.R.China 3The University of Chicago Cancer Research Center, and The Ludwig Center for Metastasis Research, The University of Chicago, Chicago, IL 60637 USA 4The Institute for Genomics and Systems Biology, and the Computational Institute, The University of Chicago, Chicago, IL 60637 USA §Corresponding author Email addresses: XY: [email protected] YH: [email protected] JC: [email protected] JX: [email protected] XS: [email protected] YL: [email protected] - 1 - Abstract Background Current outcome predictors based on “molecular profiling” rely on gene lists selected without consideration for their molecular mechanisms. This study was designed to demonstrate that we could learn about genes related to a specific mechanism and further use this knowledge to predict outcome in patients – a paradigm shift towards accurate “mechanism-anchored profiling”. We propose a novel algorithm, PGnet, which predicts a tripartite mechanism-anchored network associated to epigenetic regulation consisting of phenotypes, genes and mechanisms. Genes termed as GEMs in this network meet all of the following criteria: (i) they are co-expressed with genes known to be involved in the biological mechanism of interest, (ii) they are also differentially expressed between distinct phenotypes relevant to the study, and (iii) as a biomodule, genes correlate with both the mechanism and the phenotype.
    [Show full text]
  • TRANSCRIPTIONAL REGULATION of Hur in RENAL STRESS
    TRANSCRIPTIONAL REGULATION OF HuR IN RENAL STRESS DISSERTATION Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the Graduate School of The Ohio State University By Sudha Suman Govindaraju Graduate Program in Biochemistry The Ohio State University 2014 Dissertation Committee: Dr. Beth S. Lee, Ph.D., Advisor Dr. Kathleen Boris-Lawrie, Ph.D. Dr. Sissy M. Jhiang, Ph.D. Dr. Arthur R. Strauch, Ph.D Abstract HuR is a ubiquitously expressed RNA-binding protein that affects the post- transcriptional life of thousands of cellular mRNAs by regulating transcript stability and translation. HuR can post-transcriptionally regulate gene expression and modulate cellular responses to stress, differentiation, proliferation, apoptosis, senescence, inflammation, and the immune response. It is an important mediator of survival during cellular stress, but when inappropriately expressed, can promote oncogenic transformation. Not surprisingly, the expression of HuR itself is tightly regulated at multiple transcriptional and post-transcriptional levels. Previous studies demonstrated the existence of two alternate HuR transcripts that differ in their 5’ untranslated regions and have markedly different translatabilities. These forms were also found to be reciprocally expressed following cellular stress in kidney proximal tubule cell lines, and the shorter, more readily translatable variant was shown to be regulated by Smad 1/5/8 pathway and bone morphogenetic protein-7 (BMP-7) signaling. In this study, the factors that promote transcription of the longer alternate form were identified. NF-κB was shown to be important for expression of the long HuR mRNA, as was a newly identified region with potential for binding the Sp/KLF families of transcription factors.
    [Show full text]
  • Expression Gene Network Analyses Reveal Molecular Mechanisms And
    www.nature.com/scientificreports OPEN Diferential expression and co- expression gene network analyses reveal molecular mechanisms and candidate biomarkers involved in breast muscle myopathies in chicken Eva Pampouille1,2, Christelle Hennequet-Antier1, Christophe Praud1, Amélie Juanchich1, Aurélien Brionne1, Estelle Godet1, Thierry Bordeau1, Fréderic Fagnoul2, Elisabeth Le Bihan-Duval1 & Cécile Berri1* The broiler industry is facing an increasing prevalence of breast myopathies, such as white striping (WS) and wooden breast (WB), and the precise aetiology of these occurrences remains poorly understood. To progress our understanding of the structural changes and molecular pathways involved in these myopathies, a transcriptomic analysis was performed using an 8 × 60 K Agilent chicken microarray and histological study. The study used pectoralis major muscles from three groups: slow-growing animals (n = 8), fast-growing animals visually free from defects (n = 8), or severely afected by both WS and WB (n = 8). In addition, a weighted correlation network analysis was performed to investigate the relationship between modules of co-expressed genes and histological traits. Functional analysis suggested that selection for fast growing and breast meat yield has progressively led to conditions favouring metabolic shifts towards alternative catabolic pathways to produce energy, leading to an adaptive response to oxidative stress and the frst signs of infammatory, regeneration and fbrosis processes. All these processes are intensifed in muscles afected by severe myopathies, in which new mechanisms related to cellular defences and remodelling seem also activated. Furthermore, our study opens new perspectives for myopathy diagnosis by highlighting fne histological phenotypes and genes whose expression was strongly correlated with defects. Te poultry industry relies on the production of fast-growing chickens, which are slaughtered at high weights and intended for cutting and processing.
    [Show full text]
  • Synergistic Genetic Interactions Between Pkhd1 and Pkd1 Result in an ARPKD-Like Phenotype in Murine Models
    BASIC RESEARCH www.jasn.org Synergistic Genetic Interactions between Pkhd1 and Pkd1 Result in an ARPKD-Like Phenotype in Murine Models Rory J. Olson,1 Katharina Hopp ,2 Harrison Wells,3 Jessica M. Smith,3 Jessica Furtado,1,4 Megan M. Constans,3 Diana L. Escobar,3 Aron M. Geurts,5 Vicente E. Torres,3 and Peter C. Harris 1,3 Due to the number of contributing authors, the affiliations are listed at the end of this article. ABSTRACT Background Autosomal recessive polycystic kidney disease (ARPKD) and autosomal dominant polycystic kidney disease (ADPKD) are genetically distinct, with ADPKD usually caused by the genes PKD1 or PKD2 (encoding polycystin-1 and polycystin-2, respectively) and ARPKD caused by PKHD1 (encoding fibrocys- tin/polyductin [FPC]). Primary cilia have been considered central to PKD pathogenesis due to protein localization and common cystic phenotypes in syndromic ciliopathies, but their relevance is questioned in the simple PKDs. ARPKD’s mild phenotype in murine models versus in humans has hampered investi- gating its pathogenesis. Methods To study the interaction between Pkhd1 and Pkd1, including dosage effects on the phenotype, we generated digenic mouse and rat models and characterized and compared digenic, monogenic, and wild-type phenotypes. Results The genetic interaction was synergistic in both species, with digenic animals exhibiting pheno- types of rapidly progressive PKD and early lethality resembling classic ARPKD. Genetic interaction be- tween Pkhd1 and Pkd1 depended on dosage in the digenic murine models, with no significant enhancement of the monogenic phenotype until a threshold of reduced expression at the second locus was breached.
    [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]
  • Transcriptomic Analysis of the Aquaporin (AQP) Gene Family
    Pancreatology 19 (2019) 436e442 Contents lists available at ScienceDirect Pancreatology journal homepage: www.elsevier.com/locate/pan Transcriptomic analysis of the Aquaporin (AQP) gene family interactome identifies a molecular panel of four prognostic markers in patients with pancreatic ductal adenocarcinoma Dimitrios E. Magouliotis a, b, Vasiliki S. Tasiopoulou c, Konstantinos Dimas d, * Nikos Sakellaridis d, Konstantina A. Svokos e, Alexis A. Svokos f, Dimitris Zacharoulis b, a Division of Surgery and Interventional Science, Faculty of Medical Sciences, UCL, London, UK b Department of Surgery, University of Thessaly, Biopolis, Larissa, Greece c Faculty of Medicine, School of Health Sciences, University of Thessaly, Biopolis, Larissa, Greece d Department of Pharmacology, Faculty of Medicine, School of Health Sciences, University of Thessaly, Biopolis, Larissa, Greece e The Warren Alpert Medical School of Brown University, Providence, RI, USA f Riverside Regional Medical Center, Newport News, VA, USA article info abstract Article history: Background: This study aimed to assess the differential gene expression of aquaporin (AQP) gene family Received 14 October 2018 interactome in pancreatic ductal adenocarcinoma (PDAC) using data mining techniques to identify novel Received in revised form candidate genes intervening in the pathogenicity of PDAC. 29 January 2019 Method: Transcriptome data mining techniques were used in order to construct the interactome of the Accepted 9 February 2019 AQP gene family and to determine which genes members are differentially expressed in PDAC as Available online 11 February 2019 compared to controls. The same techniques were used in order to evaluate the potential prognostic role of the differentially expressed genes. Keywords: PDAC Results: Transcriptome microarray data of four GEO datasets were incorporated, including 142 primary Aquaporin tumor samples and 104 normal pancreatic tissue samples.
    [Show full text]
  • Deficiency in Class III PI3-Kinase Confers Postnatal Lethality with IBD
    ARTICLE DOI: 10.1038/s41467-018-05105-8 OPEN Deficiency in class III PI3-kinase confers postnatal lethality with IBD-like features in zebrafish Shaoyang Zhao1,2,3, Jianhong Xia2,3, Xiuhua Wu2,3, Leilei Zhang2,3, Pengtao Wang2,3, Haiyun Wang2,3, Heying Li2, Xiaoshan Wang 2, Yan Chen 2, Jean Agnetti2, Yinxiong Li 2, Duanqing Pei2,3 & Xiaodong Shu2,3 The class III PI3-kinase (PIK3C3) is an enzyme responsible for the generation of phospha- tidylinositol 3-phosphate (PI3P), a critical component of vesicular membrane. Here, we report 1234567890():,; that PIK3C3 deficiency in zebrafish results in intestinal injury and inflammation. In pik3c3 mutants, gut tube forms but fails to be maintained. Gene expression analysis reveals that barrier-function-related inflammatory bowel disease (IBD) susceptibility genes (e-cadherin, hnf4a, ttc7a) are suppressed, while inflammatory response genes are stimulated in the mutants. Histological analysis shows neutrophil infiltration into mutant intestinal epithelium and the clearance of gut microbiota. Yet, gut microorganisms appear dispensable as mutants cultured under germ-free condition have similar intestinal defects. Mechanistically, we show that PIK3C3 deficiency suppresses the formation of PI3P and disrupts the polarized dis- tribution of cell-junction proteins in intestinal epithelial cells. These results not only reveal a role of PIK3C3 in gut homeostasis, but also provide a zebrafish IBD model. 1 School of Life Sciences, University of Science and Technology of China, 230027 Hefei, Anhui, China. 2 CAS Key Laboratory of Regenerative Biology, Guangzhou Institute of Biomedicine and Health-Guangzhou Medical University Joint School of Biological Sciences, South China Institute for Stem Cell Biology and Regenerative Medicine, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, 510530 Guangzhou, China.
    [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]