檳榔隱癥病毒1型TaqMan實時熒光定量PCR檢測方法的建立

林兆威 牛曉慶 唐慶華 宋薇薇 孟秀利 覃偉權

摘 ?要：為了建立一種快速、精準且能定量分析檳榔隱癥病毒1型（Areca palm velarivirus 1， APV1）的檢測方法。本研究參照GenBank中已登錄的APV1全基因組序列，在p21蛋白基因保守區域和p60蛋白基因保守區域之間設計1對特異性RT-PCR引物，并構建陽性重組質粒;在p60蛋白基因保守區域部分設計熒光定量PCR特異性引物和TaqMan-MGB探針，利用陽性重組質粒，構建標準曲線，并對該檢測方法的敏感性、特異性、穩定性及其應用效果進行驗證。結果顯示：該方法對陽性質粒標準品檢測的敏感性達3.0×101 copies/μL級別，是常規RT-PCR敏感性的100倍;標準曲線顯示，Ct值與拷貝數的對數呈線性關系，其標準曲線方程為y=-3.2748x+42.957，擴增效率為102%，相關系數R2為0.9992;該方法對APV1的檢測具有較好的特異性，其他檳榔病原物及內生菌不會對本方法產生非特異性干擾;批內與批間重復性試驗結果顯示，該方法對APV1的檢測具有較好的重復性和穩定性;利用該方法對海南萬寧、瓊海、定安、文昌等4個市（縣）采集的181份疑似樣品進行檢測，陽性率分別為91.95%、86.95%、89.65%、73.68%，其陽性的病毒拷貝數最高達1.59×107 copies/μL，表明海南部分地區APV1病毒載量較高、流行情況比較嚴重。

關鍵詞：檳榔;檳榔隱癥病毒1型;TaqMan實時熒光定量PCR

中圖分類號：S763.7 ? ? ?文獻標識碼：A

Abstract： It is aimed to develop a rapid， precision and quantitative detection method for Areca palm velarivirus 1 （APV1）. A pair of specificity RT-PCR primers were designed based on the APV1 complete genome sequence registered in GenBank to clone the sequence between the conservative region of p21 and the conservative region of p60， and a positive recombinant plasmid was constructed. Fluorescence quantitative PCR specificity primers and TaqMan-MGB probes were designed in the conservative region of p60 and a positive recombinant plasmid was constructed to establish a standard curve to verify the sensitivity， specificity， stability and application effect of the detection method. The sensitivity of this method to the detection of positive plasmid standards was up to 3.0×101 copies /μL， which was 100 times the sensitivity of conventional RT-PCR. The standard curve showed that the Ct value had a linear relationship with the logarithm of the copy number. The standard curve equation was y=-3.2748x+42.957， the amplification efficiency was 102%， and the correlation coefficient R2 was 0.9992. The method has good specificity for the detection of APV1， and other areca pathogens and endophytes will not be affected by this method. The intra batch repeat test and inter batch repeat test results showed that this method had good repeatability and stability in the detection of APV1. This method was used to detect 181 suspected samples collected from 4 cities and counties in Hainan， including Wanning， Qionghai， Ding’an， and Wenchang. The positive rate was 91.95%， 86.95%， 89.65%， and 73.68%， respectively， and the number of positive virus copies was up to 1.59×107 copies/μL. It shows that APV1 virus load is relatively high and the epidemic situation is relatively serious in some areas of Hainan.

Keywords： areca; Areca palm velarivirus 1; TaqMan RT-qPCR

DOI： 10.3969/j.issn.1000-2561.2021.11.006

