侵染花葉青木的病毒種類鑒定及其基因序列分析

梁海文，蘭平秀，柳勤海，譚冠林，陳小姣，趙雁，李凡

侵染花葉青木的病毒種類鑒定及其基因序列分析

梁海文1,2，蘭平秀1，柳勤海1，譚冠林3，陳小姣1，趙雁2，李凡1

1云南農業大學植物保護學院，昆明 650201；2云南農業大學園林園藝學院，昆明 650201；3云南農業大學現代教育技術中心，昆明 650201

【目的】探明采自云南和海南表現為黃化、皺縮等表型的花葉青木病毒種類?！痉椒ā渴紫?，利用轉錄組測序技術對來自云南省昆明市的花葉青木樣品進行檢測；然后，根據轉錄組測序結果，采用RT-PCR對采自云南省昆明市和海南省?？谑械?3份花葉青木樣品開展轉錄組測序結果中發現的病毒種類檢測驗證，并擴增獲得相關病毒基因片段，對這些病毒的基因序列進行比對分析。【結果】轉錄組測序結果發現，送測的花葉青木樣品中含有蠶豆萎蔫病毒2號（broad bean wilt virus 2，BBWV2）、菜豆黃花葉病毒（bean yellow mosaic virus，BYMV）、黃瓜花葉病毒（cucumber mosaic virus，CMV）、豌豆耳突花葉病毒1號（pea enation mosaic virus 1，PEMV-1）和豌豆耳突花葉病毒2號（pea enation mosaic virus 2，PEMV-2）5種病毒的侵染。但隨后對來自云南和海南的73份花葉青木樣品的RT-PCR驗證中，只檢測到BBWV2、BYMV、CMV和PEMV-2 4種病毒。其中，云南和海南的樣品中均能檢測到BBWV2，總檢出率為11.0%；而BYMV、CMV和PEMV-2僅在云南的樣品中被檢測到，檢出率分別為5.5%、1.4%和1.4%。從兩地采集的花葉青木樣品中既未發現PEMV-1的侵染，也未發現BBWV2、BYMV、CMV和PEMV-2的復合侵染。為了進一步分析BBWV2、BYMV、CMV和PEMV-2花葉青木分離物與來自其他寄主植物相關病毒分離物的分子變異，選取檢測到上述4種病毒的樣品，分別對BBWV2、BYMV和CMV 3種病毒的及PEMV-2的序列進行擴增，并對擴增獲得的BBWV2 490 nt的small（GenBank登錄號：OQ137565、OQ137566）、BYMV 499 nt的核心區域（GenBank登錄號：OQ137568）、CMV 880 nt的（GenBank登錄號：OQ137567）和PEMV-2 800 nt的核心區域序列（GenBank登錄號：OQ137569）進行分析。序列比對分析結果表明，BBWV2云南與海南花葉青木分離物間存在85.9%的核苷酸序列一致性，二者與GenBank中其他BBWV2分離物分別存在79.8%—94.5%和80.2%—92.0%的核苷酸序列一致性；BYMV、CMV和PEMV-2云南花葉青木分離物與GenBank中相應病毒的核苷酸序列一致性分別為79.8%—99.0%、71.5%—99.4%和84.9%—98.0%。利用以上病毒相應基因核苷酸序列構建系統發育樹，發現BBWV2云南和海南花葉青木分離物分別屬于I組中的a亞組和b亞組；BYMV云南花葉青木分離物與伊朗（GenBank登錄號：MN241051，MN241060）和伊拉克（GenBank登錄號：JQ026005）蠶豆分離物具有較近的親緣關系；CMV云南花葉青木分離物屬于I組中的IB亞組，但與IB亞組的其他成員具有較遠的親緣關系；PEMV-2云南花葉青木分離物與云南豌豆分離物親緣關系最近并聚為一組，但與其他分離物存在較遠的親緣關系?！窘Y論】首次報道了花葉青木被病毒感染，同時也是BBWV2、BYMV、CMV和PEMV-2在花葉青木乃至桃葉珊瑚屬植物上的首次報道，BBWV2、BYMV、CMV和PEMV-2花葉青木分離物與來自其他植物的相關病毒分離物存在較大的分子變異。

花葉青木；病毒鑒定；蠶豆萎蔫病毒2號；菜豆黃花葉病毒；黃瓜花葉病毒；豌豆耳突花葉病毒2號

0 引言

