Tutorials › Statistical Methods for Single-Cell RNA Sequencing › Distinguishing Cell-State Uncertainty From Fate Probability

Branching and Lineage Models · Tutorial 427 of 700

Distinguishing Cell-State Uncertainty From Fate Probability

Distinguish what a fitted model cannot confidently assign from the biological chance of a cell reaching a particular fate.

Advanced 10 min read

What You'll Learn

  • Define posterior branch-assignment probabilities and biological fate probabilities as distinct quantities.
  • Identify when a branch label describes a cell’s current trajectory state rather than its terminal outcome.
  • Work through ambiguous and transitional-cell examples with explicit probability calculations.
  • Use competing transition hazards to derive eventual fate probabilities under stated assumptions.
  • Diagnose when a soft assignment is being misread as a fate prediction.
  • Compare hard labels, posterior assignments, and fate-probability estimates.

Two Questions That Can Have Different Answers

A cell near a branchpoint may be difficult to place on either branch from its measured expression. Separately, a cell in a transitional state may have a genuine chance of reaching either terminal fate. These are different questions: how uncertain is the model’s assignment of the observed cell? and what is the probability of each biological outcome, given the cell’s state and context? Confusing them can turn uncertainty about measurement or model fit into an unsupported claim about cell fate.

The distinction depends partly on what a model’s branch label means. A latent label may identify a cell’s current trajectory segment, or it may represent the cell’s eventual terminal lineage. Those definitions are not interchangeable. The interpretation of any reported probability therefore begins by naming the random quantity and the information on which the probability is conditioned.

Definition: Let \(Y_i\) denote the measured expression features for cell \(i\), \(D\) the data used to fit a branching model \(M\), and \(Z_i\) a latent branch or trajectory label. The posterior branch-assignment probability for branch \(b\) is \(P(Z_i=b\mid Y_i,D,M)\): the model’s conditional probability that the cell has label \(b\), given the observed data and model. Let \(F_i\) denote the cell’s eventual biological fate, \(S_i\) its current state, and \(E_i\) its specified biological context. A fate probability is \(P(F_i=b\mid S_i,E_i)\), under a stated biological or statistical model. If fate probabilities depend on unknown parameters \(\theta\), uncertainty about them is represented by a posterior distribution for \(\theta\); the posterior predictive fate probability averages over that uncertainty.

A posterior assignment is conditional inference about a latent label. Its uncertainty can arise because the cell’s measured features are compatible with multiple branches, because model parameters are uncertain, or because the model itself is only an approximation. A fate probability describes the distribution of a specified outcome under specified conditions. It may reflect stochastic biological processes, or it may be an operational conditional probability for a population of comparable cells. In either case, it is not automatically equal to uncertainty about a latent label.

The earlier material on soft lineage assignment and Bayesian inference for lineage fate probabilities provides the relevant modeling ideas: branch labels can be latent, and fate probabilities can be estimated probabilistically. Here the essential additional discipline is to state whether the latent label means current state or eventual fate, and to keep assignment and outcome probabilities distinct unless the model explicitly equates them.

Define the Random Quantity Before Interpreting a Probability

Suppose a fitted model assigns a cell probabilities of 0.55 to branch A and 0.45 to branch B. This statement is incomplete without the label definition. If \(Z_i\) denotes the cell’s current trajectory segment, the probabilities express uncertainty about where the cell belongs in the fitted representation. They do not say that the cell has a 55% biological chance of ending in fate A. If \(Z_i\) instead explicitly denotes the eventual fate \(F_i\), and the model is calibrated for the relevant state and context, the assignment probability can also serve as a model-based fate prediction. That interpretation requires those assumptions; it does not follow from the numbers alone.

In Bayesian notation, a fate probability may depend on unknown model parameters. If \(q_{ib}(\theta)=P(F_i=b\mid S_i,E_i,\theta)\), then the probability of fate \(b\) conditional on the available data is the posterior predictive probability:

$$ P(F_i=b\mid S_i,E_i,D) = \int P(F_i=b\mid S_i,E_i,\theta)\,p(\theta\mid D)\,d\theta $$

This averaging gives a predictive probability for the outcome while accounting for posterior uncertainty in \(\theta\). The posterior distribution \(p(\theta\mid D)\) and the predictive probability are related but not identical: the first describes uncertainty about model parameters, while the second describes the model’s probability for a future fate outcome after that uncertainty is integrated out. A posterior probability for \(Z_i=b\) answers a different question unless \(Z_i\) has been defined as that outcome and the relevant conditioning is the same.

Worked Example: An Ambiguous Cell With a Separate Fate Estimate

Suppose a two-branch model gives a cell equal prior probabilities of belonging to branches A and B. For the cell’s observed expression pattern \(Y_i\), the model assigns likelihood 0.30 under A and 0.28 under B. Calculate the posterior branch-assignment probabilities. Separately, suppose a validated model for cells in the same specified state and context estimates a 0.72 probability of eventual fate A.