檳榔（Areca cathecu L.）是棕櫚科檳榔屬多年生熱帶喬木，原產于馬來西亞，是我國的“四大南藥”（檳榔、砂仁、益智、巴戟天）之首，現主要分布在印度、中國、緬甸、孟加拉國及印度尼西亞等熱帶及亞熱帶國家與地區;在我國，檳榔主要栽培于海南、臺灣及云南等地，其中海南省檳榔種植面積占全國的95%以上，產量位居世界第二位，是海南省第二大熱帶經濟作物、海南省政府重點發展的“三棵樹”（橡膠、檳榔、椰子）之一[1-2]。近年來，國內檳榔消費市場迅速擴張，效益日益增長，據《海南省統計年鑒（2020）》統計，海南檳榔種植面積逐年擴大，2019年種植面積達11.52萬hm2。檳榔已成為海南省中東部地區230多萬農民的主要收入來源。

隨著檳榔種植面積的不斷擴大，植保問題日趨嚴重，檳榔黃化現象倍受關注。Wang等[3]研究指出，檳榔隱癥病毒1型（Areca palm velarivirus 1， APV1）是造成檳榔黃化現象的病害因子之一，且該病毒引發的癥狀與檳榔黃化病（arecanut yellow leaf disease， AYLD）的癥狀極其相似。隨后，本研究團隊對海南省檳榔黃化現象進行調查，并對APV1進行RT-PCR檢測分析，發現APV1在海南檳榔種植地區危害范圍廣，且為害面積較大，對檳榔產業健康發展構成潛在的威脅，引起本研究團隊的高度重視（該結果待發表）。然而，快速、精準的檢測技術對APV1的診斷尤為重要，當前對于該病毒的檢測主要采用RT-PCR技術[4]。由于TaqMan實時熒光定量PCR檢測技術具有較高的特異性、靈敏性及穩定性，并可對檢測的病原物進行定量分析，目前該方法廣泛應用于動植物的病原檢測[5-9]。因此，本研究根據APV1基因組序列，在p60蛋白基因保守序列中設計TaqMan- MGB探針，建立TaqMan RT-qPCR檢測方法，并對該方法的敏感性、特異性、穩定性及其應用效果進行了驗證，以期為APV1精準診斷、病原監測及分子流行病學的研究提供可靠的技術手段。

1 ?材料與方法

1.1 ?材料

1.1.1 ?供試樣品 ?檳榔隱癥病毒1型（Areca palm velarivirus 1， APV1）、檳榔壞死環斑病毒（Areca palm necrotic ringspot virus， ANRSV）及檳榔黃化植原體（arecanut yellow leaf phytoplasma， AYL）由中國熱帶農業科學院椰子研究所提供，并分別采用RT-PCR、巢氏PCR技術檢測為陽性;檳榔炭疽菌、檳榔葉片分離的內生或致病細菌：Burkholderia andropogonis、Pantoea ananatis、Sphingomonas yantingensis、Curtobacterium al-bidum、Bacillus cereus、Frondihabitans australicus由中國熱帶農業科學院椰子研究所提供;健康檳榔葉片樣品采自中國熱帶農業科學院椰子研究所檳榔苗圃，經RT-PCR對APV1檢測為陰性;疑似APV1樣品采自海南萬寧87份、瓊海46份、定安29份及文昌19份等4個市（縣）檳榔園，共181份。

1.1.2 ?主要試劑與儀器 ?主要試劑：RNAprep Pure多酚多糖植物總RNA提取試劑盒、反轉錄cDNA第一鏈合成試劑盒、2×Taq PCR mix、瓊脂糖凝膠回收試劑盒、質粒小提試劑盒及實時熒光定量PCR（探針法）試劑盒購于天根生化科技（北京）有限公司，pMDTM18-T Vector購自TaKaRa公司;DH5α感受態細胞購自北京全式金生物技術有限公司。主要儀器：Bio-Rad T100 PCR儀器，QuanStudio 6 Flex Real-time定量PCR儀，QuantStudio 3 Real-time定量PCR儀，超微量分光光度計 Thermo ND2000。

1.2 ?方法

