楊芷翊,王新凱,2,史玉婷,付思源,張鈺炘,曹琛福,賈偉新*
(1.華南農(nóng)業(yè)大學(xué)獸醫(yī)學(xué)院 國家禽流感專業(yè)實驗室(廣州)/廣東省人獸共患病防控制劑工程實驗室/人獸共患病防控制劑國家地方聯(lián)合工程實驗室/農(nóng)業(yè)農(nóng)村部人畜共患病重點實驗室/廣東省動物源性人獸共患病防控重點實驗室,廣州 510642;2.深圳海關(guān)動植物檢驗檢疫技術(shù)中心,深圳 518045)
禽流感病毒(avian influenza virus, AIV)是禽流行性感冒(禽流感)的病原,含有8個單鏈負鏈RNA片段,包含在獨立的病毒核糖核蛋白復(fù)合物(vRNPs)中[1]。根據(jù)禽流感病毒對雞的致病性,可將其分為高致病性禽流感病毒(HPAIV)和低致病性禽流感病毒(LPAIV)。HPAIV主要包括H5和H7兩種亞型,我國H5亞型HPAIV主要有H5N1、H5N6和H5N8三種亞型[2]。近年來,高致病性H5亞型禽流感疫情規(guī)模之大和地理分布之廣前所未有,不僅在亞洲、歐洲和非洲的家禽中流行,也引起了人的偶發(fā)感染[3],嚴重威脅家禽業(yè)和人類健康[4]。對H5亞型HPAIV的早期精準監(jiān)測是控制H5亞型禽流感的關(guān)鍵措施。
目前,對于H5亞型HPAIV的檢測主要包括病毒分離培養(yǎng)、抗原抗體檢測和分子生物學(xué)檢測。病毒分離培養(yǎng)法是診斷禽流感病毒感染的“金標準”[5],其檢測準確率高、檢測結(jié)果完整,但耗時長且工作量大,不適合H5亞型禽流感爆發(fā)時的大規(guī)模檢測?;诳乖贵w的AIV H5亞型檢測方法種類繁多,研究范圍廣泛,但是不可避免地都會出現(xiàn)抗原抗體檢測法本身存在的缺陷,如檢測過程中會經(jīng)常出現(xiàn)假陽性反應(yīng)從而降低檢測結(jié)果的可信度,商品化的抗原或抗體可能存在價效差異,制備抗體血清的時間較長等。相比于抗原抗體檢測,分子生物學(xué)檢測(核酸檢測)將對病原體檢測的研究提高到了基因分子水平,有著耗時短、靈敏度高、特異性強等優(yōu)勢,能夠有效地應(yīng)用于致病菌的快速檢測。
有規(guī)律成簇間隔短回文重復(fù)序列(clustered regularly interspaced short palindromic repeats, CRISPR)及其相關(guān)蛋白在RNA干擾和基因編輯等領(lǐng)域有著巨大的潛力[6]。其中CRISPR相關(guān)蛋白Cas13a和Cas12a能在crRNA的引導(dǎo)下特異性切割靶標RNA,并在切割完成后仍保持活性,繼續(xù)切割其他非靶標 RNA,即具有“附帶切割”能力[7]。基于Cas12和Cas13的旁切活性,研究團隊開發(fā)出各種核酸檢測工具。張鋒實驗室的 Gootenberg 等[8]開發(fā)了基于CRISPR-Cas13的特異性高靈敏度酶解報告基因解鎖(specific high-sensitivity enzymatic reporter unlocking,SHERLOCK)檢測系統(tǒng)。在反應(yīng)體系中加入大量信號報告RNA,該RNA一端連有熒光基團,另一端連有猝滅基團,正常情況下報告RNA不發(fā)出熒光。當(dāng)Cas13蛋白-crRNA二元復(fù)合體識別底物RNA后,Cas13蛋白被激活,除了特異性地切割底物RNA,還會非特異性地切割環(huán)境內(nèi)的報告分子,釋放信號。Jennifer Doudna團隊的Chen等[9]則基于Cas12a開發(fā)了DETECTR(DNA endonuclease targeted CRISPR trans reporter)技術(shù)。與SHERLOCK類似,Cas12a/crRNA識別檢測序列后,激活Cas12a的旁切活性,接著切割單鏈DNA報告探針并發(fā)出熒光。有研究比較了Cas13a和Cas12a酶的核酸檢測特性和報告酶切割活性,結(jié)果顯示靶向RNA的Cas13a在低活化劑濃度下的檢測速度明顯快于靶向DNA的Cas12a[10]。
但僅依靠 Cas 蛋白的切割活性檢測實際樣本存在敏感度不足、耗時長等問題[11]。目前主要的解決方法是與核酸擴增技術(shù)聯(lián)合,而等溫擴增技術(shù)(ITA)的反應(yīng)過程始終維持在恒定的溫度下,通過添加不同活性的酶和各自特異性引物來達到快速擴增核酸的目的。與其他的核酸擴增技術(shù)相比,等溫擴增有快速、高效、特異的優(yōu)點且無需專用的設(shè)備。逆轉(zhuǎn)錄重組酶輔助擴增(reverse transcription-recombinase aided amplification, RT-RAA)是一種核酸恒溫擴增技術(shù),具有靈敏度高、特異性強、檢測速度快、操作簡便等優(yōu)點。在室溫下,重組酶可以與引物DNA緊密結(jié)合形成聚合物。當(dāng)引物識別出與之完美匹配的互補序列模板DNA時,在單鏈DNA結(jié)合蛋白的幫助下,DNA聚合酶可以打開模板DNA的雙鏈結(jié)構(gòu),形成新的cDNA鏈,擴增產(chǎn)物呈指數(shù)級增長[12]。采用依賴于體溫的RAA核酸擴增技術(shù)來提高樣品中的DNA或RNA水平,從而增加CRISPR工具的靈敏度,且這種工具能夠在幾乎任何環(huán)境下使用[8]。目前,基于RT-RAA的AIV H5亞型核酸CRISPR-Cas13a檢測方法未見報道。
本研究擬運用CRISPR-Cas13a基因編輯術(shù)結(jié)合逆轉(zhuǎn)錄重組酶輔助擴增技術(shù)(RT-RAA)建立檢測 AIV H5亞型的新平臺,該類檢測平臺具有新穎、快速、方便等特點,既可以在實驗室完成大批量臨床樣本的檢測,也可以在養(yǎng)殖場等室外環(huán)境進行小規(guī)模的快速檢測。這將為H5亞型禽流感病毒的臨床檢測提供新的技術(shù)手段,也為其他RNA病毒的快速準確檢測提供新思路,具有較廣闊的應(yīng)用前景。
含有目的基因的pUC57-H5質(zhì)粒標準品購自生工生物工程(上海)股份有限公司。H3、H5、H6、H7、H9、H10亞型禽流感病毒(avian influenza virus, AIV),雞新城疫病毒(Newcastle disease virus, NDV)、雞傳染性支氣管炎病毒(infectious bronchitis virus, IBV)、雞傳染性法氏囊病病毒(infectious bursal disease virus, IBDV)和鴨坦布蘇病毒(duck tembusu virus, DTMUV)由華南農(nóng)業(yè)大學(xué)國家禽流感專業(yè)實驗室(廣州)保存。
RNA提取試劑盒購自上海飛捷生物技術(shù)有限公司,RT-RAA核酸擴增試劑(基礎(chǔ)型)購自杭州眾測生物科技有限公司,EZ-10柱式DNA純化試劑盒購自生工生物工程(上海)股份有限公司,DL1000 DNA Marker購自寶日醫(yī)生物技術(shù)(北京)有限公司,T7 RNA Polymerase和NTP Buffer Mix購自New England Biolabs,LwaCas13a蛋白購自廣州美格生物科技有限公司,Murine RNase inhibitor(MRI)購自南京諾唯贊生物科技股份有限公司。實時熒光PCR儀和Bio-rad恒溫水浴鍋購自Biod-rad公司。
1.3.1 RT-RAA引物與crRNA序列的設(shè)計與合成 從 GISAID數(shù)據(jù)庫中下載并使用342條AIV H5亞型的HA基因序列,使用DNAStar軟件對比分析序列的保守性,在保守區(qū)域內(nèi)選擇28 bp作為crRNA的目的基因靶向序列。根據(jù)RAA引物設(shè)計的原則,通過Oligo 7.0軟件設(shè)計了RT-RAA引物,設(shè)計的引物擴增區(qū)域內(nèi)包含crRNA 的目的基因靶向序列。在RT-RAA上游引物的5′端添加25 bp的T7啟動子序列,序列如表1所示。crRNA由廣州博徠斯生物科技股份有限公司合成,引物由生工生物工程(上海)股份有限公司合成。
表1 基因序列信息Table 1 Gene sequence information
1.3.2 RNA提取與RT-RAA擴增 用RNA提取試劑盒對H5病毒株樣本核酸進行提取后,按照眾測?RT-RAA核酸擴增試劑(基礎(chǔ)型)使用說明書對模板進行RT-RAA擴增。向單個干粉管中加入29.4 μL A Buffer,T7-RAA-F(10 μmol·L-1) 和RAA-R(10 μmol·L-1) 各2 μL,再加入模板和DNase/RNase-Free Water共14.1 μL,向每個反應(yīng)管蓋內(nèi)加入 2.5 μL B Buffer,同時用DNase/RNase-Free Water設(shè)置陰性對照(NC),蓋上管蓋后上下顛倒混勻7~8次,瞬間離心后,在水浴鍋中42 ℃孵育30 min。將獲得的RT-RAA擴增產(chǎn)物經(jīng)過DNA純化試劑盒純化,并經(jīng)2%瓊脂糖凝膠電泳后送生工生物工程(上海)股份有限公司進行測序驗證。
1.3.3 CRISPR-Cas13a檢測體系的建立與優(yōu)化 首先配制CRISPR-Cas13a反應(yīng)體系:包括crRNA(100 μmol·L-1)1.6 μL、TaqMan探針(50 μmol·L-1)0.64 μL、T7 RNA Polymerase 0.25 μL、NTP Buffer Mix 1.6 μL、MRI(40 U·μL-1) 1 μL、10×T7 RNA Polymerase Buffer 2.0 μL,向反應(yīng)混合物內(nèi)分別加入 0.5、1、2、4 μL LwaCas13a蛋白(60 μg·mL-1),加入0.6 μL cDNA(100 ng),用DNase/RNase-Free Water補至20 μL,37 ℃反應(yīng)20 min,以CRISPR-Cas13a檢測最短時間出現(xiàn)最強熒光信號為準則,確定最佳 LwaCas13a基因編輯蛋白用量。以最佳LwaCas13a蛋白用量配制反應(yīng)體系,在反應(yīng)體系中分別加入 0.2、0.4、0.6、0.8、1.6、3.2 μL crRNA(100 μmol·L-1),確定最佳crRNA濃度。運用順序控制變量的方法對CRISPR-Cas13a檢測體系進行優(yōu)化,其他物質(zhì)濃度設(shè)置如下:TaqMan探針(50 μmol·L-1)0.16、0.32、0.64、1.28 μL;T7 RNA Polymerase 0.25、0.5、1、2 μL;NTP Buffer Mix 0.8、1.2、1.6、2、2.4 μL,最終確定CRISPR-Cas13a的最佳反應(yīng)體系。
1.3.4 RT-RAA-CRISPR特異性試驗 對AIV H3、H5、H6、H7、H9、H10亞型以及感染家禽的其他常見病毒NDV、IBV、IBDV、DTMUV進行核酸抽提,將抽提的核酸進行RT-RAA擴增并進行純化。使用優(yōu)化后的CRISPR-Cas13a反應(yīng)體系進行反應(yīng),將DNase/RNase-Free Water作為陰性對照,于37 ℃反應(yīng)20 min,實時檢測熒光信號,驗證其特異性。
1.3.5 RT-RAA-CRISPR靈敏度試驗 使用pUC57-H5質(zhì)粒作為標準品,用DNase/RNase-Free Water進行10倍梯度稀釋,設(shè)置不同拷貝數(shù)(105~1 copies·μL-1)的質(zhì)粒作為模板進行RT-RAA擴增和CRISPR-Cas13a反應(yīng),并用DNase/RNase-Free Water設(shè)置陰性對照,檢測該反應(yīng)的檢出限。
1.3.6 RT-RAA-CRISPR符合性試驗 選取已知的20份陽性和36份陰性臨床樣本,對其進行核酸提取后,運用國標GB/T 18936—2020公布的引物(RT-qPCR-H5F和RT-qPCR-H5R)和探針(RT-qPCR-H5Taq)以及依據(jù)標準中提供的反應(yīng)體系和反應(yīng)程序進行熒光定量RT-PCR檢測,引物和探針的序列具體見表1。反應(yīng)體系為:2×RT緩沖液10 μL,酶混合液0.5 μL,上游引物(10 μmol·L-1)0.8 μL,下游引物(10 μmol·L-1)0.8 μL,探針(10 μmol·L-1)0.4 μL,無核酶滅菌水2.5 μL。反應(yīng)程序為:45 ℃ 反轉(zhuǎn)錄15 min;95 ℃ 預(yù)變性2 min;退火延伸95 ℃ 15 s,60 ℃ 60 s,40個循環(huán),在每次循環(huán)退火延伸時收集熒光信號,按照標準中提供的依據(jù)進行結(jié)果判定。同時將56份樣品用本研究建立的方法進行 RT-RAA-CRISPR 檢測,比較兩者之間的符合性。
為了獲得適用于檢測AIV H5亞型的高效特異性crRNA目的基因的靶向序列,經(jīng)過對GISAID數(shù)據(jù)庫中342條AIV H5亞型的HA基因序列進行篩選,根據(jù)crRNA設(shè)計原則找到合適的保守區(qū)域為第1 409位堿基到第1 440位堿基(圖1),該區(qū)域保守性一致率為98.1%:第1 416位堿基G占比97.6%,A占比2.4%;第1 429位堿基G占比96.7%,A占比3.3%;其他堿基一致率均>98.8%。因此,在該保守區(qū)域內(nèi)選擇了28 bp作為crRNA目的基因的靶向序列(5′-GCTTAGGGATAATGCAAAGGAGCTGGGT-3′),用于 AIV H5 亞型 CRISPR-Cas13a 核酸檢測靶點。該靶向序列的反向互補序列與36 bp的LwaCas13a 蛋白結(jié)合的錨定序列(5′-GAUUUAGACUACCCCAAAAACGAAGGGGACUAAAAC-3′)共同構(gòu)成crRNA,詳見表1。
圖1 crRNA保守性分析Fig.1 crRNA conserved analysis
按照RAA引物設(shè)計的原則通過Oligo 7.0軟件設(shè)計了RT-RAA引物,引物擴增產(chǎn)物長度為198 bp,擴增序列包含crRNA的靶向序列且與之不重合,應(yīng)用Primer-BLAST進一步檢驗引物的特異性。經(jīng)過篩選后,在RAA上游引物的5′端添加了25 bp的T7啟動子序列(5′-TAATACGACTCACTATAGGG-3′),以便后續(xù)對 RT-RAA 擴增產(chǎn)物進行DNA轉(zhuǎn)錄,序列如表1所示。
根據(jù)RT-RAA試劑盒說明書將抽提好的核酸擴增為cDNA并純化,凝膠電泳結(jié)果如圖2所示,所設(shè)計的針對H5病毒株的RT-RAA引物在198 bp的位置擴增出單一特異性條帶,與預(yù)期片段大小相符,經(jīng)測序驗證為正確序列。結(jié)果表明,本研究成功設(shè)計了可以擴增H5病毒株的RT-RAA引物,可利用擴增產(chǎn)物繼續(xù)進行后續(xù)試驗。
M. DL1000 DNA Marker;1. H5;2. 陰性對照M. DL1000 DNA Marker; 1. H5; 2. Negative control圖2 RT-RAA擴增結(jié)果Fig.2 RT-RAA amplification results
按照順序控制變量逐一優(yōu)化CRISPR-Cas13a檢測體系,以 CRISPR-Cas13a 檢測最短時間出現(xiàn)最強熒光信號為準則,確定最優(yōu)反應(yīng)體系。熒光檢測結(jié)果如圖3所示,表明CRISPR-Cas13a的最佳檢測體系為:LwaCas13a蛋白(60 μg·mL-1)2.0 μL、crRNA(100 μmol·L-1)3.2 μL、TaqMan探針(50 μmol·L-1)1.28 μL、T7 RNA Polymerase 0.25 μL、NTP Buffer Mix 2.0 μL、Murine RNase Inhibitor(40 U·μL-1)1 μL、10×Cas13a Reaction Buffer 2.0 μL、10×T7 RNA Polymerase Buffer 2.0 μL,加入0.6 μL cDNA(100 ng)后用DNase/RNase-Free Water補至20 μL,于 37 ℃ 反應(yīng)20 min,實時檢測 CRISPR-Cas13a 附帶切割釋放的熒光信號。
a. LwaCas13a蛋白用量的優(yōu)化;b. crRNA 濃度的優(yōu)化;c. TaqMan探針用量的優(yōu)化;d. T7 RNA Polymerase用量的優(yōu)化;e. NTP Buffer Mix用量的優(yōu)化a. Optimization of LwaCas13a protein dosage; b. Optimization of crRNA concentration; c. Optimization of TaqMan probe concentration; d. Optimization of the amount of T7 RNA Polymerase; e. Optimization of NTP Buffer Mix consumption圖3 CRISPR最佳檢測體系的優(yōu)化結(jié)果Fig.3 Optimized results of CRISPR optimal detection system
為了評價RT-RAA-CRISPR檢測方法的靈敏度,將pUC57-H5標準質(zhì)粒用DNase/RNase-Free Water進行10倍梯度稀釋,設(shè)置105、104、103、102、10、1 copies·μL-16個檢測濃度梯度,同時設(shè)置陰性對照(NC)。熒光檢測結(jié)果如圖4所示,反應(yīng)20 min后,從105copise·μL-1到1 copy·μL-1這6個檢測梯度的熒光值分別為14 553 189、11 586 768、11 375 226、12 802 810、6 765 644.5、841 933.25 a.u.,結(jié)果均與陰性對照的熒光值差異具有統(tǒng)計學(xué)意義(P<0.05)。結(jié)果表明,本研究建立的 RT-RAA-CRISPR 方法檢測AIV H5亞型的靈敏度為 1 copy·μL-1。
圖4 RT-RAA-Cas13a 靈敏度試驗熒光檢測結(jié)果Fig.4 Fluorescence detection results of RT-RAA-Cas13a sensitivity test
為評價 RT-RAA-CRISPR 檢測方法的特異性,分別用該方法檢測了AIV H3、H5、H6、H7、H9、H10亞型以及NDV、IBV、IBDV、DTMUV的RAA擴增產(chǎn)物。經(jīng)過CRISPR-Cas13a體系反應(yīng)后,各樣品熒光強度信號結(jié)果如圖5所示。結(jié)果表明,H5 病毒株的檢測熒光強度顯著高于 H3、H6、H7、H9、H10、NDV、IBV、IBDV、DTMUV和陰性對照,且與陰性對照有顯著性差異(P<0.01)。說明本研究建立的檢測AIV H5亞型方法特異性良好,與禽流感其他亞型及其他禽病均無交叉反應(yīng),可用于特異性檢測AIV H5亞型的病毒。
圖5 RT-RAA-Cas13a 特異性試驗熒光檢測結(jié)果Fig.5 Fluorescence results of RT-RAA-Cas13a specific assay
為進一步評估RT-RAA-CRISPR核酸檢測方法與國家標準規(guī)定熒光定量 RT-PCR方法的一致性,選取了20份AIV H5亞型的陽性臨床樣本和其他36份陰性臨床樣本,進行CRISPR核酸檢測和GB/T 18936—2020熒光定量RT-PCR檢測。結(jié)果如表2所示,其中19份熒光RT-PCR方法檢測陽性的樣本全部檢出,特異性為100%,敏感性為95%,符合率為98.2%。
表2 基于RT-RAA的CRISPR臨床樣品檢測結(jié)果Table 2 Results of CRISPR clinical samples based on RT-RAA
H5亞型HPAIV傳播快、致死率高,可引起家禽急性死亡,且具有重要的公共衛(wèi)生意義。從2003年1月至2022年10月,全球21個國家共報告了868例人感染H5N1禽流感病毒的病例,其中456例死亡,致死率高達53%[13]。因此,H5亞型HPAIV的快速、準確分型鑒別診斷對于養(yǎng)禽業(yè)和人類健康尤為重要。目前臨床對于AIV H5亞型的檢測方法主要包括病毒分離培養(yǎng)、抗原抗體檢測和分子生物學(xué)檢測。病毒分離培養(yǎng)和抗原抗體檢測雖然較為可靠,但同時也受到假陽性率高、檢測周期長、成本高、檢測過程較為復(fù)雜等缺點的影響。分子生物學(xué)研究可以同時檢測大量樣本,檢測效率高,適合于流感大面積暴發(fā)時使用。常用的分子檢測法有普通PCR、多重PCR(Multiplex PCR)、實時熒光定量PCR(Real-time PCR,RT-qPCR)、數(shù)字 PCR(digital PCR,dPCR)等。普通PCR和多重PCR雖然引物設(shè)計簡單、操作方便,但只能根據(jù)凝膠電泳進行定性試驗。實時熒光定量PCR和數(shù)字PCR雖然可以進行定量試驗,由于儀器成本高,使用復(fù)雜,且需要相關(guān)專業(yè)人員進行操作,在資源有限的地區(qū)無法大范圍推廣。因此,需要找到更合適的方法用于H5亞型禽流感病毒的即時檢驗(point-of-care testing, POCT),即在采樣現(xiàn)場即刻進行分析,省去樣本在實驗室檢驗時的復(fù)雜處理程序,快速得到檢驗結(jié)果。
核酸的快速檢測是臨床診斷和生物技術(shù)應(yīng)用中不可或缺的一部分。近幾年來,利用CRISPR-Cas酶的附帶切割活性,CRISPR-Cas13實現(xiàn)了對核酸靶標片段的檢測且被迅速開發(fā)、應(yīng)用于核酸診斷。但其單獨使用時靈敏度較低[8],因此,CRISPR-Cas13一般與其他核酸擴增技術(shù)一同使用。有研究利用PCR技術(shù)穩(wěn)定性好、實用性強及高效擴增靶序列的特點,建立了基于PCR和CRISPR-Cas13a系統(tǒng)的結(jié)核分枝桿菌脫氧核糖核酸(MTB DNA)檢測方法——PCR-CRISPR,其最高靈敏度為10 copies·μL-1[14],但PCR擴增過程所需時間較長,在檢測時間上幾乎沒有優(yōu)勢。也有研究利用CRISPR/Cas系統(tǒng)結(jié)合環(huán)介導(dǎo)等溫擴增技術(shù)(LAMP)實現(xiàn)了對非洲豬瘟病毒CD2v基因的快速現(xiàn)場化檢測,其檢測靈敏度達到8.88 copies·μL-1[15],但LAMP對引物的要求極高,需要設(shè)計4~6對引物,不利于該方法的推廣?;谥亟M酶介導(dǎo)等溫核酸擴增技術(shù)(RAA)與CRISPR-Cas13a技術(shù)相結(jié)合建立的RT-RAA-CRISPR具有引物設(shè)計簡單、反應(yīng)溫度低、操作簡易、檢測時間短等特點。本研究中對樣本RNA進行提取后,通過RT-RAA進行30 min的核酸反轉(zhuǎn)錄以及擴增過程,對獲得的cDNA進行純化,再通過CRISPR-Cas13a進行熒光檢測,在37~42 ℃條件下于1 h 20 min內(nèi)即可完成對RNA病毒樣本的檢測,可通過熒光值直觀判斷出病毒含量,短時間內(nèi)實現(xiàn)將待測靶標信號放大的同時對樣品進行定量檢測,大幅度提高了核酸檢測的效率。除具有檢測的快速性以外,該方法還具有良好的特異性和靈敏度,其最低檢測限度為1 copy·μL-1,比PCR-CRISPR和LAMP-CRISPR的檢測靈敏度高[14-15]。與其他檢測平臺(如熒光定量PCR)相比,RT-RAA-CRISPR方法的另一個優(yōu)勢是其使用成本低,該方法中熒光信號也可以通過熒光恒溫擴增儀、酶標儀進行檢測,相較于熒光定量PCR儀價格更加低廉,若采用便攜式熒光檢測儀替代,對檢測環(huán)境的要求更低、適應(yīng)性更好。
本研究成功建立了一種針對AIV H5亞型的基于RT-RAA的CRISPR快速檢測方法,可滿足生產(chǎn)實踐和臨床檢測的需要,為新型H5亞型禽流感病毒檢測方法的開發(fā)提供了技術(shù)支撐,在H5亞型禽流感病毒的防控、毒株的鑒別診斷以及流行病學(xué)調(diào)查方面均具有較好的應(yīng)用價值。