【研究意義】花葉青木（var.）為絲纓花科（Garryaceae）桃葉珊瑚屬（）常綠灌木，其葉片革質光亮且具有金黃色斑點，是青木（）最常見的栽培變種，在國內外園林綠化中廣泛應用。2021年，筆者發現昆明市園林綠化中的花葉青木葉片大多具有黃化斑點等表型，有的甚至整個植株葉片都布滿黃色斑紋。這些花葉青木葉片的黃色斑紋表型與其原種青木葉片的表型差別很大，同時還觀察到花葉青木“色彩斑斕”的葉片表型與植物病毒侵染引起的癥狀較為相似，而花葉青木主要通過扦插方式進行繁殖，極易導致病毒等病原物通過繁殖材料傳播擴散。目前尚不清楚昆明市園林綠化中這些表現為黃化斑點的表型是花葉青木固有特性，還是因病毒感染所致。為此，通過對這些表現為黃化斑點等癥狀的花葉青木樣品開展病毒種類鑒定及其分子變異分析，以期為花葉青木葉片表型與病毒侵染的相關性研究，以及彩葉觀賞植物的品種選育等提供理論依據，同時也將為觀賞植物病毒病的防控工作提供參考。【前人研究進展】多年生的園林觀賞植物很容易遭受各種病原物的侵害，對其觀賞價值等造成嚴重影響，但有些病毒的侵染對于觀賞植物幾乎沒有危害，反而使其葉片和花朵的色彩變得更加豐富，進而增加其觀賞價值和經濟價值[1]。當今，彩葉植物的觀賞價值遠高于同源的綠葉觀賞植物，在園林綠化中使用頻率高，經濟價值也隨之提升。因此在觀賞植物育種過程中，人們有意識地選擇具有花葉或雜色花性狀的品種并廣泛引種栽培。有研究表明觀賞植物豐富表型的產生與病毒侵染相關，如段永嘉等[2]在香石竹、唐菖蒲、紫羅蘭及虞美人的雜色花品種上分別檢測出香石竹斑駁病毒（carnation mottle virus，CarMV）、菜豆黃花葉病毒（bean yellow mosaic virus，BYMV）、蕪箐花葉病毒（turnip mosaic virus，TuMV），探明了這4種植物的花瓣雜色均為病毒侵染所致。吳康承[3]從具有花葉癥狀的彩葉扶桑上分離得到的番茄花葉病毒（tomato mosaic virus，ToMV），將其接種普通朱槿后可產生全株花葉的朱槿。因此，有些植物病毒的侵染與觀賞植物的彩色葉片或花瓣有關，有的病毒在提高觀賞植物的觀賞價值方面是有益的。截至目前，與花葉青木同源的青木上僅見日本學者報道過花椰菜花葉病毒科（）桿狀DNA病毒屬（）的桃葉珊瑚環斑病毒（aucuba ringspot virus，AuRV）[4]和伴生豇豆病毒科（）線蟲傳多面體病毒屬（）的蘇鐵壞死矮化病毒（cycas necrotic stunt virus，CNSV）[5]兩種病毒的侵染，尚未見有關花葉青木病毒方面的報道。深度測序（deep sequencing）也稱下一代測序（next-generation sequencing，NGS）和高通量測序（high-throughput sequencing，HTS），是一種可在短時間內以較低成本對幾十萬到幾百萬條DNA分子進行序列測定的便捷測序技術，目前已在植物病毒鑒定和檢測方面廣泛應用，可鑒定已知病毒、發現新病毒及分析病毒間的進化關系等[6]。目前宏基因組測序（metagenome sequencing）、轉錄組測序（RNA sequencing，RNA-seq）和小RNA測序（small RNA sequencing，sRNA-seq）是深度測序技術中主要的病毒檢測鑒定方法，其中轉錄組測序在植物病毒研究中應用最為廣泛。Lan等通過去除rRNA的轉錄組測序，分別在表現為壞死、褪綠癥狀的山藥樣品中發現柘橙病毒屬（）的一種新病毒——山藥褪綠壞死病毒（yam chlorotic necrosis virus，YCNV）[7]，在諾麗上發現馬鈴薯Y病毒屬（）的一種新病毒——諾麗花葉病毒（noni mosaic virus，NoMV）[8]；李梁利用轉錄組測序技術發現在表現為皺縮癥狀的藜上存在蠶豆萎蔫病毒2號（broad bean wilt virus 2，BBWV2）的侵染[9]。轉錄組測序也被用于分析病毒侵染后對觀賞植物葉色變化的影響，如王煒等利用轉錄組測序探明滇山茶花葉呈色可能與病毒侵染后誘導了多個代謝途徑的基因沉默有關[10]。【本研究切入點】彩葉植物的葉片呈色機理受多方面的影響，有的可能為品種本身遺傳特性決定，有的或是被病毒感染所致。然而，病毒感染對園林植物觀賞性方面的影響迄今鮮有報道，也尚未見有關花葉青木被病毒感染的報道。【擬解決的關鍵問題】利用轉錄組測序結合RT-PCR技術，對采自云南和海南表現為黃化、皺縮和褪綠斑等癥狀的花葉青木樣品進行病毒種類檢測鑒定，并對其中發現的病毒進行分子變異和系統發育分析。

1 材料與方法

1.1 樣品采集

2021年8月，從云南省昆明市部分住宅小區景觀花園中采集表現為黃化（圖1-a）和皺縮、褪綠斑（圖1-b）的花葉青木葉片樣品58份，挑選具有典型癥狀的7份樣品等量混合后于液氮中磨碎，送上海元莘生物醫藥科技有限公司進行轉錄組測序（RNA-seq），其余樣品（包括轉錄組測序剩余的7份）于-20 ℃冰箱保存備用。

2022年2月，從海南省?？谑胁杉憩F為不規則黃斑和黃化（圖1-c）的花葉青木樣品15份，-20 ℃冰箱保存備用。

1.2 材料與試劑

DH5菌株為本實驗室保存，CTAB為國產分析純，pMD19-T Vector和Reverse Trancsriptase M-MLV（RNase H-）購自寶生物工程（大連）有限公司，DNA Polymerase和DNA純化回收試劑盒（離心柱型）購自天根生化科技（北京）有限公司。

1.3 轉錄組測序

