The Delta Delta Ct Method Demystified: Calculating Relative Gene Expression Step by Step

The delta delta Ct method turns two Cq values per sample into a fold change: normalize the target to a reference gene (ΔCt), compare treated with control (ΔΔCt), and report 2^−ΔΔCt. The arithmetic takes a minute. The two mistakes it hides take longer to spot: amplification efficiency that is not quite 100%, and statistics run on fold changes instead of Cq values.

This SOP runs from raw Cq values to a fold change with a correct confidence interval, starting with the efficiency check that decides whether ΔΔCt is the right formula at all.

Ask Sophie about your qPCR analysis

Why the delta delta Ct method works, and what it assumes

At 100% efficiency, every cycle doubles the product. A sample that crosses threshold one cycle earlier started with twice as much template, so Cq differences convert to fold differences through powers of 2:

ΔCt = Cq(target) − Cq(reference) ΔΔCt = ΔCt(treated) − ΔCt(control) Fold change = 2^−ΔΔCt

Subtracting the reference gene removes differences in how much RNA went in. Subtracting the control removes everything the two groups share. What remains is the treatment effect.

The "2" is the assumption. It holds only when both assays double their product every cycle. Equal efficiencies are not enough: if both assays run at 90%, the product grows 1.9-fold per cycle, and a true 6.9-fold change is reported as 8-fold. At a ΔΔCt of −5 the gap widens to 24.8-fold reported as 32. The bias grows with the effect size, so it lands hardest on your biggest result.

Before you start

• Cq values for your target and at least one validated, stable reference gene

• Technical triplicates for every sample and assay

• At least three biological replicates per group (independent cultures or animals)

• A pooled cDNA sample for a dilution series

• A spreadsheet, or statistics software that runs a t-test and reports confidence intervals

Checking amplification efficiency

Step 1. Prepare a five-point dilution series of pooled cDNA spanning at least a 100-fold range. Run it with both the target and reference assays.

Step 2. For each assay, plot Cq against log₁₀ of input and fit a straight line. Calculate efficiency as E = 10^(−1/slope). A slope of −3.32 gives E = 2.0, which is 100%.

Step 3. Confirm the two efficiencies are equal. Plot ΔCt (target minus reference) against log₁₀ input. The absolute value of the slope should be below 0.1; a flat line means the two assays amplify at the same rate.

Step 4. Choose the formula:

• Both assays close to 100% and the Step 3 slope below 0.1: use ΔΔCt (Steps 5–11).

• Either assay clearly below 100%, or the Step 3 slope above 0.1: use the Pfaffl ratio (Step 12), or redesign the assay.

Caution: A passing Step 3 slope proves the efficiencies match. It does not prove they are 100%, and ΔΔCt needs both.

Averaging technical replicates

Step 5. For each sample, average the three technical Cq values for the target, and separately for the reference.

Note: If one well sits clearly apart from the other two, check the plate record for a pipetting error or bubble before averaging. Record any exclusion.

Step 6. Carry one averaged Cq per assay per biological sample into the calculation.

Caution: Technical replicates measure pipetting precision and are never your n. Treating 3 wells from one culture as n = 3 is pseudoreplication.

Running the calculation

Step 7. Calculate ΔCt for every sample: Cq(target) − Cq(reference). With several reference genes, use the mean of their Cq values, which is the geometric mean of their quantities.

Step 8. Calculate the calibrator as the mean ΔCt of the control group, then ΔΔCt for every sample as its ΔCt minus that mean.

Note: Using one control sample as the calibrator builds that sample's noise into every fold change. The group mean does not, and it lets the control samples scatter around 1 as they should.

Step 9. Calculate each sample's fold change as 2^−ΔΔCt. For the group, calculate 2^−(mean ΔΔCt).

Caution: Never average fold changes arithmetically. A 2-fold increase and a 2-fold decrease average to 1.25, though the net change is zero. Averaging on the ΔΔCt scale gives the correct 1.00.

Running the statistics on ΔCt

Step 10. Test the ΔCt values of the treated group against the control group with a t-test (or ANOVA for three or more groups). Cq is the scale on which the data behave; fold changes are skewed.

Step 11. Take the 95% confidence interval of the ΔΔCt and convert both ends: the fold-change interval runs from 2^−(upper limit) to 2^−(lower limit). It will be asymmetric around the fold change, and that asymmetry is correct.

Tip: Paste your Cq table into Soφ with a one-line description ("target and GAPDH, three controls, three treated, triplicate wells") and have it return the fold change with its interval.

