Sparse Counts and Measurement Noise

Why Sparse Counts Complicate Velocity

RNA velocity is inferred from the balance of unspliced and spliced RNA, so a few observed molecules can have a large effect on a cell's estimated direction of change. In a sparse gene-by-cell count matrix, zero may mean that no molecule was captured, not that the cell contained none. And if intronic reads are less likely to be observed than spliced reads, the two counts cannot be compared as though they had equal coverage.

The earlier material on the dynamical RNA velocity model defines the kinetic quantities used here. In brief, let $u$ and $s$ denote latent unspliced and spliced abundance for a gene in a cell. With transcription rate $\alpha$, splicing rate $\beta$, and degradation rate $\gamma$, the model is

\frac{du}{dt}=\alpha-\beta u,\qquad \frac{ds}{dt}=\beta u-\gamma s.

The spliced-RNA velocity is $v_s=\beta u-\gamma s$. This tutorial focuses on how noisy counts affect estimates of $u$, $s$, and $v_s$, not on re-estimating the kinetic model itself.

A Count Model for Unequal Detection

Definition: For one gene in one cell, let $X_U$ and $X_S$ be the observed unspliced and spliced counts, and let $u$ and $s$ be their latent abundances on a common scale. Write $a_U$ and $a_S$ for their respective detection exposures: the expected observed count per unit of latent abundance. Under a simple Poisson measurement model, $$ X_U\mid u\sim\operatorname{Poisson}(a_Uu),\qquad X_S\mid s\sim\operatorname{Poisson}(a_Ss). $$

The exposures may reflect sequencing depth, capture efficiency, and the fraction of reads available to detect each RNA class. They need not be equal: intronic reads may have lower effective coverage than spliced reads, giving $a_U<a_S$. The basic model assumes that the exposures are known or estimated well enough to use as offsets, counts are conditionally Poisson given latent abundance and exposure, and the two counts are conditionally independent. These are assumptions about measurement, not guarantees that biological variation is Poisson.

An exposure-adjusted estimate is $\hat u=X_U/a_U$ and $\hat s=X_S/a_S$. Under the model, each is unbiased conditional on its latent abundance and exposure. An exposure of 1 gives the familiar count-as-abundance estimate; a smaller exposure means that the same observed count corresponds to a larger estimated latent abundance. The velocity estimate is

\hat v_s=\beta\hat u-\gamma\hat s.

How Count Noise Enters the Velocity

For a Poisson random variable, the conditional variance equals its mean. Consequently, even if the latent abundance is fixed, repeated sequencing or capture would produce variable counts. The following derivation gives the uncertainty this induces in the adjusted velocity when the kinetic rates and exposures are treated as fixed.

\operatorname{Var}(\hat u\mid u) =\frac{\operatorname{Var}(X_U\mid u)}{a_U^2} =\frac{a_Uu}{a_U^2} =\frac{u}{a_U}.

Replacing the unknown $u$ by its count-based estimate gives the plug-in variance estimate $X_U/a_U^2$; likewise, the estimate for $\hat s$ is $X_S/a_S^2$. If the two counts are conditionally independent, the variance of the velocity estimate is the sum of the two scaled variances:

\operatorname{Var}(\hat v_s\mid u,s) =\beta^2\frac{u}{a_U}+\gamma^2\frac{s}{a_S}, \qquad \widehat{\operatorname{Var}}(\hat v_s) =\beta^2\frac{X_U}{a_U^2}+\gamma^2\frac{X_S}{a_S^2}.

If the count errors are correlated, a covariance term must be included: the variance becomes $\beta^2\operatorname{Var}(\hat u)+\gamma^2\operatorname{Var}(\hat s)-2\beta\gamma\operatorname{Cov}(\hat u,\hat s)$. These formulas describe count-sampling uncertainty only. Uncertainty in $\beta$ and $\gamma$, uncertainty in exposure estimates, and cell-to-cell biological variability add further uncertainty.

Low exposure increases variance even after adjustment. Dividing by $a_U$ corrects the expected scale of the unspliced count, but it cannot recover information that was not observed. When $a_U$ is small, each observed intronic read has greater influence on the adjusted estimate.

Worked Example: What a Zero Count Can Mean

Suppose the expected observed unspliced count for a cell is $a_Uu=3$. Under the Poisson model, the chance of observing no unspliced reads is $$ P(X_U=0\mid a_Uu=3)=e^{-3}=0.0498. $$

Thus a zero is possible even when the expected count is three. If coverage is lower, say $a_U=0.5$ with latent abundance $u=3$, the expected count is $1.5$, and $$ P(X_U=0)=e^{-1.5}=0.223. $$

In the second case, more than one in five repeated measurements would yield zero under the stated model. Calling that cell's true unspliced abundance zero would confuse an observed count with a latent quantity. The observation does provide information—it is less likely at high abundance—but a single zero is not proof of absence. Moreover, substituting $X_U=0$ into the plug-in variance estimate gives zero, which cannot sensibly mean that the latent abundance is known with certainty. It means the plug-in estimate is inadequate at this boundary; uncertainty should be assessed using a count likelihood or information pooled across suitable observations.

Worked Velocity Comparisons

Worked Example: A Difference Smaller Than Count Noise

Consider two hypothetical cells with equal exposures $a_U=a_S=1$. Let $\beta=0.5$ and $\gamma=0.2$, in compatible rate units. Their observed counts are $X_U=8, X_S=20$ in cell A and $X_U=16, X_S=20$ in cell B. Treat the kinetic rates and exposures as fixed for this illustration.

The estimated velocities are $$ \hat v_{s,A}=0.5(8)-0.2(20)=0,\qquad \hat v_{s,B}=0.5(16)-0.2(20)=4. $$

