Dynamic Expression of M6a Regulators During Multiple Human Tissue Development and Cancers

Total Page:16

File Type:pdf, Size:1020Kb

Dynamic Expression of M6a Regulators During Multiple Human Tissue Development and Cancers ORIGINAL RESEARCH published: 26 January 2021 doi: 10.3389/fcell.2020.629030 Dynamic Expression of m6A Regulators During Multiple Human Tissue Development and Cancers Ya Zhang 1,2†, Sicong Xu 1,2†, Gang Xu 1,2, Yueying Gao 1,2, Si Li 1,2, Ke Zhang 1,2, Zhanyu Tian 1,2, Jing Guo 1,2, Xia Li 1,2*, Juan Xu 1,2* and Yongsheng Li 1,2* 1 Key Laboratory of Tropical Translational Medicine of Ministry of Education, Hainan Medical University, Haikou, China, 2 College of Biomedical Information and Engineering, Hainan Medical University, Haikou, China N6-methyladenosine (m6A) plays critical roles in human development and cancer Edited by: progression. However, our knowledge regarding the dynamic expression of m6A Meng Zhou, Wenzhou Medical University, China regulators during human tissue development is still lacking. Here, we comprehensively 6 Reviewed by: analyzed the dynamic expression alterations of m A regulators during seven tissue Jian Huang, development and eight cancer types. We found that m6A regulators globally exhibited University of Electronic Science and Technology of China, China decreased expression during development. In addition, IGF2BP1/2/3 (insulinlike growth Lixin Cheng, factor 2 MRNA-binding protein 1/2/3) exhibited reverse expression pattern in cancer Jinan University, China progression, suggesting an oncofetal reprogramming in cancer. The expressions Kuixi Zhu, University of Arizona, United States of IGF2BP1/2/3 were regulated by genome alterations, particularly copy number *Correspondence: amplification in cancer. Clinical association analysis revealed that higher expressions of Yongsheng Li IGF2BP1/2/3 were associated with worse survival of cancer patients. Finally, we found [email protected] Juan Xu that genes significantly correlated with IGF2BP1/2/3 were significantly enriched in cancer [email protected] hallmark-related pathways. In summary, dynamic expression analysis will guide both Xia Li mechanistic and therapeutic roles of m6A regulators during tissue development and [email protected] cancer progression. †These authors have contributed 6 equally to this work Keywords: m A regulators, tissue development, cancer, gene experession, survival Specialty section: This article was submitted to INTRODUCTION Epigenomics and Epigenetics, 6 a section of the journal N6-methyladenosine (m A) is the most abundant RNA modification and has been shown to play Frontiers in Cell and Developmental important roles in development and cancer (Dominissini et al., 2012; Meyer et al., 2012). RNA Biology methylation has been shown to be catalyzed by a multicomponent methyltransferase complex, Received: 13 November 2020 including METTL13/14 and WTAP (Liu et al., 2014). RNA methylation has been found in Accepted: 21 December 2020 thousands of coding and non-coding genes. However, we still lack knowledge about the functions Published: 26 January 2021 of m6A regulators in development and cancer. Citation: Numerous studies have shown that m6A regulators were widely perturbed in various types Zhang Y, Xu S, Xu G, Gao Y, Li S, of cancer. METTL3 and IGF2BP2 (insulinlike growth factor 2 MRNA-binding protein 2) were Zhang K, Tian Z, Guo J, Li X, Xu J and overexpressed in colorectal cancer and promote the cancer progression (Li et al., 2019). WTAP was Li Y (2021) Dynamic Expression of 6 also shown to lead to suppression of ETS1 and contribute to the proliferation of liver cancer (Chen m A Regulators During Multiple 6 Human Tissue Development and et al., 2019b). Xu et al. found that 19 m A regulators were highly expressed in glioma tissues, and Cancers. the expression of regulators was associated with prognoses and grade (Xu et al., 2020b). One recent Front. Cell Dev. Biol. 8:629030. study has revealed an essential role of YTHDF1 in the progression of hepatocellular carcinoma doi: 10.3389/fcell.2020.629030 (Liu et al., 2020). RNA methyltransferase METTL3 has been found to facilitate colorectal cancer by Frontiers in Cell and Developmental Biology | www.frontiersin.org 1 January 2021 | Volume 8 | Article 629030 Zhang et al. m6A Regulators in Development and Cancer activating m6A-GLUT1-mTORC1 axis and is a therapeutic preprocessed with RMA normalization, merged, and batch target (Chen et al., 2020). Moreover, we recently performed effect–corrected via Combat method (Leek et al., 2012). comprehensive characterization of m6A regulators in 33 cancer types and found that there were widespread expression 6 6 Collection of m A Regulators perturbation of m A regulators in cancer (Li et al., 2019a). The gene list of m6A regulators was obtained from one of our However, there has been no research that comprehensively 6 studies (Li et al., 2019b), which curated a catalog of 20 genes that explored the expression landscape of m A regulators during function mainly as regulators of RNA methylation. In total, there human tissue development. were 11 readers, 7 writers, and 2 erasers. All these gene symbols In this study, we performed a comprehensive evaluation were converted into Ensemble gene IDs and HGNC symbols of the expression of m6A regulators across multiple human 6 based on annotation from GeneCards (https://www.genecards. tissues development and cancer. We found that the m A org/). regulators exhibited decreased expression after born. Moreover, m6A regulators were up-regulated in cancer, which suggests an oncofetal reprogramming mechanism in cancer. Moreover, Differential Expression Analysis 6 we investigated the association between the expression of m6A To evaluate whether the expression of m A regulators exhibits regulators and patient survival. Our comprehensive analysis differential expression between normal and cancer, we used the of the dynamic expression of m6A regulators provided novel Wilcoxon rank sum test to compare the expression levels. Here, insights into their function in development and cancer. only cancer types with more than five normal samples were analyzed. The gene expression of m6A regulators across normal or cancer samples was shown by box plots. MATERIALS AND METHODS Transcriptome During Human Tissue Genome Alteration of m6A Regulators 6 Development To explore the genome alteration of m A regulators, we used the cBioPortal tool (Gao et al., 2013), which provided a Web resource Gene expression data across human tissue development were for exploring, visualizing, and analyzing multidimensional downloaded from ArrayExpress with the accession code E- cancer genomic data. The gene symbols of m6A regulators were MTAB-6814 (Cardoso-Moreira et al., 2019). The samples were used as input, and TCGA pan-cancer datasets were used. The started from prenatal development at 4 weeks postconception genome alteration frequency of m6A regulators was calculated (WPC) to 20 WPC. Moreover, postnatal samples were also and plotted as bar plot. Here, genetic mutations and copy number sampled, including infants, toddlers, school, teenagers, and adults variation were considered. from each decade until 65 years of age. There were seven human tissues analyzed, including brain, cerebellum, heart, kidney, liver, ovary, and testis. In total, 297 human samples were sequenced Survival Analysis by RNA-Seq. Gene expression levels were calculated as reads per To explore whether the expression of m6A regulators was kilobase of exon model per million mapped reads. associated with patient survival, we divided all the patients into two groups based on the median expression of each m6A Gene Expression in Human Cancers regulator. The log-rank test was used to test the difference We downloaded the gene expression across human cancer survival rates between two groups. This process was performed by the survival package in R program (https://cran.r-project.org/ types from The Cancer Genome Atlas (TCGA) project (Cancer < Genome Atlas Research et al., 2013). According to the tissues in web/packages/survival/index.html). The p 0.05 was considered human development, we considered eight cancer types, including as significant. The odds ratios and 95% confidence levels were kidney renal clear cell carcinoma (KIRC), kidney renal papillary also calculated and shown by ggplot2 package in R program cell carcinoma (KIRP), kidney chromophobe (KICH), brain (Wickham, 2009). lower grade glioma, glioblastoma multiforme, ovarian serous cystadenocarcinoma, liver hepatocellular carcinoma (LIHC), Function Prediction of Regulators and testicular germ cell tumors. Gene expression levels were To predict the function of m6A regulators in human development calculated as fragments per kilobase of exon model per million and cancer, we first calculated the Spearman correlation mapped fragments. Moreover, we downloaded the clinical coefficient (SPCC) between the expression of m6A regulators information of patients, including the survival time, survival and other protein coding genes. All the genes were ranked by status from TCGA project. SPCC and subjected to preranked Gene Set Enrichment Analysis To validate the expression of IGF2BPs regulators across cancer (GSEA) (Mootha et al., 2003; Subramanian et al., 2005). The types, we collected gene expression data across ∼7,400 samples cancer hallmark-related functions were considered (Liberzon representing 11 cancers. Gene expression data were collected et al., 2011, 2015). Gene sets with more than 5 genes and fewer from Gene Expression Omnibus. To minimize interplatform than 1,500 genes were considered. We performed 1,000 times variation, only datasets generated from the Affymetrix Human randomization. The normalized enrichment scores of GSEA were
Recommended publications
  • IGF2BP1 (D-9): Sc-166344
    SANTA CRUZ BIOTECHNOLOGY, INC. IGF2BP1 (D-9): sc-166344 BACKGROUND STORAGE Insulin like growth factor 2 mRNA binding proteins (IGF2BPs) bind RNA and Store at 4° C, **DO NOT FREEZE**. Stable for one year from the date of influence RNA synthesis and metabolism. IGF2BP1, also known as coding shipment. Non-hazardous. No MSDS required. region determinant-binding protein/Insulin-like growth factor II mRNA-binding protein (CRD-BP), IMP1 or VICKZ1; IGF2BP2 (IMP2, VICKZ2, p62); and IGF2BP2 DATA (IMP3, KOC1, VICKZ3) contain a unique combination of RNA recognition motifs A B and four hnRNP K homology domains. IGF2BP1 is abundant in embryonal ABC DEF tissues and is expressed in 81% of colon cancers, 73% of sarcomas and 132 K – 58.5% of breast cancers. It recognizes c-Myc, IGF-II and t mRNAs, and H19 90 K – RNA, and plays a major role in proliferation of K-562 cells by an IGF-II-depen- 55 K – IGF2BP1 dent mechanism. IGF2BP2 binds the 5' UTR of IGF-II mRNA and influences tumor cell growth, in which IGF2BP2 is associated with apoptosis induced 43 K – by tretinoin. IGF2BP3 knockdown by RNA interference decreases levels of IGF-II protein without affecting IGF-II, c-Myc, or b Actin mRNA and H19 RNA levels. IGF2BP3 is a marker for carcinomas and high-grade dysplastic lesions IGF2BP1 (D-9): sc-166344. Western blot analysis IGF2BP1 (D-9): sc-166344. Immunofluorescence staining ofIGF2BP1 expressioninHEK293 (A), NIH/3T3 (B), of methanol-fixed HeLa cells showing nucleolar and of pancreatic ductal epithelium. K-562 (C), U-698-M (D) and Raji (E)wholecell lysates cytoplasmic localization (A).
    [Show full text]
  • Exome Sequencing Revealed DNA Variants in NCOR1, IGF2BP1, SGLT2 and NEK11 As Potential Novel Causes of Ketotic Hypoglycemia in C
    www.nature.com/scientificreports OPEN Exome sequencing revealed DNA variants in NCOR1, IGF2BP1, SGLT2 and NEK11 as potential novel causes of ketotic hypoglycemia in children Yazeid Alhaidan1,2,3,4*, Martin J. Larsen 1,2, Anders Jørgen Schou5, Maria H. Stenlid6, Mohammed A. Al Balwi 3,4, Henrik Thybo Christesen2,5,7,9 & Klaus Brusgaard1,2,8,9 Unexplained or idiopathic ketotic hypoglycemia (KH) is the most common type of hypoglycemia in children. The diagnosis is based on the exclusion of routine hormonal and metabolic causes of hypoglycemia. We aimed to identify novel genes that cause KH, as this may lead to a more targeted treatment. Deep phenotyping of ten preschool age at onset KH patients (boys, n = 5; girls, n = 5) was performed followed by trio exome sequencing and comprehensive bioinformatics analysis. Data analysis revealed four novel candidate genes: (1) NCOR1 in a patient with KH, iron defciency and loose stools; (2) IGF2BP1 in a proband with KH, short stature and delayed bone age; (3) SLC5A2 in a proband with KH, intermittent glucosuria and extremely elevated p-GLP-1; and (4) NEK11 in a proband with ketotic hypoglycemia and liver afiction. These genes are associated with diferent metabolic processes, such as gluconeogenesis, translational regulation, and glucose transport. In conclusion, WES identifed DNA variants in four diferent genes as potential novel causes of IKH, suggesting that IKH is a heterogeneous disorder that can be split into several novel diseases: NCOR1-KH, IGF2BP1-KH, SGLT2-KH or familial renal glucosuria KH, and NEK11-KH. Precision medicine treatment based on exome sequencing may lead to advances in the management of IKH.
    [Show full text]
  • The M A-Related Gene Signature for Predicting the Prognosis of Breast Cancer
    The m6A-related gene signature for predicting the prognosis of breast cancer Shanliang Zhong1,*, Zhenzhong Lin2,*, Huanwen Chen3, Ling Mao4, Jifeng Feng5 and Siying Zhou6 1 Center of Clinical Laboratory Science, The Affiliated Cancer Hospital of Nanjing Medical University & Jiangsu Cancer Hospital & Jiangsu Institute of Cancer Research, Nanjing, China 2 Department of Pathology, The Affiliated Cancer Hospital of Nanjing Medical University & Jiangsu Cancer Hospital & Jiangsu Institute of Cancer Research, Nanjing, China 3 Xinglin laboratory, The First Affiliated Hospital of Xiamen University, Nanjing, China 4 Department of Thyroid Breast Surgery, The Affiliated Huai'an Hospital of Xuzhou Medical University, Huai'an, China 5 Department of Medical Oncology, The Affiliated Cancer Hospital of Nanjing Medical University & Jiangsu Cancer Hospital & Jiangsu Institute of Cancer Research, Nanjing, China 6 Department of General Surgery, The First Affiliated Hospital of Soochow University, Suzhou, China * These authors contributed equally to this work. ABSTRACT N6-methyladenosine (m6A) modification has been shown to participate in tumorigen- esis and metastasis of human cancers. The present study aimed to investigate the roles of m6A RNA methylation regulators in breast cancer. We used LASSO regression to identify m6A-related gene signature predicting breast cancer survival with the datasets downloaded from Gene Expression Omnibus and The Cancer Genome Atlas (TCGA). RNA-Seq data of 3409 breast cancer patients from GSE96058 and 1097 from TCGA were used in present study. A 10 m6A-related gene signature associated with prognosis was identified from 22 m6A RNA methylation regulators. The signature divided patients into low- and high-risk group. High-risk patients had a worse prognosis than the low-risk group.
    [Show full text]
  • IGF2BP1 Is a Targetable SRC/MAPK-Dependent Driver Of
    bioRxiv preprint doi: https://doi.org/10.1101/2020.06.19.159905; this version posted June 20, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. 1 IGF2BP1 is a targetable SRC/MAPK-dependent driver of invasive growth in ovarian 2 cancer 3 4 Authors: Nadine Bley1*$#, Annekatrin Schott1*, Simon Müller1, Danny Misiak1, Marcell 5 Lederer1, Tommy Fuchs1, Chris Aßmann1, Markus Glaß1, Christian Ihling2, Andrea Sinz2, 6 Nikolaos Pazaitis3, Claudia Wickenhauser3, Martina Vetter4, Olga Ungurs4, Hans-Georg 7 Strauss4, Christoph Thomssen4, Stefan Hüttelmaier1$. 8 9 10 Addresses: 11 1 Sect. Molecular Cell Biology, Inst. of Molecular Medicine, Charles Tanford Protein Center, 12 Medical Faculty, Martin Luther University Halle-Wittenberg, Kurt-Mothes-Str. 3A, 06120 Halle, 13 Germany 14 2 Dept. of Pharmaceutical Chemistry & Bioanalytics, Inst. of Pharmacy, Charles Tanford 15 Protein Center, Martin Luther University Halle-Wittenberg, Kurt-Mothes Str. 3A, 06120 Halle, 16 Germany 17 3 Inst. of Pathology, Medical Faculty, Martin Luther University Halle-Wittenberg, Magdeburger 18 Str. 14, 06112 Halle, Germany 19 4 Clinics for Gynecology, Medical Faculty, Martin Luther University Halle-Wittenberg, Ernst- 20 Grube-Str. 40, 06120 Halle, Germany 21 22 23 * These authors contributed equally to this work 24 $ shared corresponding authorship 25 # Correspondence: should be addressed to Nadine Bley 26 Running title: IGF2BP1-directed AJ disassembly depends on SRC activation 27 Disclosure of Potential Conflicts of Interests: the authors declare no conflicts of interest.
    [Show full text]
  • Proteomic Profiling of Retinoblastoma by High Resolution Mass Spectrometry Ravikanth Danda1,5, Kalaivani Ganapathy1, Gajanan Sathe4, Anil K
    Danda et al. Clin Proteom (2016) 13:29 DOI 10.1186/s12014-016-9128-7 Clinical Proteomics RESEARCH Open Access Proteomic profiling of retinoblastoma by high resolution mass spectrometry Ravikanth Danda1,5, Kalaivani Ganapathy1, Gajanan Sathe4, Anil K. Madugundu4, Sharavan Ramachandran1, Uma Maheswari Krishnan5, Vikas Khetan3, Pukhraj Rishi3, T. S. Keshava Prasad4, Akhilesh Pandey6,7,8,9, Subramanian Krishnakumar1*, Harsha Gowda4* and Sailaja V. Elchuri2* Abstract Background: Retinoblastoma is an ocular neoplastic cancer caused primarily due to the mutation/deletion of RB1 gene. Due to the rarity of the disease very limited information is available on molecular changes in primary retino- blastoma. High throughput analysis of retinoblastoma transcriptome is available however the proteomic landscape of retinoblastoma remains unexplored. In the present study we used high resolution mass spectrometry-based quantita- tive proteomics to identify proteins associated with pathogenesis of retinoblastoma. Methods: We used five pooled normal retina and five pooled retinoblastoma tissues to prepare tissue lysates. Equivalent amount of proteins from each group was trypsin digested and labeled with iTRAQ tags. The samples were analyzed on Orbitrap Velos mass spectrometer. We further validated few of the differentially expressed proteins by immunohistochemistry on primary tumors. Results: We identified and quantified a total of 3587 proteins in retinoblastoma when compared with normal adult retina. In total, we identified 899 proteins that were differentially expressed in retinoblastoma with a fold change of 2 of which 402 proteins were upregulated and 497 were down regulated. Insulin growth factor 2 mRNA binding protein≥ 1 (IGF2BP1), chromogranin A, fetuin A (ASHG), Rac GTPase-activating protein 1 and midkine that were found to be overexpressed in retinoblastoma were further confirmed by immunohistochemistry by staining 15 independent retinoblastoma tissue sections.
    [Show full text]
  • 'Next- Generation' Sequencing Data Analysis
    Novel Algorithm Development for ‘Next- Generation’ Sequencing Data Analysis Agne Antanaviciute Submitted in accordance with the requirements for the degree of Doctor of Philosophy University of Leeds School of Medicine Leeds Institute of Biomedical and Clinical Sciences 12/2017 ii The candidate confirms that the work submitted is her own, except where work which has formed part of jointly-authored publications has been included. The contribution of the candidate and the other authors to this work has been explicitly given within the thesis where reference has been made to the work of others. This copy has been supplied on the understanding that it is copyright material and that no quotation from the thesis may be published without proper acknowledgement ©2017 The University of Leeds and Agne Antanaviciute The right of Agne Antanaviciute to be identified as Author of this work has been asserted by her in accordance with the Copyright, Designs and Patents Act 1988. Acknowledgements I would like to thank all the people who have contributed to this work. First and foremost, my supervisors Dr Ian Carr, Professor David Bonthron and Dr Christopher Watson, who have provided guidance, support and motivation. I could not have asked for a better supervisory team. I would also like to thank my collaborators Dr Belinda Baquero and Professor Adrian Whitehouse for opening new, interesting research avenues. A special thanks to Dr Belinda Baquero for all the hard wet lab work without which at least half of this thesis would not exist. Thanks to everyone at the NGS Facility – Carolina Lascelles, Catherine Daley, Sally Harrison, Ummey Hany and Laura Crinnion – for the generation of NGS data used in this work and creating a supportive and stimulating work environment.
    [Show full text]
  • A Novel Hypoxic Long Noncoding RNA KB-1980E6.3 Maintains Breast Cancer Stem Cell Stemness Via Interacting with IGF2BP1 to Facilitate C-Myc Mrna Stability
    Oncogene (2021) 40:1609–1627 https://doi.org/10.1038/s41388-020-01638-9 ARTICLE A novel hypoxic long noncoding RNA KB-1980E6.3 maintains breast cancer stem cell stemness via interacting with IGF2BP1 to facilitate c-Myc mRNA stability 1 2 3 4 1 1 1 1 1 Pengpeng Zhu ● Fang He ● Yixuan Hou ● Gang Tu ● Qiao Li ● Ting Jin ● Huan Zeng ● Yilu Qin ● Xueying Wan ● 1 1 5 1 Yina Qiao ● Yuxiang Qiu ● Yong Teng ● Manran Liu Received: 15 June 2020 / Revised: 13 November 2020 / Accepted: 18 December 2020 / Published online: 19 January 2021 © The Author(s), under exclusive licence to Springer Nature Limited 2021. This article is published with open access Abstract The hostile hypoxic microenvironment takes primary responsibility for the rapid expansion of breast cancer tumors. However, the underlying mechanism is not fully understood. Here, using RNA sequencing (RNA-seq) analysis, we identified a hypoxia-induced long noncoding RNA (lncRNA) KB-1980E6.3, which is aberrantly upregulated in clinical breast cancer tissues and closely correlated with poor prognosis of breast cancer patients. The enhanced lncRNA KB- 1980E6.3 facilitates breast cancer stem cells (BCSCs) self-renewal and tumorigenesis under hypoxic microenvironment both 1234567890();,: 1234567890();,: in vitro and in vivo. Mechanistically, lncRNA KB-1980E6.3 recruited insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) to form a lncRNA KB-1980E6.3/IGF2BP1/c-Myc signaling axis that retained the stability of c-Myc mRNA through increasing binding of IGF2BP1 with m6A-modified c-Myc coding region instability determinant (CRD) mRNA. In conclusion, we confirm that lncRNA KB-1980E6.3 maintains the stemness of BCSCs through lncRNA KB-1980E6.3/ IGF2BP1/c-Myc axis and suggest that disrupting this axis might provide a new therapeutic target for refractory hypoxic tumors.
    [Show full text]
  • Novel Regulators of the IGF System in Cancer
    biomolecules Review Novel Regulators of the IGF System in Cancer Caterina Mancarella 1, Andrea Morrione 2 and Katia Scotlandi 1,* 1 IRCCS Istituto Ortopedico Rizzoli, Laboratory of Experimental Oncology, 40136 Bologna, Italy; [email protected] 2 Department of Biology, Sbarro Institute for Cancer Research and Molecular Medicine and Center for Biotechnology, College of Science and Technology, Temple University, Philadelphia, PA 19122, USA; [email protected] * Correspondence: [email protected]; Tel.: +39-051-6366-760 Abstract: The insulin-like growth factor (IGF) system is a dynamic network of proteins, which includes cognate ligands, membrane receptors, ligand binding proteins and functional downstream effectors. It plays a critical role in regulating several important physiological processes including cell growth, metabolism and differentiation. Importantly, alterations in expression levels or activa- tion of components of the IGF network are implicated in many pathological conditions including diabetes, obesity and cancer initiation and progression. In this review we will initially cover some general aspects of IGF action and regulation in cancer and then focus in particular on the role of transcriptional regulators and novel interacting proteins, which functionally contribute in fine tuning IGF1R signaling in several cancer models. A deeper understanding of the biological relevance of this network of IGF1R modulators might provide novel therapeutic opportunities to block this system in neoplasia. Keywords: IGF system; cancer; transcriptional regulators; functional regulation; circular RNAs; IGF2BPs; ADAR; DDR1; E-cadherin; decorin Citation: Mancarella, C.; Morrione, A.; Scotlandi, K. Novel Regulators of the IGF System in Cancer. 1. Introduction Biomolecules 2021, 11, 273. https:// doi.org/10.3390/biom11020273 The insulin-like growth factor (IGF) system is a network of ligands, binding proteins and receptors regulating crucial physiological and pathological biological processes.
    [Show full text]
  • Large-Scale Functional Prediction of Individual M6a RNA Methylation Sites from an RNA Co-Methylation Network
    Wu et al. BMC Bioinformatics (2019) 20:223 https://doi.org/10.1186/s12859-019-2840-3 RESEARCHARTICLE Open Access m6Acomet: large-scale functional prediction of individual m6A RNA methylation sites from an RNA co- methylation network Xiangyu Wu1,2, Zhen Wei1,2, Kunqi Chen1,2, Qing Zhang1,4, Jionglong Su5, Hui Liu6, Lin Zhang6* and Jia Meng1,3* Abstract Background: Over one hundred different types of post-transcriptional RNA modifications have been identified in human. Researchers discovered that RNA modifications can regulate various biological processes, and RNA methylation, especially N6-methyladenosine, has become one of the most researched topics in epigenetics. Results: To date, the study of epitranscriptome layer gene regulation is mostly focused on the function of mediator proteins of RNA methylation, i.e., the readers, writers and erasers. There is limited investigation of the functional relevance of individual m6A RNA methylation site. To address this, we annotated human m6A sites in large-scale based on the guilt-by-association principle from an RNA co-methylation network. It is constructed based on public human MeRIP-Seq datasets profiling the m6A epitranscriptome under 32 independent experimental conditions. By systematically examining the network characteristics obtained from the RNA methylation profiles, a total of 339,158 putative gene ontology functions associated with 1446 human m6A sites were identified. These are biological functions that may be regulated at epitranscriptome layer via reversible m6A RNA methylation. The results were further validated on a soft benchmark by comparing to a random predictor. Conclusions: An online web server m6Acomet was constructed to support direct query for the predicted biological functions of m6A sites as well as the sites exhibiting co-methylated patterns at the epitranscriptome layer.
    [Show full text]
  • The Devil Is in the Domain: Understanding Protein Recognition of Multiple RNA Targets
    The devil is in the domain: understanding protein recognition of multiple RNA targets Glen R. Gronland1 and Andres Ramos1 1. Institute of Structural and Molecular Biology, University College London, London WC1E 6XA, UK Correspondence to Andres Ramos: [email protected] Abstract: RNA regulation provides a finely-tuned and highly-coordinated control of gene expression. Regulation is mediated by hundreds to thousands of multi-functional RNA-binding proteins which often interact with large sets of RNAs. In this brief review, we focus on recent work that highlights how the proteins use multiple RNA-binding domains to interact selectively with the different RNA targets. De-convoluting the molecular complexity of the RNA regulatory network is essential to understanding cell differentiation and function, and requires accurate models for protein-RNA recognition and protein target selectivity. We discuss that the structural and molecular understanding of the key determinant of recognition, together with the availability of methods to examine protein-RNA interactions at the transcriptome level, may provide an avenue to establish these models. Abbreviations list: CSD, Cold shock domain dsRBD, Double-stranded RNA-binding domain KH, hnRNP K-homology mRNA, Messenger RNA ncRNA, Non-coding RNA RBD, RNA-binding domain RBP, RNA-binding protein RRM, RNA-recognition motif ZnF, Zinc finger Introduction The combined regulation of the various steps in the metabolism and transport of messenger RNAs (mRNAs) and non-coding RNAs (ncRNAs) multiplies genomic potential, and allows cellular differentiation and the development of complex organisms. In the cell, functional RNAs are associated with a fluctuating assortment of multi-domain RNA-binding proteins (RBPs), forming integrated complexes, whose composition directs the fate of the transcript1,2.
    [Show full text]
  • B01P) Gene Alias: IME4, M6A, MGC4336, MT-A70, Spo8
    METTL3 purified MaxPab mouse Gene Symbol: METTL3 polyclonal antibody (B01P) Gene Alias: IME4, M6A, MGC4336, MT-A70, Spo8 Catalog Number: H00056339-B01P Gene Summary: This gene encodes the 70 kDa subunit of MT-A which is part of Regulatory Status: For research use only (RUO) N6-adenosine-methyltransferase. This enzyme is involved in the posttranscriptional methylation of internal Product Description: Mouse polyclonal antibody raised adenosine residues in eukaryotic mRNAs, forming against a full-length human METTL3 protein. N6-methyladenosine. [provided by RefSeq] Immunogen: METTL3 (NP_062826.2, 1 a.a. ~ 580 a.a) References: full-length human protein. 1. Synaptic N6-methyladenosine (m6A) epitranscriptome Sequence: reveals functional partitioning of localized transcripts. MSDTWSSIQAHKKQLDSLRERLQRRRKQDSGHLDLR Merkurjev D, Hong WT, Iida K, Oomoto I, Goldie BJ, NPEAALSPTFRSDSPVPTAPTSGGPKPSTASAVPELAT Yamaguti H, Ohara T, Kawaguchi SY, Hirano T, Martin DPELEKKLLHHLSDLALTLPTDAVSICLAISTPDAPATQ KC, Pellegrini M, Wang DO. Nat Neurosci. 2018 DGVESLLQKFAAQELIEVKRGLLQDDAHPTLVTYADHS Jul;21(7):1004-1014. doi: 10.1038/s41593-018-0173-6. KLSAMMGAVAEKKGPGEVAGTVTGQKRRAEQDSTTV Epub 2018 Jun 27. AAFASSLVSGLNSSASEPAKEPAKKSRKHAASDVDLEI 2. The RNA Methyltransferase Complex of WTAP, ESLLNQQSTKEQQSKKVSQEILELLNTTTAKEQSIVEKF METTL3, and METTL14 Regulates Mitotic Clonal RSRGRAQVQEFCDYGTKEECMKASDADRPCRKLHFR Expansion in Adipogenesis. Kobayashi M, Ohsugi M, RIINKHTDESLGDCSFLNTCFHMDTCKYVHYEIDACMD Sasako T, Awazawa M, Umehara T, Iwane A, Kobayashi SEAPGSKDHTPSQELALTQSVGGDSSADRLFPPQWIC
    [Show full text]
  • Zipcode Binding Protein 1 (ZBP1; IGF2BP1): a Model for Sequence-Specific RNA Regulation
    Zipcode Binding Protein 1 (ZBP1; IGF2BP1): A Model for Sequence-Specific RNA Regulation 1 1,5 1,5 1,2 JEETAYU BISWAS, LETI NUNEZ, SULAGNA DAS, YOUNG J. YOON, 1,3 1,2,4 CAROLINA ELISCOVICH, AND ROBERT H. SINGER 1Department of Anatomy and Structural Biology, 2Department of Neuroscience, 3Department of Medicine, Albert Einstein College of Medicine, Bronx, New York 10461, USA 4Howard Hughes Medical Institute, Janelia Research Campus, Ashburn, Virginia 20147, USA Correspondence: [email protected] The fate of an RNA, from its localization, translation, and ultimate decay, is dictated by interactions with RNA binding proteins (RBPs). β-actin mRNA has functioned as the classic example of RNA localization in eukaryotic cells. Studies of β-actin mRNA over the past three decades have allowed understanding of how RBPs, such as ZBP1 (IGF2BP1), can control both RNA localization and translational status. Here, we summarize studies of β-actin mRNA and focus on how ZBP1 serves as a model for understanding interactions between RNA and their binding protein(s). Central to the study of RNA and RBPs were technological developments that occurred along the way. We conclude with a future outlook highlighting new technologies that may be used to address still unanswered questions about RBP-mediated regulation of mRNA during its life cycle, within the cell. Changes in protein synthesis define all aspects of a RNAs tested (actin, vimentin, and tubulin), actin had cell’s life. Although most cells carry identical copy num- the most distinct localization pattern. Approximately bers of their DNA, differences in transcription and trans- 95% of myoblasts had a nonrandom distribution of β-actin lation define both cellular fate and moment to moment mRNA, in which RNAs were concentrated at the cell cellular decisions.
    [Show full text]