When efficiency is not 100%: the Pfaffl ratio

Step 12. Use each assay's measured efficiency from Step 2:

Ratio = E_target^ΔCq_target / E_reference^ΔCq_reference

where each ΔCq is the control-group mean Cq minus the treated-group mean Cq for that assay. With both efficiencies at exactly 2.0, this reduces to 2^−ΔΔCt.

Worked example

Illustrative data. Each Cq is the mean of a technical triplicate.

Sample	Target Cq	Reference Cq	ΔCt	ΔΔCt	Fold change
Control 1	26.10	18.00	8.10	0.00	1.00
Control 2	26.60	18.30	8.30	+0.20	0.87
Control 3	25.80	17.90	7.90	−0.20	1.15
Treated 1	24.40	18.20	6.20	−1.90	3.73
Treated 2	24.00	17.80	6.20	−1.90	3.73
Treated 3	24.90	18.10	6.80	−1.30	2.46

The control mean ΔCt (the calibrator) is 8.10 and the treated mean is 6.40, so ΔΔCt = −1.70 and the fold change is 3.25.

A t-test on the ΔCt values gives t(4) = 7.36, p = 0.0018. The 95% CI of the ΔΔCt is −2.34 to −1.06, so the fold-change CI is 2.08 to 5.07, wider above 3.25 than below it. Running the same t-test on the fold-change values instead gives p = 0.0059: a different answer from the same data, on the wrong scale.

Had the target assay run at 90% efficiency (E = 1.90) with the reference at 100%, the Pfaffl ratio gives 2.97-fold. The 2^−ΔΔCt figure of 3.25 would overstate it by 9%.

Reporting the result

Report the fold change with its asymmetric 95% CI, the test on ΔCt, each assay's efficiency, the reference genes, and the formula used. MIQE 2.0 (2025) asks for efficiency-corrected quantities, so state the efficiencies even when they justify ΔΔCt.

Troubleshooting

Symptom	Cause	Fix
Fold changes in control samples far from 1	Single-sample calibrator, or unstable reference gene	Use the control-group mean ΔCt; revalidate the reference gene
Step 3 slope above 0.1	Target and reference efficiencies differ	Use the Pfaffl ratio or redesign the weaker assay
Error bars symmetric around the fold change	Error calculated on the fold-change scale	Recalculate the CI on ΔΔCt and back-transform
Technical triplicate spread widely	Pipetting error, bubble, low template	Check the plate record; rerun if the cause is unclear
Implausibly large fold change	Efficiency below 100% inflating 2^−ΔΔCt	Measure efficiency; apply Pfaffl

The delta delta Ct method is a minute of arithmetic and two decisions

The subtraction is the easy part. The number becomes trustworthy when you have checked that both assays really double every cycle, and run the statistics on ΔCt, where the data behave well. Get those two decisions right and 2^−ΔΔCt is a sound, fast way to report relative expression.

Setting up the qPCR run itself?

Find the bench protocol in the CLYTE SOP library.

Frequently asked questions

What does a negative ΔΔCt mean? Higher expression in the treated group. The target crossed threshold earlier relative to control, so 2^−ΔΔCt is greater than 1.

Can I use ΔΔCt if my efficiency is 95%? You can, if you accept the error: at 95% efficiency on both assays, a ΔΔCt of −3 is overstated by about 8%. For small fold changes that may be tolerable; for publication, the Pfaffl ratio removes it.

Should I run statistics on ΔCt or on fold change? On ΔCt. Fold changes are skewed, so tests and error bars calculated on them are distorted. Convert only the final estimate and its confidence limits.

Which sample should be the calibrator? The mean ΔCt of the whole control group, not a single control sample.

References

1. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2^−ΔΔCT method. Methods, 2001;25(4):402–408. PubMed

2. Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Research, 2001;29(9):e45. DOI

3. Yuan JS, Reed A, Chen F, Stewart CN. Statistical analysis of real-time PCR data. BMC Bioinformatics, 2006;7:85. DOI

4. Applied Biosystems / Thermo Fisher Scientific. Guide to Performing Relative Quantitation of Gene Expression Using Real-Time Quantitative PCR.

5. Bustin SA, Ruijter JM, van den Hoff MJB, et al. MIQE 2.0: revision of the minimum information for publication of quantitative real-time PCR experiments guidelines. Clinical Chemistry, 2025;71(6):634–651. PubMed