论文里的 qPCR 怎么复现:方法部分常漏写的参数、Ct 值和熔解曲线排查
Vara 团队 · 更新于 2026-10-21
这页讲怎么复现论文里的 qPCR 结果:按作者的方法重建反应体系、扩增程序和数据分析。多数论文的方法部分只给了引物名称和循环数,扩增效率、内参验证和对照设置常常不写。下面列出要找的信息、关键参数之间怎么联动,以及常见问题怎么排查。文中的 Ct 值在 MIQE 指南里叫 Cq,是同一个东西。
方法部分常漏写什么
- 每对引物的扩增效率
2^-ΔΔCt 法假设目的基因和内参的扩增效率都接近 100% 且彼此相近。没有作者的标准曲线,就无法判断他们的倍数变化换成你的引物后是否还成立。
- 引物序列和引物终浓度
只写基因名不能确定用的是哪个检测。同一个基因的不同引物对,扩增效率、特异性和覆盖的转录本都可能不一样。
- 内参基因怎么选的,有没有验证稳定性
如果内参随处理变化,所有归一化后的数值都会跟着变。论文经常直接写 GAPDH 或 ACTB,却没有证明它在实验条件下是稳定的。
- RNA 质量和反转录体系
RNA 投入量、RNA 完整性、反转录引物(oligo-dT、随机六聚体或两者混用)和反转录酶,在 PCR 开始之前就已经影响了 Ct 值。
- 对照:无模板对照(NTC)和无反转录对照(no-RT)
没有这两个对照,就无法排除污染或基因组 DNA 扩增,而这两者正是 Ct 值偏晚、复孔间波动大的常见来源。
- n 指的是什么
同一份 RNA 做 3 个复孔是技术重复,不是 3 个生物学重复。很多论文写 n=3,却不说明是哪一种。
关键参数和联动关系
| 参数 | 常见取值 | 会跟着变的 |
|---|---|---|
| 反应体系 | 每孔 10 到 20 µL | 每孔 Master Mix、引物和模板的加样体积,以及整板需要的 Master Mix 总量。 |
| 引物终浓度 | 常见每条 200 到 500 nM,每个检测单独优化 | 引物储存液的稀释和每孔加样体积;扩增效率和引物二聚体的形成。 |
| 模板投入量 | 由固定量 RNA 反转录得到的 cDNA,使用前常先稀释 | 预期的 Ct 值范围,必须落在标准曲线的线性范围内。 |
| 退火和延伸温度 | 两步法常用 60 °C | 引物 Tm 值和特异性;熔解曲线的形状。 |
| 反应数 | 样品数 × 复孔数 × 基因数,再加上对照和移液损耗的富余量 | 每种组分的总用量。样品数一变,所有总量都要重算。 |
排查表
| 现象 | 可能原因 | 先查什么 |
|---|---|---|
| 扩增效率不在 90% 到 110% 之间(标准曲线斜率不在约 -3.6 到 -3.1 之间) | 引物设计不好、RNA 提取时残留了抑制物,或稀释梯度超出了线性范围。 | 重做 5 个点的梯度稀释,去掉两端的点,看 A260/230 判断有没有残留,斜率还是不对就重新设计引物。 |
| 同一个样品不同批次跑出来 Ct 值漂移 | 模板反复冻融后降解、小体积加样误差,或换了 Master Mix 批号。 | cDNA 分装保存,每块板都放一个校准样品,把模板稀释后用更大的体积加样。 |
| 熔解曲线出现双峰 | 引物二聚体、另一个非特异产物,偶尔也可能是单一扩增子有不止一个解链区;跑胶可以区分。 | 把产物跑胶看,降低引物浓度,或提高退火温度。 |
| 无模板对照(NTC)有信号 | 试剂污染或引物二聚体。 | NTC 的熔解峰和产物一致,就换试剂;如果是 Tm 偏低的峰,多半是二聚体。 |
| 无反转录对照(no-RT)有信号 | 基因组 DNA 污染。 | 加一步 DNase 处理,并使用跨外显子连接处的引物。 |
拆解示例:从一句方法原文到能照着做的步骤
Total RNA was reverse transcribed and target genes were quantified by SYBR Green qPCR using GAPDH as the internal control. Relative expression was calculated by the 2^-ΔΔCt method.
每一步按原文写得清不清楚打标:原文写明、部分写明、推断、原文缺失。
- 1.总 RNA 反转录成 cDNA原文写明
- 2.RNA 投入量和反转录试剂盒原文缺失
没写;要对上 Ct 值范围需要这个信息。
- 3.对每个目的基因和 GAPDH 做 SYBR Green qPCR原文写明
- 4.引物序列和终浓度原文缺失
先查补充材料,别急着自己猜。
- 5.扩增程序推断
用标准两步法起步是合理的,但这是一个假设。
- 6.确认每对引物的扩增效率和 GAPDH 接近原文缺失
这是 2^-ΔΔCt 法成立的前提。
- 7.用 2^-ΔΔCt 法计算相对表达量原文写明
- 8.生物学重复数原文缺失
方法部分没写;去图注里找,并确认 n 是生物学重复还是技术重复。
Vara 能做什么
- 粘贴方法段落或上传 PDF,Vara 把它拆成带参数的步骤,每一步标出原文写明、部分写明、推断或原文缺失,论文没写的地方一眼就能看到。
- 扩增效率、内参验证、n 指什么这类缺失项会列成复现风险,不会被悄悄补上。
- 反应体系、反应数和各组分用量用公式关联。改了样品数,Master Mix 总量会重新算。
- 每次实验都记在对应的方案版本下,Ct 值漂移时可以对比前后改了什么。
这类实验用得上的免费速算器
常见问题
用 2^-ΔΔCt 法,扩增效率要达到多少?+
目的基因和内参的扩增效率都要接近 100%,而且彼此相差在几个百分点以内。常用的接受范围是 90% 到 110%。如果两者效率差别较大,改用校正效率的模型,例如 Pfaffl 法。
我的 Ct 值为什么和论文对不上?+
Ct 值的绝对数值受 RNA 投入量、反转录、仪器和 Master Mix 影响,不同实验室之间很少能一致。应该比较相对表达量,并确认扩增效率和内参稳定性都合格。
GAPDH 做内参一定可靠吗?+
不一定。GAPDH 的表达会随处理、代谢状态和细胞类型变化。在你的实验条件下测几个候选内参,用 geNorm 这类稳定性指标排序,再用最稳定的两到三个基因的几何平均值做归一化。
论文没给引物序列怎么办?+
先查补充材料,以及同一课题组之前被引用的论文。还是找不到,就联系通讯作者。换成别的引物也能做,但那样就不再是在复现同一个检测了。
熔解曲线出现两个峰是什么原因?+
常见的是引物二聚体或另一个非特异产物,偶尔也可能是单一扩增子有不止一个解链区。把产物跑胶就能区分。确认是二聚体或非特异产物的话,可以降低引物浓度或提高退火温度。
参考资料
- Bustin 等 2009,MIQE 指南,Clinical Chemistry
- Livak 与 Schmittgen 2001,2^-ΔΔCt 法,Methods
- Pfaffl 2001,校正扩增效率的相对定量,Nucleic Acids Research
- Vandesompele 等 2002,内参基因归一化(geNorm),Genome Biology