1
|
Coscujuela Tarrero L, Famà V, D'Andrea G, Maestri S, de Polo A, Biffo S, Furlan M, Pelizzola M. Nanodynamo quantifies subcellular RNA dynamics revealing extensive coupling between steps of the RNA life cycle. Nat Commun 2024; 15:7725. [PMID: 39231948 PMCID: PMC11375098 DOI: 10.1038/s41467-024-51917-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2023] [Accepted: 08/20/2024] [Indexed: 09/06/2024] Open
Abstract
The coordinated action of transcriptional and post-transcriptional machineries shapes gene expression programs at steady state and determines their concerted response to perturbations. We have developed Nanodynamo, an experimental and computational workflow for quantifying the kinetic rates of nuclear and cytoplasmic steps of the RNA life cycle. Nanodynamo is based on mathematical modelling following sequencing of native RNA from cellular fractions and polysomes. We have applied this workflow to triple-negative breast cancer cells, revealing widespread post-transcriptional RNA processing that is mutually exclusive with its co-transcriptional counterpart. We used Nanodynamo to unravel the coupling between transcription, processing, export, decay and translation machineries. We have identified a number of coupling interactions within and between the nucleus and cytoplasm that largely contribute to coordinating how cells respond to perturbations that affect gene expression programs. Nanodynamo will be instrumental in unravelling the determinants and regulatory processes involved in the coordination of gene expression responses.
Collapse
Affiliation(s)
| | - Valeria Famà
- Center for Genomic Science of IIT@SEMM, Istituto Italiano di Tecnologia (IIT), Milan, Italy
- Department of Oncology and Emato-Oncology, University of Milan, Milan, Italy
| | - Giacomo D'Andrea
- Istituto Nazionale Genetica Molecolare "Romeo ed Enrica Invernizzi", Milano, Italy
- Department of Biosciences, University of Milan, Milan, Italy
| | - Simone Maestri
- Center for Genomic Science of IIT@SEMM, Istituto Italiano di Tecnologia (IIT), Milan, Italy
| | - Anna de Polo
- Center for Genomic Science of IIT@SEMM, Istituto Italiano di Tecnologia (IIT), Milan, Italy
| | - Stefano Biffo
- Istituto Nazionale Genetica Molecolare "Romeo ed Enrica Invernizzi", Milano, Italy
- Department of Biosciences, University of Milan, Milan, Italy
| | - Mattia Furlan
- Center for Genomic Science of IIT@SEMM, Istituto Italiano di Tecnologia (IIT), Milan, Italy.
| | - Mattia Pelizzola
- Center for Genomic Science of IIT@SEMM, Istituto Italiano di Tecnologia (IIT), Milan, Italy.
- Department of Biotechnology and Biosciences, University of Milano-Bicocca, Milan, Italy.
| |
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
|
Rummel T, Sakellaridi L, Erhard F. grandR: a comprehensive package for nucleotide conversion RNA-seq data analysis. Nat Commun 2023; 14:3559. [PMID: 37321987 PMCID: PMC10272207 DOI: 10.1038/s41467-023-39163-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Accepted: 06/01/2023] [Indexed: 06/17/2023] Open
Abstract
Metabolic labeling of RNA is a powerful technique for studying the temporal dynamics of gene expression. Nucleotide conversion approaches greatly facilitate the generation of data but introduce challenges for their analysis. Here we present grandR, a comprehensive package for quality control, differential gene expression analysis, kinetic modeling, and visualization of such data. We compare several existing methods for inference of RNA synthesis rates and half-lives using progressive labeling time courses. We demonstrate the need for recalibration of effective labeling times and introduce a Bayesian approach to study the temporal dynamics of RNA using snapshot experiments.
Collapse
Affiliation(s)
- Teresa Rummel
- Institute for Virology and Immunobiology, University of Würzburg, Versbacher Str. 7, 97078, Würzburg, Germany
| | - Lygeri Sakellaridi
- Institute for Virology and Immunobiology, University of Würzburg, Versbacher Str. 7, 97078, Würzburg, Germany
| | - Florian Erhard
- Institute for Virology and Immunobiology, University of Würzburg, Versbacher Str. 7, 97078, Würzburg, Germany.
- Faculty for Informatics and Data Science, University of Regensburg, Bajuwarenstr. 4, 93053, Regensburg, Germany.
| |
Collapse
|
4
|
Forbes Beadle L, Love JC, Shapovalova Y, Artemev A, Rattray M, Ashe HL. Combined modelling of mRNA decay dynamics and single-molecule imaging in the Drosophila embryo uncovers a role for P-bodies in 5' to 3' degradation. PLoS Biol 2023; 21:e3001956. [PMID: 36649329 PMCID: PMC9882958 DOI: 10.1371/journal.pbio.3001956] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2022] [Revised: 01/27/2023] [Accepted: 12/13/2022] [Indexed: 01/18/2023] Open
Abstract
Regulation of mRNA degradation is critical for a diverse array of cellular processes and developmental cell fate decisions. Many methods for determining mRNA half-lives rely on transcriptional inhibition or metabolic labelling. Here, we use a non-invasive method for estimating half-lives for hundreds of mRNAs in the early Drosophila embryo. This approach uses the intronic and exonic reads from a total RNA-seq time series and Gaussian process regression to model the dynamics of premature and mature mRNAs. We show how regulation of mRNA stability is used to establish a range of mature mRNA dynamics during embryogenesis, despite shared transcription profiles. Using single-molecule imaging, we provide evidence that, for the mRNAs tested, there is a correlation between short half-life and mRNA association with P-bodies. Moreover, we detect an enrichment of mRNA 3' ends in P-bodies in the early embryo, consistent with 5' to 3' degradation occurring in P-bodies for at least a subset of mRNAs. We discuss our findings in relation to recently published data suggesting that the primary function of P-bodies in other biological contexts is mRNA storage.
Collapse
Affiliation(s)
- Lauren Forbes Beadle
- Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
| | - Jennifer C. Love
- Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
| | - Yuliya Shapovalova
- Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
| | - Artem Artemev
- Department of Computing, Imperial College London, London, United Kingdom
| | - Magnus Rattray
- Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
- * E-mail: (MR); (HLA)
| | - Hilary L. Ashe
- Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
- * E-mail: (MR); (HLA)
| |
Collapse
|
5
|
Hersch M, Biasini A, Marques AC, Bergmann S. Estimating RNA dynamics using one time point for one sample in a single-pulse metabolic labeling experiment. BMC Bioinformatics 2022; 23:147. [PMID: 35459101 PMCID: PMC9034570 DOI: 10.1186/s12859-022-04672-4] [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: 07/16/2021] [Accepted: 04/04/2022] [Indexed: 11/05/2022] Open
Abstract
Background Over the past decade, experimental procedures such as metabolic labeling for determining RNA turnover rates at the transcriptome-wide scale have been widely adopted and are now turning to single cell measurements. Several computational methods to estimate RNA synthesis, processing and degradation rates from such experiments have been suggested, but they all require several RNA sequencing samples. Here we present a method that can estimate those three rates from a single sample. Methods Our method relies on the analytical solution to the Zeisel model of RNA dynamics. It was validated on metabolic labeling experiments performed on mouse embryonic stem cells. Resulting degradation rates were compared both to previously published rates on the same system and to a state-of-the-art method applied to the same data. Results Our method is computationally efficient and outputs rates that correlate well with previously published data sets. Using it on a single sample, we were able to reproduce the observation that dynamic biological processes tend to involve genes with higher metabolic rates, while stable processes involve genes with lower rates. This supports the hypothesis that cells control not only the mRNA steady-state abundance, but also its responsiveness, i.e., how fast steady state is reached. Moreover, degradation rates obtained with our method compare favourably with the other tested method. Conclusions In addition to saving experimental work and computational time, estimating rates for a single sample has several advantages. It does not require an error-prone normalization across samples and enables the use of replicates to estimate uncertainty and assess sample quality. Finally the method and theoretical results described here are general enough to be useful in other contexts such as nucleotide conversion methods and single cell metabolic labeling experiments. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-022-04672-4.
Collapse
Affiliation(s)
- Micha Hersch
- Department of Computational Biology, University of Lausanne, Lausanne, Switzerland. .,Swiss Institute of Bioinformatics, 1015, Lausanne, CH, Switzerland.
| | - Adriano Biasini
- Department of Computational Biology, University of Lausanne, Lausanne, Switzerland.,RNA Therapeutics Institute, University of Massachusetts Medical School, Worcester, MA, USA
| | - Ana C Marques
- Department of Computational Biology, University of Lausanne, Lausanne, Switzerland
| | - Sven Bergmann
- Department of Computational Biology, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, 1015, Lausanne, CH, Switzerland
| |
Collapse
|
6
|
Furlan M, de Pretis S, Pelizzola M. Dynamics of transcriptional and post-transcriptional regulation. Brief Bioinform 2021; 22:bbaa389. [PMID: 33348360 PMCID: PMC8294512 DOI: 10.1093/bib/bbaa389] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Revised: 11/12/2020] [Accepted: 11/27/2020] [Indexed: 02/07/2023] Open
Abstract
Despite gene expression programs being notoriously complex, RNA abundance is usually assumed as a proxy for transcriptional activity. Recently developed approaches, able to disentangle transcriptional and post-transcriptional regulatory processes, have revealed a more complex scenario. It is now possible to work out how synthesis, processing and degradation kinetic rates collectively determine the abundance of each gene's RNA. It has become clear that the same transcriptional output can correspond to different combinations of the kinetic rates. This underscores the fact that markedly different modes of gene expression regulation exist, each with profound effects on a gene's ability to modulate its own expression. This review describes the development of the experimental and computational approaches, including RNA metabolic labeling and mathematical modeling, that have been disclosing the mechanisms underlying complex transcriptional programs. Current limitations and future perspectives in the field are also discussed.
Collapse
Affiliation(s)
- Mattia Furlan
- Center for Genomic Science of IIT@SEMM, Fondazione Istituto Italiano di Tecnologia (IIT), 20139 Milan, Italy
| | - Stefano de Pretis
- Center for Genomic Science of IIT@SEMM, Fondazione Istituto Italiano di Tecnologia (IIT), 20139 Milan, Italy
| | - Mattia Pelizzola
- Center for Genomic Science of IIT@SEMM, Fondazione Istituto Italiano di Tecnologia (IIT), 20139 Milan, Italy
| |
Collapse
|
7
|
Furlan M, Galeota E, Gaudio ND, Dassi E, Caselle M, de Pretis S, Pelizzola M. Genome-wide dynamics of RNA synthesis, processing, and degradation without RNA metabolic labeling. Genome Res 2020; 30:1492-1507. [PMID: 32978246 PMCID: PMC7605262 DOI: 10.1101/gr.260984.120] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2020] [Accepted: 08/21/2020] [Indexed: 12/13/2022]
Abstract
The quantification of the kinetic rates of RNA synthesis, processing, and degradation are largely based on the integrative analysis of total and nascent transcription, the latter being quantified through RNA metabolic labeling. We developed INSPEcT−, a computational method based on the mathematical modeling of premature and mature RNA expression that is able to quantify kinetic rates from steady-state or time course total RNA-seq data without requiring any information on nascent transcripts. Our approach outperforms available solutions, closely recapitulates the kinetic rates obtained through RNA metabolic labeling, improves the ability to detect changes in transcript half-lives, reduces the cost and complexity of the experiments, and can be adopted to study experimental conditions in which nascent transcription cannot be readily profiled. Finally, we applied INSPEcT− to the characterization of post-transcriptional regulation landscapes in dozens of physiological and disease conditions. This approach was included in the INSPEcT Bioconductor package, which can now unveil RNA dynamics from steady-state or time course data, with or without the profiling of nascent RNA.
Collapse
Affiliation(s)
- Mattia Furlan
- Center for Genomic Science, Fondazione Istituto Italiano di Tecnologia, 20139 Milan, Italy.,Physics Department and INFN, University of Turin, 10125 Turin, Italy
| | - Eugenia Galeota
- Center for Genomic Science, Fondazione Istituto Italiano di Tecnologia, 20139 Milan, Italy
| | - Nunzio Del Gaudio
- Center for Genomic Science, Fondazione Istituto Italiano di Tecnologia, 20139 Milan, Italy
| | - Erik Dassi
- Centre for Integrative Biology, University of Trento, 38123 Trento, Italy
| | - Michele Caselle
- Physics Department and INFN, University of Turin, 10125 Turin, Italy
| | - Stefano de Pretis
- Center for Genomic Science, Fondazione Istituto Italiano di Tecnologia, 20139 Milan, Italy
| | - Mattia Pelizzola
- Center for Genomic Science, Fondazione Istituto Italiano di Tecnologia, 20139 Milan, Italy
| |
Collapse
|
8
|
Biasini A, Marques AC. A Protocol for Transcriptome-Wide Inference of RNA Metabolic Rates in Mouse Embryonic Stem Cells. Front Cell Dev Biol 2020; 8:97. [PMID: 32175319 PMCID: PMC7056730 DOI: 10.3389/fcell.2020.00097] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2019] [Accepted: 02/07/2020] [Indexed: 12/16/2022] Open
Abstract
The relative ease of mouse Embryonic Stem Cells (mESCs) culture and the potential of these cells to differentiate into any of the three primary germ layers: ectoderm, endoderm and mesoderm (pluripotency), makes them an ideal and frequently used ex vivo system to dissect how gene expression changes impact cell state and differentiation. These efforts are further supported by the large number of constitutive and inducible mESC mutants established with the aim of assessing the contributions of different pathways and genes to cell homeostasis and gene regulation. Gene product abundance is controlled by the modulation of the rates of RNA synthesis, processing, and degradation. The ability to determine the relative contribution of these different RNA metabolic rates to gene expression control using standard RNA-sequencing approaches, which only capture steady state abundance of transcripts, is limited. In contrast, metabolic labeling of RNA with 4-thiouridine (4sU) coupled with RNA-sequencing, allows simultaneous and reproducible inference of transcriptome wide synthesis, processing, and degradation rates. Here we describe, a detailed protocol for 4sU metabolic labeling in mESCs that requires short 4sU labeling times at low concentration and minimally impacts cellular homeostasis. This approach presents a versatile method for in-depth characterization of the gene regulatory strategies governing gene steady state abundance in mESC.
Collapse
Affiliation(s)
- Adriano Biasini
- Department of Computational Biology, University of Lausanne, Lausanne, Switzerland
| | - Ana Claudia Marques
- Department of Computational Biology, University of Lausanne, Lausanne, Switzerland
| |
Collapse
|
9
|
Desgranges E, Caldelari I, Marzi S, Lalaouna D. Navigation through the twists and turns of RNA sequencing technologies: Application to bacterial regulatory RNAs. BIOCHIMICA ET BIOPHYSICA ACTA-GENE REGULATORY MECHANISMS 2020; 1863:194506. [PMID: 32068131 DOI: 10.1016/j.bbagrm.2020.194506] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/14/2019] [Revised: 02/11/2020] [Accepted: 02/13/2020] [Indexed: 12/20/2022]
Abstract
Discovered in the 1980s, small regulatory RNAs (sRNAs) are now considered key actors in virtually all aspects of bacterial physiology and virulence. Together with transcriptional and translational regulatory proteins, they integrate and often are hubs of complex regulatory networks, responsible for bacterial response/adaptation to various perceived stimuli. The recent development of powerful RNA sequencing technologies has facilitated the identification and characterization of sRNAs (length, structure and expression conditions) and their RNA targets in several bacteria. Nevertheless, it could be very difficult for non-experts to understand the advantages and drawbacks related to each offered option and, consequently, to make an informed choice. Therefore, the main goal of this review is to provide a guide to navigate through the twists and turns of high-throughput RNA sequencing technologies, with a specific focus on those applied to the study of sRNAs. This article is part of a Special Issue entitled: RNA and gene control in bacteria edited by Dr. M. Guillier and F. Repoila.
Collapse
Affiliation(s)
- Emma Desgranges
- Université de Strasbourg, CNRS, ARN UPR 9002, F-67000 Strasbourg, France
| | - Isabelle Caldelari
- Université de Strasbourg, CNRS, ARN UPR 9002, F-67000 Strasbourg, France
| | - Stefano Marzi
- Université de Strasbourg, CNRS, ARN UPR 9002, F-67000 Strasbourg, France
| | - David Lalaouna
- Université de Strasbourg, CNRS, ARN UPR 9002, F-67000 Strasbourg, France.
| |
Collapse
|