轉錄組測序樣品的總核酸采用Invitrogen公司的Trizol試劑進行提取后，運用NanoDrop 2000（Wilmington，DE，USA）分光光度計對核酸質量進行檢測，將OD260/280在1.6—1.8的核酸利用TruSeq RNA Sample Prep Kit除去rRNA，使用Illumina HiSeq X-ten系統進行測序。所得序列經CLC Genomic Workbench 9.5（QIAGEN）分析后，組裝成200 bp以上的contigs，利用NCBI中的Blastn和Blastx進行比對分析并獲取與病毒相關的contigs。數據分析由美國農業部農業研究服務署（USDA-ARS）國家種質資源實驗室的李如慧博士幫助完成。

1.4 引物設計

根據轉錄組測序獲取的contigs設計檢測BYMV、PEMV-1和PEMV-2的特異引物，BBWV2和CMV的檢測引物參考文獻[11]，引物詳細信息見表1。

表1 花葉青木病毒檢測引物

1.5 總核酸提取、RT-PCR及PCR產物克隆

采用Li等[12]的CTAB法提取樣品總核酸。以總核酸為模板，按照Reverse Transcriptase M-MLV（RNase H-）試劑盒說明書進行反轉錄，獲取病毒的cDNA。使用DNA Polymerase試劑盒進行PCR擴增，10 μL反應體系：10×Buffer（Mg2+plus）1 μL、dNTP mix（2.5 mmol·L-1）0.8 μL、上下游引物各0.2 μL、DNA Polymerase 0.2 μL、ddH2O 6.6 μL、cDNA 1 μL，輕彈管壁混勻后進行以下反應：94 ℃預變性2 min，94 ℃變性30 s，57 ℃退火30 s，72 ℃延伸若干30 s，擴增35次循環，72 ℃延伸10 min后，取全部PCR產物用1%瓊脂糖凝膠于0.5×TBE緩沖液中電泳，35 min后于凝膠成像系統中觀察擴增結果。

紫光燈下切下與引物設計相符的擴增片段，經SanPrep柱式膠回收試劑盒純化后，連接至pMD19-T Vector，10 h后轉化至大腸桿菌DH5，在含Amp的LB固體培養基上培養12 h后，挑取單菌落擴繁10 h并進行PCR驗證，挑選陽性克隆送至生工生物工程（上海）股份有限公司進行測序。

1.6 序列分析

利用NCBI網站中的Nucleotide BLAST（https:// blast.ncbi.nlm.nih.gov/Blast.cgi）對測序結果進行比對分析，從GenBank數據庫下載BBWV2、BYMV、CMV和PEMV-2的序列，利用Clustal Omega（https://www. ebi.ac.uk/Tools/msa/clustalo/）進行多重比對并截取與本研究獲取序列一致的區域，刪去相互間序列變異小于1%的分離物，利用MEGA7.0軟件中的最大似然法（maximum likelihood，ML）分別構建BBWV2、BYMV、CMV和PEMV-2等病毒相關基因的系統發育樹，其中構建系統發育樹所使用的最適替代模型分別為TN93（Tamura-Nei）、T92（Tamura 3- parameter）、HKY（Hasegawa-Kishino-Yano）和K2（Kimura 2-parameter）。

2 結果

2.1 花葉青木樣品轉錄組測序

通過對7份花葉青木樣品混合物的轉錄組測序，共獲得67 563 492的reads（讀長），去除reads中的接頭序列及由于接頭自連等原因導致沒有插入片段的reads，組裝后獲得長度大于200 bp的187 602個contigs（重疊群）。Blastn分析顯示，長度532—5 297 bp的幾個contigs分別與BBWV2、BYMV、CMV、PEMV-1和PEMV-2等病毒擁有94%—99%最高的核苷酸序列一致性（表2），暗示花葉青木樣品中可能存在以上5種病毒的侵染。

表2 利用轉錄組測序從花葉青木混合樣中獲取的病毒相關數據統計

2.2 花葉青木單一樣品中轉錄組測序相關病毒的RT-PCR驗證

為驗證轉錄組測序獲得的BBWV2、BYMV、CMV、PEMV-1和PEMV-2這5種病毒在花葉青木樣品中的真實性，利用上述5種病毒的特異性檢測引物（表1），對送轉錄組測序的7份采自昆明的花葉青木單一樣品進行相關病毒的RT-PCR檢測，其中4份樣品能分別擴增到約490、499、880和800 nt與引物設計相符的病毒基因片段（圖2）。將擴增產物進行克隆和測序，Blastn分析顯示這些基因片段分別與GenBank中的BBWV2（OP722593）、BYMV（JQ026005）、CMV（MT786689）和PEMV-2（MN399704）擁有95%、99%、99%和96%的核苷酸序列一致性，證明花葉青木中存在BBWV2、BYMV、CMV和PEMV-2這4種病毒的感染，7份花葉青木單一樣品均未檢測到PEMV-1的感染。

2.3 云南和海南花葉青木樣品中BBWV2、BYMV、CMV、PEMV-1和PEMV-2的檢測結果

