葡萄雙生病毒A的序列分析及檢測

中圖分類號：S665.1;S436.65 文獻標志碼：A文章編號：1001-4330（2025）04-0887-07

摘要：【目的】研究新疆吐魯番市葡萄雙生病毒A（Grapevine geminivirus A）的田間侵染狀況，檢測采集71份疑似葡萄病毒樣本，為明確葡萄雙生病毒A的發生、分布及其分子檢測提供理論指導?！痉椒ā恳罁咄繙y序結果，驗證樣本中葡萄雙生病毒A，測定GGVA的全基因序列，以葡萄雙生病毒A的全基因序列為靶標設計特異性引物進行疑似葡萄病毒樣本的檢測?！窘Y果】測定GGVA的全基因組序列，獲得GGVA-QKT分離株，其核苷酸同源性為 94.3%～99.3% ，與已公布的分離株具有高度同源性，其中與來自以色列（KX618694.1）、（NC_031340.1）的2個GGVA分離株親緣關系最近。【結論】71份不同葡萄樣本均有感染雙生病毒A，其中，新疆鄯善縣檢出率最高，檢出率為 30.55% ，新疆托克遜縣檢出率最低，檢出率為 9.10% 。該病毒在新疆吐魯番市葡萄種植區普遍發生。

0 引言

