利用实时定量PCR和2－△△CT法分析基因相对表达量
METHODS 25, 402–408 (2001)
An经营理念
alysis of Relative Gene Expression Data Using Real-Time Quantitati ve PCR and the 2－△宠物品种
△CT Method
Kenneth J. Livak* and Thomas D. Schmittgen,1
*Applied Biosystems, Foster City, California 94404; and Department of Pharmaceutical Sciences, College of Pharmacy,
Washington State University, Pullman, Washington 99金山旅游
164-6534
摘要：
现在最常用的两种分析实时定量PCR 实验数据的方法是绝对定量和相对定量。绝对定量通过标准曲线计算起始模板的拷贝数；相对定量方法则是比较经过处理的样品和未经处理的样品目标转录本之间的表达差异。2－△△CT方法是实时定量P CR 实验中分析基因表达相对变化的一种简便方法，即相对定量的一种简便方法。本文介绍了该方法的推导，假设及其应用。另外，在本文中我们还介绍了两种2－△△CT衍生方法的推导和应用，它们在实时定量 PCR 数据分析中可能会被用到。
关键词：反转录PCR 定量PCR 相对定量实时PCR Taqman
反转录 PCR （RT-PCR ）是基因表达定量非常有用的一种方法（1 - 3象征美好的词语
）。实时PCR 技术和RT-PCR 的结合产生了反转录定量 PCR 技术（4 ,5 ）。实时定量PCR 的数据分析方法有两种：绝对定量和相对定量。绝对定量一般通过定量标准曲线来确定我们所感兴趣的转录本的拷贝数；相对定量方法则是用来确定经过不同处理的样品目标转录本之间的表达差异或是目标转录本在不同时相的表达差异。
绝对定量通常在需要确定转录本绝对拷贝数的条件下使用。通过实时 PCR 进行绝对定量已有多篇报道（6 - 9 ），包括已发表的两篇研究论文（10,11 ）。在有些情况下，并不需要对转录本进行绝对定量，只需要给出相对基因表达差异即可。显然，我们说 X 基因在经过某种处理後表达量增加 2.5 倍比说该基因的表达从1000 拷贝/ 细胞增加到2500 拷贝/ 细胞更加直观。
用实时PCR 对基因表达进行相对定量分析需要特殊的公式、假设以及对这些假设的验证。2－△△CT方法可用于定量PCR 实验来计算基因表达的相对变化：2－△△CT 公式的推导，以及实验设计，有效性评估在Applied Biosystems Ur Bulleti n No.2（P/N4303859）中有介绍。用2－△△CT方法分析基因表达数据在文献中也有
报道（5,6）。本文介绍了该方法的推导、假设以及应用。另外，本文还介绍了2－△△CT两种衍生方法的推导和应用，它们在实时定量PCR 数据分析中都可能被用到。
1. 2－△△CT方法
1.1. 2－△△CT方法的推导
PCR 指数扩增的公式是：
这里，X
n 是第 n 个循环後目标分子数，X
是初始目标分子数，E
x
是目标分子扩
增效率，n 是循环数，C
T
代表目标扩增产物达到设定阈值所经历的循环数。
因此：
X
T
是目标分子达到设定的阈值时的分子数。C T,X是目标分子扩增达到阈值时的
循环数。K
x
是一个常数。
对于内参反应而言，也有同样的公式：
用X
T 除以R
T
得到：
对于使用 Taqman 探针的实时扩增而言， X
T 和 R
T
的值由一系列因素决定：包
括探针所带的荧光报导基团、探针序列对探针荧光特性的影响、探针的水解效率和纯度以及荧光阈值的设定。因此常数K 并不一定等于1 。
假设目标序列与内参序列扩增效率相同：
或：
X N 代表经过均一化处理过的初始目标分子量；△C T 表示目标基因和内标基因C T 值的差异（C T,X －C T,R ）整理上式得：
最后用任一样本q 的X N 除以参照因子（calibrator, cb ）的X N 得到：
在这里
对于一个少于150bp 的扩增片断而言，如果 Mg 2+ 浓度、引物都进行了适当的优化，扩增效率接近于1 。因此目标序列的量通过内均一化处理之后相对于参照因子而言就是
1.2. 2－△△CT 方法的假设和应用 要使△△CT 计算方法有效，目标序列和内参序列的扩增效率必须相等。看两个反应是否具有相同的扩增效率的方法是看他们模板浓度梯度稀释後扩增产物△CT 如何变化。
图1 显示的是 cDNA 样品在 100 倍稀释范围内的实验结果。对于每一个稀释样本，都用GAPDH 和c-myc 特异的荧光探针及引物进行扩增。计算出c-myc 和 GAPDH 的平均C T值以及△C T值，通过cDNA 浓度梯度的log值对△C T值作图，如果所得直线斜率绝对值接近于0，说明目标基因和内标基因的扩增效率相同，就可以通过△△CT方法进行相对定量。在图1中，直线斜率是0.047，因而假设成立，△△C T方法可以用来分析数据。如果两个扩增反应效率不同，则需要通过定量标准曲线和绝对定量的方法来
进行相对定量；或者也可以重新设计引物，优化反应条件使得目标序列和内参序列具有相同的扩增效率。
1.3. 2－△△CT内标和参照因子的选择
使用内标基因的目的是为了对加入到反转录反应中的RNA 进行均一化处理。标准的看家基因一般都可被用作内标基因。适合于实时 PCR 反应内标基因包括GAPDH， -actin，2-microglobulin 以及rRNA 。当然，其它的措组词语
看家基因也同样能被用作内标。我们推荐在应用某一基因作为内标之前首先确证该基因的表达不会受实验处理的影响。验证实验处理是否对内标基因表达产生影响的方法在2.2 部分有描述。
2－△△CT方法中参照因子的选择决定于基因表达定量实验的类型。最简单的设计就是把未经处理的爱国名人
样品作为参照因子（calibrator ）。经内标基因均一化处理後，通过方法计算，目标基因表达差异通过经过处理的样本相对于未经处理的样本的倍数来表示。对于未经处理的参照样，△△C T=0，而20=1。所以根据定义，未处理样本的倍数变化讲座主持稿
为1 。而对于那些经过处理的样本，相对于参考因子基因表达的倍数为2－△△CT。同样的分析也可用于不同时相的基因表达差异，在这种情况下，一般选0 时刻的样本作为参照因子。
有些情况下，并不是比较不同处理样本基因表达差异。例如，有的是想看某一器官中特定 mRNA 的表
达。在这种情况下，参照因子可能是另一器官中该 mRNA 的表达。表1显示了大脑和肾脏总RNA 中c-myc 和GAPDH 转录本的CT 值。在这一个例子中，大脑被人为的选择为参照因子，通过计算得到肾脏c-myc 表达量经GAPDH 校正後相对于大脑的表达量的结果。尽管相对定量方法可用于这种组织之间的比较，但结果的生物学解释是相当复杂的。不同种类细胞中目标和参照转录本单一的相对量变化可能在任何特定的组织中都存在。
1.4. 2－△△CT方法的数据分析
实时定量PCR所得到CT值可以很容易的输出到表格程序如M自动化英文
icrosoft Excel中去。为了显示数据分析过程，我们在这里给出了一个基因表达定量的实验数据和样本列表。通过-actin 均一化处理，我们对目标基因fos-glo-myc 的表达变化进行了监测。在8h 的时间范围内，在每一时间点都取3 个重复样本，每一样本在cDNA 合成之後都做定量PCR ，数据分析用到了公式9，即：
Time x 表示任意时间点，Time 0表示经 -actin 校正后1 倍量的目标基因表达。
0 时刻目标基因和内标基因的平均C T（见图2 第8 栏）被用于公式9 中。通过公式9 计算出每一个样本目标基因表达通过 -actin 均一化处理後相对于0 时刻的倍数（见图2 第9 栏）。平均SD，CV 由每一个时间点所取的三个重复样求得。用这种分析方法，在0 时刻的平均倍数变化接近于1。我们发现通过检测在 0 时刻平均倍数变化是否为1 可以很方便的验证三个重复样品之间是否有错误或者误差。如果得到的结果与1偏差很大，则表明存在计算错误或者是很高的实验误差。
在前面的例子中，在每一时间点上分别取了三个独立的 RNA 样本进行了分析。因此对每一个样本分别处理，通过计算後取结果的平均值就非常重要。如果是同一样本进行 PCR 扩增的重复，这就需要首先求出平均CT，然後再进行计算。怎么样计算平均值就要看目标基因和内参基因是在同一个管子中扩增还是在不同的管子中扩增。表 1 给出了目标基因（c- myc ）和内参基因（GAPDH ）在不同管中扩增的实验数据。在这里不应该把任一单个的 c-myc 管子和 GAPDH 管子作比较，而应该分别计算出 c-myc 和 GAPDH 的平均C T来计算△C T。重复实验中C T值的估计偏差通过标准的指数计算转化成最後结果中相对量的变化。但其中的一个难点是 C T值与相应的拷贝数成指数关系（见第4 部分），因此，在最後的计算中，的误差通过△△C T 加上标准偏差和△△C T减去标准偏差来评估，这就使得求得的数值相对于平均值呈不对称分布。不对称分布是因为结果经指数处理後转化成量的线性比较造成的。
通过不同荧光染料标记的探针，我火灾逃生演练
们可以在同一管中同时扩增目标序列和内标序列。表2 给出了目标基因（c- myc）和内标基因（GAPD H）在同一管中扩增的实验数据。对于任意一个管子，目标基因（c- myc）和内参基因（GAPDH）扩增时加入的cDNA 量都是一样多的，所以可以分别对每个管子计算△C T值，这些值取平均后再进行计算。在这里估计误差值也是一个不对称的范围，反映了误差经指数处理转化为线性差异。