利用BBWV2、BYMV、CMV、PEMV-1和PEMV-2的特異檢測引物，對采自云南（58份，含送轉錄組測的7份）以及海南（15份）花葉青木樣品進行上述5種病毒的RT-PCR檢測，共有14份樣品檢測到BBWV2、BYMV、CMV和PEMV-2，病毒總檢出率為19.2%（表3）。其中檢出BBWV2的樣品數為8份，檢出率最高，為11.0%；其次為BYMV，樣品檢出數為4份，檢出率為5.5%；CMV和PEMV-2的樣品檢出率均較低，分別僅有1份樣品檢測到陽性條帶，檢出率為1.4%。BBWV2在云南和海南的樣品中均被檢測到，而BYMV、CMV和PEMV-2僅在云南的樣品中被檢測到。73份花葉青木樣品中未發現PEMV-1的侵染，也未發現BBWV2、BYMV、CMV和PEMV-2等病毒之間的復合侵染。

表3 云南和海南73份花葉青木樣品BBWV2、BYMV、CMV、PEMV-1和PEMV-2的RT-PCR檢測結果

-：未檢出病毒no virus was detected

1—7：花葉青木樣品A. japonica var. variegata sample；M：Marker；CK：陰性對照Negative control

2.4 BBWV2花葉青木分離物small cp的遺傳變異及系統發育分析

結合RT-PCR檢測結果，從云南和海南的BBWV2陽性花葉青木樣品中各選擇一份樣品進行BBWV2擴增產物的測序，序列分析顯示獲取的490 nt片段為BBWV2的small序列。序列比對結果顯示BBWV2花葉青木云南分離物（BBWV2-YNHYQM，GenBank登錄號：OQ137565）與海南分離物（BBWV2-HNHYQM，GenBank 登錄號：OQ137566）間具有85.9%的核苷酸序列一致性，且這兩個BBWV2花葉青木分離物與GenBank中的其他BBWV2分離物分別具有79.8%—94.5%和80.2%—92.0%的核苷酸序列一致性。

從GenBank中獲取BBWV2不同分離物的small基因序列，去除相互間序列差異小于1%的分離物，對來自中國、日本、韓國、英國、德國等不同國家和地區的番茄、辣椒、白術、山藥、菊三七、車前草、本氏煙、蠶豆等不同寄主植物22個BBWV2分離物構建系統發育樹，結果顯示這些BBWV2分離物可分為I、II兩組，I組又可分為a、b兩個亞組，其中本文獲得的兩個BBWV2花葉青木分離物YNHYQM和HNHYQM分別分布于I組中的a亞組和b亞組，分別與來自德國的本氏煙分離物（OK058514）和來自云南的辣椒分離物（MW271032）具有較近的親緣關系（圖3）。

2.5 BYMV花葉青木分離物cp核心區的遺傳變異及系統發育分析

根據RT-PCR檢測結果，對唯一一個含有BYMV感染的花葉青木樣品進行病毒擴增序列分析，顯示獲取的499 nt片段為BYMV核心區（BYMV-YNHYQM，GenBank登錄號：OQ137568）。從GenBank下載BYMV其他分離物相應區域序列進行多重比對，發現BYMV花葉青木分離物YNHYQM與其他分離物在核心區域具有79.8%—99.0%的核苷酸序列一致性，分子變異較大。去除相互間序列差異小于1%的分離物，利用來自中國、日本、韓國、印度、英國、美國、阿根廷、伊朗、伊拉克等不同國家和地區豆科、蘭科等不同作物上的20個BYMV分離物構建系統發育樹，顯示這些BYMV分離物可分為6個組，其中YNHYQM分離物與來自伊朗（MN241051，MN241060）和伊拉克（JQ026005）的蠶豆分離物具有較近的親緣關系，聚集在一起歸為V組（圖4）。各組中，來自豆科、蘭科等相近寄主的分離物分別顯示了相對較近的親緣關系。

分離物GenBank登錄號/寄主植物/地理來源GenBank accession number of different isolates/host plant/geographical origin；p：本研究獲得的分離物The isolate obtained in this study；置信值：分支上的數值，用來評估該分支的可信度，置信值一般＞70%時為可信The value on the branch is used to evaluate the credibility of the branch. The confidence value is generally greater than 70%；距離標尺：序列間差異數值的單位長度，發育樹的比例尺The unit length of the value of the difference between sequences, the scale of the evolutionary tree。下同The same as below

2.6 CMV花葉青木分離物cp的遺傳變異及系統發育分析

根據RT-PCR檢測結果，對本文獲得的唯一一個CMV花葉青木分離物880 nt序列進行分析，所獲得的為CMV的（CMV-YNHYQM，GenBank登錄號：OQ137567）。從GenBank下載CMV其他分離物相應區域序列進行多重比對，CMV花葉青木分離物YNHYQM顯示了與其他分離物有71.5%—99.4%的核苷酸序列一致性。去除相互間序列差異小于1%的分離物，利用分布于IA、IB和II組的26個CMV分離物構建發育樹，YNHYQM分離物屬于IB亞組，但單獨成一簇，與IB亞組的其他成員顯示了較遠的親緣關系（圖5）。

2.7 PEMV-2花葉青木分離物RdRp核心區的遺傳變異及系統發育分析

