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.

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.

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.