Tutorials › Statistical Methods for Single-Cell RNA Sequencing › Transcription, Splicing, and RNA Degradation

RNA Velocity · Tutorial 384 of 700

Transcription, Splicing, and RNA Degradation

Learn how transcription, splicing, and degradation move RNA through distinct processes, and how to interpret their rates without confusing a rate constant with an RNA count or flux.

Advanced 10 min read

What You'll Learn

  • Distinguish transcriptional production from splicing and RNA removal.
  • Interpret transcription rates, per-molecule rate constants, and molecular fluxes with appropriate units.
  • Calculate splicing flux from an unspliced pool and a first-order splicing rate.
  • Derive and use the half-life associated with a first-order degradation rate.
  • Recognize why observed spliced and unspliced counts do not directly identify kinetic rates.
  • Identify biological and measurement conditions that can make simple rate interpretations misleading.

Three Processes, Three Different Roles

Spliced and unspliced RNA counts describe RNA molecules at different processing stages, but a count is not itself a rate. To interpret RNA dynamics, distinguish the processes that produce, convert, and remove RNA. Transcription produces a precursor transcript; splicing removes intronic sequence and converts eligible unspliced RNA into spliced RNA; and degradation removes RNA molecules. These processes can affect different RNA pools and need not proceed at the same pace.

The introductory material on spliced and unspliced RNA counts establishes that sequencing-based measurements can be separated into these two categories. Here, the focus is the underlying biological interpretation of each process. For clarity, let \(U\) denote the number of unspliced molecules of a gene in a cell at a particular time, and let \(S\) denote its number of spliced molecules. These are molecular abundances in the biological description; observed sequencing counts are measurements related to those abundances, not guaranteed exact molecule counts.

Definition: A transcription rate describes the production of RNA molecules per unit time. A splicing rate constant describes the rate of conversion per unspliced molecule, and a degradation rate constant describes the rate of removal per molecule in a specified RNA pool. Multiplying a per-molecule rate constant by the number of molecules in that pool gives an instantaneous molecular flux, measured in molecules per unit time.

A useful distinction is therefore between a rate constant and a flux. A transcription rate such as 12 molecules per minute already describes a production flux. By contrast, a splicing rate constant such as \(0.08\) per minute must be combined with the amount of eligible unspliced RNA to obtain a splicing flux. If there are \(U\) eligible molecules and each has the same first-order splicing rate constant \(\beta\), the expected instantaneous splicing flux is \(\beta U\) molecules per unit time.

Transcription Produces Precursor RNA

Transcription is the process by which RNA is synthesized from a DNA template. In the simplified picture used to interpret spliced and unspliced RNA, the immediate product is a precursor that contains intronic sequence and can contribute to the unspliced pool. The effective transcription rate, often denoted \(\alpha\), is expressed in molecules per unit time. It is a production rate, not a per-molecule rate constant: it does not need to be multiplied by the existing RNA abundance to determine the production flux.

This effective rate should not automatically be equated with promoter initiation frequency. A transcript must be initiated and synthesized, and experimental observations may reflect delays or incomplete transcripts. In a simplified kinetic description, \(\alpha\) summarizes the rate at which completed precursor molecules enter the RNA pool being considered. The interpretation depends on what counts as entry into that pool and on whether transcription is approximately stable over the time interval of interest.

Transcription can vary across genes, cells, and time. A larger unspliced count does not, by itself, prove that transcription is faster: the count could also be larger because splicing or degradation is slower, or because the cell was sampled during a transient change. Counts are snapshots of abundance, while transcription is a process rate.

Worked Example: Interpreting a Transcription Rate

In a hypothetical cell, suppose the effective transcription rate for a gene is estimated as 18 completed precursor molecules per minute. Interpret this number and compare it with an unspliced abundance of 100 molecules.

The rate \(\alpha=18\) has units of molecules per minute. Under the stated interpretation, it means that the process produces 18 precursor molecules per minute on average over the interval represented by the estimate. It is already a production flux.

$$ \alpha = 18\ \text{molecules/minute} $$

The unspliced abundance \(U=100\) molecules has units of molecules, not molecules per minute. It is therefore not valid to compare 18 and 100 as if they were two rates, or to infer that transcription is \(18/100\) per minute. Dividing the production rate by the abundance would give \(0.18\) per minute, but that ratio is not the transcription rate constant defined here; it is only a rate-to-abundance ratio and requires a separate model and interpretation.