For branch assignment, Bayes’ rule gives:

$$ P(Z_i=A\mid Y_i,D,M) = \frac{P(Y_i\mid Z_i=A,M)P(Z_i=A)} {P(Y_i\mid Z_i=A,M)P(Z_i=A)+P(Y_i\mid Z_i=B,M)P(Z_i=B)} $$

Substituting the equal prior probabilities:

$$ P(Z_i=A\mid Y_i,D,M) = \frac{0.30(0.50)}{0.30(0.50)+0.28(0.50)} = \frac{0.15}{0.29} = 0.5172 $$

The corresponding probability for B is \(1-0.5172=0.4828\). The expression pattern barely distinguishes the branch labels under this model. The separate fate estimate is 0.72 for A, conditional on the specified state and context. It does not contradict the assignment result: one concerns uncertainty about a model label from this cell’s expression, while the other is a fate prediction from a model intended to predict outcomes. The 0.72 estimate is not a guarantee that this particular cell will reach A, and it should not be substituted for the assignment probability.

Transitional Cells and Competing Outcomes

A transitional cell makes the distinction especially important. It may be confidently located in a shared trunk or intermediate state, even though its eventual destination is not yet determined. Thus, high confidence about a cell’s present position does not imply high confidence about its future fate. Conversely, a cell can be hard to assign to a current branch while its fate probability, under a separate validated model, is quite unequal between outcomes.

One simple fate model treats transitions to two terminal fates as competing events. Assume that, conditional on the current state and context, the instantaneous transition rates are constant: \(\lambda_A\) per day to fate A and \(\lambda_B\) per day to fate B. Assume also that these are the only possible exits, that neither rate changes over time, and that the cell remains at risk until one event occurs. The probability of no transition by time \(t\) is \(e^{-(\lambda_A+\lambda_B)t}\). The probability of transitioning to A in a short interval after time \(t\) is approximately \(\lambda_A\,dt\), conditional on still being in the state. Integrating over all possible transition times yields the eventual probability of A:

$$ P(F_i=A\mid S_i,E_i) = \int_0^\infty \lambda_A e^{-(\lambda_A+\lambda_B)t}\,dt = \frac{\lambda_A}{\lambda_A+\lambda_B} $$

The final equality follows because the integral of \(e^{-ct}\) from zero to infinity is \(1/c\) for \(c>0\). The probability for B is correspondingly \(\lambda_B/(\lambda_A+\lambda_B)\). These are fate probabilities under the stated constant-rate, two-outcome model, not posterior probabilities that a cell currently belongs to a branch.

Worked Example: A Transitional Cell With Unequal Exit Rates

A transitional state has estimated constant transition rates of \(\lambda_A=0.06\) per day and \(\lambda_B=0.04\) per day. An expression-based model assigns a cell a 0.90 posterior probability of being in the shared trunk state. Find the model-based eventual fate probabilities and interpret both results.

Under the two-outcome constant-rate assumptions, the eventual probability of fate A is:

$$ P(F_i=A\mid S_i,E_i) = \frac{0.06}{0.06+0.04} = \frac{0.06}{0.10} = 0.60 $$

The probability of fate B is \(0.04/0.10=0.40\). Separately, the posterior assignment probability of 0.90 concerns the cell’s current label, the shared trunk. It does not mean that the cell has a 90% chance of remaining in the trunk or reaching either terminal fate. Under the specified transition model, A is more likely than B, but B remains a substantial possibility. Both statements are conditional on their respective models and assumptions.

Derivation: Why Better Assignment Need Not Change Fate Propensity

The source of assignment uncertainty can be seen directly from Bayes’ rule. For two possible labels A and B, posterior odds equal prior odds multiplied by the likelihood ratio:

$$ \frac{P(Z_i=A\mid Y_i,D,M)}{P(Z_i=B\mid Y_i,D,M)} = \frac{P(Y_i\mid Z_i=A,M)}{P(Y_i\mid Z_i=B,M)} \cdot \frac{P(Z_i=A)}{P(Z_i=B)} $$

If the observed features have nearly equal likelihood under A and B, the likelihood ratio is near one. The data then provide little information for distinguishing the labels, so the posterior assignment remains close to the prior odds. This is a statement about the information in the observed features under the model. It is not a derivation of the cell’s biological transition propensity.

A fate probability changes when information relevant to the fate model changes, or when the fate model’s parameters are updated. Improving measurement of a current label may sharpen a posterior assignment without changing the biological state, context, or fate model. If the improved measurement also provides validated information about future fate, then the fate prediction should be updated too. Whether it does so is a modeling question, not something implied by a more confident branch assignment.

