differential_expression
Read-only
Differential expression from raw RNA-seq read counts, by one of the three standard count models. DESeq2 (Love, Huber & Anders 2014): median-of-ratios size factors, negative-binomial GLM with empirical-Bayes shrunken dispersions, Wald or likelihood-ratio test, Cook's-distance outlier handling, independent filtering, optional LFC shrinkage. edgeR quasi-likelihood (Chen et al. 2025): filterByExpr, TMM normalization, a negative-binomial GLM with bias-adjusted QL dispersions squeezed toward an abundance trend, QL F-test. limma-voom (Law et al. 2014): filterByExpr, TMM, log-CPM with mean-variance precision weights, weighted linear models, moderated t. Takes a gene-by-sample count table and a sample table with the design's variables. Every gene of every validated fit is bit-identical to DESeq2 1.52.0, edgeR 4.10.6 and limma 3.68.5 running in R. Hands its significant genes, tested genes and ranked list on for gene-set enrichment.
Input Schema
TableJSON Schema
| Name | Required | Description | Default |
|---|---|---|---|
| name | No | one fitted coefficient to report, by its DESeq2 name (e.g. 'condition_treated_vs_untreated') or its model-matrix column ('conditiontreated'); empty is the last. | |
| test | No | Wald: one coefficient or contrast. LRT: the full design against `reduced`, testing every removed term at once. Only applies when method is DESeq2. | Wald |
| alpha | No | the target FDR: the summary counts genes below this adjusted p, and DESeq2's independent filtering maximizes rejections at it. | |
| _files | No | The project: record references shaped like _subject. Their records become library documents a tool names in sequences/parts/document; a trace or alignment is found by file name. | |
| counts | Yes | raw read counts (never normalized values or TPM): genes in rows with the gene ID first, samples in columns, tab- or comma-separated, a header row naming the samples. | |
| design | No | the design formula over the sample table's columns, the variable of interest last (e.g. '~ batch + condition'). | ~ condition |
| method | No | the count model: DESeq2 (negative binomial, Wald or LRT), edgeR-QL (negative binomial with quasi-likelihood F-test) or limma-voom (precision-weighted linear model on log-CPM, moderated t). | DESeq2 |
| sfType | No | size factors: median of ratios over genes with no zero (ratio), or with zeros allowed (poscounts, for sparse data). Only applies when method is DESeq2. | ratio |
| shrink | No | shrink the reported coefficient's log2 fold changes toward zero where the counts cannot support them (DESeq2 lfcShrink), for ranking and plotting by effect size; p-values stay the unshrunken test's. apeglm: adaptive Cauchy prior, posterior mode. ashr: mixture-of-normals prior, posterior mean. normal: DESeq2's original zero-centred Normal prior (Wald only). Only applies when method is DESeq2. | none |
| fitType | No | the dispersion trend: parametric (asymptDisp + extraPois / mean, DESeq2's default) or the mean of gene-wise estimates. Only applies when method is DESeq2. | parametric |
| reduced | No | the reduced design for the likelihood-ratio test. Only applies when method is DESeq2. | |
| samples | Yes | the sample table: a header row, then one row per sample with its name (as in the count table's header) first and the design's variables after; the run uses this table's sample order. | |
| _subject | No | The record to act on: exactly one of sequence, accession (+source), sample (an id from list_samples) or file ({name, text} or {name, base64}); record_index picks one record of several. | |
| contrast | No | 'factor, numerator, denominator': the log2 fold change of numerator over denominator. Empty is the design's last coefficient. Under edgeR and voom the denominator must be the factor's reference level. | |
| minCount | No | filterByExpr: the count, at the median library size, a gene needs in the smallest group's worth of samples (taken as CPM). Only applies when method is edgeR-QL or limma-voom. | |
| reference | No | reference level per factor, 'variable=level' separated by commas; a factor not named takes its alphabetically first level, as R does. | |
| cooksCutoff | No | the Cook's distance cutoff; empty is DESeq2's F(0.99; p, m - p). Only applies when method is DESeq2. | |
| cooksFilter | No | set the p-value of genes with a Cook's-distance outlier to NA. Only applies when method is DESeq2. | |
| filterByExpr | No | keep only genes edgeR's filterByExpr finds expressed (the edgeR and limma workflows' first step); off fits every gene. Only applies when method is edgeR-QL or limma-voom. | |
| lfcThreshold | No | test |log2 fold change| against this threshold instead of 0. Only applies when method is DESeq2. | |
| altHypothesis | No | which side of the threshold is the alternative. Only applies when method is DESeq2. | greaterAbs |
| minTotalCount | No | filterByExpr: the total count a gene needs across all samples. Only applies when method is edgeR-QL or limma-voom. | |
| pAdjustMethod | No | multiple-testing adjustment across genes (R `p.adjust`). | BH |
| independentFiltering | No | drop low-mean genes from the adjustment when that gains rejections at `alpha`. Only applies when method is DESeq2. | |
| minReplicatesForReplace | No | samples per group from which outlier counts are replaced and the gene refitted, instead of flagged. Only applies when method is DESeq2. |
Output Schema
TableJSON Schema
| Name | Required | Description | Default |
|---|---|---|---|
| genes | Yes | one row per gene: log2FoldChange, stat (DESeq2 Wald or LRT, edgeR F, voom moderated t), pvalue, padj, why padj is NA, the normalized counts per sample; DESeq2 adds baseMean, lfcSE, the gene-wise, fitted and final dispersions, whether it was refitted and with `shrink` the shrunken fold change, its SE and s-value; edgeR and voom add logCPM, the raw, prior and posterior variance (edgeR's QL dispersion) and the total df, and voom the B-statistic | |
| maPlot | Yes | log2 fold change (shrunken, with `shrink`) against mean normalized count (DESeq2) or average log2 CPM (edgeR, voom), significant genes apart | |
| summary | Yes | genes up and down at alpha and tested, with the method's own fit: DESeq2's outliers, low counts, filter threshold, size factors, dispersion trend and shrinkage prior; edgeR's common and average QL dispersion and prior df; voom's span and prior variance; the filterByExpr cutoffs and library sizes | |
| upGenes | Yes | significant genes with positive log2 fold change, one per line | |
| volcano | Yes | -log10 p against log2 fold change, significant genes apart | |
| downGenes | Yes | significant genes with negative log2 fold change, one per line | |
| rankedList | Yes | tested genes by their test statistic, decreasing: DESeq2's Wald statistic, voom's moderated t, or edgeR's signed square root of F (a t with the same df): the ranked list for GSEA | |
| testedGenes | Yes | every gene with an adjusted p-value, one per line: the background (universe) for over-representation analysis | |
| dispersionPlot | Yes | DESeq2: gene-wise, fitted and final dispersions against mean count. edgeR: quarter-root raw, squeezed and trend QL dispersions against average log2 CPM (plotQLDisp). voom: sqrt standard deviation against log2 count size, with the lowess trend that sets the weights | |
| significantGenes | Yes | genes with padj below alpha, one per line: the gene list for over-representation analysis |