1.2.1 ?引物與探針設計與合成 ?根據APV1基因組序列（GenBank登錄號：MK956940.2），利用DNAMAN軟件，在p21蛋白基因保守區域和p60蛋白基因保守區域設計上下游引物APV1-F/ APV1-R，該擴增片段用于構建陽性質粒標準品;并利用Primer Express 3.0軟件在p60蛋白基因保守的區域部分設計實時熒光定量PCR特異性引物APV1-F1/APV1-R1和探針APV1-Probe（表1）。引物與探針由生工生物工程（上海）股份有限公司合成。

1.2.2 ?RNA提取及反轉錄cDNA合成 ?利用RNAprep Pure多酚多糖植物總RNA提取試劑盒，參照說明書，提取樣品總RNA，并利用反轉錄cDNA第一鏈合成試劑盒，將提取的總RNA進行反轉錄，獲得cDNA，?20 ℃保存。

1.2.3 ?實時熒光定量PCR引物驗證 ?使用引物APV1-F1/APV1-R1和探針APV1-Probe，以3份陽性樣品cDNA為模板，健康檳榔樣品cDNA陰性對照，ddH2O為空白對照進行TaqMan RT- qPCR，參照試劑盒說明書，20 μL反應體系：2×SuperReal PreMix （Probe） 10μL，正反向引物（10 μmol/L）各0.6 μL，熒光探針（10 μmol/L）0.4 μL，cDNA模板2 μL，50×ROX Reference Dye3 0.2 μL，RNase-Free ddH2O 6.5 μL。使用QuantStudio 3 Real-time定量PCR儀，設置TaqMan RT-qPCR反應程序：95 ℃預變性15 min;95 ℃變性3 s，60 ℃退火/延伸30 s，40個循環。

1.2.4 ?陽性質粒標準品的制備 ?使用引物APV1-F/APV1-R，以APV1 cDNA為模板，健康檳榔樣品cDNA陰性對照，ddH2O為空白對照進行RT-PCR，PCR產物進行瓊脂糖凝膠電泳檢測，切下目的條帶，利用瓊脂糖凝膠回收試劑盒進行目的條帶純化，純化后產物連接到pMDTM18-T Vector，并轉化進DH5α感受態細胞中，篩選并提取陽性克隆子的質粒，送生工生物工程（上海）股份有限公司進行測序。將測序正確的質粒使用超微量分光光度計測定濃度，置于?20 ℃保存備用。運用Shirima等[10]的方法計算重組質粒中的拷貝數，重組質粒拷貝數（copies/μL）=（阿伏伽德羅常數×C）/（109×重組質粒堿基對數×660 dalton/bp）。其中阿伏伽德羅常數為6.02×1023;C為重組質粒的濃度（ng/μL）。最終計算得到重組質粒的拷貝數為3.0×1010 copies/μL。RT-PCR反應體系：2×Taq PCR mix 12.5 μL，上下游引物（10 μmol/L）各1.0 μL，cDNA模板2.0 μL，ddH2O 8.5 μL。目的片段RT-PCR反應程序：95 ℃預變性5 min;95 ℃變性45 s，53 ℃退火45 s，72 ℃延伸40 s，35個循環;72 ℃延伸10 min。

1.2.5 ?敏感性試驗與標準曲線的建立 ?將1.2.4獲得的陽性質粒標準品為初始模板，按10倍梯度稀釋出10個標準品濃度，拷貝數依次為3.0×109、3.0×108、3.0×107、3.0×106、3.0×105、3.0×104、3.0×103、3.0×102、3.0×101、3.0×100 copies/μL。以上10個梯度濃度陽性質粒標準品為模板，健康檳榔樣品為陰性對照，進行TaqMan RT-qPCR。使用QuanStudio 6 Flex Real-time定量PCR儀，反應體系和程序參照1.2.3，得到動力學擴增曲線和PCR擴增循環閾值（Ct值）。運用Excel 2020軟件，以Ct值為縱坐標，10為底質粒拷貝數的對數值為橫坐標作標準曲線，獲得回歸直線方程，并計算該組體系的擴增效率：擴增效率（%）= 10（?1/斜率）?1。并將上述稀釋的10個濃度梯度的陽性質粒標準品為模板，健康檳榔樣品為陰性對照，使用APV1-F1/APV1-R1引物進行RT-PCR，PCR產物經凝膠電泳后驗證，與TaqMan RT-qPCR敏感性的結果比較。