根據RT-PCR檢測結果，對獲得的唯一一個PEMV-2花葉青木分離物800 nt的序列進行分析，顯示為的核心區域（PEMV-2-YNHYQM，GenBank登錄號：OQ137569）。從GenBank下載PEMV-2其他分離物相應區域進行多重序列比對，顯示PEMV-2花葉青木分離物YNHYQM與其他PEMV-2分離物擁有84.9%—98.0%的核苷酸序列一致性。去除相互間序列差異小于1%的PEMV-2分離物構建系統發育樹，發現這19個PEMV-2分離物可分為3個組，大多數的分離物聚集在I組；而II組又可分為a亞組和b亞組，其中德國豌豆分離物（MW961148）單獨形成II組的a亞組，PEMV-2云南花葉青木分離物與云南的4個豌豆分離物聚集在II組的b亞組，親緣關系較近；美國豌豆分離物（JF723436）與其他分離物親緣關系較遠，單獨形成亞組III（圖6）。

豆科 Leguminosae；蘭科 Orchidaceae

3 討論

3.1 云南和海南花葉青木的病毒種類

花葉青木是常見的園林植物，國內外報道過花葉青木上有炭疽病、白粉病、根腐病、褐斑病等病害的發生[13-17]，但尚未見病毒病的相關報道。本研究利用轉錄組測序從采自昆明的疑似病毒病癥狀花葉青木樣品中發現可能有BBWV2、BYMV、CMV、PEMV-1和PEMV-2感染。隨后采用RT-PCR從供試的73份花葉青木樣品中只檢測到BBWV2、BYMV、CMV和PEMV-2，表明花葉青木中存在上述4種病毒的侵染。但4種病毒中只有BBWV2同時在云南和海南的樣品中被檢測到，而BYMV、CMV和PEMV-2僅在云南的樣品中被檢測到。自然界中，PEMV-1與PEMV-2通常復合侵染引起豌豆耳突花葉病，但植株中也存在僅檢測到PEMV-2而檢測不到PEMV-1的現象[18-19]。本研究轉錄組測序結果顯示獲得了PEMV-1的近全長基因組序列，但利用RT-PCR技術在73份花葉青木樣品中均未檢測到PEMV-1，原因有待進一步探究。此前日本科學家曾報道在與花葉青木同源的青木上發現有AuRV和CNSV侵染[4-5]。為此，本研究也利用這兩種病毒的特異引物對供試的73份花葉青木樣品進行RT-PCR檢測，但未擴增到相關病毒的目標條帶，因此供試的花葉青木樣品中暫未發現含有AuRV和CNSV這兩種病毒。

3.2 4種病毒花葉青木分離物與其他分離物的系統發育關系

BBWV2在全世界范圍內均有發生和分布，在44個科186屬328種植物上曾有報道[20-22]。本研究中BBWV2是唯一在云南與海南花葉青木樣品上均檢測到的病毒，但來自這兩個地區的BBWV2陽性花葉青木樣品上觀察到的癥狀表現不一致，云南的陽性樣品表現為邊緣不清晰的淡黃色斑塊，而海南的陽性樣品嫩葉有邊緣明顯的黃色斑點，老葉黃化。兩個BBWV2花葉青木分離物在small區域具有85.9%的核苷酸序列一致性，系統發育上分屬于不同亞組，親緣關系較遠。尚不清楚這些差異是否是造成云南和海南兩地花葉青木表型差異的主要原因。另外，本實驗室也在花葉青木BBWV2陽性樣品中檢測到1種新病毒（結果未發表），因此云南和海南兩地BBWV2陽性樣品葉片癥狀不一的原因尚無法確定。BYMV是蠶豆和豌豆上發生最普遍和最嚴重的病毒之一，可侵染12個科的114種植物[23-24]，本研究中獲得BYMV花葉青木分離物與伊朗（MN241051和MN241060）和伊拉克（JQ026005）蠶豆分離物的親緣關系較近，但更翔實的數據還有待BYMV花葉青木分離物全基因組序列的獲取及進一步分析。CMV被稱為全球第四大植物病毒，可侵染100多個科1 000多種植物并造成危害和嚴重經濟損失[25-26]。根據癥狀特點、寄主范圍、致病性及核苷酸序列同源性的差異，CMV可分為I和II兩個亞群[27]，亞群I進一步分為IA和IB兩個亞組，亞群II和IA亞組全球分布，而IB亞組主要在亞洲發生[28]。本研究在花葉青木上發現的CMV分離物屬于IB亞組，與IB亞組的印度和烏干達分離物最近，但與國內其他CMV分離物不在同一組，暗示CMV花葉青木分離物可能是從國外傳入。另外，PEMV-2花葉青木分離物在系統發育中與本實驗室獲得的4個PEMV-2云南豌豆分離物（數據未發表）單獨聚成一組，與其他國家PEMV-2豌豆分離物的親緣關系較遠，推測本研究獲得的PEMV-2花葉青木分離物與云南豌豆中PEMV-2分離物可能有著相同的起源。不過，本研究中花葉青木4種病毒的來源還需進行更多深入的探究。

圖5 基于CMV不同分離物cp構建的系統發育樹

圖6 基于PEMV-2不同分離物RdRp核心區構建的系統發育樹

3.3 花葉青木葉片表型與病毒種類的相關性

