DNA Methylation Biomarkers for Hepatocellular Carcinoma Guorun Fan1†, Yaqin Tu1†, Cai Chen3, Haiying Sun1, Chidan Wan2* and Xiong Cai2*

Total Page:16

File Type:pdf, Size:1020Kb

DNA Methylation Biomarkers for Hepatocellular Carcinoma Guorun Fan1†, Yaqin Tu1†, Cai Chen3, Haiying Sun1, Chidan Wan2* and Xiong Cai2* Fan et al. Cancer Cell Int (2018) 18:140 https://doi.org/10.1186/s12935-018-0629-5 Cancer Cell International PRIMARY RESEARCH Open Access DNA methylation biomarkers for hepatocellular carcinoma Guorun Fan1†, Yaqin Tu1†, Cai Chen3, Haiying Sun1, Chidan Wan2* and Xiong Cai2* Abstract Background: Aberrant methylation of DNA is a key driver of hepatocellular carcinoma (HCC). In this study, we sought to integrate four cohorts profle datasets to identify such abnormally methylated genes and pathways associated with HCC. Methods: To this end, we downloaded microarray datasets examining gene expression (GSE84402, GSE46408) and gene methylation (GSE73003, GSE57956) from the GEO database. Abnormally methylated diferentially expressed genes (DEGs) were sorted and pathways were analyzed. The String database was then used to perform enrichment and functional analysis of identifed pathways and genes. Cytoscape software was used to create a protein–protein interaction network, and MCODE was used for module analysis. Finally, overall survival analysis of hub genes was performed by the OncoLnc online tool. Results: In total, we identifed 19 hypomethylated highly expressed genes and 14 hypermethylated lowly expressed genes at the screening step, and fnally found six mostly changed hub genes including MAD2L1, CDC20, CCNB1, CCND1, AR and ESR1. Pathway analysis showed that aberrantly methylated-DEGs mainly associated with the cell cycle process, p53 signaling, and MAPK signaling in HCC. After validation in TCGA database, the methylation and expres- sion status of hub genes was signifcantly altered and same with our results. Patients with high expression of MAD2L1, CDC20 and CCNB1 and low expression of CCND1, AR, and ESR1 was associated with shorter overall survival. Conclusions: Taken together, we have identifed novel aberrantly methylated genes and pathways linked to HCC, potentially ofering novel insights into the molecular mechanisms governing HCC progression and serving as novel biomarkers for precision diagnosis and disease treatment. Keywords: Hepatocellular carcinoma, Methylation, Hub genes Background with hepatitis C virus (HCV) [3, 4]. Other risk factors Hepatocellular carcinoma (HCC), an infammation- for developing HCC include exposure to afatoxin, alco- driven disease, is the third deadliest cancer worldwide, hol intake, smoking, and diabetes [5]. Te best curative and HCC prevalence is predicted to continue to rise in treatments to date in early stage HCC patients involve coming years, serving as a major economic burden [1, surgical resection, tumor ablation, and potentially liver 2]. Most cases of HCC occur in developing countries, transplantation [6, 7]. However, the prognosis after cura- such as China, and the leading cause of HCC is infection tive therapy for HCC remains unsatisfactory because with hepatitis B virus (HBV); in contrast, the main cause of a high postoperative recurrence rate. An improved in developed countries, such as the USA, is infection understanding of the basic biology of HCC is needed to enhance prognostic predictions and to enhance thera- peutic efcacy against this deadly disease. *Correspondence: [email protected]; [email protected] Te term epigenetics refers to heritable gene expression †Guorun Fan and Yaqin Tu contributed equally to this work 2 Department of Hepatobiliary Surgery, Union Hospital, Tongji Medical alterations no associated with DNA sequence changes College, Huazhong University of Science and Technology, Wuhan 430022, [8]. Te DNA methylation is closely related to embry- China onic development [9], regulation of gene expression [10], Full list of author information is available at the end of the article © The Author(s) 2018. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creat​iveco​mmons​.org/licen​ses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creat​iveco​mmons​.org/ publi​cdoma​in/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Fan et al. Cancer Cell Int (2018) 18:140 Page 2 of 13 X-chromosome inactivation [11], genomic imprinting threshold set for up and down-regulated genes was a | [12], and genomic stability [13]. Altered DNA methyla- log FC | ≥ 2 and P ≤ 0.05. For gene methylation profling tion such as tumor suppressor gene hypermethylation dataset (GSE73003 and GSE57956), the GEO2R software or oncogene hypomethylation is thought to promote was used to analyze the raw data and to identify difer- tumorigenesis. Genes including P15, P16, Ras asso- entially methylated genes (DMGs). GEO2R is an inter- ciation domain family 1 isoform A (RASSF1A), and Ret- active web tool which allows users to compare diferent inoblastoma 1 are inactivated in HCC due to promoter groups of samples in a GEO series to screen genes that hypermethylation of these genes [14–17]. Given that are diferentially expressed in experimental conditions. methylation is potentially reversible, detection of such Te adjusted P-values (adj. P) and Benjamini and Hoch- aberrant DNA methylation of tumor suppressors and berg false discovery rate were applied to provide a bal- oncogenes in HCC could be useful as a therapeutic ance between discovery of statistically signifcant genes target. and limitations of false-positives. Te DMGs cutof crite- While altered methylation of many genes has been ria were P < 0.05 and |t| > 2. Additionally, one-sided tests demonstrated to date in the context of HCC, a com- were used to categorize the upregulated or downregu- plete interaction network documenting the relation- lated DMGs. Finally, a Venn diagram was used to identify ship between said genes remains to be produced. Te hypomethylated highly-expressed genes and hypermeth- comprehensive analysis of multiple datasets ofers the ylated lowly-expressed genes. power needed to properly identify and assess pertinent pathways and genes mediating the biological processes Functional and pathway enrichment analysis associated with HCC. To this end, we used datasets from microarrays examining gene expression (GSE84402, Gene ontology (GO) analyses and Kyoto Encyclopedia GSE46408) and gene methylation (GSE73003, GSE57956) of Genes and Genomes (KEGG) pathway enrichment to assess HCC gene and epigenetic signatures, allowing analyses were conducted on identifed genes with altered for identifcation of genes and pathways that were both methylation/expression using the Search Tool for the abnormally methylated and diferentially expressed. Retrieval of Interacting Genes (STRING, https ://strin Using a protein–protein interaction network we were g-db.org/). Tis allowed for protein–protein interaction also able to identify key so-called “hub” genes central to (PPI) network generation and functional annotation of these signaling events. Trough this analysis, we believe genes of interest. P < 0.05 was the threshold for statistical it is possible to identify novel diferentially methylated signifcance. genes associated with HCC, ofering key insights into the molecular mechanisms governing HCC development and Generation and analysis of a protein–protein interaction progression. (PPI) network In order to interpret the drivers of carcinogenesis in a Methods meaningful manner, a functional PPI analysis was neces- Microarray data sary. Te STRING database allowed us to generate PPI Gene expression profling datasets GSE84402 and networks for both hypomethylated/and highly-expressed GSE46408, and gene methylation profling datasets genes as well as for hypermethylated and lowly-expressed GSE73003 and GSE57956, were downloaded from the genes. A 0.4 interaction score served as the cutof prior to gene expression omnibus (GEO, https ://www.ncbi.nlm. visualization. Modules were then screened using Molec- nih.gov/geo/) [18, 19]. 14 HCC and 14 normal specimens ular Complex Detection (MCODE) in the Cytoscape were obtained in GSE84402, while 6 HCC and 6 normal software package, with MCODE score > 4 and num- samples were obtained in GSE46408. GSE57956 included ber of nodes > 5. Te top 3 Hub genes were selected by a total of 59 primary HCC tumor samples and 59 adja- CytoHubba app in Cytoscape software. To determine the cent normal tissue samples, while GSE73003 included 40 expression pattern of six hub genes in HCC, we used the paired normal and HCC samples from 20 patients. datasets in the Oncomine (https ://www.oncom ine.org) and UALCAN (http://ualca n.path.uab.edu/) database. Data acquisition and processing Oncomine and UALCAN are an online database consist- Raw gene expression profling datasets GSE84402 and ing of previously published and open-access microarray GSE46408 were downloaded from GEO public reposito- data. Te analysis enables multiple comparisons of gene ries. Data processing was performed using robust multi- expression between diferent studies; the signifcance of array average (RMA) in GeneSpring GX 11.5 (Agilent the gene expression across the available studies was also Technologies Pty Ltd), including background adjustment, taken into account [20, 21]. Te results were fltered by normalization and log transformation of the values. Te selecting hepatocellular carcinoma vs. normal tissue. Fan et al. Cancer Cell Int (2018) 18:140 Page 3 of 13 Hub gene validation two sets of genes, a total of 19 hypomethylated,
Recommended publications
  • Profiling of Cell Cycle Genes of Breast Cells Exposed to Etodolac
    1383-1391.qxd 29/3/2010 11:48 Ì ™ÂÏ›‰·1383 ONCOLOGY REPORTS 23: 1383-1391, 2010 Profiling of cell cycle genes of breast cells exposed to etodolac D. ROY1,2*, G.A. ARASON2*, B. CHOWDHURY3,4, A. MITRA5 and G.M. CALAF6,7 1Department of Natural Sciences, Hostos College of the City University of New York, Bronx, NY; 2Graduate Biotechnology Program, Manhattan College, Riverdale, NY; 3Human Retrovirus Section, Vaccine Branch, National Cancer Institute, National Institutes of Health, Frederick, MD; 4Department of Cellular Injury, Military Casualty Research, Walter Reed Army Institute of Research, Silver Spring, MD; 5Department of Chemistry, The College of Mount Saint Vincent, Riverdale, NY; 6Center for Radiological Research, Columbia University Medical Center, New York, NY, USA; 7Instituto de Alta Investigación, Universidad de Tarapaca, Arica, Chile Received September 29, 2009; Accepted November 6, 2009 DOI: 10.3892/or_00000775 Abstract. Breast cancer represents the second leading cause Introduction of cancer-related deaths in the world. There is increasing evidence that perturbation of cell cycle regulation is an Breast cancer represents the second leading cause of cancer- important contributing factor to various cancer progression related deaths in the United States and other Western countries; stages. There are key checkpoints in the cell cycle involving accounting for about 30-40% of all newly diagnosed cancers various regulatory proteins. The relationship between these (1,2). Hereditary breast cancer accounts for just 5-10% of cell cycle regulatory proteins and cell cycle arrest by cyclo- cases. It is generally believed that a family history of breast oxygenase (COX) inhibitors during neoplastic progression cancer and hormonal imbalances also contributes to the remains largely unknown.
    [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]
  • 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]
  • The Mitotic Checkpoint Protein MAD2 Delivers Monoamine Transporters to Endocytosis
    bioRxiv preprint doi: https://doi.org/10.1101/2021.06.09.447721; this version posted June 10, 2021. 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. 1 The mitotic checkpoint protein MAD2 delivers monoamine transporters to endocytosis. Florian Koban1*, Michael Freissmuth1 https://orcid.org/0000-0003-0912-3114, https://orcid.org/0000-0001-9398-1765 1Institute of Pharmacology and the Gaston H. Glock Research Laboratories for Exploratory Drug Development, Center of Physiology and Pharmacology, Medical University of Vienna, Vienna, Austria. *Corresponding Author: Florian Koban, PhD Institute of Pharmacology and the Gaston H. Glock Research Laboratories for Exploratory Drug Development, Center of Physiology and Pharmacology, Medical University of Vienna. Waehringer Strasse 13A, 1090 Vienna, Austria E-mail: [email protected] Tel.: +43 1 40160 31328 Author contributions: F. K. conceptualized the study; F. K. conducted research; F. K. and M. F. analyzed data; F.K. and M.F. wrote the paper. bioRxiv preprint doi: https://doi.org/10.1101/2021.06.09.447721; this version posted June 10, 2021. 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. 2 ABSTRACT Monoamine transporters retrieve serotonin (SERT), dopamine (DAT) and norepinephrine (NET) from the synaptic cleft. Surface levels of transporters are also regulated by internalization. Clathrin-mediated endocytosis of cargo proteins requires adaptor protein 2 (AP2), which recruits cargo to the nascent clathrin- cage. The transporter C-terminus is required for internalization but lacks an AP2- binding site.
    [Show full text]
  • Genetic and Genomic Analysis of Hyperlipidemia, Obesity and Diabetes Using (C57BL/6J × TALLYHO/Jngj) F2 Mice
    University of Tennessee, Knoxville TRACE: Tennessee Research and Creative Exchange Nutrition Publications and Other Works Nutrition 12-19-2010 Genetic and genomic analysis of hyperlipidemia, obesity and diabetes using (C57BL/6J × TALLYHO/JngJ) F2 mice Taryn P. Stewart Marshall University Hyoung Y. Kim University of Tennessee - Knoxville, [email protected] Arnold M. Saxton University of Tennessee - Knoxville, [email protected] Jung H. Kim Marshall University Follow this and additional works at: https://trace.tennessee.edu/utk_nutrpubs Part of the Animal Sciences Commons, and the Nutrition Commons Recommended Citation BMC Genomics 2010, 11:713 doi:10.1186/1471-2164-11-713 This Article is brought to you for free and open access by the Nutrition at TRACE: Tennessee Research and Creative Exchange. It has been accepted for inclusion in Nutrition Publications and Other Works by an authorized administrator of TRACE: Tennessee Research and Creative Exchange. For more information, please contact [email protected]. Stewart et al. BMC Genomics 2010, 11:713 http://www.biomedcentral.com/1471-2164/11/713 RESEARCH ARTICLE Open Access Genetic and genomic analysis of hyperlipidemia, obesity and diabetes using (C57BL/6J × TALLYHO/JngJ) F2 mice Taryn P Stewart1, Hyoung Yon Kim2, Arnold M Saxton3, Jung Han Kim1* Abstract Background: Type 2 diabetes (T2D) is the most common form of diabetes in humans and is closely associated with dyslipidemia and obesity that magnifies the mortality and morbidity related to T2D. The genetic contribution to human T2D and related metabolic disorders is evident, and mostly follows polygenic inheritance. The TALLYHO/ JngJ (TH) mice are a polygenic model for T2D characterized by obesity, hyperinsulinemia, impaired glucose uptake and tolerance, hyperlipidemia, and hyperglycemia.
    [Show full text]
  • Intrinsic Disorder in PTEN and Its Interactome Confers Structural Plasticity and 63
    OPEN Intrinsic Disorder in PTEN and its SUBJECT AREAS: Interactome Confers Structural Plasticity CELLULAR SIGNALLING NETWORKS and Functional Versatility CANCER GENOMICS Prerna Malaney1*, Ravi Ramesh Pathak1*, Bin Xue3, Vladimir N. Uversky3,4,5 & Vrushank Dave´1,2 SYSTEMS ANALYSIS GENE REGULATORY NETWORKS 1Morsani College of Medicine, Department of Pathology and Cell Biology, 2Department of Molecular Oncology, H. Lee Moffitt Cancer Center and Research Institute, 3Department of Molecular Medicine, 4USF Health Byrd Alzheimer’s Research Institute, Received University of South Florida, Tampa, FL, 33612, USA, 5Institute for Biological Instrumentation, Russian Academy of Sciences, 142290 11 March 2013 Pushchino, Moscow Region, Russia. Accepted 3 June 2013 IDPs, while structurally poor, are functionally rich by virtue of their flexibility and modularity. However, how mutations in IDPs elicit diseases, remain elusive. Herein, we have identified tumor suppressor PTEN as Published an intrinsically disordered protein (IDP) and elucidated the molecular principles by which its intrinsically 20 June 2013 disordered region (IDR) at the carboxyl-terminus (C-tail) executes its functions. Post-translational modifications, conserved eukaryotic linear motifs and molecular recognition features present in the C-tail IDR enhance PTEN’s protein-protein interactions that are required for its myriad cellular functions. PTEN primary and secondary interactomes are also enriched in IDPs, most being cancer related, revealing that Correspondence and PTEN functions emanate from and are nucleated by the C-tail IDR, which form pliable network-hubs. requests for materials Together, PTEN higher order functional networks operate via multiple IDP-IDP interactions facilitated by should be addressed to its C-tail IDR. Targeting PTEN IDR and its interaction hubs emerges as a new paradigm for treatment of V.D.
    [Show full text]
  • Gene Section Review
    Atlas of Genetics and Cytogenetics in Oncology and Haematology OPEN ACCESS JOURNAL AT INIST-CNRS Gene Section Review MAD2L1 (mitotic arrest deficient 2, yeast, human homolog like-1) Elizabeth M. Petty, Kenute Myrie Division of Medical Genetics Departments of Human Genetics and Internal Medicine University of Michigan Medical School 1150 West Medical Center Drive, 4301 MSRB III, Ann Arbor, Michigan 48109- 0638, USA (EMP, KM) Published in Atlas Database: March 2001 Online updated version : http://AtlasGeneticsOncology.org/Genes/MAD2L1ID304.html DOI: 10.4267/2042/37729 This work is licensed under a Creative Commons Attribution-Noncommercial-No Derivative Works 2.0 France Licence. © 2001 Atlas of Genetics and Cytogenetics in Oncology and Haematology Identity Other names: HsMAD2; MAD2; MAD2A HGNC (Hugo): MAD2L1 Shadded boxes (1-5) depict the 5 exons of MAD2L1. The black Location : 4q27 triangle indicates a del A mutation that was found in the CAL51 Local order: As noted on the GM99-GB4 breast cancer cell line. Open triangkes depict the locations of Chromosome 4 map: Position: 548.24 (cR3000) Lod identified sequence variants. Figure is not drawn to scale. score: 1.16 Reference Interval: D4S2945-D4S430 Transcription (115.1-125.1 cM) It is located within the NCBI BAC MAD2L1 has 5 coding exons. No alternative splicing genomic contig: NT_006302.2 which is part of the has been described. Regulation of its transcription in homo sapiens chromosome 4 sequence segment. human cells is currently poorly understood. Note: MAD2L1 was intially (and errounously) mapped by fluorescence in situ hybridization (FISH) to 5q23- Protein q31. Subsequent comprehensive mapping studies using somatic cell hybrid analysis, radiation hybrid (RH) Description mapping, and FISH localized it to 4q27.
    [Show full text]
  • WO 2019/079361 Al 25 April 2019 (25.04.2019) W 1P O PCT
    (12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) (19) World Intellectual Property Organization I International Bureau (10) International Publication Number (43) International Publication Date WO 2019/079361 Al 25 April 2019 (25.04.2019) W 1P O PCT (51) International Patent Classification: CA, CH, CL, CN, CO, CR, CU, CZ, DE, DJ, DK, DM, DO, C12Q 1/68 (2018.01) A61P 31/18 (2006.01) DZ, EC, EE, EG, ES, FI, GB, GD, GE, GH, GM, GT, HN, C12Q 1/70 (2006.01) HR, HU, ID, IL, IN, IR, IS, JO, JP, KE, KG, KH, KN, KP, KR, KW, KZ, LA, LC, LK, LR, LS, LU, LY, MA, MD, ME, (21) International Application Number: MG, MK, MN, MW, MX, MY, MZ, NA, NG, NI, NO, NZ, PCT/US2018/056167 OM, PA, PE, PG, PH, PL, PT, QA, RO, RS, RU, RW, SA, (22) International Filing Date: SC, SD, SE, SG, SK, SL, SM, ST, SV, SY, TH, TJ, TM, TN, 16 October 2018 (16. 10.2018) TR, TT, TZ, UA, UG, US, UZ, VC, VN, ZA, ZM, ZW. (25) Filing Language: English (84) Designated States (unless otherwise indicated, for every kind of regional protection available): ARIPO (BW, GH, (26) Publication Language: English GM, KE, LR, LS, MW, MZ, NA, RW, SD, SL, ST, SZ, TZ, (30) Priority Data: UG, ZM, ZW), Eurasian (AM, AZ, BY, KG, KZ, RU, TJ, 62/573,025 16 October 2017 (16. 10.2017) US TM), European (AL, AT, BE, BG, CH, CY, CZ, DE, DK, EE, ES, FI, FR, GB, GR, HR, HU, ΓΕ , IS, IT, LT, LU, LV, (71) Applicant: MASSACHUSETTS INSTITUTE OF MC, MK, MT, NL, NO, PL, PT, RO, RS, SE, SI, SK, SM, TECHNOLOGY [US/US]; 77 Massachusetts Avenue, TR), OAPI (BF, BJ, CF, CG, CI, CM, GA, GN, GQ, GW, Cambridge, Massachusetts 02139 (US).
    [Show full text]
  • Protein Interactions in the Cancer Proteome† Cite This: Mol
    Molecular BioSystems View Article Online PAPER View Journal | View Issue Small-molecule binding sites to explore protein– protein interactions in the cancer proteome† Cite this: Mol. BioSyst., 2016, 12,3067 David Xu,ab Shadia I. Jalal,c George W. Sledge Jr.d and Samy O. Meroueh*aef The Cancer Genome Atlas (TCGA) offers an unprecedented opportunity to identify small-molecule binding sites on proteins with overexpressed mRNA levels that correlate with poor survival. Here, we analyze RNA-seq and clinical data for 10 tumor types to identify genes that are both overexpressed and correlate with patient survival. Protein products of these genes were scanned for binding sites that possess shape and physicochemical properties that can accommodate small-molecule probes or therapeutic agents (druggable). These binding sites were classified as enzyme active sites (ENZ), protein–protein interaction sites (PPI), or other sites whose function is unknown (OTH). Interestingly, the overwhelming majority of binding sites were classified as OTH. We find that ENZ, PPI, and OTH binding sites often occurred on the same structure suggesting that many of these OTH cavities can be used for allosteric modulation of Creative Commons Attribution 3.0 Unported Licence. enzyme activity or protein–protein interactions with small molecules. We discovered several ENZ (PYCR1, QPRT,andHSPA6)andPPI(CASC5, ZBTB32,andCSAD) binding sites on proteins that have been seldom explored in cancer. We also found proteins that have been extensively studied in cancer that have not been previously explored with small molecules that harbor ENZ (PKMYT1, STEAP3,andNNMT) and PPI (HNF4A, MEF2B,andCBX2) binding sites. All binding sites were classified by the signaling pathways to Received 29th March 2016, which the protein that harbors them belongs using KEGG.
    [Show full text]
  • Identi Cation of Novel Hub Genes Associated with Gastric Cancer
    Identication of Novel Hub Genes Associated With Gastric Cancer Using Integrated Bioinformatics Analysis Xiao-Qing Lu Department of Breast Surgery, the second hospital of Shanxi Medical University Jia-qian Zhang Shanxi Medical University Second Aliated Hospital https://orcid.org/0000-0001-9042-0305 Jun Qiao Department of Rheumatology, the Second Hospital of Shanxi Medical University Sheng-Xiao Zhang Department of Rheumatology, the Second Hospital of Shanxi Medical University Meng-Ting Qiu Department of Rheumatology, the Second Hospital of Shanxi Medical University Xiang-Rong Liu Department of Breast Surgery, the Second Hospital of Shanxi Medical University Xiao-Xia Chen Department of Breast Surgery, the Second Hospital of Shanxi Medical University Chong Gao Department of Pathology, Brigham and Women's Hospital, Harvard Medical School Huan-Hu Zhang ( [email protected] ) Department of Gastroenterology Shanxi Cancer Hospital Taiyuan, Shanxi 030001 China Research article Keywords: Gastric cancer, Bioinformatics analysis, Differentially expressed genes Posted Date: September 10th, 2020 DOI: https://doi.org/10.21203/rs.3.rs-58756/v1 License: This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License Page 1/25 Abstract Background: Gastric cancer (GC) is one of the most common solid malignant tumors worldwide with a high- recurrence-rate. Identifying the molecular signatures and specic biomarkers of GC might provide novel clues for GC prognosis and targeted therapy. Methods: Gene expression proles were obtained from the ArrayExpress and Gene Expression Omnibus database. Differentially expressed genes (DEGs) were picked out by R software. The hub genes were screened by cytoHubba plugin. Their prognostic values were assessed by Kaplan–Meier survival analyses and the gene expression proling interactive analysis (GEPIA).
    [Show full text]
  • Proteomic Analysis of Ubiquitin Ligase KEAP1 Reveals Associated Proteins That Inhibit NRF2 Ubiquitination
    Published OnlineFirst February 4, 2013; DOI: 10.1158/0008-5472.CAN-12-4400 Cancer Molecular and Cellular Pathobiology Research Proteomic Analysis of Ubiquitin Ligase KEAP1 Reveals Associated Proteins That Inhibit NRF2 Ubiquitination Bridgid E. Hast1, Dennis Goldfarb2, Kathleen M. Mulvaney1, Michael A. Hast4, Priscila F. Siesser1, Feng Yan1, D. Neil Hayes3, and Michael B. Major1,2 Abstract Somatic mutations in the KEAP1 ubiquitin ligase or its substrate NRF2 (NFE2L2) commonly occur in human cancer, resulting in constitutive NRF2-mediated transcription of cytoprotective genes. However, many tumors display high NRF2 activity in the absence of mutation, supporting the hypothesis that alternative mechanisms of pathway activation exist. Previously, we and others discovered that via a competitive binding mechanism, the proteins WTX (AMER1), PALB2, and SQSTM1 bind KEAP1 to activate NRF2. Proteomic analysis of the KEAP1 protein interaction network revealed a significant enrichment of associated proteins containing an ETGE amino acid motif, which matches the KEAP1 interaction motif found in NRF2. Like WTX, PALB2, and SQSTM1, we found that the dipeptidyl peptidase 3 (DPP3) protein binds KEAP1 via an "ETGE" motif to displace NRF2, thus inhibiting NRF2 ubiquitination and driving NRF2-dependent transcription. Comparing the spectrum of KEAP1-interacting proteins with the genomic profile of 178 squamous cell lung carcinomas characterized by The Cancer Genome Atlas revealed amplification and mRNA overexpression of the DPP3 gene in tumors with high NRF2 activity but lacking NRF2 stabilizing mutations. We further show that tumor-derived mutations in KEAP1 are hypomorphic with respect to NRF2 inhibition and that DPP3 overexpression in the presence of these mutants further promotes NRF2 activation.
    [Show full text]
  • Expression and Mutational Analyses of the Human MAD2L1 Gene in Breast Cancer Cells
    GENES,CHROMOSOMES&CANCER29:356–362(2000) BRIEFCOMMUNICATION ExpressionandMutationalAnalysesoftheHuman MAD2L1GeneinBreastCancerCells MelanieJ.Percy,1 KenuteA.Myrie,2 ChristopherK.Neeley,1 JamesN.Azim,1 StephenP.Ethier,3 and ElizabethM.Petty1,2* 1DepartmentofInternalMedicine,UniversityofMichiganMedicalCenter,AnnArbor,Michigan 2DepartmentofHumanGenetics,UniversityofMichiganMedicalCenter,AnnArbor,Michigan 3DepartmentofRadiationOncology,UniversityofMichiganComprehensiveCancerCenter,AnnArbor,Michigan Breastcancerisaheterogeneousdisorderinwhichmosttumorsdisplaysomedegreeofaneuploidy,especiallythoseatlater stagesofthedisease.Aneuploidyandassociatedchromosomeinstabilitymaybeimportantintheprogressionofmammary tumorigenesis.Aneuploidyispreventedduringnormalcelldivisioninpartthroughregulationofamitoticspindlecheckpoint wheremitoticarrestpreventssegregationofmisalignedchromosomesintodaughtercellsatanaphase.Mitoticarrestgenes, includingtheMADfamily,whichwasoriginallycharacterizedinyeast,helpregulatenormalfunctionofthemitoticspindle checkpoint.DecreasedexpressionofthehumangeneMAD2L1waspreviouslyreportedinabreastcancercelllineexhibiting chromosomeinstabilityandaneuploidy.ToexplorefurtherthepotentialroleofMAD2L1inbreastcancer,weanalyzed MAD2L1geneexpressionin13minimallytogrosslyaneuploidhumanbreastcancercelllinesandfoundsignificantdifferences ofexpressioninthreelines.SequenceanalysisofMAD2L1cDNAintheseaswellasnineadditionalaneuploidbreastcancer andfiveimmortalizednormalhumanmammaryepithelialcelllinesrevealedoneheterozygousframeshift(572delA)mutation inacancercelllinethatdemonstratedahighleveloftranscriptexpression.Inaddition,two3ЈUTRsequencevariantswere
    [Show full text]