1.2.6 ?特異性試驗 ?以檳榔APV1的cDNA為陽性對照、健康檳榔cDNA為陰性對照，ANRSV、檳榔黃化植原體、Burkholderia andropogonis、Pantoea ananatis、Sphingomonas yantingensis、Curtobacterium albidum、Bacillus cereus、Frondihabitans australicus的cDNA為模板進行TaqMan RT-qPCR，檢測其特異性，使用QuantStudio 3 Real-time定量PCR儀，反應體系和設置程序參照1.2.3。

1.2.7 ?重復性試驗 ?選取3.0×109、3.0×108、3.0×107、3.0×106 copies/μL 4個濃度的陽性質粒標準品為模板，進行TaqMan RT-qPCR檢測，每個濃度重復3次，共檢測3次，進行批內重復性試驗，計算Ct值的變異系數。同時，將上述濃度的陽性質粒標準品保存在?20 ℃ 3、6、9 d后，以相同反應條件重復進行批間試驗，計算Ct值的變異系數，驗證該方法的重復性和穩定性。使用QuanStudio 6 Flex Real-time定量PCR儀，反應體系和設置程序參照1.2.3。

1.2.8 ?田間樣本檢測 ?從海南萬寧、瓊海、定安、文昌等4個市（縣）檳榔園中分別采集87、46、29、19份疑似APV1病樣，并將這181份疑似樣品進行總RNA提取，反轉錄成cDNA，運用本研究建立的TaqMan RT-qPCR方法進行檢測，健康檳榔樣品為陰性對照，獲得該4個市（縣）檳榔APV1的檢出率，并根據1.2.5的標準曲線回歸方程，計算田間樣品病毒的最高、最低拷貝數。使用QuanStudio 6 Flex Real-time定量PCR儀，反應體系和設置程序參照1.2.3。

2 ?結果與分析

2.1 ?熒光定量PCR引物的驗證

使用引物APV1-F1/APV1-R1和探針APV1- Probe對3份APV1陽性樣品和1份陰性樣品進行TaqMan RT-qPCR。如圖1所示，3份陽性樣品均出現“S”形擴增曲線，而陰性樣品和空白對照未出現擴增曲線，說明該TaqMan RT-qPCR的引物/探針組合能檢測到APV1，可進一步實驗。

2.2 ?陽性標準品的制備

以3份APV1陽性樣品為模板，使用引物APV1-F/APV1-R進行RT-PCR，經凝膠電泳，結果如圖2，獲得約490 bp的擴增片段。經純化后的產物進行載體連接、感受態細胞轉化，并將重組質粒測序后進行NCBI blast比對，與報道的APV1序列（GenBank登錄號：MK956940.2）一致性達100%，說明已獲得陽性標準品。

2.3 ?敏感性試驗與標準曲線的建立

將3.0×109～3.0×100 copies/μL 10個濃度梯度的陽性質粒標準品為模板，健康檳榔樣品為陰性對照，進行TaqMan RT-qPCR和RT-PCR，并利用Ct值與拷貝數的對數構建標準曲線。結果表明，該引物/探針對10個濃度梯度的陽性質粒標準品最低能檢測到3.0×101 copies/μL（圖3），而RT-PCR最低能檢測到3.0×103 copies/μL（圖4），TaqMan RT-qPCR檢測APV1的敏感性是常規RT-PCR的100倍。標準曲線顯示，Ct值與拷貝數的對數呈線性關系，其標準曲線方程為y=?3.2748x+42.957，擴增效率為102%，相關系數R2為0.9992（圖5）。

2.4 ?特異性試驗

