qRT-PCR实验及数据分析

qRT-PCR实验原理

RT-qPCR由普通PCR技术发展而来,它是在传统PCR反应体系中加入荧光化学物质(荧光染料或者荧光探针),根据各自不同的发光机制实时检测PCR退火、延伸过程中荧光信号变化来计算PCR每个循环中产物的变化量,目前最常见的方法为荧光染料法和探针法。

荧光染料法:

一些荧光染料如SYBR Green Ⅰ,PicoGreen,BEBO等,它们本身不发光,但结合于dsDNA的小沟后会发出荧光。所以当PCR反应刚开始时机器并不能检测到荧光信号,当反应进行到退火-延伸(二步法)或者延伸阶段(三步法),此时双链打开,在DNA聚合酶的作用下新链合成,荧光分子就结合于dsDNA的小沟中并发出荧光,随着PCR循环数的增加越来越多的染料与dsDNA结合,荧光信号也不断的增强。以SYBR Green Ⅰ为例。

探针法:

Taqman探针是最为常用的一种水解探针,在探针的5’端存在一个荧光基团,通常为FAM,探针本身则为一段与目的基因互补的序列,在探针的3’端有一个荧光猝灭基团,根据荧光共振能量转移原理(Förster resonance energy transfer, FRET),当报告荧光基团(供体荧光分子)和猝灭荧光基团(受体荧光分子)激发光谱重叠且距离很近时(7-10nm),供体分子的激发可以诱发受体分子发荧光,而自身荧光减弱。所以PCR反应开始,探针游离于体系中完整存在时,报告荧光基团并不会发出荧光,当退火时,引物和探针结合于模板,在延伸阶段,聚合酶不断的合成新链,由于DNA聚合酶具有5’-3’核酸外切酶活性,到达探针时,DNA聚合酶就会将探针从模板上水解下来,报告荧光基团和猝灭荧光基团分开,释放荧光信号。由于探针和模板存在一对一的关系,所以在试验的精度和灵敏度上,探针法都要优于染料法。

Fig 1 qRT-PCR原理



引物设计

原则:

引物应用核酸系列保守区内设计并具有特异性。最好用cDNA序列,mRNA序列也可,如果没有则找出DNA序列的cds区域设计。

做荧光定量产物长度80-150bp最好,最长是300bp

引物长度一般在17-25碱基之间,上下游引物不能相差太大。

G+C含量在40%~60%之间,45-55%最佳。

TM值在58-62度之间。

尽量避免引物二聚体和自二聚体,(不要出现超过4对连续互补的碱基)发卡结构,如果不可避免,使ΔG<4.5kJ/mol*

如果不能确保反转录过程中gDNA已经除干净,最好设计出夸内含子的引物*

3′端不可修饰,而且要避开AT,GC rich的区域,避开T/C,A/G连续结构(2-3个)

引物与非特异性扩增序列的同源性最好小于70%或者有8个互补碱基同源。

数据库:

CottonFGD search by keywords

引物设计:

IDT-qPCR primer design


Fig2 IDT在线引物设计工具页面


Fig3 结果页面展示

ncRNA引物设计:

lncRNA: 和mRNA的步骤一致。

miRNA:参考这个方法:茎环法miRNA引物设计-丁香园 茎环法原理:由于miRNA全部为23nt左右的短序列,无法进行直接PCR检测,所以使用了茎环序列这个工具,茎环序列是一段50nt左右单链DNA,可以自身形成发卡结构,3’末端可以设计成与miRNA部分片段互补的序列,那么在反转录时就可以将目的miRNA连接到茎环序列上,那么总长度就达到70bp,就符合qPCR确定的扩增产物的长度。加尾法miRNA引物设计



天根反转录引物设计


扩增特异性检测:

Online blast database:CottonFGD blast by sequence similarity

Local blast:参照我之前写的用Blast+做本地blast,linux和macos可以直接建立本地数据库,win10系统需要安装ubuntu bash后也可以做。本地blast数据库建立及本地blastwin10开启ubuntu bash