Splicing Converts RNA Between Pools

Splicing removes intronic regions from a precursor transcript. In the two-pool description, it is a conversion: an eligible molecule leaves the unspliced pool and contributes to the spliced pool. It does not create a new RNA molecule from nothing, nor does it itself mean that the RNA has been degraded. Real transcripts may undergo multiple intron-removal events, so a single effective splicing rate can summarize a more complicated sequence of processing steps.

Under a first-order approximation, each molecule in the eligible unspliced pool has the same instantaneous chance of being spliced per unit time. If the rate constant is \(\beta\), with units of inverse time, and the pool contains \(U\) eligible molecules, then the expected instantaneous conversion flux is

$$ \text{splicing flux} = \beta U. $$

The units confirm the interpretation: \((1/\text{time})\times(\text{molecules})=\text{molecules}/\text{time}\). This calculation depends on the first-order approximation and on \(U\) referring to molecules that are eligible for the modeled conversion. If processing is delayed, depends on transcript features, or occurs in multiple stages, one effective \(\beta\) may not capture all relevant biology.

Worked Example: Converting a Splicing Rate Constant to a Flux

Suppose a hypothetical cell has \(U=120\) eligible unspliced molecules of a gene. Assume a first-order splicing rate constant of \(\beta=0.08\) per minute. Calculate the expected instantaneous splicing flux.

The formula is the rate constant multiplied by the number of eligible molecules:

$$ \beta U = (0.08\ \text{minute}^{-1})(120\ \text{molecules}) = 9.6\ \text{molecules/minute}. $$

Thus, under the stated approximation, the instantaneous expected conversion flux is 9.6 molecules per minute. This is not a claim that exactly 9.6 molecules will be observed to splice during every one-minute interval; actual molecular events are discrete and variable. It is also not the splicing rate constant: that remains \(0.08\) per minute.

Degradation Removes RNA

Degradation is the removal of RNA molecules. It can act on unspliced RNA, spliced RNA, or both, and the rates need not be equal. For example, a precursor that is incorrectly processed may be removed, while a mature transcript may have a different lifetime determined by its sequence, cellular context, and interactions with RNA-binding molecules. It is therefore often more biologically meaningful to distinguish degradation rates for different pools than to speak of one universal RNA degradation rate.

Under a first-order degradation assumption, each molecule in a specified pool has a constant instantaneous removal hazard. If \(\gamma\) is the degradation rate constant, its units are inverse time, and the expected instantaneous degradation flux from a pool of \(N\) molecules is \(\gamma N\) molecules per unit time. This assumption means the per-molecule removal rate is constant over the interval being modeled; it does not mean the pool size itself is constant.

A consequence of a constant first-order hazard is an exponential survival curve. If \(P(T>t)\) is the probability that an individual RNA molecule remains undegraded beyond time \(t\), then over a short interval \(\Delta t\), the probability of surviving is approximately multiplied by \(1-\gamma\Delta t\). Taking the continuous-time limit gives the survival function below. The half-life is the time at which this survival probability equals one half.

$$ P(T>t)=e^{-\gamma t}, \qquad \frac{1}{2}=e^{-\gamma t_{1/2}}, \qquad t_{1/2}=\frac{\ln 2}{\gamma}. $$

This derivation relies on a constant per-molecule hazard and a homogeneous interpretation of the pool. If molecules have different hazards, or the cellular environment changes during the measurement, a single exponential survival curve and a single half-life may be inadequate summaries.

Worked Example: Comparing Degradation Rates and Half-Lives

For a hypothetical gene, suppose the first-order degradation rate constants are \(\gamma_u=0.20\) per hour for unspliced RNA and \(\gamma_s=0.035\) per hour for spliced RNA. Calculate both half-lives and interpret their difference.

For unspliced RNA, the first-order half-life formula gives

$$ t_{1/2,u} = \frac{\ln 2}{\gamma_u} = \frac{0.6931}{0.20\ \text{hour}^{-1}} = 3.466\ \text{hours}. $$

For spliced RNA, the calculation is

$$ t_{1/2,s} = \frac{\ln 2}{\gamma_s} = \frac{0.6931}{0.035\ \text{hour}^{-1}} = 19.80\ \text{hours}. $$