使用檳榔APV1及葉片中主要的病原菌和內生細菌進行TaqMan RT-qPCR特異性驗證，結果如圖6所示，僅APV1出現擴增曲線，健康檳榔陰性對照、檳榔植原體、ANRSV、檳榔炭疽菌、Burkholderia andropogonis、Pantoea ananatis、Sph-i?ngomonas yantingensis、Curtobacterium albid?um、Bacillus cereus、Frondihabitans australicus均未擴增出，說明該TaqMan RT-qPCR方法對APV1的檢測具有較好的特異性，其他檳榔病原物和內生菌不會對本檢測方法產生非特異性干擾。

2.5 ?重復性試驗

選取3.0×109、3.0×108、3.0×107及3.0×106 copies/L 4個濃度的陽性標準品為模板，進行TaqMan RT-qPCR檢測，分別驗證該方法的批內和批間的重復性和穩定性。結果如表2所示，批內變異系數小于1%，批間變異系數小于2%，表明該方法對APV1的檢測具有較好的重復性和穩定性。

2.6 ?田間樣本檢測

運用建立的TaqMan qRT-PCR方法對海南萬寧、瓊海、定安、文昌等4個市（縣）的檳榔園中采集的181份疑似樣品進行檢測。結果顯示，萬寧、瓊海、定安、文昌的陽性率分別為91.95%、86.95%、89.65%、73.68%，其陽性的病毒拷貝數在11.23 copies/μL～1.59×107 copies/μL之間，表明海南部分地區APV1病毒載量較高、流行情況比較嚴重。

3 ?討論

3.1 ?TaqMan RT-qPCR檢測APV1方法的建立

目前對APV1的檢測主要有RT-PCR技術和酶聯免疫檢測技術[4， 11]。相較于TaqMan RT-qPCR技術，傳統的RT-PCR技術敏感性較低，且耗時久;雖然酶聯免疫檢測具有便捷、快速及直觀的特點，但該技術敏感性未能達到更高的檢測要求，并且準確性往往依賴于抗血清的品質[12]。TaqMan RT-qPCR技術具有較高的特異性和敏感性，該技術不僅可用于定性檢測，也可用于定量分析。在本研究中，該TaqMan RT-qPCR檢測方法的敏感性可達3.0×101 copies/μL級別，比RT-PCR的敏感性高100倍，在實際應用中，可檢測含量最低至11.23 copies/μL的病毒，并且具有較高的穩定性和特異性。因此，該TaqMan RT-qPCR檢測方法的建立，對APV1的病原監測及分子流行病學的研究提供了可靠的技術手段。

3.2 ?APV1對檳榔產業發展具有潛在的威脅性

隱癥病毒屬可為害甜菜、葡萄、柑橘、小麥及胡蘿卜等多種經濟作物，屬于正義ssRNA病毒[13]。由隱癥病毒屬引起的甜菜黃化病毒（Beet yellows virus， BYV）一般年份的發病率達50%～60%，塊根減產20%～25%，含糖率下降2%～3%，染病種子減產30%[14];葡萄卷葉伴隨病毒2（Grapevine leafroll-associated virus 2， GLRaV-2）是我國葡萄卷葉病的主要病原之一，該病毒嚴重影響葡萄產量和質量[15];柑橘衰退病毒（Citrus tristeza virus， CTV）則嚴重為害世界柑橘生產的重要病原物之一[16]。檳榔隱癥病毒1型（APV1）為長線形病毒科、隱癥病毒屬病毒，2015年Yu等[17]首次報道在國內發生流行，并對APV1-HN株進行全基因組測序，此后關于APV1的報道較少。至2020年Wang等[3]對APV1-WNY株進行全基因組測序，并闡述APV1與檳榔黃化病存在相關性和APV1的潛在傳播媒介為粉蚧。隨后，本研究團隊于2020年對海南省檳榔黃化現象進行調查，并檢測APV1，發現APV1在全省檳榔中分布廣，并且發病率較高。在本研究中，應用TaqMan RT-qPCR檢測方法對海南萬寧、瓊海、定安、文昌等4個市（縣）的檳榔園中采集的181份疑似樣品進行檢測，其陽性率分別為91.95%、86.95%、89.65%、73.68%，表明APV1病毒在海南檳榔主栽區發病率高，對海南檳榔產業發展具有潛在的威脅性，需提高重視，加強該病害的監測、檢測及其流行病學的研究。

