DOCSLIB.ORG

IRBIS: a Systematic Search for Conserved Complementarity

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
IRBIS: a Systematic Search for Conserved Complementarity Downloaded from rnajournal.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press BIOINFORMATICS IRBIS: a systematic search for conserved complementarity DMITRI D. PERVOUCHINE1,2 1Centre for Genomic Regulation and UPF, Barcelona 08003, Spain 2Faculty of Bioengineering and Bioinformatics, Moscow State University, 119992 Moscow, Russia ABSTRACT IRBIS is a computational pipeline for detecting conserved complementary regions in unaligned orthologous sequences. Unlike other methods, it follows the “first-fold-then-align” principle in which all possible combinations of complementary k-mers are searched for simultaneous conservation. The novel trimming procedure reduces the size of the search space and improves the performance to the point where large-scale analyses of intra- and intermolecular RNA–RNA interactions become possible. In this article, I provide a rigorous description of the method, benchmarking on simulated and real data, and a set of stringent predictions of intramolecular RNA structure in placental mammals, drosophilids, and nematodes. I discuss two particular cases of long-range RNA structures that are likely to have a causal effect on single- and multiple-exon skipping, one in the mammalian gene Dystonin and the other in the insect gene Ca-α1D.InDystonin, one of the two complementary boxes contains a binding site of Rbfox protein similar to one recently described in Enah gene. I also report that snoRNAs and long noncoding RNAs (lncRNAs) have a high capacity of base-pairing to introns of protein-coding genes, suggesting possible involvement of these transcripts in splicing regulation. I also find that conserved sequences that occur equally likely on both strands of DNA (e.g., transcription factor binding sites) contribute strongly to the false-discovery rate and, therefore, would confound every such analysis. IRBIS is an open-source software that is available at http://genome.crg.es/~dmitri/irbis/. Keywords: RNA–RNA interaction; evolutionary conservation; long-range RNA structure; exon skipping; alternative splicing; Ca-α1D; Dystonin; snoRNA; lncRNA INTRODUCTION ular interactions are driven by the same molecular forces, this distinction is crucial for algorithms because RSS is historically RNA–RNA interactions (RRIs) received increasing attention assumed to be nested (i.e., unknotted; see discussed below), in recent years, especially in the light of growing evidence for while simultaneous prediction of intra- and intermolecular abundant expression of noncoding RNAs (Ponting et al. base-pairings is equivalent to RNA folding with pseudoknots 2009; Derrien et al. 2012). One current hypothesis is that (Pervouchine 2004; Alkan et al. 2006; Huang et al. 2009). Here RRI could specifically guide some of the regulatory programs I discuss the two problems jointly without assuming that RSS in the RNA processing pathway, similar to what small RNAs is nested and, in particular, ascribe long-range intramolecular do in the post-transcriptional gene silencing and translational base-pairings also to RRI. attenuation. RRI plays a fundamental role in the functioning Both RRI and RSS predictions comprise a broad range of of the spliceosome, where small nuclear RNAs (snRNAs) methods that admit single-sequence (de novo) and multi- interact with each other and with the pre-mRNA by form- ple-sequence (comparative) formulations. Most of the de ing hetero-duplexes (Will and Luhrmann 2011). Not only novo methods are based on thermodynamic energy model, snRNAs do this; for instance, the C/D box snoRNA HBII- which assumes additive contributions to the free energy func- 52 contains a sequence that is complementary to HT R 2C tion from elementary structural units (Mathews et al. 1999). mRNA and affects alternative splicing in this disease-associ- However, the optimization method that is used to find the ated gene (Kishore and Stamm 2006). minimum free energy (dynamic programming) is computa- The problem of RRI prediction is technically very similar to tionally efficient only for nested RSS: For arbitrary pseudo- RNA secondary structure (RSS) prediction, with the major knots, it is NP complete (Lyngsø and Pedersen 2000), and difference being that base pairs both within and between even for the most generic type of pseudoknots, the required RNA molecules are allowed. Although intra- and intermolec- time is O(n6) (Rivas and Eddy 1999). Besides this technical Corresponding author: [email protected] Article published online ahead of print. Article and publication date are at © 2014 Pervouchine This article, published in RNA, is available under a http://www.rnajournal.org/cgi/doi/10.1261/rna.045088.114. Freely available Creative Commons License (Attribution 4.0 International), as described at online through the RNA Open Access option. http://creativecommons.org/licenses/by/4.0/. RNA 20:1519–1531; Published by Cold Spring Harbor Laboratory Press for the RNA Society 1519 Downloaded from rnajournal.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press Pervouchine limitation, there is a fundamental problem that the additive for long, continuous helices occurring in long-range eukary- model is insufficient to describe entropy contribution of loops otic RSS. in molecules with pseudoknots and that important steric and This line of reasoning was elaborated in our recent reports topological limitations also need to be taken into account on conserved long-range RSS in introns of mammalian (Pervouchine 2004). Consequently, RRI prediction methods and insect genes (Raker et al. 2009; Pervouchine et al. avoid intramolecular interactions to be computationally effi- 2012). The strategy, which shares some technical ideas with cient (Mückstein et al. 2006; Wenzel et al. 2012). The best GUUGle (Gerlach and Giegerich 2006), was to convert se- current trade-off approach uses precomputed accessibility quences into hash tables that store the location of each k- profiles in addition to free energy scoring of exposed binding mer and to apply set-theoretic intersection (1) with the re- sites (Mückstein et al. 2006; Tafer et al. 2011). Some methods verse complement and (2) across orthologs to detect instances gain additional speed by simplifications to the free energy of simultaneous complementarity and conservation. One im- model, which makes them practicable on a genome-wide scale portant advantage of this method is that poorly conserved re- as, for instance, microRNA target finders, although elimina- gions do not need to be aligned and, on the contrary, the lack tion of the internal RNA structure results in a dramatic in- of conservation becomes useful when assigning statistical sig- crease of false-positive predictions (Rehmsmeier et al. 2004; nificance to conserved “islands” immersed in a nonconserved John et al. 2006; Ragan et al. 2009). intronic background (Raker et al. 2009). By construction, In contrast, comparative methods take advantage of the there are neither constraints on the distance between base evolutionary information to reduce the false-positive rate pairs nor limitations on pseudoknots. This technique, in (Gardner and Giegerich 2004). Simultaneous alignment and fact, applies to a broader range of settings than RSSs near folding, known as the Sankoff algorithm, is computationally splice sites. Figure 1 illustrates a general formulation in which overexpensive (Sankoff 1985). Instead, most existing methods the input is organized as a collection of unaligned orthologous take a so-called “first-align-then-fold” route in which a mul- sequence segments. In particular, such collection could con- tiple sequence alignment (MSA) is analyzed, for instance, as a sist of intronic windows adjacent to orthologous donor and profile by a single-sequence algorithm (Seemann et al. 2010, acceptor splice sites or of orthologous miRNA precursors 2011; Li et al. 2011) or by a probabilistic model (Knudsen and 3′-UTRs in a miRNA target search, etc. The aim is to iden- and Hein 2003; Pedersen et al. 2006; Rivas et al. 2012). This tify all pairs of complementary k-mers that occur in sufficient- approach strongly depends on the accuracy of MSA, and al- ly many orthologous segments—no matter where, since the though some improvement can be achieved by considering positions in unaligned sequences are not matched. suboptimal sequence alignments (Will et al. 2013), the major This setup is implemented here as a computational pipe- limitation remains that MSA does not always exist. The oppo- line called IRBIS (intermolecular RNA interaction search, site, “first-fold-then-align” route has not been systematical- where “B” acknowledges a BLAT-like algorithm) (Kent ly investigated because optimal single-sequence predictions 2002), which is designed to search for conserved, long-range are not accurate enough to build a consistent alignment RRI. The pipeline is fully automated and contains all necessary (Shapiro and Zhang 1990; Hochsmann et al. 2003; Gardner preprocessing steps, including genomic sequence download, and Giegerich 2004). identification of orthologous segments, selection of unique Recent studies on RSSs in eukaryotic genes revealed wide- orthologs, etc. (Supplemental Material). Compared to the spread occurrence of long-range RRI with diverse functions previous analyses (Raker et al. 2009; Pervouchine et al. such as riboswitches (Li and Breaker 2013) and mediators 2012), which were restricted to short segments and to a rela- of exon skipping (Lovci et al. 2013), mutually exclusive tively small number of their combinations, IRBIS
Load more