Worked Example: Sharper Assignment, Same Stated Fate Probability

A current-state classifier begins with equal prior probabilities for labels A and B. A new expression feature has likelihood 0.90 under label A and 0.10 under label B. Assume the feature is used only to identify the current label, and does not alter the separately specified fate model, which gives a 0.65 probability of eventual fate A for the cell’s fixed state and context.

The posterior assignment probability is:

$$ P(Z_i=A\mid Y_i) = \frac{0.90(0.50)}{0.90(0.50)+0.10(0.50)} = \frac{0.45}{0.50} = 0.90 $$

The classifier’s assignment is now much more certain than it was under an uninformative feature, which would have left the probability at 0.50. By the example’s explicit assumption, the fate model remains at 0.65 for A and 0.35 for B. The result illustrates that certainty about a current label can change without changing the specified fate propensity. If the feature were known to predict fate, the assumption of no effect on the fate model would be inappropriate and the fate prediction would need to be reconsidered.

Diagnostics, Alternatives, and Limitations

Before reporting a probability, check the label definition, the conditioning information, and the model’s biological interpretation. A posterior predictive check can reveal that a branching model fails to reproduce relevant features of the data, but satisfactory checks do not by themselves prove that its latent labels equal biological outcomes. Evidence about model fit and evidence about fate semantics are separate requirements.

  • Check what the label represents. If it denotes a current trajectory segment, report it as a segment-assignment probability, not as the probability of terminal fate.
  • Check what information conditions the probability. State whether the quantity is conditional on expression, current state, context, fitted parameters, or an outcome model. Probabilities with different conditioning information need not agree.
  • Check calibration and outcome definitions. A fate probability needs a defined endpoint and a model supported by relevant outcome information. A soft clustering output alone does not establish calibration for future fate.
  • Check assumptions about transitions. Constant competing rates are convenient but may fail if rates vary with time, cells can enter a third fate, or the model omits cell death or another exit. Then the two-rate formula is not justified.
  • Inspect ambiguity rather than hiding it. If several labels have similar posterior probabilities, a hard label conceals uncertainty. If one label dominates, that still does not establish a biological fate probability unless the label and model support that interpretation.

Hard assignment replaces a probability vector with the most probable label. It is simple to display, but discards uncertainty and can make borderline cells look certain. Soft assignment preserves the posterior distribution over labels, but it remains a statement about those labels as defined by the model. Direct fate prediction instead targets eventual outcomes and should be based on an explicitly defined outcome model; a posterior predictive probability can account for uncertainty in that model’s parameters. None of these approaches removes the need to justify the meaning of the target quantity.

The principal limitation is that observational expression data may not identify either a unique latent branch assignment or an individual cell’s eventual fate. A high posterior probability can still be conditional on a misspecified model, while an estimated fate probability can be poorly calibrated or apply only to a different state or context. Finally, probability does not imply a fractional biological outcome: the cell either realizes a particular outcome under the defined endpoint or it does not. The probability describes uncertainty or outcome frequency under stated conditions, not a literal fraction of the cell.

Key takeaway: Always name the event before interpreting its probability. \(P(Z_i=b\mid Y_i,D,M)\) is a posterior assignment probability for a model-defined label; \(P(F_i=b\mid S_i,E_i)\) is a probability for an eventual fate under specified conditions. They coincide only when the model’s label, outcome, and conditioning are explicitly aligned and the model is appropriate.

Check Your Understanding

For each question, identify what the probability refers to before interpreting it.

  1. A model assigns a cell probabilities 0.80 to the trunk and 0.20 to branch A. What additional definition is needed before interpreting 0.20 as a fate probability?
  2. Under a two-fate constant-rate model, \(\lambda_A=0.03\) per day and \(\lambda_B=0.02\) per day. What are the eventual fate probabilities?
  3. A cell is assigned to its current branch with posterior probability 0.95. Does that imply a 95% probability of reaching that branch’s terminal fate? Explain.
  4. A new measurement makes a current-label posterior more concentrated. Give one condition under which the fate probability should also be updated.

Solutions. (1) Establish whether the label denotes the eventual fate rather than the cell’s current state, and whether the model is valid and conditioned on the information relevant to that fate. Without this, 0.20 is only a posterior probability for the model-defined label. (2) Under the stated two-outcome, constant-rate assumptions, \(P(F_i=A)=0.03/(0.03+0.02)=0.60\) and \(P(F_i=B)=0.02/0.05=0.40\). (3) No. The 0.95 concerns assignment to the current branch; a terminal-fate probability requires a model connecting current state and context to eventual outcomes. (4) Update the fate probability if the measurement carries validated information about the eventual outcome, or changes relevant information or parameters in the fate model. A sharper estimate of a current label alone does not require a changed fate probability.

Clinical Biostats