Single-Cell RNA-seq Quality Control: Cells, Genes, and Mitochondrial Reads

By Dr. Zubair Khalid, DVM, MS, PhD · 2026-07-12

Key Takeaways

Context-aware filtering of single-cell RNA sequencing (scRNA-seq) data is paramount, moving beyond universal thresholds to account for tissue of origin, dissociation methods, and specific biological questions.
Mitochondrial read fraction is a critical quality metric, but its interpretation is highly context-dependent; elevated levels (potentially exceeding 20%) can be biologically relevant in metabolically active or hypoxic tissues, rather than solely indicative of cell death.
Gene and UMI (Unique Molecular Identifier) counts per cell require careful evaluation using knee plots to identify inflection points, distinguishing between empty droplets and potentially doublets, with typical ranges varying significantly by protocol and sequencing depth.
Doublet detection is a distinct and essential step, as gene and UMI filters alone are insufficient to remove cells originating from two or more cells, which can confound downstream analyses.
Iterative filtering and validation are necessary, involving re-examination of QC metrics per cluster after initial analysis to ensure that biologically relevant cell populations are retained and not inadvertently discarded.
Transparent documentation of all filtering thresholds, justifications, software versions, and doublet removal strategies is crucial for ensuring the reproducibility of scRNA-seq analyses.

Single cell RNA seq quality control is not about applying universal thresholds. It is about making context aware filtering decisions based on your tissue of origin, dissociation method, and biological question. This guide is for bench scientists and bioinformaticians who want to filter cells and genes with confidence and document their choices transparently. You will learn why mitochondrial read percentage varies across cell types, how to set gene and UMI cutoffs that match your data, and how to report your decisions so that others can reproduce your analysis.

At a Glance

Metric	Typical Filtering Range	Context Sensitive Considerations
Number of genes per cell	200 to 5000 detected genes	Neurons and stem cells often express more genes. Fibroblasts may express fewer. Check your UMI saturation curve.
Number of UMIs per cell	500 to 20000 unique molecules	Low UMI cells may be empty droplets or damaged. High UMI cells could be doublets. Vary by protocol and sequencing depth.
Mitochondrial read fraction	5% to 20% maximum	High in stressed or dying cells. Hypoxic tissues such as keloid (source) or COPD lung (source) naturally have higher mito reads.
Doublet rate	0.5% to 8% per 1000 cells	Depends on loading concentration. Use software like DoubletFinder. Do not rely on a single QC filter to remove doublets.

Decision Criteria

Why Universal Thresholds Fail

Standard cutoffs such as “keep cells with fewer than 20% mitochondrial reads” appear in many tutorials from sources like the Galaxy Training Network. These numbers work for common cultured cell lines. They break for tissues with high metabolic activity. For example, a study on human ovarian aging (source) found that granulosa cells naturally express more mitochondrial transcripts because they support follicular development. Applying a 10% mito cutoff would discard healthy cells in that context.

Tissue Specific Considerations

Your filtering thresholds must reflect the biology of your sample. Use these questions to decide:

Is your tissue hypoxic or metabolically active? Sources and describe conditions where mitochondrial reads are elevated even in viable cells.
Did you use enzymatic dissociation? Harsh dissociation increases stress and mito reads. Compare with single nucleus protocols such as the minimized bias profiling in source.
What is your expected cell size and complexity? Immune cells are small and have fewer transcripts. Epithelial cells are larger and may have more genes.

Gene and UMI Cutoffs

Do not set gene cutoffs by guesswork. Plot the number of genes per cell against the number of UMIs per cell. Cells with very low gene counts are likely empty droplets. Cells with very high gene counts may be doublets. Use the knee plot of cumulative UMI counts to identify the inflection point. The NCBI Bookshelf provides background on UMI based quantification that supports this approach. Document the inflection point value and the number of cells retained.

Mitochondrial Read Fraction

Mitochondrial reads measure cell health, but they also reflect cell type. Cells that rely on oxidative phosphorylation have higher baseline mito fractions. For example, a study on keloid infiltrating zones (source) reported stromal cells with higher mitochondrial activity than immune cells. Your cutoff should be a percentile based decision, not a fixed number. Examine the bimodal distribution of mito fraction. Set the threshold at the upper tail of the main population, not at an arbitrary constant.

Practical Workflow

Step 1: Generate Raw Counts

Begin with filtered feature barcode matrices from your alignment software. Ensure that you have removed ambient RNA using tools such as EmptyDrops or SoupX. The EMBL-EBI Training resource includes practical modules on preprocessing scRNA seq data that explain ambient RNA correction.

Step 2: Calculate QC Metrics