Notice:陆地棉海岛棉为四倍体作物,所以blast结果往往会是两个或多个匹配,以往经验以NAU的cds作为database进行blast很可能找到两个只有几个SNP差异的同源基因,通常通过引物设计是分不开这两个同源基因的,所以就将他们当做同一个来做,如果有明显的indel,通常把引物设计在indel上,但这样就可能导致引物的二级结构的自由能变高,导致扩增效率的下降,但这是不可避免的。

引物二级结构的检测:

步骤:  打开oligo 7→输入模版序列→关闭子窗口→保存→将引物定位到模版上,按ctrl+D设定引物长度→分析各种二级结构,例如自二聚体,异二聚体,发卡,错配等。图4中最后两张图片为引物测试结果,前引物结果较好,没有明显的二聚体和发卡结构,没有连续的互补碱基,且自由能绝对值也小于4.5,而后引物则出现了连续6碱基互补,自由能为8.8;另外在3’端也出现了较为严重的二聚体话,出现了连续4碱基的二聚体,虽然自由能不高,但3’的二聚体化可以严重的影响扩增特异性和扩增效率。此外还需要对发卡及异二聚体,错配进行检查。


Fig3 oligo7 检测结果

扩增效率检测:

PCR反应的扩增效率严重的影响到PCR的结果,同样在qRT-PCR中,扩增效率对定量的结果尤为重要,除去反应液(reaction buffer)中其他物质及机器和protocol带来的影响,引物的质量对于qRT-PCR的扩增效率影响也非常大,为了保证结果的准确性,无论是相对荧光定量还是绝对荧光定量均需要对引物的扩增效率进行检测,而公认的有效的qRT-PCR的扩增效率为85%-115%之间有两种方法:

1.标准曲线法: 

a.将cDNA混合

b.梯度稀释

c.qPCR

d.线性回归方程计算扩增效率

2.LinRegPCR

LinRegPCR is a program for the analysis of real time RT-PCR Data,also called quantitative PCR(qPCR)data based on SYBR Green or similar chemistry. The program uses non-baseline corrected data,performs a baseline correction on each sample separately,determines a window-of-linearity and then uses linear regression analysis to fit a straight line through the PCR data set. From the slope of this line the PCR efficiency of each individual sample is calculated. The mean PCR efficiency per amplicon and the Ct value per sample are used to calculate a starting concentration per sample,expressed in arbitrary fluorescence units. Data input and output are through an Excel spreadsheet.

只需要混样,不需要做梯度

步骤: (以伯乐CFX96为例,不太清楚ABI的机器)

实验:为标准qPCR实验。

qPCR 数据输出 : LinRegPCR 可以识别两种形式的输出文件形式:RDML或者quantification Amplification result.实际上就是机器对循环数和荧光信号的实时检测值,通过对线性区段的荧光变化值分析得出扩增效率。

数据选择:理论上RDML值应该是可以用的,估计是我电脑的问题软件并不能识别RDML,所以我已excel输出值作为原始数据,建议先对数据进行粗筛选,例如加样失败的等造成的点可以在输出的数据中删掉(当然不删也可以,LinRegPCR后期会对这些点进行忽略)

Fig5 qPCR data export


Fig6 selection of candidate samples

数据输入:打开 qualification amplification results. xls, →打开 LinRegPCR →文件→ read from excel → 如图7选择参数→ OK → 点击determine baselines


Fig7 steps of linRegPCR data input

结果:如果没有做重复则不需要分组,如果做了重复在sample grouping中可以编辑分组,indentifier 中键入基因的名称,然后相同的会自动分为一组。最后,点击文件,输出excel,查看结果,每个孔的扩增效率及R2结果都会显示,其次如果你分了组,会显示校正过的扩增效率均值。保证每个引物的扩增效率位于85%-115%之间,过大或者过小都说明引物扩增效率较差。


Fig8 Result and data output



实验过程:

RNA的质量要求:

纯度:1.7 <OD260/OD280<2.0(<1.7时表明有蛋白质或酚污染;>2.0时表明可能有异硫氰酸残存。 干净的核酸A260/A230应该在2左右。如果有在230 nm 处的强吸收提示污染有酚盐离子等有机化合物。具体吸光值和纯度的关系看这篇文章:OD260/OD280比值的问题-丁香通;此外可以用1.5%的琼脂糖凝胶电泳检测,具体方法:RNA琼脂糖凝胶电泳-丁香通,要注意的是如果只检测纯度是否准确用非变性胶检测,不要点marker,因为ssRNA没有变性的情况下RNA的涌动率和分子量对数不程线性关系,无法正确表示分子量。

浓度:理论上不低于100ng/ul,如果浓度太低,一般情况下纯度也不高

Fig9 RNA gel

此外,如果样品比较珍贵,且RNA浓度较高时建议提完后可以分装,将RNA稀释到终浓度100-300ng/ul做反转录。

反转录过程

在进行mRNA转录时,采用可以特异结合到polyA尾巴的oligo(dt)引物反转录,而lncRNA,circRNA则采用随即六聚体(Random 6 mer)引物对总RNA进行反转录,对于miRNA则使用miRNA特异的颈环引物进行反转录,很多公司现在推出了专门的加尾法试剂盒,对于茎环法而言,加尾法更加便捷,高通量,节省试剂,但是对于同家族的miRNA区分效果应该没有茎环法好。对于基因特异性引物(茎环)的浓度各反转录试剂盒都有要求。miRNA现在用的内参为U6,在茎环法反转的过程中要单独反转录一管U6,直接加入U6的前后引物。而circRNA和lncRNA均可以用HKGs作为内参。


cDNA检测

一般情况下RNA如果没有问题,cDNA也应该没有问题,但如果追求实验的完美性,最好用一个可以区分gDNA和cds的内参基因(Reference gene,RG)一般RG为管家基因(housekeeping gene, HKG)如图10;当时在做大豆贮藏蛋白,用了包含内含子的actin7作为内参,该引物在gDNA中的扩增片段大小为452bp,cDNA为模板的话为142bp,那么检测结果发现有部分cDNA其实受到了gDNA的污染,同时也证明了反转录的结果没有问题,完全可以作为PCR的模板。直接用cDNA跑琼脂糖凝胶电泳没有用,跑出来是弥散的带,说服力不强。


Fig 10 cDNA检测


qPCR条件的确定

一般来说按照试剂盒的protocol来没有问题,主要是在tm值这一步,如果在引物设计时部分引物非常设计,导致tm值与理论的60℃上下差异较大,建议将cDNA样品混合后用引物跑一个梯度PCR,尽量避免将没有条带的温度设为TM值。



数据分析

常规的相对荧光定量PCR的处理方法基本是按照2-ΔΔCT.数据处理模板

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 206,126评论 6 481
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 88,254评论 2 382
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 152,445评论 0 341
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 55,185评论 1 278
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 64,178评论 5 371
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,970评论 1 284
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,276评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,927评论 0 259
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 43,400评论 1 300
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,883评论 2 323
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,997评论 1 333
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,646评论 4 322
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,213评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,204评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,423评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,423评论 2 352
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,722评论 2 345

推荐阅读更多精彩内容

  • PCR技术完全讲座 PCR技术概论 聚合酶链反应(Polymerase Chain Reaction,PCR)是8...
    三豆同学阅读 2,832评论 0 14
  • 基因表达 什么是基因表达,如下是来自于维基百科的解释: 研究方法 定量PCR 这部分我不太懂,所以就放几段百度百科...
    xuzhougeng阅读 120,865评论 15 242
  • 通常来讲,real-time qPCR的反应程序不需要像常规的PCR那样,要变性、退火、延伸3步。由于其产物长度在...
    toneker阅读 1,783评论 0 5
  • 每一个结束,其实就是下一个开始。所以,哪有什么结束,只是另外一个开始罢了。 5月始,去衡山看日。 层层叠叠,掩掩映...
    serendipityw阅读 385评论 5 0
  • 相遇总是那么猝不及防,缘分出现的刚刚好最美人心。 有时候,童话般的美好也会出现在现实中,那么,我是幸运的。你,...
    囧囧囧景江style阅读 1,281评论 5 4