【研究意義】葡萄（Vitisspp.）是世界四大水果之一，由于多年的無性繁殖、苗木流轉等多導致葡萄中存在多種病原體，葡萄病毒普遍發生[1]。目前，世界已報道超過86種病毒[2，我國已報道22 種葡萄病毒[3]。在葡萄藤上發現多種DNA 病毒[4]，如葡萄清靜脈病毒（GVCV）[5]、葡萄葉變色相關病毒（GRLDaV）、葡萄紅斑病毒（ GRBaV ）[7，8.9]和葡萄雙生病毒A（GrapevinegeminivirusA，GGVA）[10]。而葡萄雙生病毒從最初被檢測到，陸續在多個國家和地區報道[11-14]在我國云南、山東、遼寧、四川、新疆等?。▍^）均有發現[1.5.6]。雙生病毒（GGVA）屬于雙生病毒科Maldovirus屬的成員，基因組為環狀單鏈DNA（ssDNA），被包裹在一對不完全二十面體中，可通過感染寄主細胞中的雙鏈DNA（dsDNA）中間體復制[7]。其全基因組序列長度為 2903～2907 bp，具有一個保守的莖環序列“TAATATTAC”，正義鏈上有2個開放閱讀框（ORF）（V1，外殼蛋白；V2，推定運動蛋白），負補鏈上有四個開放閱讀框（ORF），（C1，復制相關蛋白；C2，轉錄激活蛋白；C3，復制增強子;C4，宿主激活蛋白）[10]。葡萄雙生病毒A是近幾年報道的葡萄DNA病毒，研究發現，該病毒屬于潛隱性病毒，多數在無癥狀的葡萄葉片檢測出?！厩叭搜芯窟M展】范旭東等[\"]在引進的栽培葡萄品種‘夏季黑’上發現該病毒，且發生普遍;孫蘇偉等[15]在云南葡萄種植區發現雙生病毒（GGVA）并進行基因組序列分析;趙小春[16]在新疆葡萄種植區發現葡萄雙生病毒（GGVA）；田沂民等[18]對葡萄雙生病毒進行傳毒介體檢測鑒定相關研究?！颈狙芯壳腥朦c】關于吐魯番市葡萄雙生病毒A的檢測相關報道較少，而在新疆吐魯蕃市采集疑似病毒樣本中，通過混樣測序發現葡萄雙生病毒A，需明確吐魯番市葡萄雙生病毒A田間侵染情況?！緮M解決的關鍵問題】采集吐魯番市主要葡萄種植區疑似葡萄病毒樣品材料并低溫保存，對采集的71份疑似病毒病樣品進行PCR檢測，并進行序列鑒定及系統進化分析。提取葡萄植株中總DNA并進行PCR檢測，對擴增產物進行序列測定并在GeneBank中進行序列比對，分析吐魯番市葡萄雙生病毒A田間發生情況，為明確我國葡萄雙生病毒A的發生、分布及其分子檢測提供理論依據。

材料與方法

1.1 材料

1.1.1 供試葡萄

于2022年6月以吐魯番市葡萄種植區為調查區域，對3個葡萄種植區進行疑似病毒樣本的采集，共71份，將采集的樣品進行表面清洗并放在攜帶的液氮罐里，注明采集日期、品種名稱及編號，樣本帶回實驗室置于 -80% 冰箱保存、備用。表1，圖1

Notes：A、B、C：Suspected grape virus disease symptoms leaf；D： Healthyleaf

1. 1.2 主要儀器與試劑

儀器：超低溫冰箱（青島海爾公司）、超凈工作臺、PCR儀（美國伯樂公司）電泳儀、凝膠成像系統（上海嘉鵬公司）、恒溫水浴鍋、臺式智能精密搖床、高速臺式冷凍離心機等。

試劑：總DNA提取試劑盒DNAsecure新型植物基因組DNA提取試劑盒購自天根公司； 2×Es.

TaqMasterMix（Dye）購自康為世紀公司；Marker（DNAMarkerII、DL200O）購自北京全式金公司、pUCm-T克隆載體、DH5a感受態、凝膠回收試劑盒等購自上海生工公司。

1.1. 3 供試引物

依據高通量測序結果，通過NCBI對比分析后，參照相似性最高的GGVA全基因組序列，以其

為靶標，利用Primer5.0軟件設計（ F^2* 和 R^2* ）GG-VA-F（5’-ACCATGTCTGAGGGTAAGGC-3'和GGVA -R（5’ -GGAGTCAATTGGAACGTGCG -

3）特異性檢測引物和擴增全長基因組引物，并在NCBI數據庫中進行特異性引物比對，預計產物片段大小為 329bp ，引物均有由上海生工合成。表2

表2供試引物

1.2 方法

1. 2.1 總DNA提取

參考DNAsecure新型植物基因組DNA提取試劑盒說明書，提取71份樣本總DNA，進行電泳檢測其完整性，并保存于 -80^°C 冰箱備用。

1.2.2 高通量測序

將采集的疑似葡萄病毒樣本混合，送往公司在IlluminaHiseq高通量測序平臺檢測，確認RNA的純度（ 0D_260/280 ）、濃度、核酸吸收峰與RNA的完整性是否正常；樣品檢測合格后，通過試劑盒去除rRNA，進行文庫的構建；文庫構建后，檢測文庫質量，通過不同文庫按照目標下機數據量進行poo-ling，用IlluminaHiSeq平臺進行測序分析;測序最終得到原始數據，處理原始數據，處理后進行測序數據質量的控制，得到有效數據。最后，進行有效數據轉錄本拼接和組裝，拼接完成后得到Nr數據庫和CDS預測序列;將Nr數據庫進行Taxonomic分析，得到未注釋序列和病毒序列信息；基于病毒序列進行組裝分析，獲得完整的病毒種類信息，對于未注釋的序列信息，進行dsRNA預測分析，獲得RNA病毒和類病毒的序列信息[16，19，20]

1.2.3 PCR擴增

根據（ F^2* 和 R^2* ）的引物進行PCR擴增，PCR反應體系為 25μL 2×Es Taq MasterMix（Dye） 12.5μL ，正反引物（ 10μmol/L ）各 1μL cDNA模板 2μL ， RNase free ddH₂O 補平 25μL PCR擴增條件： 94°C 5min;94°C 3min，55°C 45s ，72°C 1min，35 個循環： 72%10min 。PCR擴增產物在 1.5% 瓊脂糖 1×TAE 緩沖系統中電泳，每孔加樣 5μL ，用凝膠成像系統觀察并攝影記錄。

1.2.4PCR產物測定及序列

將PCR擴增產物進行回收，具體步驟參考瓊脂糖凝膠DNA回收試劑盒說明書。經過連接、純化DNA，連接至 克隆載體，轉化大腸桿菌DH5a，獲得陽性克隆，測序由上海生工完成。測序結果利用DNAman軟件進行拼接；通過登錄NCBI利用BLAST檢索功能進行序列同源性分析，使用SnapGene軟件進行序列分析，采用MEGA11軟件構建系統發育樹（neighbor joining，NJ，Bootstrap重復次為1000）。

1. 2.5 GGVA全長基因組序列鑒定

通過NCBI對比分析后，參照相似性最高的韓國分離物jaok-vf-1全基因組序列（登錄號：MF163263）設計引物，分段擴增GGVA全長基因組序列。PCR反應體系為（ 25μL ）： 2×Es Taq MasterMix（Dye） 12.5μL ，正反引物（ 10μmol/L ））各 1μL ，DNA 模板 2μL ，RNase free ddH₂O 補平25μL ;PCR反應條件： 94°C 5min，94°C 30s，55°C 電 30s，72%1min，72%10min，35 個循環。 1.5% 瓊脂糖凝膠電泳，凝膠成像，經過連接、純化DNA，連接至 克隆載體，轉化大腸桿菌DH5a，獲得陽性克隆，測序由上海生工公司完成，為確保送測的準確性，至少送測3個陽性克隆樣品進行測序。

2 結果與分析

2. 1 高通量測序結果

研究表明，共得到5種葡萄病毒序列，其中在一個未知屬中有4個未知contigs，包含3588條Reads，將得到的contig通過NCBI進行序列比對；有1個contigs與已報道的葡萄雙生病毒A有較高的匹配度，其長度為 586～2058bp ，與已公布的雙生病毒（GGVA）的序列同源性有 94.87% ＼～99.73% 。表3

2.2 引物的擴增驗證

研究表明，所有目的條帶與空白對照均正常，特異性引物檢測條帶約為 329bp ，而擴增條帶與目的帶存在偏差，后期測序結果證實，特異性引物擴增的 500bp 條帶也為目的基因序列。與美國分離株（ KX570608.1，KX618694.1） 、印度分離株（OQ427641.1）相似性為 95.1%～99% ，所測序列與公布序列具有高度同源性，所檢測病毒及類病毒種類與公布的種類相同。圖2

注：A：病毒的部分PCR檢測；B：PCR的擴增；M：marker;CK：陰性對照Notes：A：Part of the virus PCR detection；B：PCR amplification；M：marker；CK：Negative control

2.3 GGVA分離物全基因組序列

研究表明，該病毒序列含有6個ORF，分別編碼6種蛋白，ORF1編碼1個長度為 312nt （313～624bp 的 Pre-CP/V2 ;ORF2編碼1個長度為771nt（446～1216bp） 的CP/V1;ORF3編碼1個長度為 429nt（1213～1641bp） 的 Ren/C3 ; ORF4編碼1個長度為 的TrAP/C2;ORF5編碼1個長度為 1221nt（1683 ～2903bp 的 Rep/Cl 。ORF6編碼1個長度為 的HAP/C4。GGVA-QKT分離株與已報道的分離株全基因組核苷酸同源性為 94.3%～99.3% ，所測序列與公布的種類相同。獲得新疆分離株與已公布分離株具有高度同源性，且沒有顯著差異，該病毒不存在地域分布規律。圖3

2.4 GGVA分離物的系統進化

研究表明，供試的全基因組序列分離物共分為3個大組，獲得的（中國新疆吐魯番市鄯善縣）七克臺分離株與日本（KX570616.1）、（KX570615.1）、（KX570612.1）、（KX570613.1）、（LC424098.1）、匈牙利（KX570618.1）、新西蘭（MK690474.1）、中國（KX950822.1）、（MT344710.1）、（MT344704.1）、（MT344706.1）、（MT344712.1）、以色列（KX618694.1）、（NC_031340.1）、韓國（KX570607.1）等的GGVA分離株聚為一支，且與來自以色列（KX618694.1）、（NC_031340.1）的2個GGVA分離株親緣關系最近；來自中國（MT344703.1）、（MT344713.1）、（KX570611.1）、日本（KX570610.1）、韓國（MK336184.1）、中國新疆（ON868755.1）的GG-VA分離株聚為一支；與來自中國（KX574323.1）、韓國（KX570609.1）、（KX570614.1）的GGVA分離物聚為一支。圖3

2.5 葡萄樣本的PCR檢測

研究表明，71份不同地區樣本均有感染葡萄雙生病毒（GGVA），其中，善縣檢出率最高，檢出率為 30.55% ;其次是吐魯番市高昌區，檢出率為 20.83% ；托克遜縣檢出率最低，為 9.10% ，該病毒在吐魯番葡萄種植區發生普遍。表4

3討論

3.1該病毒是在鮮食葡萄上通過深度測序技術篩查發現的葡萄新病毒[10]；Jo等[12]在多個葡萄品種中首次檢出該病毒;其后，范旭東等[\"在多地葡萄園進行調查時，發現普通栽培品種和引進栽培品種夏季黑攜帶GGVA，且發生普遍，而在本地栽培品種未發現GGVA。

3.2研究測定了GGVA的全基因組序列，并進行序列相似性分析，雙生病毒（GGVA-QKT）分離株核苷酸同源性為 94.3%～99.3% ，與已公布分離株具有高度同源性，其中與來自以色列（KX618694.1）（NC_031340.1）的2個GGVA分離物親緣關系最近。通過對采集的71份疑似病毒樣本進行PCR檢測。結果表明，71份不同地區樣品感染均有感染雙生病毒（GGVA），其中，善縣檢出率最高，檢出率為 30.55% ；其次是高昌區，檢出率為 20.83% ；托克遜縣檢出率最低，為9.10% ，從檢測結果表明，該病毒在吐魯番葡萄種植區發生非常普遍，可見葡萄雙生病毒（GGVA）在新疆已廣泛分布。

3.3葡萄雙生病毒A可侵染木本植物和草本植物，如已報道的柑橘褪綠矮縮相關病毒（CCDaV）[21]、麻風樹花葉病毒（JaMV）[22]、桑樹花葉矮縮相關病毒病毒（MMDaV）[23]等。雙生病毒（GGVA）為典型的蟲傳病毒，主要傳毒介體為粉虱、葉蟬等，田沂民等利用PCR或RT-PCR方法在粉虱、薊馬、狹領紋唇盲蝽、菊方翅網蝽、蚜蟲和寡脈蠅等6種昆蟲體內均檢出GGVA。而薊馬、狹領紋唇盲蝽、菊方翅網蝽、寡脈蠅未見相關研報道證明是否傳毒，幾種昆蟲是否傳毒以及傳毒特點需要進一步研究。

4結論

對吐魯番市采集的71份樣品進行PCR檢測，共檢測出17份陽性樣品，檢出率為 23.94% ;以鄯善縣、高昌區檢出率較高，托克遜縣檢出率較低；測定GGVA的全基因組序列，其核苷酸同源性為94.3%～99.3% ，與來自以色列（KX618694.1、NC_031340.1）的2個GGVA分離物親緣關系最近。

參考文獻（References）

[1]梁巧玲，乾義柯，張娜，等．新疆三種葡萄病毒RT-PCR 檢 測及序列分析[J].植物保護學報，2015，42（3）：376-381. LIANGQiaoling，QIAN Yike，ZHANG Na，et al.RT-PCR detection and sequence analysis on three grapevine viruses from XinjiangUyghur Autonomous Region[J].Journal of Plant Protection，2015，42（3）：376-381.

[2]Fuchs M. Grapevine viruses：a multitude of diverse species with simple but overall poorly adopted management solutions in the vineyard[J]. Journal of Plant Pathology，2020，102（3）：643- 653.

[3］張亞萍．新疆6個釀酒葡萄品種病毒及類病毒的分子檢測 [D]．石河子：石河子大學，2023. ZHANG Yaping. Molecular detection of viruses and viroids in 6 winegrape varieties in Xinjiang[D]. Shihezi；Shihezi University， 2023.

[4]GPMartelli. Directory of virus and virus -like diseases of the grapevine and their agents：introduction[J].Journal of Plant Pa 5 thology，2014， 96（Suppl.1）：1-4

[5]GuoQ，Honesty S，XuML，etal.Genetic diversity and tissue and host specificity of Grapevine vein clearing virus[J]. Phytopathology，2014， 104（5）：539-547 ：

[6]Maliogka VI， Olmos A，Pappi PG，et al.A novel grapevine badnavirus is associated with the Roditis leaf discoloration disease [J].VirusResearch，2015，203：47-55.

[7]AlRwahnihM，DaveA，AndersonMM，etal.Associationofa DNA virus with grapevines affected byred blotch disease in California[J].Phytopathology，2013，103（10）： 1069-1076.

[8]AlRwahnihM，RowhaniA，GolinoD A ，et al.Detection and genetic diversity of Grapevine red blotch-associated virus isolates in table grape accessions in the National Clonal Germplasm Repositoryin California[J]. Canadian Journalof PlantPathology， 2015，37（1）：130 -135.

[9]Blanco-UlateB，HopferH，Figueroa-BalderasR，et al．Red blotch disease alters grape berry development and metabolism by interfering with the transcriptional and hormonal regulation of ripening[J].Journal of Experimental Botany，2017，68（5）：1225 -1238.

[10]AlRwahnihM，AlabiOJ，Westrick NM，et al.Description of a novel monopartite geminivirus and its defective subviral genome in grapevine[J].Phytopathology，2017，107（2）：240-251.

[11]FanXD，ZhangZP，RenF，etal.First report of grapevine geminivirus A from grapevines in China[J].Plant Disease，

[12]JoY，ChoiH，Song MK，etal.First report of grapevine geminivirusA indiverse Vitis speciesin Korea[J].Plant Disease， 2018，102（1）：255.

[13]Veerakone S，Blouin AG，Barrero R A，et al. First report of grapevine geminivirus A in Vitis in New Zealand[J]. Plant Disease，2020，104（2）：600.

[14] Kishan G， Kumar R， Sharma S K ， et al. Development and application of crude sap-based recombinase polymerase amplification assay for the detection and occurrence of grapevine geminivirusAin Indian grapevine cultivars[J].Frontiers in Plant Science， 2023，14： 1151471.

[15]SunSW，HuY，JiangGZ，etal．Molecular characterization and genomic function of grapevine geminivirus A[J]. Frontiers in Microbiology，2020，11：555194.

[16］趙小春．新疆葡萄病毒病的病原鑒定及檢測[D]．石河子： 石河子大學，2022. ZHAO Xiaochun. Pathogen identification and detection of grape virus disease in Xinjiang[D]. Shihezi： Shihezi University，2022.

[17]KingA.Virus taxonomy ninth report of the International Committee on TaxonomyofViruses[M].Oxford：Elsevier，2011.

[18]田沂民，朱雅君，羅金燕，等．上海地區葡萄病毒種類及傳 毒介體檢測鑒定[J]．上海農業學報，2021，37（04）：52- 56. TIAN Yimin，ZHU Yajun，LUN Jinyan，et al. Detection and identification of grape virus species and vectors in Shanghai[J]. Journal of Shanghai Agricultural Sciences，2019， 37（O4）： 52 - 56.

[19]王哲彥.細胞珠蛋白互作蛋白的篩選、鑒定及過表達、敲除 前后人肝細胞轉錄組分析比較[D].南方醫科大學，2020. WANG Zheyan. Cell pearl protein interactions proteinbefore and after the screening，identification and expression，knock out person livercelltranscriptome analysis[D]. Southern Medical University，2020.

[20］吳越．鎘脅迫日本沼蝦肝胰腺轉錄組學及樁蛋白基因的表 達與功能分析[D]．杭州：浙江大學，2020. WU Yue. Transcriptomics and expression and function analysis of hepatopancreas and post protein genes of Macrobrachium nipponense under cadmium stress[D].Hangzhou：Zhejiang University， 2020.

[21]Loconsole G，SaldarelliP，Doddapaneni H，etal. Identification of a single-stranded DNA virus associated with Citrus chlorotic dwarf disease，a new member in the family Geminiviridae[J].Virology，2012，432（1）：162-172.

[22]Wang G，Sun Y W，Xu R R， et al.DNA -A of a highly pathogenic Indian cassava mosaic virus isolated from Jatropha curcas causes symptoms in Nicotiana benthamiana[J]. Virus Genes， 2014，48（2）：402-405.

[23]MaYX，NavarroB，Zhang ZX，et al.Identification and molecular characterization of a novel monopartite geminivirus associated with mulberry mosaic dwarf disease[J].The Journal of General Virology，2015，96（8）：2421-2434.

Abstract：【Objective】 This project aims to determine the field infection status of Grapevine Geminivirus A in Turpan City，Xinjiang.【Methods】71 suspected grape virus samples were collected and detected and according to the results of high -throughput sequencing，the existence of grapevine geminivirus A in the sample was verified，and the whole gene sequence of GGVA was determined. The whole gene sequence of grapevine geminivirus A was used as the target to design specific primers for the detection of suspected grape virus samples.【Results】 The whole genome sequence of GGVA was determined，and GGVA -QKT isolates were obtained，with nucleotide identity of 94.3%-99.3% ，and published isolates were highly homologous，among which，they were closely related to two GGVA isolates from Israel（KX618694.1） and（NC - 031340.1） . 【Conclusion】71 diffrent parts of the grape samples were Grapevine geminivirus A，71 grape samples from diferent regions were infected with Grapevine geminivirus A.among them， Shanshan County，Xinjiang had the highest detection rate 30.55% ，while Tokson County，Xinjiang had the lowest detection rate 9.10% . The virus was very common in the grape growing area of Turpan County，Xinjiang.In thisstudy，the field infection status of Grapevine geminivirus A in Turpan City is detected.

Key words：grape virus；DNA virus;sequence analysis;detection