在花葉青木上發現的4種病毒中，海南的樣品只檢測到BBWV2一種病毒，而云南的樣品BBWV2、BYMV、CMV和PEMV-2這4種病毒均被檢測到，可看出云南和海南的花葉青木樣品上的病毒種類及數量差異較大，云南的病毒種類多樣性更豐富。暫不清楚造成上述差異的原因，推測可能是兩地栽培的花葉青木種質來源不同或是其周圍植物及其攜帶的病毒種類各異、氣候條件差異等，導致兩地花葉青木上的病毒種類不一，致使兩地的花葉青木葉片斑紋也不盡相同。本研究雖未在花葉青木樣品中檢測到AuRV和CNSV，但根據日本學者報道，青木受這兩種病毒侵染后均產生黃色環斑，這與本研究供試的花葉青木樣品表現的不規則黃色斑點表型相似，因此花葉青木葉片黃色斑紋表型的產生很有可能受病毒侵染影響，且病毒種類不一樣，在花葉青木上造成的表型可能相同，也可能不相同。除了本文中所報道的BBWV2、BYMV、CMV和PEMV-2這4種病毒外，本實驗室還在花葉青木上發現了多種新病毒（結果另文報道）。由于目前較難找到不含病毒的青木植株，且尚未分離出單一病毒分離物在花葉青木上進行回接試驗并驗證柯赫式法則，因此暫不能斷定花葉青木上葉片斑紋的變化是否與病毒侵染有關；若與病毒侵染有關，則病毒種類與花葉青木上葉片斑紋的相關性及其形成機制等都值得深入研究。

4 結論

采用高通量測序并結合RT-PCR技術，從來自云南和海南的73份花葉青木樣品中檢測到BBWV2、BYMV、CMV和PEMV-2 4種病毒，證明花葉青木中存在上述4種病毒的侵染，這是BBWV2、BYMV、CMV和PEMV-2在花葉青木乃至桃葉珊瑚屬植物上的首次報道。研究結果可為病毒侵染與花葉青木葉片表型相關性研究提供理論依據。

[1] 侯元凱. 彩葉植物研究進展. 世界林業研究, 2010, 23(6): 24-28.

HOU Y K. Study advances in color-leaf plants. World Forestry Research, 2010, 23(6): 24-28. (in Chinese)

[2] 段永嘉, 蔡紅. 花卉雜色花的研究. 植物病理學報, 1998, 28(1): 85-89.

DUAN Y J, CAI H. Studies on the variegated flowers of flowering plants. Acta Phytopathologica Sinica, 1998, 28(1): 85-89. (in Chinese)

[3] 吳康承. 病毒造就的繽紛植物. 生命世界, 2010(7): 52-55.

WU K C. Colorful plants created by viruses. Life World, 2010(7): 52-55. (in Chinese)

[4] UKE A, PINILI M S, NATSUAKI K T, GEERING A D W. Complete genome sequence of aucuba ringspot virus. Archives of Virology, 2021, 166(4): 1227-1230.

[5] KUSUNOKI M, YAMASHITA S, CHANG M U, DOI Y, YORA K. Isolation of new, cycas necrotic stunt virus from,and. Japan Journal Phytopathology, 1979, 45: 571-572.

[6] MATSUMURA E E, COLETTA-FILHO H D, NOURI S, FALK B W,

NERVA L, OLIVEIRA T S, DORTA S O, MACHADO M A. Deep sequencing analysis of RNAs from citrus plants grown in a citrus sudden death-affected area reveals diverse known and putative novel viruses. Viruses, 2017, 9(4): 92.

[7] LAN P X, MENG Y, SHEN P, LI R H, MA Y, TAN S T, CHEN H R, CAO M J, LI F. Complete genome sequence of yam chlorotic necrosis virus, a novel macluravirus infecting yam. Archives of Virology, 2018, 163(8): 2275-2278.

[8] LAN P X, HE P, ZHANG Y K, ZHANG S, ZHANG Z B, CHEN X J, TAN S T , LUO H M, CAO M J, LI F. Molecular characterization of a novel potyvirus infecting noni. Archives of Virology, 2019, 164(12): 3099-3102.

[9] 李梁. 衛矛與藜上病毒的檢測與鑒定[D]. 沈陽: 沈陽大學, 2020.

LI L. Detection and identification of viruses on theMaxin andL. plants[D]. Shenyang: Shenyang University, 2020. (in Chinese)

[10] 王煒, 鄭偉, 徐曉丹, 陳菁, 王談笑. 基于轉錄組測序的滇山茶花葉呈色機理分析. 西北植物學報, 2017, 37(9): 1720-1727.

WANG W, ZHENG W, XU X D, CHEN J, WANG T X. Coloring mechanism analysis of mosaic leaves inLindl. based on sequencing of transcriptome. Acta Botanica Boreali- Occidentalia Sinica, 2017, 37(9): 1720-1727. (in Chinese)

[11] 劉勇, 李凡, 李月月, 張松柏, 高希武, 謝艷, 燕飛, 張安盛, 戴良英, 程兆榜, 等. 侵染我國主要蔬菜作物的病毒種類、分布與發生趨勢. 中國農業科學, 2019, 52(2): 239-261. doi: 10.3864/j.issn. 0578-1752.2019.02.005.

LIU Y, LI F, LI Y Y, ZHANG S B, GAO X W, XIE Y, YAN F, ZHANG A S, DAI L Y, CHENG Z B,. Identification, distribution and occurrence of viruses in the main vegetables of China. Scientia Agricultura Sinica, 2019, 52(2): 239-261. doi: 10.3864/j.issn.0578- 1752.2019.02.005. (in Chinese)

