Tutorials › Statistical Methods for Single-Cell RNA Sequencing › Reporting Cell-Type Labels With Uncertainty

Cell-Type Annotation and Its Uncertainty · Tutorial 352 of 700

Reporting Cell-Type Labels With Uncertainty

Report cell-type annotations so readers can see what was assigned, how uncertain it was, which reference informed it, and why labels were reported at a particular level of detail.

Intermediate 12 min read

What You'll Learn

  • Distinguish assigned labels, model scores, estimated probabilities, and prediction sets.
  • Report abstentions and label proportions with explicit denominators.
  • State the reference, annotation method, and relevant provenance.
  • Document reproducible rules for choosing coarse or fine labels.
  • Recognize summaries that can conceal uncertainty or dependence among cells.

Why Reporting Is Part of Annotation

A cell-type annotation is not just a name attached to a barcode. A useful report shows what label was assigned, what evidence or uncertainty accompanies it, how unresolved cells were handled, and which reference and decision rules shaped the result. Without those details, readers cannot tell whether a confident-looking label came from strong evidence, a forced choice among unsuitable alternatives, or a decision to report only a broad parent category.

This tutorial focuses on communicating annotation results rather than choosing an annotation method. The earlier material on selective classification, prediction sets, reference mismatch, and adjudicating disagreements provides ways to assess or resolve uncertainty. Here the task is to make the resulting choices interpretable and auditable. In particular, a reported score should not be called a probability unless its meaning and calibration support that description.

Definition: For cell \(i\), let \(\hat y_i\) denote the reported label, \(q_{ik}\) a reported score for candidate label \(k\), \(A_i\) a reported set of plausible labels when the method returns a set, and \(R_i=1\) indicate that the workflow abstained from assigning a cell-type label. These are distinct report fields: an assigned label is a decision, a score summarizes method output, a set records multiple candidates, and an abstention records that no label was assigned under the stated rule.

A report should define the interpretation of every uncertainty quantity. For example, \(q_{ik}\) might be a normalized classifier output, a similarity score, or a calibrated estimate of the probability of label \(k\) conditional on the model and its inputs. These quantities are not interchangeable. A normalized vector sums to one, but that alone does not establish that its entries are calibrated probabilities. If the method returns scores rather than calibrated probabilities, call them scores and state how they were produced. If it returns calibrated probabilities, state the calibration procedure and the validation setting relevant to interpreting them.

What a Reproducible Report Contains

At minimum, report the unit being annotated, the label vocabulary and its hierarchy, the method and reference used, the rule that converts method output into reported labels, and the handling of unresolved cases. For a cell-level data release, these details can be paired with a table containing one row per cell. For a manuscript or summary, provide aggregate counts as well, while making clear which cells form each denominator.

  • Label: Give the assigned label, or an explicit value such as “abstain” or “unresolved.” Do not silently drop cells without a label.
  • Uncertainty: Specify whether the field is a score, calibrated probability, prediction set, or other summary. State the candidate labels and any threshold or decision rule used.
  • Reference provenance: Identify the reference dataset or atlas, its version or release, relevant tissue and species scope, and any filtering or mapping that affects the candidate labels.
  • Granularity: State the hierarchy and the rule for reporting a parent category instead of a child label. Explain whether labels were assigned directly at that level or rolled up from finer predictions.
  • Scope: Report cell counts, the number of abstentions, and the denominator for each proportion. If relevant, distinguish a cell-level description from population-level inference.

Provenance should be detailed enough that another analyst can identify the reference and reconstruct the relevant label vocabulary. “Annotated using a public atlas” is not enough if multiple releases, tissues, or label mappings exist. Conversely, provenance should not imply that the reference is a ground truth for every query cell: a reference defines the comparison context, not an independent guarantee that every possible cell type is represented or correctly labeled.

Worked Example: Reporting Coverage and Label Proportions

Suppose a hypothetical dataset contains 100 cells. A stated rule assigns 82 cells a label and abstains on 18. Among the assigned cells, 50 are labeled “T cell,” 20 “B cell,” and 12 “myeloid cell.” Show both the label proportions among assigned cells and the proportions among all cells.

Let \(n=100\) be the total number of cells and \(n_A=82\) the number assigned a label. Coverage is the fraction receiving a label:

$$ \widehat{\text{coverage}}=\frac{n_A}{n}=\frac{82}{100}=0.82. $$

Thus, 82% of cells received a label under the stated rule, and 18% were abstentions. Among assigned cells, the T-cell proportion is \(50/82=0.6098\), the B-cell proportion is \(20/82=0.2439\), and the myeloid-cell proportion is \(12/82=0.1463\). These conditional proportions sum to one, up to rounding.