Compute per cell metrics: number of UMIs, number of genes, and fraction of reads mapped to mitochondrial genes. Also compute per gene metrics: total counts across cells and number of cells expressing the gene. Use libraries from Bioconductor such as scater or scran.

Step 3: Plot and Inspect

Generate violin plots for each metric. Look for outlier populations. Create scatter plots of UMI count vs gene count and color by mito fraction. Identify clusters of cells that deviate from the main cloud. Do not rely on automated outliers alone. Manual inspection is essential. Refer to the workflow described in source for a protocol that includes quality inspection steps.

Step 4: Apply Contextual Thresholds

Set lower bounds for genes and UMIs based on the knee point. Set upper bounds based on the 95th or 99th percentile, not on a universal value. For mito fraction, set the threshold at the inflection point where the density curve drops. Alternatively, use a MAD (median absolute deviation) approach. For example, keep cells within 3 MADs of the median for each metric. Document every cutoff value.

Step 5: Filter and Validate

Remove cells that fail any filter. Then check whether expected cell types remain. For example, if you lose all T cells after filtering mito reads, your threshold may be too strict. Reexamine the biological basis. Source discusses cell type specific stress responses that can mimic poor quality. Iterate until you retain a diverse population consistent with your tissue.

Step 6: Document Your Choices

Write a QC report that includes the following:

Number of cells before and after filtering.
Threshold values and the justification for each.
Number of genes retained (typically genes expressed in at least 3 to 10 cells).
Software versions and parameters.
Any doublet removal steps.

This documentation is critical for reproducibility. The Galaxy Training Network offers templates for tracking analysis steps.

Common Mistakes

Overfiltering Healthy Cells

The most frequent error is applying a strict mitochondrial filter to tissues with high oxidative metabolism. Cells from sources and show that mito reads can reach 30% in viable cells. Remove these cells and you remove the biological signal you want to study.

Ignoring Doublet Detection

Gene and UMI filters alone do not remove doublets. Doublets look like large healthy cells with high gene counts. You must use computational doublet detection. Source describes cell interaction studies that require clean single cell data. Doublets would confound those analyses.

Using the Same Thresholds for Different Libraries

Batch effects are real. Each sequencing run may have different quality distributions. Apply QC per sample or per batch, not globally. Source on the Galaxy Training Network emphasizes per sample quality checks.

Not Checking Ambient RNA

Empty droplets contain ambient RNA that looks like real cells with low UMIs. Prefiltering by UMI count removes many empties, but some can persist. Use dedicated tools before QC filtering.

Limits and Uncertainty

Technical Noise from Droplet Based Methods

Droplet protocols produce variable capture efficiency. Cells with low UMI counts may be real cells with low expression or technical failures. You cannot distinguish them with certainty. Use a conservative lower bound to retain small cells such as lymphocytes.

Biological Heterogeneity Within the Same Sample

Tissues are heterogeneous. A single threshold for gene count may discard a rare cell type with low transcriptional activity. For example, quiescent stem cells express fewer genes than proliferating ones. Evaluate whether your filters disproportionately remove one cluster.

Uncertainty in Mitochondrial Read Origin

High mitochondrial reads can indicate cell stress, but they can also indicate contamination from platelet aggregates or cell fragments. Source in plant single cell studies shows that mitochondrial reads vary by cell type even in healthy tissue. Your filter may inadvertently select for one cell type over another.

Iterative Filtering Is Required

Do not expect to set filters once. After initial filtering and clustering, reexamine the QC metrics per cluster. If one cluster has unusually high mito reads, decide whether it is a dying cell population or a genuine metabolic state. This iterative approach is described in the context of liver and adipose tissue profiling in source.

Frequently Asked Questions

Why do some tutorials recommend a 5% mitochondrial cutoff while others say 20%?

The recommended cutoff depends on the tissue and protocol. Cultured cell lines and immune cells from blood often have low mitochondrial reads. Solid tissues with high energy demand contain cells with naturally higher mito fractions. Always inspect your own data distribution before setting a cutoff.

Should I remove cells with very high gene counts?

Not automatically. High gene counts can indicate doublets, but they can also indicate large cells such as macrophages or megakaryocytes. Use doublet detection software to distinguish. If doublet scores are low, keep those cells.

What do I do if my data shows a bimodal distribution in mitochondrial reads?

Examine which cell types fall into each mode. The high mode may represent a genuine cell population such as stressed or metabolically active cells. Consult tissue specific references like source or to see if that pattern is expected.

Is it better to filter cells by gene count or by UMI count?

Both metrics are needed. Gene count reflects complexity. UMI count reflects capture efficiency. A cell with many UMIs but few genes may be a technical artifact. A cell with many genes but few UMIs may be a damaged cell with high ambient contamination. Use both in combination.