[12] LI R, MOCK R, HUANG Q, ABAD J, HARTUNG J, KINARD G. A reliable and inexpensive method of nucleic acid extraction for the PCR-based detection of diverse plant pathogens. Journal of Virological Methods, 2008, 154: 48-55.

[13] 李翰文, 劉歡歡, 張敏, 蔡琴, 劉丹, 李沛利. 成都市灑金珊瑚炭疽病病原鑒定及潛在侵染源初探. 植物保護, 2021, 47(4): 38-45.

LI H W, LIU H H, ZHANG M, CAI Q, LIU D, LI P L. Identification of the pathogen causing anthracnose and preliminary investigation of the potential infection source ofin Chengdu. Plant Protection, 2021, 47(4): 38-45. (in Chinese)

[14] LI P L, LIU D, GONG G S, Chen S R, Yang X X. First report ofcausing anthracnose onin Sichuan province of China. Plant Disease, 2016, 100(5): 1019.

[15] CHO S E, ZHAO T, LEE S H, LEE S Y, SHIN H D. First report of powdery mildew caused byonin Korea. Plant Disease, 2017, 101(6): 1033.

[16] 馮福娟. 灑金東瀛珊瑚菌核性根腐病初步研究. 浙江林業科技, 2006, 26(3): 47-49, 58.

FENG F J. A preliminary study on the root rot ofvar.D’Om-Brain. Journal of Zhejiang Forestry Science and Technology, 2006, 26(3): 47-49, 58. (in Chinese)

[17] 呂傳紅, 楊銀虎. 灑金珊瑚特性及其培育技術. 現代園藝, 2015(7): 43-44.

Lü C H, YANG Y H. Characteristics and cultivation techniques ofvar.. Contemporary Horticulture, 2015(7): 43-44. (in Chinese)

[18] GAAFA Y Z A, HERZ K, HARTRICK J, FLETCHER J, BLOUIN A G, MACDIARMID R, ZIEBELL H. Investigating the pea virome in Germany—old friends and new players in the field(s). Frontiers in Microbiology, 2020, 11: 583242.

[19] FOWKES A R, MCGREIG S, PUFAL H, DUFFY S, HOWARD B, ADAMS I P, MACARTHUR R, WEEKES R, FOX A. Integrating high throughput sequencing into survey design reveals turnip yellows virus and soybean dwarf virus in pea () in the United Kingdom. Viruses, 2021, 13(12): 2530.

[20] EDWARDSON J R, CHRISTIC R G. Fabaviruses//CRC handbook of viruses infecting legumes. Boca Raton: CRC Press, 1991: 263-274.

[21] 周雪平, 余永杰, 劉勇, 李德葆. 侵染豌豆和蠶豆的蠶豆萎蔫病毒研究. 浙江農業大學學報, 1995, 21(3): 221-226.

ZHOU X P, YU Y J, LIU Y, LI D B. Studies on two isolates of broad bean wilt virus fromand. Journal of Zhejiang Agricultural University, 1995, 21(3): 221-226. (in Chinese)

[22] 青玲, 吳建祥, 戚益軍, 周雪平, 李德葆. 蠶豆萎蔫病毒單克隆抗體制備及檢測應用. 微生物學報, 2000, 40(2): 166-173.

QING L, WU J X, QI Y J, ZHOU X P, LI D B. Preparation of monoclonal antibodies to broad bean wilt virus and application in virus detection. Acta Microbiologica Sinica, 2000, 40(2): 166-173. (in Chinese)

[23] SOFY A R, HMED A, ALNAGGAR A, DAWOUD R A, ELSHAARAWY R, SOFY M R. Mitigating effects of bean yellow mosaic virus infection in faba bean using new carboxymethyl chitosan-titania nanobiocomposites.International Journal of Biological Macromolecules, 2020, 163: 1261-1275.

[24] 涂麗琴, 吳淑華, 干射香, 崔曉艷, 趙文浩, 程兆榜, 陳新, 周益軍, 季英華, 朱月林. 江蘇省蠶豆上菜豆黃花葉病毒的分子鑒定. 江蘇農業學報, 2019, 35(4): 804-810.

TU L Q, WU S H, GAN S X, CUI X Y, ZHAO W H, CHENG Z B, CHEN X, ZHOU Y J, JI Y H, ZHU Y L. Molecular identification of bean yellow mosaic virus infectingfrom Jiangsu province. Jiangsu Journal of Agricultural Sciences, 2019, 35(4): 804-810. (in Chinese)

[25] SCHOLTHOF K G, ADKINS S, CZOSNEK H, PALUKAITIS P, JACQUOT E, HOHN T, HOHN B, SAUNDERS K, CANDRESSE T, AHLQUIST P, HEMENWAY C, FOSTER G D. Top 10 plant viruses in molecular plant pathology. Molecular Plant Pathology, 2011, 12(9): 938-954.

[26] 宋嘉偉, 侯盼盼, 費莉玢, 龍奎, 劉穎, 高詩琪, 蘇秀, 馬良進. 黃瓜花葉病毒浙江白術分離物全基因組序列測定與分析. 植物病理學報, 2020, 50(2): 251-254.

SONG J W, HOU P P, FEI L B, LONG K, LIU Y, GAO S Q, SU X, MA L J. Complete genomic sequence analysis of cucumber mosaic virus isolated fromKoidz. in Zhejiang. Acta Phytopathologica Sinica, 2020, 50(2): 251-254. (in Chinese)

