RNA-seq Differential Expression — Transcriptomic DE Analysis
RNA Sequencing Differential Expression Analysis · Also known as: RNA-seq DE analysis, transcriptomic differential expression, bulk RNA-seq DE, DEA
RNA-seq differential expression (DE) analysis identifies genes whose transcript abundance differs significantly between two or more biological conditions — for example, treated versus control, or diseased versus healthy tissue. Starting from raw sequencing reads, the pipeline moves through alignment, count-based normalization, statistical modeling of count dispersion, hypothesis testing, and multiple-testing correction to produce a ranked list of differentially expressed genes accompanied by fold-change estimates and adjusted p-values.
Read the full method
Sign in with a free account to read this section.
Method map
The neighbourhood of related methods — select a node to explore.
+50 more
When to use it
Use RNA-seq DE analysis when you have bulk RNA sequencing data from two or more biological conditions and want to identify genes whose expression levels change in a statistically robust way. The method requires a minimum of 3 biological replicates per condition (6 total) to estimate dispersion reliably; 3–5 replicates is the practical minimum, with 5–8 providing substantially better power. Do not use this pipeline for single-cell RNA-seq data, which has different sparsity characteristics and requires scRNA-seq-specific tools (Seurat, Scanpy). Also avoid applying it when replicates are actually technical rather than biological, when total RNA input is very low (high ERCC spike-in proportion), or when the experimental design confounds batch with treatment — correct for batch effects explicitly before testing.
Strengths & limitations
- Genome-wide, unbiased survey of transcript abundance without prior hypothesis about which genes are involved.
- Well-validated statistical frameworks (DESeq2, edgeR, limma-voom) with extensive benchmarking literature.
- Handles complex multi-factor designs including batch correction, paired samples, and interaction terms within a GLM framework.
- Produces quantitative fold-change estimates alongside significance calls, enabling biological prioritization.
- Seamlessly integrates with downstream enrichment and pathway analysis tools via standardized output formats.
- Supported by mature, actively maintained Bioconductor packages with thorough documentation.
- Requires biological replicates — pseudo-replication (e.g., technical replicates only) inflates false positives drastically.
- Bulk RNA-seq averages expression across all cells in a sample; cell-type-specific effects are obscured without deconvolution.
- Low-count genes have high variability and low statistical power; very lowly expressed transcripts are often filtered before testing.
- Results are sensitive to alignment choice, quantification strategy, and normalization method — pipelines must be documented and justified.
- Multiple-testing correction at FDR 0.05 still permits a fraction of false discoveries in large gene lists.
Frequently asked
DESeq2 or edgeR — which should I use?
Both are valid and well-benchmarked. DESeq2 is often preferred for smaller sample sizes due to its shrinkage estimators, and its default output includes shrunken log2 fold changes suitable for ranking. edgeR's quasi-likelihood pipeline performs comparably and is flexible for complex designs. Many analyses run both and report consistent findings. The choice matters less than having adequate replication and a well-documented pipeline.
How many replicates do I need?
A minimum of 3 biological replicates per condition is necessary to estimate dispersion; 3 is workable but underpowered for lowly expressed genes. Power analyses (e.g., using the RnaSeqSampleSize or RNASeqPower R packages) should guide design. For rare-event detection or subtle fold changes, 5–8 replicates per condition provide substantially better sensitivity.
What is the difference between biological and technical replicates?
Biological replicates are independent samples from different individuals or culture wells — they capture true biological variability. Technical replicates are repeated measurements of the same biological material. DE analysis requires biological replicates because it aims to generalize beyond a single sample. Using technical replicates as if they were biological replicates severely underestimates variance and inflates false positives.
My gene list is huge (thousands of DE genes) — is something wrong?
Not necessarily, but it warrants scrutiny. A very large DE gene list can result from a strong biological effect (e.g., a drastic treatment), but also from batch confounding, failure to apply fold-change filtering, or overly liberal FDR thresholds. Inspect your PCA plot: if PC1 separates conditions cleanly and accounts for a large fraction of variance, a large DE list may be genuine. Also apply a minimum absolute log2 fold change threshold (e.g., |log2FC| > 1) alongside the FDR cutoff.
Can I use RNA-seq DE analysis for single-cell data?
Standard bulk tools (DESeq2, edgeR) should not be applied directly to raw single-cell count matrices because single-cell data are far sparser and zero-inflated. Pseudo-bulk approaches — aggregating single-cell counts to the sample level and then running bulk DE tools — are currently recommended for differential expression between conditions in scRNA-seq experiments. Within-dataset cluster marker detection uses different methods (Wilcoxon rank-sum, MAST, etc.) that are designed for single-cell sparsity.
Sources
- Love, M. I., Huber, W., & Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biology, 15(12), 550. DOI: 10.1186/s13059-014-0550-8 ↗
- Robinson, M. D., McCarthy, D. J., & Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics, 26(1), 139–140. DOI: 10.1093/bioinformatics/btp616 ↗
How to cite this page
ScholarGate. (2026, June 3). RNA Sequencing Differential Expression Analysis. ScholarGate. https://scholargate.app/en/bioinformatics/rna-seq-differential-expression
Which method?
Set this method beside its closest kin and read them side by side — the library lays the books on the table; the choice is yours.
- ChIP-seq Peak CallingBioinformatics↔ compare
- Gene Set Enrichment AnalysisBioinformatics↔ compare
- Pathway Enrichment AnalysisBioinformatics↔ compare
- Sequence AlignmentBioinformatics↔ compare
- Single-cell RNA-seq analysisBioinformatics↔ compare
- Variant CallingBioinformatics↔ compare