Among all 100 cells, the corresponding labeled proportions are \(50/100=0.50\), \(20/100=0.20\), and \(12/100=0.12\), with an abstention proportion of \(18/100=0.18\). The identity

$$ \frac{n_k}{n} = \frac{n_A}{n}\times\frac{n_k}{n_A} $$

shows why the two denominators answer different questions: the all-cell proportion for a label equals coverage multiplied by that label’s proportion among assigned cells. A report that gives only the assigned-cell proportions describes the composition of the labeled subset, not the full dataset. A suitable summary therefore states both the denominator and the abstention count.

Reporting Uncertainty Without Overstatement

The form of the uncertainty report should match the method output. For a classifier that returns a score for each candidate label, a cell-level record might include the assigned label and the leading scores, along with a note that they are uncalibrated scores if that is their status. For a method returning a prediction set \(A_i\), report the set or a clearly defined summary of its size; do not present one member as uniquely established unless a separate decision rule selects it. If an abstention rule is used, report its criterion and the number of cells affected.

A compact report can summarize distributions of uncertainty, but it should not replace the cell-level information needed to inspect ambiguous cases. For example, report the median and range of a defined uncertainty score across cells, or the fraction of cells whose prediction set contains more than one label. State whether those summaries include abstentions and how missing or undefined scores were handled. A mean score alone can conceal a small group of highly uncertain cells, so useful summaries often include a distribution or thresholded counts.

Where probability estimates are reported, the interpretation is conditional on the model, reference, and validation context. A probability of 0.8 should not automatically be paraphrased as “there is an 80% chance this cell is truly this type.” That reading is justified only under an appropriate probabilistic interpretation and calibration evidence for a sufficiently relevant setting. Even calibrated performance in one validation population may not carry over unchanged under tissue, donor, batch, or reference distribution differences.

Worked Example: Writing a Cell-Level Uncertainty Record

Consider a hypothetical cell for which a classifier reports normalized outputs of 0.72 for “memory B cell,” 0.23 for “naive B cell,” and 0.05 for “T cell.” Suppose the workflow has not established calibration, and its reporting rule assigns the highest-scoring label only when the leading score is at least 0.70. What should be reported?

The largest output is \(0.72\), so the rule assigns “memory B cell.” The cell-level record should retain the assigned label, the relevant scores, and the threshold rule. Because calibration has not been established, the values should be described as normalized model scores, not as probabilities of true cell identity. For example: “Assigned label: memory B cell; normalized classifier score: 0.72; next score: 0.23; rule: assign the highest-scoring label when its score is at least 0.70.”

For contrast, if the leading score were 0.68 with the same rule, the cell would not receive a label because \(0.68<0.70\). The record should say “abstain under the 0.70 threshold,” rather than assigning “memory B cell” simply because it ranked first. The threshold decision is reproducible, but it does not by itself make the score a probability or establish that the label is biologically correct.

Making Label Granularity Explicit

Cell-type labels often form a hierarchy: a broad category can have several more specific child labels. Reporting at fine granularity can be informative when evidence supports it, but it can imply distinctions the data or reference do not resolve. Reporting only broad categories can avoid unsupported specificity while hiding distinctions that are well supported. A transparent report states the available hierarchy and the rule that determines the level shown.

A granularity rule might require a minimum score, a minimum margin between the top two child labels, or a specified validation criterion before reporting a child. These rules should be fixed and described rather than chosen after inspecting which labels produce the preferred narrative. If a child-level assignment fails the criterion, a workflow may report its parent if parent-level evidence supports that category; otherwise it may abstain. The report should distinguish a parent label selected directly from evidence from a parent label produced by rolling up child predictions.

Worked Example: Applying a Parent-versus-Child Rule

Suppose a hypothetical label hierarchy has “B cell” as a parent and “naive B cell” and “memory B cell” as children. The reporting protocol assigns a child label only if its score is at least 0.75 and exceeds the other child’s score by at least 0.20. If neither child qualifies, it reports “B cell” only when the combined B-cell score is at least 0.80; otherwise it abstains. For one cell, the child scores are 0.54 for naive B cell and 0.31 for memory B cell, and the combined B-cell score is 0.85.

The leading child score is 0.54, which is below 0.75. Its margin over the other child is \(0.54-0.31=0.23\), which meets the margin condition but cannot compensate for failing the score threshold. Therefore, neither child label is reported. The combined B-cell score is \(0.85\), and \(0.85\geq0.80\), so the rule reports the parent label “B cell.” The report should identify this as a parent-level assignment under the stated rule, not as a resolved naive-versus-memory distinction.