[27] JACQUEMOND M. Cucumber mosaic virus. Advances in Virus Research, 2012, 84: 439-504.

[28] 瞿青云, 沈文濤, 庹德財, 言普, 陳萍, 周鵬. 黃瓜花葉病毒香蕉分離物基因組全長克隆及侵染性克隆構建. 分子植物育種, https://kns.cnki.net/kcms/detail//46.1068.S.20221214.1136.002. html.

QU Q Y, SHEN W T, TUO D C, YAN P, CHEN P, ZHOU P. The genome full-length cloning and invasive cloning construction of cucumber mosaic virus banana isolate. Molecular Plant Breeding, https://kns.cnki.net/kcms/detail//46.1068.S.20221214.1136.002.html. (in Chinese)

Viruses Identification and Their Gene Sequences Analysis Infectingvar.

LIANG HaiWen1,2, LAN PingXiu1, LIU QinHai1, TAN GuanLin3, CHEN XiaoJiao1, ZHAO Yan2, LI Fan1

1College of Plant Protection, Yunnan Agricultural University, Kunming 650201;2College of Horticulture and Landscape, Yunnan Agricultural University, Kunming 650201;3Modern Education Technology Center, Yunnan Agricultural University, Kunming 650201

【Objective】The objective of this study is to identify the viruses infectingvar.plants with symptoms of yellowing and shrinking which collected from Yunnan and Hainan provinces.【Method】Firstly, RNA sequencing was used todetect the viruses infectingvar.plants from Kunming, Yunnan province. Then, based on the RNA sequencing results, RT-PCR was used to detect these viruses in 73var.field samplescollected in Kunming of Yunnan province and Haikou of Hainan province, and the related virus gene fragments were amplified, sequenced and analyzed.【Result】RNA sequencing results showed that 5 viruses were detected from the pooled diseased plants, such as broad bean wilt virus 2 (BBWV2), bean yellow mosaic virus (BYMV), cucumber mosaic virus (CMV), pea enation mosaic virus 1 (PEMV-1), and pea enation mosaic virus 2 (PEMV-2). However, only BBWV2, BYMV, CMV and PEMV-2 were detected in the 73var.samplesfrom Yunnan and Hainan. Among them, BBWV2 was detected in thevar.samples both from Yunnan and Hainan, with total detection rate of 11.0%, while BYMV, CMV and PEMV-2 were only detected in the Yunnan samples with detection rate of 5.5%, 1.4% and 1.4%, respectively. However, neither PEMV-1 nor the co-infections among BBWV2, BYMV, CMV and PEMV-2 was detected in these 73var.samples from Yunnan or Hainan. In order to further analyze the molecular variation of BBWV2, BYMV, CMV and PEMV-2 aucuba isolates with the related virus isolates from other host plants, the 490 nt of smallof BBWV2 (accession number OQ137565, OQ137566), 499 nt of the coreregion of BYMV (accession number OQ137568), 880 nt of theof CMV (accession number OQ137567), and 800 nt of the coreof PEMV-2 (accession number OQ137569) were amplified from thevar.samples infected with the above 4 viruses, and then the amplified gene sequences were analyzed, respectively. Sequence comparison analysis results showed that BBWV2 Yunnan aucuba isolate shared 85.9% nt identity with Hainan aucuba isolate, the BBWV2 Yunnan and Hainan aucuba isolates shared 79.8%-94.5% and 80.2%-92.0% nt identity with other BBWV2 isolates retrieved from GenBank, respectively. BYMV, CMV, and PEMV-2 aucuba isolates shared 79.8%-99.0%, 71.5%-99.4% and 84.9%-98.0% nt identity with other corresponding viruses in GenBank, respectively. Phylogenetic trees based on the gene nucleotide sequences of the above viruses revealed that the BBWV2 Yunnan and Hainan aucuba isolates were grouped into subgroup a and subgroup b in BBWV2 group I, respectively, and the BYMV Yunnan aucuba isolate showed a closest relationship with Iranian (accession number MN241051, MN241060) and Iraqi (accession number JQ026005) broad bean isolates. Moreover, the CMV Yunnan aucuba isolate was clustered into CMV IB subgroup of group I, but showed a distant relationship with other CMV isolates in the same subgroup. While PEMV-2 Yunnan aucuba isolate showed a closest relationship with Yunnan pea isolates and clustered into the same group, and had a distant relationship with other PEMV-2 isolates.【Conclusion】To our knowledge, this is the first record of virus infectingvar., and also the first report of BBWV2, BYMV, CMV and PEMV-2 infectingvar.as well as the plants in genus. Significant molecular variation was detected among the BBWV2, BYMV, CMV and PEMV-2 aucuba isolates with the related virus isolates from other plants worldwide.

var.; virus identification; broad bean wilt virus 2; bean yellow mosaic virus; cucumber mosaic virus; pea enation mosaic virus 2

10.3864/j.issn.0578-1752.2023.10.007

2023-02-01；

2023-03-14

國家重點研發計劃（2022YFD1400700）、云南省專家工作站（202005AF150040）

梁海文，E-mail：13627546002@163.com。蘭平秀，E-mail：lanpx22@126.com。梁海文和蘭平秀為同等貢獻作者。通信作者趙雁，E-mail：744234036@qq.com。通信作者李凡，E-mail：fanli@ynau.edu.cn

（責任編輯 岳梅）