Why Captured Proportions Need Not Match Source Proportions
A cell-type proportion calculated from a single-cell assay describes the cells that entered and remained in the analyzed data, not automatically the cells present in the source sample. Cell types can differ in dissociation sensitivity, viability, sorting, loading, or recovery. If some types are less likely to be captured, the observed composition is shifted even when the source composition is known and unchanged. Statistical adjustment for sampling variation does not, by itself, remove this measurement distortion.
The distinction matters for differential-abundance comparisons. A difference in the captured proportions between conditions may reflect a biological change in source composition, a condition-related change in capture or survival, or both. The data alone do not generally separate those explanations. The appropriate target and interpretation therefore depend on whether capture probabilities are known, stable, or plausibly different across conditions.
A Model for Capture-Weighted Proportions
This model assumes that each source cell has a type, that capture and retention can be represented by a type- and condition-specific probability, and that the analyzed labels correctly identify the captured types. It also treats the source composition and capture probabilities as fixed within each condition for the comparison. Conditional on the total number \(n_g\) of captured cells, a basic sampling model is \(Y_{g1},\ldots,Y_{gK}\sim\operatorname{Multinomial}(n_g;r_{g1},\ldots,r_{gK})\), where \(Y_{gk}\) is the captured count of type \(k\). This multinomial model describes sampling variation around the captured composition; it does not assert that \(r_{gk}=\pi_{gk}\).
The expression for \(r_{gk}\) follows by considering the probability that a randomly selected source cell is both type \(k\) and captured. That joint probability is \(\pi_{gk}q_{gk}\). The probability that a source cell is captured, regardless of type, is \(\sum_j\pi_{gj}q_{gj}\). Conditioning the joint probability on capture gives the observed proportion:
The denominator is essential: observed proportions sum to one, so a low capture probability for one type changes the observed proportions of the other types as well. Thus capture bias is compositional, not merely a type-specific downward correction.
Worked Example: Unequal Capture Changes the Observed Composition
Suppose a hypothetical source sample contains 40% type A and 60% type B. Assume capture probabilities are 0.50 for type A and 0.90 for type B, with no other types.
The unnormalized captured contributions are \(0.40(0.50)=0.20\) for A and \(0.60(0.90)=0.54\) for B. The total capture probability is \(0.20+0.54=0.74\). Therefore,
Although type A makes up 40% of the source sample, it makes up about 27.0% of captured cells. Type B is correspondingly overrepresented at about 73.0%. These are expected captured proportions under the assumptions, not guaranteed percentages in a finite data set.
What Can Be Identified Without Calibration?
Observed proportions do not generally identify source proportions and capture probabilities separately. For each type, the observed data depend on the product \(\pi_{gk}q_{gk}\), followed by normalization. Different combinations of source composition and capture probabilities can therefore yield the same captured composition. More cells can reduce uncertainty about \(r_{gk}\), but cannot resolve this structural ambiguity.
For two types, the relationship is especially clear in odds form:
The observed odds combine source odds with relative capture probability. Unless the relative capture probability is known or constrained, the source odds cannot be recovered from the observed odds alone.
Worked Example: The Same Observed Data Allow Different Source Compositions
Suppose a hypothetical assay yields observed proportions \(r_A=r_B=0.50\). One explanation is equal source proportions and equal capture probabilities. Another explanation has \(q_A=0.25\) and \(q_B=0.75\), with source proportions \(\pi_A=0.75\) and \(\pi_B=0.25\).
Under the second explanation, the unnormalized captured contributions are \(0.75(0.25)=0.1875\) for A and \(0.25(0.75)=0.1875\) for B. Normalizing gives \(r_A=r_B=0.50\), just as under the first explanation. Thus the observed 50:50 composition cannot distinguish a 50:50 source with equal capture from a 75:25 source with unequal capture.
The example is not evidence that either source composition is correct. It demonstrates non-identifiability: captured proportions alone do not supply the missing calibration information.
The correction assumes that the calibration applies to the samples being analyzed, that the relevant cell types are represented in the calibration, and that capture probabilities are positive. If a type has zero probability of capture, its source proportion cannot be recovered from captured cells by inverse weighting. If the probabilities are estimated rather than known, uncertainty in those estimates should also be reflected in uncertainty about the corrected proportions.
Worked Example: Correcting With Calibrated Capture Probabilities
In a hypothetical sample, the observed proportions are \(r_A=0.30\) and \(r_B=0.70\). An external calibration applicable to this sample gives \(q_A=0.50\) and \(q_B=0.875\). Calculate the corrected source proportions.
Divide each observed proportion by its capture probability to obtain unnormalized values: \(0.30/0.50=0.60\) for A and \(0.70/0.875=0.80\) for B. Their sum is \(1.40\), so the normalized estimates are
Under the stated calibration and model, the corrected estimates are about 42.9% A and 57.1% B. The correction is not justified merely because these values seem more plausible; it depends on the calibration being relevant and sufficiently accurate.
Comparing Conditions: Stable and Changing Capture
Whether capture bias cancels in a comparison depends on the contrast and on how capture probabilities behave across conditions. If \(q_{gk}=q_k\) is the same in every condition, pairwise observed log-ratios have a useful property:
Comparing this log-ratio between two conditions cancels the fixed capture term. Therefore, the change in the pairwise log-ratio equals the change in the source pairwise log-ratio, under the model and stable relative capture. This does not mean that the observed percentages themselves equal the source percentages, or that every summary of them is unbiased. If capture probabilities vary by condition, an additional capture term remains in the contrast.
Worked Example: Stable Capture Preserves a Pairwise Log-Ratio Change
Consider two hypothetical conditions with source proportions \((\pi_A,\pi_B)=(0.40,0.60)\) and \((0.60,0.40)\). Suppose the capture probabilities are unchanged across conditions: \(q_A=0.50\) and \(q_B=0.90\).
In the first condition, the captured contributions are 0.20 and 0.54, giving \(r_A=0.20/0.74=0.2703\) and \(r_B=0.54/0.74=0.7297\). In the second condition, they are 0.30 and 0.36, giving \(r_A=0.30/0.66=0.4545\) and \(r_B=0.36/0.66=0.5455\).
The source odds for A versus B change from \(0.40/0.60=0.6667\) to \(0.60/0.40=1.500\), a ratio of odds of \(1.500/0.6667=2.250\). The observed odds change from \(0.2703/0.7297=0.3704\) to \(0.4545/0.5455=0.8333\), also a ratio of odds of \(0.8333/0.3704=2.250\). Hence the pairwise log-ratio change is preserved. But the observed A proportions, 27.0% and 45.5%, are not the source proportions, 40.0% and 60.0%; the percentage-point changes are not interchangeable.
Diagnostics, Alternatives, and Limits
Capture bias is difficult to diagnose from the captured composition alone, because that composition is the quantity jointly shaped by source proportions and capture. Useful evidence must therefore come from design information or independent measurements. Relevant checks include whether sample handling, dissociation, viability, sorting, and loading procedures were comparable across conditions; whether cell-type recovery differs in technical controls; and whether calibration material provides type-specific capture estimates that apply to the study samples. These checks can reveal concerns, but a lack of visible imbalance does not prove equal capture.
There are three practical approaches, each answering a different question:
| Approach | What it supports | Main limitation |
|---|---|---|
| Analyze observed proportions as captured | A comparison of cell composition among cells captured under the specified assay and processing workflow | Does not by itself estimate source composition when capture differs by type |
| Correct using calibrated relative capture probabilities | An estimate of source proportions under the capture model | Depends on transportable, sufficiently precise calibration and positive capture probabilities |
| Sensitivity analysis over plausible capture probabilities | An assessment of how conclusions change under specified capture scenarios | Does not identify which scenario is true without external evidence |
A differential-abundance test can quantify sampling uncertainty in observed proportions, but a small p-value does not establish that a difference is biological rather than induced by condition-specific capture. Similarly, a large number of captured cells cannot compensate for missing calibration when the target is the source composition. If capture probabilities are believed stable, that belief should be justified by the protocol and supporting evidence; it is an assumption, not a consequence of observing similar sample processing.
A common mistake is to reason, “Type A is rare among captured cells, so it must be rare in the tissue,” or, “Both conditions used the same assay, so capture bias cancels.” The first ignores type-specific capture. The second is too broad: stable capture can preserve certain pairwise log-ratio contrasts, but does not make observed proportions equal to source proportions; condition-specific differences in survival or recovery can also invalidate the cancellation. State the target explicitly and describe conclusions as applying to the captured assay composition when source-level calibration is unavailable.
Check Your Understanding
Answer each question, then compare with the solutions below.
- A source contains 30% type X and 70% type Y. Their capture probabilities are 0.80 and 0.40. What are the expected observed proportions?
- Why can increasing the number of captured cells not, by itself, identify source proportions when capture probabilities are unknown?
- For two types, what condition allows the capture term to cancel from a between-condition change in pairwise log-ratio?
- Observed proportions are 20% type X and 80% type Y. Calibration gives relative capture probabilities \(q_X=0.50\) and \(q_Y=1.00\). Compute the corrected proportions.
Solutions. (1) The unnormalized contributions are \(0.30(0.80)=0.24\) and \(0.70(0.40)=0.28\); the total is 0.52. Thus \(r_X=0.24/0.52=0.4615\) and \(r_Y=0.28/0.52=0.5385\). (2) More captured cells improve estimation of the captured proportions \(r_k\), but \(r_k\) still depends on both \(\pi_k\) and \(q_k\). Without calibration or assumptions, different source compositions and capture probabilities can produce the same \(r_k\). (3) Relative capture probabilities must be stable across the compared conditions, so the term \(\log(q_X/q_Y)\) is the same and cancels in the contrast. (4) The unnormalized corrected values are \(0.20/0.50=0.40\) and \(0.80/1.00=0.80\), with sum 1.20. The corrected proportions are \(0.40/1.20=0.3333\) for X and \(0.80/1.20=0.6667\) for Y, conditional on the calibration being applicable.