1 Supplementary Materials and Methods Gene Set Enrichment Analysis for Gene Expression Data We Compared the Gene Expression Leve

Total Page:16

File Type:pdf, Size:1020Kb

1 Supplementary Materials and Methods Gene Set Enrichment Analysis for Gene Expression Data We Compared the Gene Expression Leve Supplementary Materials and Methods Gene set enrichment analysis for gene expression data We compared the gene expression levels from 2 different phenotypes (i.e., drug sensitive versus resistant) and picked up the genes which had significant different expression for Gene set enrichment analysis (GSEA) by using Molecular Signatures Database(V3.0). Gene set enrichment analysis was carried out by computing overlaps with canonical pathways(CP) and gene ontology (GO) gene sets(C5), obtained from the Broad Institute [1]. Genes in Gene Set (K), Genes in Overlap (k), k/K and P value were used to rank the pathways enriched in each phenotype. We used with 10363 genes in 2334 pathways as the gene set in this study. Supplementary References 1.Subramanian A, Tamayo P, Mootha VK, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 2005;102:15545– 50. 1 Supplementary Results Other drug response-genotype associations in this study Cells with wild genotype were more sensitive to irinotecan than those with TGFβR2 mutation and deletion together (P=0.03) (Supplementary Table S6). PC cells with wild genotype in ADAM gene family were more sensitive to gemcitabine and docetaxel than those with mutation in this family (p<0.05) (Supplementary Table S6). PC cells with wild genotype were more sensitive to docetaxel than those with DEPDC or TTN mutation (p<0.05) (Supplementary Table S6). PC cells with BAI3 mutation were more sensitive to cisplatin than those with wild genotype (p=0.02) (Supplementary Table S6). We also examined some important members of TGFβ pathway (SMAD3, SMAD4, TGFβR1 and TGFβR2) individually and together for additional alterations. Cells were grouped by whether they had a mutation or deletion in any of these four genes. PC cells with the deletion of SMAD Pathway genes (SMAD3, SMAD4, TGFβR2 or TGFβR3) were less sensitive to triptolide (P=0.03) (Supplementary Table S6). Logistic regression models were set up to analyze the factors related to drug response. Different gene statuses were considered as sensitive or resistant genetic factors for each drug. We found AHNAK or DEPDC mutation were resistant factors for cisplatin, the mutation of DEPDC, TTN, FMN2, BAI3 or MLL were resistant factors for docetaxel, ODZ4, OVCH1 or SMADPath gene inactivation were resistant factors for triptolide, and BAI3 inactivation was a resistant factor for gemcitabine (Supplementary Table S7). CDK.SMAD.P53 inactivation was a sensitive factor for artemisinin (Supplementary Table S7). We compared the IC50 values of each drug between familial and sporadic pancreatic cancer cell lines and didn’t find differences of drug response for anticancer drugs (p>0.05) (Supplementary Table S17). Different gene expression profiles were correlated with drug responses of broad classes of anticancer agents in human PC cells We did the Spearman correlation analysis between IC50 values of each drug and the gene expression levels for 19 cell lines in the discovery screen. According to this rule, we discovered a series of related genes, as shown in Supplementary Table S10. We found that expression levels of some genes were closely correlated to anticancer drug responses in human PC cell lines. NIT2, AKT1S1, STK17A, MAP3K12, MAP4K2 and MAPK11 were closely related to drug response of cisplatin. MYH9, ABL2, GPR125, RCC2, STK17B and TTK were related to chemosensitivity of MMC. RAB6B, RHOC, SNIP1, TUSC1, MED23, MYB and WNT11 were related to drug response of triptolide. FBXL12, LDB1, SIN3B, TSC22D3, ABCG2, CDGAP, 2 MAP3K10, MED16, NOTCH3 and STAT6 were related to drug response of gemcitabine. CDC37, PIM3, TSPAN5 and USP5 were correlated to drug response of Parp1inhibitor. We also compared gene expression levels of 20,661 protein-coding genes from sensitive and resistant cell lines of each drug. The rules to define sensitive versus resistant by IC50 value were as follows: <5 nM versus >20 nM for gemcitabine, <1 nM versus >5 nM for docetaxel and triptolide, <10 nM versus >50 nM for MMC, <200 nM versus >10000 nM for irinotecan, cisplatin and Parp1 inhibitor and <300 nM versus >10000 nM for artemisinin. We selected the genes which had significant different expression levels for Gene set enrichment analysis and results were shown in Supplementary Table S18. We found that genes with significant different gene expression gathered in cell cycle, cell cycle checkpoint, DNA replication, DNA repair, apoptosis and MAPK Pathway for gemcitabine; genes for MMC gathered in DNA damage,P38 signaling, protein deacetylase activity and apoptosis pathway; genes for triptolide were involved in MTA3 and JUN Pathway; genes for Parp1 inhibitor gathered in mitochondrial metabolism; genes for docetaxel were involved in serine/threonine/tyrosine kinase activity, chromatin assembly and disaseembly, p38 and JUN Signaling and nucleosome assembly; genes for artemisinin were involved in PI3K Pathway. We listed discovered genes of cell cycle, apoptosis, autophagy, DNA repair, tumor suppressor gene and anticancer drug metabolism pathway between sensitive and resistant human pancreatic cancer cells in Supplementary Figure S4 and Supplementary Table S16. Genes were selected for sensitive versus resistant based on P<0.05 of the median IC50s. Cell cycle genes (CDC2L1,CDC37L1,CDC73,E2F1), apaoptosis and autophagy genes(GADD45B,ATG10, ATG4B), DNA repair genes(ERCC3, PARP2,PLK2, RAD17,RAD54L2,TP53BP2,TSC2) and some MAPK pathway genes were upregulated in triptolide sensitive PC cell lines (Supplementary Figure S4A and Supplementary Table S16). It was also observed that cell cycle genes (CDC20,CDC25A,CDC27,CDC45L,CDC6, CDK5R1,CDKN2B), apaoptosis and autophagy genes (ATG10,ATG4C,CASP5), and DNA repair genes (AURKA, AURKB, BLM,BRCA1MSH2,RAD21,RAD23A, RAD51, RAD51AP1,SMC2,XRCC5,XRCC6) were upregulated in gemcitabine sensitive PC cell lines (Supplementary Figure S4C,S4D and Supplementary Table S16). Apaoptosis and autophagy genes (ATG9A,BCL3) were upregulated, but cell cycle genes (CDC2,CDC40) and DNA repair genes (MDM2, MLH1,MSH3,BRCA1) were downregulated in MMC sensitive pancreatic cell lines (Supplementary Figure S4E and Supplementary Table S16). Different gene expression profiles were correlated with DPC4/SMAD4, TP53 or P16/CDKN2A inactivations 3 In order to assess the gene expression profile under DPC4 deficiency, we compared gene expression levels of 20,661 protein-coding genes from 2 different genotypes (DPC4.md versus DPC4.wt) and then discovered the genes in cell cycle, apoptosis and DNA repair pathway which had statistically significantly different expression levels. Compared to cells with DPC4.wt, cells with DPC4.md showed these genes were mainly involved in pathways related to cell cycle checkpoint (CDKN1B, CDC2L6, CDCA4), DNA repair (FANCD2, TP73, XRCC3), apoptosis (BCL2L1,BCL2L13,BCL7B) and MAPK signaling (Supplementary Figure S5A and Supplementary Table S19). Overall, these results suggest that under DPC4.md genetic background, cellular decreased cell cycle checkpoint arrest, apoptosis and DNA repair pathway change may confer sensitivity to cisplatin and irinotecan and nonresponse to gemcitabine. In order to know the gene expression profile under TP53 deficiency, we compared gene expression levels from 2 different genotypes (TP53.md versus TP53.wt) and then discovered the genes with P value less than 0.05. Compared to cells with TP53.wt, cells with TP53.md showed upregulated genes gathered in pathways related to cell cycle checkpoint (CDKN1C,CDKN2A,CDKN2B, CDKN2C,CDC25B), DNA repair (RAD51,RAD54L,BLM,MLF1), apoptosis (CASP7,CASP8AP2,CASP9,BCL2L11,CYCS), autophagy (ATG12, ATG4B,ATG4D,ATG7) and MAPK signaling (Supplementary Figure S5B,S5C and Supplementary Table S20). These results suggest that increased cell cycle G1/S checkpoint arrest, increased exacerbated spontaneous HR observed in p53 defective cells, upregulation of apoptosis and autophagy genes may contribute to the sensitivity of triptolide and Parp1 inhibitor. We compared gene expression levels from 2 different genotypes (P16.md deletion versus P16.wt) and then found the genes with P value less than 0.05. Compared to cells with P16 wild type, cells with P16 deletion showed changed genes gathered in pathways related to cell cycle checkpoint (AKT1,CHEK2,CDC25C,CCNB1IP1,CDC14A, CDCA2,CDK10,CDK2AP1), apoptosis (FADD,FAIM,FAS,BAG5, BMF),DNA repair(PLK1,RAD21,RAD9B,WRN,XRCC2, PCNA) and autophagy (ATG9A) pathway (Supplementary Figure S5D and Supplementary Table S21). Overall, these results suggest that underregulation of DNA repair genes and cell cycle checkpoint change may confer nonresponse to gemcitabine and MMC under P16.md. 4 Supplementary tables: Table S1 Characteristics of human pancreatic cancer cell lines used in this study Table S2 Characteristics of human DPC4/SMAD4 isogenic pancreatic cancer cell lines Table S3 Characteristics of human TP53 isogenic colon cancer cell lines Table S4 IC50 values (nM) of Gemcitabine, Triptolide, Docetaxel, MMC, Cisplatin, Irinotecan, Parp1 inhibitor and Artemisinin for human PC cell lines Table S5 IC50 values (nM) of Gemcitabine, Triptolide, Docetaxel, MMC, Cisplatin, Irinotecan, Parp1 inhibitor and Artemisinin for human pancreatic cancer cell lines Table S6 Correlation analysis between chemosensitivity and somatic mutations alone, homozygous deletion alone or both in pancreatic cancer cells identified in the discovery and prevalence screen arranged by drugs Table S7 Logistic regression analysis for drug response to genetic status Table
Recommended publications
  • Uwaterloo Latex Thesis Template
    New Algorithms for Predicting Conformational Polymorphism and Inferring Direct Couplings for Side Chains of Proteins by Laleh Soltan Ghoraie A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Doctor of Philosophy in Computer Science Waterloo, Ontario, Canada, 2015 c Laleh Soltan Ghoraie 2015 I hereby declare that I am the sole author of this thesis. This is a true copy of the thesis, including any required final revisions, as accepted by my examiners. I understand that my thesis may be made electronically available to the public. ii Abstract Protein crystals populate diverse conformational ensembles. Despite much evidence that there is widespread conformational polymorphism in protein side chains, most of the x- ray crystallography data are modelled by single conformations in the Protein Data Bank. The ability to extract or to predict these conformational polymorphisms is of crucial im- portance, as it facilitates deeper understanding of protein dynamics and functionality. This dissertation describes a computational strategy capable of predicting side-chain poly- morphisms. The applied approach extends a particular class of algorithms for side-chain prediction by modelling the side-chain dihedral angles more appropriately as continuous rather than discrete variables. Employing a new inferential technique known as particle belief propagation (PBP), we predict residue-specific distributions that encode informa- tion about side-chain polymorphisms. The predicted polymorphisms are in relatively close agreement with results from a state-of-the-art approach based on x-ray crystallography data. This approach characterizes the conformational polymorphisms of side chains us- ing electron density information, and has successfully discovered previously unmodelled conformations.
    [Show full text]
  • High-Throughput Discovery of Novel Developmental Phenotypes
    High-throughput discovery of novel developmental phenotypes The Harvard community has made this article openly available. Please share how this access benefits you. Your story matters Citation Dickinson, M. E., A. M. Flenniken, X. Ji, L. Teboul, M. D. Wong, J. K. White, T. F. Meehan, et al. 2016. “High-throughput discovery of novel developmental phenotypes.” Nature 537 (7621): 508-514. doi:10.1038/nature19356. http://dx.doi.org/10.1038/nature19356. Published Version doi:10.1038/nature19356 Citable link http://nrs.harvard.edu/urn-3:HUL.InstRepos:32071918 Terms of Use This article was downloaded from Harvard University’s DASH repository, and is made available under the terms and conditions applicable to Other Posted Material, as set forth at http:// nrs.harvard.edu/urn-3:HUL.InstRepos:dash.current.terms-of- use#LAA HHS Public Access Author manuscript Author ManuscriptAuthor Manuscript Author Nature. Manuscript Author Author manuscript; Manuscript Author available in PMC 2017 March 14. Published in final edited form as: Nature. 2016 September 22; 537(7621): 508–514. doi:10.1038/nature19356. High-throughput discovery of novel developmental phenotypes A full list of authors and affiliations appears at the end of the article. Abstract Approximately one third of all mammalian genes are essential for life. Phenotypes resulting from mouse knockouts of these genes have provided tremendous insight into gene function and congenital disorders. As part of the International Mouse Phenotyping Consortium effort to generate and phenotypically characterize 5000 knockout mouse lines, we have identified 410 Users may view, print, copy, and download text and data-mine the content in such documents, for the purposes of academic research, subject always to the full Conditions of use:http://www.nature.com/authors/editorial_policies/license.html#terms #Corresponding author: [email protected].
    [Show full text]
  • RAB6B Polyclonal Antibody
    For Research Use Only RAB6B Polyclonal antibody Catalog Number:10340-1-AP 2 Publications www.ptgcn.com Catalog Number: GenBank Accession Number: Recommended Dilutions: Basic Information 10340-1-AP BC002510 WB 1:500-1:2000 Size: GeneID (NCBI): IP 0.5-4.0 ug for IP and 1:500-1:2000 450 μg/ml 51560 for WB IHC 1:20-1:200 Source: Full Name: IF 1:10-1:100 Rabbit RAB6B, member RAS oncogene family Isotype: Calculated MW: IgG 23 kDa Purification Method: Observed MW: Antigen affinity purification 24 kDa Immunogen Catalog Number: AG0322 Applications Tested Applications: Positive Controls: IF, IHC, IP, WB,ELISA WB : mouse brain tissue; C6 cells, rat brain tissue Cited Applications: IP : mouse brain tissue; IF, WB IHC : human gliomas tissue; Species Specificity: human, mouse, rat IF : C6 cells; Cited Species: mouse, rat Note-IHC: suggested antigen retrieval with TE buffer pH 9.0; (*) Alternatively, antigen retrieval may be performed with citrate buffer pH 6.0 The human RAB genes share structural and biochemical properties with the Ras gene superfamily. Accumulating Background Information data suggests an important role for RAB proteins either in endocytosis or in biosynthetic protein transport. The transport of newly synthesized proteins from endoplasmic reticulum to the Golgi complex and to secretory vesicles involves the movement of carrier vesicles, a process that appears to involve RAB protein function. Rab6A has been shown to be a regulator of membrane traffic from the Golgi apparatus towards the endoplasmic reticulum (ER). Rab6B is encoded by an independent gene which is located on chromosome 3 region q21-q23.
    [Show full text]
  • Prognostic Significance of Autophagy-Relevant Gene Markers in Colorectal Cancer
    ORIGINAL RESEARCH published: 15 April 2021 doi: 10.3389/fonc.2021.566539 Prognostic Significance of Autophagy-Relevant Gene Markers in Colorectal Cancer Qinglian He 1, Ziqi Li 1, Jinbao Yin 1, Yuling Li 2, Yuting Yin 1, Xue Lei 1 and Wei Zhu 1* 1 Department of Pathology, Guangdong Medical University, Dongguan, China, 2 Department of Pathology, Dongguan People’s Hospital, Southern Medical University, Dongguan, China Background: Colorectal cancer (CRC) is a common malignant solid tumor with an extremely low survival rate after relapse. Previous investigations have shown that autophagy possesses a crucial function in tumors. However, there is no consensus on the value of autophagy-associated genes in predicting the prognosis of CRC patients. Edited by: This work screens autophagy-related markers and signaling pathways that may Fenglin Liu, Fudan University, China participate in the development of CRC, and establishes a prognostic model of CRC Reviewed by: based on autophagy-associated genes. Brian M. Olson, Emory University, United States Methods: Gene transcripts from the TCGA database and autophagy-associated gene Zhengzhi Zou, data from the GeneCards database were used to obtain expression levels of autophagy- South China Normal University, China associated genes, followed by Wilcox tests to screen for autophagy-related differentially Faqing Tian, Longgang District People's expressed genes. Then, 11 key autophagy-associated genes were identified through Hospital of Shenzhen, China univariate and multivariate Cox proportional hazard regression analysis and used to Yibing Chen, Zhengzhou University, China establish prognostic models. Additionally, immunohistochemical and CRC cell line data Jian Tu, were used to evaluate the results of our three autophagy-associated genes EPHB2, University of South China, China NOL3, and SNAI1 in TCGA.
    [Show full text]
  • Inhibition of Mitochondrial Complex II in Neuronal Cells Triggers Unique
    www.nature.com/scientificreports OPEN Inhibition of mitochondrial complex II in neuronal cells triggers unique pathways culminating in autophagy with implications for neurodegeneration Sathyanarayanan Ranganayaki1, Neema Jamshidi2, Mohamad Aiyaz3, Santhosh‑Kumar Rashmi4, Narayanappa Gayathri4, Pulleri Kandi Harsha5, Balasundaram Padmanabhan6 & Muchukunte Mukunda Srinivas Bharath7* Mitochondrial dysfunction and neurodegeneration underlie movement disorders such as Parkinson’s disease, Huntington’s disease and Manganism among others. As a corollary, inhibition of mitochondrial complex I (CI) and complex II (CII) by toxins 1‑methyl‑4‑phenylpyridinium (MPP+) and 3‑nitropropionic acid (3‑NPA) respectively, induced degenerative changes noted in such neurodegenerative diseases. We aimed to unravel the down‑stream pathways associated with CII inhibition and compared with CI inhibition and the Manganese (Mn) neurotoxicity. Genome‑wide transcriptomics of N27 neuronal cells exposed to 3‑NPA, compared with MPP+ and Mn revealed varied transcriptomic profle. Along with mitochondrial and synaptic pathways, Autophagy was the predominant pathway diferentially regulated in the 3‑NPA model with implications for neuronal survival. This pathway was unique to 3‑NPA, as substantiated by in silico modelling of the three toxins. Morphological and biochemical validation of autophagy markers in the cell model of 3‑NPA revealed incomplete autophagy mediated by mechanistic Target of Rapamycin Complex 2 (mTORC2) pathway. Interestingly, Brain Derived Neurotrophic Factor
    [Show full text]
  • Viewed and Published Immediately Upon Acceptance Cited in Pubmed and Archived on Pubmed Central Yours — You Keep the Copyright
    BMC Genomics BioMed Central Research article Open Access Differential gene expression in ADAM10 and mutant ADAM10 transgenic mice Claudia Prinzen1, Dietrich Trümbach2, Wolfgang Wurst2, Kristina Endres1, Rolf Postina1 and Falk Fahrenholz*1 Address: 1Johannes Gutenberg-University, Institute of Biochemistry, Mainz, Johann-Joachim-Becherweg 30, 55128 Mainz, Germany and 2Helmholtz Zentrum München – German Research Center for Environmental Health, Institute for Developmental Genetics, Ingolstädter Landstraße 1, 85764 Neuherberg, Germany Email: Claudia Prinzen - [email protected]; Dietrich Trümbach - [email protected]; Wolfgang Wurst - [email protected]; Kristina Endres - [email protected]; Rolf Postina - [email protected]; Falk Fahrenholz* - [email protected] * Corresponding author Published: 5 February 2009 Received: 19 June 2008 Accepted: 5 February 2009 BMC Genomics 2009, 10:66 doi:10.1186/1471-2164-10-66 This article is available from: http://www.biomedcentral.com/1471-2164/10/66 © 2009 Prinzen et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract Background: In a transgenic mouse model of Alzheimer disease (AD), cleavage of the amyloid precursor protein (APP) by the α-secretase ADAM10 prevented amyloid plaque formation, and alleviated cognitive deficits. Furthermore, ADAM10 overexpression increased the cortical synaptogenesis. These results suggest that upregulation of ADAM10 in the brain has beneficial effects on AD pathology. Results: To assess the influence of ADAM10 on the gene expression profile in the brain, we performed a microarray analysis using RNA isolated from brains of five months old mice overexpressing either the α-secretase ADAM10, or a dominant-negative mutant (dn) of this enzyme.
    [Show full text]
  • Download Tool
    by Submitted in partial satisfaction of the requirements for degree of in in the GRADUATE DIVISION of the UNIVERSITY OF CALIFORNIA, SAN FRANCISCO Approved: ______________________________________________________________________________ Chair ______________________________________________________________________________ ______________________________________________________________________________ ______________________________________________________________________________ ______________________________________________________________________________ Committee Members Copyright 2019 by Adolfo Cuesta ii Acknowledgements For me, completing a doctoral dissertation was a huge undertaking that was only possible with the support of many people along the way. First, I would like to thank my PhD advisor, Jack Taunton. He always gave me the space to pursue my own ideas and interests, while providing thoughtful guidance. Nearly every aspect of this project required a technique that was completely new to me. He trusted that I was up to the challenge, supported me throughout, helped me find outside resources when necessary. I remain impressed with his voracious appetite for the literature, and ability to recall some of the most subtle, yet most important details in a paper. Most of all, I am thankful that Jack has always been so generous with his time, both in person, and remotely. I’ve enjoyed our many conversations and hope that they will continue. I’d also like to thank my thesis committee, Kevan Shokat and David Agard for their valuable support, insight, and encouragement throughout this project. My lab mates in the Taunton lab made this such a pleasant experience, even on the days when things weren’t working well. I worked very closely with Tangpo Yang on the mass spectrometry aspects of this project. Xiaobo Wan taught me almost everything I know about protein crystallography. Thank you as well to Geoff Smith, Jordan Carelli, Pat Sharp, Yazmin Carassco, Keely Oltion, Nicole Wenzell, Haoyuan Wang, Steve Sethofer, and Shyam Krishnan, Shawn Ouyang and Qian Zhao.
    [Show full text]
  • Presynaptic Mechanisms for Cargo Recruitment
    Presynaptic Mechanisms for Cargo Recruitment The Harvard community has made this article openly available. Please share how this access benefits you. Your story matters Citation Nyitrai, Hajnalka. 2019. Presynaptic Mechanisms for Cargo Recruitment. Doctoral dissertation, Harvard University, Graduate School of Arts & Sciences. Citable link http://nrs.harvard.edu/urn-3:HUL.InstRepos:42013129 Terms of Use This article was downloaded from Harvard University’s DASH repository, and is made available under the terms and conditions applicable to Other Posted Material, as set forth at http:// nrs.harvard.edu/urn-3:HUL.InstRepos:dash.current.terms-of- use#LAA Presynaptic Mechanisms for Cargo Recruitment A dissertation presented by Hajnalka Nyitrai to The Division of Medical Sciences in partial fulfillment of the requirements for the degree of Doctor of Philosophy in the subject of Biological and Biomedical Sciences Harvard University Cambridge, Massachusetts August 2019 © 2019 Hajnalka Nyitrai All rights reserved. Dissertation Advisor: Dr. Pascal S. Kaeser Hajnalka Nyitrai Presynaptic Mechanisms for Cargo Recruitment Abstract Delivery and assembly of presynaptic material over long axonal distances is a daunting challenge for neurons. Many presynaptic proteins are thought to be delivered from the Golgi apparatus via precursor vesicles, the recruitment of which at nerve terminals is crucial for the maintenance of presynaptic function. But how do these transport vesicles know where to stop along the axon? In most cellular secretory pathways, Golgi-derived cargos are labeled with Rab GTPases that are recognized by tethering complexes at their fusion sites. Surprisingly, we know very little about what presynaptic tethers and cargo labels facilitate presynaptic cargo capture.
    [Show full text]
  • Computational Proteogenomic Identification and Functional
    bioRxiv preprint doi: https://doi.org/10.1101/168377; this version posted July 25, 2017. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Computational proteogenomic identification and functional interpretation of translated fusions and micro structural variations in cancer Yen Yi Lin1;y, Alexander Gawronski1;y, Faraz Hach1;3;y, Sujun Li2, Ibrahim Numanagic,´ 1, Iman Sarrafi1;3, Swati Mishra5, Andrew McPherson1, Colin Collins3;4, Milan Radovich5, Haixu Tang2, S. Cenk Sahinalp 1;2;3; ∗ 1School of Computing Science, Simon Fraser University, Burnaby, BC, Canada, 2School of Informatics and Computing, Indiana University, Bloomington, IN, USA, 3Vancouver Prostate Centre, Vancouver, BC, Canada, 4Dept. of Urologic Sciences, University of British Columbia, Vancouver, BC, Canada 5Department of Surgery, Indiana University, School of Medicine, Indianapolis, IN, USA Motivation: Rapid advancement in high throughput genome and transcriptome sequencing (HTS) and mass spectrometry (MS) technologies has enabled the acquisition of the genomic, transcriptomic and proteomic data from the same tissue sample. In this paper we in- troduce a novel computational framework which can integratively analyze all three types of omics data to obtain a complete molecular profile of a tissue sample, in normal and disease conditions. Our framework includes MiStrVar, an algorithmic method we developed to identify micro structural variants (microSVs) on genomic HTS data. Coupled with deFuse, a popular gene fusion detection method we devel- oped earlier, MiStrVar can provide an accurate profile of structurally aberrant tran- scripts in cancer samples. Given the breakpoints obtained by MiStrVar and deFuse, our framework can then identify all relevant peptides that span the breakpoint junc- tions and match them with unique proteomic signatures in the respective proteomics data sets.
    [Show full text]
  • 39UTR Shortening Identifies High-Risk Cancers with Targeted Dysregulation
    OPEN 39UTR shortening identifies high-risk SUBJECT AREAS: cancers with targeted dysregulation of GENE REGULATORY NETWORKS the ceRNA network REGULATORY NETWORKS Li Li1*, Duolin Wang2*, Mengzhu Xue1*, Xianqiang Mi1, Yanchun Liang2 & Peng Wang1,3 Received 1 2 15 April 2014 Key Laboratory of Systems Biology, Shanghai Advanced Research Institute, Chinese Academy of Sciences, College of Computer Science and Technology, Jilin University, 3School of Life Science and Technology, ShanghaiTech University. Accepted 3 June 2014 Competing endogenous RNA (ceRNA) interactions form a multilayered network that regulates gene Published expression in various biological pathways. Recent studies have demonstrated novel roles of ceRNA 23 June 2014 interactions in tumorigenesis, but the dynamics of the ceRNA network in cancer remain unexplored. Here, we examine ceRNA network dynamics in prostate cancer from the perspective of alternative cleavage and polyadenylation (APA) and reveal the principles of such changes. Analysis of exon array data revealed that both shortened and lengthened 39UTRs are abundant. Consensus clustering with APA data stratified Correspondence and cancers into groups with differing risks of biochemical relapse and revealed that a ceRNA subnetwork requests for materials enriched with cancer genes was specifically dysregulated in high-risk cancers. The novel connection between should be addressed to 39UTR shortening and ceRNA network dysregulation was supported by the unusually high number of P.W. (wangpeng@ microRNA response elements (MREs) shared by the dysregulated ceRNA interactions and the significantly sari.ac.cn) altered 39UTRs. The dysregulation followed a fundamental principle in that ceRNA interactions connecting genes that show opposite trends in expression change are preferentially dysregulated. This targeted dysregulation is responsible for the majority of the observed expression changes in genes with significant * These authors ceRNA dysregulation and represents a novel mechanism underlying aberrant oncogenic expression.
    [Show full text]
  • A Systematic Genome-Wide Association Analysis for Inflammatory Bowel Diseases (IBD)
    A systematic genome-wide association analysis for inflammatory bowel diseases (IBD) Dissertation zur Erlangung des Doktorgrades der Mathematisch-Naturwissenschaftlichen Fakultät der Christian-Albrechts-Universität zu Kiel vorgelegt von Dipl.-Biol. ANDRE FRANKE Kiel, im September 2006 Referent: Prof. Dr. Dr. h.c. Thomas C.G. Bosch Koreferent: Prof. Dr. Stefan Schreiber Tag der mündlichen Prüfung: Zum Druck genehmigt: der Dekan “After great pain a formal feeling comes.” (Emily Dickinson) To my wife and family ii Table of contents Abbreviations, units, symbols, and acronyms . vi List of figures . xiii List of tables . .xv 1 Introduction . .1 1.1 Inflammatory bowel diseases, a complex disorder . 1 1.1.1 Pathogenesis and pathophysiology. .2 1.2 Genetics basis of inflammatory bowel diseases . 6 1.2.1 Genetic evidence from family and twin studies . .6 1.2.2 Single nucleotide polymorphisms (SNPs) . .7 1.2.3 Linkage studies . .8 1.2.4 Association studies . 10 1.2.5 Known susceptibility genes . 12 1.2.5.1 CARD15. .12 1.2.5.2 CARD4. .15 1.2.5.3 TNF-α . .15 1.2.5.4 5q31 haplotype . .16 1.2.5.5 DLG5 . .17 1.2.5.6 TNFSF15 . .18 1.2.5.7 HLA/MHC on chromosome 6 . .19 1.2.5.8 Other proposed IBD susceptibility genes . .20 1.2.6 Animal models. 21 1.3 Aims of this study . 23 2 Methods . .24 2.1 Laboratory information management system (LIMS) . 24 2.2 Recruitment. 25 2.3 Sample preparation. 27 2.3.1 DNA extraction from blood. 27 2.3.2 Plate design .
    [Show full text]
  • Polyglutamine Expanded Huntingtin Dramatically Alters the Genome-Wide Binding of HSF1
    Polyglutamine Expanded Huntingtin Dramatically Alters the Genome-Wide Binding of HSF1 The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation Riva, Laura, et al. "Polyglutamine Expanded Huntingtin Dramatically Alters the Genome-Wide Binding of HSF1." Journal of Huntington's Disease, Vol. 1, No. 1 (2012). As Published http://dx.doi.org/10.3233/JHD-2012-120020 Publisher IOS Press Version Author's final manuscript Citable link http://hdl.handle.net/1721.1/88924 Terms of Use Creative Commons Attribution-Noncommercial-Share Alike Detailed Terms http://creativecommons.org/licenses/by-nc-sa/4.0/ NIH Public Access Author Manuscript J Huntingtons Dis. Author manuscript; available in PMC 2013 January 04. Published in final edited form as: J Huntingtons Dis. 2012 ; 1(1): 33–45. Poly-glutamine expanded huntingtin dramatically alters the genome wide binding of HSF1 $watermark-text $watermark-text $watermark-text Laura Rivaa,d,†, Martina Koevaa,b,†, Ferah Yildirima, Leila Pirhajia, Deepika Dinesha, Tali Mazora, Martin L. Duennwaldc, and Ernest Fraenkela,* aDepartment of Biological Engineering, Massachusetts Institute of Technology, 77 Mass Ave., Cambridge, MA, USA bThe Whitehead Institute for Biomedical Research, Cambridge, MA, USA cBoston Biomedical Research Institute, Watertown, MA, USA Abstract In Huntington’s disease (HD), polyglutamine expansions in the huntingtin (Htt) protein cause subtle changes in cellular functions that, over-time, lead to neurodegeneration and death. Studies have indicated that activation of the heat shock response can reduce many of the effects of mutant Htt in disease models, suggesting that the heat shock response is impaired in the disease.
    [Show full text]