Recommended publications
  • Histone Isoform H2A1H Promotes Attainment of Distinct Physiological
    Bhattacharya et al. Epigenetics & Chromatin (2017) 10:48 DOI 10.1186/s13072-017-0155-z Epigenetics & Chromatin RESEARCH Open Access Histone isoform H2A1H promotes attainment of distinct physiological states by altering chromatin dynamics Saikat Bhattacharya1,4,6, Divya Reddy1,4, Vinod Jani5†, Nikhil Gadewal3†, Sanket Shah1,4, Raja Reddy2,4, Kakoli Bose2,4, Uddhavesh Sonavane5, Rajendra Joshi5 and Sanjay Gupta1,4* Abstract Background: The distinct functional efects of the replication-dependent histone H2A isoforms have been dem- onstrated; however, the mechanistic basis of the non-redundancy remains unclear. Here, we have investigated the specifc functional contribution of the histone H2A isoform H2A1H, which difers from another isoform H2A2A3 in the identity of only three amino acids. Results: H2A1H exhibits varied expression levels in diferent normal tissues and human cancer cell lines (H2A1C in humans). It also promotes cell proliferation in a context-dependent manner when exogenously overexpressed. To uncover the molecular basis of the non-redundancy, equilibrium unfolding of recombinant H2A1H-H2B dimer was performed. We found that the M51L alteration at the H2A–H2B dimer interface decreases the temperature of melting of H2A1H-H2B by ~ 3 °C as compared to the H2A2A3-H2B dimer. This diference in the dimer stability is also refected in the chromatin dynamics as H2A1H-containing nucleosomes are more stable owing to M51L and K99R substitu- tions. Molecular dynamic simulations suggest that these substitutions increase the number of hydrogen bonds and hydrophobic interactions of H2A1H, enabling it to form more stable nucleosomes. Conclusion: We show that the M51L and K99R substitutions, besides altering the stability of histone–histone and histone–DNA complexes, have the most prominent efect on cell proliferation, suggesting that the nucleosome sta- bility is intimately linked with the physiological efects observed.
    [Show full text]
  • Genome-Wide Analysis of 5-Hmc in the Peripheral Blood of Systemic Lupus Erythematosus Patients Using an Hmedip-Chip
    INTERNATIONAL JOURNAL OF MOLECULAR MEDICINE 35: 1467-1479, 2015 Genome-wide analysis of 5-hmC in the peripheral blood of systemic lupus erythematosus patients using an hMeDIP-chip WEIGUO SUI1*, QIUPEI TAN1*, MING YANG1, QIANG YAN1, HUA LIN1, MINGLIN OU1, WEN XUE1, JIEJING CHEN1, TONGXIANG ZOU1, HUANYUN JING1, LI GUO1, CUIHUI CAO1, YUFENG SUN1, ZHENZHEN CUI1 and YONG DAI2 1Guangxi Key Laboratory of Metabolic Diseases Research, Central Laboratory of Guilin 181st Hospital, Guilin, Guangxi 541002; 2Clinical Medical Research Center, the Second Clinical Medical College of Jinan University (Shenzhen People's Hospital), Shenzhen, Guangdong 518020, P.R. China Received July 9, 2014; Accepted February 27, 2015 DOI: 10.3892/ijmm.2015.2149 Abstract. Systemic lupus erythematosus (SLE) is a chronic, Introduction potentially fatal systemic autoimmune disease characterized by the production of autoantibodies against a wide range Systemic lupus erythematosus (SLE) is a typical systemic auto- of self-antigens. To investigate the role of the 5-hmC DNA immune disease, involving diffuse connective tissues (1) and modification with regard to the onset of SLE, we compared is characterized by immune inflammation. SLE has a complex the levels 5-hmC between SLE patients and normal controls. pathogenesis (2), involving genetic, immunologic and envi- Whole blood was obtained from patients, and genomic DNA ronmental factors. Thus, it may result in damage to multiple was extracted. Using the hMeDIP-chip analysis and valida- tissues and organs, especially the kidneys (3). SLE arises from tion by quantitative RT-PCR (RT-qPCR), we identified the a combination of heritable and environmental influences. differentially hydroxymethylated regions that are associated Epigenetics, the study of changes in gene expression with SLE.
    [Show full text]
  • Environmental Influences on Endothelial Gene Expression
    ENDOTHELIAL CELL GENE EXPRESSION John Matthew Jeff Herbert Supervisors: Prof. Roy Bicknell and Dr. Victoria Heath PhD thesis University of Birmingham August 2012 University of Birmingham Research Archive e-theses repository This unpublished thesis/dissertation is copyright of the author and/or third parties. The intellectual property rights of the author or third parties in respect of this work are as defined by The Copyright Designs and Patents Act 1988 or as modified by any successor legislation. Any use made of information contained in this thesis/dissertation must be in accordance with that legislation and must be properly acknowledged. Further distribution or reproduction in any format is prohibited without the permission of the copyright holder. ABSTRACT Tumour angiogenesis is a vital process in the pathology of tumour development and metastasis. Targeting markers of tumour endothelium provide a means of targeted destruction of a tumours oxygen and nutrient supply via destruction of tumour vasculature, which in turn ultimately leads to beneficial consequences to patients. Although current anti -angiogenic and vascular targeting strategies help patients, more potently in combination with chemo therapy, there is still a need for more tumour endothelial marker discoveries as current treatments have cardiovascular and other side effects. For the first time, the analyses of in-vivo biotinylation of an embryonic system is performed to obtain putative vascular targets. Also for the first time, deep sequencing is applied to freshly isolated tumour and normal endothelial cells from lung, colon and bladder tissues for the identification of pan-vascular-targets. Integration of the proteomic, deep sequencing, public cDNA libraries and microarrays, delivers 5,892 putative vascular targets to the science community.
    [Show full text]
  • Repertoire of Morphable Proteins in an Organism
    Repertoire of morphable proteins in an organism Keisuke Izumi, Eitaro Saho, Ayuka Kutomi, Fumiaki Tomoike and Tetsuji Okada Department of Life Science, Gakushuin University, Tokyo, Japan ABSTRACT All living organisms have evolved to contain a set of proteins with variable physical and chemical properties. Efforts in the field of structural biology have contributed to uncovering the shape and the variability of each component. However, quantification of the variability has been performed mostly by multiple pair-wise comparisons. A set of experimental coordinates for a given protein can be used to define the “morphness/unmorphness”. To understand the evolved repertoire in an organism, here we show the results of global analysis of more than a thousand Escherichia coli proteins, by the recently introduced method, distance scoring analysis (DSA). By collecting a new index “UnMorphness Factor” (UMF), proposed in this study and determined from DSA for each of the proteins, the lowest and the highest boundaries of the experimentally observable structural variation are comprehensibly defined. The distribution plot of UMFs obtained for E. coli represents the first view of a substantial fraction of non-redundant proteome set of an organism, demonstrating how rigid and flexible components are balanced. The present analysis extends to evaluate the growing data from single particle cryo-electron microscopy, providing valuable information on effective interpretation to structural changes of proteins and the supramolecular complexes. Subjects Biochemistry, Bioinformatics, Biophysics, Molecular Biology Keywords Protein, Structure, Crystallography, cryo-EM, E. coli, Human, PDB, Coordinates Submitted 21 October 2019 Accepted 20 January 2020 INTRODUCTION Published 11 February 2020 Self-replicating (living) species are defined by the presence of a genome that is used to Corresponding author Tetsuji Okada, produce a set of proteins and RNAs.
    [Show full text]
  • HIST1H2AK (HIST1H2AG) Human Shrna Lentiviral Particle (Locus ID 8969) Product Data
    OriGene Technologies, Inc. 9620 Medical Center Drive, Ste 200 Rockville, MD 20850, US Phone: +1-888-267-4436 [email protected] EU: [email protected] CN: [email protected] Product datasheet for TL304103V HIST1H2AK (HIST1H2AG) Human shRNA Lentiviral Particle (Locus ID 8969) Product data: Product Type: shRNA Lentiviral Particles Product Name: HIST1H2AK (HIST1H2AG) Human shRNA Lentiviral Particle (Locus ID 8969) Locus ID: 8969 Synonyms: H2A.1b; H2A/p; H2AC13; H2AC15; H2AC16; H2AC17; H2AFP; H2AG; HIST1H2AG; pH2A/f Vector: pGFP-C-shLenti (TR30023) Format: Lentiviral particles RefSeq: NM_021064, NM_021064.1, NM_021064.2, NM_021064.3, NM_021064.4, BC016677, BC016677.1, BC067782, NM_021064.5 Summary: Histones are basic nuclear proteins that are responsible for the nucleosome structure of the chromosomal fiber in eukaryotes. Two molecules of each of the four core histones (H2A, H2B, H3, and H4) form an octamer, around which approximately 146 bp of DNA is wrapped in repeating units, called nucleosomes. The linker histone, H1, interacts with linker DNA between nucleosomes and functions in the compaction of chromatin into higher order structures. This gene is intronless and encodes a replication-dependent histone that is a member of the histone H2A family. Transcripts from this gene lack polyA tails but instead contain a palindromic termination element. This gene is found in the histone microcluster on chromosome 6p21.33. [provided by RefSeq, Aug 2015] shRNA Design: These shRNA constructs were designed against multiple splice variants at this gene locus. To be certain that your variant of interest is targeted, please contact [email protected]. If you need a special design or shRNA sequence, please utilize our custom shRNA service.
    [Show full text]
  • Download Download
    Supplementary Figure S1. Results of flow cytometry analysis, performed to estimate CD34 positivity, after immunomagnetic separation in two different experiments. As monoclonal antibody for labeling the sample, the fluorescein isothiocyanate (FITC)- conjugated mouse anti-human CD34 MoAb (Mylteni) was used. Briefly, cell samples were incubated in the presence of the indicated MoAbs, at the proper dilution, in PBS containing 5% FCS and 1% Fc receptor (FcR) blocking reagent (Miltenyi) for 30 min at 4 C. Cells were then washed twice, resuspended with PBS and analyzed by a Coulter Epics XL (Coulter Electronics Inc., Hialeah, FL, USA) flow cytometer. only use Non-commercial 1 Supplementary Table S1. Complete list of the datasets used in this study and their sources. GEO Total samples Geo selected GEO accession of used Platform Reference series in series samples samples GSM142565 GSM142566 GSM142567 GSM142568 GSE6146 HG-U133A 14 8 - GSM142569 GSM142571 GSM142572 GSM142574 GSM51391 GSM51392 GSE2666 HG-U133A 36 4 1 GSM51393 GSM51394 only GSM321583 GSE12803 HG-U133A 20 3 GSM321584 2 GSM321585 use Promyelocytes_1 Promyelocytes_2 Promyelocytes_3 Promyelocytes_4 HG-U133A 8 8 3 GSE64282 Promyelocytes_5 Promyelocytes_6 Promyelocytes_7 Promyelocytes_8 Non-commercial 2 Supplementary Table S2. Chromosomal regions up-regulated in CD34+ samples as identified by the LAP procedure with the two-class statistics coded in the PREDA R package and an FDR threshold of 0.5. Functional enrichment analysis has been performed using DAVID (http://david.abcc.ncifcrf.gov/)
    [Show full text]
  • For Reference Only
    Datasheet Version: 1.0.0 Revision date: 05 Nov 2020 Histone H3R2me2a Antibody Catalogue No.:abx000029 Western blot analysis of extracts of various cell lines, using Asymmetric Dimethyl-Histone H3- R2 antibody (abx000029). Dot-blot analysis of various methylation peptides using Asymmetric Dimethyl-Histone H3-R2 antibody (abx000029). Immunofluorescence analysis of 293T cells using Asymmetric Dimethyl-Histone H3-R2 antibody (abx000029). Blue: DAPI for nuclear staining. Histone H3R2me2a Antibody is a Rabbit Polyclonal antibody against Histone H3R2me2a. Histones are basic nuclear proteins that are responsible for the nucleosome structure of the chromosomal fiber in eukaryotes. Nucleosomes consist of approximately 146 bp of DNA wrapped around a histone octamer composed of pairs of each of the four core histones (H2A, H2B, H3, and H4). The chromatin fiber is further compacted through the interaction of a linker histone, H1, with the DNA between the nucleosomes to form higher order chromatin structures. This gene is intronless and encodes a member of the histone H3 family. TranscriptsFor from this Reference gene lack polyA tails; instead, they contain a palindromic Only termination element. This gene is located separately from the other H3 genes that are in the histone gene cluster on chromosome 6p22-p21.3. Target: Histone H3R2me2a Clonality: Polyclonal Reactivity: Human, Mouse, Rat v1.0.0 Abbexa Ltd, Cambridge, UK · Phone: +44 1223 755950 · Fax: +44 1223 755951 1 Abbexa LLC, Houston, TX, USA · Phone: +1 832 327 7413 www.abbexa.com · Email: [email protected] Datasheet Version: 1.0.0 Revision date: 05 Nov 2020 Tested Applications: WB, IHC, IF/ICC, IP, ChIP Host: Rabbit Recommended dilutions: WB: 1/500 - 1/2000, IHC: 1/50 - 1/200, IF/ICC: 1/50 - 1/200, IP: 1/50 - 1/200, ChIP: 1/20 - 1/100, ChIPseq: 1/20 - 1/100.
    [Show full text]
  • Anti-Histone H3 (Mono Methyl K4) Polyclonal Antibody (CABT-LM00401) This Product Is for Research Use Only and Is Not Intended for Diagnostic Use
    Anti-Histone H3 (mono methyl K4) polyclonal antibody (CABT-LM00401) This product is for research use only and is not intended for diagnostic use. PRODUCT INFORMATION Antigen Description Modulation of chromatin structure plays an important role in the regulation of transcription in eukaryotes. The nucleosome, made up of DNA wound around eight core histone proteins (two each of H2A, H2B, H3, and H4), is the primary building block of chroma Target Histone H3K4me2 Immunogen A synthetic methylated peptide corresponding to residues surrounding K4 of human histone H3 Isotype IgG Source/Host Rabbit Species Reactivity Human Purification Antibodies were produced by immunizing rabbits and were purified by antigen affinity- chromatography. Conjugate Unconjugated Applications WB,IHC,IF,IP,ChIP Molecular Weight 15kDa Format Buffer: PBS with 0.02% sodium azide, 50% glycerol, pH7.3. Preservative 0.02% Sodium Azide Storage Store at -20°C or -80°C. Avoid freeze / thaw cycles. GENE INFORMATION Gene Name HIST3H3 histone cluster 3, H3 [ Homo sapiens (human) ] Official Symbol HIST3H3 Synonyms HIST3H3; histone cluster 3, H3; H3 histone family, member T , H3FT, histone 3, H3; histone H3.1t; H3/g; H3t; H3/t; histone 3, H3; H3 histone family, member T; H3.4; H3FT; MGC126886; MGC126888; 45-1 Ramsey Road, Shirley, NY 11967, USA Email: [email protected] Tel: 1-631-624-4882 Fax: 1-631-938-8221 1 © Creative Diagnostics All Rights Reserved Entrez Gene ID 8290 Protein Refseq NP_003484 UniProt ID Q16695 Chromosome Location 1q42 Pathway Alcoholism; Cell Cycle; Cellular Senescence; Cellular responses to stress; Chromosome Maintenance; DNA Damage/Telomere Stress Induced Senescence; EGFR1 Signaling Pathway; FSH signaling pathway; Meiosis; Meiotic Recombination; Meiotic Synapsis; Packaging O Function DNA binding; protein binding; protein heterodimerization activity References 1.
    [Show full text]
  • Genome-Wide Screen of Cell-Cycle Regulators in Normal and Tumor Cells
    bioRxiv preprint doi: https://doi.org/10.1101/060350; this version posted June 23, 2016. 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. Genome-wide screen of cell-cycle regulators in normal and tumor cells identifies a differential response to nucleosome depletion Maria Sokolova1, Mikko Turunen1, Oliver Mortusewicz3, Teemu Kivioja1, Patrick Herr3, Anna Vähärautio1, Mikael Björklund1, Minna Taipale2, Thomas Helleday3 and Jussi Taipale1,2,* 1Genome-Scale Biology Program, P.O. Box 63, FI-00014 University of Helsinki, Finland. 2Science for Life laboratory, Department of Biosciences and Nutrition, Karolinska Institutet, SE- 141 83 Stockholm, Sweden. 3Science for Life laboratory, Division of Translational Medicine and Chemical Biology, Department of Medical Biochemistry and Biophysics, Karolinska Institutet, S-171 21 Stockholm, Sweden To identify cell cycle regulators that enable cancer cells to replicate DNA and divide in an unrestricted manner, we performed a parallel genome-wide RNAi screen in normal and cancer cell lines. In addition to many shared regulators, we found that tumor and normal cells are differentially sensitive to loss of the histone genes transcriptional regulator CASP8AP2. In cancer cells, loss of CASP8AP2 leads to a failure to synthesize sufficient amount of histones in the S-phase of the cell cycle, resulting in slowing of individual replication forks. Despite this, DNA replication fails to arrest, and tumor cells progress in an elongated S-phase that lasts several days, finally resulting in death of most of the affected cells.
    [Show full text]
  • Histone-Related Genes Are Hypermethylated in Lung Cancer
    Published OnlineFirst October 1, 2019; DOI: 10.1158/0008-5472.CAN-19-1019 Cancer Genome and Epigenome Research Histone-Related Genes Are Hypermethylated in Lung Cancer and Hypermethylated HIST1H4F Could Serve as a Pan-Cancer Biomarker Shihua Dong1,Wei Li1, Lin Wang2, Jie Hu3,Yuanlin Song3, Baolong Zhang1, Xiaoguang Ren1, Shimeng Ji3, Jin Li1, Peng Xu1, Ying Liang1, Gang Chen4, Jia-Tao Lou2, and Wenqiang Yu1 Abstract Lung cancer is the leading cause of cancer-related deaths lated in all 17 tumor types from TCGA datasets (n ¼ 7,344), worldwide. Cytologic examination is the current "gold stan- which was further validated in nine different types of cancer dard" for lung cancer diagnosis, however, this has low sensi- (n ¼ 243). These results demonstrate that HIST1H4F can tivity. Here, we identified a typical methylation signature of function as a universal-cancer-only methylation (UCOM) histone genes in lung cancer by whole-genome DNA methyl- marker, which may aid in understanding general tumorigen- ation analysis, which was validated by The Cancer Genome esis and improve screening for early cancer diagnosis. Atlas (TCGA) lung cancer cohort (n ¼ 907) and was further confirmed in 265 bronchoalveolar lavage fluid samples with Significance: These findings identify a new biomarker for specificity and sensitivity of 96.7% and 87.0%, respectively. cancer detection and show that hypermethylation of histone- More importantly, HIST1H4F was universally hypermethy- related genes seems to persist across cancers. Introduction to its low specificity, LDCT is far from satisfactory as a screening tool for clinical application, similar to other currently used cancer Lung cancer is one of the most common malignant tumors and biomarkers, such as carcinoembryonic antigen (CEA), neuron- the leading cause of cancer-related deaths worldwide (1, 2).
    [Show full text]
  • A Yeast Phenomic Model for the Influence of Warburg Metabolism on Genetic Buffering of Doxorubicin Sean M
    Santos and Hartman Cancer & Metabolism (2019) 7:9 https://doi.org/10.1186/s40170-019-0201-3 RESEARCH Open Access A yeast phenomic model for the influence of Warburg metabolism on genetic buffering of doxorubicin Sean M. Santos and John L. Hartman IV* Abstract Background: The influence of the Warburg phenomenon on chemotherapy response is unknown. Saccharomyces cerevisiae mimics the Warburg effect, repressing respiration in the presence of adequate glucose. Yeast phenomic experiments were conducted to assess potential influences of Warburg metabolism on gene-drug interaction underlying the cellular response to doxorubicin. Homologous genes from yeast phenomic and cancer pharmacogenomics data were analyzed to infer evolutionary conservation of gene-drug interaction and predict therapeutic relevance. Methods: Cell proliferation phenotypes (CPPs) of the yeast gene knockout/knockdown library were measured by quantitative high-throughput cell array phenotyping (Q-HTCP), treating with escalating doxorubicin concentrations under conditions of respiratory or glycolytic metabolism. Doxorubicin-gene interaction was quantified by departure of CPPs observed for the doxorubicin-treated mutant strain from that expected based on an interaction model. Recursive expectation-maximization clustering (REMc) and Gene Ontology (GO)-based analyses of interactions identified functional biological modules that differentially buffer or promote doxorubicin cytotoxicity with respect to Warburg metabolism. Yeast phenomic and cancer pharmacogenomics data were integrated to predict differential gene expression causally influencing doxorubicin anti-tumor efficacy. Results: Yeast compromised for genes functioning in chromatin organization, and several other cellular processes are more resistant to doxorubicin under glycolytic conditions. Thus, the Warburg transition appears to alleviate requirements for cellular functions that buffer doxorubicin cytotoxicity in a respiratory context.
    [Show full text]
  • Genome-Wide DNA Methylation Analysis Reveals Molecular Subtypes of Pancreatic Cancer
    www.impactjournals.com/oncotarget/ Oncotarget, 2017, Vol. 8, (No. 17), pp: 28990-29012 Research Paper Genome-wide DNA methylation analysis reveals molecular subtypes of pancreatic cancer Nitish Kumar Mishra1 and Chittibabu Guda1,2,3,4 1Department of Genetics, Cell Biology and Anatomy, University of Nebraska Medical Center, Omaha, NE, 68198, USA 2Bioinformatics and Systems Biology Core, University of Nebraska Medical Center, Omaha, NE, 68198, USA 3Department of Biochemistry and Molecular Biology, University of Nebraska Medical Center, Omaha, NE, 68198, USA 4Fred and Pamela Buffet Cancer Center, University of Nebraska Medical Center, Omaha, NE, 68198, USA Correspondence to: Chittibabu Guda, email: [email protected] Keywords: TCGA, pancreatic cancer, differential methylation, integrative analysis, molecular subtypes Received: October 20, 2016 Accepted: February 12, 2017 Published: March 07, 2017 Copyright: Mishra et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC-BY), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. ABSTRACT Pancreatic cancer (PC) is the fourth leading cause of cancer deaths in the United States with a five-year patient survival rate of only 6%. Early detection and treatment of this disease is hampered due to lack of reliable diagnostic and prognostic markers. Recent studies have shown that dynamic changes in the global DNA methylation and gene expression patterns play key roles in the PC development; hence, provide valuable insights for better understanding the initiation and progression of PC. In the current study, we used DNA methylation, gene expression, copy number, mutational and clinical data from pancreatic patients.
    [Show full text]