參考文獻

[1] 楊連珍， 劉小香， 李增平. 世界檳榔生產現狀及生產技術研究[J]. 世界農業， 2018（7）： 121-128.

[2] 代文婷， 吉建邦， 康效寧， 等. 檳榔干制技術研究進展[J]. 食品工業， 2021， 42（3）： 258-262.

[3] Wang H X， Zhao R B， Zhang H W， et al. Prevalence of Yel-low Leaf Disease （YLD） and its associated Areca palm velarivirus 1 （APV1） in betel palm （Areca catechu） plantations in Hainan， China[J]. Plant disease， 2020， 104（10）： 1-33.

[4] 黃 ?惜， 王洪星， 趙瑞白， 等. 一種檳榔黃化病相關病毒及其檢測方法： CN111187756A[P]. 2020-05-22.

[5] 韓 ?玉， 王 ?濤， 潘 ?力， 等. 基于非洲豬瘟病毒CD2v基因TaqMan熒光定量PCR檢測方法的建立[J]. 中國獸醫學報， 2021， 41（5）： 847-852.

[6] 王 ?琛， 馬躍宇， 黃思超， 等. 蜜蜂殘翼病毒TaqMan- MGB探針熒光定量RT-PCR檢測方法的建立[J]. 病毒學報， 2021， 37（3）： 686-694.

[7] 王 ?濤， 韓 ?玉， 潘 ?力， 等. 針對非洲豬瘟病毒MGF360-13L基因的TaqMan熒光定量PCR的建立[J]. 中國農業科學， 2021， 54（5）： 1073-1080.

[8] 王柱華， 王文鵬， 袁恩平， 等. 喜樹叢枝植原體的分子鑒定及TaqMan探針實時熒光定量PCR檢測方法的建立[J].植物病理學報， 2021， 51（3）： 429-440.

[9] 姜自紅， 殷培峰. 3種菊花病毒/類病毒實時熒光定量RT-PCR檢測體系的構建與準確性評價[J]. 病毒學報， 2021， 37（1）： 169-180.

[10] Shirima R R， Maeda D G， Kanju E， et al. Absolute quantifi-cation of Cassava brown streak virus mRNA by real-time qPCR[J]. Journal of Virological Methods， 2017， 245： 5-13.

[11] 陳 ?陽， 王洪星， 黃 ?惜. 檳榔APV1病毒多克隆抗體制備及酶聯免疫檢測[J/OL]. 分子植物育種. [2021-06-23]. http：//kns.cnki.net/kcms/detail/46.1068.S.20210521.1404.008.html.

[12] 張 ?偉， 徐 ?碩， 陳德鑫. 常用植物病毒病檢測技術比較[J]. 南方農業， 2013， 7（12）： 24-28.

[13] 洪雪玉， 陳劍平. 長線形病毒科成員編碼蛋白的功能[J]. 浙江農業學報， 2004（1）： 51-56.

[14] 劉長兵. 甜菜黃化病毒發生與防治[J]. 新疆農業科技， 2009（4）： 41.

[15] Martin R R， Eastwell K C， Wagner A， et al. Survey for vi-ruses of grapevine in Oregon and Washington[J]. Plant Dis-ease， 2005， 89（7）： 763-766.

[16] 宋 ?震， 周常勇. 柑橘衰退病毒分子生物學研究進展[J]. 西南農業學報， 2005（2）： 213-218.

[17] Yu H M， Qi S S， Chang Z X， et al. Complete genome se-quence of a novel velarivirus infecting areca palm in Chi-na[J]. Archives of virology， 2015， 160（9）： 2367-2370.

責任編輯：謝龍蓮