利用实时定量PCR和2
－△△CT
法分析基因相对表达量
METHODS25,402–408(2001)
AnalysisofRelativeGeneExpressionDataUsingReal-TimeQuantitativePCRandthe2
－△△
CT
Method
KennethJ.LivakandThomasD.Schmittgen
?,1
AppliedBiosystems,FosterCity,California94404;and?DepartmentofPharmaceuticalSciences,
CollegeofPharmacy,
WashingtonStateUniversity,Pullman,Washington99164-6534
摘要：
现在最常用的两种分析实时定量PCR实验数据的方法是绝对定量和相对定量。绝对定量通过标准曲
线计算起始模板的拷贝数；相对定量方法则是比较经过处理的样品和未经处理的样品目标转录本之间的表
达差异。2
－△△CT
方法是实时定量PCR实验中分析基因表达相对变化的一种简便方法，即相对定量的一种
简便方法。本文介绍了该方法的推导，假设及其应用。另外，在本文中我们还介绍了两种2
－△△CT
衍生方法
的推导和应用，它们在实时定量PCR数据分析中可能会被用到。
关键词：反转录PCR定量PCR相对定量实时PCRTaqman
反转录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
公式的推导，以及实验设计，有效性评估
在AppliedBiosystemsUserBulletinNo.2（P/N4303859）中有介绍。用2
－△△CT
方法分析基因表达数据在
文献中也有报道（5,6）。本文介绍了该方法的推导、假设以及应用。另外，本文还介绍了2
－△△CT
两种衍
生方法的推导和应用，它们在实时定量PCR数据分析中都可能被用到。
1.2
－△△CT
方法
1.1.2
－△△CT
方法的推导
PCR指数扩增的公式是：
这里，X
n
是第n个循环後目标分子数，X
0
是初始目标分子数，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，因而假设成立，△
△CT方法可以用来分析数据。如果两个扩增反应效率不同，则需要通过定量标准曲线和绝对定量的方法来
进行相对定量；或者也可以重新设计引物，优化反应条件使得目标序列和内参序列具有相同的扩增效率。
1.3.2
－△△CT
内标和参照因子的选择
使用内标基因的目的是为了对加入到反转录反应中的RNA进行均一化处理。标准的看家基因一般都
可被用作内标基因。适合于实时PCR反应内标基因包括GAPDH，β-actin，β
2
-microglobulin以及rRNA。
当然，其它的看家基因也同样能被用作内标。我们推荐在应用某一基因作为内标之前首先确证该基因的表
达不会受实验处理的影响。验证实验处理是否对内标基因表达产生影响的方法在2.2部分有描述。
2
－△△CT
方法中参照因子的选择决定于基因表达定量实验的类型。最简单的设计就是把未经处理的样
品作为参照因子（calibrator）。经内标基因均一化处理後，通过方法计算，目标基因表达差异通过经过处
理的样本相对于未经处理的样本的倍数来表示。对于未经处理的参照样，△△C
T
=0，而2
0
=1。所以根据定
义，未处理样本的倍数变化为1。而对于那些经过处理的样本，相对于参考因子基因表达的倍数为2
－△△CT
。
同样的分析也可用于不同时相的基因表达差异，在这种情况下，一般选0时刻的样本作为参照因子。
有些情况下，并不是比较不同处理样本基因表达差异。例如，有的是想看某一器官中特定mRNA的
表达。在这种情况下，参照因子可能是另一器官中该mRNA的表达。表1显示了大脑和肾脏总RNA中
c-myc和GAPDH转录本的CT值。在这一个例子中，大脑被人为的选择为参照因子，通过计算得到肾脏
c-myc表达量经GAPDH校正後相对于大脑的表达量的结果。尽管相对定量方法可用于这种组织之间的比
较，但结果的生物学解释是相当复杂的。不同种类细胞中目标和参照转录本单一的相对量变化可能在任何
特定的组织中都存在。
1.4.2
－△△CT
方法的数据分析
实时定量PCR所得到CT值可以很容易的输出到表格程序如MicrosoftExcel中去。为了显示数据分
析过程，我们在这里给出了一个基因表达定量的实验数据和样本列表。通过β-actin均一化处理，我们对
目标基因fos-glo-myc的表达变化进行了监测。在8h的时间范围内，在每一时间点都取3个重复样本，
每一样本在cDNA合成之後都做定量PCR，数据分析用到了公式9，即：
Timex表示任意时间点，Time0表示经β-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）和内标基因（GAPDH）在同一管中扩增的实验数据。对于任意一个管子，目标基因（c-
myc）和内参基因（GAPDH）扩增时加入的cDNA量都是一样多的，所以可以分别对每个管子计算△C
T
值，
这些值取平均后再进行计算。在这里估计误差值也是一个不对称的范围，反映了误差经指数处理转化为线
性差异。
在表1和表2中，估计误差在从△C
T
到△△C
T
的计算中未见有增加，这是因为我们把参照基因和
检测基因的误差都显示出来了。我们把△CT，cb当作一个人为设定的常数来减去，得到△△CT。这样
得到的结果就与图2所显示的在求平均之前对不同重复样本分别通过各自的CT值求实所得结果相
当。另一种方法是将参照基因当作没有任何误差的1倍的量，在这种情况下，平均△C
T,cb
的误差值被引入
到每一样本的△△C
T
中。在表1中，肾脏中△△C
T
变成-2.50±0.20而经过校正的c-myc量是5.6倍，
范围从4.9到5.6。而在大脑中的结果是没有误差的1倍。
2.2
－△Cf
方法
2.1.2
－△Cf
方法的推导
通过内标RNA可以对加入RNA的量的差异进行校正。2
－△△CT
方法的数据分析的一个特点就是能
够利用实时PCR实验的一部分数据来完成这种校正。在其它的方法不能定量初始RNA量的时候：例如，
在能得到的RNA量非常有限的时候或者需要处理高通量的样品的时候，这一方法的优势就格外明显。当
然我们也可以利用PCR实验以外的方法来完成这种校正。最常用的一种方法就是用紫外吸收来确定用于
cDNA合成的RNA量，然後将相同的RNA反转录产生的cDNA用于PCR定量反应。这种外标法校正
的一个应用例子就是研究某种实验处理是否影响内标基因的表达。在这里，目标基因和内标合二为一。在
这个例子中，公式[2]不被公式[3]除，公式[5]变成：
整理得：
任一样品X
0,q
除以参照品X
0,cb
得：
在这里△C''
T
=C
T
,
q
-C
T,cb
。△C’
T
与前面计算中用的△C
T
（用目标基因C
T
值减去参照基因C
T
值）相
互区别。
就象在1.1部分所描述的，如果条件优化较好，效率接近于1，内标相对于参照因子为：
2.2.2
－△Cf
方法的应用
2
－△CT''
方法的一个应用就是确定实验处理对某一候选内标基因的影响。为了显示这一过程，我们做了
血清饥饿/诱导实验（7）。血清饥饿/诱导是研究某些mRNA降解的常用方法（8）。然而，血清可能
影响一些基因的表达包括标准的看家基因的表达。
在24-h血清饥饿培养之後，在NIH3T3细胞中加入15%血清诱导基因表达。从细胞中提取Poly(A)
+
RNA，并将之反转录成cDNA。利用SYBRGreen通过实时定量PCR检测GAPDH，β2-microglobulin
cDNA的量。GAPDH和β
2
-microglobulin各自的相对量通过2
－△CT''
公式求得。细胞处理对于GAPDH的
基因表达有明显影响，但对β
2
-microglobulin没有什么影响。因此β
2
-microglobulin很适合做血清刺激定
量实验的内标，而GAPDH并不适合。这一例子向大家展示了在只研究一个基因的时候怎么用2
－△CT''
的方
法分析基因相对表达数据。
3.实时PCR数据的统计学分析
实时PCR最终分析的是阈值循环或C
T
。
C
T
值通过PCR信号的对数值和循环数来确定。因此C
T
值是一个指数而非线性概念。因此，在任何
统计分析中都不要用原始的C
T
值来表示结果。正如我们在前文中所描述的一样，PCR相对量通常和内标
和参照样本一起计算而很少直接用C
T
值来表示，除非我们想检验重复样本之间的差别。为了向大家显示这
一点，我们用SYBRGreen通过real-timePCR来检测相同cDNA的96个重复反应。所有反应组分在
同一管中混好後分装到96个管中，做实时PCR分析，得到了每一个样本的CT值。为了比较样品间
变化，计算了96个样本的平均±SD，如果通过原始C
T
值计算，平均±SD是20.00±0.193，CV为0.971%。
但是如果把原始C
T
值用2-CT转化成线性形式，平均±SD是9.08×10
-7
±1.33×10
-7
，CV为13.5%。
从这个简单的例子我们可以看出，通过原始C
T
值来反映变化是错误的，应该避免。用2
－CT
将单个数据转
化成线性形式来说明重复样本之间的变化和差异更准确可靠。
4.结论
实时定量PCR实验设计和数据分析可以采用相对定量和绝对定量两种方法，研究人员在设计实时定
量PCR实验分析基因表达的时候首先要问的一个问题就是：数据最後会以一个什么样的形式得到。如果
需要知道绝对的拷贝数，就必须用绝对定量的方法，否则只需要给出基因表达相对量就足够了。相对定量
可能比绝对定量要更容易一些，因为它不需要作标准曲线。
本文所给出的公式对于每个用相对定量的方法分析基因表达差异的研究人员都足够了。下面，我们
总结一下实验设计和评估中的一些重要步骤：
选择一个内标基因。
确定内标的有效性，确保它不会受到实验处理的影响。
通过PCR扩增目标基因和内标基因RNA或cDNA的一系列梯度稀释模板确保它们的扩增效率相
同。
最後通过2
－△CT
计算将统计数据转化成线性形式而不是原始C
T
值。
参考文献
1.Murphy,L.D.,Herzog,C.E.,Rudick,J.B.,Fojo,A.T.,andBates,S.E.(1990)Biochemistry29,
10351–10356.
2.Noonan,K.E.,Beck,C.,Holzmayer,T.A.,Chin,J.E.,Wunder,J.S.,Andrulis,I.L.,Gazdar,A.
F.,Willman,C.L.,Griffith,B.,Von-Hoff,D.D.,andRobinson,I.B.(1990)Proc.Natl.Acad.Sci.toalinear
formusing22CTmoreaccuratelydepictstheUSA87,7160–7164.
3.Horikoshi,T.,etal.(1992)CancerRes.52,108–116.
4.Heid,C.A.,Stevens,J.,Livak,K.J.,andWilliams,P.M.(1996)GenomeRes.6,986–994.
5.Winer,J.,Jung,C.K.,Shackel,I.,andWilliams,P.M.(1999)Anal.Biochem.270,41–9.
6.Schmittgen,T.D.,Zakrajsek,B.A.,Mills,A.G.,Gorn,V.,Singer,M.J.,andReed,M.W.(2000)
Anal.Biochem.285,194–204.
7.Schmittgen,T.D.,andZakrajsek,B.A.(2000)J.Biochem.Biophys.Methods46,69–81.
8.Chen,C.Y.,andShyu,A.B.(1994)Mol.Cell.Biol.14,8471–8482.
9.Iyer,V.R.,etal.(1999)Science283,83–87.
10.Giulietti,A.,Overbergh,L.,Valckx,D.,Decallone,B.,Bouillon,R.,andMathieu,C.(2001)
Methods25,386–401.
11.Niesters,H.G.M.(2001)Methods25,419–429.
|
|
