1
|
Vock IW, Simon MD. bakR: uncovering differential RNA synthesis and degradation kinetics transcriptome-wide with Bayesian hierarchical modeling. RNA 2023; 29:958-976. [PMID: 37028916 DOI: 10.1261/rna.079451.122] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/16/2022] [Accepted: 03/14/2023] [Indexed: 06/18/2023]
Abstract
Differential expression analysis of RNA sequencing (RNA-seq) data can identify changes in cellular RNA levels, but provides limited information about the kinetic mechanisms underlying such changes. Nucleotide recoding RNA-seq methods (NR-seq; e.g., TimeLapse-seq, SLAM-seq, etc.) address this shortcoming and are widely used approaches to identify changes in RNA synthesis and degradation kinetics. While advanced statistical models implemented in user-friendly software (e.g., DESeq2) have ensured the statistical rigor of differential expression analyses, no such tools that facilitate differential kinetic analysis with NR-seq exist. Here, we report the development of Bayesian analysis of the kinetics of RNA (bakR; https:// github.com/simonlabcode/bakR), an R package to address this need. bakR relies on Bayesian hierarchical modeling of NR-seq data to increase statistical power by sharing information across transcripts. Analyses of simulated data confirmed that bakR implementations of the hierarchical model outperform attempts to analyze differential kinetics with existing models. bakR also uncovers biological signals in real NR-seq data sets and provides improved analyses of existing data sets. This work establishes bakR as an important tool for identifying differential RNA synthesis and degradation kinetics.
Collapse
Affiliation(s)
- Isaac W Vock
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut 06536, USA
- Institute of Biomolecular Design and Discovery, Yale University, West Haven, Connecticut 06477, USA
| | - Matthew D Simon
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut 06536, USA
- Institute of Biomolecular Design and Discovery, Yale University, West Haven, Connecticut 06477, USA
| |
Collapse
|
2
|
Petkovic S, Graff S, Feller N, Berghaus J, Ruppert VP, Dülfer J, Sczakiel G. Circular versus linear RNA topology: different modes of RNA-RNA interactions in vitro and in human cells. RNA Biol 2021; 18:674-683. [PMID: 34839802 PMCID: PMC8782184 DOI: 10.1080/15476286.2021.1978214] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
Circular RNA is progressively reported to occur in various species including mammals where it is thought to be involved in the post-transcriptional regulation of gene expression, partly via interactions with microRNA. Here, we asked whether the circular topology causes functional differences to linear forms when interacting with short RNA strands in vitro and in human cells. Kinetic studies with human bladder cancer-derived synthetic circular RNA versus linear transcripts, respectively, with short oligoribonucleotides showed similar association rates for both topologies. Conversely, a substantial topology-related difference was measured for the activation entropy and the activation enthalpy of RNA–RNA annealing. This finding strongly indicates a significant difference of the mechanism of RNA–RNA interactions. To investigate whether these characteristics of circular RNA are biologically meaningful we performed transient transfection experiments with a microRNA-regulated expression system for luciferase in bladder cancer-derived cells. We co-transfected linear or circular RNA containing one microRNA binding site for the target-suppressing microRNA mlet7a. Here, the circular isoform showed a strongly increased competition with microRNA function versus linear versions. In summary, this study suggests novel topology-related characteristics of RNA–RNA interactions involving circRNA in vitro and in living cells.
Collapse
Affiliation(s)
- Sonja Petkovic
- Institut für Molekulare Medizin, Universität zu Lübeck and UKSH, Campus Lübeck, Lübeck, Germany
| | - Sarah Graff
- Institut für Molekulare Medizin, Universität zu Lübeck and UKSH, Campus Lübeck, Lübeck, Germany
| | - Nina Feller
- Institut für Molekulare Medizin, Universität zu Lübeck and UKSH, Campus Lübeck, Lübeck, Germany
| | - Julia Berghaus
- Institut für Molekulare Medizin, Universität zu Lübeck and UKSH, Campus Lübeck, Lübeck, Germany
| | | | - Jasmin Dülfer
- Institut für Molekulare Medizin, Universität zu Lübeck and UKSH, Campus Lübeck, Lübeck, Germany
| | - Georg Sczakiel
- Institut für Molekulare Medizin, Universität zu Lübeck and UKSH, Campus Lübeck, Lübeck, Germany
| |
Collapse
|
3
|
Findeiß S, Hammer S, Wolfinger MT, Kühnl F, Flamm C, Hofacker IL. In silico design of ligand triggered RNA switches. Methods 2018; 143:90-101. [PMID: 29660485 DOI: 10.1016/j.ymeth.2018.04.003] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2017] [Revised: 03/06/2018] [Accepted: 04/06/2018] [Indexed: 02/06/2023] Open
Abstract
This contribution sketches a work flow to design an RNA switch that is able to adapt two structural conformations in a ligand-dependent way. A well characterized RNA aptamer, i.e., knowing its Kd and adaptive structural features, is an essential ingredient of the described design process. We exemplify the principles using the well-known theophylline aptamer throughout this work. The aptamer in its ligand-binding competent structure represents one structural conformation of the switch while an alternative fold that disrupts the binding-competent structure forms the other conformation. To keep it simple we do not incorporate any regulatory mechanism to control transcription or translation. We elucidate a commonly used design process by explicitly dissecting and explaining the necessary steps in detail. We developed a novel objective function which specifies the mechanistics of this simple, ligand-triggered riboswitch and describe an extensive in silico analysis pipeline to evaluate important kinetic properties of the designed sequences. This protocol and the developed software can be easily extended or adapted to fit novel design scenarios and thus can serve as a template for future needs.
Collapse
Affiliation(s)
- Sven Findeiß
- Bioinformatics, Institute of Computer Science, and Interdisciplinary Center for Bioinformatics, Leipzig University, Härtelstraße 16-18, 04107 Leipzig, Germany; University of Vienna, Faculty of Computer Science, Research Group Bioinformatics and Computational Biology, Währingerstraße 29, 1090 Vienna, Austria; University of Vienna, Faculty of Chemistry, Department of Theoretical Chemistry, Währingerstraße 17, 1090 Vienna, Austria.
| | - Stefan Hammer
- Bioinformatics, Institute of Computer Science, and Interdisciplinary Center for Bioinformatics, Leipzig University, Härtelstraße 16-18, 04107 Leipzig, Germany; University of Vienna, Faculty of Computer Science, Research Group Bioinformatics and Computational Biology, Währingerstraße 29, 1090 Vienna, Austria; University of Vienna, Faculty of Chemistry, Department of Theoretical Chemistry, Währingerstraße 17, 1090 Vienna, Austria
| | - Michael T Wolfinger
- University of Vienna, Faculty of Chemistry, Department of Theoretical Chemistry, Währingerstraße 17, 1090 Vienna, Austria; Medical University of Vienna, Center for Anatomy and Cell Biology, Währingerstraße 13, 1090 Vienna, Austria
| | - Felix Kühnl
- Bioinformatics, Institute of Computer Science, and Interdisciplinary Center for Bioinformatics, Leipzig University, Härtelstraße 16-18, 04107 Leipzig, Germany
| | - Christoph Flamm
- University of Vienna, Faculty of Chemistry, Department of Theoretical Chemistry, Währingerstraße 17, 1090 Vienna, Austria
| | - Ivo L Hofacker
- University of Vienna, Faculty of Computer Science, Research Group Bioinformatics and Computational Biology, Währingerstraße 29, 1090 Vienna, Austria; University of Vienna, Faculty of Chemistry, Department of Theoretical Chemistry, Währingerstraße 17, 1090 Vienna, Austria
| |
Collapse
|