Using the Poisson plug-in variances, cell A has $$ \widehat{\operatorname{Var}}(\hat v_{s,A})=0.5^2(8)+0.2^2(20)=2.8, $$ so its standard error is $\sqrt{2.8}=1.673$. Cell B has variance $0.5^2(16)+0.2^2(20)=4.8$, and standard error $\sqrt{4.8}=2.191$. If the two measurement errors are independent, the estimated standard error of their difference is $$ \sqrt{2.8+4.8}=\sqrt{7.6}=2.757. $$

The observed difference is 4, or about $4/2.757=1.45$ estimated standard errors. This calculation does not establish a biological difference or constitute a formal test: the cells may differ biologically, and kinetic-rate uncertainty is omitted. It does show why a difference in point estimates should not automatically be read as a clear kinetic contrast. Here the count-noise scale is substantial relative to the observed difference.

Worked Example: Unequal Intronic Coverage Can Reverse a Sign

Suppose a hypothetical gene has latent abundances $u=10$ and $s=20$, with $\beta=0.4$ and $\gamma=0.1$. Its latent spliced velocity is $$ v_s=0.4(10)-0.1(20)=2. $$

Assume spliced exposure $a_S=1$, but lower intronic coverage gives $a_U=0.25$. A possible observed pair is $X_U=2$, $X_S=20$. Comparing the raw counts as if they had the same exposure gives $0.4(2)-0.1(20)=-1.2$, the opposite sign from the latent velocity.

Adjusting for exposure gives $\hat u=2/0.25=8$, $\hat s=20/1=20$, and $$ \hat v_s=0.4(8)-0.1(20)=1.2. $$

The adjustment corrects the expected scale, but this single estimate remains noisy. Its plug-in variance is $$ 0.4^2\frac{2}{0.25^2}+0.1^2\frac{20}{1^2} =0.16(32)+0.2=5.32, $$ giving a standard error of $\sqrt{5.32}=2.307$. The estimate of 1.2 is compatible with substantial measurement variation around the true value of 2. In fact, under the assumed exposures the expected observed intronic count is $0.25(10)=2.5$; observing 2 is ordinary sampling variation. Exposure correction prevents a systematic scale mismatch, but it does not eliminate count noise.

Diagnostics, Alternatives, and Limits

Several checks help distinguish a plausible kinetic pattern from a pattern dominated by measurement:

Inspect counts and exposures together. A zero with very low exposure is weak evidence for low abundance. Compare coverage across cells and between unspliced and spliced counts.
Compare signal with its uncertainty. A large-looking velocity based on one or two reads may still have a large standard error. Report or inspect uncertainty alongside the point estimate.
Check for variance beyond the Poisson model. If replicate counts vary more than their means predict, unmodeled heterogeneity, technical variation, or other sources of overdispersion may be present. A Poisson standard error may then be too small; a model allowing variance greater than the mean or an empirical variance estimate may be more appropriate.
Examine whether a result depends on a few cells. Strong changes in the estimated direction when individual low-count cells are removed suggest that the signal is fragile.
Separate measurement problems from kinetic misfit. Noise can make an estimate imprecise; systematic disagreement between the kinetic model and the observed pattern is a different problem, addressed by the earlier material on diagnosing kinetic model misfit.

The main alternatives have different trade-offs. Using raw counts is simple but can confound abundance with depth and unequal coverage. Exposure adjustment targets the correct expected scale, but low exposure still produces high variance and uncertain exposures can introduce bias. Pooling counts across cells can stabilize estimates when the cells share a meaningful state, but it can conceal genuine cell-to-cell changes and blur transitions. Smoothing or borrowing information across cells can also reduce random variation, yet may create an overly coherent pattern if the chosen cells are not comparable. None of these approaches turns sparse measurements into direct observations of latent abundance.

The Poisson assumption is a useful starting point, not a universal description. Ambient RNA, variable capture, dependence between count classes, and biological heterogeneity can all make the actual measurement process more complicated. In addition, the exposure-adjusted equations assume that the two RNA classes are measured on a comparable abundance scale and that the relative coverage is known. If $a_U$ is unknown, lower unspliced abundance and lower intronic detection can be difficult to distinguish from counts alone. An apparent velocity direction should therefore be treated cautiously when it depends on an uncertain coverage correction.

Key takeaway: A count is a noisy observation of abundance, not abundance itself. Account for exposure before comparing unspliced with spliced counts, and judge velocity estimates against the uncertainty that remains after adjustment.

Check Your Understanding

Try each question, then compare your reasoning with the brief solutions.

If $X_U\mid u\sim\operatorname{Poisson}(a_Uu)$, what is the variance of $X_U/a_U$ conditional on $u$?
With $a_U=0.2$, $X_U=3$, $\beta=0.5$, $a_S=1$, $X_S=10$, and $\gamma=0.1$, calculate the exposure-adjusted velocity estimate.
Why can the plug-in Poisson variance estimate be zero when the observed count is zero without implying certainty about abundance?
What is one limitation of pooling counts across cells to reduce sampling noise?

Solutions. (1) The variance is $u/a_U$, since $\operatorname{Var}(X_U/a_U)=a_Uu/a_U^2$. (2) $\hat u=3/0.2=15$ and $\hat s=10/1=10$, so $\hat v_s=0.5(15)-0.1(10)=6.5$ in the model's velocity units. (3) The plug-in substitutes the observed zero for the unknown Poisson mean; that boundary estimate does not rule out positive abundance, since a positive mean can still produce a zero count. (4) Pooling may combine cells with genuinely different states and obscure real cell-to-cell or transitional variation.