2^-ΔΔCt 的公式只有三行,但两个细节决定了结果对不对:校准用什么、检验做在哪个量上。
一、三步算法(务必按顺序)
- ΔCt = Ct(目的基因)− Ct(内参基因),逐样本计算。
- ΔΔCt = ΔCt(处理组)− mean ΔCt(对照组)。注意是对照组 ΔCt 的**均值**。
- 相对表达量 = 2^(−ΔΔCt)。对照组均值因此被归一为 1,处理组是相对倍数。
二、错误一:用单个对照样本做校准
- 很多教程写 "ΔΔCt = ΔCt(处理) − ΔCt(对照样本1)"。这会让整批数据随这一个样本的随机误差整体偏移。
- 正确做法是用对照组所有重复的 ΔCt 均值;对照组的离散度同时反映技术重复的好坏。
- 如果对照组 ΔCt 的 SD > 0.5,先查加样与内参稳定性,不要急着解释生物学差异。
三、错误二:对 2^-ΔΔCt 做 t 检验
- 倍数经过指数变换后分布明显右偏,直接做 t 检验会低估 P 值(更容易"显著")。
- 规范做法:检验做在 ΔCt(近似正态)上,再用 ΔCt 的均值差换算倍数。
- 组内方差为 0 时(重复完全一致),Welch t 检验会给出虚高的 P 值,此时应改用非参数检验或直接说明无法检验。
四、报告怎么写
- 写清内参基因、重复数、校准组,以及"以对照组 ΔCt 均值校准"。
- 图:柱状图 + 单点散点(不要只画均值柱),误差棒注明是 SD 还是 SEM。
- 给出统计方法与检验变量(ΔCt 还是倍数),这一句能挡掉大部分审稿质疑。
五、用接口直接算(可复制)
四组 Ct 值(对照组目的/内参、处理组目的/内参)用逗号分隔,接口直接返回倍数、P 值、所用的检验方法与异常值提示。
REST:2^-ΔΔCt
curl -X POST https://sanyu.wentangkeji.com/api/open/v1/run/act_qpcr_ddct \
-H "Authorization: Bearer sk-sanyu-你的密钥" \
-H "Content-Type: application/json" \
-d '{"gene":"NLRP3","ref_gene":"GAPDH",\
"ctrl_ct":"22.5,22.8,22.6","ctrl_ref_ct":"18.3,18.1,18.4",\
"exp_ct":"20.1,20.4,20.2","exp_ref_ct":"18.2,18.5,18.3"}'
常见问题
- 2^-ΔΔCt 和 2^-ΔCt 有什么区别?
- 2^-ΔCt 只是相对内参归一后的表达量,没有跨组比较;2^-ΔΔCt 是相对"对照组均值"的倍数,对照组被归一到 1。写论文报告组间差异时用的是后者,且必须说明校准组。
- 内参基因怎么选?
- 选在该组织/处理下表达稳定的管家基因(GAPDH、ACTB、18S、HPRT1 等),并最好用两个内参做几何平均。注意:缺氧、增殖、药物处理常会改变 GAPDH 或 ACTB 的表达,这时用它们做内参本身就会引入偏差。
- 重复之间 Ct 差多少算正常?
- 技术重复的 SD 通常应在 0.3 个 Ct 以内(约 20% 定量误差);超过 0.5 提示加样或扩增效率问题,应先排查再解释生物学差异。