牛春艷,楊佳怡,王 晶,李 晶,劉文軍
(1.中國(guó)計(jì)量科學(xué)研究院,北京 100029;2.中國(guó)科學(xué)院微生物研究所,北京 100101)
豬繁殖與呼吸綜合征(porcine reproductive and respiratory syndrome,PRRS),又稱豬藍(lán)耳病,是近年來(lái)造成養(yǎng)豬產(chǎn)業(yè)嚴(yán)重經(jīng)濟(jì)損失的高度傳染性疾病[1],最初在1980年左右于北美及歐洲被報(bào)道[2]。該疾病的病原體豬繁殖與呼吸綜合征病毒(Porcine reproductive and respiratory syndrome virus,PRRSV)為單股正鏈RNA病毒,隸屬套式病毒目動(dòng)脈炎病毒科動(dòng)脈炎病毒屬。根據(jù)血清學(xué)和遺傳特性的差異,PRRSV可分為2個(gè)型,即歐洲型(LV株為代表株)和美洲型(VR2332株為代表株)[3-4]。我國(guó)分離得到的毒株多為美洲型。2006年,我國(guó)多地區(qū)出現(xiàn)高致病性豬繁殖與呼吸綜合征病毒(Highly pathogenic PRRSV,HP-PRRSV)[5-6],給養(yǎng)豬業(yè)造成巨大經(jīng)濟(jì)損失?;蚪M測(cè)序結(jié)果顯示,相比于典型豬繁殖與呼吸綜合征病毒,HP-PRRSV(變異株)在病毒GP5及Nsp2蛋白區(qū)域存在氨基酸突變及缺失[7]。
PRRSV的檢測(cè)方法最初為病毒培養(yǎng)法,但是此方法耗時(shí)長(zhǎng)且診斷敏感性低[8]。近年來(lái)發(fā)展的巢式PCR、實(shí)時(shí)熒光定量PCR(quantitative polymerase chain reaction,qPCR)[9-10]、環(huán)介導(dǎo)等溫?cái)U(kuò)增(loop-mediated isothermal amplification,LAMP)技術(shù)等方法極大縮短了檢測(cè)時(shí)間,提高了檢測(cè)敏感性。但是,qPCR對(duì)樣本的定量分析依賴于標(biāo)準(zhǔn)曲線。數(shù)字PCR方法是最近發(fā)展起來(lái)的不依賴于外標(biāo)對(duì)核酸分子進(jìn)行定量的技術(shù)。由于數(shù)字PCR具有更高的靈敏度和準(zhǔn)確性,已在微量核酸樣本檢測(cè)、稀有突變檢測(cè)等方面有廣闊應(yīng)用前景,并受到越來(lái)越多的關(guān)注。對(duì)于RNA分子的定量檢測(cè),通常需要在逆轉(zhuǎn)錄酶的作用下,將RNA分子逆轉(zhuǎn)錄得到cDNA,從而進(jìn)行擴(kuò)增定量。因此,逆轉(zhuǎn)錄過(guò)程是RNA分子定量檢測(cè)的關(guān)鍵步驟。本研究利用建立的逆轉(zhuǎn)錄數(shù)字PCR定量檢測(cè)方法,分析了不同逆轉(zhuǎn)錄試劑盒產(chǎn)品對(duì)PRRSV定量結(jié)果的影響。
1.1 病毒和試劑 HP-PRRSV毒株HuN4及PRRSV毒株CH-1a為中國(guó)農(nóng)業(yè)科學(xué)院哈爾濱獸醫(yī)研究所提供。病毒RNA提取試劑盒(QIAamp Viral RNA Mini Kit)及QuantiTect Reverse Transcription Kit購(gòu)自QIAGEN公司;ProtoScript Ⅱ cDNA 第一鏈合成試劑盒購(gòu)自NEB公司;Reverse Transcription System購(gòu)自Promega公司;PrimeScript RT reagent Kit with gDNA Eraser購(gòu)自Takara公司; SuperScriptTMIV First-Strand Synthesis System 購(gòu)自Invitrogen公司;TaqManTMGene Expression Master Mix購(gòu)自Applied Biosystem(AB)公司;微滴式數(shù)字PCR反應(yīng)液SuperMix for Probes、微滴發(fā)生專用油、微滴分析專用油購(gòu)自Bio-Rad公司。
1.2 引物設(shè)計(jì) 將GenBank數(shù)據(jù)庫(kù)中PRRSV基因序列進(jìn)行比對(duì),選擇序列相對(duì)保守區(qū)域ORF6,設(shè)計(jì)引物及探針。F1:5'-CTAGGCCGCAAGTACATTCT-3';R1:5'- GACGACAAATGCGTGGTTATC-3';P1:5'- FAM-ATTTGCCGCAATCGGATGAAAGC C-BHQ1-3'。所有引物及探針由生工生物工程(上海)股份有限公司合成。
1.3 RNA提取及逆轉(zhuǎn)錄 按照QIAamp Viral RNA Mini Kit(QIAGEN)試劑盒使用說(shuō)明書,進(jìn)行病毒基因組RNA提取。RNA保存于-80℃。按逆轉(zhuǎn)錄試劑盒說(shuō)明書,將6 μL所提取的RNA在20 μL的反應(yīng)體系中反轉(zhuǎn)錄形成cDNA。cDNA保存于-20℃。
1.4 熒光定量PCR實(shí)驗(yàn) 以逆轉(zhuǎn)錄得到的cDNA為模板進(jìn)行熒光定量PCR實(shí)驗(yàn)。實(shí)驗(yàn)反應(yīng)體系:TaqManTMGene Expression Master Mix(AB)12.5 μL,上、下游引物各1 μL(終濃度200 nmol/L),探針1 μL(終濃度400 nmol/L),cDNA模板2 μL,加ddH2O補(bǔ)足終體積至25 μL。擴(kuò)增條件為:95℃預(yù)變性10 min;95℃ 變性30 s,60℃退火1 min,40個(gè)循環(huán);98℃變性10 min。用包含PRRSV ORF6及ORF7基因序列的線性化質(zhì)粒(中國(guó)計(jì)量科學(xué)研究院制備并定值)繪制標(biāo)準(zhǔn)曲線,對(duì)不同濃度的PRRSV cDNA進(jìn)行定量檢測(cè)。
1.5 數(shù)字PCR實(shí)驗(yàn) 以逆轉(zhuǎn)錄得到的cDNA為模板進(jìn)行數(shù)字PCR實(shí)驗(yàn)。反應(yīng)體系為:SuperMix for Probes(Bio-rad)10 μL,上、下游引物各0.8 μL(終濃度400 nmol/L),探針0.8 μL(終濃度240 nmol/L),cDNA模板4 μL,加ddH2O補(bǔ)足終體積至20 μL。將20 μL反應(yīng)混合液分裝至對(duì)應(yīng)微滴生成卡的sample孔位中,在微滴生成卡的Oil孔位中加入70 μL 微滴發(fā)生專用油,進(jìn)行微滴生成。將生成的微滴(40 μL)轉(zhuǎn)入96孔深孔板中進(jìn)行PCR反應(yīng)。擴(kuò)增條件為:95℃預(yù)變性10 min;95℃ 變性30 s,54℃退火1 min,40個(gè)循環(huán);98℃ 變性10 min。數(shù)字PCR反應(yīng)擴(kuò)增結(jié)束后,利用Droplet Reader進(jìn)行結(jié)果的讀取。選擇反應(yīng)中產(chǎn)生的可接受的平均微滴數(shù)目大于10 000的結(jié)果用于分析,保證定量分析的準(zhǔn)確性。
1.6 數(shù)字PCR反應(yīng)線性關(guān)系及重復(fù)性研究 將逆轉(zhuǎn)錄得到的cDNA進(jìn)行梯度稀釋,每個(gè)濃度梯度設(shè)置6次重復(fù),進(jìn)行數(shù)字PCR反應(yīng)的線性關(guān)系研究。以6次重復(fù)實(shí)驗(yàn)的結(jié)果計(jì)算變異系數(shù)(CV),研究數(shù)字PCR反應(yīng)重復(fù)性。
1.7 逆轉(zhuǎn)錄反應(yīng)重復(fù)性及再現(xiàn)性研究 利用NEB公司逆轉(zhuǎn)錄試劑盒對(duì)RNA進(jìn)行逆轉(zhuǎn)錄,設(shè)置3次獨(dú)立重復(fù)逆轉(zhuǎn)錄實(shí)驗(yàn),將得到的cDNA稀釋后作為模板,進(jìn)行數(shù)字PCR擴(kuò)增,并在不同時(shí)間內(nèi)分別進(jìn)行獨(dú)立重復(fù)逆轉(zhuǎn)錄實(shí)驗(yàn),研究逆轉(zhuǎn)錄反應(yīng)重復(fù)性及再現(xiàn)性。
1.8 不同逆轉(zhuǎn)錄試劑盒的擴(kuò)增結(jié)果比較 使用5種逆轉(zhuǎn)錄試劑盒將RNA模板逆轉(zhuǎn)錄,設(shè)置3次獨(dú)立重復(fù)逆轉(zhuǎn)錄實(shí)驗(yàn)。將得到的cDNA稀釋后作為模板,進(jìn)行數(shù)字PCR定量分析。比較不同逆轉(zhuǎn)錄試劑盒對(duì)定量結(jié)果的影響。
2.1 數(shù)字PCR方法的參數(shù)優(yōu)化 采用熒光定量PCR對(duì)引物及探針的濃度進(jìn)行優(yōu)化,設(shè)置引物濃度為100、200、300、400、800 nmol/L,設(shè)置探針濃度為120、240、400 nmol/L。綜合考慮Cq值及熒光擴(kuò)增信號(hào),并選擇在數(shù)字PCR反應(yīng)中能將陰性擴(kuò)增信號(hào)與陽(yáng)性擴(kuò)增信號(hào)明顯分開的條件作為最優(yōu)條件。最終選擇的體系為:引物濃度400 nmol/L,探針濃度240 nmol/L。設(shè)置退火溫度梯度為50℃~60℃。根據(jù)陰性擴(kuò)增信號(hào)與陽(yáng)性擴(kuò)增信號(hào)差異最大的原則,選擇54℃作為退火條件(圖1)。
2.2 數(shù)字PCR方法的線性關(guān)系及重復(fù)性 在利用數(shù)字PCR方法檢測(cè)核酸濃度時(shí),需要將核酸濃度稀釋到數(shù)字PCR的檢測(cè)范圍內(nèi)進(jìn)行檢測(cè)。將cDNA梯度稀釋后進(jìn)行數(shù)字PCR擴(kuò)增反應(yīng),結(jié)果表明,實(shí)驗(yàn)所建立的數(shù)字PCR方法至少在5個(gè)數(shù)量級(jí)范圍內(nèi)線性關(guān)系良好,HuN4及CH-1a對(duì)應(yīng)的R2分別為0.9995及0.9999(圖2)。以6次重復(fù)實(shí)驗(yàn)的結(jié)果計(jì)算變異系數(shù),結(jié)果如表1所示,實(shí)驗(yàn)所建立的數(shù)字PCR方法在4個(gè)數(shù)量級(jí)范圍內(nèi)變異系數(shù)小于10%(1.01%~9.61%,最小變異系數(shù)為1.01%),實(shí)驗(yàn)在此范圍內(nèi)重復(fù)性良好。
圖1 數(shù)字PCR反應(yīng)退火溫度優(yōu)化Fig.1 Optimization of annealing temperature for ddPCR
2.3 逆轉(zhuǎn)錄反應(yīng)的重復(fù)性及再現(xiàn)性 由于實(shí)驗(yàn)擬考察不同逆轉(zhuǎn)錄產(chǎn)品對(duì)于數(shù)字PCR定值結(jié)果的影響,因此必須首先保證實(shí)驗(yàn)方法的可重復(fù)性與可再現(xiàn)性。從結(jié)果2.2可以看到,對(duì)于同一逆轉(zhuǎn)錄產(chǎn)物(cDNA)進(jìn)行數(shù)字PCR定量分析,最小變異系數(shù)為1.01%,證明基于數(shù)字PCR定量方法的重復(fù)性較好。通過(guò)設(shè)置不同逆轉(zhuǎn)錄獨(dú)立重復(fù)實(shí)驗(yàn),并在不同時(shí)間內(nèi)分別進(jìn)行反應(yīng),計(jì)算實(shí)驗(yàn)結(jié)果變異系數(shù),反映逆轉(zhuǎn)錄反應(yīng)重復(fù)性及再現(xiàn)性。如表2及圖3所示,三次平行反應(yīng)間的變異系數(shù)為0.69%~3.00%,三次不同時(shí)間獨(dú)立重復(fù)實(shí)驗(yàn)間的變異系數(shù)分別為2.38%及2.47%。
2.4 數(shù)字PCR與熒光定量PCR檢測(cè)結(jié)果的比較 利用10倍梯度稀釋的線性化質(zhì)粒繪制標(biāo)準(zhǔn)曲線,標(biāo)曲斜率為-3.3927,對(duì)應(yīng)擴(kuò)增效率為97%,R2為0.9993。將10倍梯度稀釋的HuN4與Ch-1a兩種病毒cDNA分子分別利用數(shù)字PCR與熒光定量PCR進(jìn)行定量檢測(cè),并將二者的檢測(cè)結(jié)果進(jìn)行相關(guān)性分析,如圖4所示,發(fā)現(xiàn)數(shù)字PCR與熒光定量PCR方法之間存在良好的線性相關(guān)性(R2分別為0.998及0.994)。
圖2 數(shù)字PCR線性關(guān)系Fig.2 Linearity of ddPCR
表1 數(shù)字PCR方法的重復(fù)性Table 1 Repeatability of ddPCR
2.5 不同逆轉(zhuǎn)錄試劑盒對(duì)PRRSV定量檢測(cè)結(jié)果的影響 利用5種不同逆轉(zhuǎn)錄試劑盒,對(duì)HuN4及CH-1a兩種病毒RNA進(jìn)行逆轉(zhuǎn)錄后,用逆轉(zhuǎn)錄數(shù)字PCR方法進(jìn)行定量檢測(cè)。從圖4可以看到,對(duì)于同一種RNA,利用5種不同逆轉(zhuǎn)錄試劑盒,最后的定量結(jié)果有所不同。利用試劑盒A進(jìn)行逆轉(zhuǎn)錄,HuN4及CH-1a兩種樣品的定量結(jié)果均為最低;利用試劑盒B、C、D、E進(jìn)行逆轉(zhuǎn)錄,將HuN4及CH-1a兩種樣品得到的定量結(jié)果分別進(jìn)行統(tǒng)計(jì)學(xué)分析,均發(fā)現(xiàn)p<0.01,因此利用這4個(gè)試劑盒進(jìn)行逆轉(zhuǎn)錄,得到的結(jié)果存在顯著性,具有統(tǒng)計(jì)學(xué)意義。本研究采用5種不同逆轉(zhuǎn)錄試劑盒對(duì)同種病毒RNA進(jìn)行逆轉(zhuǎn)錄后,采用同種逆轉(zhuǎn)錄數(shù)字PCR方法進(jìn)行定量檢測(cè),且逆轉(zhuǎn)錄數(shù)字PCR方法的線性、重復(fù)性及再現(xiàn)性均經(jīng)過(guò)確認(rèn)和驗(yàn)證,由此可見不同逆轉(zhuǎn)錄試劑盒獲得的不同結(jié)果是由逆轉(zhuǎn)錄試劑盒本身所造成的。
表2 逆轉(zhuǎn)錄反應(yīng)的重復(fù)性及再現(xiàn)性Table 2 Repeatability and reproducibility of reverse transcription process
圖3 逆轉(zhuǎn)錄反應(yīng)的重復(fù)性及再現(xiàn)性Fig.3 Repeatability and reproducibility of reverse transcription process
數(shù)字PCR的原理是將一個(gè)PCR反應(yīng)體系分配到大量微小的反應(yīng)單元中,在每個(gè)微反應(yīng)器中包含或不包含1個(gè)或多個(gè)拷貝的目標(biāo)核酸分子,進(jìn)行“單分子模板”PCR擴(kuò)增[11]。擴(kuò)增結(jié)束后,通過(guò)陽(yáng)性反應(yīng)單元(通過(guò)終點(diǎn)熒光信號(hào)判斷)的數(shù)目和統(tǒng)計(jì)學(xué)方法計(jì)算原始樣本中目標(biāo)基因的拷貝數(shù)。因此,數(shù)字PCR可以實(shí)現(xiàn)核酸分子的絕對(duì)定量。對(duì)于RNA定量,首先需要將RNA進(jìn)行逆轉(zhuǎn)錄得到cDNA,才能進(jìn)行擴(kuò)增定量。如果逆轉(zhuǎn)錄過(guò)程中沒有將RNA分子全部轉(zhuǎn)換為cDNA,則無(wú)法被準(zhǔn)確定量。因此逆轉(zhuǎn)錄的過(guò)程直接影響最后的定量結(jié)果。逆轉(zhuǎn)錄酶已有商品化試劑,或是整合了逆轉(zhuǎn)錄酶及反應(yīng)所需的緩沖液、核苷酸、引物等,形成的逆轉(zhuǎn)錄試劑盒。實(shí)驗(yàn)室在通常情況下,都是選擇一種逆轉(zhuǎn)錄酶產(chǎn)品進(jìn)行方法研究。本文通過(guò)建立的PRRSV逆轉(zhuǎn)錄數(shù)字PCR定量檢測(cè)方法,探討了不同逆轉(zhuǎn)錄試劑盒產(chǎn)品對(duì)定量結(jié)果的影響。
圖4 PRRSV數(shù)字PCR及熒光定量PCR檢測(cè)Fig. 4 Quantification of PRRSV by ddPCR and qPCR
圖5 不同逆轉(zhuǎn)錄試劑盒對(duì)PRRSV定量檢測(cè)結(jié)果的影響Fig.5 The effect of different reverse transcription kits on the quantitative results of PRRSV
針對(duì)PRRSV RNA分子的檢測(cè),本實(shí)驗(yàn)建立了逆轉(zhuǎn)錄數(shù)字PCR定量方法。首先利用逆轉(zhuǎn)錄試劑盒產(chǎn)品將RNA分子逆轉(zhuǎn)錄,將得到的cDNA作為模板進(jìn)行數(shù)字PCR擴(kuò)增??梢钥吹剑瑢?duì)于同一種逆轉(zhuǎn)錄酶得到的cDNA模板,建立的數(shù)字PCR方法線性及重復(fù)性良好,并且與建立的熒光定量PCR方法的檢測(cè)結(jié)果具有良好的線性相關(guān)性,因此可用于比較逆轉(zhuǎn)錄酶產(chǎn)品對(duì)定量結(jié)果的影響[9]。本研究共選擇了5種不同逆轉(zhuǎn)錄試劑盒進(jìn)行比較,其中利用試劑盒A進(jìn)行逆轉(zhuǎn)錄,對(duì)HuN4及CH-1a兩種菌株的RNA分子定量結(jié)果均明顯低于其他試劑盒[6]。利用其他4種試劑盒進(jìn)行逆轉(zhuǎn)錄,得到的定量結(jié)果通過(guò)統(tǒng)計(jì)分析組內(nèi)差異與組間差異,發(fā)現(xiàn)存在統(tǒng)計(jì)學(xué)差異。而且,對(duì)于HuN4及CH-1a兩種模板,不同試劑盒定量結(jié)果差異的程度有所不同。分析數(shù)字PCR定量結(jié)果存在差異原因,可能有兩個(gè)方面:一方面,可能由于不同廠家所生產(chǎn)的逆轉(zhuǎn)錄酶活力、純度等有所不同,造成目標(biāo)RNA轉(zhuǎn)錄為cDNA的數(shù)量有所不同,未被逆轉(zhuǎn)錄的目標(biāo)RNA則不會(huì)被擴(kuò)增;另一方面,可能由于不同廠家試劑盒生產(chǎn)試劑盒中加入的試劑不同,對(duì)后續(xù)PCR反應(yīng)有不同程度的抑制作用存在,造成逆轉(zhuǎn)錄得到的cDNA不能被完全擴(kuò)增定量,從而影響了后續(xù)基于數(shù)字PCR的定量檢測(cè)結(jié)果。
因此,在對(duì)RNA分子進(jìn)行定量檢測(cè)時(shí),無(wú)論擴(kuò)增反應(yīng)優(yōu)化的程度如何,逆轉(zhuǎn)錄的過(guò)程都將直接影響定量結(jié)果。在涉及逆轉(zhuǎn)錄過(guò)程的定量檢測(cè)方法建立的過(guò)程中,應(yīng)對(duì)逆轉(zhuǎn)錄產(chǎn)品進(jìn)行考察和選擇。