1
|
Fonseca A, Riveras E, Moyano TC, Alvarez JM, Rosa S, Gutiérrez RA. Dynamic changes in mRNA nucleocytoplasmic localization in the nitrate response of Arabidopsis roots. PLANT, CELL & ENVIRONMENT 2024. [PMID: 38950037 DOI: 10.1111/pce.15018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Revised: 05/23/2024] [Accepted: 06/14/2024] [Indexed: 07/03/2024]
Abstract
Nitrate is a nutrient and signal that regulates gene expression. The nitrate response has been extensively characterized at the organism, organ, and cell-type-specific levels, but intracellular mRNA dynamics remain unexplored. To characterize nuclear and cytoplasmic transcriptome dynamics in response to nitrate, we performed a time-course expression analysis after nitrate treatment in isolated nuclei, cytoplasm, and whole roots. We identified 402 differentially localized transcripts (DLTs) in response to nitrate treatment. Induced DLT genes showed rapid and transient recruitment of the RNA polymerase II, together with an increase in the mRNA turnover rates. DLTs code for genes involved in metabolic processes, localization, and response to stimulus indicating DLTs include genes with relevant functions for the nitrate response that have not been previously identified. Using single-molecule RNA FISH, we observed early nuclear accumulation of the NITRATE REDUCTASE 1 (NIA1) transcripts in their transcription sites. We found that transcription of NIA1, a gene showing delayed cytoplasmic accumulation, is rapidly and transiently activated; however, its transcripts become unstable when they reach the cytoplasm. Our study reveals the dynamic localization of mRNAs between the nucleus and cytoplasm as an emerging feature in the temporal control of gene expression in response to nitrate treatment in Arabidopsis roots.
Collapse
Affiliation(s)
- Alejandro Fonseca
- Millennium Institute for Integrative Biology (iBio), Santiago, Chile
- Center for Genome Regulation, Millennium Institute Center for Genome Regulation (CRG), Santiago, Chile
- Departamento de Genética Molecular y Microbiología, Facultad de Ciencias Biológicas, Pontificia Universidad Católica de Chile, Santiago, Chile
- Department of Plant Biology, Swedish University of Agricultural Sciences (SLU), Uppsala, Sweden
| | - Eleodoro Riveras
- Millennium Institute for Integrative Biology (iBio), Santiago, Chile
- Center for Genome Regulation, Millennium Institute Center for Genome Regulation (CRG), Santiago, Chile
- Departamento de Genética Molecular y Microbiología, Facultad de Ciencias Biológicas, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Tomás C Moyano
- Millennium Institute for Integrative Biology (iBio), Santiago, Chile
- Center for Genome Regulation, Millennium Institute Center for Genome Regulation (CRG), Santiago, Chile
- Departamento de Genética Molecular y Microbiología, Facultad de Ciencias Biológicas, Pontificia Universidad Católica de Chile, Santiago, Chile
- Centro de Biotecnología Vegetal, Facultad de Ciencias de la Vida, Universidad Andrés Bello, Santiago, Chile
| | - José M Alvarez
- Millennium Institute for Integrative Biology (iBio), Santiago, Chile
- Centro de Biotecnología Vegetal, Facultad de Ciencias de la Vida, Universidad Andrés Bello, Santiago, Chile
| | - Stefanie Rosa
- Department of Plant Biology, Swedish University of Agricultural Sciences (SLU), Uppsala, Sweden
| | - Rodrigo A Gutiérrez
- Millennium Institute for Integrative Biology (iBio), Santiago, Chile
- Center for Genome Regulation, Millennium Institute Center for Genome Regulation (CRG), Santiago, Chile
- Departamento de Genética Molecular y Microbiología, Facultad de Ciencias Biológicas, Pontificia Universidad Católica de Chile, Santiago, Chile
| |
Collapse
|
2
|
Tomaz da Silva P, Zhang Y, Theodorakis E, Martens LD, Yépez VA, Pelechano V, Gagneur J. Cellular energy regulates mRNA degradation in a codon-specific manner. Mol Syst Biol 2024; 20:506-520. [PMID: 38491213 PMCID: PMC11066088 DOI: 10.1038/s44320-024-00026-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2023] [Revised: 02/19/2024] [Accepted: 02/20/2024] [Indexed: 03/18/2024] Open
Abstract
Codon optimality is a major determinant of mRNA translation and degradation rates. However, whether and through which mechanisms its effects are regulated remains poorly understood. Here we show that codon optimality associates with up to 2-fold change in mRNA stability variations between human tissues, and that its effect is attenuated in tissues with high energy metabolism and amplifies with age. Mathematical modeling and perturbation data through oxygen deprivation and ATP synthesis inhibition reveal that cellular energy variations non-uniformly alter the effect of codon usage. This new mode of codon effect regulation, independent of tRNA regulation, provides a fundamental mechanistic link between cellular energy metabolism and eukaryotic gene expression.
Collapse
Affiliation(s)
- Pedro Tomaz da Silva
- School of Computation, Information and Technology, Technical University of Munich, Garching, Germany
- Munich Center for Machine Learning, Munich, Germany
| | - Yujie Zhang
- Scilifelab, Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden
| | - Evangelos Theodorakis
- School of Computation, Information and Technology, Technical University of Munich, Garching, Germany
| | - Laura D Martens
- School of Computation, Information and Technology, Technical University of Munich, Garching, Germany
- Computational Health Center, Helmholtz Center Munich, Neuherberg, Germany
| | - Vicente A Yépez
- School of Computation, Information and Technology, Technical University of Munich, Garching, Germany
| | - Vicent Pelechano
- Scilifelab, Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden
| | - Julien Gagneur
- School of Computation, Information and Technology, Technical University of Munich, Garching, Germany.
- Computational Health Center, Helmholtz Center Munich, Neuherberg, Germany.
- Institute of Human Genetics, School of Medicine, Technical University of Munich, Munich, Germany.
| |
Collapse
|
3
|
Liu H, Arsiè R, Schwabe D, Schilling M, Minia I, Alles J, Boltengagen A, Kocks C, Falcke M, Friedman N, Landthaler M, Rajewsky N. SLAM-Drop-seq reveals mRNA kinetic rates throughout the cell cycle. Mol Syst Biol 2023; 19:1-23. [PMID: 38778223 PMCID: PMC10568207 DOI: 10.15252/msb.202211427] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Revised: 07/24/2023] [Accepted: 08/04/2023] [Indexed: 05/25/2024] Open
Abstract
RNA abundance is tightly regulated in eukaryotic cells by modulating the kinetic rates of RNA production, processing, and degradation. To date, little is known about time‐dependent kinetic rates during dynamic processes. Here, we present SLAM‐Drop‐seq, a method that combines RNA metabolic labeling and alkylation of modified nucleotides in methanol‐fixed cells with droplet‐based sequencing to detect newly synthesized and preexisting mRNAs in single cells. As a first application, we sequenced 7280 HEK293 cells and calculated gene‐specific kinetic rates during the cell cycle using the novel package Eskrate. Of the 377 robust‐cycling genes that we identified, only a minor fraction is regulated solely by either dynamic transcription or degradation (6 and 4%, respectively). By contrast, the vast majority (89%) exhibit dynamically regulated transcription and degradation rates during the cell cycle. Our study thus shows that temporally regulated mRNA degradation is fundamental for the correct expression of a majority of cycling genes. SLAM‐Drop‐seq, combined with Eskrate, is a powerful approach to understanding the underlying mRNA kinetics of single‐cell gene expression dynamics in continuous biological processes.
Collapse
Affiliation(s)
- Haiyue Liu
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany
| | - Roberto Arsiè
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany
| | - Daniel Schwabe
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany
| | - Marcel Schilling
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany
- Lübeck Interdisciplinary Platform for Genome Analytics (LIGA), University of Lübeck, Lübeck, Germany
| | - Igor Minia
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany
| | - Jonathan Alles
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany
| | - Anastasiya Boltengagen
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany
| | - Christine Kocks
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany
| | - Martin Falcke
- Mathematical Cell Physiology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany
- Department of Physics, Humboldt University Berlin, Berlin, Germany
| | - Nir Friedman
- The Rachel and Selim Benin School of Computer Science and Engineering, The Hebrew University of Jerusalem, Jerusalem, Israel
- The Lautenberg Center for Immunology and Cancer Research, Institute of Medical Research Israel-Canada (IMRIC), Faculty of Medicine, The Hebrew University of Jerusalem, Jerusalem, Israel
- The Center for Computational Medicine, Faculty of Medicine, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Markus Landthaler
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany.
- Institut für Biologie, Humboldt Universität zu Berlin, Berlin, Germany.
| | - Nikolaus Rajewsky
- Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin, Germany.
- Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany.
- German Center for Cardiovascular Research (DZHK), Berlin, Germany.
- NeuroCure Cluster of Excellence, Berlin, Germany.
- German Cancer Consortium (DKTK), Berlin, Germany.
- National Center for Tumor Diseases (NCT), Berlin, Germany.
| |
Collapse
|
4
|
Sarkar S, Rammohan J. Nearly maximal information gain due to time integration in central dogma reactions. iScience 2023; 26:106767. [PMID: 37235057 PMCID: PMC10206154 DOI: 10.1016/j.isci.2023.106767] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2022] [Revised: 02/21/2023] [Accepted: 04/24/2023] [Indexed: 05/28/2023] Open
Abstract
Living cells process information about their environment through the central dogma processes of transcription and translation, which drive the cellular response to stimuli. Here, we study the transfer of information from environmental input to the transcript and protein expression levels. Evaluation of both experimental and analogous simulation data reveals that transcription and translation are not two simple information channels connected in series. Instead, we demonstrate that the central dogma reactions often create a time-integrating information channel, where the translation channel receives and integrates multiple outputs from the transcription channel. This information channel model of the central dogma provides new information-theoretic selection criteria for the central dogma rate constants. Using the data for four well-studied species we show that their central dogma rate constants achieve information gain because of time integration while also keeping the loss because of stochasticity in translation relatively low (<0.5 bits).
Collapse
Affiliation(s)
- Swarnavo Sarkar
- National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
| | - Jayan Rammohan
- National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
| |
Collapse
|
5
|
Aubé S, Nielly-Thibault L, Landry CR. Evolutionary trade-off and mutational bias could favor transcriptional over translational divergence within paralog pairs. PLoS Genet 2023; 19:e1010756. [PMID: 37235586 DOI: 10.1371/journal.pgen.1010756] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2023] [Accepted: 04/21/2023] [Indexed: 05/28/2023] Open
Abstract
How changes in the different steps of protein synthesis-transcription, translation and degradation-contribute to differences of protein abundance among genes is not fully understood. There is however accumulating evidence that transcriptional divergence might have a prominent role. Here, we show that yeast paralogous genes are more divergent in transcription than in translation. We explore two causal mechanisms for this predominance of transcriptional divergence: an evolutionary trade-off between the precision and economy of gene expression and a larger mutational target size for transcription. Performing simulations within a minimal model of post-duplication evolution, we find that both mechanisms are consistent with the observed divergence patterns. We also investigate how additional properties of the effects of mutations on gene expression, such as their asymmetry and correlation across levels of regulation, can shape the evolution of paralogs. Our results highlight the importance of fully characterizing the distributions of mutational effects on transcription and translation. They also show how general trade-offs in cellular processes and mutation bias can have far-reaching evolutionary impacts.
Collapse
Affiliation(s)
- Simon Aubé
- Département de biochimie, de microbiologie et de bio-informatique, Faculté des sciences et de génie, Université Laval, Québec, Québec, Canada
- Institut de Biologie Intégrative et des Systèmes, Université Laval, Québec, Québec, Canada
- PROTEO, Le regroupement québécois de recherche sur la fonction, l'ingénierie et les applications des protéines, Université Laval, Québec, Québec, Canada
- Centre de Recherche en Données Massives, Université Laval, Québec, Québec, Canada
| | - Lou Nielly-Thibault
- Institut de Biologie Intégrative et des Systèmes, Université Laval, Québec, Québec, Canada
- PROTEO, Le regroupement québécois de recherche sur la fonction, l'ingénierie et les applications des protéines, Université Laval, Québec, Québec, Canada
- Centre de Recherche en Données Massives, Université Laval, Québec, Québec, Canada
- Département de biologie, Faculté des sciences et de génie, Université Laval, Québec, Québec, Canada
| | - Christian R Landry
- Département de biochimie, de microbiologie et de bio-informatique, Faculté des sciences et de génie, Université Laval, Québec, Québec, Canada
- Institut de Biologie Intégrative et des Systèmes, Université Laval, Québec, Québec, Canada
- PROTEO, Le regroupement québécois de recherche sur la fonction, l'ingénierie et les applications des protéines, Université Laval, Québec, Québec, Canada
- Centre de Recherche en Données Massives, Université Laval, Québec, Québec, Canada
- Département de biologie, Faculté des sciences et de génie, Université Laval, Québec, Québec, Canada
| |
Collapse
|
6
|
Nicholas BA, Purohit R, Woods AD, Kannan K, Srinivasa G, Bridge JA, Kim JA, Keller C. BCR-ABL is enriched in S- and G 2-cell cycle phases. Leuk Res 2023; 126:107036. [PMID: 36764024 DOI: 10.1016/j.leukres.2023.107036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2022] [Revised: 02/03/2023] [Accepted: 02/06/2023] [Indexed: 02/11/2023]
Affiliation(s)
- Bradley A Nicholas
- Children's Cancer Therapy Development Institute, Beaverton, OR 97005 USA
| | - Reshma Purohit
- Children's Cancer Therapy Development Institute, Beaverton, OR 97005 USA
| | - Andrew D Woods
- Children's Cancer Therapy Development Institute, Beaverton, OR 97005 USA
| | | | | | | | - Jin-Ah Kim
- Children's Cancer Therapy Development Institute, Beaverton, OR 97005 USA
| | - Charles Keller
- Children's Cancer Therapy Development Institute, Beaverton, OR 97005 USA.
| |
Collapse
|
7
|
Ofir R, Kriecherbauer T, Grüne L, Margaliot M. On the gain of entrainment in the n-dimensional ribosome flow model. J R Soc Interface 2023; 20:20220763. [PMID: 36751928 PMCID: PMC9905980 DOI: 10.1098/rsif.2022.0763] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2022] [Accepted: 01/11/2023] [Indexed: 02/09/2023] Open
Abstract
The ribosome flow model (RFM) is a phenomenological model for the flow of particles along a one-dimensional chain of n sites. It has been extensively used to study ribosome flow along the mRNA molecule during translation. When the transition rates along the chain are time-varying and jointly T-periodic the RFM entrains, i.e. every trajectory of the RFM converges to a unique T-periodic solution that depends on the transition rates, but not on the initial condition. In general, entrainment to periodic excitations like the 24 h solar day or the 50 Hz frequency of the electric grid is important in numerous natural and artificial systems. An interesting question, called the gain of entrainment (GOE) in the RFM, is whether proper coordination of the periodic translation rates along the mRNA can lead to a larger average protein production rate. Analysing the GOE in the RFM is non-trivial and partial results exist only for the RFM with dimensions n = 1, 2. We use a new approach to derive several results on the GOE in the general n-dimensional RFM. Perhaps surprisingly, we rigorously characterize several cases where there is no GOE, so to maximize the average production rate in these cases, the best choice is to use constant transition rates all along the chain.
Collapse
Affiliation(s)
- Ron Ofir
- Andrew and Erna Viterbi Faculty of Electrical and Computer Engineering, Technion—Israel Institute of Technology, Haifa 3200003, Israel
| | | | - Lars Grüne
- Mathematical Institute, University of Bayreuth, Bayreuth, Germany
| | - Michael Margaliot
- School of Electrical Engineering, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
8
|
A dynamical stochastic model of yeast translation across the cell cycle. Heliyon 2023; 9:e13101. [PMID: 36793957 PMCID: PMC9922973 DOI: 10.1016/j.heliyon.2023.e13101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2022] [Revised: 01/04/2023] [Accepted: 01/16/2023] [Indexed: 01/27/2023] Open
Abstract
Translation is a central step in gene expression, however its quantitative and time-resolved regulation is poorly understood. We developed a discrete, stochastic model for protein translation in S. cerevisiae in a whole-transcriptome, single-cell context. A "base case" scenario representing an average cell highlights translation initiation rates as the main co-translational regulatory parameters. Codon usage bias emerges as a secondary regulatory mechanism through ribosome stalling. Demand for anticodons with low abundancy is shown to cause above-average ribosome dwelling times. Codon usage bias correlates strongly both with protein synthesis rates and elongation rates. Applying the model to a time-resolved transcriptome estimated by combining data from FISH and RNA-Seq experiments, it could be shown that increased total transcript abundance during the cell cycle decreases translation efficiency at single transcript level. Translation efficiency grouped by gene function shows highest values for ribosomal and glycolytic genes. Ribosomal proteins peak in S phase while glycolytic proteins rank highest in later cell cycle phases.
Collapse
|
9
|
Conte F, Papa F, Paci P, Farina L. StaRTrEK:in silico estimation of RNA half-lives from genome-wide time-course experiments without transcriptional inhibition. BMC Bioinformatics 2022; 23:190. [PMID: 35596139 PMCID: PMC9123730 DOI: 10.1186/s12859-022-04730-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2021] [Accepted: 05/05/2022] [Indexed: 11/26/2022] Open
Abstract
Background Gene expression is the result of the balance between transcription and degradation. Recent experimental findings have shown fine and specific regulation of RNA degradation and the presence of various molecular machinery purposely devoted to this task, such as RNA binding proteins, non-coding RNAs, etc. A biological process can be studied by measuring time-courses of RNA abundance in response of internal and/or external stimuli, using recent technologies, such as the microarrays or the Next Generation Sequencing devices. Unfortunately, the picture provided by looking only at the transcriptome abundance may not gain insight into its dynamic regulation. By contrast, independent simultaneous measurement of RNA expression and half-lives could provide such valuable additional insight. A computational approach to the estimation of RNAs half-lives from RNA expression time profiles data, can be a low-cost alternative to its experimental measurement which may be also affected by various artifacts. Results Here we present a computational methodology, called StaRTrEK (STAbility Rates ThRough Expression Kinetics), able to estimate half-life values basing only on genome-wide gene expression time series without transcriptional inhibition. The StaRTrEK algorithm makes use of a simple first order kinetic model and of a \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$l_1$$\end{document}l1-norm regularized least square optimization approach to find its parameter values. Estimates provided by StaRTrEK are validated using simulated data and three independent experimental datasets of two short (6 samples) and one long (48 samples) time-courses. Conclusions We believe that our algorithm can be used as a fast valuable computational complement to time-course experimental gene expression studies by adding a relevant kinetic property, i.e. the RNA half-life, with a strong biological interpretation, thus providing a dynamic picture of what is going in a cell during the biological process under study. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-022-04730-x.
Collapse
Affiliation(s)
- Federica Conte
- Institute for Systems Analysis and Computer Science "Antonio Ruberti", National Research Council, Rome, Italy
| | - Federico Papa
- Institute for Systems Analysis and Computer Science "Antonio Ruberti", National Research Council, Rome, Italy.,SysBio Centre for Systems Biology, Milan, Italy
| | - Paola Paci
- Department of Computer, Control and Management Engineering, Sapienza University of Rome, Rome, Italy
| | - Lorenzo Farina
- Department of Computer, Control and Management Engineering, Sapienza University of Rome, Rome, Italy.
| |
Collapse
|
10
|
Monti M, Fiorentino J, Milanetti E, Gosti G, Tartaglia GG. Prediction of Time Series Gene Expression and Structural Analysis of Gene Regulatory Networks Using Recurrent Neural Networks. ENTROPY (BASEL, SWITZERLAND) 2022; 24:141. [PMID: 35205437 PMCID: PMC8871363 DOI: 10.3390/e24020141] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Revised: 01/14/2022] [Accepted: 01/15/2022] [Indexed: 11/17/2022]
Abstract
Methods for time series prediction and classification of gene regulatory networks (GRNs) from gene expression data have been treated separately so far. The recent emergence of attention-based recurrent neural network (RNN) models boosted the interpretability of RNN parameters, making them appealing for the understanding of gene interactions. In this work, we generated synthetic time series gene expression data from a range of archetypal GRNs and we relied on a dual attention RNN to predict the gene temporal dynamics. We show that the prediction is extremely accurate for GRNs with different architectures. Next, we focused on the attention mechanism of the RNN and, using tools from graph theory, we found that its graph properties allow one to hierarchically distinguish different architectures of the GRN. We show that the GRN responded differently to the addition of noise in the prediction by the RNN and we related the noise response to the analysis of the attention mechanism. In conclusion, this work provides a way to understand and exploit the attention mechanism of RNNs and it paves the way to RNN-based methods for time series prediction and inference of GRNs from gene expression data.
Collapse
Affiliation(s)
- Michele Monti
- RNA System Biology Lab, Department of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
- Centre for Genomic Regulation (CRG), The Barcelona Institute for Science and Technology, Dr. Aiguader 88, 08003 Barcelona, Spain
| | - Jonathan Fiorentino
- Center for Life Nano- & Neuro-Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy; (J.F.); (E.M.); (G.G.)
| | - Edoardo Milanetti
- Center for Life Nano- & Neuro-Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy; (J.F.); (E.M.); (G.G.)
- Department of Physics, Sapienza University of Rome, 00185 Rome, Italy
| | - Giorgio Gosti
- Center for Life Nano- & Neuro-Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy; (J.F.); (E.M.); (G.G.)
- Department of Physics, Sapienza University of Rome, 00185 Rome, Italy
| | - Gian Gaetano Tartaglia
- RNA System Biology Lab, Department of Neuroscience and Brain Technologies, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
- Centre for Genomic Regulation (CRG), The Barcelona Institute for Science and Technology, Dr. Aiguader 88, 08003 Barcelona, Spain
- Center for Life Nano- & Neuro-Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy; (J.F.); (E.M.); (G.G.)
- Department of Biology and Biotechnology Charles Darwin, Sapienza University of Rome, 00185 Rome, Italy
| |
Collapse
|
11
|
Ziane R, Camasses A, Radman-Livaja M. The asymmetric distribution of RNA polymerase II and nucleosomes on replicated daughter genomes is caused by differences in replication timing between the lagging and the leading strand. Genome Res 2022; 32:337-356. [PMID: 35042724 PMCID: PMC8805712 DOI: 10.1101/gr.275387.121] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2021] [Accepted: 12/19/2021] [Indexed: 11/24/2022]
Abstract
Chromatin features are thought to have a role in the epigenetic transmission of transcription states from one cell generation to the next. It is unclear how chromatin structure survives disruptions caused by genomic replication or whether chromatin features are instructive of the transcription state of the underlying gene. We developed a method to monitor budding yeast replication, transcription, and chromatin maturation dynamics on each daughter genome in parallel, with which we identified clusters of secondary origins surrounding known origins. We found a difference in the timing of lagging and leading strand replication on the order of minutes at most yeast genes. We propose a model in which the majority of old histones and RNA polymerase II (RNAPII) bind to the gene copy that replicated first, while newly synthesized nucleosomes are assembled on the copy that replicated second. RNAPII enrichment then shifts to the sister copy that replicated second. The order of replication is largely determined by genic orientation: If transcription and replication are codirectional, the leading strand replicates first; if they are counterdirectional, the lagging strand replicates first. A mutation in the Mcm2 subunit of the replicative helicase Mcm2-7 that impairs Mcm2 interactions with histone H3 slows down replication forks but does not qualitatively change the asymmetry in nucleosome distribution observed in the WT. We propose that active transcription states are inherited simultaneously and independently of their underlying chromatin states through the recycling of the transcription machinery and old histones, respectively. Transcription thus actively contributes to the reestablishment of the active chromatin state.
Collapse
Affiliation(s)
- Rahima Ziane
- Institut de Génétique Moléculaire de Montpellier, UMR 5535, CNRS, 34293 Montpellier Cedex 5, France.,Université de Montpellier, 34090 Montpellier, France
| | - Alain Camasses
- Institut de Génétique Moléculaire de Montpellier, UMR 5535, CNRS, 34293 Montpellier Cedex 5, France.,Université de Montpellier, 34090 Montpellier, France
| | - Marta Radman-Livaja
- Institut de Génétique Moléculaire de Montpellier, UMR 5535, CNRS, 34293 Montpellier Cedex 5, France.,Université de Montpellier, 34090 Montpellier, France
| |
Collapse
|
12
|
Alkallas R, Najafabadi HS. Analysis of mRNA Dynamics Using RNA Sequencing Data. Methods Mol Biol 2022; 2515:129-150. [PMID: 35776350 DOI: 10.1007/978-1-0716-2409-8_9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The RNA abundance of each gene is determined by its rates of transcription and RNA decay. Biochemical experiments that measure these rates, including transcription inhibition and metabolic labelling, are challenging to perform and are largely limited to in vitro settings. Most transcriptomic studies have focused on analyzing changes in RNA abundances without attributing those changes to transcriptional or posttranscriptional regulation. Estimating differential transcription and decay rates of RNA molecules would enable the identification of regulatory factors, such as transcription factors, RNA binding proteins, and microRNAs, that govern large-scale shifts in RNA expression. Here, we describe a protocol for estimating differential stability of RNA molecules between conditions using standard RNA-sequencing data, without the need for transcription inhibition or metabolic labeling. We apply this protocol to in vivo RNA-seq data from individuals with Alzheimer's disease and demonstrate how estimates of differential stability can be leveraged to infer the regulatory factors underlying them.
Collapse
Affiliation(s)
- Rached Alkallas
- Department of Human Genetics, McGill University, Montréal, QC, Canada
- McGill Genome Centre, Montréal, QC, Canada
- Rosalind and Morris Goodman Cancer Research Centre, Montréal, QC, Canada
| | - Hamed S Najafabadi
- Department of Human Genetics, McGill University, Montréal, QC, Canada.
- McGill Genome Centre, Montréal, QC, Canada.
| |
Collapse
|
13
|
Krenning L, Sonneveld S, Tanenbaum M. Time-resolved single-cell sequencing identifies multiple waves of mRNA decay during the mitosis-to-G1 phase transition. eLife 2022; 11:71356. [PMID: 35103592 PMCID: PMC8806192 DOI: 10.7554/elife.71356] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Accepted: 01/17/2022] [Indexed: 01/20/2023] Open
Abstract
Accurate control of the cell cycle is critical for development and tissue homeostasis, and requires precisely timed expression of many genes. Cell cycle gene expression is regulated through transcriptional and translational control, as well as through regulated protein degradation. Here, we show that widespread and temporally controlled mRNA decay acts as an additional mechanism for gene expression regulation during the cell cycle in human cells. We find that two waves of mRNA decay occur sequentially during the mitosis-to-G1 phase transition, and we identify the deadenylase CNOT1 as a factor that contributes to mRNA decay during this cell cycle transition. Collectively, our data show that, akin to protein degradation, scheduled mRNA decay helps to reshape cell cycle gene expression as cells move from mitosis into G1 phase.
Collapse
Affiliation(s)
- Lenno Krenning
- Oncode Institute, Hubrecht Institute – KNAW and University Medical Center UtrechtUtrechtNetherlands
| | - Stijn Sonneveld
- Oncode Institute, Hubrecht Institute – KNAW and University Medical Center UtrechtUtrechtNetherlands
| | - Marvin Tanenbaum
- Oncode Institute, Hubrecht Institute – KNAW and University Medical Center UtrechtUtrechtNetherlands
| |
Collapse
|
14
|
Mahadevan S, Cornwell JA, Chami B, Kelly E, Zoellner H. Cell-Projection Pumping of Fibroblast Contents into Osteosarcoma SAOS-2 Cells Correlates with Increased SAOS-2 Proliferation and Migration, as well as Altered Morphology. Biomolecules 2021; 11:1875. [PMID: 34944519 PMCID: PMC8699393 DOI: 10.3390/biom11121875] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2021] [Revised: 12/06/2021] [Accepted: 12/11/2021] [Indexed: 01/16/2023] Open
Abstract
We earlier reported that cell-projection pumping transfers fibroblast contents to cancer cells and this alters the cancer cell phenotype. Here, we report on single-cell tracking of time lapse recordings from co-cultured fluorescent fibroblasts and SAOS-2 osteosarcoma cells, tracking 5201 cells across 7 experiments. The fluorescent lipophilic marker DiD was used to label fibroblast organelles and to trace the transfer of fibroblast cytoplasm into SAOS-2 cells. We related SAOS-2 phenotypic change to levels of fluorescence transfer from fibroblasts to SAOS-2 cells, as well as what we term 'compensated fluorescence', that numerically projects mother cell fluorescence post-mitosis into daughter cells. The comparison of absolute with compensated fluorescence allowed us to deduct if the phenotypic effects in mother SAOS-2 cells were inherited by their daughters. SAOS-2 receipt of fibroblast fluorescence correlated by Kendall's tau with cell-profile area and without evidence of persistence in daughter cells (median tau = 0.51, p < 0.016); negatively and weakly with cell circularity and with evidence of persistence (median tau = -0.19, p < 0.05); and very weakly with cell migration velocity and without evidence of persistence (median tau = 0.01, p < 0.016). In addition, mitotic SAOS-2 cells had higher rates of prior fluorescence uptake (median = 64.9 units/day) than non-dividing cells (median = 35.6 units/day, p < 0.016) and there was no evidence of persistence post-mitosis. We conclude that there was an appreciable impact of cell-projection pumping on cancer cell phenotype relevant to cancer histopathological diagnosis, clinical spread and growth, with most effects being 'reset' by cancer cell mitosis.
Collapse
Affiliation(s)
- Swarna Mahadevan
- The Cellular and Molecular Pathology Research Unit, Oral Pathology and Oral Medicine, School of Dentistry, Faculty of Medicine and Health, The University of Sydney, Westmead Hospital, Westmead, NSW 2145, Australia; (S.M.); (J.A.C.); (B.C.); (E.K.)
| | - James A Cornwell
- The Cellular and Molecular Pathology Research Unit, Oral Pathology and Oral Medicine, School of Dentistry, Faculty of Medicine and Health, The University of Sydney, Westmead Hospital, Westmead, NSW 2145, Australia; (S.M.); (J.A.C.); (B.C.); (E.K.)
- Laboratory of Cancer Biology and Genetics, Center for Cancer Research, National Cancer Institute, Bethesda, MD 20892, USA
| | - Belal Chami
- The Cellular and Molecular Pathology Research Unit, Oral Pathology and Oral Medicine, School of Dentistry, Faculty of Medicine and Health, The University of Sydney, Westmead Hospital, Westmead, NSW 2145, Australia; (S.M.); (J.A.C.); (B.C.); (E.K.)
- Molecular Biomedicine, Charles Perkins Centre, The University of Sydney, Sydney, NSW 2006, Australia
| | - Elizabeth Kelly
- The Cellular and Molecular Pathology Research Unit, Oral Pathology and Oral Medicine, School of Dentistry, Faculty of Medicine and Health, The University of Sydney, Westmead Hospital, Westmead, NSW 2145, Australia; (S.M.); (J.A.C.); (B.C.); (E.K.)
| | - Hans Zoellner
- The Cellular and Molecular Pathology Research Unit, Oral Pathology and Oral Medicine, School of Dentistry, Faculty of Medicine and Health, The University of Sydney, Westmead Hospital, Westmead, NSW 2145, Australia; (S.M.); (J.A.C.); (B.C.); (E.K.)
- Biomedical Engineering, Faculty of Engineering, The University of Sydney, Sydney, NSW 2006, Australia
- Graduate School of Biomedical Engineering, University of NSW, Sydney, NSW 2052, Australia
- Strongarch Pty Ltd., Sydney, NSW 2000, Australia
| |
Collapse
|
15
|
Mok CH, MacLeod JN. Kinetics of Gene Expression Changes in Equine Fetal Interzone and Anlagen Cells Over 14 Days of Induced Chondrogenesis. Front Vet Sci 2021; 8:722324. [PMID: 34434986 PMCID: PMC8380811 DOI: 10.3389/fvets.2021.722324] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2021] [Accepted: 07/13/2021] [Indexed: 11/13/2022] Open
Abstract
Within developing synovial joints, interzone and anlagen cells progress through divergent chondrogenic pathways to generate stable articular cartilage and transient hypertrophic anlagen cartilage, respectively. Understanding the comparative cell biology between interzone and anlagen cells may provide novel insights into emergent cell-based therapies to support articular cartilage regeneration. The aim of this study was to assess the kinetics of gene expression profiles in these skeletal cell lines after inducing chondrogenesis in culture. Interzone and anlagen cells from seven equine fetuses were isolated and grown in a TGF-β1 chondrogenic inductive medium. Total RNA was isolated at ten time points (0, 1.5, 3, 6, 12, 24, 48, 96, 168, and 336 h), and gene expression for 93 targeted gene loci was measured in a microfluidic RT-qPCR system. Differential transcriptional responses were observed as early as 1.5 h after the initiation of chondrogenesis. Genes with functional annotations that include transcription regulation responded to the chondrogenic stimulation earlier (1.5–96 h) than genes involved in signal transduction (1.5–336 h) and the extracellular matrix biology (3–336 h). Between interzone and anlagen cell cultures, expression levels of 73 out of the 93 targeted genes were not initially different at 0 h, but 47 out of the 73 genes became differentially expressed under the chondrogenic stimulation. While interzone and anlagen cells are both chondrogenic, they display clear differences in response to the same TGF-β1 chondrogenic stimulation. This study provides new molecular insight into a timed sequence of the divergent developmental fates of interzone and anlagen cells in culture over 14 days.
Collapse
Affiliation(s)
- Chan Hee Mok
- Department of Veterinary Science, Gluck Equine Research Center, University of Kentucky, Lexington, KY, United States
| | - James N MacLeod
- Department of Veterinary Science, Gluck Equine Research Center, University of Kentucky, Lexington, KY, United States
| |
Collapse
|
16
|
Abstract
The protein-coding regions of mRNAs have the information to make proteins and hence have been at the center of attention for understanding altered protein functions in disease states, including cancer. Indeed, the discovery of genomic alterations and driver mutations that change protein levels and/or activity has been pivotal in our understanding of cancer biology. However, to better understand complex molecular mechanisms that are deregulated in cancers, we also need to look at non-coding parts of mRNAs, including 3'UTRs (untranslated regions), which control mRNA stability, localization, and translation efficiency. Recently, these rather overlooked regions of mRNAs are gaining attention as mounting evidence provides functional links between 3'UTRs, protein functions, and cancer-related molecular mechanisms. Here, roles of 3'UTRs in cancer biology and mechanisms that result in cancer-specific 3'-end isoform variants will be reviewed. An increased appreciation of 3'UTRs may help the discovery of new ways to explain as of yet unknown oncogene activation and tumor suppressor inactivation cases in cancers, and provide new avenues for diagnostic and therapeutic applications.
Collapse
Affiliation(s)
- Ayse Elif Erson-Bensan
- Department of Biological Sciences and Cancer Systems Biology Laboratory, Middle East Technical University (METU, ODTU), Dumlupinar Blv No: 1, Universiteler Mah, 06800, Ankara, Turkey.
| |
Collapse
|
17
|
Abstract
Circadian rhythms are constituted by a complex dynamical system with intertwined feedback loops, molecular switches, and self-sustained oscillations. Mathematical modeling supports understanding available heterogeneous kinetic data, highlights basic mechanisms, and can guide experimental research. Here, we introduce the basic steps from a biological question to simple models providing insight into gene-regulatory mechanisms. We illustrate the general approach by three examples: modeling decay processes, clock-controlled genes, and self-sustained oscillations.
Collapse
Affiliation(s)
- J Patrick Pett
- Institute for Theoretical Biology, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Pål O Westermark
- Leibniz Institute for Farm Animal Biology, Institute of Genetics and Biometry, Dummerstorf, Germany
| | - Hanspeter Herzel
- Institute for Theoretical Biology, Charité-Universitätsmedizin Berlin, Berlin, Germany.
| |
Collapse
|
18
|
Abstract
The identification and characterization of rhythmically expressed mRNAs have been an active area of research over the past 20 years, as these mRNAs are believed to produce the daily rhythms in a wide range of biological processes. Circadian transcriptome studies have used mature mRNA as a primary readout and focused largely on rhythmic RNA synthesis as a regulatory mechanism underlying rhythmic mRNA expression. However, RNA synthesis, RNA degradation, or a combination of both must be rhythmic to drive rhythmic RNA profiles, and it is still unclear to what extent rhythmic synthesis leads to rhythmic RNA profiles. In addition, circadian RNA expression is also often tissue specific. Although a handful of genes cycle in all or most tissues, others are rhythmic only in certain tissues, even though the same core clock mechanism is believed to control the rhythmic RNA profiles in all tissues. This review focuses on the dynamics of rhythmic RNA synthesis and degradation and discusses how these steps collectively determine the rhythmicity, phase, and amplitude of RNA accumulation. In particular, we highlight a possible role of RNA degradation in driving tissue-specific RNA rhythms. By unifying findings from experimental and theoretical studies, we will provide a comprehensive overview of how rhythmic gene expression can be achieved and how each regulatory step contributes to tissue-specific circadian transcriptome output in mammals.
Collapse
Affiliation(s)
| | - Shihoko Kojima
- To whom all correspondence should be addressed: Shihoko Kojima, Department of Biological Sciences, Fralin Life Sciences Institute, Virginia Tech, 1015 Life Science Circle, Blacksburg, VA, 24061, USA; .
| |
Collapse
|
19
|
de Pretis S, Furlan M, Pelizzola M. INSPEcT-GUI Reveals the Impact of the Kinetic Rates of RNA Synthesis, Processing, and Degradation, on Premature and Mature RNA Species. Front Genet 2020; 11:759. [PMID: 32765590 PMCID: PMC7379887 DOI: 10.3389/fgene.2020.00759] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2020] [Accepted: 06/26/2020] [Indexed: 12/23/2022] Open
Abstract
The abundance of RNA species and their response to perturbations are set by the kinetics rates of RNA synthesis, processing, and degradation. However, the visualization, interpretation, and manipulation of these data require familiarity with mathematical modeling and command line tools. INSPEcT-GUI is an R-Shiny interface that allows researchers without specific training to effortlessly explore how the fine kinetic regulation of the RNA life cycle can shape gene expression programs. In particular, it allows to: (i) interactively visualize gene-level RNA dynamics; (ii) refine the model fit of experimental data; (iii) test alternative regulatory models; (iv) explore, independently from the availability of data, how the combined action of the RNA kinetic rates impacts on premature and mature RNA. INSPEcT-GUI is freely available within the R/Bioconductor package INSPEcT at http://bioconductor.org/packages/INSPEcT/. An HTML vignette including documentation on the tool startup and usage, executable examples, and a video demonstration, are available at: http://bioconductor.org/packages/release/bioc/vignettes/INSPEcT/inst/doc/INSPEcT_GUI.html.
Collapse
Affiliation(s)
- Stefano de Pretis
- Center for Genomic Science of IIT@SEMM, Fondazione Istituto Italiano di Tecnologia, Milan, Italy
| | - Mattia Furlan
- Center for Genomic Science of IIT@SEMM, Fondazione Istituto Italiano di Tecnologia, Milan, Italy
| | - Mattia Pelizzola
- Center for Genomic Science of IIT@SEMM, Fondazione Istituto Italiano di Tecnologia, Milan, Italy
| |
Collapse
|
20
|
Fei Q, Zou Z, Roundtree IA, Sun HL, He C. YTHDF2 promotes mitotic entry and is regulated by cell cycle mediators. PLoS Biol 2020; 18:e3000664. [PMID: 32267835 PMCID: PMC7170294 DOI: 10.1371/journal.pbio.3000664] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2019] [Revised: 04/20/2020] [Accepted: 03/20/2020] [Indexed: 11/19/2022] Open
Abstract
The N6-methyladenosine (m6A) modification regulates mRNA stability and translation. Here, we show that transcriptomic m6A modification can be dynamic and the m6A reader protein YTH N6-methyladenosine RNA binding protein 2 (YTHDF2) promotes mRNA decay during cell cycle. Depletion of YTHDF2 in HeLa cells leads to the delay of mitotic entry due to overaccumulation of negative regulators of cell cycle such as Wee1-like protein kinase (WEE1). We demonstrate that WEE1 transcripts contain m6A modification, which promotes their decay through YTHDF2. Moreover, we found that YTHDF2 protein stability is dependent on cyclin-dependent kinase 1 (CDK1) activity. Thus, CDK1, YTHDF2, and WEE1 form a feedforward regulatory loop to promote mitotic entry. We further identified Cullin 1 (CUL1), Cullin 4A (CUL4A), damaged DNA-binding protein 1 (DDB1), and S-phase kinase-associated protein 2 (SKP2) as components of E3 ubiquitin ligase complexes that mediate YTHDF2 proteolysis. Our study provides insights into how cell cycle mediators modulate transcriptomic m6A modification, which in turn regulates the cell cycle.
Collapse
Affiliation(s)
- Qili Fei
- Department of Chemistry, The University of Chicago, Chicago, Illinois, United States of America
- Howard Hughes Medical Institute, The University of Chicago, Chicago, Illinois, United States of America
- Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, Guangdong, China
| | - Zhongyu Zou
- Department of Chemistry, The University of Chicago, Chicago, Illinois, United States of America
- Howard Hughes Medical Institute, The University of Chicago, Chicago, Illinois, United States of America
| | - Ian A. Roundtree
- Department of Chemistry, The University of Chicago, Chicago, Illinois, United States of America
- Howard Hughes Medical Institute, The University of Chicago, Chicago, Illinois, United States of America
- Medical Scientist Training Program, The University of Chicago, Chicago, Illinois, United States of America
- Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois, United States of America
- Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois, United States of America
| | - Hui-Lung Sun
- Department of Chemistry, The University of Chicago, Chicago, Illinois, United States of America
- Howard Hughes Medical Institute, The University of Chicago, Chicago, Illinois, United States of America
| | - Chuan He
- Department of Chemistry, The University of Chicago, Chicago, Illinois, United States of America
- Howard Hughes Medical Institute, The University of Chicago, Chicago, Illinois, United States of America
- Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois, United States of America
- Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois, United States of America
- * E-mail:
| |
Collapse
|
21
|
Caglar HO, Biray Avci C. Alterations of cell cycle genes in cancer: unmasking the role of cancer stem cells. Mol Biol Rep 2020; 47:3065-3076. [DOI: 10.1007/s11033-020-05341-6] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2019] [Accepted: 02/22/2020] [Indexed: 02/07/2023]
|
22
|
De los Santos H, Collins EJ, Mann C, Sagan AW, Jankowski MS, Bennett KP, Hurley JM. ECHO: an application for detection and analysis of oscillators identifies metabolic regulation on genome-wide circadian output. Bioinformatics 2020; 36:773-781. [PMID: 31384918 PMCID: PMC7523678 DOI: 10.1093/bioinformatics/btz617] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2018] [Revised: 07/24/2019] [Accepted: 08/02/2019] [Indexed: 01/07/2023] Open
Abstract
MOTIVATION Time courses utilizing genome scale data are a common approach to identifying the biological pathways that are controlled by the circadian clock, an important regulator of organismal fitness. However, the methods used to detect circadian oscillations in these datasets are not able to accommodate changes in the amplitude of the oscillations over time, leading to an underestimation of the impact of the clock on biological systems. RESULTS We have created a program to efficaciously identify oscillations in large-scale datasets, called the Extended Circadian Harmonic Oscillator application, or ECHO. ECHO utilizes an extended solution of the fixed amplitude oscillator that incorporates the amplitude change coefficient. Employing synthetic datasets, we determined that ECHO outperforms existing methods in detecting rhythms with decreasing oscillation amplitudes and in recovering phase shift. Rhythms with changing amplitudes identified from published biological datasets revealed distinct functions from those oscillations that were harmonic, suggesting purposeful biologic regulation to create this subtype of circadian rhythms. AVAILABILITY AND IMPLEMENTATION ECHO's full interface is available at https://github.com/delosh653/ECHO. An R package for this functionality, echo.find, can be downloaded at https://CRAN.R-project.org/package=echo.find. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Hannah De los Santos
- Department of Computer Science, Troy, NY 12180, USA,Institute for Data Exploration and Applications, Troy, NY 12180, USA
| | | | | | - April W Sagan
- Department of Mathematical Sciences, Troy, NY 12180, USA
| | | | - Kristin P Bennett
- Department of Computer Science, Troy, NY 12180, USA,Institute for Data Exploration and Applications, Troy, NY 12180, USA,Department of Mathematical Sciences, Troy, NY 12180, USA
| | - Jennifer M Hurley
- Department of Biological Sciences, Troy, NY 12180, USA,Center for Biotechnology and Interdisciplinary Sciences, Rensselaer Polytechnic Institute, Troy, NY 12180, USA,To whom correspondence should be addressed.
| |
Collapse
|
23
|
Catala M, Abou Elela S. Promoter-dependent nuclear RNA degradation ensures cell cycle-specific gene expression. Commun Biol 2019; 2:211. [PMID: 31240249 PMCID: PMC6572803 DOI: 10.1038/s42003-019-0441-3] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2018] [Accepted: 04/29/2019] [Indexed: 11/09/2022] Open
Abstract
Cell cycle progression depends on phase-specific gene expression. Here we show that the nuclear RNA degradation machinery plays a lead role in promoting cell cycle-dependent gene expression by triggering promoter-dependent co-transcriptional RNA degradation. Single molecule quantification of RNA abundance in different phases of the cell cycle indicates that relative curtailment of gene expression in certain phases is attained even when transcription is not completely inhibited. When nuclear ribonucleases are deleted, transcription of the Saccharomyces cerevisiae G1-specific axial budding gene AXL2 is detected throughout the cell cycle and its phase-specific expression is lost. Promoter replacement abolished cell cycle-dependent RNA degradation and rendered the RNA insensitive to the deletion of nuclear ribonucleases. Together the data reveal a model of gene regulation whereby RNA abundance is controlled by promoter-dependent induction of RNA degradation.
Collapse
Affiliation(s)
- Mathieu Catala
- Département de microbiologie et d’infectiologie, Faculté de médecine et des sciences de la santé, Université de Sherbrooke, Sherbrooke, QC J1E 4K8 Canada
| | - Sherif Abou Elela
- Département de microbiologie et d’infectiologie, Faculté de médecine et des sciences de la santé, Université de Sherbrooke, Sherbrooke, QC J1E 4K8 Canada
| |
Collapse
|
24
|
Dronamraju R, Hepperla AJ, Shibata Y, Adams AT, Magnuson T, Davis IJ, Strahl BD. Spt6 Association with RNA Polymerase II Directs mRNA Turnover During Transcription. Mol Cell 2019; 70:1054-1066.e4. [PMID: 29932900 DOI: 10.1016/j.molcel.2018.05.020] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2017] [Revised: 03/23/2018] [Accepted: 05/17/2018] [Indexed: 10/28/2022]
Abstract
Spt6 is an essential histone chaperone that mediates nucleosome reassembly during gene transcription. Spt6 also associates with RNA polymerase II (RNAPII) via a tandem Src2 homology domain. However, the significance of Spt6-RNAPII interaction is not well understood. Here, we show that Spt6 recruitment to genes and the nucleosome reassembly functions of Spt6 can still occur in the absence of its association with RNAPII. Surprisingly, we found that Spt6-RNAPII association is required for efficient recruitment of the Ccr4-Not de-adenylation complex to transcribed genes for essential degradation of a range of mRNAs, including mRNAs required for cell-cycle progression. These findings reveal an unexpected control mechanism for mRNA turnover during transcription facilitated by a histone chaperone.
Collapse
Affiliation(s)
- Raghuvar Dronamraju
- Department of Biochemistry & Biophysics, University of North Carolina School of Medicine, Chapel Hill, NC 27599, USA; Lineberger Comprehensive Cancer Center, University of North Carolina School of Medicine, Chapel Hill, NC 27599, USA
| | - Austin J Hepperla
- Curriculum in Genetics and Molecular Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Yoichiro Shibata
- Lineberger Comprehensive Cancer Center, University of North Carolina School of Medicine, Chapel Hill, NC 27599, USA; Curriculum in Genetics and Molecular Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Alexander T Adams
- Department of Biochemistry & Biophysics, University of North Carolina School of Medicine, Chapel Hill, NC 27599, USA
| | - Terry Magnuson
- Lineberger Comprehensive Cancer Center, University of North Carolina School of Medicine, Chapel Hill, NC 27599, USA; Curriculum in Genetics and Molecular Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA; Department of Genetics, The Carolina Center for Genome Sciences, University of North Carolina, Chapel Hill, NC 27599, USA
| | - Ian J Davis
- Lineberger Comprehensive Cancer Center, University of North Carolina School of Medicine, Chapel Hill, NC 27599, USA; Curriculum in Genetics and Molecular Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA; Department of Genetics, The Carolina Center for Genome Sciences, University of North Carolina, Chapel Hill, NC 27599, USA; Departments of Pediatrics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Brian D Strahl
- Department of Biochemistry & Biophysics, University of North Carolina School of Medicine, Chapel Hill, NC 27599, USA; Lineberger Comprehensive Cancer Center, University of North Carolina School of Medicine, Chapel Hill, NC 27599, USA; Curriculum in Genetics and Molecular Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA.
| |
Collapse
|
25
|
A transcriptome-wide analysis deciphers distinct roles of G1 cyclins in temporal organization of the yeast cell cycle. Sci Rep 2019; 9:3343. [PMID: 30833602 PMCID: PMC6399449 DOI: 10.1038/s41598-019-39850-7] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2018] [Accepted: 01/30/2019] [Indexed: 12/11/2022] Open
Abstract
Oscillating gene expression is crucial for correct timing and progression through cell cycle. In Saccharomyces cerevisiae, G1 cyclins Cln1-3 are essential drivers of the cell cycle and have an important role for temporal fine-tuning. We measured time-resolved transcriptome-wide gene expression for wild type and cyclin single and double knockouts over cell cycle with and without osmotic stress. Clustering of expression profiles, peak time detection of oscillating genes, integration with transcription factor network dynamics, and assignment to cell cycle phases allowed us to quantify the effect of genetic or stress perturbations on the duration of cell cycle phases. Cln1 and Cln2 showed functional differences, especially affecting later phases. Deletion of Cln3 led to a delay of START followed by normal progression through later phases. Our data and network analysis suggest mutual effects of cyclins with the transcriptional regulators SBF and MBF.
Collapse
|
26
|
Hausser J, Mayo A, Keren L, Alon U. Central dogma rates and the trade-off between precision and economy in gene expression. Nat Commun 2019; 10:68. [PMID: 30622246 PMCID: PMC6325141 DOI: 10.1038/s41467-018-07391-8] [Citation(s) in RCA: 106] [Impact Index Per Article: 21.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2018] [Accepted: 10/18/2018] [Indexed: 12/31/2022] Open
Abstract
Steady-state protein abundance is set by four rates: transcription, translation, mRNA decay and protein decay. A given protein abundance can be obtained from infinitely many combinations of these rates. This raises the question of whether the natural rates for each gene result from historical accidents, or are there rules that give certain combinations a selective advantage? We address this question using high-throughput measurements in rapidly growing cells from diverse organisms to find that about half of the rate combinations do not exist: genes that combine high transcription with low translation are strongly depleted. This depletion is due to a trade-off between precision and economy: high transcription decreases stochastic fluctuations but increases transcription costs. Our theory quantitatively explains which rate combinations are missing, and predicts the curvature of the fitness function for each gene. It may guide the design of gene circuits with desired expression levels and noise.
Collapse
Affiliation(s)
- Jean Hausser
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, 76100, Israel.
| | - Avi Mayo
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, 76100, Israel
| | - Leeat Keren
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, 76100, Israel
| | - Uri Alon
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, 76100, Israel.
| |
Collapse
|
27
|
Yamada T, Akimitsu N. Contributions of regulated transcription and mRNA decay to the dynamics of gene expression. WILEY INTERDISCIPLINARY REVIEWS-RNA 2018; 10:e1508. [PMID: 30276972 DOI: 10.1002/wrna.1508] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/29/2018] [Revised: 08/06/2018] [Accepted: 08/27/2018] [Indexed: 12/21/2022]
Abstract
Organisms have acquired sophisticated regulatory networks that control gene expression in response to cellular perturbations. Understanding of the mechanisms underlying the coordinated changes in gene expression in response to external and internal stimuli is a fundamental issue in biology. Recent advances in high-throughput technologies have enabled the measurement of diverse biological information, including gene expression levels, kinetics of gene expression, and interactions among gene expression regulatory molecules. By coupling these technologies with quantitative modeling, we can now uncover the biological roles and mechanisms of gene regulation at the system level. This review consists of two parts. First, we focus on the methods using uridine analogs that measure synthesis and decay rates of RNAs, which demonstrate how cells dynamically change the regulation of gene expression in response to both internal and external cues. Second, we discuss the underlying mechanisms of these changes in kinetics, including the functions of transcription factors and RNA-binding proteins. Overall, this review will help to clarify a system-level view of gene expression programs in cells. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Turnover and Surveillance > Regulation of RNA Stability RNA Methods > RNA Analyses in vitro and In Silico.
Collapse
Affiliation(s)
- Toshimichi Yamada
- Department of Molecular and Cellular Biochemistry, Meiji Pharmaceutical University, Tokyo, Japan
| | | |
Collapse
|
28
|
Farina L, Paci P. A feature-based integrated scoring scheme for cell cycle-regulated genes prioritization. J Theor Biol 2018; 459:130-141. [PMID: 30261169 DOI: 10.1016/j.jtbi.2018.09.025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2018] [Revised: 08/03/2018] [Accepted: 09/23/2018] [Indexed: 10/28/2022]
Abstract
Prioritization of cell cycle-regulated genes from expression time-profiles is still an open problem. The point at issue is the surprisingly poor overlap among ranked lists obtained from different experimental protocols. Instead of developing a general-purpose computational methodology for detecting periodic signals, we focus on the budding yeast mitotic cell cycle. The reason being that the current availability of a total of 12 datasets, produced by 6 independent groups using 4 different synchronization methods, permits a re-analysis and re-consideration of this problem in a more reliable and extensive data domain. Notably, budding yeast is a model organism for studying cancer and testing new drugs. Here we propose a novel multi-feature score (called PERLA, PERiodicity, Regulation and Lag-Autocorrelation) that integrates different features of cell cycle-regulated gene expression time-profiles. We obtained increased performances on a wide range of benchmarks and, most importantly, a substantially increased overlap of the top ranking genes among different datasets, thus proving the effectiveness of the proposed prioritization algorithm. Examples on how to use PERLA to gain new insight into the biology of the cell cycle, are provided in a final dedicated section.
Collapse
Affiliation(s)
- Lorenzo Farina
- Department of Computer, Control and Management Engineering "A. Ruberti", Sapienza University of Rome, Italy; Institute for Systems Analysis and Computer Science "A. Ruberti", National Research Council, Rome, Italy.
| | - Paola Paci
- Institute for Systems Analysis and Computer Science "A. Ruberti", National Research Council, Rome, Italy; SysBio Centre for Systems Biology, Rome, Italy.
| |
Collapse
|
29
|
Wittmann G, Lechan RM. Prss56 expression in the rodent hypothalamus: Inverse correlation with pro-opiomelanocortin suggests oscillatory gene expression in adult rat tanycytes. J Comp Neurol 2018; 526:2444-2461. [PMID: 30242838 DOI: 10.1002/cne.24504] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2018] [Revised: 06/27/2018] [Accepted: 07/12/2018] [Indexed: 11/11/2022]
Abstract
We recently reported that the number of hypothalamic tanycytes expressing pro-opiomelanocortin (Pomc) is highly variable among brains of adult rats. While its cause and significance remain unknown, identifying other variably expressed genes in tanycytes may help understand this curious phenomenon. In this in situ hybridization study, we report that the Prss56 gene, which encodes a trypsin-like serine protease and is expressed in neural stem/progenitor cells, shows a similarly variable mRNA expression in tanycytes of adult rats and correlates inversely with tanycyte Pomc mRNA. Prss56 was expressed in α1, β1, subsets of α2, and some median eminence γ tanycytes, but virtually absent from β2 tanycytes. Prss56 was also expressed in vimentin positive tanycyte-like cells in the parenchyma of the ventromedial and arcuate nuclei, and in thyrotropin beta subunit-expressing cells of the pars tuberalis of the pituitary. In contrast to adults, Prss56 expression was uniformly high in tanycytes in adolescent rats. In mice, Prss56-expressing tanycytes and parenchymal cells were also observed but fewer in number and without significant variations. The results identify Prss56 as a second gene that is expressed variably in tanycytes of adult rats. We propose that the variable, inversely correlating expression of Prss56 and Pomc reflect periodically oscillating gene expression in tanycytes rather than stable expression levels that vary between individual rats. A possible functional link between Prss56 and POMC, and Prss56 as a potential marker for migrating tanycytes are discussed.
Collapse
Affiliation(s)
- Gábor Wittmann
- Department of Medicine, Division of Endocrinology, Diabetes and Metabolism, Tupper Research Institute, Tufts Medical Center, Boston, Massachusetts
| | - Ronald M Lechan
- Department of Medicine, Division of Endocrinology, Diabetes and Metabolism, Tupper Research Institute, Tufts Medical Center, Boston, Massachusetts.,Department of Neuroscience, Tufts University School of Medicine, Boston, Massachusetts
| |
Collapse
|
30
|
Arani BMS, Mahmoudi M, Lahti L, González J, Wit EC. Stability estimation of autoregulated genes under Michaelis-Menten-type kinetics. Phys Rev E 2018; 97:062407. [PMID: 30011543 DOI: 10.1103/physreve.97.062407] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2017] [Indexed: 11/07/2022]
Abstract
Feedback loops are typical motifs appearing in gene regulatory networks. In some well-studied model organisms, including Escherichia coli, autoregulated genes, i.e., genes that activate or repress themselves through their protein products, are the only feedback interactions. For these types of interactions, the Michaelis-Menten (MM) formulation is a suitable and widely used approach, which always leads to stable steady-state solutions representative of homeostatic regulation. However, in many other biological phenomena, such as cell differentiation, cancer progression, and catastrophes in ecosystems, one might expect to observe bistable switchlike dynamics in the case of strong positive autoregulation. To capture this complex behavior we use the generalized family of MM kinetic models. We give a full analysis regarding the stability of autoregulated genes. We show that the autoregulation mechanism has the capability to exhibit diverse cellular dynamics including hysteresis, a typical characteristic of bistable systems, as well as irreversible transitions between bistable states. We also introduce a statistical framework to estimate the kinetics parameters and probability of different stability regimes given observational data. Empirical data for the autoregulated gene SCO3217 in the SOS system in Streptomyces coelicolor are analyzed. The coupling of a statistical framework and the mathematical model can give further insight into understanding the evolutionary mechanisms toward different cell fates in various systems.
Collapse
Affiliation(s)
- Babak M S Arani
- Department of Aquatic Ecology and Water Quality Management, Wageningen University, P.O. Box 47, NL-6700 AA Wageningen, The Netherlands
| | - Mahdi Mahmoudi
- Faculty of Mathematics, Statistics and Computer Science, Semnan University, P.O. Box 35195-363, Semnan, Iran
| | - Leo Lahti
- Department of Mathematics and Statistics, University of Turku, FI-20014 Turku, Finland
| | - Javier González
- Department of Mathematics and Statistics, Lancaster University, Lancaster LA1 4YF, United Kingdom and Amazon Research Cambridge, Cambridge, United Kingdom
| | - Ernst C Wit
- Institute of Computational Science, USI, Via G. Buffi 13, Lugano 6900, Switzerland
| |
Collapse
|
31
|
Transcriptional timing and noise of yeast cell cycle regulators-a single cell and single molecule approach. NPJ Syst Biol Appl 2018; 4:17. [PMID: 29844922 PMCID: PMC5962571 DOI: 10.1038/s41540-018-0053-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2017] [Revised: 04/05/2018] [Accepted: 04/19/2018] [Indexed: 12/12/2022] Open
Abstract
Gene expression is a stochastic process and its appropriate regulation is critical for cell cycle progression. Cellular stress response necessitates expression reprogramming and cell cycle arrest. While previous studies are mostly based on bulk experiments influenced by synchronization effects or lack temporal distribution, time-resolved methods on single cells are needed to understand eukaryotic cell cycle in context of noisy gene expression and external perturbations. Using smFISH, microscopy and morphological markers, we monitored mRNA abundances over cell cycle phases and calculated transcriptional noise for SIC1, CLN2, and CLB5, the main G1/S transition regulators in budding yeast. We employed mathematical modeling for in silico synchronization and for derivation of time-courses from single cell data. This approach disclosed detailed quantitative insights into transcriptional regulation with and without stress, not available from bulk experiments before. First, besides the main peak in G1 we found an upshift of CLN2 and CLB5 expression in late mitosis. Second, all three genes showed basal expression throughout cell cycle enlightening that transcription is not divided in on and off but rather in high and low phases. Finally, exposing cells to osmotic stress revealed different periods of transcriptional inhibition for CLN2 and CLB5 and the impact of stress on cell cycle phase duration. Combining experimental and computational approaches allowed us to precisely assess cell cycle progression timing, as well as gene expression dynamics.
Collapse
|
32
|
A Multi-Parameter Analysis of Cellular Coordination of Major Transcriptome Regulation Mechanisms. Sci Rep 2018; 8:5742. [PMID: 29636505 PMCID: PMC5893539 DOI: 10.1038/s41598-018-24039-1] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2017] [Accepted: 03/21/2018] [Indexed: 01/06/2023] Open
Abstract
To understand cellular coordination of multiple transcriptome regulation mechanisms, we simultaneously measured transcription rate (TR), mRNA abundance (RA) and translation activity (TA). This revealed multiple insights. First, the three parameters displayed systematic statistical differences. Sequentially more genes exhibited extreme (low or high) expression values from TR to RA, and then to TA; that is, cellular coordination of multiple transcriptome regulatory mechanisms leads to sequentially enhanced gene expression selectivity as the genetic information flow from the genome to the proteome. Second, contribution of the stabilization-by-translation regulatory mechanism to the cellular coordination process was assessed. The data enabled an estimation of mRNA stability, revealing a moderate but significant positive correlation between mRNA stability and translation activity. Third, the proportion of mRNA occupied by un-translated regions (UTR) exhibited a negative relationship with the level of this correlation, and was thus a major determinant of the mode of regulation of the mRNA. High-UTR-proportion mRNAs tend to defy the stabilization-by-translation regulatory mechanism, staying out of the polysome but remaining stable; mRNAs with little UTRs largely followed this regulation. In summary, we quantitatively delineated the relationship among multiple transcriptome regulation parameters, i.e., cellular coordination of corresponding regulatory mechanisms.
Collapse
|
33
|
Genome-Wide Mapping of Decay Factor-mRNA Interactions in Yeast Identifies Nutrient-Responsive Transcripts as Targets of the Deadenylase Ccr4. G3-GENES GENOMES GENETICS 2018; 8:315-330. [PMID: 29158339 PMCID: PMC5765359 DOI: 10.1534/g3.117.300415] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
The Ccr4 (carbon catabolite repression 4)-Not complex is a major regulator of stress responses that controls gene expression at multiple levels, from transcription to mRNA decay. Ccr4, a “core” subunit of the complex, is the main cytoplasmic deadenylase in Saccharomyces cerevisiae; however, its mRNA targets have not been mapped on a genome-wide scale. Here, we describe a genome-wide approach, RNA immunoprecipitation (RIP) high-throughput sequencing (RIP-seq), to identify the RNAs bound to Ccr4, and two proteins that associate with it, Dhh1 and Puf5. All three proteins were preferentially bound to lowly abundant mRNAs, most often at the 3′ end of the transcript. Furthermore, Ccr4, Dhh1, and Puf5 are recruited to mRNAs that are targeted by other RNA-binding proteins that promote decay and mRNA transport, and inhibit translation. Although Ccr4-Not regulates mRNA transcription and decay, Ccr4 recruitment to mRNAs correlates better with decay rates, suggesting it imparts greater control over transcript abundance through decay. Ccr4-enriched mRNAs are refractory to control by the other deadenylase complex in yeast, Pan2/3, suggesting a division of labor between these deadenylation complexes. Finally, Ccr4 and Dhh1 associate with mRNAs whose abundance increases during nutrient starvation, and those that fluctuate during metabolic and oxygen consumption cycles, which explains the known genetic connections between these factors and nutrient utilization and stress pathways.
Collapse
|
34
|
Mena A, Medina DA, García-Martínez J, Begley V, Singh A, Chávez S, Muñoz-Centeno MC, Pérez-Ortín JE. Asymmetric cell division requires specific mechanisms for adjusting global transcription. Nucleic Acids Res 2017; 45:12401-12412. [PMID: 29069448 PMCID: PMC5716168 DOI: 10.1093/nar/gkx974] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2017] [Accepted: 10/10/2017] [Indexed: 12/19/2022] Open
Abstract
Most cells divide symmetrically into two approximately identical cells. There are many examples, however, of asymmetric cell division that can generate sibling cell size differences. Whereas physical asymmetric division mechanisms and cell fate consequences have been investigated, the specific problem caused by asymmetric division at the transcription level has not yet been addressed. In symmetrically dividing cells the nascent transcription rate increases in parallel to cell volume to compensate it by keeping the actual mRNA synthesis rate constant. This cannot apply to the yeast Saccharomyces cerevisiae, where this mechanism would provoke a never-ending increasing mRNA synthesis rate in smaller daughter cells. We show here that, contrarily to other eukaryotes with symmetric division, budding yeast keeps the nascent transcription rates of its RNA polymerases constant and increases mRNA stability. This control on RNA pol II-dependent transcription rate is obtained by controlling the cellular concentration of this enzyme.
Collapse
Affiliation(s)
- Adriana Mena
- Departamento de Bioquímica y Biología Molecular and E.R.I. Biotecmed, Universitat de València, Dr. Moliner, 50, Burjassot 46100, Valencia, Spain
| | - Daniel A Medina
- Departamento de Bioquímica y Biología Molecular and E.R.I. Biotecmed, Universitat de València, Dr. Moliner, 50, Burjassot 46100, Valencia, Spain
| | - José García-Martínez
- Departamento de Genética and E.R.I. Biotecmed, Universitat de València, Dr. Moliner, 50, Burjassot 46100, Valencia, Spain
| | - Victoria Begley
- Departamento de Genética, Universidad de Sevilla and Instituto de Biomedicina de Sevilla (IBiS), Hospital Virgen del Rocío-CSIC-Universidad de Sevilla, 41013 Sevilla, Spain
| | - Abhyudai Singh
- Department of Electrical and Computer Engineering, University of Delaware, Newark, DE 19716, USA
| | - Sebastián Chávez
- Departamento de Genética, Universidad de Sevilla and Instituto de Biomedicina de Sevilla (IBiS), Hospital Virgen del Rocío-CSIC-Universidad de Sevilla, 41013 Sevilla, Spain
| | - Mari C Muñoz-Centeno
- Departamento de Genética, Universidad de Sevilla and Instituto de Biomedicina de Sevilla (IBiS), Hospital Virgen del Rocío-CSIC-Universidad de Sevilla, 41013 Sevilla, Spain
| | - José E Pérez-Ortín
- Departamento de Bioquímica y Biología Molecular and E.R.I. Biotecmed, Universitat de València, Dr. Moliner, 50, Burjassot 46100, Valencia, Spain
| |
Collapse
|
35
|
Chávez S, Eastman G, Smircich P, Becco LL, Oliveira-Rizzo C, Fort R, Potenza M, Garat B, Sotelo-Silveira JR, Duhagon MA. Transcriptome-wide analysis of the Trypanosoma cruzi proliferative cycle identifies the periodically expressed mRNAs and their multiple levels of control. PLoS One 2017; 12:e0188441. [PMID: 29182646 PMCID: PMC5705152 DOI: 10.1371/journal.pone.0188441] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2017] [Accepted: 11/07/2017] [Indexed: 12/02/2022] Open
Abstract
Trypanosoma cruzi is the protozoan parasite causing American trypanosomiasis or Chagas disease, a neglected parasitosis with important human health impact in Latin America. The efficacy of current therapy is limited, and its toxicity is high. Since parasite proliferation is a fundamental target for rational drug design, we sought to progress into its understanding by applying a genome-wide approach. Treating a TcI linage strain with hydroxyurea, we isolated epimastigotes in late G1, S and G2/M cell cycle stages at 70% purity. The sequencing of each phase identified 305 stage-specific transcripts (1.5-fold change, p≤0.01), coding for conserved cell cycle regulated proteins and numerous proteins whose cell cycle dependence has not been recognized before. Comparisons with the parasite T. brucei and the human host reveal important differences. The meta-analysis of T. cruzi transcriptomic and ribonomic data indicates that cell cycle regulated mRNAs are subject to sub-cellular compartmentalization. Compositional and structural biases of these genes- including CAI, GC content, UTR length, and polycistron position- may contribute to their regulation. To discover nucleotide motifs responsible for the co-regulation of cell cycle regulated genes, we looked for overrepresented motifs at their UTRs and found a variant of the cell cycle sequence motif at the 3' UTR of most of the S and G2 stage genes. We additionally identified hairpin structures at the 5' UTRs of a high proportion of the transcripts, suggesting that periodic gene expression might also rely on translation initiation in T. cruzi. In summary, we report a comprehensive list of T. cruzi cell cycle regulated genes, including many previously unstudied proteins, we show evidence favoring a multi-step control of their expression, and we identify mRNA motifs that may mediate their regulation. Our results provide novel information of the T. cruzi proliferative proteins and the integrated levels of their gene expression control.
Collapse
Affiliation(s)
- Santiago Chávez
- Laboratory of Molecular Interactions, School of Sciences, Universidad de la República, Montevideo, Uruguay
- Department of Genetics, School of Medicine, Universidad de la República, Montevideo, Uruguay
| | - Guillermo Eastman
- Department of Genomics, Instituto de Investigaciones Biológicas Clemente Estable, Montevideo, Uruguay
| | - Pablo Smircich
- Laboratory of Molecular Interactions, School of Sciences, Universidad de la República, Montevideo, Uruguay
- Department of Genetics, School of Medicine, Universidad de la República, Montevideo, Uruguay
| | - Lorena Lourdes Becco
- Laboratory of Molecular Interactions, School of Sciences, Universidad de la República, Montevideo, Uruguay
| | - Carolina Oliveira-Rizzo
- Laboratory of Molecular Interactions, School of Sciences, Universidad de la República, Montevideo, Uruguay
- Department of Genetics, School of Medicine, Universidad de la República, Montevideo, Uruguay
| | - Rafael Fort
- Laboratory of Molecular Interactions, School of Sciences, Universidad de la República, Montevideo, Uruguay
- Department of Genetics, School of Medicine, Universidad de la República, Montevideo, Uruguay
| | - Mariana Potenza
- Institute for Research in Genetic Engineering and Molecular Biology 'Dr. N.H. Torres', Buenos Aires, Argentina
| | - Beatriz Garat
- Laboratory of Molecular Interactions, School of Sciences, Universidad de la República, Montevideo, Uruguay
| | - José Roberto Sotelo-Silveira
- Department of Genomics, Instituto de Investigaciones Biológicas Clemente Estable, Montevideo, Uruguay
- Department of Cell and Molecular Biology, School of Sciences, Universidad de la República, Montevideo, Uruguay
| | - María Ana Duhagon
- Laboratory of Molecular Interactions, School of Sciences, Universidad de la República, Montevideo, Uruguay
- Department of Genetics, School of Medicine, Universidad de la República, Montevideo, Uruguay
- * E-mail:
| |
Collapse
|
36
|
Golan-Lavi R, Giacomelli C, Fuks G, Zeisel A, Sonntag J, Sinha S, Köstler W, Wiemann S, Korf U, Yarden Y, Domany E. Coordinated Pulses of mRNA and of Protein Translation or Degradation Produce EGF-Induced Protein Bursts. Cell Rep 2017; 18:3129-3142. [PMID: 28355565 DOI: 10.1016/j.celrep.2017.03.014] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2014] [Revised: 01/16/2017] [Accepted: 03/01/2017] [Indexed: 11/25/2022] Open
Abstract
Protein responses to extracellular cues are governed by gene transcription, mRNA degradation and translation, and protein degradation. In order to understand how these time-dependent processes cooperate to generate dynamic responses, we analyzed the response of human mammary cells to the epidermal growth factor (EGF). Integrating time-dependent transcript and protein data into a mathematical model, we inferred for several proteins their pre-and post-stimulus translation and degradation coefficients and found that they exhibit complex, time-dependent variation. Specifically, we identified strategies of protein production and degradation acting in concert to generate rapid, transient protein bursts in response to EGF. Remarkably, for some proteins, for which the response necessitates rapidly decreased abundance, cells exhibit a transient increase in the corresponding degradation coefficient. Our model and analysis allow inference of the kinetics of mRNA translation and protein degradation, without perturbing cells, and open a way to understanding the fundamental processes governing time-dependent protein abundance profiles.
Collapse
Affiliation(s)
- Roni Golan-Lavi
- Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 7610001, Israel; Department of Biological Regulation, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Chiara Giacomelli
- Division of Molecular Genome Analysis, German Cancer Research Center (DKFZ), Heidelberg 69120, Germany
| | - Garold Fuks
- Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Amit Zeisel
- Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Johanna Sonntag
- Division of Molecular Genome Analysis, German Cancer Research Center (DKFZ), Heidelberg 69120, Germany
| | - Sanchari Sinha
- Department of Biological Regulation, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Wolfgang Köstler
- Department of Biological Regulation, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Stefan Wiemann
- Division of Molecular Genome Analysis, German Cancer Research Center (DKFZ), Heidelberg 69120, Germany
| | - Ulrike Korf
- Division of Molecular Genome Analysis, German Cancer Research Center (DKFZ), Heidelberg 69120, Germany
| | - Yosef Yarden
- Department of Biological Regulation, Weizmann Institute of Science, Rehovot 7610001, Israel.
| | - Eytan Domany
- Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 7610001, Israel.
| |
Collapse
|
37
|
Szedlak A, Sims S, Smith N, Paternostro G, Piermarocchi C. Cell cycle time series gene expression data encoded as cyclic attractors in Hopfield systems. PLoS Comput Biol 2017; 13:e1005849. [PMID: 29149186 PMCID: PMC5711035 DOI: 10.1371/journal.pcbi.1005849] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2017] [Revised: 12/01/2017] [Accepted: 10/25/2017] [Indexed: 12/18/2022] Open
Abstract
Modern time series gene expression and other omics data sets have enabled unprecedented resolution of the dynamics of cellular processes such as cell cycle and response to pharmaceutical compounds. In anticipation of the proliferation of time series data sets in the near future, we use the Hopfield model, a recurrent neural network based on spin glasses, to model the dynamics of cell cycle in HeLa (human cervical cancer) and S. cerevisiae cells. We study some of the rich dynamical properties of these cyclic Hopfield systems, including the ability of populations of simulated cells to recreate experimental expression data and the effects of noise on the dynamics. Next, we use a genetic algorithm to identify sets of genes which, when selectively inhibited by local external fields representing gene silencing compounds such as kinase inhibitors, disrupt the encoded cell cycle. We find, for example, that inhibiting the set of four kinases AURKB, NEK1, TTK, and WEE1 causes simulated HeLa cells to accumulate in the M phase. Finally, we suggest possible improvements and extensions to our model. Cell cycle—the process in which a parent cell replicates its DNA and divides into two daughter cells—is an upregulated process in many forms of cancer. Identifying gene inhibition targets to regulate cell cycle is important to the development of effective therapies. Although modern high throughput techniques offer unprecedented resolution of the molecular details of biological processes like cell cycle, analyzing the vast quantities of the resulting experimental data and extracting actionable information remains a formidable task. Here, we create a dynamical model of the process of cell cycle using the Hopfield model (a type of recurrent neural network) and gene expression data from human cervical cancer cells and yeast cells. We find that the model recreates the oscillations observed in experimental data. Tuning the level of noise (representing the inherent randomness in gene expression and regulation) to the “edge of chaos” is crucial for the proper behavior of the system. We then use this model to identify potential gene targets for disrupting the process of cell cycle. This method could be applied to other time series data sets and used to predict the effects of untested targeted perturbations.
Collapse
Affiliation(s)
- Anthony Szedlak
- Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan, United States of America
| | - Spencer Sims
- Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan, United States of America
| | - Nicholas Smith
- Salgomed Inc., Del Mar, California, United States of America
| | - Giovanni Paternostro
- Sanford Burnham Prebys Medical Discovery Institute, La Jolla, California, United States of America
| | - Carlo Piermarocchi
- Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan, United States of America
- * E-mail:
| |
Collapse
|
38
|
An improved MS2 system for accurate reporting of the mRNA life cycle. Nat Methods 2017; 15:81-89. [PMID: 29131164 PMCID: PMC5843578 DOI: 10.1038/nmeth.4502] [Citation(s) in RCA: 206] [Impact Index Per Article: 29.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2017] [Accepted: 10/04/2017] [Indexed: 12/19/2022]
Abstract
The MS2-MCP system enables researchers to image multiple steps of the mRNA life cycle with high temporal and spatial resolution. However, for short-lived mRNAs, the tight binding of the MS2 coat protein (MCP) to the MS2 binding sites (MBS) protects the RNA from being efficiently degraded, and this confounds the study of mRNA regulation. Here, we describe a reporter system (MBSV6) with reduced affinity for the MCP, which allows mRNA degradation while preserving single-molecule detection determined by single-molecule FISH (smFISH) or live imaging. Constitutive mRNAs (MDN1 and DOA1) and highly-regulated mRNAs (GAL1 and ASH1) endogenously tagged with MBSV6 in Saccharomyces cerevisiae degrade normally. As a result, short-lived mRNAs were imaged throughout their complete life cycle. The MBSV6 reporter revealed that, in contrast to previous findings, coordinated recruitment of mRNAs at specialized structures such as P-bodies during stress did not occur, and mRNA degradation was heterogeneously distributed in the cytoplasm.
Collapse
|
39
|
Cho CY, Motta FC, Kelliher CM, Deckard A, Haase SB. Reconciling conflicting models for global control of cell-cycle transcription. Cell Cycle 2017; 16:1965-1978. [PMID: 28934013 PMCID: PMC5638368 DOI: 10.1080/15384101.2017.1367073] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2017] [Accepted: 08/07/2017] [Indexed: 10/18/2022] Open
Abstract
Models for the control of global cell-cycle transcription have advanced from a CDK-APC/C oscillator, a transcription factor (TF) network, to coupled CDK-APC/C and TF networks. Nonetheless, current models were challenged by a recent study that concluded that the cell-cycle transcriptional program is primarily controlled by a CDK-APC/C oscillator in budding yeast. Here we report an analysis of the transcriptome dynamics in cyclin mutant cells that were not queried in the previous study. We find that B-cyclin oscillation is not essential for control of phase-specific transcription. Using a mathematical model, we demonstrate that the function of network TFs can be retained in the face of significant reductions in transcript levels. Finally, we show that cells arrested at mitotic exit with non-oscillating levels of B-cyclins continue to cycle transcriptionally. Taken together, these findings support a critical role of a TF network and a requirement for CDK activities that need not be periodic.
Collapse
Affiliation(s)
- Chun-Yi Cho
- Department of Biology, Duke University, Durham, NC, USA
| | | | | | | | | |
Collapse
|
40
|
Heffner KM, Hizal DB, Yerganian GS, Kumar A, Can Ö, O’Meally R, Cole R, Chaerkady R, Wu H, Bowen MA, Betenbaugh MJ. Lessons from the Hamster: Cricetulus griseus Tissue and CHO Cell Line Proteome Comparison. J Proteome Res 2017; 16:3672-3687. [DOI: 10.1021/acs.jproteome.7b00382] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
| | | | | | - Amit Kumar
- Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Özge Can
- Acibadem University, Medical Biochemistry, Istanbul, Maltepe, Turkey
| | - Robert O’Meally
- Johns Hopkins Medical Institute, Baltimore, Maryland 21205, United States
| | - Robert Cole
- Johns Hopkins Medical Institute, Baltimore, Maryland 21205, United States
| | | | - Herren Wu
- MedImmune, Gaithersburg, Maryland 20878, United States
| | | | | |
Collapse
|
41
|
Chen T, van Steensel B. Comprehensive analysis of nucleocytoplasmic dynamics of mRNA in Drosophila cells. PLoS Genet 2017; 13:e1006929. [PMID: 28771467 PMCID: PMC5557608 DOI: 10.1371/journal.pgen.1006929] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2016] [Revised: 08/15/2017] [Accepted: 07/17/2017] [Indexed: 01/14/2023] Open
Abstract
Eukaryotic mRNAs undergo a cycle of transcription, nuclear export, and degradation. A major challenge is to obtain a global, quantitative view of these processes. Here we measured the genome-wide nucleocytoplasmic dynamics of mRNA in Drosophila cells by metabolic labeling in combination with cellular fractionation. By mathematical modeling of these data we determined rates of transcription, export and cytoplasmic decay for 5420 genes. We characterized these kinetic rates and investigated links with mRNA features, RNA-binding proteins (RBPs) and chromatin states. We found prominent correlations between mRNA decay rate and transcript size, while nuclear export rates are linked to the size of the 3'UTR. Transcription, export and decay rates are each associated with distinct spectra of RBPs. Specific classes of genes, such as those encoding cytoplasmic ribosomal proteins, exhibit characteristic combinations of rate constants, suggesting modular control. Binding of splicing factors is associated with faster rates of export, and our data suggest coordinated regulation of nuclear export of specific functional classes of genes. Finally, correlations between rate constants suggest global coordination between the three processes. Our approach provides insights into the genome-wide nucleocytoplasmic kinetics of mRNA and should be generally applicable to other cell systems. All mRNAs start from production in the nucleus, undergo exportation through nuclear pores and finally are degraded in the cytoplasm. A comprehensive characterization of the kinetic rates of all mRNAs is an important prerequisite for a global understanding of the regulation of the transcriptome and the cell. By conducting a time-series experiment and building a mathematical model, we trace the dynamics of mRNAs from the nucleus to the cytoplasm and determine the rates at each kinetic step at transcriptome-wide level. This information allows us to associate mRNA kinetic rates with a wealth of biological features and made some intriguing discoveries. We show mRNA decay is positively linked to transcript length while mRNA export is negatively linked to the length of the 3' UTR. We show binding of splicing factors is associated with faster rates of mRNA export. We provide evidence for global coordination between nuclear export an decay of mRNA. We show genes sharing specific functions tend to have similar nucleoplasmic kinetics, in which ribosomal proteins possessing special kinetic features exclusively stand out. Altogether, our integrated approach to quantitatively determine the rates of kinetic steps on a gene-by-gene basis provides a blueprint to obtain the global understanding of RNA regulation.
Collapse
Affiliation(s)
- Tao Chen
- Division of Gene Regulation, Netherlands Cancer Institute, Amsterdam, The Netherlands
| | - Bas van Steensel
- Division of Gene Regulation, Netherlands Cancer Institute, Amsterdam, The Netherlands
- * E-mail:
| |
Collapse
|
42
|
Gibbs DL, Shmulevich I. Solving the influence maximization problem reveals regulatory organization of the yeast cell cycle. PLoS Comput Biol 2017; 13:e1005591. [PMID: 28628618 PMCID: PMC5495484 DOI: 10.1371/journal.pcbi.1005591] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2016] [Revised: 07/03/2017] [Accepted: 05/24/2017] [Indexed: 02/07/2023] Open
Abstract
The Influence Maximization Problem (IMP) aims to discover the set of nodes with the greatest influence on network dynamics. The problem has previously been applied in epidemiology and social network analysis. Here, we demonstrate the application to cell cycle regulatory network analysis for Saccharomyces cerevisiae. Fundamentally, gene regulation is linked to the flow of information. Therefore, our implementation of the IMP was framed as an information theoretic problem using network diffusion. Utilizing more than 26,000 regulatory edges from YeastMine, gene expression dynamics were encoded as edge weights using time lagged transfer entropy, a method for quantifying information transfer between variables. By picking a set of source nodes, a diffusion process covers a portion of the network. The size of the network cover relates to the influence of the source nodes. The set of nodes that maximizes influence is the solution to the IMP. By solving the IMP over different numbers of source nodes, an influence ranking on genes was produced. The influence ranking was compared to other metrics of network centrality. Although the top genes from each centrality ranking contained well-known cell cycle regulators, there was little agreement and no clear winner. However, it was found that influential genes tend to directly regulate or sit upstream of genes ranked by other centrality measures. The influential nodes act as critical sources of information flow, potentially having a large impact on the state of the network. Biological events that affect influential nodes and thereby affect information flow could have a strong effect on network dynamics, potentially leading to disease. Code and data can be found at: https://github.com/gibbsdavidl/miergolf.
Collapse
Affiliation(s)
- David L. Gibbs
- Institute for Systems Biology, Seattle, Washington, United States of America
| | - Ilya Shmulevich
- Institute for Systems Biology, Seattle, Washington, United States of America
| |
Collapse
|
43
|
Inoue M, Horimoto K. Relationship between regulatory pattern of gene expression level and gene function. PLoS One 2017; 12:e0177430. [PMID: 28494005 PMCID: PMC5426767 DOI: 10.1371/journal.pone.0177430] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2016] [Accepted: 04/27/2017] [Indexed: 01/03/2023] Open
Abstract
Regulation of gene expression levels is essential for all living systems and transcription factors (TFs) are the main regulators of gene expression through their ability to repress or induce transcription. A balance between synthesis and degradation rates controls gene expression levels. To determine which rate is dominant, we analyzed the correlation between expression levels of a TF and its regulated gene based on a mathematical model. We selected about 280,000 expression patterns of 355 TFs and 647 regulated genes using DNA microarray data from the Gene Expression Omnibus (GEO) data repository. Based on our model, correlation between the expressions of TF–regulated gene pairs corresponds to tuning of the synthesis rate, whereas no correlation indicates excessive synthesis and requires tuning of the degradation rate. The gene expression relationships between TF–regulated gene pairs were classified into four types that correspond to different gene regulatory mechanisms. It was surprising that fewer than 20% of these genes were governed by the familiar regulatory mechanism, i.e., through the synthesis rate. Moreover, we performed pathway analysis and found that each classification type corresponded to distinct gene functions: cellular regulation pathways were dominant in the type with synthesis rate regulation and terms associated with diseases such as cancer, Parkinson’s disease, and Alzheimer’s disease were dominant in the type with degradation rate regulation. Interestingly, these diseases are caused by the accumulation of proteins. These results indicated that gene expression is regulated structurally, not arbitrarily, according to the gene function. This funding is indicative of a systematic control of transcription processes at the whole-cell level.
Collapse
Affiliation(s)
- Masayo Inoue
- Molecular Profiling Research Center for Drug Discovery, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-ku, Tokyo 135-0064, Japan
- * E-mail:
| | - Katsuhisa Horimoto
- Molecular Profiling Research Center for Drug Discovery, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-ku, Tokyo 135-0064, Japan
| |
Collapse
|
44
|
Delineating functional principles of the bow tie structure of a kinase-phosphatase network in the budding yeast. BMC SYSTEMS BIOLOGY 2017; 11:38. [PMID: 28298210 PMCID: PMC5353956 DOI: 10.1186/s12918-017-0418-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/23/2016] [Accepted: 03/08/2017] [Indexed: 11/10/2022]
Abstract
Background Kinases and phosphatases (KP) form complex self-regulating networks essential for cellular signal processing. In spite of having a wealth of data about interactions among KPs and their substrates, we have very limited models of the structures of the directed networks they form and consequently our ability to formulate hypotheses about how their structure determines the flow of information in these networks is restricted. Results We assembled and studied the largest bona fide kinase-phosphatase network (KP-Net) known to date for the yeast Saccharomyces cerevisiae. Application of the vertex sort (VS) algorithm on the KP-Net allowed us to elucidate its hierarchical structure in which nodes are sorted into top, core and bottom layers, forming a bow tie structure with a strongly connected core layer. Surprisingly, phosphatases tend to sort into the top layer, implying they are less regulated by phosphorylation than kinases. Superposition of the widest range of KP biological properties over the KP-Net hierarchy shows that core layer KPs: (i), receive the largest number of inputs; (ii), form bottlenecks implicated in multiple pathways and in decision-making; (iii), and are among the most regulated KPs both temporally and spatially. Moreover, top layer KPs are more abundant and less noisy than those in the bottom layer. Finally, we showed that the VS algorithm depends on node degrees without biasing the biological results of the sorted network. The VS algorithm is available as an R package (https://cran.r-project.org/web/packages/VertexSort/index.html). Conclusions The KP-Net model we propose possesses a bow tie hierarchical structure in which the top layer appears to ensure highest fidelity and the core layer appears to mediate signal integration and cell state-dependent signal interpretation. Our model of the yeast KP-Net provides both functional insight into its organization as we understand today and a framework for future investigation of information processing in yeast and eukaryotes in general. Electronic supplementary material The online version of this article (doi:10.1186/s12918-017-0418-0) contains supplementary material, which is available to authorized users.
Collapse
|
45
|
Kelliher CM, Haase SB. Connecting virulence pathways to cell-cycle progression in the fungal pathogen Cryptococcus neoformans. Curr Genet 2017; 63:803-811. [PMID: 28265742 PMCID: PMC5605583 DOI: 10.1007/s00294-017-0688-5] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2017] [Revised: 02/22/2017] [Accepted: 02/22/2017] [Indexed: 11/01/2022]
Abstract
Proliferation and host evasion are critical processes to understand at a basic biological level for improving infectious disease treatment options. The human fungal pathogen Cryptococcus neoformans causes fungal meningitis in immunocompromised individuals by proliferating in cerebrospinal fluid. Current antifungal drugs target "virulence factors" for disease, such as components of the cell wall and polysaccharide capsule in C. neoformans. However, mechanistic links between virulence pathways and the cell cycle are not as well studied. Recently, cell-cycle synchronized C. neoformans cells were profiled over time to identify gene expression dynamics (Kelliher et al., PLoS Genet 12(12):e1006453, 2016). Almost 20% of all genes in the C. neoformans genome were periodically activated during the cell cycle in rich media, including 40 genes that have previously been implicated in virulence pathways. Here, we review important findings about cell-cycle-regulated genes in C. neoformans and provide two examples of virulence pathways-chitin synthesis and G-protein coupled receptor signaling-with their putative connections to cell division. We propose that a "comparative functional genomics" approach, leveraging gene expression timing during the cell cycle, orthology to genes in other fungal species, and previous experimental findings, can lead to mechanistic hypotheses connecting the cell cycle to fungal virulence.
Collapse
Affiliation(s)
- Christina M Kelliher
- Department of Biology, Duke University, Box 90338, 130 Science Drive, Durham, NC, 27708-0338, USA
| | - Steven B Haase
- Department of Biology, Duke University, Box 90338, 130 Science Drive, Durham, NC, 27708-0338, USA.
| |
Collapse
|
46
|
Histone H3K4 and H3K36 Methylation Independently Recruit the NuA3 Histone Acetyltransferase in Saccharomyces cerevisiae. Genetics 2017; 205:1113-1123. [PMID: 28108585 DOI: 10.1534/genetics.116.199422] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2016] [Accepted: 12/23/2016] [Indexed: 11/18/2022] Open
Abstract
Histone post-translational modifications (PTMs) alter chromatin structure by promoting the interaction of chromatin-modifying complexes with nucleosomes. The majority of chromatin-modifying complexes contain multiple domains that preferentially interact with modified histones, leading to speculation that these domains function in concert to target nucleosomes with distinct combinations of histone PTMs. In Saccharomyces cerevisiae, the NuA3 histone acetyltransferase complex contains three domains, the PHD finger in Yng1, the PWWP domain in Pdp3, and the YEATS domain in Taf14; which in vitro bind to H3K4 methylation, H3K36 methylation, and acetylated and crotonylated H3K9, respectively. While the in vitro binding has been well characterized, the relative in vivo contributions of these histone PTMs in targeting NuA3 is unknown. Here, through genome-wide colocalization and by mutational interrogation, we demonstrate that the PHD finger of Yng1, and the PWWP domain of Pdp3 independently target NuA3 to H3K4 and H3K36 methylated chromatin, respectively. In contrast, we find no evidence to support the YEATS domain of Taf14 functioning in NuA3 recruitment. Collectively our results suggest that the presence of multiple histone PTM binding domains within NuA3, rather than restricting it to nucleosomes containing distinct combinations of histone PTMs, can serve to increase the range of nucleosomes bound by the complex. Interestingly, however, the simple presence of NuA3 is insufficient to ensure acetylation of the associated nucleosomes, suggesting a secondary level of acetylation regulation that does not involve control of HAT-nucleosome interactions.
Collapse
|
47
|
Investigating Conservation of the Cell-Cycle-Regulated Transcriptional Program in the Fungal Pathogen, Cryptococcus neoformans. PLoS Genet 2016; 12:e1006453. [PMID: 27918582 PMCID: PMC5137879 DOI: 10.1371/journal.pgen.1006453] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2016] [Accepted: 11/01/2016] [Indexed: 12/24/2022] Open
Abstract
The pathogenic yeast Cryptococcus neoformans causes fungal meningitis in immune-compromised patients. Cell proliferation in the budding yeast form is required for C. neoformans to infect human hosts, and virulence factors such as capsule formation and melanin production are affected by cell-cycle perturbation. Thus, understanding cell-cycle regulation is critical for a full understanding of virulence factors for disease. Our group and others have demonstrated that a large fraction of genes in Saccharomyces cerevisiae is expressed periodically during the cell cycle, and that proper regulation of this transcriptional program is important for proper cell division. Despite the evolutionary divergence of the two budding yeasts, we found that a similar percentage of all genes (~20%) is periodically expressed during the cell cycle in both yeasts. However, the temporal ordering of periodic expression has diverged for some orthologous cell-cycle genes, especially those related to bud emergence and bud growth. Genes regulating DNA replication and mitosis exhibited a conserved ordering in both yeasts, suggesting that essential cell-cycle processes are conserved in periodicity and in timing of expression (i.e. duplication before division). In S. cerevisiae cells, we have proposed that an interconnected network of periodic transcription factors (TFs) controls the bulk of the cell-cycle transcriptional program. We found that temporal ordering of orthologous network TFs was not always maintained; however, the TF network topology at cell-cycle commitment appears to be conserved in C. neoformans. During the C. neoformans cell cycle, DNA replication genes, mitosis genes, and 40 genes involved in virulence are periodically expressed. Future work toward understanding the gene regulatory network that controls cell-cycle genes is critical for developing novel antifungals to inhibit pathogen proliferation.
Collapse
|
48
|
Abstract
We give an overview of experimental and computational methods to estimate RNA metabolism rates genome-wide. We then advocate a local definition of RNA metabolism rate at the level of individual phosphodiester bonds. Rates of formation and disappearance of individual bonds are unambiguously defined, in contrast to rates of complete transcripts. We show that over previous approaches, the recently developed transient transcriptome sequencing (TT-seq) protocol allows for estimation of metabolism rates of individual bonds with least positional bias.
Collapse
Affiliation(s)
- Leonhard Wachutka
- a Department of Informatics , Technical University of Munich, Garching bei München , Germany
| | - Julien Gagneur
- a Department of Informatics , Technical University of Munich, Garching bei München , Germany
| |
Collapse
|
49
|
Hillenbrand P, Maier KC, Cramer P, Gerland U. Inference of gene regulation functions from dynamic transcriptome data. eLife 2016; 5. [PMID: 27652904 PMCID: PMC5072840 DOI: 10.7554/elife.12188] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2015] [Accepted: 09/20/2016] [Indexed: 11/17/2022] Open
Abstract
To quantify gene regulation, a function is required that relates transcription factor binding to DNA (input) to the rate of mRNA synthesis from a target gene (output). Such a ‘gene regulation function’ (GRF) generally cannot be measured because the experimental titration of inputs and simultaneous readout of outputs is difficult. Here we show that GRFs may instead be inferred from natural changes in cellular gene expression, as exemplified for the cell cycle in the yeast S. cerevisiae. We develop this inference approach based on a time series of mRNA synthesis rates from a synchronized population of cells observed over three cell cycles. We first estimate the functional form of how input transcription factors determine mRNA output and then derive GRFs for target genes in the CLB2 gene cluster that are expressed during G2/M phase. Systematic analysis of additional GRFs suggests a network architecture that rationalizes transcriptional cell cycle oscillations. We find that a transcription factor network alone can produce oscillations in mRNA expression, but that additional input from cyclin oscillations is required to arrive at the native behaviour of the cell cycle oscillator. DOI:http://dx.doi.org/10.7554/eLife.12188.001 Living cells rely on networks of genes to control their behavior, including how they grow, develop and respond to stress. Genes encode instructions needed to make proteins and other molecules, and much of the control is exerted at the first stage of protein production, known as transcription. During this process, a gene is copied to make molecules known as transcripts that may later be used as templates to make proteins. Many genes encode proteins that act to regulate transcription. Therefore, an individual gene may receive inputs from other genes, and these inputs affect how much transcript the gene produces, which can be considered as the gene’s output. While these inputs and outputs can often be wired together to form a network, it is less clear exactly how all the different inputs at a gene interact to determine its output. These interactions are known as “gene regulation functions”, and knowing them would be an important step towards understanding gene networks, which would help us to predict how cells will behave in different situations. Gene regulation functions are difficult to measure directly, so researchers would like to find other ways to assess them indirectly. A recently developed experimental technique called “dynamic transcriptome analysis” seemed promising as it measures both the inputs and outputs of all genes in a cell over time. Hillenbrand et al. used this technique to infer gene regulation functions with one or two inputs in yeast cells. Comparing these estimates with experimental data from previous studies showed that these inferred gene regulation functions could successfully predict the output of a gene based on its inputs. Hillenbrand et al. then used these estimates to search and model a well-known genetic network that is thought to be part of the molecular clockwork that controls the timing of events that cause a cell to divide. Currently, the approach used by Hillenbrand et al. treats gene regulation functions like “black boxes”. This means that, while an output can be predicted if the inputs are known, it cannot reveal all of the detailed mechanisms behind it. Gaining insights into the inner workings of these black boxes will require taking more data into account, such as how abundant the proteins that regulate transcription are, where they are located within cells or whether they are active or not. Therefore, the next challenge is to incorporate these kinds of data to gain a fuller picture of how gene networks operate within cells. DOI:http://dx.doi.org/10.7554/eLife.12188.002
Collapse
Affiliation(s)
- Patrick Hillenbrand
- Lehrstuhl für Theorie komplexer Biosysteme, Physik-Department, Technische Universität München, Garching, Germany
| | - Kerstin C Maier
- Max-Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Patrick Cramer
- Max-Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Ulrich Gerland
- Lehrstuhl für Theorie komplexer Biosysteme, Physik-Department, Technische Universität München, Garching, Germany
| |
Collapse
|
50
|
Tyanova S, Temu T, Sinitcyn P, Carlson A, Hein MY, Geiger T, Mann M, Cox J. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat Methods 2016; 13:731-40. [DOI: 10.1038/nmeth.3901] [Citation(s) in RCA: 4028] [Impact Index Per Article: 503.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2016] [Accepted: 05/10/2016] [Indexed: 02/06/2023]
|