If the combined score had instead been 0.76, it would fail the parent threshold as well; the rule would abstain. Showing the hierarchy and thresholds makes these outcomes auditable and prevents broad labels from being mistaken for fine-grained evidence.

Derivation: Why Denominators Must Be Shown

The distinction between all-cell proportions and assigned-cell proportions follows directly from the counts. Let \(n_k\) be the number of cells assigned label \(k\), and let \(n_A\) be the total number of labeled cells. When \(n_A>0\), the assigned-cell proportion is \(n_k/n_A\), while coverage is \(n_A/n\). Multiplying gives:

$$ \frac{n_A}{n}\frac{n_k}{n_A}=\frac{n_k}{n}. $$

The cancellation is valid because \(n_A>0\). This identity does not make the two proportions interchangeable: \(n_k/n_A\) conditions on receiving a label, whereas \(n_k/n\) describes the share of the entire dataset assigned that label. If no cells are labeled, the assigned-cell proportions are undefined; the report should give zero coverage and the abstention count, not a proportion with a zero denominator.

Alternatives, Diagnostics, and Limitations

Reporting formWhat it communicatesWhat it can conceal
One label per cellA concise, usable annotation decisionCompeting candidates and the degree of uncertainty
Label plus scores or probabilitiesRelative model outputs for candidate labelsScore meaning or calibration, if not documented
Prediction setMultiple labels retained as plausible by a defined procedureWhich label is preferred, unless a selection rule is also reported
Abstention indicatorCells for which the workflow withheld a label under a ruleWhy the evidence was insufficient, unless the rule and reason are documented

These forms can be combined: a report may give an assigned label for some cells, prediction sets for others, and abstentions where the protocol calls for them. A single-label table is easy to use downstream, but should not be the only record when uncertainty matters. Conversely, a large collection of scores is not automatically informative if its scale, interpretation, and decision rule are absent.

Several checks help detect reporting problems. Confirm that category counts plus abstentions equal the total number of cells; that proportions use stated denominators; and that per-cell labels match the declared thresholds and hierarchy. Check for missing reference versions, undocumented label mappings, and labels finer than the reference or evidence can support. If uncertainty summaries are unexpectedly uniform, verify that the reported quantity was not inadvertently reduced to the chosen label alone. Compare summaries across relevant batches or donors when those groupings could reveal systematic differences.

Cell-level counts are descriptive of the analyzed dataset. Cells from the same donor are not generally independent biological replicates, so treating thousands of cells as thousands of independent replicates can overstate precision in population-level claims. If uncertainty intervals or inferential comparisons are added, their sampling and resampling units must reflect the study design, often donors rather than individual cells. Annotation reporting alone does not establish that a cell type is present in a broader population or that an observed difference is causal.

Common Reporting Mistakes

  • Calling every normalized output a probability. Values summing to one do not guarantee calibration. Name the quantity according to what the method actually provides.
  • Hiding abstentions. Removing unlabeled cells changes the denominator and can make coverage appear complete. Report their count and the rule that produced them.
  • Reporting a broad label as a fine-grained result. A parent assignment does not resolve which child label applies. State the level and how it was chosen.
  • Omitting reference provenance. A label depends partly on the candidate vocabulary and reference context. Give enough identifying information to reproduce that context.
  • Using proportions without denominators. “Twenty percent B cells” is ambiguous if it is unclear whether abstained cells are included.

The practical standard is simple: let readers distinguish the assigned decision from the supporting uncertainty information, and let them reconstruct how the reference, thresholds, abstentions, and hierarchy shaped that decision.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions below.

  1. In a dataset of 250 cells, 200 receive labels. What is coverage, and what denominator should be used for the proportion of one label among assigned cells?
  2. A classifier’s normalized outputs sum to one. Is that sufficient to call them calibrated probabilities? Explain.
  3. A cell fails a child-label threshold but passes the stated parent-label threshold. What should the report say about its granularity?
  4. Why should a report identify the reference version and label mapping?
  5. What is wrong with reporting a label proportion after silently removing all abstentions?

Solutions. (1) Coverage is \(200/250=0.80\), or 80%. For a label proportion among assigned cells, use 200 as the denominator and state that choice. (2) No. Summing to one describes normalization, not calibration; report them as scores unless calibration and the probability interpretation are supported. (3) Report the parent label under the stated rule and make clear that the child distinction was not resolved. (4) Reference releases and mappings can change the candidate vocabulary and label interpretation, so provenance is needed to interpret and reproduce annotations. (5) The proportion then describes only the labeled subset, may differ from the all-cell proportion, and hides coverage; report abstentions and the denominator explicitly.

Clinical Biostats