Under the assumed constant-hazard model, half of an initially considered group of unspliced molecules would remain undegraded after about 3.47 hours, while half of the spliced molecules would remain after about 19.8 hours. These are model-based lifetime interpretations, not claims that every molecule has exactly that lifetime. The distinct rate constants imply different removal dynamics for the two pools; they do not, by themselves, establish why the rates differ biologically.

Rate Assumptions, Measurements, and Failure Modes

These rate interpretations rely on specific assumptions. A first-order splicing or degradation constant assumes a proportional relationship between pool size and instantaneous flux, with the same effective per-molecule rate for molecules in the pool. A transcription rate expressed in molecules per time assumes a defined interval and a consistent definition of the precursor product. Each rate is local to its biological context: it may change with cell state, gene, time, or experimental condition.

Sequencing counts do not directly reveal these rate constants. Counts are affected by capture efficiency, sequencing depth, transcript features, and the probability that a molecule is assigned to the unspliced or spliced category. A count of zero can reflect low abundance, incomplete capture, or sampling variability; it is not proof that the corresponding biological process has stopped. Likewise, a large unspliced count is compatible with high transcription, slow splicing, slow degradation, or combinations of these mechanisms.

A common mistake is to interpret a large spliced-to-unspliced count ratio as direct evidence that splicing is fast. The ratio compares two abundances, not a conversion rate. The pools may have different degradation rates, and the cell may not be at a stable state. Rate inference requires additional assumptions and information; it cannot be obtained from a single abundance ratio alone.

A second diagnostic is dimensional consistency. A number reported in molecules per minute is a flux or production rate; a number reported in inverse minutes is a per-molecule rate constant; and a molecule count is an abundance. Multiplying an inverse-time constant by a pool abundance gives a flux. Confusing these units can lead to incorrect biological conclusions even when the arithmetic is correct.

Experimental strategies can provide information beyond a static snapshot. For example, pulse-labeling or pulse-chase designs track newly synthesized RNA over time and can help distinguish production and removal processes. Such measurements still require assumptions about labeling, processing, and sample comparability. A simple first-order summary may be useful, but it can obscure multiple introns, processing delays, transcript-specific stability, or cell-to-cell heterogeneity.

Comparing the Quantities

QuantityBiological roleTypical unitsExample interpretation
Transcription rate \(\alpha\)Production of precursor RNAMolecules per time18 precursor molecules produced per minute
Splicing rate constant \(\beta\)Per-molecule conversion from unspliced RNAInverse timeEach eligible molecule has an effective splicing rate of 0.08 per minute
Splicing flux \(\beta U\)Conversion events per time in the current poolMolecules per time9.6 molecules per minute for \(U=120\) and \(\beta=0.08\) per minute
Degradation rate constant \(\gamma\)Per-molecule removal from a specified poolInverse timeA constant hazard of 0.035 per hour
Abundance \(U\) or \(S\)Number of molecules in a poolMolecules120 unspliced molecules at the observation time
Key takeaway: Transcription supplies precursor RNA, splicing converts eligible unspliced RNA into spliced RNA, and degradation removes molecules from specified pools. Keep production rates, per-molecule rate constants, fluxes, and measured abundances distinct; they answer different biological questions.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. A gene has 80 eligible unspliced molecules and a first-order splicing rate constant of 0.05 per minute. What is the expected instantaneous splicing flux?
  2. A spliced RNA has a first-order degradation rate constant of 0.10 per hour. What is its half-life, to three significant figures?
  3. Why is a transcription rate of 15 molecules per minute not interchangeable with a splicing rate constant of 0.15 per minute?
  4. A cell has more unspliced RNA than another cell. Does this establish that transcription is faster in the first cell? Give two alternative explanations.

Solutions. (1) The flux is \(\beta U=(0.05\ \text{minute}^{-1})(80\ \text{molecules})=4.0\) molecules per minute. (2) \(t_{1/2}=\ln(2)/\gamma=0.6931/(0.10\ \text{hour}^{-1})=6.93\) hours, under the constant-hazard assumption. (3) The transcription rate is a production flux with units of molecules per minute; the splicing constant is a per-molecule rate with units of inverse minutes and must be multiplied by an eligible pool size to obtain a flux. (4) No. Slower splicing or slower degradation of unspliced RNA could also increase its abundance; a transient change in transcription is another possibility.

Clinical Biostats